A whole process damage constitutive model for layered sandstone under uniaxial compression based on Logistic function

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.

A nanoparticle formation model considering layered motion based on an electrical explosion experiment with AI wires

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.

Finite-difference modeling of Maxwell viscoelastic media developed from perfectly matched layer

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.

Quantitative prediction model for the depth limit of oil accumulation in the deep carbonate rocks:A case study of Lower Ordovician in Tazhong area of Tarim Basin

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.

Assessing the effects of model parameter assumptions on surface-wave inversion results

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.

Evaluation of Nonbreaking Wave-Induced Mixing Parameterization Schemes Based on a One-Dimensional Ocean Model

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.

Dynamic and electrical responses of a curved sandwich beam with glass reinforced laminate layers and a pliable core in the presence of a piezoelectric layer under low-velocity impact

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.

2-D Modeling and Calculations of Stratospheric Ozone and Influences of Convection, Diffusion, and Time

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 DO3 and active atomic oxygen DO 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 DO was found to be more sensitive and important than that of ozone DO3 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. .

Determination of Natural Logarithm of Diffusion Coefficient and Activation Energy of Thin Layer Drying Process of Ginger Rhizome Slices

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.

Lower Bounds of Decay Rates for Solution to the Single-Layer Quasi-Geostrophic Model

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 .

Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

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.

LayeredModel:一个面向室内空间的移动对象数据模型

基于室内空间的移动对象管理或称室内移动对象管理,是一个崭新而又富有挑战的研究领域.如何建立语义完备并且支持多种应用的室内移动对象模型是一个重要的基础性问题.提出了一个分层的室内移动对象数据模型:LayeredModel.该模型提出了室内空间距离的概念,并使用不同的层次关系表达室内元素、传感器以及移动对象之间的关系.从而能支持多种基于室内空间的应用,例如移动对象的跟踪、监控、导航以及室内最近邻查询等.LayeredModel模型也为后续的室内移动对象的索引、查询的研究奠定了基础.

Experimental Study of Air Layer Drag Reduction with Bottom Cavity for A Bulk Carrier Ship Model

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 Improvement Made by a Modified TLM in 4DVAR with a Geophysical Boundary Layer Model

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.

Physical model test and numerical simulation on the failure mechanism of the roadway in layered soft rocks

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.

Performance and Analysis of a Model for Describing Layered Leaf Area Index of Rice

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 for the effect of a soft middle layer on Love waves propagating in layered piezoelectric systems

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.

A Simplified Scheme of the Generalized Layered Radiative Transfer Model

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.

Breakdown voltage model and structure realization of a thin silicon layer with linear variable doping on a silicon on insulator high voltage device with multiple step field plates

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.

ALGEBRAIC TURBUL ENCE MODEL WITH MEM ORY FOR COMPUTATION OF 3-D TURBULENT BOUNDARY LAYERS WITH VALIDATION

Additional equations were found based on experiments for an algebraic turbulence model to improve the prediction of the behavior of three dimensional turbulent boundary layers by taking account of the effects of pressure gradient and the historical variation of eddy viscosity, so the model is with memory. Numerical calculation by solving boundary layer equations was carried out for the five pressure driven three dimensional turbulent boundary layers developed on flat plates, swept wing, and prolate spheroid in symmetrical plane. Comparing the computational results with the experimental data, it is obvious that the prediction will be more accurate if the proposed closure equations are used, especially for the turbulent shear stresses.