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.

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.

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.

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.

Two-grid method for characteristic mixed finite-element solutions of nonlinear convection-diffusion equations

A two-grid method for solving nonlinear convection-dominated diffusion equations is presented. The method use discretizations based on a characteristic mixed finite-element method and give the linearization for nonlinear systems by two steps. The error analysis shows that the two-grid scheme combined with the characteristic mixed finite-element method can decrease numerical oscillation caused by dominated convections and solve nonlinear advection-dominated diffusion problems efficiently.

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.

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.

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.

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.

A Finite Difference Scheme on a Priori Adapted Meshes for a Singularly Perturbed Parabolic Convection-Diffusion Equation

A boundary value problem is considered for a singularly perturbed parabolic convection-diffusion equation; we construct a finite difference scheme on α priori （sequentially） adapted meshes and study its convergence. The scheme on α priori adapted meshes is constructed using a majorant function for the singular component of the discrete solution, which allows us to find α priori a subdomain where the computed solution requires a further improvement. This subdomain is defined by the perturbation parameter ε, the step-size of a uniform mesh in χ, and also by the required accuracy of the discrete solution and the prescribed number of refinement iterations K for improving the solution. To solve the discrete problems aimed at the improvement of the solution, we use uniform meshes on the subdomains. The error of the numerical solution depends weakly on the parameter ε. The scheme converges almost ε-uniformly, precisely, under the condition N^-1 = o （ε^v）, where N denotes the number of nodes in the spatial mesh, and the value v = v（K） can be chosen arbitrarily small for suitable K.

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.

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.

UPWIND SPLITTING SCHEME FOR CONVECTION-DIFFUSION EQUATIONS

WT5,5”BX] A new class of numerical schemes is proposed to solve convection diffusion equations by combining the upwind technique and the method of operator splitting. For every time step, the multi dimensional approximation is performed in several independent directions alternatively, while the upwind technique is applied to treat the convection term in every individual direction. This scheme possesses maximum principle. Stability and convergence are analysed by energy method.[WT5,5”HZ]

Incremental Unknowns Method for Solving Three-Dimensional Convection-Diffusion Equations

We use the incremental unknowns method in conjunction with the iterative methods to approximate the solution of the nonsymmetric and positive-definite linear systems generated from a multilevel discretization of three-dimensional convection-diffusion equations. The condition numbers of incremental unknowns matrices associated with the convection-diffusion equations and the number of iterations needed to attain an acceptable accuracy are estimated. Numerical results are presented with two-level approximations, which demonstrate that the incremental unknowns method when combined with some iter- ative methods is very effcient.

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.

An Interpolation-Free Cell-Centered Finite Volume Scheme for 3D Anisotropic Convection-Diffusion Equations on Arbitrary Polyhedral Meshes

Most existing cell-centered finite volume schemes need to introduce auxiliary unknowns in order to maintain the second-order accuracy when the mesh is distorted or the problem is discontinuous,so interpolation algorithms of auxiliary unknowns are required.Interpolation algorithms are not only difficult to construct,but also bring extra computation.In this paper,an interpolation-free cell-centered finite volume scheme is proposed for the heterogeneous and anisotropic convectiondiffusion problems on arbitrary polyhedral meshes.We propose a new interpolationfree discretization method for diffusion term,and two new second-order upwind algorithms for convection term.Most interestingly,the scheme can be adapted to any mesh topology and can handle any discontinuity strictly.Numerical experiments show that this new scheme is robust,possesses a small stencil,and has approximately secondorder accuracy for both diffusion-dominated and convection-dominated problems.

Quantum Implementation of Numerical Methods for Convection-Diffusion Equations:Toward Computational Fluid Dynamics

We present quantum numerical methods for the typical initial boundary value problems(IBVPs)of convection-diffusion equations in fluid dynamics.The IBVP is discretized into a series of linear systems via finite difference methods and explicit time marching schemes.To solve these discrete systems in quantum computers,we design a series of quantum circuits,including four stages of encoding,amplification,adding source terms,and incorporating boundary conditions.In the encoding stage,the initial condition is encoded in the amplitudes of quantum registers as a state vector to take advantage of quantum algorithms in space complexity.In the following three stages,the discrete differential operators in classical computing are converted into unitary evolutions to satisfy the postulate in quantum systems.The related arithmetic calculations in quantum amplitudes are also realized to sum up the increments from these stages.The proposed quantum algorithm is implemented within the open-source quantum computing framework Qiskit[2].By simulating one-dimensional transient problems,including the Helmholtz equation,the Burgers’equation,and Navier-Stokes equations,we demonstrate the capability of quantum computers in fluid dynamics.

FINITE PROXIMATE METHOD FOR CONVECTION-DIFFUSION EQUATION

A finite proximate method was presented to solve the convection-diffusion equation in curvilinear grids. The method has characteristics of automatic upwind effect and the good stability. It was verified through exact solution and other calculation results of two-dimensional dam-break flow in a frictionless, horizontal channel. The calculation results are in good agreement with the exact solution and other calculation results, which show that the finite proximate method can be applied to solve the convection-diffusion equation directly not only in the rectangular grids, but also in the curvilinear grids.

A PERTURBATIONAL FOURTH-ORDER UPWIND FINITE DIFFERENCE SCHEME FOR THE CONVECTION-DIFFUSION EQUATION

In this study,a compact fourth-order upwind finite difference scheme for the con- vection-diffusion equation is developed,by the scheme perturbation technique and the compact second-order upwind scheme proposed by the authors.The basic fourth-order scheme,which like the classical upwind scheme is free of cell Reynolds-number limitation in terms of spurious oscil- lation and involves only immediate neighbouring nodal points,is presented for the one-dimen- sional equation,and subsequently generalized to multi-dimensional cases.Numerical examples including one-to three-dimensional model equations,with available analytical solutions,of fluid flow and a problem,with benchmark solutions,of natural convective heat transfer are given to illustrate the excellent behavior in such aspects as accuracy,resolution to‘shock wave’-and ‘boundary layer’-effects in convection dominant cases,of the present scheme.Besides,the fourth-order accuracy is specially verified using double precision arithmetic.