Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userEvolving surface finite element methods for random advection-diffusion equations
70
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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-dimensional (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 define new unfitted finite element methods for numerically approximating the solution of surface partial differential equations using bulk finite elements. The key idea is that the $n$-dimensional hypersurface, $Gamma subset mathbb{R}^{n+1}$, is embedded in a polyhedral domain in $mathbb R^{n+1}$ consisting of a union, $mathcal{T}_h$, of $(n+1)$-simplices. The finite element approximating space is based on continuous piece-wise linear finite element functions on $mathcal{T}_h$. Our first method is a sharp interface method, emph{SIF}, which uses the bulk finite element space in an approximating weak formulation obtained from integration on a polygonal approximation, $Gamma_{h}$, of $Gamma$. The full gradient is used rather than the projected tangential gradient and it is this which distinguishes emph{SIF} from the method of [42]. The second method, emph{NBM}, is a narrow band method in which the region of integration is a narrow band of width $O(h)$. emph{NBM} is similar to the method of [13]. but again the full gradient is used in the discrete weak formulation. The a priori error analysis in this paper shows that the methods are of optimal order in the surface $L^{2}$ and $H^{1}$ norms and have the advantage that the normal derivative of the discrete solution is small and converges to zero. Our third method combines bulk finite elements, discrete sharp interfaces and narrow bands in order to give an unfitted finite element method for parabolic equations on evolving surfaces. We show that our method is conservative so that it preserves mass in the case of an advection diffusion conservation law. Numerical results are given which illustrate the rates of convergence.
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.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
