Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery...Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery.We rst derive the projection formulas for a vector onto the feasible sets.The centralized circumcentered-reection method is designed to solve the convex feasibility problem.Some numerical experiments demonstrate the feasibility and e ectiveness of the proposed algorithm,showing superior performance compared to conventional alternating projection methods.展开更多
To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fract...To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.展开更多
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 order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error...In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.展开更多
Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout a...Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.展开更多
In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least s...In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least squares(IMLS)method,which enhances the efficiency and stability of the numerical solution.The numerical solution formulas are derived using the improved element-free Galerkin(IEFG)method.We introduce the solid isotropic microstructures with penalization(SIMP)model to formulate a mathematical model for topology opti-mization,which effectively penalizes intermediate densities.The optimization problem is defined with the numerical solution formula and volume fraction as constraints.The objective function,which is the minimum value of flexibility,is optimized iteratively using the optimization criterion method to update the design variables efficiently and converge to an optimal solution.Sensitivity analysis is performed using the adjoint method,which provides accurate and efficient gradient information for the optimization algorithm.We validate the proposed framework through a series of numerical examples,including clamped beam,cantilever beam,and simply supported beam made of orthotropic materials.The convergence of the objective function is demonstrated by increasing the number of iterations.Additionally,the stability of the iterative process is analyzed by examining the fluctuation law of the volume fraction.By adjusting the parameters to an appropriate range,we achieve the final optimization results of the IEFG method without the checkerboard phenomenon.Comparative studies between the Element-Free Galerkin(EFG)and IEFG methods reveal that both methods yield consistent optimization results under identical parameter settings.However,the IEFG method significantly reduces computational time,highlighting its efficiency and suitability for orthotropic materials.展开更多
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.展开更多
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.展开更多
In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this sche...In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this scheme seamlessly combines the material point method(MPM)for accurately tracking solid-phase deformation and heat transport,and the Eulerian finite element method(FEM)for effectively capturing fluid flow and heat advection-diffusion behavior.The proposed approach circumvents the substantial challenges posed by large nonlinear equation systems with the monolithic solution scheme.The staggered solution process strategically separates each physical field through explicit or implicit integration.The characteristic-based method is used to stabilize advection-dominated heat flows for efficient numerical implementation.Furthermore,a fractional step approach is employed to decompose fluid velocity and pressure,thereby suppressing pore pressure oscillation on the linear background grid.The fracturing initiation and propagation are simulated by a rate-dependent phase field model.Through a series of quasi-static and transient simulations,the exceptional performance and promising potential of the proposed model in addressing THM fracturing problems in poro-elastic media is demonstrated.展开更多
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.展开更多
This paper presents a novel element differential method for modeling cracks in piezoelectric materials,aiming to simulate fracture behaviors and predict the fracture parameter known as the J-integral accurately.The me...This paper presents a novel element differential method for modeling cracks in piezoelectric materials,aiming to simulate fracture behaviors and predict the fracture parameter known as the J-integral accurately.The method leverages an efficient collocation technique to satisfy traction and electric charge equilibrium on the crack surface,aligning internal nodes with piezoelectric governing equations without needing integration or variational principles.It combines the strengths of the strong form collocation and finite element methods.The J-integral is derived analytically using the equivalent domain integral method,employing Green's formula and Gauss's divergence theorem to transform line integrals into area integrals for solving two-dimensional piezoelectric material problems.The accuracy of the method is validated through comparison with three typical examples,and it offers fracture prevention strategies for engineering piezoelectric structures under different electrical loading patterns.展开更多
In essence,the negotiation of license fees on standard essential patent(SEP)belongs to a kind of market be⁃havior,and the pricing right should be given to the market subjects under the requirements of patent law.In re...In essence,the negotiation of license fees on standard essential patent(SEP)belongs to a kind of market be⁃havior,and the pricing right should be given to the market subjects under the requirements of patent law.In recent years,the frequent disputes on SEP license fees witnessed in the industrial and academic worlds,together with the lack of systematic supporting functions like FRAND,make SEP pricing excessively reliant on judicial judgment in practice.Fortunately,a variety of pricing methods have been proposed by theoretical research and practiced in judicial cases,which provide possible solutions for the license fee pricing of SEP from the operational level.In this paper,by focusing on the characteristics of the existing SEP pricing methods in the academic fields and judicial system,the dispute caused by license fees of SEP is clarified firstly,then by combining and interpreting twelve existing pricing methods of license fee of SEP with academic literature and judicial cases,four categories of methods are composed based on the application stages and calculation logic.Thirdly,the application barriers and dilemmas caused by the inherent limita⁃tions of the four categories of methods are analyzed,and the possible ways to put these methods into practice are ex⁃plored.Lastly,suggestions are presented from the aspects of preconditions for application,pricing stages,dispute reso⁃lution mechanisms,and comprehensive applications.The purpose of this paper is to provide enlightenment for getting back on track with the pricing right and further optimization of the pricing mechanism of license fees of SEP.展开更多
This paper investigates the frictional adhesive contact of a rigid,electrically/magnetically conductive spherical indenter sliding past a multiferroic coating deposed onto a rigid substrate,based on the hybrid element...This paper investigates the frictional adhesive contact of a rigid,electrically/magnetically conductive spherical indenter sliding past a multiferroic coating deposed onto a rigid substrate,based on the hybrid element method.The adhesion behavior is described based on the Maugis-Dugdale model.The adhesion-driven conjugate gradient method is employed to calculate the distribution of unknown pressures,while the discrete convolution-fast Fourier transform is utilized to compute the deformations,surface electric and magnetic potentials as well as the subsurface stresses,electric displacements,and magnetic inductions.The goal of this study is to investigate the influences of adhesion parameter,friction coefficient,coating thickness,and surface electric and magnetic charge densities on contact behaviors,such as contact area and pressures,electric and magnetic potentials,and subsurface stresses.展开更多
This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solv...This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solving a specific numerical problem under the scope of the linear finite element method(LFEM),so the method is termed computational method for analytical solutions with finite elements(CMAS-FE).The primary objective of the CMAS-FE is to construct analytical expressions for displacements and reaction forces at nodes,as well as for strains and stresses at elemental quadrature points,all of which are formulated as infinite series solutions of various orders of Poisson’s ratios.Like the conventional LFEM,the CMAS-FE forms global sparse linear equations,but the Young’s modulus and Poisson’s ratio remain variables(or symbols).By employing a direct inverse method to solve these symbolic linear systems,an analytical expression of the displacement field can be constructed.The CMAS-FE is validated via patch and bending tests,which demonstrate convergence with mesh and term refine-ment.Furthermore,the CMAS-FE is applied to obtain the bending stiffness of a beam structure and to estimate an approximate stress intensity factor for a straight crack within a square-shaped plate.展开更多
The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRE...The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRESP),is an important factor affecting the homogeneity.This study investigated the influence of the size combination and MRESP on the distribution homogeneity of particles in aggregate blend using discrete element method(DEM).An indicator quantifying the distribution homogeneity was established according to the coefficient of variation(CV)for particle number.Two-size,three-size,and four-size aggregate blends with various compositions were designed.Laboratory tests show the DEM simulation is feasible.The particle distribution homogeneity in various blends was analyzed.The results showed the distribution homogeneity of each size particles in a blend is closely related to their mass fraction.The higher the mass fraction of the particles,the more homogeneous the distribution of them.The MRESP has no significant influence on the homogeneity of the blend composed of only coarse aggregates.However,the homogeneity of the blend composed of coarse and fine aggregates improves gradually with the increase of the mass fraction of fine aggregates.The smaller the maximum particle size in a blend,the better the homogeneity.It is suggested that the mass fraction of fine aggregates should be between 33%and 50%for achieving good homogeneity of aggregate blends.The research results can provide a reference for gradation design of asphalt mixture.展开更多
This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic...This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic simulations and its high computational cost in broadband problems are examined.To improve efficiency,a Taylor series expansion is applied to decouple frequency-dependent terms in BEM.Additionally,the SecondOrder Arnoldi(SOAR)model order reduction method is integrated to reduce computational costs and enhance numerical stability.Furthermore,an isogeometric sensitivity boundary integral equation is formulated using the direct differentiation method,incorporating Cauchy principal value integrals and Hadamard finite part integrals to handle singularities.The proposed method improves the computational efficiency,and the acoustic sensitivity analysis provides theoretical support for further acoustic structure optimization.展开更多
We present a hybrid smoothed particle magnetohydrodynamics(SPMHD)code integrating smoothed particle hydrodynamics(SPH)and finite element methods(FEM)to simulate coupled fluid-electromagnetic phenomena.The framework em...We present a hybrid smoothed particle magnetohydrodynamics(SPMHD)code integrating smoothed particle hydrodynamics(SPH)and finite element methods(FEM)to simulate coupled fluid-electromagnetic phenomena.The framework employs SPH for fluid dynamics,addressing large deformations,shocks,and plasma behavior,while FEM resolves electromagnetic fields via Maxwell's equations for magnetic vector and electric scalar potentials,ensuring divergence-free conditions and global current density calculations in conductive region.Operator splitting method couples these modules,enabling real-time integration of magnetic,electric,thermal,and fluid fields.Benchmark tests validate the code against analytical solutions and existing models,including blow-by instability simulations that demonstrate the method's accuracy in capturing fluid-magnetic interactions.Designed for 3D applications,SPMHD offers robust scalability across multiprocessor architectures,establishing it as a versatile tool for plasma physics research.展开更多
Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the chall...Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale reservoirs.展开更多
A parallel finite element scheme for 3D resistivity method forward modeling is introduced in this article.The domain decomposition algorithm,along with a message passing interface,is used to implement parallelism.The ...A parallel finite element scheme for 3D resistivity method forward modeling is introduced in this article.The domain decomposition algorithm,along with a message passing interface,is used to implement parallelism.The computational domain is divided into subdomains,and mesh partitioning is combined with load balancing.Unstructured meshes and local mesh refinement strategies are used to realize high precision for complex topography models.Furthermore,an improved linear solver for multi-electrode resistivity method modeling is adopted.Recycling preconditioned conjugate gradient,which is a linear solver,is based on the similarity of linear systems between point sources.The multiple right-hand-side linear systems corresponding to different point source positions are constructed,and the accelerated convergence is obtained through recycling subspace using the linear solver.The computational accuracy and efficiency of the forward scheme for complex topography models are verified using the numerical test results.展开更多
Highly efficient removal of impurity Bi element from scrap brass can facilitate the recycling process of brass.The effects of melting temperature,holding time and Mg-Ca alloy content on the removal effect of impurity ...Highly efficient removal of impurity Bi element from scrap brass can facilitate the recycling process of brass.The effects of melting temperature,holding time and Mg-Ca alloy content on the removal effect of impurity Bi element were investigated by compound-separation method.The mechanism of the compound-separation method was revealed for removing the Bi element from a thermodynamic point of view.The results showed that the Bi content was decreased from 1.95 wt.%to 0.178 wt.%at the optimum process parameters of melting temperature of 980°C,holding time of 20 min,and Mg-Ca alloy content of 6 wt.%,achieving a removal rate of 90.9%.A small amount of Ca-Bi compound remained in the brass matrix after refining.NaF flux can effectively wet and adsorb Ca-Bi compounds due to its low viscosity and the function of lowering the surface tension,which facilitate the agglomeration and flotation of Ca-Bi compounds to the melt surface,thereby ensuring the sufficient removal of Bi element.展开更多
基金Supported by the Natural Science Foundation of Guangxi Province(Grant Nos.2023GXNSFAA026067,2024GXN SFAA010521)the National Natural Science Foundation of China(Nos.12361079,12201149,12261026).
文摘Convex feasibility problems are widely used in image reconstruction,sparse signal recovery,and other areas.This paper is devoted to considering a class of convex feasibility problem arising from sparse signal recovery.We rst derive the projection formulas for a vector onto the feasible sets.The centralized circumcentered-reection method is designed to solve the convex feasibility problem.Some numerical experiments demonstrate the feasibility and e ectiveness of the proposed algorithm,showing superior performance compared to conventional alternating projection methods.
基金funded by the project of the Major Scientific and Technological Projects of CNOOC in the 14th Five-Year Plan(No.KJGG2022-0701)the CNOOC Research Institute(No.2020PFS-03).
文摘To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.
基金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.
文摘In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.
基金supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.
基金supported by the Graduate Student Scientific Research Innovation Project through Research Innovation Fund for Graduate Students in Shanxi Province(Project No.2024KY648).
文摘In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least squares(IMLS)method,which enhances the efficiency and stability of the numerical solution.The numerical solution formulas are derived using the improved element-free Galerkin(IEFG)method.We introduce the solid isotropic microstructures with penalization(SIMP)model to formulate a mathematical model for topology opti-mization,which effectively penalizes intermediate densities.The optimization problem is defined with the numerical solution formula and volume fraction as constraints.The objective function,which is the minimum value of flexibility,is optimized iteratively using the optimization criterion method to update the design variables efficiently and converge to an optimal solution.Sensitivity analysis is performed using the adjoint method,which provides accurate and efficient gradient information for the optimization algorithm.We validate the proposed framework through a series of numerical examples,including clamped beam,cantilever beam,and simply supported beam made of orthotropic materials.The convergence of the objective function is demonstrated by increasing the number of iterations.Additionally,the stability of the iterative process is analyzed by examining the fluctuation law of the volume fraction.By adjusting the parameters to an appropriate range,we achieve the final optimization results of the IEFG method without the checkerboard phenomenon.Comparative studies between the Element-Free Galerkin(EFG)and IEFG methods reveal that both methods yield consistent optimization results under identical parameter settings.However,the IEFG method significantly reduces computational time,highlighting its efficiency and suitability for orthotropic materials.
基金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.
基金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.
基金supported by National Natural Science Foundation of China(Grant No.42377149)the Research Grants Council of Hong Kong(General Research Fund Project No.17202423).
文摘In this study,a powerful thermo-hydro-mechanical(THM)coupling solution scheme for saturated poroelastic media involving brittle fracturing is developed.Under the local thermal non-equilibrium(LTNE)assumption,this scheme seamlessly combines the material point method(MPM)for accurately tracking solid-phase deformation and heat transport,and the Eulerian finite element method(FEM)for effectively capturing fluid flow and heat advection-diffusion behavior.The proposed approach circumvents the substantial challenges posed by large nonlinear equation systems with the monolithic solution scheme.The staggered solution process strategically separates each physical field through explicit or implicit integration.The characteristic-based method is used to stabilize advection-dominated heat flows for efficient numerical implementation.Furthermore,a fractional step approach is employed to decompose fluid velocity and pressure,thereby suppressing pore pressure oscillation on the linear background grid.The fracturing initiation and propagation are simulated by a rate-dependent phase field model.Through a series of quasi-static and transient simulations,the exceptional performance and promising potential of the proposed model in addressing THM fracturing problems in poro-elastic media is demonstrated.
基金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.
基金Financial support of this work by the Technology Development program of China(Grant No.2022204B003)National Natural Science Foundation of China(12272083 and 12172078)the Fundamental Research Funds for the Central Universities(DUT24YJ136)is gratefully acknowledged.
文摘This paper presents a novel element differential method for modeling cracks in piezoelectric materials,aiming to simulate fracture behaviors and predict the fracture parameter known as the J-integral accurately.The method leverages an efficient collocation technique to satisfy traction and electric charge equilibrium on the crack surface,aligning internal nodes with piezoelectric governing equations without needing integration or variational principles.It combines the strengths of the strong form collocation and finite element methods.The J-integral is derived analytically using the equivalent domain integral method,employing Green's formula and Gauss's divergence theorem to transform line integrals into area integrals for solving two-dimensional piezoelectric material problems.The accuracy of the method is validated through comparison with three typical examples,and it offers fracture prevention strategies for engineering piezoelectric structures under different electrical loading patterns.
基金Hierarchical Identification and Cross-Layer Correlation of Key Core Technologies from the Perspective of Industrial Chain Structure (National Social Science Fund of China, 24BTQ067)Chongqing Education Commission (CEC) Funding:Research on the Co-governance Mechanism of Patent Quality Based on the Dual-Filter Perspective(24SKGH213)Chongqing Graduate Education and Teaching Funding:Research on the Interdisciplinary Law of Intellectual Property and Optimization of Graduate Talent Training Mode(yjg213122)。
文摘In essence,the negotiation of license fees on standard essential patent(SEP)belongs to a kind of market be⁃havior,and the pricing right should be given to the market subjects under the requirements of patent law.In recent years,the frequent disputes on SEP license fees witnessed in the industrial and academic worlds,together with the lack of systematic supporting functions like FRAND,make SEP pricing excessively reliant on judicial judgment in practice.Fortunately,a variety of pricing methods have been proposed by theoretical research and practiced in judicial cases,which provide possible solutions for the license fee pricing of SEP from the operational level.In this paper,by focusing on the characteristics of the existing SEP pricing methods in the academic fields and judicial system,the dispute caused by license fees of SEP is clarified firstly,then by combining and interpreting twelve existing pricing methods of license fee of SEP with academic literature and judicial cases,four categories of methods are composed based on the application stages and calculation logic.Thirdly,the application barriers and dilemmas caused by the inherent limita⁃tions of the four categories of methods are analyzed,and the possible ways to put these methods into practice are ex⁃plored.Lastly,suggestions are presented from the aspects of preconditions for application,pricing stages,dispute reso⁃lution mechanisms,and comprehensive applications.The purpose of this paper is to provide enlightenment for getting back on track with the pricing right and further optimization of the pricing mechanism of license fees of SEP.
基金support from the National Natural Science Foundation of China(12102085)the Postdoctoral Science Foundation of China(2023M730504)+2 种基金the Sichuan Province Regional Innovation and Cooperation Project(2024YFHZ0210)supported by the European Union-NextGenerationEU through the Italian Ministry of University and Research under the following programs:(NM)PRIN2022(Projects of Relevant National Interest)grant no.2022SJ8HTC-Electroactive Gripper for Micro-Object Manipulation(ELFIN)(NM)PRIN2022 PNRR(Projects of Relevant National Interest)grant no.P2022MAZHX-Tribological Modeling for Sustainable Design of Industrial Frictional Interfaces(TRIBOSCORE).
文摘This paper investigates the frictional adhesive contact of a rigid,electrically/magnetically conductive spherical indenter sliding past a multiferroic coating deposed onto a rigid substrate,based on the hybrid element method.The adhesion behavior is described based on the Maugis-Dugdale model.The adhesion-driven conjugate gradient method is employed to calculate the distribution of unknown pressures,while the discrete convolution-fast Fourier transform is utilized to compute the deformations,surface electric and magnetic potentials as well as the subsurface stresses,electric displacements,and magnetic inductions.The goal of this study is to investigate the influences of adhesion parameter,friction coefficient,coating thickness,and surface electric and magnetic charge densities on contact behaviors,such as contact area and pressures,electric and magnetic potentials,and subsurface stresses.
基金supported by the National Natural Science Foundation of China Excellence Research Group Program for“Multiscale Problems in Nonlinear Mechanics”(Grant No.12588201)the National Key R&D Program of China(Grant No.2023YFA1008901)+1 种基金the National Nat-ural Science Foundation of China(Grant No.12172009)supported by“The Fundamental Research Funds for the Central Universities,Peking University”.
文摘This study presents a novel methodology to obtain an approximate analytical solution for an isotropic homo-geneous elastic medium with displacement and traction boundary conditions.The solution is derived through solving a specific numerical problem under the scope of the linear finite element method(LFEM),so the method is termed computational method for analytical solutions with finite elements(CMAS-FE).The primary objective of the CMAS-FE is to construct analytical expressions for displacements and reaction forces at nodes,as well as for strains and stresses at elemental quadrature points,all of which are formulated as infinite series solutions of various orders of Poisson’s ratios.Like the conventional LFEM,the CMAS-FE forms global sparse linear equations,but the Young’s modulus and Poisson’s ratio remain variables(or symbols).By employing a direct inverse method to solve these symbolic linear systems,an analytical expression of the displacement field can be constructed.The CMAS-FE is validated via patch and bending tests,which demonstrate convergence with mesh and term refine-ment.Furthermore,the CMAS-FE is applied to obtain the bending stiffness of a beam structure and to estimate an approximate stress intensity factor for a straight crack within a square-shaped plate.
基金funded by the National Natural Science Foundation of China(No.51978048).
文摘The homogeneity of aggregate blend has a significant influence on the performance of asphalt mixture.The composition of aggregate blend,including the size combination and the mass ratio between each size particles(MRESP),is an important factor affecting the homogeneity.This study investigated the influence of the size combination and MRESP on the distribution homogeneity of particles in aggregate blend using discrete element method(DEM).An indicator quantifying the distribution homogeneity was established according to the coefficient of variation(CV)for particle number.Two-size,three-size,and four-size aggregate blends with various compositions were designed.Laboratory tests show the DEM simulation is feasible.The particle distribution homogeneity in various blends was analyzed.The results showed the distribution homogeneity of each size particles in a blend is closely related to their mass fraction.The higher the mass fraction of the particles,the more homogeneous the distribution of them.The MRESP has no significant influence on the homogeneity of the blend composed of only coarse aggregates.However,the homogeneity of the blend composed of coarse and fine aggregates improves gradually with the increase of the mass fraction of fine aggregates.The smaller the maximum particle size in a blend,the better the homogeneity.It is suggested that the mass fraction of fine aggregates should be between 33%and 50%for achieving good homogeneity of aggregate blends.The research results can provide a reference for gradation design of asphalt mixture.
基金supported by the Shanxi Scholarship Council of China(Grant No.2023-036)the Natural Science Foundation of Shanxi Province(Grant No.202303021222020).
文摘This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic simulations and its high computational cost in broadband problems are examined.To improve efficiency,a Taylor series expansion is applied to decouple frequency-dependent terms in BEM.Additionally,the SecondOrder Arnoldi(SOAR)model order reduction method is integrated to reduce computational costs and enhance numerical stability.Furthermore,an isogeometric sensitivity boundary integral equation is formulated using the direct differentiation method,incorporating Cauchy principal value integrals and Hadamard finite part integrals to handle singularities.The proposed method improves the computational efficiency,and the acoustic sensitivity analysis provides theoretical support for further acoustic structure optimization.
基金supported by the Major National Science and Technology Infrastructure(No.2208-000000-04-01249628)the Shanghai Science and Technology Commission(No.21DZ1206500)。
文摘We present a hybrid smoothed particle magnetohydrodynamics(SPMHD)code integrating smoothed particle hydrodynamics(SPH)and finite element methods(FEM)to simulate coupled fluid-electromagnetic phenomena.The framework employs SPH for fluid dynamics,addressing large deformations,shocks,and plasma behavior,while FEM resolves electromagnetic fields via Maxwell's equations for magnetic vector and electric scalar potentials,ensuring divergence-free conditions and global current density calculations in conductive region.Operator splitting method couples these modules,enabling real-time integration of magnetic,electric,thermal,and fluid fields.Benchmark tests validate the code against analytical solutions and existing models,including blow-by instability simulations that demonstrate the method's accuracy in capturing fluid-magnetic interactions.Designed for 3D applications,SPMHD offers robust scalability across multiprocessor architectures,establishing it as a versatile tool for plasma physics research.
文摘Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale reservoirs.
基金supported by the National Natural Science Foundation of China(No:42274182)Guangxi Natural Science Foundation(No:2020 GXNSFAA297079)。
文摘A parallel finite element scheme for 3D resistivity method forward modeling is introduced in this article.The domain decomposition algorithm,along with a message passing interface,is used to implement parallelism.The computational domain is divided into subdomains,and mesh partitioning is combined with load balancing.Unstructured meshes and local mesh refinement strategies are used to realize high precision for complex topography models.Furthermore,an improved linear solver for multi-electrode resistivity method modeling is adopted.Recycling preconditioned conjugate gradient,which is a linear solver,is based on the similarity of linear systems between point sources.The multiple right-hand-side linear systems corresponding to different point source positions are constructed,and the accelerated convergence is obtained through recycling subspace using the linear solver.The computational accuracy and efficiency of the forward scheme for complex topography models are verified using the numerical test results.
基金the financial supplies supported by the National Natural Science Foundation of China(Nos.U2202255,52371038)the Science and Technology Innovation Program of Hunan Province,China(No.2023RC1019)。
文摘Highly efficient removal of impurity Bi element from scrap brass can facilitate the recycling process of brass.The effects of melting temperature,holding time and Mg-Ca alloy content on the removal effect of impurity Bi element were investigated by compound-separation method.The mechanism of the compound-separation method was revealed for removing the Bi element from a thermodynamic point of view.The results showed that the Bi content was decreased from 1.95 wt.%to 0.178 wt.%at the optimum process parameters of melting temperature of 980°C,holding time of 20 min,and Mg-Ca alloy content of 6 wt.%,achieving a removal rate of 90.9%.A small amount of Ca-Bi compound remained in the brass matrix after refining.NaF flux can effectively wet and adsorb Ca-Bi compounds due to its low viscosity and the function of lowering the surface tension,which facilitate the agglomeration and flotation of Ca-Bi compounds to the melt surface,thereby ensuring the sufficient removal of Bi element.
