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 propose...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.展开更多
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.展开更多
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.展开更多
This paper presents a finite element procedure for solving tran-sient, multidimensional convection-diffusion equations. The procedure is based onthe characteristic Galerkin method with an implicit algorithm using prec...This paper presents a finite element procedure for solving tran-sient, multidimensional convection-diffusion equations. The procedure is based onthe characteristic Galerkin method with an implicit algorithm using precise integra-tion method. With the operator splitting procedure, the precise integration methodis introduced to determine the material derivative in the convection-diffusion equa-tion, consequently, the physical quantities of material points. An implicit algorithmwith a combination of both the precise and the traditional numerical integration pro-cedures in time domain in the Lagrange coordinates for the characteristic Galerkinmethod is formulated. The stability analysis of the algorithm shows that the uncondi-tional stability of present implicit algorithm is enhanced as compared with that of thetraditional implicit numerical integration procedure. The numerical results validatethe presented method in solving convection-diffusion equations. As compared withSUPG method and explicit characteristic Galerkin method, the present method givesthe results with higher accuracy and better stability.展开更多
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 nnique solution of a nonlinear boundary-initial value problem for t...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 nnique 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.展开更多
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.展开更多
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.展开更多
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 conve...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.展开更多
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)functio...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.展开更多
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.展开更多
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 condi...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.展开更多
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 R...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.展开更多
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 ...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.展开更多
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 replac...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.展开更多
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 compactsecond-order upwind scheme proposed by the a...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 compactsecond-order upwind scheme proposed by the authors.The basic fourth-order scheme,which likethe 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 examplesincluding one-to three-dimensional model equations,with available analytical solutions,of fluidflow and a problem,with benchmark solutions,of natural convective heat transfer are given toillustrate 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,thefourth-order accuracy is specially verified using double precision arithmetic.展开更多
The Perturbational Finite Difference (PFD) method is a kind of high-order-accurate compact difference method, But its idea is different from the normal compact method and the multi-nodes method. This method can get a ...The Perturbational Finite Difference (PFD) method is a kind of high-order-accurate compact difference method, But its idea is different from the normal compact method and the multi-nodes method. This method can get a Perturbational Exact Numerical Solution (PENS) scheme for locally linearlized Convection-Diffusion (CD) equation. The PENS scheme is similar to the Finite Analytical (FA) scheme and Exact Difference Solution (EDS) scheme, which are all exponential schemes, but PENS scheme is simpler and uses only 3, 5 and 7 nodes for 1-, 2- and 3-dimensional problems, respectively. The various approximate schemes of PENS scheme are also called Perturbational-High-order-accurate Difference (PHD) scheme. The PHD schemes can be got by expanding the exponential terms in the PENS scheme into power series of grid Renold number, and they are all upwind schemes and remain the concise structure form of first-order upwind scheme. For 1-dimensional (1-D) CD equation and 2-D incompressible Navier-Stokes equation, their PENS and PHD schemes were constituted in this paper, they all gave highly accurate results for the numerical examples of three 1-D CD equations and an incompressible 2-D flow in a square cavity.展开更多
We consider an optimal control problem with an 1D singularly perturbed differential state equation.For solving such problems one uses the enhanced system of the state equation and its adjoint form.Thus,we obtain a sys...We consider an optimal control problem with an 1D singularly perturbed differential state equation.For solving such problems one uses the enhanced system of the state equation and its adjoint form.Thus,we obtain a system of two convectiondiffusion equations.Using linear finite elements on adapted grids we treat the effects of two layers arising at different boundaries of the domain.We proof uniform error estimates for this method on meshes of Shishkin type.We present numerical results supporting our analysis.展开更多
An operator-splitting algorithm for the three-dimensional convection-diffusion equa-tion is presented.The flow region is discretized into tetrahedronal elements which are fixed in time.The transport equation is split ...An operator-splitting algorithm for the three-dimensional convection-diffusion equa-tion is presented.The flow region is discretized into tetrahedronal elements which are fixed in time.The transport equation is split into two successive initial value problems:a pure convection problemand a pure diffusion problem.For the pure convection problem,solutions are found by the method ofcharacteristiCS.The solution algorithm involves tracing the characteristic lines backwards in time froma vertex of an element to an interior point.A cubic polynomial is used to interpolate the concentrationand its derivatives within each element.For the diffusion problem,an explicit finite element algorithmis employed.Numerical examples are given which agree well with the analytical solutions.展开更多
This paper presents an exponential compact higher order scheme for Convection-Diffusion Equations(CDE)with variable and nonlinear convection coeffi-cients.The scheme is O(h4)for one-dimensional problems and produces a...This paper presents an exponential compact higher order scheme for Convection-Diffusion Equations(CDE)with variable and nonlinear convection coeffi-cients.The scheme is O(h4)for one-dimensional problems and produces a tri-diagonal system of equations which can be solved efficiently using Thomas algorithm.For twodimensional problems,the scheme produces an O(h4+k4)accuracy over a compact nine point stencil which can be solved using any line iterative approach with alternate direction implicit procedure.The convergence of the iterative procedure is guaranteed as the coefficient matrix of the developed scheme satisfies the conditions required to be positive.Wave number analysis has been carried out to establish that the scheme is comparable in accuracy with spectral methods.The higher order accuracy and better rate of convergence of the developed scheme have been demonstrated by solving numerous model problems for one and two-dimensional CDE,where the solutions have the sharp gradient at the solution boundary.展开更多
基金supported by the Natural Science Foundation of Zhejiang Province,China(Grant Nos.LY20A010021,LY19A010002,LY20G030025)the Natural Science Founda-tion of Ningbo City,China(Grant Nos.2021J147,2021J235).
文摘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.
基金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.
基金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.
文摘This paper presents a finite element procedure for solving tran-sient, multidimensional convection-diffusion equations. The procedure is based onthe characteristic Galerkin method with an implicit algorithm using precise integra-tion method. With the operator splitting procedure, the precise integration methodis introduced to determine the material derivative in the convection-diffusion equa-tion, consequently, the physical quantities of material points. An implicit algorithmwith a combination of both the precise and the traditional numerical integration pro-cedures in time domain in the Lagrange coordinates for the characteristic Galerkinmethod is formulated. The stability analysis of the algorithm shows that the uncondi-tional stability of present implicit algorithm is enhanced as compared with that of thetraditional implicit numerical integration procedure. The numerical results validatethe presented method in solving convection-diffusion equations. As compared withSUPG method and explicit characteristic Galerkin method, the present method givesthe results with higher accuracy and better stability.
基金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 nnique 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.
文摘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.
文摘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.
基金the NSFC grant 11871428the Nature Science Research Program for Colleges and Universities of Jiangsu Province grant 20KJB110011Qiang Zhang:Research supported by the NSFC grant 11671199。
文摘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.
文摘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.
文摘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.
文摘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.
基金Research sponsored by NSFC grants 11871428 and 12071214Nature Science Research Program for Colleges and Universities of Jiangsu Province grant 20KJB110011+1 种基金Research is supported in part by NSFC grants U1930402the fellowship of China Postdoctoral Science Foundation(No.2020TQ0030).
文摘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.
基金supported by National Natural Science Foundation of China(Grant Nos.11971132 and 11971131)Natural Science Foundation of Heilongjiang Province(Grant No.YQ2021A002)Guangdong Basic and Applied Basic Research Foundation(Grant No.2020B1515310006)。
文摘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.
基金supported in part by National Natural Science Foundation of China (No.11871038).
文摘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.
文摘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 compactsecond-order upwind scheme proposed by the authors.The basic fourth-order scheme,which likethe 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 examplesincluding one-to three-dimensional model equations,with available analytical solutions,of fluidflow and a problem,with benchmark solutions,of natural convective heat transfer are given toillustrate 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,thefourth-order accuracy is specially verified using double precision arithmetic.
文摘The Perturbational Finite Difference (PFD) method is a kind of high-order-accurate compact difference method, But its idea is different from the normal compact method and the multi-nodes method. This method can get a Perturbational Exact Numerical Solution (PENS) scheme for locally linearlized Convection-Diffusion (CD) equation. The PENS scheme is similar to the Finite Analytical (FA) scheme and Exact Difference Solution (EDS) scheme, which are all exponential schemes, but PENS scheme is simpler and uses only 3, 5 and 7 nodes for 1-, 2- and 3-dimensional problems, respectively. The various approximate schemes of PENS scheme are also called Perturbational-High-order-accurate Difference (PHD) scheme. The PHD schemes can be got by expanding the exponential terms in the PENS scheme into power series of grid Renold number, and they are all upwind schemes and remain the concise structure form of first-order upwind scheme. For 1-dimensional (1-D) CD equation and 2-D incompressible Navier-Stokes equation, their PENS and PHD schemes were constituted in this paper, they all gave highly accurate results for the numerical examples of three 1-D CD equations and an incompressible 2-D flow in a square cavity.
文摘We consider an optimal control problem with an 1D singularly perturbed differential state equation.For solving such problems one uses the enhanced system of the state equation and its adjoint form.Thus,we obtain a system of two convectiondiffusion equations.Using linear finite elements on adapted grids we treat the effects of two layers arising at different boundaries of the domain.We proof uniform error estimates for this method on meshes of Shishkin type.We present numerical results supporting our analysis.
文摘An operator-splitting algorithm for the three-dimensional convection-diffusion equa-tion is presented.The flow region is discretized into tetrahedronal elements which are fixed in time.The transport equation is split into two successive initial value problems:a pure convection problemand a pure diffusion problem.For the pure convection problem,solutions are found by the method ofcharacteristiCS.The solution algorithm involves tracing the characteristic lines backwards in time froma vertex of an element to an interior point.A cubic polynomial is used to interpolate the concentrationand its derivatives within each element.For the diffusion problem,an explicit finite element algorithmis employed.Numerical examples are given which agree well with the analytical solutions.
基金Nachiketa Mishra is greatly indebted to the Council of Scientific and Industrial Research for the financial support 09/084(0389)/2006-EMR-I.
文摘This paper presents an exponential compact higher order scheme for Convection-Diffusion Equations(CDE)with variable and nonlinear convection coeffi-cients.The scheme is O(h4)for one-dimensional problems and produces a tri-diagonal system of equations which can be solved efficiently using Thomas algorithm.For twodimensional problems,the scheme produces an O(h4+k4)accuracy over a compact nine point stencil which can be solved using any line iterative approach with alternate direction implicit procedure.The convergence of the iterative procedure is guaranteed as the coefficient matrix of the developed scheme satisfies the conditions required to be positive.Wave number analysis has been carried out to establish that the scheme is comparable in accuracy with spectral methods.The higher order accuracy and better rate of convergence of the developed scheme have been demonstrated by solving numerous model problems for one and two-dimensional CDE,where the solutions have the sharp gradient at the solution boundary.
