A Dimension-Splitting Variational Multiscale Element-Free Galerkin Method for Three-Dimensional Singularly Perturbed Convection-Diffusion Problems

By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.

Computational analysis for fractional characterization of coupled convection-diffusion equations arising in MHD fows

The work is devoted to the fractional characterization of time-dependent coupled convection-diffusion systems arising in magnetohydrodynamics(MHD)flows.The time derivative is expressed by means of Caputo’s fractional derivative concept,while the model is solved via the full-spectral method(FSM)and the semi-spectral scheme(SSS).The FSM is based on the operational matrices of derivatives constructed by using higher-order orthogonal polynomials and collocation techniques.The SSS is developed by discretizing the time variable,and the space domain is collocated by using equal points.A detailed comparative analysis is made through graphs for various parameters and tables with existing literature.The contour graphs are made to show the behaviors of the velocity and magnetic fields.The proposed methods are reasonably efficient in examining the behavior of convection-diffusion equations arising in MHD flows,and the concept may be extended for variable order models arising in MHD flows.

THE FINITE DIFFERENCE STREAMLINE DIFFUSION METHODS FOR TIME-DEPENDENT CONVECTION-DIFFUSION EQUATIONS

In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for considered schemes are derived in the norm stronger than L^2-norm.

CHARACTERISTIC GALERKIN METHOD FOR CONVECTION-DIFFUSION EQUATIONS AND IMPLICIT ALGORITHM USING PRECISE INTEGRATION

This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability.

Multidomain pseudospectral methods for nonlinear convection-diffusion equations

Multidomain pseudospectral approximations to nonlinear convection-diffusion equations are considered. The schemes are formulated with the Legendre-Galerkin method but the nonlinear term is collocated at the Legendre/Chebyshev-Gauss-Lobatto points inside each subinterval. Appropriate base functions are introduced so that the matrix of the system is sparse, and the method can be implemented efficiently and in parallel. The stability and the optimal rate of convergence of the methods are proved. Numerical results are given for both the single domain and the multidomain methods to make a comparison.

A NONLINEAR TRANSFORMATION AND A BOUNDARY-INITIAL VALUE PROBLEM FOR A CLASS OF NONLINEAR CONVECTION-DIFFUSION EQUATIONS

With the aid of a nonlinear transformation, a class of nonlinear convection-diffusion PDE in one space dimension is converted into a linear one, the unique solution of a nonlinear boundary-initial value problem for the nonlinear PDE can be exactly expressed by the nonlinear transformation, and several illustrative examples are given.

Analysis of an Il’in Scheme for a System of Singularly Perturbed Convection-Diffusion Equations

In this paper, a numerical solution for a system of singularly perturbed convection-diffusion equations is studied. The system is discretized by the Il’in scheme on a uniform mesh. It is proved that the numerical scheme has first order accuracy, which is uniform with respect to the perturbation parameters. We show that the condition number of the discrete linear system obtained from applying the Il’in scheme for a system of singularly perturbed convection-diffusion equations is O(N) and the relevant coefficient matrix is well conditioned in comparison with the matrices obtained from applying upwind finite difference schemes on this problem. Numerical results confirm the theory of the method.

The Direct Discontinuous Galerkin Methods with Implicit-Explicit Runge-Kutta Time Marching for Linear Convection-Diffusion Problems

In this paper,a fully discrete stability analysis is carried out for the direct discontinuous Galerkin(DDG)methods coupled with Runge-Kutta-type implicit-explicit time marching,for solving one-dimensional linear convection-diffusion problems.In the spatial discretization,both the original DDG methods and the refined DDG methods with interface corrections are considered.In the time discretization,the convection term is treated explicitly and the diffusion term implicitly.By the energy method,we show that the corresponding fully discrete schemes are unconditionally stable,in the sense that the time-step is only required to be upper bounded by a constant which is independent of the mesh size h.Opti-mal error estimate is also obtained by the aid of a special global projection.Numerical experiments are given to verify the stability and accuracy of the proposed schemes.

Exponential B-Spline Solution of Convection-Diffusion Equations

We present an exponential B-spline collocation method for solving convection-diffusion equation with Dirichlet’s type boundary conditions. The method is based on the Crank-Nicolson formulation for time integration and exponential B-spline functions for space integration. Using the Von Neumann method, the proposed method is shown to be unconditionally stable. Numerical experiments have been conducted to demonstrate the accuracy of the current algorithm with relatively minimal computational effort. The results showed that use of the present approach in the simulation is very applicable for the solution of convection-diffusion equation. The current results are also seen to be more accurate than some results given in the literature. The proposed algorithm is seen to be very good alternatives to existing approaches for such physical applications.

Numerical simulation of effect of convection-diffusion on oxygen transport in microcirculation

The entire process of oxygen transport in microcirculation by developing a3 D porous media model is calculated numerically with coupled solid deformation-fluid seepage-convection and diffusion. The principal novelty of the model is that it takes into account volumetric deformation of both capillary and tissues resulting from capillary fluctuation. How solid deformation, fluid seepage, and convection-diffusion combine to affect oxygen transport is examined quantitatively：（1） Solid deformation is more significant in the middle of capillary, where the maximum value of volumetric deformation reaches about 0.5%.（2） Solid deformation has positive influence on the tissue fluid so that it flows more uniformly and causes oxygen to be transported more uniformly, and eventually impacts oxygen concentration by 0.1%–0.5%.（3） Convection-diffusion coupled deformation and seepage has a maximum（16%） and average（3%） increase in oxygen concentration,compared with pure molecular diffusion. Its more significant role is to allow oxygen to be transported more evenly.

NUMERICAL SIMULATION FOR CONVECTION-DIFFUSION PROBLEM WITH PERIODIC MICRO-STRUCTURE

In this article, the convection dominated convection-diffusion problems with the periodic micro-structure are discussed. A two-scale finite element scheme based on the homogenization technique for this kind of problems is provided. The error estimates between the exact solution and the approximation solution, of the homogenized equation or the two-scale finite element scheme are analyzed. It is shown that the scheme provided in this article is convergent for any fixed diffusion coefficient 5, and it may be convergent independent of δ under some conditions. The numerical results demonstrating the theoretical results are presented in this article.

A Comparative Numerical Study of Parabolic Partial Integro-Differential Equation Arising from Convection-Diffusion

This article studies the development of two numerical techniques for solving convection-diffusion type partial integro-differential equation(PIDE)with a weakly singular kernel.Cubic trigonometric B-spline(CTBS)functions are used for interpolation in both methods.The first method is CTBS based collocation method which reduces the PIDE to an algebraic tridiagonal system of linear equations.The other method is CTBS based differential quadrature method which converts the PIDE to a system of ODEs by computing spatial derivatives as weighted sum of function values.An efficient tridiagonal solver is used for the solution of the linear system obtained in the first method as well as for determination of weighting coefficients in the second method.An explicit scheme is employed as time integrator to solve the system of ODEs obtained in the second method.The methods are tested with three nonhomogeneous problems for their validation.Stability,computational efficiency and numerical convergence of the methods are analyzed.Comparison of errors in approximations produced by the present methods versus different values of discretization parameters and convection-diffusion coefficients are made.Convection and diffusion dominant cases are discussed in terms of Peclet number.The results are also compared with cubic B-spline collocation method.

H^1 space-time discontinuous finite element method for convection-diffusion equations

An H1 space-time discontinuous Galerkin （STDG） scheme for convection- diffusion equations in one spatial dimension is constructed and analyzed. This method is formulated by combining the H1 Galerkin method and the space-time discontinuous finite element method that is discontinuous in time and continuous in space. The existence and the uniqueness of the approximate solution are proved. The convergence of the scheme is analyzed by using the techniques in the finite difference and finite element methods. An optimal a-priori error estimate in the L∞ （H1） norm is derived. The numerical exper- iments are presented to verify the theoretical results.

An approach for choosing discretization schemes and grid size based on the convection-diffusion equation

A new approach for selecting proper discretization schemes and grid size is presented. This method is based on the convection-diffusion equation and can provide insight for the Navier-Stokes equation. The approach mainly addresses two aspects, i.e., the practical accuracy of diffusion term discretization and the behavior of high wavenum- ber disturbances. Two criteria are included in this approach. First, numerical diffusion should not affect the theoretical diffusion accuracy near the length scales of interest. This is achieved by requiring numerical diffusion to be smaller than the diffusion discretization error. Second, high wavenumber modes that are.much smaller than the length scales of interest should not be amplified. These two criteria provide a range of suitable scheme combinations for convective flux and diffusive flux and an ideal interval for grid spacing. The effects of time discretization on these criteria are briefly discussed.

A Spectral Method for Convection-Diffusion Equations

In the practical problems such as nuclear waste pollution and seawater intrusion etc., many problems are reduced to solving the convection-diffusion equation, so the research of convection-diffusion equation is of great value. In this work, a spectral method is presented for solving one and two dimensional convection-diffusion equation with source term. The finite difference method is also used to solve the convection diffusion equation. The numerical experiments show that the spectral method is more efficient than other methods for solving the convection-diffusion equation.

Comparison of Fixed Point Methods and Krylov Subspace Methods Solving Convection-Diffusion Equations

The paper first introduces two-dimensional convection-diffusion equation with boundary value condition, later uses the finite difference method to discretize the equation and analyzes positive definite, diagonally dominant and symmetric properties of the discretization matrix. Finally, the paper uses fixed point methods and Krylov subspace methods to solve the linear system and compare the convergence speed of these two methods.

Adaptive Finite Element Method for Steady Convection-Diffusion Equation

This paper examines the numerical solution of the convection-diffusion equation in 2-D. The solution of this equation possesses singularities in the form of boundary or interior layers due to non-smooth boundary conditions. To overcome such singularities arising from these critical regions, the adaptive finite element method is employed. This scheme is based on the streamline diffusion method combined with Neumann-type posteriori estimator. The effectiveness of this approach is illustrated by different examples with several numerical experiments.

ANALYSIS OF THE IMPLICIT-EXPLICIT ULTRA-WEAK DISCONTINUOUS GALERKIN METHOD FOR CONVECTION-DIFFUSION PROBLEMS

In this paper,we first present the optimal error estimates of the semi-discrete ultra-weak discontinuous Galerkin method for solving one-dimensional linear convection-diffusion equations.Then,coupling with a kind of Runge-Kutta type implicit-explicit time discretization which treats the convection term explicitly and the diffusion term implicitly,we analyze the stability and error estimates of the corresponding fully discrete schemes.The fully discrete schemes are proved to be stable if the time-stepτ≤τ0,whereτ0 is a constant independent of the mesh-size h.Furthermore,by the aid of a special projection and a careful estimate for the convection term,the optimal error estimate is also obtained for the third order fully discrete scheme.Numerical experiments are displayed to verify the theoretical results.

Analysis of the local discontinuous Galerkin method with generalized fluxes for one-dimensional nonlinear convection-diffusion systems

In this paper,we present optimal error estimates of the local discontinuous Galerkin method with generalized numerical fluxes for one-dimensional nonlinear convection-diffusion systems.The upwind-biased flux with the adjustable numerical viscosity for the convective term is chosen based on the local characteristic decomposition,which is helpful in resolving discontinuities of degenerate parabolic equations without enforcing any limiting procedure.For the diffusive term,a pair of generalized alternating fluxes is considered.By constructing and analyzing generalized Gauss-Radau projections with respect to different convective or diffusive terms,we derive optimal error estimates for nonlinear convection-diffusion systems with the symmetrizable flux Jacobian and fully nonlinear diffusive problems.Numerical experiments including long time simulations,different boundary conditions and degenerate equations with discontinuous initial data are provided to demonstrate the sharpness of theoretical results.

A MODIFIED WEAK GALERKIN FINITE ELEMENTMETHOD FOR SINGULARLY PERTURBED PARABOLIC CONVECTION-DIFFUSION-REACTION PROBLEMS

In this work,a modified weak Galerkin finite element method is proposed for solving second order linear parabolic singularly perturbed convection-diffusion equations.The key feature of the proposed method is to replace the classical gradient and divergence operators by the modified weak gradient and modified divergence operators,respectively.We apply the backward finite difference method in time and the modified weak Galerkin finite element method in space on uniform mesh.The stability analyses are presented for both semi-discrete and fully-discrete modified weak Galerkin finite element methods.Optimal order of convergences are obtained in suitable norms.We have achieved the same accuracy with the weak Galerkin method while the degrees of freedom are reduced in our method.Various numerical examples are presented to support the theoretical results.It is theoretically and numerically shown that the method is quite stable.