共查询到20条相似文献,搜索用时 15 毫秒
1.
Daniel Kessler;Ricardo H. Nochetto;Alfred Schmidt 《Mathematical Modelling and Numerical Analysis》2004,38(1):129-142
Phase-field models, the simplest of which is Allen–Cahn's problem, are characterized by a small parameter ε that dictates the interface thickness. These models naturally call for mesh adaptation techniques, which rely on aposteriori error control. However, their error analysis usually deals with the underlying non-monotone nonlinearity via a Gronwall argument which leads to an exponential dependence on ε-2. Using an energy argument combined with a topological continuation argument and a spectral estimate, we establish an aposteriori error control result with only a low order polynomial dependence in ε-1. Our result is applicable to any conforming discretization technique that allows for a posteriori residual estimation. Residual estimators for an adaptive finite element scheme are derived to illustrate the theory.https://doi.org/10.1051/m2an:2004006 相似文献
2.
Shuaichao Pei Yanren Hou Wenjing Yan 《Numerical Methods for Partial Differential Equations》2022,38(1):65-101
We consider numerical approximations for a modified phase field crystal model with a strong nonlinear vacancy potential. Based on the invariant energy quadratization approach and stabilized strategies, we develop linear, unconditionally energy stable numerical schemes using the first-order Euler method, the second-order backward differentiation formulas and the second-order Crank–Nicolson method, respectively. We rigorously prove the unconditional energy stability, the mass conservation of these three numerical schemes and carry out error estimates in time for the first-order numerical scheme. Various numerical experiments in 2D and 3D are carried out to validate the accuracy, energy stability, mass conservation, and efficiency of the proposed schemes. 相似文献
3.
We introduce a piecewise linear finite-element scheme with semi-implicittime discretization for an evolutionary phase field system modellingthe isothermal solidification process of a binary alloy. Thissystem can be written in a vectorial form as a nonlinear parabolicsystem. The convergence of the scheme with error estimate isthen proved by introducing a generalized vectorial ellipticprojector. 相似文献
4.
崔明 《高等学校计算数学学报》2002,24(3):206-211
1 引 言设Ω R2为具有光滑边界的有界区域,考虑非定常的,无量纲化的,而且带有热传导的粘性不可压缩流体力学问题:问题(Ⅰ):求u=(u1,u2),p,T满足: 相似文献
5.
Some scalar double-well problems eventually lead to a degenerate convex minimization problem with unique stress. We consider the adaptive conforming and nonconforming finite element methods for the scalar double-well problem with the reliable a posteriori error analysis. A number of experiments confirm the effective decay rates of the methods. 相似文献
6.
Averaging techniques are popular tools in adaptive finite element methods for the numerical treatment of second order partial differential equations since they provide efficient a posteriori error estimates by a simple postprocessing. In this paper, their reliablility is shown for conforming, nonconforming, and mixed low order finite element methods in a model situation: the Laplace equation with mixed boundary conditions. Emphasis is on possibly unstructured grids, nonsmoothness of exact solutions, and a wide class of averaging techniques. Theoretical and numerical evidence supports that the reliability is up to the smoothness of given right-hand sides.
7.
8.
P. Bringmann C. Carstensen C. Merdon 《Numerical Methods for Partial Differential Equations》2016,32(5):1411-1432
The pseudostress approximation of the Stokes equations rewrites the stationary Stokes equations with pure (but possibly inhomogeneous) Dirichlet boundary conditions as another (equivalent) mixed scheme based on a stress in H(div) and the velocity in L2. Any standard mixed finite element function space can be utilized for this mixed formulation, e.g., the Raviart‐Thomas discretization which is related to the Crouzeix‐Raviart nonconforming finite element scheme in the lowest‐order case. The effective and guaranteed a posteriori error control for this nonconforming velocity‐oriented discretization can be generalized to the error control of some piecewise quadratic velocity approximation that is related to the discrete pseudostress. The analysis allows for local inf‐sup constants which can be chosen in a global partition to improve the estimation. Numerical examples provide strong evidence for an effective and guaranteed error control with very small overestimation factors even for domains with large anisotropy.© 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1411–1432, 2016 相似文献
9.
本文针对线性对流占优扩散方程提出了一种新型数值模拟方法一扩展特征混合有限元法,即对对流部分沿特征线方向离散,而对扩散部分采用扩展混合有限元方法,同时高精度逼近未知函数,未知函数的梯度及伴随向量函数,通过严格的数值分析,得到其最优L^2模误差估计。 相似文献
10.
Yan Yang & Xiaoping Xie 《高等学校计算数学学报(英文版)》2022,15(1):227-250
A family of conforming mixed finite elements with mass lumping on triangular grids are presented for linear elasticity. The stress field is approximatedby symmetric $H($div) − $P_k (k ≥ 3)$ polynomial tensors enriched with higher orderbubbles so as to allow mass lumping, and the displacement field is approximated by $C^{−1}− P_{k−1}$ polynomial vectors enriched with higher order terms. For both the proposed mixed elements and their mass lumping schemes, optimal error estimates arederived for the stress and displacement in $H$(div) norm and $L^2$ norm, respectively.Numerical results confirm the theoretical analysis. 相似文献
11.
Averaging techniques are popular tools in adaptive finite element methods since they provide efficient a posteriori error estimates by a simple postprocessing. In the second paper of our analysis of their reliability, we consider conforming -FEM of higher (i.e., not of lowest) order in two or three space dimensions. In this paper, reliablility is shown for conforming higher order finite element methods in a model situation, the Laplace equation with mixed boundary conditions. Emphasis is on possibly unstructured grids, nonsmoothness of exact solutions, and a wide class of local averaging techniques. Theoretical and numerical evidence supports that the reliability is up to the smoothness of given right-hand sides.
12.
1 引 言在地下水含水层中 ,污染物随地下水运移并常常发生各种化学反应 [1 ] .描述地下水含水层中一类阳离子交换反应 m M1 +r M2 k2k1 r M2 +m M1 的数学模型[2 ] 为 : s1 t- dΔs1 =f1 , x∈Ω ,t∈ J (1.1a) s2 t- dΔs2 =f2 , x∈Ω ,t∈ J (1.1b) c1 t+ρ s1 t- DΔc1 =0 , x∈Ω ,t∈ J (1.2 a) c2 t+ρ s2 t- DΔc2 =0 , x∈Ω ,t∈ J (1.2 b)其中 Ω R2为具有光滑边界的有界区域 ,J=(0 ,T].这里 D>d>0为扩散系数 ,ρ>0为固体颗粒密度 ,均为常数 .根据 [1,2 ]应有 :f1 =m[k1 Rm1 Rr2 cm1 sr2 -… 相似文献
13.
14.
We consider an augmented mixed finite element method applied to the linear elasticity problem and derive a posteriori error estimators that are simpler and easier to implement than the ones available in the literature. In the case of homogeneous Dirichlet boundary conditions, the new a posteriori error estimator is reliable and locally efficient, whereas for non-homogeneous Dirichlet boundary conditions, we derive an a posteriori error estimator that is reliable and satisfies a quasi-efficiency bound. Numerical experiments illustrate the performance of the corresponding adaptive algorithms and support the theoretical results. 相似文献
15.
This article aims to study the partitioned method for magnetohydrodynamic flows at small magnetic Reynolds numbers. We design a partitioned second‐order method and show that this method is stable under a time step () restrict condition. Our method can decouple the magnetohydrodynamic equations so that we can solve two relatively simple subproblems separately at each time step, which is computationally economic. A complete theoretical analysis of error estimates is also given. Finally, we present numerical experiments to support our theory.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1966–1986, 2017 相似文献
16.
In this study, a fully discrete defect correction finite element method for the unsteady incompressible Magnetohydrodynamics (MHD) equations, which is leaded by combining the Back Euler time discretization with the two-step defect correction in space, is presented. It is a continuous work of our formal paper [Math Method Appl Sci. 2017. DOI:10.1002/mma.4296]. The defect correction method is an iterative improvement technique for increasing the accuracy of a numerical solution without applying a grid refinement. Firstly, the nonlinear MHD equation is solved with an artificial viscosity term. Then, the numerical solutions are improved on the same grid by a linearized defect-correction technique. Then, we introduce the numerical analysis including stability analysis and error analysis. The numerical analysis proves that our method is stable and has an optimal convergence rate. Some numerical results [see Math Method Appl Sci. 2017. DOI:10.1002/mma.4296] show that this method is highly efficient for the unsteady incompressible MHD problems. 相似文献
17.
This paper discusses convergence and complexity of arbitrary,but fixed,order adaptive mixed element methods for the Poisson equation in two and three dimensions.The two main ingredients in the analysis,namely the quasi-orthogonality and the discrete reliability,are achieved by use of a discrete Helmholtz decomposition and a discrete inf-sup condition.The adaptive algorithms are shown to be contractive for the sum of the error of flux in L2-norm and the scaled error estimator after each step of mesh refinement and to be quasi-optimal with respect to the number of elements of underlying partitions.The methods do not require a separate treatment for the data oscillation. 相似文献
18.
19.
1 IntroductionInrecentyears,theintentionaloraccidentalreleaseofthechemicalwastesonsoilshasfurtherstimulatedcurrentinterestsinthemovementofchemicals.Displacementstudieshavebecomeimportanttoolsinsoilphysics,particularlyforpredictingthemovementofpestcides… 相似文献
20.
This paper aims at a general guideline to obtain a posteriori error estimates for the finite element error control in computational partial differential equations. In the abstract setting of mixed formulations, a generalised formulation of the corresponding residuals is proposed which then allows for the unified estimation of the respective dual norms. Notably, this can be done with an approach which is applicable in the same way to conforming, nonconforming and mixed discretisations. Subsequently, the unified approach is applied to various model problems. In particular, we consider the Laplace, Stokes, Navier-Lamé, and the semi-discrete eddy current equations. 相似文献