A finite-element model of the thermosetting epoxy asphalt mixture（EAM） microstructure is developed to simulate the indirect tension test（IDT）.Image techniques are used to capture the EAM microstructure which is di...A finite-element model of the thermosetting epoxy asphalt mixture（EAM） microstructure is developed to simulate the indirect tension test（IDT）.Image techniques are used to capture the EAM microstructure which is divided into two phases：aggregates and mastic.A viscoelastic constitutive relationship,which is obtained from the results of a creep test,is used to represent the mastic phase at intermittent temperatures.Model simulation results of the stiffness modulus in IDT compare favorably with experimental data.Different loading directions and velocities are employed in order to account for their influence on the modulus and the localized stress of the microstructure model.It is pointed out that the modulus is not consistent when the loading direction changes since the heterogeneous distribution of the mixture internal structure,and the loading velocity affects the localized stress as a result of the viscoelasticity of the mastic.The study results can provide a theoretical basis for the finite-element method,which can be extended to the numerical simulations of asphalt mixture micromechanical behavior.展开更多
To predict the behavior of geogrids embedded in sand under pullout loading conditions, the two dimensional plane-stress finite element model was presented. The interactions between soil and geogrid were simulated as ...To predict the behavior of geogrids embedded in sand under pullout loading conditions, the two dimensional plane-stress finite element model was presented. The interactions between soil and geogrid were simulated as non-linear springs, and the stiffness of the springs was determined from simple tests in the specially designed pullout box. The predicted behavior of the geogrid under pullout load agrees well with the observed data including the load-displacement properties, the displacement distribution along the longitudinal direction and the mobilization of the frictional and bearing resistance. (Edited author abstract) 8 Refs.展开更多
The finite-element modeling and simulations of the intra-body communication （IBC） were investigated to provide a theoretical basis for biomedical monitoring. A finite-element model for the whole human body was devel...The finite-element modeling and simulations of the intra-body communication （IBC） were investigated to provide a theoretical basis for biomedical monitoring. A finite-element model for the whole human body was developed to simulate the IBC. The simulation of galvanic coupling IBC and electrostatic coupling IBC were implemented along with different signal transmission paths, and their attenuations were calculated. Our study showed that the position near the signal electrode had higher potential than other positions in the two types of IBC, while the potential generally decreased along the axis of the body parts. Both signal attenuations of the two types IBC increased with increasing signal transmission distance, and the electrostatic coupling IBC had comparatively higher receiving potential than the galvanic coupling IBC. The results indicated that the proposed modeling method could be used for the research of biomedical monitoring based on IBC technology.展开更多
Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and ...Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and with the nucleus removed were developed and validatedusing published experimental and clinical data. Then the models with a stainless steel nucleusprosthesis implanted and with polymer nucleus prostheses of various properties implanted were usedfor the 3D finite-element biomechanical analysis. All the above simulation and analysis were carriedout for the L4/L5 disc under a human worst--daily compression load of 2000 N. The results show thatthe polymer materials with Young's modulus of elasticity E = 0.1-100 MPa and Poisson's ratio v=0.35-0.5 are suitable to produce artificial nucleus prosthesis in view of biomechanicalconsideration.展开更多
The 1739 M8.0 Pingluo earthquake occurred around the Yinchuan Graben,bounded by the Helan Mountains to the west and the Ordos Block to the east.Seismological observations have shown that surface fault displacement rea...The 1739 M8.0 Pingluo earthquake occurred around the Yinchuan Graben,bounded by the Helan Mountains to the west and the Ordos Block to the east.Seismological observations have shown that surface fault displacement reaches about 2–3 m,mainly by dip-slip motion along the Helanshan Piedmont Fault.However,the documented seismic intensity is distributed predominantly within the basin area,exhibiting a sharp asymmetry across the Helanshan Piedmont Fault.Thus,the general pattern of earthquake faulting is still under debate.We built a three-dimensional elastodynamic finiteelement model to reappraise the fault mechanism.In the model,predictions from synthetic rupture models,based on available observations and the earthquake scaling law,were used as an input with the split-node technique,and the effect of basin sediment on elastic wave propagation was considered.The numerical results show that if an earthquake occurred on the Helanshan Piedmont Fault characterized by a high-angle(70°)normal fault,earthquake shaking,as predicted from the modeled peak ground velocity and peak ground acceleration,has difficulty fitting the observed result,even when the effect of sediment amplification is considered.To better fit the observed shaking pattern,the dip angle of the Helanshan Piedmont Fault must be less than about 35°between the depths of about 8–27 km,where the coseismic slip may reach about 6 m.This result leads us to conclude that the 1739 M8.0 great earthquake likely occurred on a listric normal fault at depth,in agreement with the geometry of the Helanshan Piedmont Fault,as recently evidenced by seismic reflection explorations.This conclusion means that in an intracontinental setting,a reduction in the fault dip angle along the subsurface could increase the width of the fault in the elastic crust,making misalignment between the surface rupture and the isoseismals and resulting in an increase in the upper bound of earthquake magnitude relative to simple high-angle faulting.展开更多
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.展开更多
Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor...Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor CSAMT forward modeling using the vector finite-element method. To verify the feasibility of the algorithm, we calculate the electric field, magnetic field, and tensor impedance of the 3D CSAMT far-zone field in layered media and compare them with theoretical solutions. In addition, a three-dimensional anomaly in half-space is also simulated, and the response characteristics of the impedance tensor and the apparent resistivity and impedance phase are analyzed. The results suggest that the vector finite-element method produces high-precision electromagnetic field and impedance tensor data, satisfies the electric field discontinuity, and does not require divergence correction using the vector finite-element method.展开更多
The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such...The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5D forward modeling of induced polarization（IP）. In the presented method, an unstructured triangulation mesh that allows for local mesh refinement and flexible description of arbitrary model geometries is used. Furthermore, the mesh refinement process is guided by dual error estimate weighting to bias the refinement towards elements that affect the solution at the receiver locations. After the final mesh is generated, the Jacobian matrix is used to obtain the IP response on 2D structure models. We validate the adaptive FE algorithm using a vertical contact model. The validation shows that the elements near the receivers are highly refined and the average relative error of the potentials converges to 0.4 % and 1.2 % for the IP response. This suggests that the numerical solution of the adaptive FE algorithm converges to an accurate solution with the refined mesh. Finally, the accuracy and flexibility of the adaptive FE procedure are also validated using more complex models.展开更多
Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficien...Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.展开更多
A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been deve...A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.展开更多
Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at diffe...Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at different angles, is assumed to follow a lithological interface and to experience either compression or shear. The stress associated with the destruction shows an asymmetrical pattern with different distances from the highest stress sites to the fault plane in the hanging and foot walls. As the dip angle decreases,the high-stress zone becomes wider in the hanging wall but its width changes negligibly in the foot wall.The length of the simulated fault and the deformation type affect only the magnitude of maximum stress,which remains asymmetrical relative to the fault plane. The Lh/Lfratio, where Lhand Lfare the widths of high-stress zones in the hanging and foot walls of the fault, respectively, is inversely proportional to the fault plane dip. The arithmetic mean of this ratio over different fault lengths in fractures subject to compression changes from 0.29 at a dip of 80°to 1.67 at 30°. In the case of shift displacement, ratios are increasing to 1.2 and 2.94, respectively.Usually they consider vertical fault planes and symmetry in a damage zone of faults. Following that assumption may cause errors in reconstructions of stress and fault patterns in areas of complex structural setting. According geological data, we know the structures are different and asymmetric in hanging and foot walls of fault. Thus, it is important to quantify zones of that asymmetry. The modeling results have to be taken into account in studies of natural faults, especially for practical applications in seismic risk mapping, engineering geology, hydrogeology, and tectonics.展开更多
In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can be...In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling （LWD）, in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.展开更多
The accumulated large amount of satellite magnetic data strengthens our capability of resolving the electrical conductivity of Earth’s mantle.To invert these satellite magnetic data,accurate and efficient forward mod...The accumulated large amount of satellite magnetic data strengthens our capability of resolving the electrical conductivity of Earth’s mantle.To invert these satellite magnetic data,accurate and efficient forward modeling solvers are needed.In this study,a new finite-element based forward modeling solver is developed to accurately and efficiently compute the induced electromagnetic field for a realistic 3D Earth.Firstly,the nodal-based finite element method with linear shape function on tetrahedral grid is used to assemble the final system of linear equations for the magnetic vector potential and electric scalar potential.The FGMRES solver with algebraic multigrid(AMG)preconditioner is used to quickly solve the final system of linear equations.The weighted moving least-square method is employed to accurately recover the electromagnetic field from the numerical solutions of magnetic vector and electric scalar potentials.Furthermore,a local mesh refinement technique is employed to improve the accuracy of the estimated electromagnetic field.At the end,two synthetic models are used to verify the accuracy and efficiency of our newly developed forward modeling solver.A realistic 3D Earth model is used to simulate the induced magnetic field at 450 and 200 km altitudes which are the planned flying altitudes of Macao’s geomagnetic satellites.The simulation indicates that(1)the amplitude of the mantle-induced magnetic field can reach 10–30 nT at 450 km altitude,which is 10–30%of the primary magnetic field.The induced magnetic field at 200 km altitude has larger amplitudes.These mantleinduced magnetic fields can be measured by Macao geomagnetic satellites;(2)the amplitude of the ocean-induced magnetic field can reach 5–30 nT at satellite altitudes,which needs to be carefully considered in the interpretation of satellite magnetic data.We are confident that our newly developed forward modeling solver will become a key tool for interpreting satellite magnetic data.展开更多
The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LW...The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.展开更多
Rural domestic sewage treatment is critical for environmental protection.This study defines the spatial pattern of villages from the perspective of rural sewage treatment and develops an integrated decision-making sys...Rural domestic sewage treatment is critical for environmental protection.This study defines the spatial pattern of villages from the perspective of rural sewage treatment and develops an integrated decision-making system to propose a sewage treatment mode and scheme suitable for local conditions.By considering the village spatial layout and terrain factors,a decision tree model of residential density and terrain type was constructed with accuracies of 76.47%and 96.00%,respectively.Combined with binary classification probability unit regression,an appropriate sewage treatment mode for the village was determined with 87.00%accuracy.The Analytic Hierarchy Process(AHP),combined with the Technique for Order Preference(TOPSIS)by Similarity to an Ideal Solution model,formed the basis for optimal treatment process selection under different emission standards.Verification was conducted in 542 villages across three counties of the Inner Mongolia Autonomous Region,focusing on the standard effluent effect(0.3773),low investment cost(0.3196),and high standard effluent effect(0.5115)to determine the best treatment process for the same emission standard under different needs.The annual environmental and carbon emission benefits of sewage treatment in these villages were estimated.This model matches village density,geographic feature,and social development level,and provides scientific support and a theoretical basis for rural sewage treatment decision-making.展开更多
We present a comprehensive description and benchmark evaluation of the global–regional chemical transport model called the Emission and Atmospheric Processes Integrated and Coupled Community(EPICC)model.The framework...We present a comprehensive description and benchmark evaluation of the global–regional chemical transport model called the Emission and Atmospheric Processes Integrated and Coupled Community(EPICC)model.The framework incorporates(1)grid configuration,(2)transport dynamics,(3)chemical mechanisms,(4)aerosol processes,(5)wet/dry deposition parameterizations,and(6)heterogeneous chemistry treatments associated with sulfate,nitrous acid(HONO)chemistry,and aerosol/cloud–photolysis interactions(APIs/CPIs).Openly shared with the atmospheric research community,the model facilitates integration of advanced physicochemical schemes to enhance simulation accuracy.Globally,the model demonstrates realistic representations of ozone(O_(3))and aerosol optical depth.The EPICC model generally demonstrates robust performance in simulating regional concentrations of O_(3) and PM_(2.5)(and its components)in China.It successfully captures vertical profiles of both global and regional O_(3).Notably,the model mitigates frequently reported sulfate underestimations in highly industrialized regions of China.The model accurately captures two regional severe pollution episodes observed in eastern China(January/June 2021).Sensitivity experiments highlight the critical roles of heterogeneous chemical mechanisms associated with sulfate,HONO chemistry,APIs,and CPIs in capturing PM_(2.5) and O_(3) concentrations in China.Improved sulfate mechanisms result in an increase of approximately 32.4%(2.8μg m^(−3))in simulated winter sulfate concentrations when observations exceed 10μg m^(−3).Enhanced HONO elevates winter O_(3) and PM_(2.5) by≤20 and≤10μg m^(−3),respectively.Overall,CPIs dominate over APIs in improving O_(3) and PM_(2.5) simulations across China.Locally,APIs mitigate PM_(2.5) and O_(3) discrepancies in the Sichuan Basin.Seasonal cloud–chemistry coupling explains the weaker impact of PM_(2.5) in summer.展开更多
Model evaluation using benchmark datasets is an important method to measure the capability of large language models(LLMs)in specific domains,and it is mainly used to assess the knowledge and reasoning abilities of LLM...Model evaluation using benchmark datasets is an important method to measure the capability of large language models(LLMs)in specific domains,and it is mainly used to assess the knowledge and reasoning abilities of LLMs.Therefore,in order to better assess the capability of LLMs in the agricultural domain,Agri-Eval was proposed as a benchmark for assessing the knowledge and reasoning ability of LLMs in agriculture.The assessment dataset used in Agri-Eval covered seven major disciplines in the agricultural domain:crop science,horticulture,plant protection,animal husbandry,forest science,aquaculture science,and grass science,and contained a total of 2283 questions.Among domestic general-purpose LLMs,DeepSeek R1 performed best with an accuracy rate of 75.49%.In the realm of international general-purpose LLMs,Gemini 2.0 pro exp 0205 standed out as the top performer,achieving an accuracy rate of 74.28%.As an LLMs in agriculture vertical,Shennong V2.0 outperformed all the LLMs in China,and the answer accuracy rate of agricultural knowledge exceeded that of all the existing general-purpose LLMs.The launch of Agri-Eval helped the LLM developers to comprehensively evaluate the model's capability in the field of agriculture through a variety of tasks and tests to promote the development of the LLMs in the field of agriculture.展开更多
In this paper,we establish and study a single-species logistic model with impulsive age-selective harvesting.First,we prove the ultimate boundedness of the solutions of the system.Then,we obtain conditions for the asy...In this paper,we establish and study a single-species logistic model with impulsive age-selective harvesting.First,we prove the ultimate boundedness of the solutions of the system.Then,we obtain conditions for the asymptotic stability of the trivial solution and the positive periodic solution.Finally,numerical simulations are presented to validate our results.Our results show that age-selective harvesting is more conducive to sustainable population survival than non-age-selective harvesting.展开更多
The proliferation of high-dimensional data and the widespread use of complex models present central challenges in contemporary statistics and data science.Dimension reduction and model checking,as two foundational pil...The proliferation of high-dimensional data and the widespread use of complex models present central challenges in contemporary statistics and data science.Dimension reduction and model checking,as two foundational pillars supporting scientific inference and data-driven decisionmaking,have evolved through the collective wisdom of generations of statisticians.This special issue,titled"Recent Developments in Dimension Reduction and Model Checking for regressions",not only aims to showcase cutting-edge advances in the field but also carries a distinct sense of academic homage to honor the groundbreaking and enduring contributions of Professor Lixing Zhu,a leading scholar whose work has profoundly shaped both areas.展开更多
In their recent paper Pereira et al.(2025)claim that validation is overlooked in mapping and modelling of ecosystem services(ES).They state that“many studies lack critical evaluation of the results and no validation ...In their recent paper Pereira et al.(2025)claim that validation is overlooked in mapping and modelling of ecosystem services(ES).They state that“many studies lack critical evaluation of the results and no validation is provided”and that“the validation step is largely overlooked”.This assertion may have been true several years ago,for example,when Ochoa and Urbina-Cardona(2017)made a similar observation.However,there has been much work on ES model validation over the last decade.展开更多
基金Program for New Century Excellent Talents in University(No. NCET-08-0118)Specialized Research Fund for the Doctoral Program of Higher Education (No. 20090092110049)
文摘A finite-element model of the thermosetting epoxy asphalt mixture（EAM） microstructure is developed to simulate the indirect tension test（IDT）.Image techniques are used to capture the EAM microstructure which is divided into two phases：aggregates and mastic.A viscoelastic constitutive relationship,which is obtained from the results of a creep test,is used to represent the mastic phase at intermittent temperatures.Model simulation results of the stiffness modulus in IDT compare favorably with experimental data.Different loading directions and velocities are employed in order to account for their influence on the modulus and the localized stress of the microstructure model.It is pointed out that the modulus is not consistent when the loading direction changes since the heterogeneous distribution of the mixture internal structure,and the loading velocity affects the localized stress as a result of the viscoelasticity of the mastic.The study results can provide a theoretical basis for the finite-element method,which can be extended to the numerical simulations of asphalt mixture micromechanical behavior.
文摘To predict the behavior of geogrids embedded in sand under pullout loading conditions, the two dimensional plane-stress finite element model was presented. The interactions between soil and geogrid were simulated as non-linear springs, and the stiffness of the springs was determined from simple tests in the specially designed pullout box. The predicted behavior of the geogrid under pullout load agrees well with the observed data including the load-displacement properties, the displacement distribution along the longitudinal direction and the mobilization of the frictional and bearing resistance. (Edited author abstract) 8 Refs.
基金Supported by the National Natural Science Foundation of China(60801050)the Excellent Talent Fund of Beijing(2011)Excellent Young Scholars Research Fund of Beijing Institute ofTechnology(2012)
文摘The finite-element modeling and simulations of the intra-body communication （IBC） were investigated to provide a theoretical basis for biomedical monitoring. A finite-element model for the whole human body was developed to simulate the IBC. The simulation of galvanic coupling IBC and electrostatic coupling IBC were implemented along with different signal transmission paths, and their attenuations were calculated. Our study showed that the position near the signal electrode had higher potential than other positions in the two types of IBC, while the potential generally decreased along the axis of the body parts. Both signal attenuations of the two types IBC increased with increasing signal transmission distance, and the electrostatic coupling IBC had comparatively higher receiving potential than the galvanic coupling IBC. The results indicated that the proposed modeling method could be used for the research of biomedical monitoring based on IBC technology.
文摘Pilot biomechanical design of biomaterials for artificial nucleus prosthesiswas carried out based on the 3D finite-element method. Two 3D models of lumbar intervertebral discrespectively with a real human nucleus and with the nucleus removed were developed and validatedusing published experimental and clinical data. Then the models with a stainless steel nucleusprosthesis implanted and with polymer nucleus prostheses of various properties implanted were usedfor the 3D finite-element biomechanical analysis. All the above simulation and analysis were carriedout for the L4/L5 disc under a human worst--daily compression load of 2000 N. The results show thatthe polymer materials with Young's modulus of elasticity E = 0.1-100 MPa and Poisson's ratio v=0.35-0.5 are suitable to produce artificial nucleus prosthesis in view of biomechanicalconsideration.
基金Natural Science Foundation of China(No.42120104004)。
文摘The 1739 M8.0 Pingluo earthquake occurred around the Yinchuan Graben,bounded by the Helan Mountains to the west and the Ordos Block to the east.Seismological observations have shown that surface fault displacement reaches about 2–3 m,mainly by dip-slip motion along the Helanshan Piedmont Fault.However,the documented seismic intensity is distributed predominantly within the basin area,exhibiting a sharp asymmetry across the Helanshan Piedmont Fault.Thus,the general pattern of earthquake faulting is still under debate.We built a three-dimensional elastodynamic finiteelement model to reappraise the fault mechanism.In the model,predictions from synthetic rupture models,based on available observations and the earthquake scaling law,were used as an input with the split-node technique,and the effect of basin sediment on elastic wave propagation was considered.The numerical results show that if an earthquake occurred on the Helanshan Piedmont Fault characterized by a high-angle(70°)normal fault,earthquake shaking,as predicted from the modeled peak ground velocity and peak ground acceleration,has difficulty fitting the observed result,even when the effect of sediment amplification is considered.To better fit the observed shaking pattern,the dip angle of the Helanshan Piedmont Fault must be less than about 35°between the depths of about 8–27 km,where the coseismic slip may reach about 6 m.This result leads us to conclude that the 1739 M8.0 great earthquake likely occurred on a listric normal fault at depth,in agreement with the geometry of the Helanshan Piedmont Fault,as recently evidenced by seismic reflection explorations.This conclusion means that in an intracontinental setting,a reduction in the fault dip angle along the subsurface could increase the width of the fault in the elastic crust,making misalignment between the surface rupture and the isoseismals and resulting in an increase in the upper bound of earthquake magnitude relative to simple high-angle faulting.
基金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 National Natural Science Foundation of China(No.41104068)the Deep Exploration in China,Sino Probe-03-05
文摘Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor CSAMT forward modeling using the vector finite-element method. To verify the feasibility of the algorithm, we calculate the electric field, magnetic field, and tensor impedance of the 3D CSAMT far-zone field in layered media and compare them with theoretical solutions. In addition, a three-dimensional anomaly in half-space is also simulated, and the response characteristics of the impedance tensor and the apparent resistivity and impedance phase are analyzed. The results suggest that the vector finite-element method produces high-precision electromagnetic field and impedance tensor data, satisfies the electric field discontinuity, and does not require divergence correction using the vector finite-element method.
基金financially supported by the National Natural Science Foundation of China(No.41204055,41164003,and 41104074)Opening Project(No.SMIL-2014-06) of Hubei Subsurface Multi-scale Imaging Lab(SMIL),China University of Geosciences(Wuhan)
文摘The conventional finite-element（FE） method often uses a structured mesh, which is designed according to the user’s experience, and it is not sufficiently accurate and flexible to accommodate complex structures such as dipping interfaces and rough topography. We present an adaptive FE method for 2.5D forward modeling of induced polarization（IP）. In the presented method, an unstructured triangulation mesh that allows for local mesh refinement and flexible description of arbitrary model geometries is used. Furthermore, the mesh refinement process is guided by dual error estimate weighting to bias the refinement towards elements that affect the solution at the receiver locations. After the final mesh is generated, the Jacobian matrix is used to obtain the IP response on 2D structure models. We validate the adaptive FE algorithm using a vertical contact model. The validation shows that the elements near the receivers are highly refined and the average relative error of the potentials converges to 0.4 % and 1.2 % for the IP response. This suggests that the numerical solution of the adaptive FE algorithm converges to an accurate solution with the refined mesh. Finally, the accuracy and flexibility of the adaptive FE procedure are also validated using more complex models.
基金National High Technology Research and Development Program(863 Program)(No.2006AA06Z105,2007AA06Z134)
文摘Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.
基金Projects(41674080,41674079)supported by the National Natural Science Foundation of China
文摘A modeling tool for simulating three-dimensional land frequency-domain controlled-source electromagnetic surveys,based on a finite-element discretization of the Helmholtz equation for the electric fields,has been developed.The main difference between our modeling method and those previous works is edge finite-element approach applied to solving the three-dimensional land frequency-domain electromagnetic responses generated by horizontal electric dipole source.Firstly,the edge finite-element equation is formulated through the Galerkin method based on Helmholtz equation of the electric fields.Secondly,in order to check the validity of the modeling code,the numerical results are compared with the analytical solutions for a homogeneous half-space model.Finally,other three models are simulated with three-dimensional electromagnetic responses.The results indicate that the method can be applied for solving three-dimensional electromagnetic responses.The algorithm has been demonstrated,which can be effective to modeling the complex geo-electrical structures.This efficient algorithm will help to study the distribution laws of3-D land frequency-domain controlled-source electromagnetic responses and to setup basis for research of three-dimensional inversion.
文摘Stresses in a block around a dipping fracture simulating a damage zone of a fault are reconstructed by finite-element modeling. A fracture corresponding to a fault of different lengths, with its plane dipping at different angles, is assumed to follow a lithological interface and to experience either compression or shear. The stress associated with the destruction shows an asymmetrical pattern with different distances from the highest stress sites to the fault plane in the hanging and foot walls. As the dip angle decreases,the high-stress zone becomes wider in the hanging wall but its width changes negligibly in the foot wall.The length of the simulated fault and the deformation type affect only the magnitude of maximum stress,which remains asymmetrical relative to the fault plane. The Lh/Lfratio, where Lhand Lfare the widths of high-stress zones in the hanging and foot walls of the fault, respectively, is inversely proportional to the fault plane dip. The arithmetic mean of this ratio over different fault lengths in fractures subject to compression changes from 0.29 at a dip of 80°to 1.67 at 30°. In the case of shift displacement, ratios are increasing to 1.2 and 2.94, respectively.Usually they consider vertical fault planes and symmetry in a damage zone of faults. Following that assumption may cause errors in reconstructions of stress and fault patterns in areas of complex structural setting. According geological data, we know the structures are different and asymmetric in hanging and foot walls of fault. Thus, it is important to quantify zones of that asymmetry. The modeling results have to be taken into account in studies of natural faults, especially for practical applications in seismic risk mapping, engineering geology, hydrogeology, and tectonics.
基金supported by the National Natural Science Foundation of China(No.41204094)Science Foundation of China University of Petroleum,Beijing(No.2462015YQ0506)
文摘In this paper, we propose a hybrid PML （H-PML） combining the normal absorption factor of convolutional PML （C-PML） with tangential absorption factor of Mutiaxial PML （M-PML）. The H-PML boundary conditions can better suppress the numerical instability in some extreme models, and the computational speed of finite-element method and the dynamic range are greatly increased using this HPML. We use the finite-element method with a hybrid PML to model the acoustic reflection of the interface when wireline and well logging while drilling （LWD）, in a formation with a reflector outside the borehole. The simulation results suggests that the PS- and SP- reflected waves arrive at the same time when the inclination between the well and the outer interface is zero, and the difference in arrival times increases with increasing dip angle. When there are fractures outside the well, the reflection signal is clearer in the subsequent reflection waves and may be used to identify the fractured zone. The difference between the dominant wavelength and the model scale shows that LWD reflection logging data are of higher resolution and quality than wireline acoustic reflection logging.
基金supported by the National Natural Science Foundation of China(Grant Nos.72088101,41922027,41830107,41811530010)Innovation-Driven Project of Central South University(Grant No.2020CX0012)+1 种基金the National Natural Science Foundation of Hunan Province of China(Grant No.2019JJ20032)Macao Foundation and the pre-research project on Civil Aerospace Technologies funded by China’s National Space Administration(Grant Nos.D020308,D020303).
文摘The accumulated large amount of satellite magnetic data strengthens our capability of resolving the electrical conductivity of Earth’s mantle.To invert these satellite magnetic data,accurate and efficient forward modeling solvers are needed.In this study,a new finite-element based forward modeling solver is developed to accurately and efficiently compute the induced electromagnetic field for a realistic 3D Earth.Firstly,the nodal-based finite element method with linear shape function on tetrahedral grid is used to assemble the final system of linear equations for the magnetic vector potential and electric scalar potential.The FGMRES solver with algebraic multigrid(AMG)preconditioner is used to quickly solve the final system of linear equations.The weighted moving least-square method is employed to accurately recover the electromagnetic field from the numerical solutions of magnetic vector and electric scalar potentials.Furthermore,a local mesh refinement technique is employed to improve the accuracy of the estimated electromagnetic field.At the end,two synthetic models are used to verify the accuracy and efficiency of our newly developed forward modeling solver.A realistic 3D Earth model is used to simulate the induced magnetic field at 450 and 200 km altitudes which are the planned flying altitudes of Macao’s geomagnetic satellites.The simulation indicates that(1)the amplitude of the mantle-induced magnetic field can reach 10–30 nT at 450 km altitude,which is 10–30%of the primary magnetic field.The induced magnetic field at 200 km altitude has larger amplitudes.These mantleinduced magnetic fields can be measured by Macao geomagnetic satellites;(2)the amplitude of the ocean-induced magnetic field can reach 5–30 nT at satellite altitudes,which needs to be carefully considered in the interpretation of satellite magnetic data.We are confident that our newly developed forward modeling solver will become a key tool for interpreting satellite magnetic data.
基金supported by the Major Research Project on Scientific Instrument Development of the National Natural Science Foundation of China(42327901)National Natural Science Foundation of China(42030806,42074120,41904104,423B2405).
文摘The electromagnetic(EM)telemetry systems,employed for real-time data transmission from the borehole and the earth surface during drilling,are widely used in measurement-while-drilling(MWD)and logging-while-drilling(LWD).Several numerical methods,including the method of moments(MoM),the electric field integral equation(EFIE)method,and the finite-element(FE)method have been developed for the simulation of EM telemetry systems.The computational process of these methods is complicated and time-consuming.To solve this problem,we introduce an axisymmetric semi-analytical FE method(SAFEM)in the cylindrical coordinate system with the virtual layering technique for rapid simulation of EM telemetry in a layered earth.The proposed method divides the computational domain into a series of homogeneous layers.For each layer,only its cross-section is discretized,and a high-precision integration method based on Riccati equations is employed for the calculation of longitudinally homogeneous sections.The block-tridiagonal structure of the global coefficient matrix enables the use of the block Thomas algorithm,facilitating the efficient simulation of EM telemetry problems in layered media.After the theoretical development,we validate the accuracy and efficiency of our algorithm through a series of numerical experiments and comparisons with the Multiphysics modeling software COMSOL.We also discussed the impact of system parameters on EM telemetry signal and demonstrated the applicability of our method by testing it on a field dataset acquired from Dezhou,Shandong Province,China.
基金supported by the Central Government Guiding Local Science and Technology Development Fund Project(No.2024SZY0343)the Joint Research Program for Ecological Conservation and High Quality Development of the Yellow River Basin(No.2022-YRUC-01-050205)+2 种基金the Higher Education Scientific Research Project of Inner Mongolia Autonomous Region(No.NJZZ23078)the project of Inner Mongolia"Prairie Talents"Engineering Innovation Entrepreneurship Talent Team,the Major Projects of Erdos Science and Technology(No.2022EEDSKJZDZX015)the Innovation Team of the Inner Mongolia Academy of Science and Technology(No.CXTD2023-01-016).
文摘Rural domestic sewage treatment is critical for environmental protection.This study defines the spatial pattern of villages from the perspective of rural sewage treatment and develops an integrated decision-making system to propose a sewage treatment mode and scheme suitable for local conditions.By considering the village spatial layout and terrain factors,a decision tree model of residential density and terrain type was constructed with accuracies of 76.47%and 96.00%,respectively.Combined with binary classification probability unit regression,an appropriate sewage treatment mode for the village was determined with 87.00%accuracy.The Analytic Hierarchy Process(AHP),combined with the Technique for Order Preference(TOPSIS)by Similarity to an Ideal Solution model,formed the basis for optimal treatment process selection under different emission standards.Verification was conducted in 542 villages across three counties of the Inner Mongolia Autonomous Region,focusing on the standard effluent effect(0.3773),low investment cost(0.3196),and high standard effluent effect(0.5115)to determine the best treatment process for the same emission standard under different needs.The annual environmental and carbon emission benefits of sewage treatment in these villages were estimated.This model matches village density,geographic feature,and social development level,and provides scientific support and a theoretical basis for rural sewage treatment decision-making.
基金National Key Scientific and Technological Infrastructure project “Earth System Science Numerical Simulator Facility” (EarthLab)supported by the National Natural Science Foundation of China (Grant No. 92044302)the National Key Research Development Program of China (Grant No. 2022YFC3700703)
文摘We present a comprehensive description and benchmark evaluation of the global–regional chemical transport model called the Emission and Atmospheric Processes Integrated and Coupled Community(EPICC)model.The framework incorporates(1)grid configuration,(2)transport dynamics,(3)chemical mechanisms,(4)aerosol processes,(5)wet/dry deposition parameterizations,and(6)heterogeneous chemistry treatments associated with sulfate,nitrous acid(HONO)chemistry,and aerosol/cloud–photolysis interactions(APIs/CPIs).Openly shared with the atmospheric research community,the model facilitates integration of advanced physicochemical schemes to enhance simulation accuracy.Globally,the model demonstrates realistic representations of ozone(O_(3))and aerosol optical depth.The EPICC model generally demonstrates robust performance in simulating regional concentrations of O_(3) and PM_(2.5)(and its components)in China.It successfully captures vertical profiles of both global and regional O_(3).Notably,the model mitigates frequently reported sulfate underestimations in highly industrialized regions of China.The model accurately captures two regional severe pollution episodes observed in eastern China(January/June 2021).Sensitivity experiments highlight the critical roles of heterogeneous chemical mechanisms associated with sulfate,HONO chemistry,APIs,and CPIs in capturing PM_(2.5) and O_(3) concentrations in China.Improved sulfate mechanisms result in an increase of approximately 32.4%(2.8μg m^(−3))in simulated winter sulfate concentrations when observations exceed 10μg m^(−3).Enhanced HONO elevates winter O_(3) and PM_(2.5) by≤20 and≤10μg m^(−3),respectively.Overall,CPIs dominate over APIs in improving O_(3) and PM_(2.5) simulations across China.Locally,APIs mitigate PM_(2.5) and O_(3) discrepancies in the Sichuan Basin.Seasonal cloud–chemistry coupling explains the weaker impact of PM_(2.5) in summer.
文摘Model evaluation using benchmark datasets is an important method to measure the capability of large language models(LLMs)in specific domains,and it is mainly used to assess the knowledge and reasoning abilities of LLMs.Therefore,in order to better assess the capability of LLMs in the agricultural domain,Agri-Eval was proposed as a benchmark for assessing the knowledge and reasoning ability of LLMs in agriculture.The assessment dataset used in Agri-Eval covered seven major disciplines in the agricultural domain:crop science,horticulture,plant protection,animal husbandry,forest science,aquaculture science,and grass science,and contained a total of 2283 questions.Among domestic general-purpose LLMs,DeepSeek R1 performed best with an accuracy rate of 75.49%.In the realm of international general-purpose LLMs,Gemini 2.0 pro exp 0205 standed out as the top performer,achieving an accuracy rate of 74.28%.As an LLMs in agriculture vertical,Shennong V2.0 outperformed all the LLMs in China,and the answer accuracy rate of agricultural knowledge exceeded that of all the existing general-purpose LLMs.The launch of Agri-Eval helped the LLM developers to comprehensively evaluate the model's capability in the field of agriculture through a variety of tasks and tests to promote the development of the LLMs in the field of agriculture.
基金Supported by the National Natural Science Foundation of China(12261018)Universities Key Laboratory of Mathematical Modeling and Data Mining in Guizhou Province(2023013)。
文摘In this paper,we establish and study a single-species logistic model with impulsive age-selective harvesting.First,we prove the ultimate boundedness of the solutions of the system.Then,we obtain conditions for the asymptotic stability of the trivial solution and the positive periodic solution.Finally,numerical simulations are presented to validate our results.Our results show that age-selective harvesting is more conducive to sustainable population survival than non-age-selective harvesting.
文摘The proliferation of high-dimensional data and the widespread use of complex models present central challenges in contemporary statistics and data science.Dimension reduction and model checking,as two foundational pillars supporting scientific inference and data-driven decisionmaking,have evolved through the collective wisdom of generations of statisticians.This special issue,titled"Recent Developments in Dimension Reduction and Model Checking for regressions",not only aims to showcase cutting-edge advances in the field but also carries a distinct sense of academic homage to honor the groundbreaking and enduring contributions of Professor Lixing Zhu,a leading scholar whose work has profoundly shaped both areas.
文摘In their recent paper Pereira et al.(2025)claim that validation is overlooked in mapping and modelling of ecosystem services(ES).They state that“many studies lack critical evaluation of the results and no validation is provided”and that“the validation step is largely overlooked”.This assertion may have been true several years ago,for example,when Ochoa and Urbina-Cardona(2017)made a similar observation.However,there has been much work on ES model validation over the last decade.
