Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference ...Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference techniques. The existence and uniqueness of the weak solution are proved without any assumptions on choice of the spacetime meshes. Basic error estimates in L-infinity (L-2) norm, that is maximum-norm in time, L-2-norm in space are obtained. The numerical results are given in the last part and the analysis between theoretic and experimental results are obtained.展开更多
A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order...A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.展开更多
Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved thro...Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.展开更多
A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FE...A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FEM is described and compared with ES-FEM.A practical modification of IES-FEM is then introduced that used the technique employed by ES-FEM for the nodal strain calculation.The differences in the strain computation among ES-FEM,IES-FEM,and FEM are then discussed.The modified IES-FEM exhibited superior performance in displacement and a slight advantage in stress compared to FEM using the same mesh according to the results obtained from both the regular and irregular elements.The robustness of the IES-FEM to severely deformed meshes was also verified.展开更多
We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to...We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to achieve exceptional computational efficiency and accuracy.The workflow is demonstrated through the modeling of wireline electromagnetic propagation resistivity logs,where the measured responses exhibit a highly nonlinear relationship with formation properties.The motivation for this research is the need for advanced modeling al-gorithms that are fast enough for use in modern quantitative interpretation tools,where thousands of simulations may be required in iterative inversion processes.The proposed algorithm achieves a remarkable enhancement in performance,being up to 3000 times faster than the finite element method alone when utilizing a GPU.While still ensuring high accuracy,this makes it well-suited for practical applications when reliable payzone assessment is needed in complex environmental scenarios.Furthermore,the algorithm’s efficiency positions it as a promising tool for stochastic Bayesian inversion,facilitating reliable uncertainty quantification in subsurface property estimation.展开更多
Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant i...Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant importance.The traditional finite element method(FEM)remains one of the primary approaches for addressing such issues.However,the application of FEM typically necessitates the use of a fine finite element mesh to accurately capture the heterogeneous properties of the materials and meet the required computational precision,which inevitably leads to a reduction in computational efficiency.To enhance the computational accuracy and efficiency of the FEM for heterogeneous multi-field coupling problems,this study presents the coupling magneto-electro-elastic multiscale finite element method(CM-MsFEM)for heterogeneous MEE structures.Unlike the conventional multiscale FEM(MsFEM),the proposed algorithm simultaneously constructs displacement,electric,and magnetic potential multiscale basis functions to address the heterogeneity of the corresponding parameters.The macroscale formulation of CM-MsFEM was derived,and the macroscale/microscale responses of the problems were obtained through up/downscaling calculations.Evaluation using numerical examples analyzing the transient behavior of heterogeneous MEE structures demonstrated that the proposed method outperforms traditional FEM in terms of both accuracy and computational efficiency,making it an appropriate choice for numerically modeling the dynamics of heterogeneous MEE structures.展开更多
The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX a...The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX are determined by the grain level deformation and sequentially coupled.In order to fully capture the microstructure and mechanical property evolution,a crystal plasticity finite element based modelling method for DRX and SRX is proposed in the current work.The grain level deformation is calculated with crystal plasticity which is coupled with the recrystallization model straightforwardly,and both the grain deformation and microstructure evolution are updated simultaneously.The proposed method is validated with discontinuous DRX experiments and the effects of initial deformation conditions are well-captured.Two controversial mechanisms for recrystallization microstructure evolution,i.e.oriented nucleation and growth selection,are discussed in the current framework with the advantages of accurate grain level deformation and interaction predictions.Furthermore,the sequentially coupled DRX and SRX are modelled seamlessly in the current work which provides a critical method for fully integrated thermo-mechanical processes analysis.展开更多
In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a ge...In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a generalized nonlinear Stokes problem of displacement vector field related to pseudo pressure and a diffusion problem of other pseudo pressure fields.Secondly,a fully discrete multiphysics finite element method is performed to solve the reformulated system numerically.Thirdly,existence and uniqueness of the weak solution of the reformulated model and stability analysis and optimal convergence order for the multiphysics finite element method are proven theoretically.Lastly,numerical tests are given to verify the theoretical results.展开更多
Skin panels on supersonic vehicles are subjected to aero-thermo-acoustic loads,resulting in a well-known multi-physics dynamic problem.The high-frequency dynamic response of these panels significantly impacts the stru...Skin panels on supersonic vehicles are subjected to aero-thermo-acoustic loads,resulting in a well-known multi-physics dynamic problem.The high-frequency dynamic response of these panels significantly impacts the structural safety of supersonic vehicles,but it has been rarely investigated.Given that existing methods are inefficient for high-frequency dynamic analysis in multi-physics fields,the present work addresses this challenge by proposing a Stochastic Energy Finite Element Method(SEFEM).SEFEM uses energy density instead of displacement to describe the dynamic response,thereby significantly enhancing its efficiency.In SEFEM,the effects of aerodynamic and thermal loads on the energy propagation characteristics are studied analytically and incorporated into the energy density governing equation.These effects are also considered when calculating the input power generated by the acoustic load,and two effective approaches named Frequency Response Function Method(FRFM)and Mechanical Impedance Method(MIM)are developed accordingly and integrated into SEFEM.The good accuracy,applicability,and high efficiency of the proposed SEFEM are demonstrated through numerical simulations performed on a two-dimensional panel under aero-thermoacoustic loads.Additionally,the effects and underlying mechanisms of aero-thermo-acoustic loads on the high-frequency response are explored.This work not only presents an efficient approach for predicting high-frequency dynamic response of panels subjected to aero-thermo-acoustic loads,but also provides insights into the high-frequency dynamic characteristics in multi-physics fields.展开更多
Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump mate...Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.展开更多
Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres wi...Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres with artificial inherent pores are an important type of lightweight proppant,enabling their transport to distant fracture extremities and enhancing fracture conductivity.However,the focus frequently gravitates towards the low-density advantage,often overlooking the pore geometry impacts on compressive strength by traditional strength evaluation.This paper numerically bypasses such limitations by using a combined finite and discrete element method(FDEM)considering experimental results.The mesh size of the model undergoes validation,followed by the calibration of cohesive element parameters via the single particle compression test.The stimulation elucidates that proppants with a smaller pore size(40μm)manifest crack propagation evolution at a more rapid pace in comparison to their larger-pore counterparts,though the influence of pore diameter on overall strength is subtle.The inception of pores not only alters the trajectory of crack progression but also,with an increase in porosity,leads to a discernible decline in proppant compressive strength.Intriguingly,upon crossing a porosity threshold of 10%,the decrement in strength becomes more gradual.A denser congregation of pores accelerates crack propagation,undermining proppant robustness,suggesting that under analogous conditions,hollow proppants might not match the strength of their porous counterparts.This exploration elucidates the underlying mechanisms of proppant failure from a microstructural perspective,furnishing pivotal insights that may guide future refinements in the architectural design of porous proppant.展开更多
This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is e...This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.展开更多
In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow e...In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.展开更多
Controlled nuclear fusion represents a significant solution for future clean energy,with ion cyclotron range of frequency(ICRF)heating emerging as one of the most promising technologies for heating the fusion plasma.T...Controlled nuclear fusion represents a significant solution for future clean energy,with ion cyclotron range of frequency(ICRF)heating emerging as one of the most promising technologies for heating the fusion plasma.This study primarily presents a self-developed 2D ion cyclotron resonance antenna electromagnetic field solver(ICRAEMS)code implemented on the MATLAB platform,which solves the electric field wave equation by using the finite element method,establishing perfectly matched layer(PML)boundary conditions,and post-processing the electromagnetic field data.This code can be utilized to facilitate the design and optimization processes of antennas for ICRF heating technology.Furthermore,this study examines the electric field distribution and power spectrum associated with various antenna phases to investigate how different antenna configurations affect the electromagnetic field propagation and coupling characteristics.展开更多
This study explores the magnetohydrodynamic(MHD)boundary layer flow of a water-based Boger nanofluid over a stretching sheet,with particular focus on the influences of nanoparticle diameter,nanolayer effects,and therm...This study explores the magnetohydrodynamic(MHD)boundary layer flow of a water-based Boger nanofluid over a stretching sheet,with particular focus on the influences of nanoparticle diameter,nanolayer effects,and thermal radiation.The primary aim is to examine how variations in nanoparticle size and nanolayer thickness affect the hydrothermal behavior of the nanofluid.The model also incorporates the contributions of viscous dissipation and Joule heating within the heat transfer equation.The governing momentum and energy equations are converted into dimensionless partial differential equations(PDEs)using appropriate similarity variables and are numerically solved using the finite element method(FEM)implemented in MATLAB.Extensive validation of this method confirms its reliability and accuracy in numerical solutions.The findings reveal that increasing the diameter of copper nanoparticles significantly enhances the velocity profile,with a more pronounced effect observed at wider inter-particle spacings.A higher solvent volume fraction leads to decreased velocity and temperature distributions,while a greater relaxation time ratio improves velocity and temperature profiles due to the increased elastic response of the fluid.Moreover,enhancements in the magnetic parameter,thermal radiation,and Eckert number lead to an elevation in temperature profiles.Furthermore,higher nanolayer thickness reduces the temperature profile,whereas particle radius yields the opposite outcome.展开更多
Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis r...Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.展开更多
Physics-informed neural networks(PINNs)have prevailed as differentiable simulators to investigate flow in porous media.Despite recent progress PINNs have achieved,practical geotechnical scenarios cannot be readily sim...Physics-informed neural networks(PINNs)have prevailed as differentiable simulators to investigate flow in porous media.Despite recent progress PINNs have achieved,practical geotechnical scenarios cannot be readily simulated because conventional PINNs fail in discontinuous heterogeneous porous media or multi-layer strata when labeled data are missing.This work aims to develop a universal network structure to encode the mass continuity equation and Darcy’s law without labeled data.The finite element approximation,which can decompose a complex heterogeneous domain into simpler ones,is adopted to build the differentiable network.Without conventional DNNs,physics-encoded finite element network(PEFEN)can avoid spectral bias and learn high-frequency functions with sharp/steep gradients.PEFEN rigorously encodes Dirichlet and Neumann boundary conditions without training.Benefiting from its discretized formulation,the discontinuous heterogeneous hydraulic conductivity is readily embedded into the network.Three typical cases are reproduced to corroborate PEFEN’s superior performance over conventional PINNs and the PINN with mixed formulation.PEFEN is sparse and demonstrated to be capable of dealing with heterogeneity with much fewer training iterations(less than 1/30)than the improved PINN with mixed formulation.Thus,PEFEN saves energy and contributes to low-carbon AI for science.The last two cases focus on common geotechnical settings of impermeable sheet pile in singlelayer and multi-layer strata.PEFEN solves these cases with high accuracy,circumventing costly labeled data,extra computational burden,and additional treatment.Thus,this study warrants the further development and application of PEFEN as a novel differentiable network in porous flow of practical geotechnical engineering.展开更多
In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hy...In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.展开更多
In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
This article presents a complete discretization of a nonlinear Sobolev equation using space-time discontinuous Galerkin method that is discontinuous in time and continuous in space. The scheme is formulated by introdu...This article presents a complete discretization of a nonlinear Sobolev equation using space-time discontinuous Galerkin method that is discontinuous in time and continuous in space. The scheme is formulated by introducing the equivalent integral equation of the primal equation. The proposed scheme does not explicitly include the jump terms in time, which represent the discontinuity characteristics of approximate solution. And then the complexity of the theoretical analysis is reduced. The existence and uniqueness of the approximate solution and the stability of the scheme are proved. The optimalorder error estimates in L2 （H1） and L2 （L2） norms are derived. These estimates are valid under weak restrictions on the space-time mesh, namely, without the condition kn ≥ ch2, which is necessary in traditional space-time discontinuous Galerkin methods. Numerical experiments are presented to verify the theoretical results.展开更多
文摘Adaptive space-time finite element method, continuous in space but discontinuous in time for semi-linear parabolic problems is discussed. The approach is based on a combination of finite element and finite difference techniques. The existence and uniqueness of the weak solution are proved without any assumptions on choice of the spacetime meshes. Basic error estimates in L-infinity (L-2) norm, that is maximum-norm in time, L-2-norm in space are obtained. The numerical results are given in the last part and the analysis between theoretic and experimental results are obtained.
基金supported by the National Natural Science Foundation of China (No. 10601022)NSF ofInner Mongolia Autonomous Region of China (No. 200607010106)513 and Science Fund of InnerMongolia University for Distinguished Young Scholars (No. ND0702)
文摘A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.
基金Project supported by the National Basic Research Program of China (973 program) (No.G1999032804)
文摘Energy conservation of nonlinear Schrodinger ordinary differential equation was proved through using continuous finite element methods of ordinary differential equation; Energy integration conservation was proved through using space-time continuous fully discrete finite element methods and the electron nearly conservation with higher order error was obtained through using time discontinuous only space continuous finite element methods of nonlinear Schrodinger partial equation. The numerical results are in accordance with the theory.
基金the National Natural Science Foundation of China(No.11672238)the 111 Project(No.BP0719007)the Shaanxi Province Natural Science Foundation(No.2020JZ-06)for the financial support.
文摘A modified inner-element edge-based smoothed finite element method(IES-FEM)is developed and integrated with ABAQUS using a user-defined element(UEL)in this study.Initially,the smoothing domain discretization of IES-FEM is described and compared with ES-FEM.A practical modification of IES-FEM is then introduced that used the technique employed by ES-FEM for the nodal strain calculation.The differences in the strain computation among ES-FEM,IES-FEM,and FEM are then discussed.The modified IES-FEM exhibited superior performance in displacement and a slight advantage in stress compared to FEM using the same mesh according to the results obtained from both the regular and irregular elements.The robustness of the IES-FEM to severely deformed meshes was also verified.
基金financially supported by the Russian federal research project No.FWZZ-2022-0026“Innovative aspects of electro-dynamics in problems of exploration and oilfield geophysics”.
文摘We propose a novel workflow for fast forward modeling of well logs in axially symmetric 2D models of the nearwellbore environment.The approach integrates the finite element method with deep residual neural networks to achieve exceptional computational efficiency and accuracy.The workflow is demonstrated through the modeling of wireline electromagnetic propagation resistivity logs,where the measured responses exhibit a highly nonlinear relationship with formation properties.The motivation for this research is the need for advanced modeling al-gorithms that are fast enough for use in modern quantitative interpretation tools,where thousands of simulations may be required in iterative inversion processes.The proposed algorithm achieves a remarkable enhancement in performance,being up to 3000 times faster than the finite element method alone when utilizing a GPU.While still ensuring high accuracy,this makes it well-suited for practical applications when reliable payzone assessment is needed in complex environmental scenarios.Furthermore,the algorithm’s efficiency positions it as a promising tool for stochastic Bayesian inversion,facilitating reliable uncertainty quantification in subsurface property estimation.
基金supported by the National Natural Science Foundation of China(Grant Nos.42102346,42172301).
文摘Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant importance.The traditional finite element method(FEM)remains one of the primary approaches for addressing such issues.However,the application of FEM typically necessitates the use of a fine finite element mesh to accurately capture the heterogeneous properties of the materials and meet the required computational precision,which inevitably leads to a reduction in computational efficiency.To enhance the computational accuracy and efficiency of the FEM for heterogeneous multi-field coupling problems,this study presents the coupling magneto-electro-elastic multiscale finite element method(CM-MsFEM)for heterogeneous MEE structures.Unlike the conventional multiscale FEM(MsFEM),the proposed algorithm simultaneously constructs displacement,electric,and magnetic potential multiscale basis functions to address the heterogeneity of the corresponding parameters.The macroscale formulation of CM-MsFEM was derived,and the macroscale/microscale responses of the problems were obtained through up/downscaling calculations.Evaluation using numerical examples analyzing the transient behavior of heterogeneous MEE structures demonstrated that the proposed method outperforms traditional FEM in terms of both accuracy and computational efficiency,making it an appropriate choice for numerically modeling the dynamics of heterogeneous MEE structures.
基金supported by the National Natural Science Foundation of China(Nos.52105384 and U2141215).
文摘The microstructure and related property evolution induced by dynamic recrystallization(DRX)and static recrystallization(SRX)in thermo-mechanical process are two critical factors for the metal forming.The DRX and SRX are determined by the grain level deformation and sequentially coupled.In order to fully capture the microstructure and mechanical property evolution,a crystal plasticity finite element based modelling method for DRX and SRX is proposed in the current work.The grain level deformation is calculated with crystal plasticity which is coupled with the recrystallization model straightforwardly,and both the grain deformation and microstructure evolution are updated simultaneously.The proposed method is validated with discontinuous DRX experiments and the effects of initial deformation conditions are well-captured.Two controversial mechanisms for recrystallization microstructure evolution,i.e.oriented nucleation and growth selection,are discussed in the current framework with the advantages of accurate grain level deformation and interaction predictions.Furthermore,the sequentially coupled DRX and SRX are modelled seamlessly in the current work which provides a critical method for fully integrated thermo-mechanical processes analysis.
基金Supported by the National Natural Science Foundation of China(Grant Nos.12371393,11971150 and 11801143)Natural Science Foundation of Henan Province(Grant No.242300421047).
文摘In this paper,we propose a multiphysics finite element method for a nonlinear poroelasticity model with nonlinear stress-strain relation.Firstly,we reformulate the original problem into a new coupled fluid system-a generalized nonlinear Stokes problem of displacement vector field related to pseudo pressure and a diffusion problem of other pseudo pressure fields.Secondly,a fully discrete multiphysics finite element method is performed to solve the reformulated system numerically.Thirdly,existence and uniqueness of the weak solution of the reformulated model and stability analysis and optimal convergence order for the multiphysics finite element method are proven theoretically.Lastly,numerical tests are given to verify the theoretical results.
基金financially supported by the National Natural Science Foundation of China(Nos.12302228 and 12372170)。
文摘Skin panels on supersonic vehicles are subjected to aero-thermo-acoustic loads,resulting in a well-known multi-physics dynamic problem.The high-frequency dynamic response of these panels significantly impacts the structural safety of supersonic vehicles,but it has been rarely investigated.Given that existing methods are inefficient for high-frequency dynamic analysis in multi-physics fields,the present work addresses this challenge by proposing a Stochastic Energy Finite Element Method(SEFEM).SEFEM uses energy density instead of displacement to describe the dynamic response,thereby significantly enhancing its efficiency.In SEFEM,the effects of aerodynamic and thermal loads on the energy propagation characteristics are studied analytically and incorporated into the energy density governing equation.These effects are also considered when calculating the input power generated by the acoustic load,and two effective approaches named Frequency Response Function Method(FRFM)and Mechanical Impedance Method(MIM)are developed accordingly and integrated into SEFEM.The good accuracy,applicability,and high efficiency of the proposed SEFEM are demonstrated through numerical simulations performed on a two-dimensional panel under aero-thermoacoustic loads.Additionally,the effects and underlying mechanisms of aero-thermo-acoustic loads on the high-frequency response are explored.This work not only presents an efficient approach for predicting high-frequency dynamic response of panels subjected to aero-thermo-acoustic loads,but also provides insights into the high-frequency dynamic characteristics in multi-physics fields.
基金the financial support provided by MHRD,Govt.of IndiaCoal India Limited for providing financial assistance for the research(Project No.CIL/R&D/01/73/2021)the partial financial support provided by the Ministry of Education,Government of India,under SPARC project(Project No.P1207)。
文摘Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.
基金the financial support provided by Tianfu Yongxing Laboratory Organized Research Project Funding(No.2023CXXM01)the ARC linkage program(No.LP200100420).
文摘Ceramic spheres,typically with a particle diameter of less than 0.8 mm,are frequently utilized as a critical proppant material in hydraulic fracturing for petroleum and natural gas extraction.Porous ceramic spheres with artificial inherent pores are an important type of lightweight proppant,enabling their transport to distant fracture extremities and enhancing fracture conductivity.However,the focus frequently gravitates towards the low-density advantage,often overlooking the pore geometry impacts on compressive strength by traditional strength evaluation.This paper numerically bypasses such limitations by using a combined finite and discrete element method(FDEM)considering experimental results.The mesh size of the model undergoes validation,followed by the calibration of cohesive element parameters via the single particle compression test.The stimulation elucidates that proppants with a smaller pore size(40μm)manifest crack propagation evolution at a more rapid pace in comparison to their larger-pore counterparts,though the influence of pore diameter on overall strength is subtle.The inception of pores not only alters the trajectory of crack progression but also,with an increase in porosity,leads to a discernible decline in proppant compressive strength.Intriguingly,upon crossing a porosity threshold of 10%,the decrement in strength becomes more gradual.A denser congregation of pores accelerates crack propagation,undermining proppant robustness,suggesting that under analogous conditions,hollow proppants might not match the strength of their porous counterparts.This exploration elucidates the underlying mechanisms of proppant failure from a microstructural perspective,furnishing pivotal insights that may guide future refinements in the architectural design of porous proppant.
基金supported by the National Natural Science Foundation of China(Grant Nos.51890912,51979025 and 52011530189).
文摘This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.
基金supported by the Natural Science Foundation of Shandong Province(ZR2021MA019)the National Natural Science Foundation of China(11871312)。
文摘In this paper,a composite numerical scheme is proposed to solve the threedimensional Darcy-Forchheimer miscible displacement problem with positive semi-definite assumptions.A mixed finite element is used for the fow equation.The velocity and pressure are computed simultaneously.The accuracy of velocity is improved one order.The concentration equation is solved by using mixed finite element,multi-step difference and upwind approximation.A multi-step method is used to approximate time derivative for improving the accuracy.The upwind approximation and an expanded mixed finite element are adopted to solve the convection and diffusion,respectively.The composite method could compute the diffusion flux and its gradient.It possibly becomes an eficient tool for solving convection-dominated diffusion problems.Firstly,the conservation of mass holds.Secondly,the multi-step method has high accuracy.Thirdly,the upwind approximation could avoid numerical dispersion.Using numerical analysis of a priori estimates and special techniques of differential equations,we give an error estimates for a positive definite problem.Numerical experiments illustrate its computational efficiency and feasibility of application.
基金Project supported by the National MCF Energy R&D Program(Grant No.2022YFE03190100)the National Natural Science Foundation of China(Grant Nos.12422513,12105035,and U21A20438)the Xiaomi Young Talents Program.
文摘Controlled nuclear fusion represents a significant solution for future clean energy,with ion cyclotron range of frequency(ICRF)heating emerging as one of the most promising technologies for heating the fusion plasma.This study primarily presents a self-developed 2D ion cyclotron resonance antenna electromagnetic field solver(ICRAEMS)code implemented on the MATLAB platform,which solves the electric field wave equation by using the finite element method,establishing perfectly matched layer(PML)boundary conditions,and post-processing the electromagnetic field data.This code can be utilized to facilitate the design and optimization processes of antennas for ICRF heating technology.Furthermore,this study examines the electric field distribution and power spectrum associated with various antenna phases to investigate how different antenna configurations affect the electromagnetic field propagation and coupling characteristics.
基金supported by the Fundamental Research Funds for the Central Universities(Grant No.D5000230061)the Guangdong Basic and Applied Basic Research Foundation(Grant No.2025A1515011192).
文摘This study explores the magnetohydrodynamic(MHD)boundary layer flow of a water-based Boger nanofluid over a stretching sheet,with particular focus on the influences of nanoparticle diameter,nanolayer effects,and thermal radiation.The primary aim is to examine how variations in nanoparticle size and nanolayer thickness affect the hydrothermal behavior of the nanofluid.The model also incorporates the contributions of viscous dissipation and Joule heating within the heat transfer equation.The governing momentum and energy equations are converted into dimensionless partial differential equations(PDEs)using appropriate similarity variables and are numerically solved using the finite element method(FEM)implemented in MATLAB.Extensive validation of this method confirms its reliability and accuracy in numerical solutions.The findings reveal that increasing the diameter of copper nanoparticles significantly enhances the velocity profile,with a more pronounced effect observed at wider inter-particle spacings.A higher solvent volume fraction leads to decreased velocity and temperature distributions,while a greater relaxation time ratio improves velocity and temperature profiles due to the increased elastic response of the fluid.Moreover,enhancements in the magnetic parameter,thermal radiation,and Eckert number lead to an elevation in temperature profiles.Furthermore,higher nanolayer thickness reduces the temperature profile,whereas particle radius yields the opposite outcome.
文摘Recent advancements in additive manufacturing(AM)have revolutionized the design and production of complex engineering microstructures.Despite these advancements,their mathematical modeling and computational analysis remain significant challenges.This research aims to develop an effective computational method for analyzing the free vibration of functionally graded(FG)microplates under high temperatures while resting on a Pasternak foundation(PF).This formulation leverages a new thirdorder shear deformation theory(new TSDT)for improved accuracy without requiring shear correction factors.Additionally,the modified couple stress theory(MCST)is incorporated to account for sizedependent effects in microplates.The PF is characterized by two parameters including spring stiffness(k_(w))and shear layer stiffness(k_(s)).To validate the proposed method,the results obtained are compared with those of the existing literature.Furthermore,numerical examples explore the influence of various factors on the high-temperature free vibration of FG microplates.These factors include the length scale parameter(l),geometric dimensions,material properties,and the presence of the elastic foundation.The findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the findings significantly enhance our comprehension of the free vibration of FG microplates in high thermal environments.In addition,the results of this research will have great potential in military and defense applications such as components of submarines,fighter aircraft,and missiles.
基金supported by the National Natural Science Foundation of China(Grant Nos.42272338 and 41827807)Department of Transportation of Zhejiang Province,China(Grant No.202213).
文摘Physics-informed neural networks(PINNs)have prevailed as differentiable simulators to investigate flow in porous media.Despite recent progress PINNs have achieved,practical geotechnical scenarios cannot be readily simulated because conventional PINNs fail in discontinuous heterogeneous porous media or multi-layer strata when labeled data are missing.This work aims to develop a universal network structure to encode the mass continuity equation and Darcy’s law without labeled data.The finite element approximation,which can decompose a complex heterogeneous domain into simpler ones,is adopted to build the differentiable network.Without conventional DNNs,physics-encoded finite element network(PEFEN)can avoid spectral bias and learn high-frequency functions with sharp/steep gradients.PEFEN rigorously encodes Dirichlet and Neumann boundary conditions without training.Benefiting from its discretized formulation,the discontinuous heterogeneous hydraulic conductivity is readily embedded into the network.Three typical cases are reproduced to corroborate PEFEN’s superior performance over conventional PINNs and the PINN with mixed formulation.PEFEN is sparse and demonstrated to be capable of dealing with heterogeneity with much fewer training iterations(less than 1/30)than the improved PINN with mixed formulation.Thus,PEFEN saves energy and contributes to low-carbon AI for science.The last two cases focus on common geotechnical settings of impermeable sheet pile in singlelayer and multi-layer strata.PEFEN solves these cases with high accuracy,circumventing costly labeled data,extra computational burden,and additional treatment.Thus,this study warrants the further development and application of PEFEN as a novel differentiable network in porous flow of practical geotechnical engineering.
文摘In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
基金This work was supported in part by the National Natural Science Foundation of China （Grant No. 11061021）, Natural Science Fund of Inner Mongolia Autonomous Region （2012MS0106, 2012MS0108）, Scientific Research Projection of Higher Schools of Inner Mongolia （NJZZ12011, N J10006, NJZY13199）, and the Program of Higherlevel talents of Inner Mongolia University （125119, 30105-125132）.
文摘This article presents a complete discretization of a nonlinear Sobolev equation using space-time discontinuous Galerkin method that is discontinuous in time and continuous in space. The scheme is formulated by introducing the equivalent integral equation of the primal equation. The proposed scheme does not explicitly include the jump terms in time, which represent the discontinuity characteristics of approximate solution. And then the complexity of the theoretical analysis is reduced. The existence and uniqueness of the approximate solution and the stability of the scheme are proved. The optimalorder error estimates in L2 （H1） and L2 （L2） norms are derived. These estimates are valid under weak restrictions on the space-time mesh, namely, without the condition kn ≥ ch2, which is necessary in traditional space-time discontinuous Galerkin methods. Numerical experiments are presented to verify the theoretical results.
