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...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.展开更多
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 nonlin...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.展开更多
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 prec...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.展开更多
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 c...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 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/C...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.展开更多
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 an...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.展开更多
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 th...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.展开更多
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 s...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.展开更多
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 sch...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.展开更多
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 gre...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 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 dom...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.展开更多
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 approx...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]展开更多
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...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.展开更多
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 algorit...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.展开更多
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 differe...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.展开更多
For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of ...For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.展开更多
The compact second-order upwind finite difference schemes free of ceil Reynolds number limitation are developed in this paper for the one-to three-dimensional steady convection- diffusion equations,using a perturbatio...The compact second-order upwind finite difference schemes free of ceil Reynolds number limitation are developed in this paper for the one-to three-dimensional steady convection- diffusion equations,using a perturbational technique applied to the classical first-order upwind schemes.The present second-order schemes take essentially the same form as those of the first- order schemes,but involve a simple modification to the diffusive coefficients.Numerical exam- ples including one-to three-dimensional model equations of fluid flow and a problem of natural convection with boundary-layer effect are given to illustrate the excellent behavior of the present schemes.展开更多
A discontinuous finite element method for convection-diffusion equations is proposed and analyzed. This scheme is designed to produce an approximate solution which is completely discontinuous. Optimal order of converg...A discontinuous finite element method for convection-diffusion equations is proposed and analyzed. This scheme is designed to produce an approximate solution which is completely discontinuous. Optimal order of convergence is obtained for model problem. This is the same convergence rate known for the classical methods. [ABSTRACT FROM AUTHOR]展开更多
Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based...Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based on these bilinear equations, some exact solutions to the three modified equations are derived.展开更多
In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled b...In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.12250410244,11872151)the Jiangsu Province Education Development Special Project-2022 for Double First-ClassSchool Talent Start-up Fund of China(No.2022r109)the Longshan Scholar Program of Jiangsu Province of China。
文摘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.
文摘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.
文摘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.
基金Project supported by National Natural Science Foundation of China and China State Key project for Basic Researchcs.
文摘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.
基金supported by the National Natural Science Foundation of China(No.60874039)the Leading Academic Discipline Project of Shanghai Municipal Education Commission(No.J50101)
文摘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.
文摘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.
基金Natural Science Foundation of Gansu Province of China
文摘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.
基金Project supported by the National Natural Science Foundation of China (No. 11061021)the Inner Mongolia College Research Project (No. NJ10006)the Natural Science Foundation of Inner Mongolia of China (No. 2012MS0106)
文摘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.
文摘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.
文摘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 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.
文摘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]
基金This work was supported by the Foundation of Gansu Natural Science (3ZS041-A25-011).
文摘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.
基金partially supported by the National Natural Science Foundation of China(Nos.11871009,12271055,12171048)the foundation of CAEP(CX20210044)the Foundation of LCP.
文摘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.
基金NSFC Basic Science Center Program for”Multiscale Problems in Nonlinear Mechanics”(Grant No.11988102)National Natural Science Foundation of China(Grant No.12202454).
文摘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.
基金supported by National Natural Science Foundation of China(11771257)the Shandong Provincial Natural Science Foundation of China(ZR2023YQ002,ZR2023MA007,ZR2021MA004)。
文摘For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.
文摘The compact second-order upwind finite difference schemes free of ceil Reynolds number limitation are developed in this paper for the one-to three-dimensional steady convection- diffusion equations,using a perturbational technique applied to the classical first-order upwind schemes.The present second-order schemes take essentially the same form as those of the first- order schemes,but involve a simple modification to the diffusive coefficients.Numerical exam- ples including one-to three-dimensional model equations of fluid flow and a problem of natural convection with boundary-layer effect are given to illustrate the excellent behavior of the present schemes.
文摘A discontinuous finite element method for convection-diffusion equations is proposed and analyzed. This scheme is designed to produce an approximate solution which is completely discontinuous. Optimal order of convergence is obtained for model problem. This is the same convergence rate known for the classical methods. [ABSTRACT FROM AUTHOR]
基金supported by the National Natural Science Foundation of China (Grant Nos. 11931017 and 12071447)。
文摘Three modified sine-Hilbert(sH)-type equations, i.e., the modified sH equation, the modified damped sH equation, and the modified nonlinear dissipative system, are proposed, and their bilinear forms are provided.Based on these bilinear equations, some exact solutions to the three modified equations are derived.
文摘In this study, the Bernstein collocation method has been expanded to Stancu collocation method for numerical solution of the charged particle motion for certain configurations of oscillating magnetic fields modelled by a class of linear integro-differential equations. As the method has been improved, the Stancu polynomials that are generalization of the Bernstein polynomials have been used. The method has been tested on a physical problem how the method can be applied. Moreover, numerical results of the method have been compared with the numerical results of the other methods to indicate the efficiency of the method.