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

Primal dual mixed finite element methods for indefinite advection--diffusion equations

134   0   0.0 ( 0 )
 نشر من قبل Cuiyu He
 تاريخ النشر 2018
  مجال البحث
والبحث باللغة English




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

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 estimates 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 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.
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 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.
139 - Zhichao Fang 2021
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.
Under some regularity assumptions, we report an a priori error analysis of a dG scheme for the Poisson and Stokes flow problem in their dual mixed formulation. Both formulations satisfy a Babuv{s}ka-Brezzi type condition within the space H(div) x L2. It is well known that the lowest order Crouzeix-Raviart element paired with piecewise constants satisfies such a condition on (broken) H1 x L2 spaces. In the present article, we use this pair. The continuity of the normal component is weakly imposed by penalizing jumps of the broken H(div) component. For the resulting methods, we prove well-posedness and convergence with constants independent of data and mesh size. We report error estimates in the methods natural norms and optimal local error estimates for the divergence error. In fact, our finite element solution shares for each triangle one DOF with the CR interpolant and the divergence is locally the best-approximation for any regularity. Numerical experiments support the findings and suggest that the other errors converge optimally even for the lowest regularity solutions and a crack-problem, as long as the crack is resolved by the mesh.
التعليقات
جاري جلب التعليقات جاري جلب التعليقات
سجل دخول لتتمكن من متابعة معايير البحث التي قمت باختيارها
mircosoft-partner

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