The phase-field model is a widely used mathematical approach for describing crack propagation in continuum damage fractures. In the context of phase field fracture simulations, adaptive finite element methods (AFEM) a...
The phase-field model is a widely used mathematical approach for describing crack propagation in continuum damage fractures. In the context of phase field fracture simulations, adaptive finite element methods (AFEM) are often employed to address the mesh size dependency of the model. However, existing AFEM approaches for this application frequently rely on heuristic adjustments and empirical parameters for mesh refinement. In this paper, we introduce an adaptive finite element method based on a recovery type posteriori error estimates approach grounded in theoretical analysis. This method transforms the gradient of the numerical solution into a smoother function space, using the difference between the recovered gradient and the original numerical gradient as an error indicator for adaptive mesh refinement. This enables the automatic capture of crack propagation directions without the need for empirical parameters. We have implemented this adaptive method for the Hybrid formulation of the phase-field model using the open-source software package FEALPy. The accuracy and efficiency of the proposed approach are demonstrated through simulations of classical 2D and 3D brittle fracture examples, validating the robustness and effectiveness of our implementation.
We propose some new weighted averaging methods for gradient recovery,and present analytical and numerical investigation on the performance of these weighted averaging *** is shown analytically that the harmonic averag...
详细信息
We propose some new weighted averaging methods for gradient recovery,and present analytical and numerical investigation on the performance of these weighted averaging *** is shown analytically that the harmonic averaging yields a superconvergent gradient for any mesh in one-dimension and the rectangular mesh in *** results indicate that these new weighted averaging methods are better recovered gradient approaches than the simple averaging and geometry averaging methods under triangular mesh.
In recent years,various kinds of cloak devices were designed by transforma-tion optics,but these cloak metamaterials are anisotropic and difficult to *** this paper,we designed the isotropic cloak metamaterials based o...
详细信息
In recent years,various kinds of cloak devices were designed by transforma-tion optics,but these cloak metamaterials are anisotropic and difficult to *** this paper,we designed the isotropic cloak metamaterials based on the numeri-cal method of the optimization theory to the inverse medium *** method has universality,and it is not limited by the shape and type of the cloak *** isotropic material is easier to manufacture in practice than anisotropic material.A large number of numerical results show the effectiveness of the method.
A new first-order optimality condition for the basis pursuit denoise (BPDN) problem is derived. This condition provides a new approach to choose the penalty param- eters adaptively for a fixed point iteration algori...
详细信息
A new first-order optimality condition for the basis pursuit denoise (BPDN) problem is derived. This condition provides a new approach to choose the penalty param- eters adaptively for a fixed point iteration algorithm. Meanwhile, the result is extended to matrix completion which is a new field on the heel of the compressed sensing. The numerical experiments of sparse vector recovery and low-rank matrix completion show validity of the theoretic results.
In this paper,we first prove the existence and uniqueness theorem of the solution of nonlinear variable-order fractional stochastic differential equations(VFSDEs).We futher constructe the Euler-Maruyama method to solv...
详细信息
In this paper,we first prove the existence and uniqueness theorem of the solution of nonlinear variable-order fractional stochastic differential equations(VFSDEs).We futher constructe the Euler-Maruyama method to solve the equations and prove the convergence in mean and the strong convergence of the *** particular,when the fractional order is no longer varying,the conclusions obtained are consistent with the relevant conclusions in the existing ***,the numerical experiments at the end of the article verify the correctness of the theoretical results obtained.
This paper mainly considers the optimal convergence analysis of the q-Maruyama method for stochastic Volterra integro-differential equations(SVIDEs)driven by Riemann-Liouville fractional Brownian motion under the glob...
详细信息
This paper mainly considers the optimal convergence analysis of the q-Maruyama method for stochastic Volterra integro-differential equations(SVIDEs)driven by Riemann-Liouville fractional Brownian motion under the global Lipschitz and linear growth ***,based on the contraction mapping principle,we prove the well-posedness of the analytical solutions of the ***,we show that the q-Maruyama method for the SVIDEs can achieve strong first-order *** particular,when the q-Maruyama method degenerates to the explicit Euler-Maruyama method,our result improves the conclusion that the convergence rate is H+1/2,H∈(0,1/2)by Yang et al.,***.,383(2021),***,the numerical experiment verifies our theoretical results.
By transforming the Caputo tempered fractional advection-diffusion equation into the Riemann–Liouville tempered fractional advection-diffusion equation,and then using the fractional-compact Grünwald–Letnikov te...
详细信息
By transforming the Caputo tempered fractional advection-diffusion equation into the Riemann–Liouville tempered fractional advection-diffusion equation,and then using the fractional-compact Grünwald–Letnikov tempered difference operator to approximate the Riemann–Liouville tempered fractional partial derivative,the fractional central difference operator to discritize the space Riesz fractional partial derivative,and the classical central difference formula to discretize the advection term,a numerical algorithm is constructed for solving the Caputo tempered fractional advection-diffusion *** stability and the convergence analysis of the numerical method are *** experiments show that the numerical method is effective.
A novel canonical Euler splitting method is proposed for nonlinear compositestiff functional differential-algebraic equations, the stability and convergence of themethod is evidenced, theoretical results are further c...
详细信息
A novel canonical Euler splitting method is proposed for nonlinear compositestiff functional differential-algebraic equations, the stability and convergence of themethod is evidenced, theoretical results are further confirmed by some numerical ***, the numerical method and its theories can be applied to specialcases, such as delay differential-algebraic equations and integral differential-algebraicequations.
Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are *** approaches are on account of two-grid skill include two major phases:find the ...
详细信息
Based on local algorithms,some parallel finite element(FE)iterative methods for stationary incompressible magnetohydrodynamics(MHD)are *** approaches are on account of two-grid skill include two major phases:find the FE solution by solving the nonlinear system on a globally coarse mesh to seize the low frequency component of the solution,and then locally solve linearized residual subproblems by one of three iterations(Stokes-type,Newton,and Oseen-type)on subdomains with fine grid in parallel to approximate the high frequency *** error estimates with regard to two mesh sizes and iterative steps of the proposed algorithms are *** numerical examples are implemented to verify the algorithm.
Four primal discontinuous Galerkin methods are applied to solve reactive transport problems, namely, Oden-BabuSka-Baumann DG (OBB-DG), non-symmetric interior penalty Galerkin (NIPG), symmetric interior penalty Gal...
详细信息
Four primal discontinuous Galerkin methods are applied to solve reactive transport problems, namely, Oden-BabuSka-Baumann DG (OBB-DG), non-symmetric interior penalty Galerkin (NIPG), symmetric interior penalty Galerkin (SIPG), and incomplete interior penalty Galerkin (IIPG). A unified a posteriori residual-type error estimation is derived explicitly for these methods. From the computed solution and given data, explicit estimators can be computed efficiently and directly, which can be used as error indicators for adaptation. Unlike in the reference [10], we obtain the error estimators in L^2 (L^2) norm by using duality techniques instead of in L^2(H^1) norm.
暂无评论