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

Finite Volume Element Methods for Two-Dimensional Time Fractional Reaction-Diffusion Equations on Triangular Grids

140   0   0.0 ( 0 )
 نشر من قبل Zhichao Fang
 تاريخ النشر 2021
  مجال البحث الهندسة المعلوماتية
والبحث باللغة English
 تأليف Zhichao Fang




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

In this paper, the time fractional reaction-diffusion equations with the Caputo fractional derivative are solved by using the classical $L1$-formula and the finite volume element (FVE) methods on triangular grids. The existence and uniqueness for the fully discrete FVE scheme are given. The stability result and optimal textit{a priori} error estimate in $L^2(Omega)$-norm are derived, but it is difficult to obtain the corresponding results in $H^1(Omega)$-norm, so another analysis technique is introduced and used to achieve our goal. Finally, two numerical examples in different spatial dimensions are given to verify the feasibility and effectiveness.



قيم البحث

اقرأ أيضاً

In this paper, we present an inverse problem of identifying the reaction coefficient for time fractional diffusion equations in two dimensional spaces by using boundary Neumann data. It is proved that the forward operator is continuous with respect t o the unknown parameter. Because the inverse problem is often ill-posed, regularization strategies are imposed on the least fit-to-data functional to overcome the stability issue. There may exist various kinds of functions to reconstruct. It is crucial to choose a suitable regularization method. We present a multi-parameter regularization $L^{2}+BV$ method for the inverse problem. This can extend the applicability for reconstructing the unknown functions. Rigorous analysis is carried out for the inverse problem. In particular, we analyze the existence and stability of regularized variational problem and the convergence. To reduce the dimension in the inversion for numerical simulation, the unknown coefficient is represented by a suitable set of basis functions based on a priori information. A few numerical examples are presented for the inverse problem in time fractional diffusion equations to confirm the theoretic analysis and the efficacy of the different regularization methods.
289 - S. Singh , S. Sircar 2019
We provide a preliminary comparison of the dispersion properties, specifically the time-amplification factor, the scaled group velocity and the error in the phase speed of four spatiotemporal discretization schemes utilized for solving the one-dimens ional (1D) linear advection diffusion reaction (ADR) equation: (a) An explicit (RK2) temporal integration combined with the Optimal Upwind Compact Scheme (or OUCS3) and the central difference scheme (CD2) for second order spatial discretization, (b) a fully implicit mid-point rule for time integration coupled with the OUCS3 and the Leles compact scheme for first and second order spatial discretization, respectively, (c) An implicit (mid-point rule)-explicit (RK2) or IMEX time integration blended with OUCS3 and Leles compact scheme (where the IMEX time integration follows the same ideology as introduced by Ascher et al.), and (d) the IMEX (mid-point/RK2) time integration melded with the New Combined Compact Difference scheme (or NCCD scheme). Analysis reveal the superior resolution features of the IMEX-NCCD scheme including an enhanced region of neutral stability (a region where the amplification factor is close to one), a diminished region of spurious propagation characteristics (or a region of negative group velocity) and a smaller region of nonzero phase speed error. The dispersion error of these numerical schemes through the role of q-waves is further investigated using the novel error propagation equation for the 1D linear ADR equation. Again, the in silico experiments divulge excellent Dispersion Relation Preservation (DRP) properties of the IMEX-NCCD scheme including minimal dissipation via implicit filtering and negligible unphysical oscillations (or Gibbs phenomena) on coarser grids.
In this paper, we introduce and analyse a surface finite element discretization of advection-diffusion equations with uncertain coefficients on evolving hypersurfaces. After stating unique solvability of the resulting semi-discrete problem, we prove optimal error bounds for the semi-discrete solution and Monte Carlo samplings of its expectation in appropriate Bochner spaces. Our theoretical findings are illustrated by numerical experiments in two and three space dimensions.
133 - Erik Burman , Cuiyu He 2018
We consider primal-dual mixed finite element methods for the advection--diffusion equation. For the primal variable we use standard continuous finite element space and for the flux we use the Raviart-Thomas space. We prove optimal a priori error esti mates in the energy- and the $L^2$-norms for the primal variable in the low Peclet regime. In the high Peclet regime we also prove optimal error estimates for the primal variable in the $H(div)$ norm for smooth solutions. Numerically we observe that the method eliminates the spurious oscillations close to interior layers that pollute the solution of the standard Galerkin method when the local Peclet number is high. This method, however, does produce spurious solutions when outflow boundary layer presents. In the last section we propose two simple strategies to remove such numerical artefacts caused by the outflow boundary layer and validate them numerically.
In this paper, we extend the class of kernel methods, the so-called diffusion maps (DM) and ghost point diffusion maps (GPDM), to solve the time-dependent advection-diffusion PDE on unknown smooth manifolds without and with boundaries. The core idea is to directly approximate the spatial components of the differential operator on the manifold with a local integral operator and combine it with the standard implicit time difference scheme. When the manifold has a boundary, a simplified version of the GPDM approach is used to overcome the bias of the integral approximation near the boundary. The Monte-Carlo discretization of the integral operator over the point cloud data gives rise to a mesh-free formulation that is natural for randomly distributed points, even when the manifold is embedded in high-dimensional ambient space. Here, we establish the convergence of the proposed solver on appropriate topologies, depending on the distribution of point cloud data and boundary type. We provide numerical results to validate the convergence results on various examples that involve simple geometry and an unknown manifold. Additionally, we also found positive results in solving the one-dimensional viscous Burgers equation where GPDM is adopted with a pseudo-spectral Galerkin framework to approximate nonlinear advection term.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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