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.展开更多
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.展开更多
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.展开更多
In this paper, firstly, a simplified version (SGRTM) of the generalized layered radiative transfer model (GRTM) within the canopy, developed by us, is presented. It reduces the information requirement of inputted ...In this paper, firstly, a simplified version (SGRTM) of the generalized layered radiative transfer model (GRTM) within the canopy, developed by us, is presented. It reduces the information requirement of inputted sky diffuse radiation, as well as of canopy morphology, and in turn saves computer resources. Results from the SGRTM agree perfectly with those of the GRTM. Secondly, by applying the linear superposition principle of the optics and by using the basic solutions of the GRTM for radiative transfer within the canopy under the condition of assumed zero soil reflectance, two sets of explicit analytical solutions of radiative transfer within the canopy with any soil reflectance magnitude are derived: one for incident diffuse, and the other for direct beam radiation. The explicit analytical solutions need two sets of basic solutions of canopy reflectance and transmittance under zero soil reflectance, run by the model for both diffuse and direct beam radiation. One set of basic solutions is the canopy reflectance αf (written as α1 for direct beam radiation) and transmittance βf (written as β1 for direction beam radiation) with zero soil reflectance for the downward radiation from above the canopy (i.e. sky), and the other set is the canopy reflectance (αb) and transmittance βb for the upward radiation from below the canopy (i.e., ground). Under the condition of the same plant architecture in the vertical layers, and the same leaf adaxial and abaxial optical properties in the canopies for the uniform diffuse radiation, the explicit solutions need only one set of basic solutions, because under this condition the two basic solutions are equal, i.e., αf = αb and βf = βb. Using the explicit analytical solutions, the fractions of any kind of incident solar radiation reflected from (defined as surface albedo, or canopy reflectance), transmitted through (defined as canopy transmittance), and absorbed by (defined as canopy absorptance) the canopy and other properties pertinent to the radiative transfer within the canopy can be estimated easily on the ground surface below the canopy (soil or snow surface) with any reflectance magnitudes. The simplified transfer model is proven to have a similar accuracy compared to the detailed model, as well as very efficient computing.展开更多
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.展开更多
Using Monte Carlo simulations, we have investigated the classical XY model on triangular lattices of ultra-thin film structures with middle ferromagnetic layers sandwiched between two antiferromagnetic layers. The int...Using Monte Carlo simulations, we have investigated the classical XY model on triangular lattices of ultra-thin film structures with middle ferromagnetic layers sandwiched between two antiferromagnetic layers. The internal energy, the specifc heat, the chirality and the chiral susceptibility are calculated in order to clarify phase transitions and critical phenomena. Prom the finite-size scaling analyses, the values of critical exponents are determined. In a range of interaction parameters, we find that the chirality steeply goes up as temperature increases in a temperature range; correspondingly the value of a critical exponent for this change is estimated.展开更多
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.展开更多
A rapid and accurate method for two-point ray tracing in horizontally layered velocity model is presented in this paper. Numerical experiments show that this method provides stable and rapid convergence with high accu...A rapid and accurate method for two-point ray tracing in horizontally layered velocity model is presented in this paper. Numerical experiments show that this method provides stable and rapid convergence with high accuracies, regardless of various 1-D velocity structures, takeoff angles and epicentral distances. This two-point ray tracing method is compared with the pseudobending technique and the method advanced by Kim and Baag (2002). It turns out that the method in this paper is much more efficient and accurate than the pseudobending technique, but is only applicable to 1-D velocity model. Kim's method is equivalent to ours for cases without large takeoff angles, but it fails to work when the takeoff angle is close to 90°. On the other hand, the method presented in this paper is applicable to cases with any takeoff angles with rapid and accurate convergence. Therefore, this method is a good choice for two-point ray tracing problems in horizontally layered velocity model and is efficient enough to be applied to a wide range of seismic problems.展开更多
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 .展开更多
In this paper, the new model of the real gas filtration problem has been presented multi-layered gas reservoir, when a gas well output and wellbore storage may be variable, and have obtained the exact solutions of pre...In this paper, the new model of the real gas filtration problem has been presented multi-layered gas reservoir, when a gas well output and wellbore storage may be variable, and have obtained the exact solutions of pressure distribution for each reservoir bed under three kinds of typical out-boundary conditions. As a special case, according to the new model have also obtained the qxact solutions of presssure distribution in homogeneous reservoir and is given important application in gas reservoir development.展开更多
基金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.
基金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.
基金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.
基金This work was supported by the National Natural Science Foundation of China under Grant Nos. 40233034, 40575043the Chinese Academy of Sciences (KZCX3_SW_229).
文摘In this paper, firstly, a simplified version (SGRTM) of the generalized layered radiative transfer model (GRTM) within the canopy, developed by us, is presented. It reduces the information requirement of inputted sky diffuse radiation, as well as of canopy morphology, and in turn saves computer resources. Results from the SGRTM agree perfectly with those of the GRTM. Secondly, by applying the linear superposition principle of the optics and by using the basic solutions of the GRTM for radiative transfer within the canopy under the condition of assumed zero soil reflectance, two sets of explicit analytical solutions of radiative transfer within the canopy with any soil reflectance magnitude are derived: one for incident diffuse, and the other for direct beam radiation. The explicit analytical solutions need two sets of basic solutions of canopy reflectance and transmittance under zero soil reflectance, run by the model for both diffuse and direct beam radiation. One set of basic solutions is the canopy reflectance αf (written as α1 for direct beam radiation) and transmittance βf (written as β1 for direction beam radiation) with zero soil reflectance for the downward radiation from above the canopy (i.e. sky), and the other set is the canopy reflectance (αb) and transmittance βb for the upward radiation from below the canopy (i.e., ground). Under the condition of the same plant architecture in the vertical layers, and the same leaf adaxial and abaxial optical properties in the canopies for the uniform diffuse radiation, the explicit solutions need only one set of basic solutions, because under this condition the two basic solutions are equal, i.e., αf = αb and βf = βb. Using the explicit analytical solutions, the fractions of any kind of incident solar radiation reflected from (defined as surface albedo, or canopy reflectance), transmitted through (defined as canopy transmittance), and absorbed by (defined as canopy absorptance) the canopy and other properties pertinent to the radiative transfer within the canopy can be estimated easily on the ground surface below the canopy (soil or snow surface) with any reflectance magnitudes. The simplified transfer model is proven to have a similar accuracy compared to the detailed model, as well as very efficient computing.
文摘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.
基金Project supported by the National Natural Science Foundation of China (Grant No 10234010), Acknowledgements The authors would like to thank Professor T. Horiguch for his useful discussion. 0ur simulations were carried out on IBM RS/6000 SP3 at the Centre for Computational Science and Engineering of Peking University.
文摘Using Monte Carlo simulations, we have investigated the classical XY model on triangular lattices of ultra-thin film structures with middle ferromagnetic layers sandwiched between two antiferromagnetic layers. The internal energy, the specifc heat, the chirality and the chiral susceptibility are calculated in order to clarify phase transitions and critical phenomena. Prom the finite-size scaling analyses, the values of critical exponents are determined. In a range of interaction parameters, we find that the chirality steeply goes up as temperature increases in a temperature range; correspondingly the value of a critical exponent for this change is estimated.
文摘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.
基金National Natural Science Foundation of China (40074008 and 40134010).
文摘A rapid and accurate method for two-point ray tracing in horizontally layered velocity model is presented in this paper. Numerical experiments show that this method provides stable and rapid convergence with high accuracies, regardless of various 1-D velocity structures, takeoff angles and epicentral distances. This two-point ray tracing method is compared with the pseudobending technique and the method advanced by Kim and Baag (2002). It turns out that the method in this paper is much more efficient and accurate than the pseudobending technique, but is only applicable to 1-D velocity model. Kim's method is equivalent to ours for cases without large takeoff angles, but it fails to work when the takeoff angle is close to 90°. On the other hand, the method presented in this paper is applicable to cases with any takeoff angles with rapid and accurate convergence. Therefore, this method is a good choice for two-point ray tracing problems in horizontally layered velocity model and is efficient enough to be applied to a wide range of seismic problems.
文摘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 .
文摘In this paper, the new model of the real gas filtration problem has been presented multi-layered gas reservoir, when a gas well output and wellbore storage may be variable, and have obtained the exact solutions of pressure distribution for each reservoir bed under three kinds of typical out-boundary conditions. As a special case, according to the new model have also obtained the qxact solutions of presssure distribution in homogeneous reservoir and is given important application in gas reservoir development.