On the implementation of linear finite element method for nonlocal diffusion model over 2D domain

Zuoqiang Shi

2024-08-28

Abstract:We propose an implementation of linear finite element method for nonlocal diffusion problem in 2D space. In the implementation, we reduce the integral from 4D to 2D which would simplify the computation significantly.

Numerical Analysis

What problem does this paper attempt to address?

A General Finite Element Method: Extension of Variational Analysis for Nonlinear Heat Conduction with Temperature-Dependent Properties and Boundary Conditions, and Its Implementation As Local Refinement.

Xin Yao,Yihe Wang,Jianxing Leng

DOI: https://doi.org/10.1016/j.camwa.2021.08.024

IF: 3.218

2021-01-01

Computers & Mathematics with Applications

Abstract:In simulation of heat conduction with temperature-independent physical properties and boundary conditions (BCs), Galerkin residual analysis and variational analysis yield equivalent finite element method (FEM), the conventional FEM. However, if the properties and BCs are temperature-dependent, it is discovered that their derivatives further induce nonlinearity of FEM which consequently generates divergence between the two analyses. A general FEM, extension of variational analysis, is derived as general form of conventional FEM modeling nonlinear heat conduction. Numerical examples demonstrate that the general FEM produces results with considerably higher accuracy and stability and also possesses higher performances on conforming with both two analyses. Since general FEM degenerates to conventional FEM if derivatives are of small-amplitude or zero and its direct implementation to the entire domain is costive, general FEM is alternatively utilized as local refinement of governing equation only to points with significant derivatives. The strategy of local refinement is optimized to enhance efficiency.
Finite Element Solution in Fluid Mechanics Using Exponential Function Based Interpolation

Shaowen Fang,Xingfei Yuan,Ruojun Qian,Shuhui Jiang

DOI: https://doi.org/10.4028/www.scientific.net/amm.580-583.3051

2014-01-01

Applied Mechanics and Materials

Abstract:The wave problem of finite element solution of one dimensional stationary convection diffusion equation is analyzed. The reason for wave is the poor continuity of linear Lagrange shape function. By using exponential function based interpolation (EFBI), the results of finite element solution are compared with that of the analytical solution. The results indicate that EFBI is effective to deal with the problem of numerical wave. The shape function of three dimensional EFBI is derived and analyzed. Morphological analysis shows that EFBI has good differentiability and adaptability.
Numerical Solution of the Nonlocal Diffusion Equation on the Real Line.

Chunxiong Zheng,Jiashun Hu,Qiang Du,Jiwei Zhang

DOI: https://doi.org/10.1137/16m1090107

IF: 2.968

2017-01-01

SIAM Journal on Scientific Computing

Abstract:Numerical computation of a nonlocal diffusion equation on the real axis is considered in this paper. We first apply an extensively studied quadrature scheme to obtain a discrete nonlocal diffusion system on an unbounded domain. Then we derive an alternative formulation of the discrete problem based on the spectral analysis of the z-transform. This new formulation can be seen as a system defined on a bounded domain with an artificial boundary condition, and it allows us to reformulate the original infinite domain problem into an equivalent bounded domain problem. To numerically implement the exact artificial boundary condition, we apply the trapezoidal quadrature rule to approximate the contour integral induced by the inverse z-transform. Numerical examples are presented to demonstrate the effectiveness of our approach.
High-Order Numerical Method for 1D Non-local Diffusive Equation

D. Do,H. Nick Zinat Matin,M. L. Delle Monache

2023-11-11

Abstract:In this paper we present a non-local numerical scheme based on the Local Discontinuous Galerkin method for a non-local diffusive partial differential equation with application to traffic flow. In this model, the velocity is determined by both the average of the traffic density as well as the changes in the traffic density at a neighborhood of each point. We discuss nonphysical behaviors that can arise when including diffusion, and our measures to prevent them in our model. The numerical results suggest that this is an accurate method for solving this type of equation and that the model can capture desired traffic flow behavior. We show that computation of the non-local convolution results in $\mathcal{O}(n^2)$ complexity, but the increased computation time can be mitigated with high-order schemes like the one proposed.

Numerical Analysis
A Finite Difference Method for a Nonlinear Diffusive Model

WANG Hai-ming

DOI: https://doi.org/10.3969/j.issn.1672-4143.2007.02.001

2007-01-01

Abstract:A finite difference method is proposed for a nonlinear diffusive model.The scheme resulting from this method is linear in the sense that it requires only to solve a single linear algebraic system at each time step.Moreover,the non-physical numerical solution does not occur,and the numerical solution converges to zero or some constant when time tends to infinity,depending on the initial value and the parameter in the model.This implies that the numerical solution possesses the same qualitative properties as that of the analytical solution.Numerical results are given to corroborate the theoretical results.
A spectral collocation method for nonlocal diffusion equations

Hao Tian,Jing Zhang

DOI: https://doi.org/10.48550/arXiv.1802.10545

2018-03-01

Abstract:Nonlocal diffusion model provides an appropriate description of the diffusion process of solute in the complex medium, which cannot be described properly by classical theory of PDE. However, the operators in the nonlocal diffusion models are nonlocal, so the resulting numerical methods generate dense or full stiffness matrices. This imposes significant computational and memory challenge for a nonlocal diffusion model. In this paper, we develop a spectral collocation method for the nonlocal diffusion model and provide a rigorous error analysis which theoretically justifies the spectral rate of convergence provided that the kernel functions and the source functions are sufficiently smooth. Compared to finite difference methods and finite element methods, because of the high order convergence rates, the numerical cost of spectral collocation methods will be greatly decreased. Numerical results confirm the exponential rate of convergence.

Numerical Analysis
Error analysis of finite element method for nonlocal diffusion model

Zuoqiang Shi

2024-09-03

Abstract:We analyze the error of finite element method for nonlocal diffusion model include both conformal and nonconformal method. We also consider the mesh with and without shape regularity. For shape regular mesh, finite element method for nonlocal diffusion model is asymptotic preserving and the error is $O(h^k+\delta)$. For shape irregular mesh, the error becomes $O(\frac{h^{k+1}}{\delta}+\delta)$.

Numerical Analysis
A local character based method for solving linear systems of radiation diffusion problems

Shuai Ye,Hengbin An,Xinhai Xu

DOI: https://doi.org/10.1016/j.jcp.2019.109218

2020-01-05

Abstract:The radiation diffusion problem is a kind of {time-dependent} nonlinear equations. For solving the radiation diffusion equations, many linear systems are obtained in the nonlinear iterations at each time step. The cost of linear equations dominates the numerical simulation of radiation diffusion applications, such as inertial confinement fusion, etc. {Usually,} iterative methods are used to solve the linear systems in {a} real application. {Moreover, the solution of the previous nonlinear iteration or the solution of the previous time step is typically used as the initial guess for solving the current linear equations.} Because of the strong local character in ICF, with the advancing of nonlinear iteration and time step, the solution of the linear system changes dramatically in some local domain, and changes mildly or even has no change in the rest domain. In this paper, a local {character-based} method is proposed to solve the linear systems of radiation diffusion problems. The proposed method consists of three steps: firstly, a local domain (algebraic domain) is constructed; secondly, the subsystem on the local domain is solved; and lastly, the whole system will be solved. Two methods are given to construct the local domain. One is based on the spatial gradient, and the other is based on the residual. Numerical tests for a two-dimensional heat conduction model problem, and two real application models, the multi-group radiation diffusion equations and the three temperature energy equations, are conducted. The test results show that the solution time for solving the linear system can be reduced dramatically by using the local character-based method.

Numerical Analysis
Finite Difference/Collocation Method for Two-Dimensional Sub-Diffusion Equation with Generalized Time Fractional Derivative

Qinwu Xu,Zhoushun Zheng

DOI: https://doi.org/10.4208/jms.v47n2.14.03

2014-01-01

Journal of Mathematical Study

Abstract:In this paper, we propose a finite difference/collocation method for two-dimensional time fractional diffusion equation with generalized fractional operator. The main purpose of this paper is to design a high order numerical scheme for the new generalized time fractional diffusion equation. First, a finite difference approximation formula is derived for the generalized time fractional derivative, which is verified with order 2 - alpha (0 < alpha < 1). Then, collocation method is introduced for the two-dimensional space approximation. Unconditional stability of the scheme is proved. To make the method more efficient, the alternating direction implicit method is introduced to reduce the computational cost. At last, numerical experiments are carried out to verify the effectiveness of the scheme.
Analysis and efficient implementation of alternating direction implicit finite volume method for Riesz space‐fractional diffusion equations in two space dimensions

Huan Liu,Xiangcheng Zheng,Hongfei Fu,Hong Wang

DOI: https://doi.org/10.1002/num.22554

2020-10-02

Numerical Methods for Partial Differential Equations

Abstract:In this article, we develop a Crank–Nicolson alternating direction implicit finite volume method for time‐dependent Riesz space‐fractional diffusion equation in two space dimensions. Norm‐based stability and convergence analysis are given to show that the developed method is unconditionally stable and of second‐order accuracy both in space and time. Furthermore, we develop a lossless matrix‐free fast conjugate gradient method for the implementation of the numerical scheme, which only has memory requirement and computational complexity per iteration with N being the total number of spatial unknowns. Several numerical experiments are presented to demonstrate the effectiveness and efficiency of the proposed scheme for large‐scale modeling and simulations.

mathematics, applied
High-order finite element methods for nonlinear convection-diffusion equation on time-varying domain

Chuwen Ma,Weiying Zheng

DOI: https://doi.org/10.48550/arXiv.2106.16026

2021-12-30

Abstract:A high-order finite element method is proposed to solve the nonlinear convection-diffusion equation on a time-varying domain whose boundary is implicitly driven by the solution of the equation. The method is semi-implicit in the sense that the boundary is traced explicitly with a high-order surface-tracking algorithm, while the convection-diffusion equation is solved implicitly with high-order backward differentiation formulas and fictitious-domain finite element methods. By two numerical experiments for severely deforming domains, we show that optimal convergence orders are obtained in energy norm for third-order and fourth-order methods.

Numerical Analysis
High accuracy error estimates of a Galerkin finite element method for nonlinear time fractional diffusion equation

Jincheng Ren,Dongyang Shi,Seakweng Vong

DOI: https://doi.org/10.1002/num.22428

2019-09-04

Numerical Methods for Partial Differential Equations

Abstract:In this work, an effective and fast finite element numerical method with high‐order accuracy is discussed for solving a nonlinear time fractional diffusion equation. A two‐level linearized finite element scheme is constructed and a temporal–spatial error splitting argument is established to split the error into two parts, that is, the temporal error and the spatial error. Based on the regularity of the time discrete system, the temporal error estimate is derived. Using the property of the Ritz projection operator, the spatial error is deduced. Unconditional superclose result in H1‐norm is obtained, with no additional regularity assumption about the exact solution of the problem considered. Then the global superconvergence error estimate is obtained through the interpolated postprocessing technique. In order to reduce storage and computation time, a fast finite element method evaluation scheme for solving the nonlinear time fractional diffusion equation is developed. To confirm the theoretical error analysis, some numerical results are provided.

mathematics, applied
Efficient L1-ADI finite difference method for the two-dimensional nonlinear time-fractional diffusion equation

Yubing Jiang,Hu Chen,Tao Sun,Chaobao Huang

DOI: https://doi.org/10.1016/j.amc.2024.128609

IF: 4.397

2024-02-22

Applied Mathematics and Computation

Abstract:In this work, we propose an efficient finite difference method for the two-dimensional nonlinear time-fractional diffusion equation with weakly singular solutions. By using backward formula for the approximation of nonlinear term, and L1 scheme on uniform mesh for discretisation of temporal Caputo fractional derivative, a linear scheme is constructed and analysed. Stability and pointwise-in-time convergence of the fully discrete scheme are rigorously established. Numerical results are provided to confirm the sharpness of theoretical analysis.

mathematics, applied
Local discontinuous Galerkin methods with explicit-implicit-null time discretizations for solving nonlinear diffusion problems

Haijin Wang,Qiang Zhang,Shiping Wang,Chi-Wang Shu

DOI: https://doi.org/10.48550/arXiv.1903.11859

2019-03-28

Abstract:In this paper we discuss the local discontinuous Galerkin methods coupled with two specific explicit-implicit-null time discretizations for solving one-dimensional nonlinear diffusion problems $U_t=(a(U)U_x)_x$. The basic idea is to add and subtract two equal terms $a_0 U_{xx}$ on the right hand side of the partial differential equation, then to treat the term $a_0 U_{xx}$ implicitly and the other terms $(a(U)U_x)_x-a_0 U_{xx}$ explicitly. We give stability analysis for the method on a simplified model by the aid of energy analysis, which gives a guidance for the choice of $a_0$, i.e, $a_0 \ge \max\{a(u)\}/2$ to ensure the unconditional stability of the first order and second order schemes. The optimal error estimate is also derived for the simplified model, and numerical experiments are given to demonstrate the stability, accuracy and performance of the schemes for nonlinear diffusion equations.

Numerical Analysis
Local modification of subdiffusion by initial Fickian diffusion: Multiscale modeling, analysis and computation

Xiangcheng Zheng,Yiqun Li,Wenlin Qiu

2024-01-30

Abstract:We propose a local modification of the standard subdiffusion model by introducing the initial Fickian diffusion, which results in a multiscale diffusion model. The developed model resolves the incompatibility between the nonlocal operators in subdiffusion and the local initial conditions and thus eliminates the initial singularity of the solutions of the subdiffusion, while retaining its heavy tail behavior away from the initial time. The well-posedness of the model and high-order regularity estimates of its solutions are analyzed by resolvent estimates, based on which the numerical discretization and analysis are performed. Numerical experiments are carried out to substantiate the theoretical findings.

Numerical Analysis
Numerical analysis of a class of penalty discontinuous Galerkin methods for nonlocal diffusion problems

Qiang Du,Lili Ju,Jianfang Lu,Xiaochuan Tian

DOI: https://doi.org/10.1051/m2an/2024064

2024-08-11

ESAIM: Mathematical Modelling and Numerical Analysis

Abstract:In this paper, we consider a class of discontinuous Galerkin (DG) methods for one-dimensional nonlocal diffusion (ND) problems. The nonlocal models are integral equations that are widely used in describing many physical phenomena with long-range interactions. The ND problem is the nonlocal analog of the classic diffusion problem, and the nonlocality disappears as the interaction radius (horizon) vanishes, then the ND problem recovers the classic diffusion problem. Under certain conditions, the exact solution to the ND problem may contain discontinuities, and this is one of the properties of the ND problem that sets it apart from the classic diffusion problem. It is natural to adopt the DG method to compute the ND problems since the DG method shows its great advantages in resolving problems with discontinuities in the computational fluid dynamics in the past several decades. Based on \cite{DuCAMC2020}, we construct the DG methods with different penalty terms, and the proposed DG methods have their local counterparts as the horizon vanishes. This indicates the proposed methods will converge to the existing DG schemes as the horizon vanishes, which is crucial to obtain the {\it asymptotic compatibility}. Rigorous proofs are provided to demonstrate the stability, error estimates, and asymptotic compatibility of the proposed DG schemes. To see the nonlocal diffusion effect, we also consider the time-dependent convection-diffusion problems with nonlocal diffusion. We conduct several numerical experiments, including accuracy tests and Burgers' equation with nonlocal diffusion, and various horizons are taken to show the good performance of the proposed algorithm and validate the theoretical findings.

mathematics, applied
Numerical study of a domain decomposition method for a two-scale linear transport equation

Xu Yang,François Golse,Zhongyi Huang,Shi Jin

DOI: https://doi.org/10.3934/nhm.2006.1.143

2006-01-01

Abstract:We perform a numerical study on a domain decomposition method proposed in [13] for the linear transport equation between a diffusive and a non-diffusive region. This method avoids iterating the diffusion and transport solutions as in a typical domain decomposition method. Our numerical results, in both one and two space dimensions, confirm the theoretical analysis of [13]. We also provide an improved second order method that provides a more accurate numerical solution than that proposed in [13].
A conservative Eulerian finite element method for transport and diffusion in moving domains

Maxim Olshanskii,Henry von Wahl

2024-04-11

Abstract:The paper introduces a finite element method for an Eulerian formulation of partial differential equations governing the transport and diffusion of a scalar quantity in a time-dependent domain. The method follows the idea from Lehrenfeld & Olshanskii [ESAIM: M2AN, 53(2): 585-614, 2019] of a solution extension to realise the Eulearian time-stepping scheme. However, a reformulation of the partial differential equation is suggested to derive a scheme which conserves the quantity under consideration exactly on the discrete level. For the spatial discretisation, the paper considers an unfitted finite element method. Ghost-penalty stabilisation is used to release the discrete solution extension and gives a scheme robust against arbitrary intersections between the mesh and geometry interface. The stability is analysed for both first- and second-order backward differentiation formula versions of the scheme. Several numerical examples in two and three spatial dimensions are included to illustrate the potential of this method.

Numerical Analysis
Stability and Error Analysis for a Second-Order Fast Approximation of the Local and Nonlocal Diffusion Equations on the Real Line

Chunxiong Zheng,Qiang Du,Xiang Ma,Jiwei Zhang

DOI: https://doi.org/10.1137/19m1285822

IF: 3.039

2020-01-01

SIAM Journal on Numerical Analysis

Abstract:The stability and error analysis of a second-order fast approximation are considered for the one-dimensional local and nonlocal diffusion equations in the unbounded spatial domain. We first use the conventional central difference scheme to discretize the local second-order spatial derivative operator and use an asymptotically compatible difference scheme to discretize the spatial nonlocal diffusion operator, and apply second-order backward differentiation formula (BDF2) to approximate the temporal derivative to achieve a fully discrete infinity system. To solve the resulting fully discrete systems, we develop a unified framework that is applicable to the discretization of both local and nonlocal problems. A key ingredient is to derive Dirichlet-to-Neumann (DtN)-type absorbing boundary conditions (ABCs). To do so, we apply the z-transform and solve an exterior problem using an iteration technique to derive a Dirichlet-to-Dirichlet (DtD)-type mapping as exact ABCs. After that, we use the Green formula to reformulate the DtD-type mapping equivalently as the DtN-type mapping. The resulting DtN-type mapping allows us to reduce the infinity discrete system into a finite discrete system in a truncated computational domain of interest, and also make it possible to present the stability and convergence analysis of the reduced problem under some open but reasonable assumptions. To efficiently implement the exact ABCs, we further develop a fast convolution algorithm based on approximation of the contour integral induced by the inverse z-transform. The stability and error analysis of the reduced finite discrete system based on the fast algorithm for exact ABCs are also established, and numerical examples are provided to demonstrate the effectiveness of our proposed approach.

On the implementation of linear finite element method for nonlocal diffusion model over 2D domain

A General Finite Element Method: Extension of Variational Analysis for Nonlinear Heat Conduction with Temperature-Dependent Properties and Boundary Conditions, and Its Implementation As Local Refinement.

Finite Element Solution in Fluid Mechanics Using Exponential Function Based Interpolation

Numerical Solution of the Nonlocal Diffusion Equation on the Real Line.

High-Order Numerical Method for 1D Non-local Diffusive Equation

A Finite Difference Method for a Nonlinear Diffusive Model

A spectral collocation method for nonlocal diffusion equations

Error analysis of finite element method for nonlocal diffusion model

A local character based method for solving linear systems of radiation diffusion problems

Finite Difference/Collocation Method for Two-Dimensional Sub-Diffusion Equation with Generalized Time Fractional Derivative

Analysis and efficient implementation of alternating direction implicit finite volume method for Riesz space‐fractional diffusion equations in two space dimensions

High-order finite element methods for nonlinear convection-diffusion equation on time-varying domain

High accuracy error estimates of a Galerkin finite element method for nonlinear time fractional diffusion equation

Efficient L1-ADI finite difference method for the two-dimensional nonlinear time-fractional diffusion equation

Local discontinuous Galerkin methods with explicit-implicit-null time discretizations for solving nonlinear diffusion problems

Local modification of subdiffusion by initial Fickian diffusion: Multiscale modeling, analysis and computation

Numerical analysis of a class of penalty discontinuous Galerkin methods for nonlocal diffusion problems

Numerical study of a domain decomposition method for a two-scale linear transport equation

A conservative Eulerian finite element method for transport and diffusion in moving domains

Stability and Error Analysis for a Second-Order Fast Approximation of the Local and Nonlocal Diffusion Equations on the Real Line