ترغب بنشر مسار تعليمي؟ اضغط هنا

A second-order accurate, operator splitting scheme for reaction-diffusion systems in an energetic variational formulation

123   0   0.0 ( 0 )
 نشر من قبل Yiwei Wang
 تاريخ النشر 2021
  مجال البحث الهندسة المعلوماتية
والبحث باللغة English




اسأل ChatGPT حول البحث

A second-order accurate in time, positivity-preserving, and unconditionally energy stable operator splitting numerical scheme is proposed and analyzed for the system of reaction-diffusion equations with detailed balance. The scheme is designed based on an energetic variational formulation, in which the reaction part is reformulated in terms of the reaction trajectory, and both the reaction and diffusion parts dissipate the same free energy. At the reaction stage, the reaction trajectory equation is approximated by a second-order Crank-Nicolson type method. The unique solvability, positivity-preserving, and energy-stability are established based on a convexity analysis. In the diffusion stage, an exact integrator is applied if the diffusion coefficients are constant, and a Crank-Nicolson type scheme is applied if the diffusion process becomes nonlinear. In either case, both the positivity-preserving property and energy stability could be theoretically established. Moreover, a combination of the numerical algorithms at both stages by the Strang splitting approach leads to a second-order accurate, structure-preserving scheme for the original reaction-diffusion system. Numerical experiments are presented, which demonstrate the accuracy of the proposed scheme.



قيم البحث

اقرأ أيضاً

A mass-conservative Lagrange--Galerkin scheme of second order in time for convection-diffusion problems is presented, and convergence with optimal error estimates is proved in the framework of $L^2$-theory. The introduced scheme maintains the advanta ges of the Lagrange--Galerkin method, i.e., CFL-free robustness for convection-dominated problems and a symmetric and positive coefficient matrix resulting from the discretization. In addition, the scheme conserves the mass on the discrete level. Unconditional stability and error estimates of second order in time are proved by employing two new key lemmas on the truncation error of the material derivative in conservative form and on a discrete Gronwall inequality for multistep methods. The mass-conservation property is achieved by the Jacobian multiplication technique introduced by Rui and Tabata in 2010, and the accuracy of second order in time is obtained based on the idea of the multistep Galerkin method along characteristics originally introduced by Ewing and Russel in 1981. For the first time step, the mass-conservative scheme of first order in time by Rui and Tabata in 2010 is employed, which is efficient and does not cause any loss of convergence order in the $ell^infty(L^2)$- and $ell^2(H^1_0)$-norms. For the time increment $Delta t$, the mesh size $h$ and a conforming finite element space of polynomial degree $k$, the convergence order is of $O(Delta t^2 + h^k)$ in the $ell^infty(L^2)cap ell^2(H^1_0)$-norm and of $O(Delta t^2 + h^{k+1})$ in the $ell^infty(L^2)$-norm if the duality argument can be employed. Error estimates of $O(Delta t^{3/2}+h^k)$ in discre
In the current work we build a difference analog of the Caputo fractional derivative with generalized memory kernel ($_lambda$L2-1$_sigma$ formula). The fundamental features of this difference operator are studied and on its ground some difference sc hemes generating approximations of the second order in time for the generalized time-fractional diffusion equation with variable coefficients are worked out. We have proved stability and convergence of the given schemes in the grid $L_2$ - norm with the rate equal to the order of the approximation error. The achieved results are supported by the numerical computations performed for some test problems.
112 - Dongdong Hu , Xuenian Cao 2019
In this paper, a second-order backward difference formula (abbr. BDF2) is used to approximate first-order time partial derivative, the Riesz fractional derivatives are approximated by fourth-order compact operators, a class of new alternating-directi on implicit difference scheme (abbr. ADI) is constructed for two-dimensional Riesz space fractional nonlinear reaction-diffusion equation. Stability and convergence of the numerical method are analyzed. Numerical experiments demonstrate that the proposed method is effective.
We formulate a well-posedness and approximation theory for a class of generalised saddle point problems with a specific form of constraints. In this way we develop an approach to a class of fourth order elliptic partial differential equations with po int constraints using the idea of splitting into coupled second order equations. An approach is formulated using a penalty method to impose the constraints. Our main motivation is to treat certain fourth order equations involving the biharmonic operator and point Dirichlet constraints for example arising in the modelling of biomembranes on curved and flat surfaces but the approach may be applied more generally. The theory for well-posedness and approximation is presented in an abstract setting. Several examples are described together with some numerical experiments.
74 - Jongho Park 2020
In this paper, we propose a novel overlapping domain decomposition method that can be applied to various problems in variational imaging such as total variation minimization. Most of recent domain decomposition methods for total variation minimizatio n adopt the Fenchel--Rockafellar duality, whereas the proposed method is based on the primal formulation. Thus, the proposed method can be applied not only to total variation minimization but also to those with complex dual problems such as higher order models. In the proposed method, an equivalent formulation of the model problem with parallel structure is constructed using a custom overlapping domain decomposition scheme with the notion of essential domains. As a solver for the constructed formulation, we propose a decoupled augmented Lagrangian method for untying the coupling of adjacent subdomains. Convergence analysis of the decoupled augmented Lagrangian method is provided. We present implementation details and numerical examples for various model problems including total variation minimizations and higher order models.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا