Extended finite element-based cohesive zone method for modeling simultaneous hydraulic fracture height growth in layered reservoirs

In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.

A coupled Legendre-Laguerre polynomial method with analytical integration for the Rayleigh waves in a quasicrystal layered half-space with an imperfect interface

The Laguerre polynomial method has been successfully used to investigate the dynamic responses of a half-space.However,it fails to obtain the correct stress at the interfaces in a layered half-space,especially when there are significant differences in material properties.Therefore,a coupled Legendre-Laguerre polynomial method with analytical integration is proposed.The Rayleigh waves in a one-dimensional(1D)hexagonal quasicrystal(QC)layered half-space with an imperfect interface are investigated.The correctness is validated by comparison with available results.Its computation efficiency is analyzed.The dispersion curves of the phase velocity,displacement distributions,and stress distributions are illustrated.The effects of the phonon-phason coupling and imperfect interface coefficients on the wave characteristics are investigated.Some novel findings reveal that the proposed method is highly efficient for addressing the Rayleigh waves in a QC layered half-space.It can save over 99%of the computation time.This method can be expanded to investigate waves in various layered half-spaces,including earth-layered media and surface acoustic wave(SAW)devices.

THE NAVIER-STOKES EQUATIONS IN STREAM LAYER AND ON STREAM SURFACE AND A DIMENSION SPLIT METHODS

In this paper,we proposal stream surface and stream layer.By using classical tensor calculus,we derive 3-D Navier-Stokes Equations(NSE)in the stream layer under semigeodesic coordinate system,Navier-Stokes equation on the stream surface and 2-D Navier-Stokes equations on a two dimensional manifold. After introducing stream function on the stream surface,a nonlinear initial-boundary value problem satisfies by stream function is obtained,existence and uniqueness of its solution are proven.Based this theory we proposal a new method called"dimension split method"to solve 3D NSE.

Application and evaluation of layering shear method in LADCP data processing

The current velocity observation of LADCP(Lowered Acoustic Doppler Current Profiler)has the advantages of a large vertical range of observation and high operability compared with traditional current measurement methods,and is being widely used in the field of ocean observation.Shear and inverse methods are now commonly used by the international marine community to process LADCP data and calculate ocean current profiles.The two methods have their advantages and shortcomings.The shear method calculates the value of current shear more accurately,while the accuracy in an absolute value of the current is lower.The inverse method calculates the absolute value of the current velocity more accurately,but the current shear is less accurate.Based on the shear method,this paper proposes a layering shear method to calculate the current velocity profile by“layering averaging”,and proposes corresponding current calculation methods according to the different types of problems in several field observation data from the western Pacific,forming an independent LADCP data processing system.The comparison results have shown that the layering shear method can achieve the same effect as the inverse method in the calculation of the absolute value of current velocity,while retaining the advantages of the shear method in the calculation of a value of the current shear.

Seasonal variability of the mixed layer depth determined using an improved maximum angle method in the Arctic basins

To investigate the spatiotemporal variations in the mixed layer depth(MLD)in the Arctic basins,a new criterion to determine the MLD,called the improved maximum angle method(IMAM),was developed.A total of 45123 potential density profiles collected using Ice-Tethered Profilers(ITPs)in the Arctic basins during 2005-2021 were used to demonstrate the method’s effectiveness.By comparing the results obtained by the fixed threshold method(FTM),percentage threshold method(PTM),and maximum gradient method(MGM)for profiles in the Canada Basin,Makarov Basin,and Eurasian Basin,we determined that the quality index(1.0 for perfect identification of the MLD)of the IMAM regarding the assessment of the MLD determination method reached 0.94,which is much greater than those of other criteria.Moreover,two types of the density profiles were identified based on the mixed layer development stage.The MLDs of the typical profiles determined using the IMAM were found to have better consistency with the original definition.By utilizing the new mixed layer criterion,the seasonal variations and regional differences in the MLD in the Arctic basins were analyzed.Spatially,the summer and winter MLDs in the Canada Basin were the shallowest(13.55 m in summer,26.76 m in winter)than those in the Makarov(29.51 m in summer,49.08 m in winter)and Eurasian(20.36 m in summer,46.81 m in winter)basins due to the stable stratification in the upper ocean and the subsequent small effects of dynamic and thermodynamic processes(wind-driven stirring and brine rejection)in the Canada Basin.Seasonally,in the three Arctic basins,the average MLD was shallowest(22.77 m)in summer;it deepened through autumn and reached a winter maximum(41.12 m).

Advances in preparation methods and mechanism analysis of layered double hydroxide for lithium-ion batteries and lithium-sulfur batteries

Lithium-ion(Li-ion) battery and lithium-sulfur(Li-S) battery have attracted significant attention as promising components for large-scale energy storage because of high theoretical capacity of Li,excellent energy density or environmental friendness for two kinds of batteries.However,there still exist some respective obstacles for commercial applications,such as limited theoretical capacity,high cost and low conductivity of Li-ion cells or shuttle effect of lithium polysulfides of Li-S cells.As typical twodimensional materials,layered double hydroxides(LDHs) exhibit excellent potential in the field of energy storage due to facile tunability of composition,structure and morphology as well as convenient composite and strong catalytic properties.Consequently,various LDHs toward novel separators or interlayers,cathodes,anodes,and interesting catalytic templates are researched to resolve these challenges.In this review,the recent progress for LDHs applied in Li-ion batteries and Li-S batteries including the synthesis methods,designs and applications is presented and reviewed.Meanwhile,the existing challenges and future perspectives associated with material designs and practical applications of LDHs for these two classes of cells are discussed.WeWe hope that the review can attract more attention and inspire more profound researches toward the LDH-based electrochemical materials for energy storage.

Quantum dynamics within curved thin layers with deviation

Combining the deviation between thin layers' adjacent surfaces with the confining potential method applied to the quantum curved systems,we derive the effective Schr?dinger equation describing the particle constrained within a curved layer,accompanied by a general geometric potential V_(gq) composed of a compression-corrected geometric potential V_(gq)~*and a novel potential V_(gq)~(**) brought by the deviation.Applying this analysis to the cylindrical layer emerges two types of deviation-induced geometric potential,resulting from the the cases of slant deviation and tangent deviation,respectively,which strongly renormalizes the purely geometric potential and contribute to the energy spectrum based on a very substantial deepening of bound states they offer.

Experimental investigation of closed-loop active control to modulate coherent structures by mu-level method

The experimental research on zero-net-mass-flux jet closed-loop active control was conducted in the wind tunnel.The mu-level method successfully detected burst events of the coherent structures. The streamwise velocity signals in the turbulent boundary layer were measured by HWA. The drag reduction rate of 16.7% is obtained comparable to that of the open-loop control and saves 75% of the input energy at the asynchronous 100 V/160 Hz control case, which reflects the advantages of the closed-loop control. The experimental findings indicate that the intensity increases in the near-wall region.The perturbation of the PZT vibrators on the skewness factor is concentrated in the region y+< 60. The generation of highspeed fluids is depressed and the downward effect of high-speed fluids weakens. The alteration of energy distribution and the discernible impact of modulation between structures of varying scales are observed. The correlation coefficient exhibits a strong positive correlation, which indicates that the large-scale structures produce modulation effect on small-scale ones.The occurrence of burst events is effectively suppressed. The disturbance has the characteristics of stable periodicity,positive and negative symmetry, low intermittency, and high pulsation strength. The conditional phase waveform shows that the fluctuation amplitude increases, indicating amplitude modulation effects on coherent structures.

Probing the electric double layer structure at nitrogen-doped graphite electrodes by constant-potential molecular dynamics simulations

Electric double layer(EDL)is a critical topic in electrochemistry and largely determines the working performance of lithium batteries.However,atomic insights into the EDL structures on heteroatom-modified graphite anodes and EDL evolution with electrode potential are very lacking.Herein,a constant-potential molecular dynamics(CPMD)method is proposed to probe the EDL structure under working conditions,taking N-doped graphite electrodes and carbonate electrolytes as an example.An interface model was developed,incorporating the electrode potential and atom electronegativities.As a result,an insightful atomic scenario for the EDL structure under varied electrode potentials has been established,which unveils the important role of doping sites in regulating both the EDL structures and the following electrochemical reactions at the atomic level.Specifically,the negatively charged N atoms repel the anions and adsorb Li~+at high and low potentials,respectively.Such preferential adsorption suggests that Ndoped graphite can promote Li~+desolvation and regulate the location of Li~+deposition.This CPMD method not only unveils the mysterious function of N-doping from the viewpoint of EDL at the atomic level but also applies to probe the interfacial structure on other complicated electrodes.

Wavelet Multi-Resolution Interpolation Galerkin Method for Linear Singularly Perturbed Boundary Value Problems

In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.

Methods for Estimating the Ultimate Bearing Capacity of Layered Foundations

The Meyerhof and Hanna′s(M-H) method to estimate the ultimate bearing capacity of layered foundations was improved. The experimental results of the load tests in Tianjin New Harbor were compared with predictions with the method recommended by the code for the foundations of harbor engineering, i.e. Hansen′s method and the improved M-H method. The results of the comparisons implied that the code and the improved M-H method could give a better prediction.

Express Methods for Measurement of Electroconductivity of Semiconductor Layered Crystal

We describe theoretically the grounded method of measuring the conductivity of anisotropic layered semiconductor materials. The suggested method implies the use of a four-probe testing device with a linear arrangement of probes. The final expressions for identifying the electrical conductivity are presented in the form of a series of analytic functions. The suggested method is experimentally verified, and practical recommendations of how to apply it are also provided.

Development and Prospect of Process Models and Simulation Methods for Atomic Layer Deposition

Thin film deposition is one of the most important processes in IC manufacturing. In this paper, several typical models and numerical simulation methods for thin film deposition and atomic layer deposition are introduced. Several modeling methods based on the characteristics of atomic layer deposition are introduced, it includes geometric method, cellular automata and multiscale simulation. The principle of each model and simulation method is explained, and their advantages and disadvantages are analyzed. Finally, the development direction of thin film deposition and atomic layer deposition modeling is prospected, and some modeling ideas are also provided.

Growth of CuI buffer layer prepared by spraying method

CuI thin films with nano-scale grains of about 35nm were deposited via spraying method with using acetonitrile as solvent. The influence of iodine doping concentration in acetonitrile solution on the structure, topographic and optical properties of CuI thin films was investigated. X-ray diffraction results showed that CuI iodine-doped films doped CuI：I2 were in γ-phase of zinc blende structure with （111） preferential plane. Scanning electron microscopy revealed that the microstructure of CuI films depended on the relative amount of doping iodine in the solution. When the iodine doping amount in acetonitrile solution was 0.025 g, the film was uniform and compact, the optical transmittance was 75.4% in the part of visible region and the energy band gap was close to 2.96 eV.

Extension of e^N method to general three-dimensional boundary layers

In order to extend the eN method to general three-dimensional boundary layers, the conservation law of the imaginary parts for the wave parameters with a fixed wave vector is deduced. The compatibility relationship （CR） and the general theory of ray tracing （RT）, which have been extensively used in conservative systems, are applied to a general three-dimensional boundary layer belonging to non-conservative systems. Two kinds of eN methods, i.e., the eN-CR method and the eN-RT method, are established. Both the two kinds of methods can wavenumber and the amplitude of the be used to predict the evolutions of the spanwise disturbances in general three-dimensional boundary layers. The reliability of the proposed methods is verified and validated by performing a direct numerical simulation （DNS） in a hypersonic general three-dimensional boundary layer over an aircraft model. The results are also compared with those obtained by other eN methods, indicating that the proposed methods have great potential applications in improving the transition prediction accuracy in general three-dimensional boundary layers.

Improvement of X-Band Polarization Radar Melting Layer Recognition by the Bayesian Method and ITS Impact on Hydrometeor Classification

Using melting layer(ML)and non-melting layer(NML)data observed with the X-band dual linear polarization Doppler weather radar(X-POL)in Shunyi,Beijing,the reflectivity(ZH),differential reflectivity(ZDR),and correlation coefficient(CC)in the ML and NML are obtained in several stable precipitation processes.The prior probability density distributions(PDDs)of the ZH,ZDR and CC are calculated first,and then the probabilities of ZH,ZDR and CC at each radar gate are determined(PBB in the ML and PNB in the NML)by the Bayesian method.When PBB>PNB the gate belongs to the ML,and when PBB<PNB the gate belongs to the NML.The ML identification results with the Bayesian method are contrasUsing melting layer(ML)and non-melting layer(NML)data observed with the X-band dual linear polarization Doppler weather radar(X-POL)in Shunyi,Beijing,the reflectivity(ZH),differential reflectivity(ZDR),and correlation coefficient(CC)in the ML and NML are obtained in several stable precipitation processes.The prior probability density distributions(PDDs)of the ZH,ZDR and CC are calculated first,and then the probabilities of ZH,ZDR and CC at each radar gate are determined(PBB in the ML and PNB in the NML)by the Bayesian method.When PBB>PNB the gate belongs to the ML,and when PBB<PNB the gate belongs to the NML.The ML identification results with the Bayesian method are contrasted under the conditions of the independent PDDs and joint PDDs of the ZH,ZDR and CC.The results suggest that MLs can be identified effectively,although there are slight differences between the two methods.Because the values of the polarization parameters are similar in light rain and dry snow,it is difficult for the polarization radar to distinguish them.After using the Bayesian method to identify the ML,light rain and dry snow can be effectively separated with the X-POL observed data.ted under the conditions of the independent PDDs and joint PDDs of the ZH,ZDR and CC.The results suggest that MLs can be identified effectively,although there are slight differences between the two methods.Because the values of the polarization parameters are similar in light rain and dry snow,it is difficult for the polarization radar to distinguish them.After using the Bayesian method to identify the ML,light rain and dry snow can be effectively separated with the X-POL observed data.

MODIFIED LAYER REMOVAL METHOD FOR MEASUREMENT OF RESIDUAL STRESS DISTRIBUTION IN THICK PRE-STRETCHED ALUMINUM PLATE

The integrated structure parts are widely used in aircraft. The distortion caused by residual stresses in thick pre-stretched aluminum plates during machining integrated parts is a common and serious problem. To predict and control the machining distortion, the residual stress distribution in the thick plate must be measured firstly. The modified removal method for measuring residual stress in thick pre-stretched aluminum plates is proposed and the stress-strain relation matrix is deduced by elasticity theory. The residual stress distribution in specimen of 7050T7451 plate is measured by using the method, and measurement results are analyzed and compared with data obtained by other methods. The method is effective to measure the residual stress.

An efficient source wavefield reconstruction scheme using single boundary layer values for the spectral element method

In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.

Intercalation Assembly Method and Intercalation Process Control of Layered Intercalated Functional Materials

Layered intercalated functional materials of layered double hydroxide type are an important class of functional materials developed in recent years. Based on long term studies on these materials in the State Key Laboratory of Chemical Resource Engineering in Beiiing University of Chemical Technology, the orinciole for the design of controlled intercalation processes in the light of tuture production processing requirements has been developed. Intercalation assembly methods and technologies have been invented to control the intercalation process for preparing layered intercalated materials with various structures and functions.

Solidification process and microstructure of transition layer of Cu-Al composite cast prepared by method of pouring molten aluminum

The Cu?Al composite casts were prepared by the method of pouring molten aluminum. The solidification process and themicrostructure of the transition layer were investigated during the recombination process of the liquid Al and the solid Cu. The results reveal that the microstructure of the transition layer in the Cu?Al composite cast consists of α(Al)+α(Al)?CuAl2 eutectic,α(Al)?CuAl2 eutectic, CuAl2+α(Al)?CuAl2 eutectic and Cu9Al4. Additionally, the pouring temperature, cooling mode of the Cu platesurface and start time of the forced cooling after pouring have no effect on the microstructure species. But the proportion of thevarious microstructures in the transition layer changes with the process parameters. The pure Al at the top of the transition layer startsto solidify first and then the α(Al) phase grows in a dendritic way, while the CuAl2 phase exhibits plane or cellular crystal growth from the two sides of the transition layer towards its interior. The stronger the cooling intensity of the Cu plate outer surface, the more developed the dendrite, and the easier it is for the CuAl2 phase to grow into a plane crystal.