Two-dimensional Diffusion Equation Research Articles

L1-FEM discretizations for two-dimensional multiterm fractional delay diffusion equations

A two-dimensional multiterm fractional delay diffusion equation is considered. The representation of the exact solution of the equation is derived and it is shown that the solution exhibits singular behaviors at multiple nodes due to the initial singularity and time delay. This results in the numerical schemes for solving the equation typically have a lower order of convergence in time. The problem is approximated in time by the L1 and Alikhanov schemes on symmetrical graded meshes, while in space the standard finite element method is applied. Numerical stability and convergence are presented for the schemes. Numerical experiments are performed to show the effectiveness of the schemes.

Numerical analysis of high vacuum packaging efficiency for narrow structured devices

High vacuum packaging is the key manufacturing technology of passive vacuum devices. The inherent characteristics of narrow and long space and large surface-to-volume ratio make it challenging to effectively remove the desorption gas inside the device. To this end, the pressure distribution of two standard packaging methods during the pump-down process is investigated using the two-dimensional diffusion equation. Experiments and simulations show that a pressure gradient persists within the slit even when pumping reaches the quasi-equilibrium state. In addition, a Monte Carlo method based on the cosine law reveals that the high-frequency collisions of gas molecules with the walls impede the pumping process. Furthermore, the results indicate that for tube pumping, the number of collisions of gas molecules with the wall during transport within the gap is inversely proportional to the gap height. For edge pumping, the number of collisions is inversely proportional to the square of the gap height. This provides a theoretical basis for efficient high vacuum packaging and long-term pressure maintenance for passive vacuum devices.

Iterated fractional Tikhonov method for recovering the source term and initial data simultaneously in a two-dimensional diffusion equation

In this paper, an inverse problem of a two-dimensional diffusion equation is considered. The purpose here includes recovering not only the source term but also the initial value from given observations at two fixed times t=T1 and t=T2. The uniqueness of the solutions for the inverse problem is given. Instead of the classical Tikhonov regularization strategy, we propose an iterated fractional Tikhonov regularization method to solve this problem, an a priori and an a posterior parameter selection rules and corresponding convergence rates are derived. For verification of the theoretical estimates, several numerical examples are constructed and compared with the standard iterated Tikhonov regularization method.

Stiff-cut leap-frog scheme for fractional Laplacian diffusion equations

In this study, we investigate the numerical discretization of two-dimensional fractional Laplacian diffusion equations. After discretizing the space derivative, we split the discrete fractional Laplacian into two parts. One is the τ preconditioner multiplying a constant, solved using the leap-frog scheme, while the other part and the remaining nonlinear term are treated with an explicit scheme. This results in what is termed the stiff-cut leap-frog scheme, which is more efficient than the traditional leap-frog scheme. The main advantage of the proposed scheme is that the coefficient matrices of the generated linear systems belong to τ algebra, which can be diagonalized by fast discrete sine transform with the complexity of O(NlogN), where N is the freedom of unknowns. Hence, no iterative method is required at each time integration. Moreover, we theoretically prove that the proposed scheme is unconditionally stable and converges in second order. Ultimately, numerical experiments are conducted to demonstrate the precision and effectiveness of the stiff-cut leap-frog scheme.

Two-dimensional counter-current capillary imbibition of a wetting phase into a partially submerged porous cylindrical matrix block

The purpose of this study is to address the two-dimensional counter-current capillary dominant imbibition of a wetting phase into a water-wet porous cylindrical matrix block partially submerged in the wetting phase. A two-dimensional unsteady-state diffusion equation is used to model the process. The governing equation is solved using a combination of the Laplace and the finite Fourier sine transforms to find and analyze the solutions for the normalized water saturation and the volume of the imbibed wetting phase. The results reveal that the volume of the imbibed wetting phase and the capillary diffusion shape factor for a partially submerged matrix block are significantly lower compared to those of a fully submerged matrix block, highlighting the overestimation of imbibed volume using available models based on full immersion in the wetting phase. It has been observed that the volume of the imbibed wetting phase increases over time until reaching a state of equilibrium. In the case of a partially submerged matrix block, the shape factor is inversely proportional to the square root of time (σ ∼ 1/t) during the early time and decreases sharply as the imbibed wetting phase reaches an equilibrium. In the case of a fully submerged matrix block, the shape factor is inversely proportional to the square root of time (σ ∼ 1/t) during the early time and later reaches a pseudo-steady-state value. The proposed model, along with the findings obtained, advances our understanding of capillary imbibition in porous media.

A general fourth-order mesoscopic multiple-relaxation-time lattice Boltzmann model and its macroscopic finite-difference scheme for two-dimensional diffusion equations

In this work, we first develop a general mesoscopic multiple-relaxation-time lattice Boltzmann (MRT-LB) model for the two-dimensional diffusion equation with the constant diffusion coefficient and source term, where the D2Q5 (five discrete velocities in two-dimensional space) lattice structure is considered. Then we exactly derive the corresponding macroscopic finite-difference scheme of the MRT-LB model. Additionally, we also propose a proper MRT-LB model for the diffusion equation with a linear source term, and obtain a macroscopic six-level finite-difference scheme. After that, we conduct the accuracy and stability analysis on the finite-difference scheme and the mesoscopic MRT-LB model, and find that at the diffusive scaling, both of them can achieve a fourth-order accuracy in space based on the Taylor expansion. The stability analysis also shows that they are both unconditionally stable. Finally, some numerical experiments are conducted, and the numerical results are also consistent with our theoretical analysis.

Application of data-driven model reduction techniques in reactor neutron field calculations

High-order harmonic techniques can be used to recreate neutron flux distributions in reactor cores using the neutron diffusion equation. However, traditional source iteration and source correction iteration techniques have sluggish convergence rates and protracted calculation periods. The correctness of the implicitly restarted Arnoldi method (IRAM) in resolving the eigenvalue problems of the one-dimensional and two-dimensional neutron diffusion equations was confirmed by computing the benchmark problems SLAB_1D_1G and two-dimensional steady-state TWIGL using IRAM. By integrating Galerkin projection with Proper Orthogonal Decomposition (POD) techniques, a POD-Galerkin reduced-order model was developed and the IRAM model was used as the full-order model. For 14 macroscopic cross-section values, the TWIGL benchmark problem was perturbed within a 20% range. We extracted 100 sample points using the Latin hypercube sampling method, and 70% of the samples were used as the testing set to assess the performance of the reduced-order model The remaining 30% were utilized as the training set to develop the reduced-order model, which was employed to rebuild the TWIGL benchmark problem. The reduced-order model demonstrates good flexibility and can efficiently and accurately forecast the effective multiplication factor and neutron flux distribution in the core. The reduced-order model predicts keff and neutron flux distribution with a high degree of agreement compared to the full-order model. Additionally, the reduced-order model's computation time is only 10.18% of that required by the full-order model.The neutron flux distribution of the steady-state TWIGL benchmark was recreated using the reduced-order model. The obtained results indicate that the reduced-order model can accurately predict the keff and neutron flux distribution of the steady-state TWIGL benchmark.Overall, the proposed technique not only has the potential to accurately project neutron flux distributions in transient settings, but is also relevant for reconstructing neutron flux distributions in steady-state conditions; thus, its applicability is bound to increase in the future.

Heavy metal transport with adsorption for instantaneous and exponential attenuation of concentration

Heavy metal pollutant is a serious problem in environmental pollution, and it is very difficult to eradicate once it enters the soil. As heavy metal adsorption has been proven to occur, the heavy metal’s behaviour can be modeled as a transport equation with adsorption. Previous adsorption term mostly due to the concentration alone, while in here, the desorption effect given by the rate of change of the concentration is also included. Also, the heavy metals are frequently considered to enter the soil after being dumped into the soil for a certain period of time. But, quick dumping onto the soil can introduce heavy metal instantaneously. Heavy metals entering the soil through leaching or when their concentration in the soil is influenced by chemical reactions, can all lead to the exponential decay of heavy metals entering the soil. Based on two-dimensional advection diffusion equation (ADE) with the new adsorption term, analytical solutions are obtained for the cases of instantaneous and exponential attenuation of heavy metals emission to soil by the method of Laplace transform. The results highlight the significant influence of emission type on the peak concentrations. If heavy metals are instantaneously enter the soil, the peak occurs in the range of 1–3 m radius from the point of emission on the first day, while for exponential attenuation the peak occurs close to the point of emission. Furthermore, there exists a correlation between retardation factors and heavy metal concentrations, where a decrease in retardation factors leads to an increase in heavy metal concentration. It is essential to investigate both types of heavy metals emission to provide valuable information for proper pollution management, effective environmental regulations and enforcement.

On numerical solution of two-dimensional variable-order fractional diffusion equation arising in transport phenomena

<abstract><p>In recent years, the application of variable-order (VO) fractional differential equations for describing complex physical phenomena ranging from biology, hydrology, mechanics and viscoelasticity to fluid dynamics has become one of the most hot topics in the context of scientific modeling. An interesting aspect of VO operators is their capability to address the behavior of scientific and engineering systems with time and spatially varying properties. The VO fractional diffusion equation is a fundamental model that allows transitions among sub-diffusive, diffusive and super-diffusive behaviors without altering the underlying governing equations. In this paper, we considered the two-dimensional fractional diffusion equation with the Caputo time VO derivative, which is essential for describing anomalous diffusion in real-world complex systems. A new Crank-Nicolson (C-N) difference scheme and an efficient explicit decoupled group (EDG) method were proposed to solve the problem under consideration. The proposed EDG method is based on a skewed difference scheme in conjunction with a grouping procedure of the solution grid points. Special attention was devoted to investigating the stability and convergence of the proposed methods. Three numerical examples with known exact analytical solutions were provided to illustrate our considerations. The proposed methods were shown to be stable and convergent theoretically as well as numerically. In addition, a comparative study was done between the EDG method and the C-N difference scheme. It was found that the proposed methods are accurate in simulating the considered problem, while the EDG method is superior to the C-N difference method in terms of Central Processing Unit (CPU) timing, verifying the efficiency of the former method in solving the VO problem.</p></abstract>

LAPLACE DECOMPOSITION METHOD FOR SOLVING THE TWO-DIMENSIONAL DIFFUSION PROBLEM IN FRACTAL HEAT TRANSFER

In this paper, the Local Fractional Laplace Decomposition Method (LFLDM) is used for solving a type of Two-Dimensional Fractional Diffusion Equation (TDFDE). In this method, first we apply the Laplace transform and its inverse to the main equation, and then the Adomian decomposition is used to obtain approximate/analytical solution. The accuracy and applicability of the LFLDM is shown through two examples. The LFLDM results are in good agreement with the exact solution of the problems.

A second-order difference scheme for two-dimensional two-sided space distributed-order fractional diffusion equations with variable coefficients

A second-order difference scheme for two-dimensional two-sided space distributed-order fractional diffusion equations with variable coefficients

Moving iso-contour method for solving partial differential equations

The numerical solution of partial differential equations is often performed on a numerical grid, where the grid points are used for estimating the partial derivatives. The grid can be fully static as in Eulerian type of solution method, or the grid points can move during the solution, which is the case in Lagrangian type of method. In the current article, a numerical solution method is presented, where the grid points are located on iso-contours of the two-dimensional field. The method calculates the local movement of the iso-contours according to an evolution equation described by the PDE, and the solution proceeds by moving the grid points towards the calculated direction. Additional stability is obtained by setting the grid points to move along the iso-contour line. To exemplify the application of the method, numerical examples are calculated for the two-dimensional diffusion equation.

On the analysis and application of a spectral collocation scheme for the nonlinear two-dimensional fractional diffusion equation

In this paper, we propose and analyze a novel spectral scheme for the numerical solution of a two-dimensional time-fractional diffusion equation. The proposed scheme approximates the unknown function and its derivatives in space by finite Lucas and Fibonacci polynomial expansion, and time derivative/variable by L1 formula with θ-weighted difference scheme. Error estimates are derived and convergence is proved. Six test problems are then solved to verify the theory numerically. The obtained results confirm the low computational cost and better accuracy than known results in the literature.

Fully Jacobi–Galerkin algorithm for two-dimensional time-dependent PDEs arising in physics

Herein, a pure shifted Jacobi–Galerkin (SJG) method is offered for handling linear two-dimensional space-time diffusion and telegraph equations but by considering their integrated forms. The semi-analytic approximate solution is expanded in spatial and temporal variables as polynomials bases built as a linear combination of the shifted Jacobi Polynomials (JPs). Expanding the exact solution regarding these polynomials satisfies the homogeneous initial and boundary conditions. The proposed Jacobi–Galerkin (JG) procedure yields exponential convergence rates when the solution is smooth enough. A careful investigation of the convergence analysis of the suggested approximate triple series expansion examines the offered numerical algorithm. Two test illustrative examples are given to prove the high accuracy of the proposed scheme.

Matrix splitting preconditioning based on sine transform for solving two-dimensional space-fractional diffusion equations

Finite difference discretization of the two-dimensional space-fractional diffusion equations derives a complicated linear system consisting of identity matrix and four scaled block-Toeplitz with Toeplitz block (BTTB) matrices resulted from the left and right Riemann–Liouville fractional derivatives in different directions. Incorporating with the diffusion coefficients and the symmetric parts of the BTTB matrices, we construct a diagonal and symmetric splitting (DSS) iteration method, which is demonstrated to be convergent conditionally when the considered space-fractional diffusion equations have sufficiently close diffusion coefficients. By further replacing the symmetric Toeplitz matrices involved in the BTTB matrices with τ matrices, an approximated DSS (ADSS) preconditioner based on two-dimensional fast sine transform is designed to accelerate the convergence rates of the Krylov subspace iteration methods. In this way, the total computational complexity of the ADSS-preconditioned GMRES method will be of O(n2logn), where n2 represents the dimension of the corresponding discrete linear system. In addition, theoretical analysis demonstrates that the eigenvalues of the ADSS-preconditioned matrix are weakly clustered around a complex disk centered at 1 with the radius less than 1. Numerical experiments show that the ADSS-preconditioned GMRES method is much more efficient than the other existing methods, and can show h-independent convergence behavior.

On τ-preconditioner for a novel fourth-order difference scheme of two-dimensional Riesz space-fractional diffusion equations

On τ-preconditioner for a novel fourth-order difference scheme of two-dimensional Riesz space-fractional diffusion equations

Numerical Solution of Two-Dimensional Fractional Order Diffusion Equation by using Elzaki Transform with Residual Power Series Method

The Elzaki transform with RPS method is an efficient and reliable approach for solution of non-integer order linear as well as non-linear differential equations. The major endeavourof current work is to find the solution of two dimensional non-integer order diffusion equation by Elzaki transform with RPS method. Firstly, Elzaki is applied on this equation. Secondly, inverse Elzaki is taken on this equation for finding the expressionof solution. Thirdly, assumed approximate solution is substituted on this equation. Then unidentified coefficient functions are obtained by using residual function is equal to zero as well as combining its initial circumstances. At last, coefficient functions are substituted in power series form for finding finite approximate logical solutions. The judgment between exact solution and approximate analytic solutions with different number of terms of this equation are determined and compared for reliability. This method reduces the size of computational works of solution of non-integer order non-linear diffusion equation

Radial Basis Function Based Partition of Unity Method for Two-Dimensional Unsteady Convection Diffusion Equations

Radial Basis Function Based Partition of Unity Method for Two-Dimensional Unsteady Convection Diffusion Equations

A Green’s Function Method based approach for effectiveness factor evaluation in two-dimensional spherical catalyst geometries with non-linear kinetics and partial wetting

The solution of the two-dimensional spherical diffusion equation inside the catalyst particle is proposed using Green’s Function Method (GFM). A Green’s Function was derived using Eigenfunction expansion in the polar direction and considering jump discontinuity in the radial direction. The advantage of using multi-precision arithmetic in the evaluation of Green’s Function has been demonstrated. A pre-computed value of Green’s Function at a selected precision level was used to calculate the concentration profile inside the fully and partially wetted spherical catalyst particle geometry. The non-linear kinetic rate expression-based formulation for concentration profile calculation gives rise to non-linear equations. These non-linear equations were solved in PETSc (Balay et al., 2017) using various combinations of linear pre-conditioners and Krylov methods. A comparison with the Finite Difference Method (FDM) formulation is presented, and the advantage of the proposed GFM-based formulation for convergence and reduction in computation time is demonstrated.

Error estimate of L1-ADI scheme for two-dimensional multi-term time fractional diffusion equation

<abstract><p>A two-dimensional multi-term time fractional diffusion equation $ {D}_{t}^{\mathit{\boldsymbol{\alpha}}} u(x, y, t)- \Delta u(x, y, t) = f(x, y, t) $ is considered in this paper, where $ {D}_{t}^{\mathit{\boldsymbol{\alpha}}} $ is the multi-term time Caputo fractional derivative. To solve the equation numerically, L1 discretisation to each fractional derivative is used on a graded temporal mesh, together with a standard finite difference method for the spatial derivatives on a uniform spatial mesh. We provide a rigorous stability and convergence analysis of a fully discrete L1-ADI scheme for solving the multi-term time fractional diffusion problem. Numerical results show that the error estimate is sharp.</p></abstract>