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.展开更多
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.展开更多
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.展开更多
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.展开更多
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 press...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.展开更多
For the purpose of investigating the influence of metastable olivine(MO) phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma...For the purpose of investigating the influence of metastable olivine(MO) phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50°dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers,which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800℃and 1200℃isotherm,respectively.The elevation of olivine transformation near~410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive,zonal,negative density anomalies,respectively.These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab:(i) slab pull force exerts above a depth of~230 km;(ii) MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km;(iii) as the subducting materials accumulated over 660 km,deepest earthquakes occur due to MO transformation;(iv) a flattened‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.展开更多
Based on the geographical circumstance, climate and the boundary layer meteorology features of the Pearl River Delta, a boundary layer concept model of the Pearl River Delta was built. The concept model consists of fo...Based on the geographical circumstance, climate and the boundary layer meteorology features of the Pearl River Delta, a boundary layer concept model of the Pearl River Delta was built. The concept model consists of four fundamental factors that affect the boundary layer meteorology of the Pearl River Delta and can convincingly explain the reason of the air quality change in the Pearl River Delta. The model can be used to the diffusion capability analysis, the air pollution potential forecasting or haze forecasting, etc.展开更多
A new mesoscale air-sea coupled model (WRF- OMLM-Noh) was constructed based on the Weather Research and Forecasting (WRF) model and an improved Mellor-Yamada ocean mixed-layer model from Noh and Kim (OMLM-Noh). Throug...A new mesoscale air-sea coupled model (WRF- OMLM-Noh) was constructed based on the Weather Research and Forecasting (WRF) model and an improved Mellor-Yamada ocean mixed-layer model from Noh and Kim (OMLM-Noh). Through off-line tests and a simulation of a real typhoon, the authors compared the performance of the WRF-OMLM-Noh with another existing ocean mixed-layer coupled model (WRF-OMLM-Pollard). In the off-line tests with Tropical Ocean Global Atmosphere Program's Coupled Ocean Atmosphere Response Experiment (TOGA-COARE) observational data, the results show that OMLM-Noh is better able to simulate sea surface temperature (SST) variational trends than OMLM -Pollard. Moreover, OMLM-Noh can sufficiently reproduce the diurnal cycle of SST. Regarding the typhoon case study, SST cooling due to wind-driven ocean mixing is underestimated in WRF-OMLM-Pollard, which artificially increases the intensity of the typhoon due to more simulated air-sea heat fluxes. Compared to the WRF- OMLM-Pollard, the performance of WRF-OMLM-Noh is superior in terms of both the spatial distribution and temporal variation of SST and air-sea heat fluxes.展开更多
The dynamic ground subsidence due to underground mining is a complicated time-dependent and rate- dependent process. Based. on the theory of rock rheology and probability integral method, this study developed the supe...The dynamic ground subsidence due to underground mining is a complicated time-dependent and rate- dependent process. Based. on the theory of rock rheology and probability integral method, this study developed the superposltlOn model for the prediction and analysis of the ground dynamic subsidence in mining area of thick !oose layer. The model consists of two parts (the prediction of overlying bedrock and the prediction of thick loose layer). The overlying bedrock is regarded as visco-elastic beam, of which the dynamic subsidence is predicted by the Kelvin visco-elastic rheological model. The thick loose layer is regarded as random medium, and the ground dynamic subsidence, is predicted by the probability integral model. At last, the two prediction models are vertically stacked in the same coordinate system, and the bedrock dynamic subsidence is regarded as a variable mining thickness input into the prediction model of ground dynamic subsidence. The prediction results obtained were compared w^th actual movement and deformation data from Zhao I and Zhao II mine, central China. The agreement of the prediction results with the. field measurements.show that the superposition model (SM) is more satisfactory and the formulae obtained are more effective than the classical single probability Integral model(SPIM), and thus can be effectively used for predicting the ground dynamic subsidence in mining area of thick loose layer.展开更多
The dynamic multichannel binocular visual image modeling is studied based on Internet of Things (IoT) Perception Layer, using mobile robot self-organizing network. By employing multigroup mobile robots with binocular ...The dynamic multichannel binocular visual image modeling is studied based on Internet of Things (IoT) Perception Layer, using mobile robot self-organizing network. By employing multigroup mobile robots with binocular visual system, the real visual images of the object will be obtained. Then through the mobile self-organizing network, a three-dimensional model is rebuilt by synthesizing the returned images. On this basis, we formalize a novel algorithm for multichannel binocular visual three-dimensional images based on fast three-dimensional modeling. Compared with the method based on single binocular visual system, the new algorithm can improve the Integrity and accuracy of the dynamic three-dimensional object modeling. The simulation results show that the new method can effectively accelerate the modeling speed, improve the similarity and not increase the data size.展开更多
Wind field simulation in the surface layer is often used to manage natural resources in terms of air quality,gene flow(through pollen drift),and plant disease transmission(spore dispersion).Although Lagrangian stochas...Wind field simulation in the surface layer is often used to manage natural resources in terms of air quality,gene flow(through pollen drift),and plant disease transmission(spore dispersion).Although Lagrangian stochastic(LS)models describe stochastic wind behaviors,such models assume that wind velocities follow Gaussian distributions.However,measured surface-layer wind velocities show a strong skewness and kurtosis.This paper presents an improved model,a non-Gaussian LS model,which incorporates controllable non-Gaussian random variables to simulate the targeted non-Gaussian velocity distribution with more accurate skewness and kurtosis.Wind velocity statistics generated by the non-Gaussian model are evaluated by using the field data from the Cooperative Atmospheric Surface Exchange Study,October 1999 experimental dataset and comparing the data with statistics from the original Gaussian model.Results show that the non-Gaussian model improves the wind trajectory simulation by stably producing precise skewness and kurtosis in simulated wind velocities without sacrificing other features of the traditional Gaussian LS model,such as the accuracy in the mean and variance of simulated velocities.This improvement also leads to better accuracy in friction velocity(i.e.,a coupling of three-dimensional velocities).The model can also accommodate various non-Gaussian wind fields and a wide range of skewness–kurtosis combinations.Moreover,improved skewness and kurtosis in the simulated velocity will result in a significantly different dispersion for wind/particle simulations.Thus,the non-Gaussian model is worth applying to wind field simulation in the surface layer.展开更多
By deriving the discrete equation of the parameterized equation for the New Medium-Range Forecast(NMRF)boundary layer scheme in the GRAPES model,the adjusted discrete equation for temperature is obviously different fr...By deriving the discrete equation of the parameterized equation for the New Medium-Range Forecast(NMRF)boundary layer scheme in the GRAPES model,the adjusted discrete equation for temperature is obviously different from the original equation under the background of hydrostatic equilibrium and adiabatic hypothesis.In the present research,three discrete equations for temperature in the NMRF boundary layer scheme are applied,namely the original(hereafter NMRF),the adjustment(hereafter NMRF-gocp),and the one in the YSU boundary-layer scheme(hereafter NMRF-TZ).The results show that the deviations of height,temperature,U and V wind in the boundary layer in the NMRF-gocp and NMRF-TZ experiments are smaller than those in the NMRF experiment and the deviations in the NMRF-gocp experiment are the smallest.The deviations of humidity are complex for the different forecasting lead time in the three experiments.Moreover,there are obvious diurnal variations of deviations from these variables,where the diurnal variations of deviations from height and temperature are similar and those from U and V wind are also similar.However,the diurnal variation of humidity is relatively complicated.The root means square errors of 2m temperature(T2m)and 10m speed(V10m)from the three experiments show that the error of NMRF-gocp is the smallest and that of NMRF is the biggest.There is also a diurnal variation of T2m and V10m,where T2m has double peaks and V10m has only one peak.Comparison of the discrete equations between NMRF and NMRF-gocp experiments shows that the deviation of temperature is likely to be caused by the calculation of vertical eddy diffusive coefficients of heating,which also leads to the deviations of other elements.展开更多
Based on a modified-Darcy-Maxwell model, two-dimensional, incompressible and heat transfer flow of two bounded layers, through electrified Maxwell fluids in porous media is performed. The driving force for the instabi...Based on a modified-Darcy-Maxwell model, two-dimensional, incompressible and heat transfer flow of two bounded layers, through electrified Maxwell fluids in porous media is performed. The driving force for the instability under an electric field, is an electrostatic force exerted on the free charges accumulated at the dividing interface. Normal mode analysis is considered to study the linear stability of the disturbances layers. The solutions of the linearized equations of motion with the boundary conditions lead to an implicit dispersion relation between the growth rate and wave number. These equations are parameterized by Weber number, Reynolds number, Marangoni number, dimensionless conductivities, and dimensionless electric potentials. The case of long waves interfaciaJ stability has been studied. The stability criteria are performed theoreticaily in which stability diagrams are obtained. In the limiting cases, some previously published results can be considered as particular cases of our results. It is found that the Reynolds number plays a destabilizing role in the stability criteria, while the damping influence is observed for the increasing of Marangoni number and Maxwell relaxation time.展开更多
基金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.
基金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.
基金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.
基金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.
基金National Natural Science F oundation of China !( No.91880 10 )National Defense Science Foundation!( 95 J13 A .1.2 )
文摘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.
基金supported by the National Natural Science Foundation of China(Nos.40574047 and 40628004)
文摘For the purpose of investigating the influence of metastable olivine(MO) phase transformations on both deep seismicity and stagnation of slabs,we constructed a 2-dimensional finite element thermal model for a 120 Ma-old 50°dipping oceanic lithosphere descending at 10 cm/yr with velocity boundary layers,which would mitigate the interference of constant velocity field for the slab. The resulting temperatures show that most of intermediate and deep earthquakes occurring within the Tonga slab are occurring inside the 800℃and 1200℃isotherm,respectively.The elevation of olivine transformation near~410 km and respective persistence of metastable olivine and spinel within the transition zone and beneath 660 km would thus result in bimodal positive,zonal,negative density anomalies,respectively.These results together with the resulting pressure anomalies may reflect the stress pattern of the Tonga slab:(i) slab pull force exerts above a depth of~230 km;(ii) MO existence changes the buoyancy force within the transition zone and facilitates slab stagnation at a depth of 660 km;(iii) as the subducting materials accumulated over 660 km,deepest earthquakes occur due to MO transformation;(iv) a flattened‘slab’ may penetrate into the lower mantle due to the density increment of Sp transformation.
基金National Natural Science Foundation of China(40645026)National Basic Research Programof China (2002CB410801)
文摘Based on the geographical circumstance, climate and the boundary layer meteorology features of the Pearl River Delta, a boundary layer concept model of the Pearl River Delta was built. The concept model consists of four fundamental factors that affect the boundary layer meteorology of the Pearl River Delta and can convincingly explain the reason of the air quality change in the Pearl River Delta. The model can be used to the diffusion capability analysis, the air pollution potential forecasting or haze forecasting, etc.
基金supported by the "Strategic Priority Research Program-Climate Change: Carbon Budget andRelated Issue" of the Chinese Academy of Sciences (Grant No.XDA-05110303)the National Basic Research Program of China(Grant Nos. 2010CB951703 and 2009CB421403)the Knowledge Innovation Program of the Chinese Academy of Sciences(Grant Nos. KZCX2-YW-Q11-01 and KZCX2-YW-BR-14)
文摘A new mesoscale air-sea coupled model (WRF- OMLM-Noh) was constructed based on the Weather Research and Forecasting (WRF) model and an improved Mellor-Yamada ocean mixed-layer model from Noh and Kim (OMLM-Noh). Through off-line tests and a simulation of a real typhoon, the authors compared the performance of the WRF-OMLM-Noh with another existing ocean mixed-layer coupled model (WRF-OMLM-Pollard). In the off-line tests with Tropical Ocean Global Atmosphere Program's Coupled Ocean Atmosphere Response Experiment (TOGA-COARE) observational data, the results show that OMLM-Noh is better able to simulate sea surface temperature (SST) variational trends than OMLM -Pollard. Moreover, OMLM-Noh can sufficiently reproduce the diurnal cycle of SST. Regarding the typhoon case study, SST cooling due to wind-driven ocean mixing is underestimated in WRF-OMLM-Pollard, which artificially increases the intensity of the typhoon due to more simulated air-sea heat fluxes. Compared to the WRF- OMLM-Pollard, the performance of WRF-OMLM-Noh is superior in terms of both the spatial distribution and temporal variation of SST and air-sea heat fluxes.
基金provided by the National Natural Science Foundation of China Youth Found of China (No.41102169)the doctoral foundation of Henan Polytechnic University of China (No. B2014-056)
文摘The dynamic ground subsidence due to underground mining is a complicated time-dependent and rate- dependent process. Based. on the theory of rock rheology and probability integral method, this study developed the superposltlOn model for the prediction and analysis of the ground dynamic subsidence in mining area of thick !oose layer. The model consists of two parts (the prediction of overlying bedrock and the prediction of thick loose layer). The overlying bedrock is regarded as visco-elastic beam, of which the dynamic subsidence is predicted by the Kelvin visco-elastic rheological model. The thick loose layer is regarded as random medium, and the ground dynamic subsidence, is predicted by the probability integral model. At last, the two prediction models are vertically stacked in the same coordinate system, and the bedrock dynamic subsidence is regarded as a variable mining thickness input into the prediction model of ground dynamic subsidence. The prediction results obtained were compared w^th actual movement and deformation data from Zhao I and Zhao II mine, central China. The agreement of the prediction results with the. field measurements.show that the superposition model (SM) is more satisfactory and the formulae obtained are more effective than the classical single probability Integral model(SPIM), and thus can be effectively used for predicting the ground dynamic subsidence in mining area of thick loose layer.
基金supported by HiTech Researchand Development Program of China under Grant No.2007AA10Z235
文摘The dynamic multichannel binocular visual image modeling is studied based on Internet of Things (IoT) Perception Layer, using mobile robot self-organizing network. By employing multigroup mobile robots with binocular visual system, the real visual images of the object will be obtained. Then through the mobile self-organizing network, a three-dimensional model is rebuilt by synthesizing the returned images. On this basis, we formalize a novel algorithm for multichannel binocular visual three-dimensional images based on fast three-dimensional modeling. Compared with the method based on single binocular visual system, the new algorithm can improve the Integrity and accuracy of the dynamic three-dimensional object modeling. The simulation results show that the new method can effectively accelerate the modeling speed, improve the similarity and not increase the data size.
基金financial support for this research from a USDA-AFRI Foundational Grant (Grant No. 2012-67013-19687)from the Illinois State Water Survey at the University of Illinois at Urbana—Champaign
文摘Wind field simulation in the surface layer is often used to manage natural resources in terms of air quality,gene flow(through pollen drift),and plant disease transmission(spore dispersion).Although Lagrangian stochastic(LS)models describe stochastic wind behaviors,such models assume that wind velocities follow Gaussian distributions.However,measured surface-layer wind velocities show a strong skewness and kurtosis.This paper presents an improved model,a non-Gaussian LS model,which incorporates controllable non-Gaussian random variables to simulate the targeted non-Gaussian velocity distribution with more accurate skewness and kurtosis.Wind velocity statistics generated by the non-Gaussian model are evaluated by using the field data from the Cooperative Atmospheric Surface Exchange Study,October 1999 experimental dataset and comparing the data with statistics from the original Gaussian model.Results show that the non-Gaussian model improves the wind trajectory simulation by stably producing precise skewness and kurtosis in simulated wind velocities without sacrificing other features of the traditional Gaussian LS model,such as the accuracy in the mean and variance of simulated velocities.This improvement also leads to better accuracy in friction velocity(i.e.,a coupling of three-dimensional velocities).The model can also accommodate various non-Gaussian wind fields and a wide range of skewness–kurtosis combinations.Moreover,improved skewness and kurtosis in the simulated velocity will result in a significantly different dispersion for wind/particle simulations.Thus,the non-Gaussian model is worth applying to wind field simulation in the surface layer.
基金National Key R&D Program of China(2018YFC1506902)National Natural Science Foundation of China(42175105,U2142213)Special Fund of China Meteorological Administration for Innovation and Development(CXFZ2021Z006)。
文摘By deriving the discrete equation of the parameterized equation for the New Medium-Range Forecast(NMRF)boundary layer scheme in the GRAPES model,the adjusted discrete equation for temperature is obviously different from the original equation under the background of hydrostatic equilibrium and adiabatic hypothesis.In the present research,three discrete equations for temperature in the NMRF boundary layer scheme are applied,namely the original(hereafter NMRF),the adjustment(hereafter NMRF-gocp),and the one in the YSU boundary-layer scheme(hereafter NMRF-TZ).The results show that the deviations of height,temperature,U and V wind in the boundary layer in the NMRF-gocp and NMRF-TZ experiments are smaller than those in the NMRF experiment and the deviations in the NMRF-gocp experiment are the smallest.The deviations of humidity are complex for the different forecasting lead time in the three experiments.Moreover,there are obvious diurnal variations of deviations from these variables,where the diurnal variations of deviations from height and temperature are similar and those from U and V wind are also similar.However,the diurnal variation of humidity is relatively complicated.The root means square errors of 2m temperature(T2m)and 10m speed(V10m)from the three experiments show that the error of NMRF-gocp is the smallest and that of NMRF is the biggest.There is also a diurnal variation of T2m and V10m,where T2m has double peaks and V10m has only one peak.Comparison of the discrete equations between NMRF and NMRF-gocp experiments shows that the deviation of temperature is likely to be caused by the calculation of vertical eddy diffusive coefficients of heating,which also leads to the deviations of other elements.
文摘Based on a modified-Darcy-Maxwell model, two-dimensional, incompressible and heat transfer flow of two bounded layers, through electrified Maxwell fluids in porous media is performed. The driving force for the instability under an electric field, is an electrostatic force exerted on the free charges accumulated at the dividing interface. Normal mode analysis is considered to study the linear stability of the disturbances layers. The solutions of the linearized equations of motion with the boundary conditions lead to an implicit dispersion relation between the growth rate and wave number. These equations are parameterized by Weber number, Reynolds number, Marangoni number, dimensionless conductivities, and dimensionless electric potentials. The case of long waves interfaciaJ stability has been studied. The stability criteria are performed theoreticaily in which stability diagrams are obtained. In the limiting cases, some previously published results can be considered as particular cases of our results. It is found that the Reynolds number plays a destabilizing role in the stability criteria, while the damping influence is observed for the increasing of Marangoni number and Maxwell relaxation time.