The similarities and differences in inherent mechanism and characteristic frequency between the onedimensional(1D)poroelastic model and the layered White model were investigated.This investigation was conducted under ...The similarities and differences in inherent mechanism and characteristic frequency between the onedimensional(1D)poroelastic model and the layered White model were investigated.This investigation was conducted under the assumption that the rock was homogenous and isotropic at the mesoscopic scale.For the inherent mechanism,both models resulted from quasi-static flow in a slow P-wave diffusion mode,and the differences between them originated from saturated fluids and boundary conditions.On the other hand,for the characteristic frequencies of the models,the characteristic frequency of the 1D poroelastic model was first modified because the elastic constant and formula for calculating it were misused and then compared to that of the layered White model.Both of them moved towards higher frequencies with increasing permeability and decreasing viscosity and diffusion length.The differences between them were due to the diffusion length.The diffusion length for the 1D poroelastic model was determined by the sample length,whereas that for the layered White model was determined by the length of the representative elementary volume(REV).Subsequently,a numerical example was presented to demonstrate the similarities and differences between the models.Finally,published experimental data were interpreted using the 1D poroelastic model combined with the Cole-Cole model.The prediction of the combined model was in good agreement with the experimental data,thereby validating the effectiveness of the 1D poroelastic model.Furthermore,the modified characteristic frequency in our study was much closer to the experimental data than the previous prediction,validating the effectiveness of our modification of the characteristic frequency of the 1D poroelastic model.The investigation provided insight into the internal relationship between wave-induced fluid flow(WIFF)models at macroscopic and mesoscopic scales and can aid in a better understanding of the elastic modulus dispersion and attenuation caused by the WIFF at different scales.展开更多
Bedding structural planes significantly influence the mechanical properties and stability of engineering rock masses.This study conducts uniaxial compression tests on layered sandstone with various bedding angles(0...Bedding structural planes significantly influence the mechanical properties and stability of engineering rock masses.This study conducts uniaxial compression tests on layered sandstone with various bedding angles(0°,15°,30°,45°,60°,75°and 90°)to explore the impact of bedding angle on the deformational mechanical response,failure mode,and damage evolution processes of rocks.It develops a damage model based on the Logistic equation derived from the modulus’s degradation considering the combined effect of the sandstone bedding dip angle and load.This model is employed to study the damage accumulation state and its evolution within the layered rock mass.This research also introduces a piecewise constitutive model that considers the initial compaction characteristics to simulate the whole deformation process of layered sandstone under uniaxial compression.The results revealed that as the bedding angle increases from 0°to 90°,the uniaxial compressive strength and elastic modulus of layered sandstone significantly decrease,slightly increase,and then decline again.The corresponding failure modes transition from splitting tensile failure to slipping shear failure and back to splitting tensile failure.As indicated by the modulus’s degradation,the damage characteristics can be categorized into four stages:initial no damage,damage initiation,damage acceleration,and damage deceleration termination.The theoretical damage model based on the Logistic equation effectively simulates and predicts the entire damage evolution process.Moreover,the theoretical constitutive model curves closely align with the actual stress−strain curves of layered sandstone under uniaxial compression.The introduced constitutive model is concise,with fewer parameters,a straightforward parameter determination process,and a clear physical interpretation.This study offers valuable insights into the theory of layered rock mechanics and holds implications for ensuring the safety of rock engineering.展开更多
To study the evolution of nanoparticles during Al wire electrical explosion,a nanoparticle formation model that considered layered motion was developed,and an experimental system was set up to carry out electrical exp...To study the evolution of nanoparticles during Al wire electrical explosion,a nanoparticle formation model that considered layered motion was developed,and an experimental system was set up to carry out electrical explosion experiments using 0.1 mm and 0.2 mm Al wires.The characteristic parameters and evolution process during the formation of nanoparticles were calculated and analyzed.The results show that the maximum velocities of the innermost and outermost layers are about 1200 m·s-1and 1600 m·s-1,and the velocity of the middle layer is about 1400 m·s-1,respectively.Most of the nanoparticles are formed in the temperature range of2600 K-2500 K.The characteristic temperature for the formation of Al nanoparticles is~2520K,which is also the characteristic temperature of other parameters.The size distribution range of the formed nanoparticles is 18 to 110 nm,and most of them are around 22 nm.The variation of saturated vapor pressure determines the temperature distribution range of particle nucleation.There is a minimum critical diameter of particles(~25 nm);particles smaller than the critical diameter can grow into larger particles during surface growth.Particle motion has an effect on the surface growth and aggregation process of particles,and also on the distribution area of larger-diameter particles.The simulation results are in good agreement with the experiments.We provide a method to estimate the size and distribution of nanoparticles,which is of great significance to understand the formation process of particles during the evolution of wire electrical explosion.展开更多
In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.Ho...In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.However,convolution operator is intractable in finite-difference time-domain(FDTD)method.A great deal of progress has been made in using time stepping instead of convolution in FDTD.To incorporate PML into viscoelastic media,more memory variables need to be introduced,which increases the code complexity and computation costs.By modifying the nonsplitting PML formulation,I propose a viscoelastic model,which can be used as a viscoelastic material and/or a PML just by adjusting the parameters.The proposed viscoelastic model is essentially equivalent to a Maxwell model.Compared with existing PML methods,the proposed method requires less memory and its implementation in existing finite-difference codes is much easier.The attenuation and phase velocity of P-and S-waves are frequency independent in the viscoelastic model if the related quality factors(Q)are greater than 10.The numerical examples show that the method is stable for materials with high absorption(Q=1),and for heterogeneous media with large contrast of acoustic impedance and large contrast of viscosity.展开更多
With continuous hydrocarbon exploration extending to deeper basins,the deepest industrial oil accumulation was discovered below 8,200 m,revealing a new exploration field.Hence,the extent to which oil exploration can b...With continuous hydrocarbon exploration extending to deeper basins,the deepest industrial oil accumulation was discovered below 8,200 m,revealing a new exploration field.Hence,the extent to which oil exploration can be extended,and the prediction of the depth limit of oil accumulation(DLOA),are issues that have attracted significant attention in petroleum geology.Since it is difficult to characterize the evolution of the physical properties of the marine carbonate reservoir with burial depth,and the deepest drilling still cannot reach the DLOA.Hence,the DLOA cannot be predicted by directly establishing the relationship between the ratio of drilling to the dry layer and the depth.In this study,by establishing the relationships between the porosity and the depth and dry layer ratio of the carbonate reservoir,the relationships between the depth and dry layer ratio were obtained collectively.The depth corresponding to a dry layer ratio of 100%is the DLOA.Based on this,a quantitative prediction model for the DLOA was finally built.The results indicate that the porosity of the carbonate reservoir,Lower Ordovician in Tazhong area of Tarim Basin,tends to decrease with burial depth,and manifests as an overall low porosity reservoir in deep layer.The critical porosity of the DLOA was 1.8%,which is the critical geological condition corresponding to a 100%dry layer ratio encountered in the reservoir.The depth of the DLOA was 9,000 m.This study provides a new method for DLOA prediction that is beneficial for a deeper understanding of oil accumulation,and is of great importance for scientific guidance on deep oil drilling.展开更多
Surface-wave inversion is a powerful tool for revealing the Earth's internal structure.However,aside from shear-wave velocity(v_(S)),other parameters can influence the inversion outcomes,yet these have not been sy...Surface-wave inversion is a powerful tool for revealing the Earth's internal structure.However,aside from shear-wave velocity(v_(S)),other parameters can influence the inversion outcomes,yet these have not been systematically discussed.This study investigates the influence of various parameter assumptions on the results of surface-wave inversion,including the compressional and shear velocity ratio(v_(P)/v_(S)),shear-wave attenuation(Q_(S)),density(ρ),Moho interface,and sedimentary layer.We constructed synthetic models to generate dispersion data and compared the obtained results with different parameter assumptions with those of the true model.The results indicate that the v_(P)/v_(S) ratio,Q_(S),and density(ρ) have minimal effects on absolute velocity values and perturbation patterns in the inversion.Conversely,assumptions about the Moho interface and sedimentary layer significantly influenced absolute velocity values and perturbation patterns.Introducing an erroneous Mohointerface depth in the initial model of the inversion significantly affected the v_(S) model near that depth,while using a smooth initial model results in relatively minor deviations.The assumption on the sedimentary layer not only affects shallow structure results but also impacts the result at greater depths.Non-linear inversion methods outperform linear inversion methods,particularly for the assumptions of the Moho interface and sedimentary layer.Joint inversion with other data types,such as receiver functions or Rayleigh wave ellipticity,and using data from a broader period range or higher-mode surface waves,can mitigate these deviations.Furthermore,incorporating more accurate prior information can improve inversion results.展开更多
Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve...Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve the simulation of the upper ocean.Thus far,several nonbreaking wave-induced mixing parameterization schemes have been proposed;however,no quantitative comparison has been performed among them.In this paper,a one-dimensional ocean model was used to compare the performances of five schemes,including those of Qiao et al.(Q),Hu and Wang(HW),Huang and Qiao(HQ),Pleskachevsky et al.(P),and Ghantous and Babanin(GB).Similar to previous studies,all of these schemes can decrease the simulated sea surface temperature(SST),increase the subsurface temperature,and deepen the mixed layer,thereby alleviating the common thermal deviation problem of the ocean model for upper ocean simulation.Among these schemes,the HQ scheme exhibited the weakest wave-induced mixing effect,and the HW scheme exhibited the strongest effect;the other three schemes exhibited roughly the same effect.In particular,the Q and P schemes exhibited nearly the same effect.In the simulation based on observations from the Ocean Weather Station Papa,the HQ scheme exhibited the best performance,followed by the Q scheme.In the experiment with the HQ scheme,the root-mean-square deviation of the simulated SST from the observations was 0.43℃,and the mixed layer depth(MLD)was 2.0 m.As a contrast,the deviations of the SST and MLD reached 1.25℃ and 8.4 m,respectively,in the experiment without wave-induced mixing.展开更多
The polyurethane foam(PU)compressible layer is a viable solution to the problem of damage to the secondary lining in squeezing tunnels.Nevertheless,the mechanical behaviour of the multi-layer yielding supports has not...The polyurethane foam(PU)compressible layer is a viable solution to the problem of damage to the secondary lining in squeezing tunnels.Nevertheless,the mechanical behaviour of the multi-layer yielding supports has not been thoroughly investigated.To fill this gap,large-scale model tests were conducted in this study.The synergistic load-bearing mechanics were analyzed using the convergenceconfinement method.Two types of multi-layer yielding supports with different thicknesses(2.5 cm,3.75 cm and 5 cm)of PU compressible layers were investigated respectively.Digital image correlation(DIC)analysis and acoustic emission(AE)techniques were used for detecting the deformation fields and damage evolution of the multi-layer yielding supports in real-time.Results indicated that the loaddisplacement relationship of the multi-layer yielding supports could be divided into the crack initiation,crack propagation,strain-hardening,and failure stages.Compared with those of the stiff support,the toughness,deformability and ultimate load of the yielding supports were increased by an average of 225%,61%and 32%,respectively.Additionally,the PU compressible layer is positioned between two primary linings to allow the yielding support to have greater mechanical properties.The analysis of the synergistic bearing effect suggested that the thickness of PU compressible layer and its location significantly affect the mechanical properties of the yielding supports.The use of yielding supports with a compressible layer positioned between the primary and secondary linings is recommended to mitigate the effects of high geo-stress in squeezing tunnels.展开更多
The dynamic responses and generated voltage in a curved sandwich beam with glass reinforced laminate(GRL)layers and a pliable core in the presence of a piezoelectric layer under low-velocity impact(LVI)are investigate...The dynamic responses and generated voltage in a curved sandwich beam with glass reinforced laminate(GRL)layers and a pliable core in the presence of a piezoelectric layer under low-velocity impact(LVI)are investigated.The current study aims to carry out a dynamic analysis on the sandwich beam when the impactor hits the top face sheet with an initial velocity.For the layer analysis,the high-order shear deformation theory(HSDT)and Frostig's second model for the displacement fields of the core layer are used.The classical non-adhesive elastic contact theory and Hunter's principle are used to calculate the dynamic responses in terms of time.In order to validate the analytical method,the outcomes of the current investigation are compared with those gained by the experimental tests carried out by other researchers for a rectangular composite plate subject to the LVI.Finite element(FE)simulations are conducted by means of the ABAQUS software.The effects of the parameters such as foam modulus,layer material,fiber angle,impactor mass,and its velocity on the generated voltage are reviewed.展开更多
An engineering system approach of 2-D cylindrical model of transient mass balance calculations of ozone and other concerned chemicals along with fourteen photolysis, ozone-generating and ozone-depleting chemical react...An engineering system approach of 2-D cylindrical model of transient mass balance calculations of ozone and other concerned chemicals along with fourteen photolysis, ozone-generating and ozone-depleting chemical reaction equations was developed, validated, and used for studying the ozone concentrations, distribution and peak of the layer, ozone depletion and total ozone abundance in the stratosphere. The calculated ozone concentrations and profile at both the Equator and a 60˚N location were found to follow closely with the measured data. The calculated average ozone concentration was within 1% of the measured average, and the deviation of ozone profiles was within 14%. The monthly evolution of stratospheric ozone concentrations and distribution above the Equator was studied with results discussed in details. The influences of slow air movement in both altitudinal and radial directions on ozone concentrations and profile in the stratosphere were explored and discussed. Parametric studies of the influences of gas diffusivities of ozone D<sub>O3</sub> and active atomic oxygen D<sub>O</sub> on ozone concentrations and distributions were also studied and delineated. Having both influences through physical diffusion and chemical reactions, the diffusivity (and diffusion) of atomic oxygen D<sub>O</sub> was found to be more sensitive and important than that of ozone D<sub>O3</sub> on ozone concentrations and distribution. The 2-D ozone model present in this paper for stratospheric ozone and its layer and depletion is shown to be robust, convenient, efficient, and executable for analyzing the complex ozone phenomena in the stratosphere. .展开更多
This study is an extension of the previous work done with ARS-680 Environmental Chamber. Drying is a complex operation that demands much energy and time. Drying is essentially important for preservation of ginger rhiz...This study is an extension of the previous work done with ARS-680 Environmental Chamber. Drying is a complex operation that demands much energy and time. Drying is essentially important for preservation of ginger rhizome. Drying of ginger was modeled, and then the effective diffusion coefficient and activation energy were determined. For this purpose, the experiments were done at six levels of varied temperatures: 10°C, 20°C, 30°C, 40°C, 50°C and 60°C. The values of effective diffusion coefficients obtained in this work for the variously treated ginger rhizomes closely agreed with the average effective diffusion coefficients of other notable authors who determined the drying kinetics and convective heat transfer coefficients of ginger slices.展开更多
In this paper, we study the long-time behavior of solutions of the single-layer quasi-geostrophic model arising from geophysical fluid dynamics. We obtain the lower bound of the decay estimate of the solution. Utilizi...In this paper, we study the long-time behavior of solutions of the single-layer quasi-geostrophic model arising from geophysical fluid dynamics. We obtain the lower bound of the decay estimate of the solution. Utilizing the Fourier splitting method, under suitable assumptions on the initial data, for any multi-index α, we show that the solution Ψ satisfies .展开更多
The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-eleme...The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.展开更多
Air lubrication by means of a bottom cavity is a promising method for ship drag reduction. The characteristics of the bottom cavity are sensitive to the flow field around the ship hull and the effect of drag reduction...Air lubrication by means of a bottom cavity is a promising method for ship drag reduction. The characteristics of the bottom cavity are sensitive to the flow field around the ship hull and the effect of drag reduction, especially the depth of the bottom cavity. In this study, a ship model experiment of a bulk carrier is conducted in a towing tank using the method of air layer drag reduction (ALDR) with different bottom cavity depths. The shape of the air layer is observed, and the changes in resistance are measured. The model experiments produce results of approximately 20% for the total drag reduction at the ship design speed for a 25-mm cavity continuously supplied with air at Cq = 0.224 in calm water, and the air layer covers the whole cavity when the air flow rate is suitable. In a regular head wave, the air layer is easily broken and reduces the drag reduction rate in short waves, particularly when λ/Lw1 is close to one;however, it still has a good drag reduction effect in the long waves.展开更多
The strong nonlinearity of boundary layer parameterizations in atmospheric and oceanic models can cause difficulty for tangent linear models in approximating nonlinear perturbations when the time integration grows lon...The strong nonlinearity of boundary layer parameterizations in atmospheric and oceanic models can cause difficulty for tangent linear models in approximating nonlinear perturbations when the time integration grows longer. Consequently, the related 4—D variational data assimilation problems could be difficult to solve. A modified tangent linear model is built on the Mellor-Yamada turbulent closure (level 2.5) for 4-D variational data assimilation. For oceanic mixed layer model settings, the modified tangent linear model produces better finite amplitude, nonlinear perturbation than the full and simplified tangent linear models when the integration time is longer than one day. The corresponding variational data assimilation performances based on the adjoint of the modified tangent linear model are also improved compared with those adjoints of the full and simplified tangent linear models.展开更多
To explore the failure mechanism of roadway in layered soft rocks,a physical model with the physically finite elemental slab assemblage(PFESA)method was established.Infrared thermography and a video camera were employ...To explore the failure mechanism of roadway in layered soft rocks,a physical model with the physically finite elemental slab assemblage(PFESA)method was established.Infrared thermography and a video camera were employed to capture thermal responses and deformation.The model results showed that layered soft roadway suffered from large deformation.A three-dimensional distinct element code(3 DEC)model with tetrahedral blocks was built to capture the characteristics of roadway deformation,stress,and cracks.The results showed two failure patterns,layer bending fracture and layer slipping after excavation.The layer bending fracture occurred at positions where the normal direction of layers pointed to the inside of the roadway and the layer slipping occurred in the ribs.Six schemes were proposed to investigate the effects of layered soft rocks.The results showed that the deformation of ribs was obviously larger than that of the roof and floor when the roadway passed through three types of strata.When the roadway was completely in a coal seam,the change of deformation in ribs was not obvious,while the deformation in the roof and floor increased obviously.These results can provide guidance for excavation and support design of roadways in layered soft rocks.展开更多
Layered leaf area index (LAIk) is one of the major determinants for rice canopy. The objective of this study is to attain rice LAI k using morphological traits especially leaf traits that affected plant type. A theo...Layered leaf area index (LAIk) is one of the major determinants for rice canopy. The objective of this study is to attain rice LAI k using morphological traits especially leaf traits that affected plant type. A theoretical model based on rice geometrical structure was established to describe LAI k of rice with leaf length (Li), width (Wi), angle (Ai), and space (Si), and plant pole height (H) at booting and heading stages. In correlation with traditional manual measurement, the model was performed by high R2-values (0.95-0.89, n=24) for four rice hybrids (Liangyoupeijiu, Liangyou E32, Liangyou Y06, and Shanyou 63) with various plant types and four densities (3 750, 2 812, 1 875, and 1 125 plants per 100 m2) of a particular hybrid (Liangyoupeijiu). The analysis of leaf length, width, angle, and space on LAI k for two hybrids (Liangyoupeijiu and Shanyou 63) showed that leaves length and space exhibited greater effects on the change of rice LAI k . The radiation intensity showed a significantly negative exponential relation to the accumulation of LAI k , which agreed to the coefficient of light extinction (K). Our results suggest that plant type regulates radiation distribution through changing LAI k . The present model would be helpful to acquire leaf distribution and judge canopy structure of rice field by computer system after a simple and less-invasive measurement of leaf length, width, angle (by photo), and space at field with non-dilapidation of plants.展开更多
A three-layer structure model is proposed for investigating the effect of a soft elastic middle layer on the propagation behavior of Love waves in piezoelectric layered systems, with "soft" implying that the bulk-sh...A three-layer structure model is proposed for investigating the effect of a soft elastic middle layer on the propagation behavior of Love waves in piezoelectric layered systems, with "soft" implying that the bulk-shear-wave velocity of the middle layer is smaller than that of the upper sensitive layer. Dispersion equations are obtained for unelectroded and traction-free upper surfaces which, in the limit, can be reduced to those for classical Love waves. Systematic parametric studies are subsequently carried out to quantify the effects of the soft middle layer upon Love wave propagation, including its thickness, mass density, dielectric constant and elastic coefficient. It is demonstrated that whilst the thickness and elastic coefficient of the middle layer affect significantly Love wave propagation, its mass density and dielectric constant have negligible influence. On condition that both the thickness and elastic coefficient of the middle layer are vanishingly small so that it degenerates into an imperfectly bonded interface, the three-layer model is also employed to investigate the influence of imperfect interfaces on Love waves propagating in piezoelectric layer/elastic sub- strate systems. Upon comparing with the predictions ob- tained by employing the traditional shear-lag model, the present three-layer structure model is found to be more ac- curate as it avoids the unrealistic displacement discontinuity across imperfectly bonded interfaces assumed by the shearlag model, especially for long waves when the piezoelectric layer is relatively thin.展开更多
Based on the theoretical and experimental investigation of a thin silicon layer(TSL) with linear variable doping(LVD) and further research on the TSL LVD with a multiple step field plate(MSFP),a breakdown voltag...Based on the theoretical and experimental investigation of a thin silicon layer(TSL) with linear variable doping(LVD) and further research on the TSL LVD with a multiple step field plate(MSFP),a breakdown voltage(BV) model is proposed and experimentally verified in this paper.With the two-dimensional Poisson equation of the silicon on insulator(SOI) device,the lateral electric field in drift region of the thin silicon layer is assumed to be constant.For the SOI device with LVD in the thin silicon layer,the dependence of the BV on impurity concentration under the drain is investigated by an enhanced dielectric layer field(ENDIF),from which the reduced surface field(RESURF) condition is deduced.The drain in the centre of the device has a good self-isolation effect,but the problem of the high voltage interconnection(HVI) line will become serious.The two step field plates including the source field plate and gate field plate can be adopted to shield the HVI adverse effect on the device.Based on this model,the TSL LVD SOI n-channel lateral double-diffused MOSFET(nLDMOS) with MSFP is realized.The experimental breakdown voltage(BV) and specific on-resistance(R on,sp) of the TSL LVD SOI device are 694 V and 21.3 ·mm 2 with a drift region length of 60 μm,buried oxide layer of 3 μm,and silicon layer of 0.15 μm,respectively.展开更多
基金supported by the National Natural Science Foundation of China (42030810,42104115)。
文摘The similarities and differences in inherent mechanism and characteristic frequency between the onedimensional(1D)poroelastic model and the layered White model were investigated.This investigation was conducted under the assumption that the rock was homogenous and isotropic at the mesoscopic scale.For the inherent mechanism,both models resulted from quasi-static flow in a slow P-wave diffusion mode,and the differences between them originated from saturated fluids and boundary conditions.On the other hand,for the characteristic frequencies of the models,the characteristic frequency of the 1D poroelastic model was first modified because the elastic constant and formula for calculating it were misused and then compared to that of the layered White model.Both of them moved towards higher frequencies with increasing permeability and decreasing viscosity and diffusion length.The differences between them were due to the diffusion length.The diffusion length for the 1D poroelastic model was determined by the sample length,whereas that for the layered White model was determined by the length of the representative elementary volume(REV).Subsequently,a numerical example was presented to demonstrate the similarities and differences between the models.Finally,published experimental data were interpreted using the 1D poroelastic model combined with the Cole-Cole model.The prediction of the combined model was in good agreement with the experimental data,thereby validating the effectiveness of the 1D poroelastic model.Furthermore,the modified characteristic frequency in our study was much closer to the experimental data than the previous prediction,validating the effectiveness of our modification of the characteristic frequency of the 1D poroelastic model.The investigation provided insight into the internal relationship between wave-induced fluid flow(WIFF)models at macroscopic and mesoscopic scales and can aid in a better understanding of the elastic modulus dispersion and attenuation caused by the WIFF at different scales.
基金Projects(52074299,41941018)supported by the National Natural Science Foundation of ChinaProject(2023JCCXSB02)supported by the Fundamental Research Funds for the Central Universities,China。
文摘Bedding structural planes significantly influence the mechanical properties and stability of engineering rock masses.This study conducts uniaxial compression tests on layered sandstone with various bedding angles(0°,15°,30°,45°,60°,75°and 90°)to explore the impact of bedding angle on the deformational mechanical response,failure mode,and damage evolution processes of rocks.It develops a damage model based on the Logistic equation derived from the modulus’s degradation considering the combined effect of the sandstone bedding dip angle and load.This model is employed to study the damage accumulation state and its evolution within the layered rock mass.This research also introduces a piecewise constitutive model that considers the initial compaction characteristics to simulate the whole deformation process of layered sandstone under uniaxial compression.The results revealed that as the bedding angle increases from 0°to 90°,the uniaxial compressive strength and elastic modulus of layered sandstone significantly decrease,slightly increase,and then decline again.The corresponding failure modes transition from splitting tensile failure to slipping shear failure and back to splitting tensile failure.As indicated by the modulus’s degradation,the damage characteristics can be categorized into four stages:initial no damage,damage initiation,damage acceleration,and damage deceleration termination.The theoretical damage model based on the Logistic equation effectively simulates and predicts the entire damage evolution process.Moreover,the theoretical constitutive model curves closely align with the actual stress−strain curves of layered sandstone under uniaxial compression.The introduced constitutive model is concise,with fewer parameters,a straightforward parameter determination process,and a clear physical interpretation.This study offers valuable insights into the theory of layered rock mechanics and holds implications for ensuring the safety of rock engineering.
文摘To study the evolution of nanoparticles during Al wire electrical explosion,a nanoparticle formation model that considered layered motion was developed,and an experimental system was set up to carry out electrical explosion experiments using 0.1 mm and 0.2 mm Al wires.The characteristic parameters and evolution process during the formation of nanoparticles were calculated and analyzed.The results show that the maximum velocities of the innermost and outermost layers are about 1200 m·s-1and 1600 m·s-1,and the velocity of the middle layer is about 1400 m·s-1,respectively.Most of the nanoparticles are formed in the temperature range of2600 K-2500 K.The characteristic temperature for the formation of Al nanoparticles is~2520K,which is also the characteristic temperature of other parameters.The size distribution range of the formed nanoparticles is 18 to 110 nm,and most of them are around 22 nm.The variation of saturated vapor pressure determines the temperature distribution range of particle nucleation.There is a minimum critical diameter of particles(~25 nm);particles smaller than the critical diameter can grow into larger particles during surface growth.Particle motion has an effect on the surface growth and aggregation process of particles,and also on the distribution area of larger-diameter particles.The simulation results are in good agreement with the experiments.We provide a method to estimate the size and distribution of nanoparticles,which is of great significance to understand the formation process of particles during the evolution of wire electrical explosion.
文摘In numerical simulation of wave propagation,both viscoelastic materials and perfectly matched layers(PMLs)attenuate waves.The wave equations for both the viscoelastic model and the PML contain convolution operators.However,convolution operator is intractable in finite-difference time-domain(FDTD)method.A great deal of progress has been made in using time stepping instead of convolution in FDTD.To incorporate PML into viscoelastic media,more memory variables need to be introduced,which increases the code complexity and computation costs.By modifying the nonsplitting PML formulation,I propose a viscoelastic model,which can be used as a viscoelastic material and/or a PML just by adjusting the parameters.The proposed viscoelastic model is essentially equivalent to a Maxwell model.Compared with existing PML methods,the proposed method requires less memory and its implementation in existing finite-difference codes is much easier.The attenuation and phase velocity of P-and S-waves are frequency independent in the viscoelastic model if the related quality factors(Q)are greater than 10.The numerical examples show that the method is stable for materials with high absorption(Q=1),and for heterogeneous media with large contrast of acoustic impedance and large contrast of viscosity.
基金This work was supported by the Beijing Nova Program[Z211100002121136]Open Fund Project of State Key Laboratory of Lithospheric Evolution[SKL-K202103]+1 种基金Joint Funds of National Natural Science Foundation of China[U19B6003-02]the National Natural Science Foundation of China[42302149].We would like to thank Prof.Zhu Rixiang from the Institute of Geology and Geophysics,Chinese Academy of Sciences.
文摘With continuous hydrocarbon exploration extending to deeper basins,the deepest industrial oil accumulation was discovered below 8,200 m,revealing a new exploration field.Hence,the extent to which oil exploration can be extended,and the prediction of the depth limit of oil accumulation(DLOA),are issues that have attracted significant attention in petroleum geology.Since it is difficult to characterize the evolution of the physical properties of the marine carbonate reservoir with burial depth,and the deepest drilling still cannot reach the DLOA.Hence,the DLOA cannot be predicted by directly establishing the relationship between the ratio of drilling to the dry layer and the depth.In this study,by establishing the relationships between the porosity and the depth and dry layer ratio of the carbonate reservoir,the relationships between the depth and dry layer ratio were obtained collectively.The depth corresponding to a dry layer ratio of 100%is the DLOA.Based on this,a quantitative prediction model for the DLOA was finally built.The results indicate that the porosity of the carbonate reservoir,Lower Ordovician in Tazhong area of Tarim Basin,tends to decrease with burial depth,and manifests as an overall low porosity reservoir in deep layer.The critical porosity of the DLOA was 1.8%,which is the critical geological condition corresponding to a 100%dry layer ratio encountered in the reservoir.The depth of the DLOA was 9,000 m.This study provides a new method for DLOA prediction that is beneficial for a deeper understanding of oil accumulation,and is of great importance for scientific guidance on deep oil drilling.
基金supported by the Special Fund of the Institute of Geophysics, China Earthquake Administration (No. DQJB21B32)the National Key R&D Program of China (No. 2022YFF0800601)。
文摘Surface-wave inversion is a powerful tool for revealing the Earth's internal structure.However,aside from shear-wave velocity(v_(S)),other parameters can influence the inversion outcomes,yet these have not been systematically discussed.This study investigates the influence of various parameter assumptions on the results of surface-wave inversion,including the compressional and shear velocity ratio(v_(P)/v_(S)),shear-wave attenuation(Q_(S)),density(ρ),Moho interface,and sedimentary layer.We constructed synthetic models to generate dispersion data and compared the obtained results with different parameter assumptions with those of the true model.The results indicate that the v_(P)/v_(S) ratio,Q_(S),and density(ρ) have minimal effects on absolute velocity values and perturbation patterns in the inversion.Conversely,assumptions about the Moho interface and sedimentary layer significantly influenced absolute velocity values and perturbation patterns.Introducing an erroneous Mohointerface depth in the initial model of the inversion significantly affected the v_(S) model near that depth,while using a smooth initial model results in relatively minor deviations.The assumption on the sedimentary layer not only affects shallow structure results but also impacts the result at greater depths.Non-linear inversion methods outperform linear inversion methods,particularly for the assumptions of the Moho interface and sedimentary layer.Joint inversion with other data types,such as receiver functions or Rayleigh wave ellipticity,and using data from a broader period range or higher-mode surface waves,can mitigate these deviations.Furthermore,incorporating more accurate prior information can improve inversion results.
基金supported by the Laoshan Laboratory(No.LSKJ202201600)the National Key Research and Development Program of China(No.2022YFC2808304).
文摘Surface waves have a considerable effect on vertical mixing in the upper ocean.In the past two decades,the vertical mixing induced through nonbreaking surface waves has been used in ocean and climate models to improve the simulation of the upper ocean.Thus far,several nonbreaking wave-induced mixing parameterization schemes have been proposed;however,no quantitative comparison has been performed among them.In this paper,a one-dimensional ocean model was used to compare the performances of five schemes,including those of Qiao et al.(Q),Hu and Wang(HW),Huang and Qiao(HQ),Pleskachevsky et al.(P),and Ghantous and Babanin(GB).Similar to previous studies,all of these schemes can decrease the simulated sea surface temperature(SST),increase the subsurface temperature,and deepen the mixed layer,thereby alleviating the common thermal deviation problem of the ocean model for upper ocean simulation.Among these schemes,the HQ scheme exhibited the weakest wave-induced mixing effect,and the HW scheme exhibited the strongest effect;the other three schemes exhibited roughly the same effect.In particular,the Q and P schemes exhibited nearly the same effect.In the simulation based on observations from the Ocean Weather Station Papa,the HQ scheme exhibited the best performance,followed by the Q scheme.In the experiment with the HQ scheme,the root-mean-square deviation of the simulated SST from the observations was 0.43℃,and the mixed layer depth(MLD)was 2.0 m.As a contrast,the deviations of the SST and MLD reached 1.25℃ and 8.4 m,respectively,in the experiment without wave-induced mixing.
基金supported by the National Key Research and Development Program of China (Grant No.2021YFB2600800)the National Key Research and Development 451 Program of China (Grant No.2021YFC3100803)the Guangdong Innovative and Entrepreneurial Research Team Program (Grant No.2016ZT06N340).
文摘The polyurethane foam(PU)compressible layer is a viable solution to the problem of damage to the secondary lining in squeezing tunnels.Nevertheless,the mechanical behaviour of the multi-layer yielding supports has not been thoroughly investigated.To fill this gap,large-scale model tests were conducted in this study.The synergistic load-bearing mechanics were analyzed using the convergenceconfinement method.Two types of multi-layer yielding supports with different thicknesses(2.5 cm,3.75 cm and 5 cm)of PU compressible layers were investigated respectively.Digital image correlation(DIC)analysis and acoustic emission(AE)techniques were used for detecting the deformation fields and damage evolution of the multi-layer yielding supports in real-time.Results indicated that the loaddisplacement relationship of the multi-layer yielding supports could be divided into the crack initiation,crack propagation,strain-hardening,and failure stages.Compared with those of the stiff support,the toughness,deformability and ultimate load of the yielding supports were increased by an average of 225%,61%and 32%,respectively.Additionally,the PU compressible layer is positioned between two primary linings to allow the yielding support to have greater mechanical properties.The analysis of the synergistic bearing effect suggested that the thickness of PU compressible layer and its location significantly affect the mechanical properties of the yielding supports.The use of yielding supports with a compressible layer positioned between the primary and secondary linings is recommended to mitigate the effects of high geo-stress in squeezing tunnels.
文摘The dynamic responses and generated voltage in a curved sandwich beam with glass reinforced laminate(GRL)layers and a pliable core in the presence of a piezoelectric layer under low-velocity impact(LVI)are investigated.The current study aims to carry out a dynamic analysis on the sandwich beam when the impactor hits the top face sheet with an initial velocity.For the layer analysis,the high-order shear deformation theory(HSDT)and Frostig's second model for the displacement fields of the core layer are used.The classical non-adhesive elastic contact theory and Hunter's principle are used to calculate the dynamic responses in terms of time.In order to validate the analytical method,the outcomes of the current investigation are compared with those gained by the experimental tests carried out by other researchers for a rectangular composite plate subject to the LVI.Finite element(FE)simulations are conducted by means of the ABAQUS software.The effects of the parameters such as foam modulus,layer material,fiber angle,impactor mass,and its velocity on the generated voltage are reviewed.
文摘An engineering system approach of 2-D cylindrical model of transient mass balance calculations of ozone and other concerned chemicals along with fourteen photolysis, ozone-generating and ozone-depleting chemical reaction equations was developed, validated, and used for studying the ozone concentrations, distribution and peak of the layer, ozone depletion and total ozone abundance in the stratosphere. The calculated ozone concentrations and profile at both the Equator and a 60˚N location were found to follow closely with the measured data. The calculated average ozone concentration was within 1% of the measured average, and the deviation of ozone profiles was within 14%. The monthly evolution of stratospheric ozone concentrations and distribution above the Equator was studied with results discussed in details. The influences of slow air movement in both altitudinal and radial directions on ozone concentrations and profile in the stratosphere were explored and discussed. Parametric studies of the influences of gas diffusivities of ozone D<sub>O3</sub> and active atomic oxygen D<sub>O</sub> on ozone concentrations and distributions were also studied and delineated. Having both influences through physical diffusion and chemical reactions, the diffusivity (and diffusion) of atomic oxygen D<sub>O</sub> was found to be more sensitive and important than that of ozone D<sub>O3</sub> on ozone concentrations and distribution. The 2-D ozone model present in this paper for stratospheric ozone and its layer and depletion is shown to be robust, convenient, efficient, and executable for analyzing the complex ozone phenomena in the stratosphere. .
文摘This study is an extension of the previous work done with ARS-680 Environmental Chamber. Drying is a complex operation that demands much energy and time. Drying is essentially important for preservation of ginger rhizome. Drying of ginger was modeled, and then the effective diffusion coefficient and activation energy were determined. For this purpose, the experiments were done at six levels of varied temperatures: 10°C, 20°C, 30°C, 40°C, 50°C and 60°C. The values of effective diffusion coefficients obtained in this work for the variously treated ginger rhizomes closely agreed with the average effective diffusion coefficients of other notable authors who determined the drying kinetics and convective heat transfer coefficients of ginger slices.
文摘In this paper, we study the long-time behavior of solutions of the single-layer quasi-geostrophic model arising from geophysical fluid dynamics. We obtain the lower bound of the decay estimate of the solution. Utilizing the Fourier splitting method, under suitable assumptions on the initial data, for any multi-index α, we show that the solution Ψ satisfies .
基金sponsored by the National Natural Science Foundation of China Research(Grant No.41274138)the Science Foundation of China University of Petroleum(Beijing)(No.KYJJ2012-05-02)
文摘The perfectly matched layer (PML) is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.
基金supported by the Ministry of Industry and High Technology Marine Scientific Research Projects(Grant No.2011530)the High Performance Marine Technology Key Laboratory of the Ministry of Education Open Foundation(Grant No.2013033102)
文摘Air lubrication by means of a bottom cavity is a promising method for ship drag reduction. The characteristics of the bottom cavity are sensitive to the flow field around the ship hull and the effect of drag reduction, especially the depth of the bottom cavity. In this study, a ship model experiment of a bulk carrier is conducted in a towing tank using the method of air layer drag reduction (ALDR) with different bottom cavity depths. The shape of the air layer is observed, and the changes in resistance are measured. The model experiments produce results of approximately 20% for the total drag reduction at the ship design speed for a 25-mm cavity continuously supplied with air at Cq = 0.224 in calm water, and the air layer covers the whole cavity when the air flow rate is suitable. In a regular head wave, the air layer is easily broken and reduces the drag reduction rate in short waves, particularly when λ/Lw1 is close to one;however, it still has a good drag reduction effect in the long waves.
基金Acknowledgments. The authors would like to thank Prof. Z. Yuan for her numerous suggestions in the writing of this paper. This work is supported by the National Natural Science Foundation of China (Grant No.40176009), the National Key Programme for Devel
文摘The strong nonlinearity of boundary layer parameterizations in atmospheric and oceanic models can cause difficulty for tangent linear models in approximating nonlinear perturbations when the time integration grows longer. Consequently, the related 4—D variational data assimilation problems could be difficult to solve. A modified tangent linear model is built on the Mellor-Yamada turbulent closure (level 2.5) for 4-D variational data assimilation. For oceanic mixed layer model settings, the modified tangent linear model produces better finite amplitude, nonlinear perturbation than the full and simplified tangent linear models when the integration time is longer than one day. The corresponding variational data assimilation performances based on the adjoint of the modified tangent linear model are also improved compared with those adjoints of the full and simplified tangent linear models.
基金supported by the National Key Research and Development Plan of China(No.2016YFC0600901)the National Natural Science Foundation of China(No.51874311)the Natural Science Foundation of China(No.51904306)。
文摘To explore the failure mechanism of roadway in layered soft rocks,a physical model with the physically finite elemental slab assemblage(PFESA)method was established.Infrared thermography and a video camera were employed to capture thermal responses and deformation.The model results showed that layered soft roadway suffered from large deformation.A three-dimensional distinct element code(3 DEC)model with tetrahedral blocks was built to capture the characteristics of roadway deformation,stress,and cracks.The results showed two failure patterns,layer bending fracture and layer slipping after excavation.The layer bending fracture occurred at positions where the normal direction of layers pointed to the inside of the roadway and the layer slipping occurred in the ribs.Six schemes were proposed to investigate the effects of layered soft rocks.The results showed that the deformation of ribs was obviously larger than that of the roof and floor when the roadway passed through three types of strata.When the roadway was completely in a coal seam,the change of deformation in ribs was not obvious,while the deformation in the roof and floor increased obviously.These results can provide guidance for excavation and support design of roadways in layered soft rocks.
基金supported by the National Natural Science Foundation of China (NSFC,30871479)
文摘Layered leaf area index (LAIk) is one of the major determinants for rice canopy. The objective of this study is to attain rice LAI k using morphological traits especially leaf traits that affected plant type. A theoretical model based on rice geometrical structure was established to describe LAI k of rice with leaf length (Li), width (Wi), angle (Ai), and space (Si), and plant pole height (H) at booting and heading stages. In correlation with traditional manual measurement, the model was performed by high R2-values (0.95-0.89, n=24) for four rice hybrids (Liangyoupeijiu, Liangyou E32, Liangyou Y06, and Shanyou 63) with various plant types and four densities (3 750, 2 812, 1 875, and 1 125 plants per 100 m2) of a particular hybrid (Liangyoupeijiu). The analysis of leaf length, width, angle, and space on LAI k for two hybrids (Liangyoupeijiu and Shanyou 63) showed that leaves length and space exhibited greater effects on the change of rice LAI k . The radiation intensity showed a significantly negative exponential relation to the accumulation of LAI k , which agreed to the coefficient of light extinction (K). Our results suggest that plant type regulates radiation distribution through changing LAI k . The present model would be helpful to acquire leaf distribution and judge canopy structure of rice field by computer system after a simple and less-invasive measurement of leaf length, width, angle (by photo), and space at field with non-dilapidation of plants.
基金supported by the National Natural Science Foundation of China(10972171)the Program for New Century Excellent Talents in Universities(NCET-08-0429)the National 111 Project(B06024)
文摘A three-layer structure model is proposed for investigating the effect of a soft elastic middle layer on the propagation behavior of Love waves in piezoelectric layered systems, with "soft" implying that the bulk-shear-wave velocity of the middle layer is smaller than that of the upper sensitive layer. Dispersion equations are obtained for unelectroded and traction-free upper surfaces which, in the limit, can be reduced to those for classical Love waves. Systematic parametric studies are subsequently carried out to quantify the effects of the soft middle layer upon Love wave propagation, including its thickness, mass density, dielectric constant and elastic coefficient. It is demonstrated that whilst the thickness and elastic coefficient of the middle layer affect significantly Love wave propagation, its mass density and dielectric constant have negligible influence. On condition that both the thickness and elastic coefficient of the middle layer are vanishingly small so that it degenerates into an imperfectly bonded interface, the three-layer model is also employed to investigate the influence of imperfect interfaces on Love waves propagating in piezoelectric layer/elastic sub- strate systems. Upon comparing with the predictions ob- tained by employing the traditional shear-lag model, the present three-layer structure model is found to be more ac- curate as it avoids the unrealistic displacement discontinuity across imperfectly bonded interfaces assumed by the shearlag model, especially for long waves when the piezoelectric layer is relatively thin.
基金Project supported partially by the National Natural Science Foundation of China (Grant Nos. 60906038 and 61076082)
文摘Based on the theoretical and experimental investigation of a thin silicon layer(TSL) with linear variable doping(LVD) and further research on the TSL LVD with a multiple step field plate(MSFP),a breakdown voltage(BV) model is proposed and experimentally verified in this paper.With the two-dimensional Poisson equation of the silicon on insulator(SOI) device,the lateral electric field in drift region of the thin silicon layer is assumed to be constant.For the SOI device with LVD in the thin silicon layer,the dependence of the BV on impurity concentration under the drain is investigated by an enhanced dielectric layer field(ENDIF),from which the reduced surface field(RESURF) condition is deduced.The drain in the centre of the device has a good self-isolation effect,but the problem of the high voltage interconnection(HVI) line will become serious.The two step field plates including the source field plate and gate field plate can be adopted to shield the HVI adverse effect on the device.Based on this model,the TSL LVD SOI n-channel lateral double-diffused MOSFET(nLDMOS) with MSFP is realized.The experimental breakdown voltage(BV) and specific on-resistance(R on,sp) of the TSL LVD SOI device are 694 V and 21.3 ·mm 2 with a drift region length of 60 μm,buried oxide layer of 3 μm,and silicon layer of 0.15 μm,respectively.