The higher-order numerical scheme of nonlinear advection-diffusion equations is studied in this article, where the space fractional derivatives are evaluated by using weighted and shifted Grünwald difference oper...The higher-order numerical scheme of nonlinear advection-diffusion equations is studied in this article, where the space fractional derivatives are evaluated by using weighted and shifted Grünwald difference operators and combining the compact technique, in the time direction is discretized by the Crank-Nicolson method. Through the energy method, the stability and convergence of the numerical scheme in the sense of L<sub>2</sub>-norm are proved, and the convergence order is . Some examples are given to show that our numerical scheme is effective.展开更多
Atmospheric air pollution turbulent fluxes can be assumed to be proportional to the mean concentration gradient. This assumption, along with the equation of continuity, leads to the advection-diffusion equation. Many ...Atmospheric air pollution turbulent fluxes can be assumed to be proportional to the mean concentration gradient. This assumption, along with the equation of continuity, leads to the advection-diffusion equation. Many models simulating air pollution dispersion are based upon the solution (numerical or analytical) of the advection-diffusion equation as- suming turbulence parameterization for realistic physical scenarios. We present the general steady three-dimensional solution of the advection-diffusion equation considering a vertically inhomogeneous atmospheric boundary layer for arbitrary vertical profiles of wind and eddy-diffusion coefficients. Numerical results and comparison with experimental data are shown.展开更多
A comparison between a non-Gaussian puff model and an advanced time-dependent model to simulate the pollutant dispersion in the Planetary Boundary Layer is presented. The puff model is based on a general technique for...A comparison between a non-Gaussian puff model and an advanced time-dependent model to simulate the pollutant dispersion in the Planetary Boundary Layer is presented. The puff model is based on a general technique for solving the K-equation, using the truncated Gram-Charlier expansion (type A) of the concentration field and finite set equations for the corresponding moments. The other model (named ADMM: Analytical Dispersion Multilayers Model) is an semi- analytical solution to the time-dependent two-dimensional advection-diffusion equation based on a discretization of the PBL in N sub-layers;in each sub-layers the advection-diffusion equation is solved by the Laplace transform technique, considering an average value for eddy diffusivity and the wind speed. A preliminary performance evaluation is shown in the case of continuous emission from an elevated source in a variable boundary layer. Both models were able to correctly reproduce the concentration field measured and so to be used as operative air pollution models.展开更多
In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equi...In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equivalent functional of two-dimensional advection-diffusion problems is formed,the variation method is used to obtain the equation system,and the penalty method is employed to impose the essential boundary conditions.The difference method for twopoint boundary value problems is used to obtain the discrete equations.Then the corresponding formulas of the ICVMM for advection-diffusion problems are presented.Two numerical examples with different node distributions are used to validate and investigate the accuracy and efficiency of the new method in this paper.It is shown that ICVMM is very effective for advection-diffusion problems,and has a good convergent character,accuracy,and computational efficiency.展开更多
Linear fractional map type (LFMT) nonlinear QCA (NLQCA), one of the simplest reversible NLQCA is studied analytically as well as numerically. Linear advection equation or Time Dependent Schrödinger Equation (...Linear fractional map type (LFMT) nonlinear QCA (NLQCA), one of the simplest reversible NLQCA is studied analytically as well as numerically. Linear advection equation or Time Dependent Schrödinger Equation (TDSE) is obtained from the continuum limit of linear QCA. Similarly it is found that some nonlinear advection-diffusion equations including inviscid Burgers equation and porous-medium equation are obtained from LFMT NLQCA.展开更多
The element of pesedospectral-multiwavelet-Galerkin method, and how tocombine it with penalty method for treating boundary conditions are given. Multiwavelet bases don'toverlap on the given scale, and possess the ...The element of pesedospectral-multiwavelet-Galerkin method, and how tocombine it with penalty method for treating boundary conditions are given. Multiwavelet bases don'toverlap on the given scale, and possess the same compact set in a group of several functions, sothey can be directly used to the numerical discretion on the finite interval. Numerical tests showthat general boundary conditions can be enforced with the penalty method, and thatpesedospectral-multiwavelet-Galerkin method can well track the solutions' development. This alsoproves that pesedospectral-multiwavelet-Galerkin method is effective.展开更多
In this paper, a posteriori error estimates for the generalized Schwartz method with Dirichlet boundary conditions on the interfaces for advection-diffusion equation with second order boundary value problems are prove...In this paper, a posteriori error estimates for the generalized Schwartz method with Dirichlet boundary conditions on the interfaces for advection-diffusion equation with second order boundary value problems are proved by using the Euler time scheme combined with Galerkin spatial method. Furthermore, an asymptotic behavior in Sobolev norm is de- duced using Benssoussau-Lions' algorithm. Finally, the results of some numerical experiments are presented to support the theory.展开更多
The complex variable reproducing kernel particle method (CVRKPM) of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape...The complex variable reproducing kernel particle method (CVRKPM) of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape function of a two-dimensional problem is formed with a one-dimensional basis function. The Galerkin weak form is employed to obtain the discretized system equation, and the penalty method is used to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional variable coefficient advection-diffusion problems are obtained. Two numerical examples are given to show that the method in this paper has greater accuracy and computational efficiency than the conventional meshless method such as reproducing the kernel particle method (RKPM) and the element- free Galerkin (EFG) method.展开更多
In this paper, two modified QUICK schemes, namely Q-QUICK and UQ-QUICK, for improving the preci-sion of convective flux approximation are verified in advection-diffusion equation of pollutants on unstruc-tured grids. ...In this paper, two modified QUICK schemes, namely Q-QUICK and UQ-QUICK, for improving the preci-sion of convective flux approximation are verified in advection-diffusion equation of pollutants on unstruc-tured grids. The constructed auxiliary nodes for Q-QUICK/UQ-QUICK are composed of two neighboring nodes plus the next upwind node, the later node is generated from intersection of the line of current neighboring nodes and their corresponding interfaces. 2D unsteady advection-diffusion equation of pollut-ants is conducted for their verifications on unstructured grids. The numerical results show that Q-QUICK and UQ-QUICK have similar computational accuracy to the central difference scheme and similar numerical stability to upwind difference scheme after applying the deferred correction method. In addition, their corre-sponding CPU times are approximately equivalent to those of traditional difference schemes and their abili-ties for adapting high grid deformation are robust.展开更多
A singularly perturbed advection-diffusion two-point Robin boundary value problem whose solution has a single boundary layer is considered. Based on the piecewise linear polynomial approximation, the finite element me...A singularly perturbed advection-diffusion two-point Robin boundary value problem whose solution has a single boundary layer is considered. Based on the piecewise linear polynomial approximation, the finite element method is applied to the problem. Estimation of the error between solution and the finite element approximation is given in energy norm on shishkin-type mesh.展开更多
This paper presents discontinuous Legendre wavelet Galerkin (DLWG) approach for solving one-dimensional advection-diffusion equation (ADE). Variational formulation of this type equation and corresponding numerical flu...This paper presents discontinuous Legendre wavelet Galerkin (DLWG) approach for solving one-dimensional advection-diffusion equation (ADE). Variational formulation of this type equation and corresponding numerical fluxes are devised by utilizing the advantages of both the Legendre wavelet bases and discontinuous Galerkin (DG) method. The distinctive features of the proposed method are its simple applicability for a variety of boundary conditions and able to effectively approximate the solution of PDEs with less storage space and execution. The results of a numerical experiment are provided to verify the efficiency of the designed new technique.展开更多
We construct an approximate Riemann solver for scalar advection-diffusion equations with piecewise polynomial initial data.The objective is to handle advection and diffusion simultaneously to reduce the inherent numer...We construct an approximate Riemann solver for scalar advection-diffusion equations with piecewise polynomial initial data.The objective is to handle advection and diffusion simultaneously to reduce the inherent numerical diffusion produced by the usual advection flux calculations.The approximate solution is based on the weak formulation of the Riemann problem and is solved within a space-time discontinuous Galerkin approach with two subregions.The novel generalized Riemann solver produces piecewise polynomial solutions of the Riemann problem.In conjunction with a recovery polynomial,the Riemann solver is then applied to define the numerical flux within a finite volume method.Numerical results for a piecewise linear and a piecewise parabolic approximation are shown.These results indicate a reduction in numerical dissipation compared with the conventional separated flux calculation of advection and diffusion.Also,it is shown that using the proposed solver only in the vicinity of discontinuities gives way to an accurate and efficient finite volume scheme.展开更多
By transforming the Caputo tempered fractional advection-diffusion equation into the Riemann–Liouville tempered fractional advection-diffusion equation,and then using the fractional-compact Grünwald–Letnikov te...By transforming the Caputo tempered fractional advection-diffusion equation into the Riemann–Liouville tempered fractional advection-diffusion equation,and then using the fractional-compact Grünwald–Letnikov tempered difference operator to approximate the Riemann–Liouville tempered fractional partial derivative,the fractional central difference operator to discritize the space Riesz fractional partial derivative,and the classical central difference formula to discretize the advection term,a numerical algorithm is constructed for solving the Caputo tempered fractional advection-diffusion equation.The stability and the convergence analysis of the numerical method are given.Numerical experiments show that the numerical method is effective.展开更多
The finite difference method such as alternating group iterative methods is useful in numerical method for evolutionary equations and this is the standard approach taken in this paper. Alternating group explicit (AGE)...The finite difference method such as alternating group iterative methods is useful in numerical method for evolutionary equations and this is the standard approach taken in this paper. Alternating group explicit (AGE) iterative methods for one-dimensional convection diffusion equations problems are given. The stability and convergence are analyzed by the linear method. Numerical results of the model problem are taken. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show that the behavior of the method with emphasis on treatment of boundary conditions is valuable.展开更多
According to Godunov theorem for numerical calculations of advection equations, there exist no high-er-order schemes with constant positive difference coefficients in a family of polynomial schemes with an accuracy ex...According to Godunov theorem for numerical calculations of advection equations, there exist no high-er-order schemes with constant positive difference coefficients in a family of polynomial schemes with an accuracy exceeding the first-order. In case of advection-diffusion equations, so far there have been not found stable schemes with positive difference coefficients in a family of numerical schemes exceeding the second-order accuracy. We propose a third-order computational scheme for numerical fluxes to guarantee the non-negative difference coefficients of resulting finite difference equations for advection-diffusion equations. The present scheme is optimized so as to minimize truncation errors for the numerical fluxes while fulfilling the positivity condition of the difference coefficients which are variable depending on the local Courant number and diffusion number. The feature of the present optimized scheme consists in keeping the third-order accuracy anywhere without any numerical flux limiter by using the same stencil number as convemtional third-order shemes such as KAWAMURA and UTOPIA schemes. We extend the present method into multi-dimensional equations. Numerical experiments for linear and nonlinear advection-diffusion equations were performed and the present scheme’s applicability to nonlinear Burger’s equation was confirmed.展开更多
Flood wave propagation modeling is of critical importance to advancing water resources management and protecting human life and property. In this study, we investigated how the advection-diffusion routing model perfor...Flood wave propagation modeling is of critical importance to advancing water resources management and protecting human life and property. In this study, we investigated how the advection-diffusion routing model performed in flood wave propagation on a 16 km long downstream section of the Big Piney River, MO. Model performance was based on gaging station data at the upstream and downstream cross sections. We demonstrated with advection-diffusion theory that for small differences in watershed drainage area between the two river cross sections, inflow along the reach mainly contributes to the downstream hydrograph's rising limb and not to the falling limb. The downstream hydrograph's falling limb is primarily determined by the propagated flood wave originating at the upstream cross section. This research suggests the parameter for the advectiondiffusion routing model can be calibrated by fitting the hydrograph falling limb. Application of the advection diffusion model to the flood wave of January 29, 2013 supports our theoretical finding that the propagated flood wave determines the downstream cross section falling limb, and the model has good performance in our test examples.展开更多
We propose and analyze a posteriori energy-norm error estimates for weighted interior penalty discontinuous Galerkin approximations of advection-diffusion-reaction equations with heterogeneous and anisotropic diffusio...We propose and analyze a posteriori energy-norm error estimates for weighted interior penalty discontinuous Galerkin approximations of advection-diffusion-reaction equations with heterogeneous and anisotropic diffusion. The weights, which play a key role in the analysis, depend on the diffusion tensor and are used to formulate the consistency terms in the discontinuous Galerkin method. The error upper bounds, in which all the constants are specified, consist of three terms: a residual estimator which depends only on the elementwise fluctuation of the discrete solution residual, a diffusive flux estimator where the weights used in the method enter explicitly, and a non-conforming estimator which is nonzero because of the use of discontinuous finite element spaces. The three estimators can be bounded locally by the approximation error. A particular attention is given to the dependency on problem parameters of the constants in the local lower error bounds. For moderate advection, it is shown that full robustness with respect to diffusion heterogeneities is achieved owing to the specific design of the weights in the discontinuous Galerkin method, while diffusion anisotropies remain purely local and impact the constants through the square root of the condition number of the diffusion tensor. For dominant advection, it is shown, in the spirit of previous work by Verfiirth on continuous finite elements, that the local lower error bounds can be written with constants involving a cut-off for the ratio of local mesh size to the reciprocal of the square root of the lowest local eignevalue of the diffusion tensor.展开更多
We develop a family of characteristic discontinuous Galerkin methods for transient advection-diffusion equations,including the characteristic NIPG,OBB,IIPG,and SIPG schemes.The derived schemes possess combined advanta...We develop a family of characteristic discontinuous Galerkin methods for transient advection-diffusion equations,including the characteristic NIPG,OBB,IIPG,and SIPG schemes.The derived schemes possess combined advantages of EulerianLagrangian methods and discontinuous Galerkin methods.An optimal-order error estimate in the L2 norm and a superconvergence estimate in a weighted energy norm are proved for the characteristic NIPG,IIPG,and SIPG scheme.Numerical experiments are presented to confirm the optimal-order spatial and temporal convergence rates of these schemes as proved in the theorems and to show that these schemes compare favorably to the standard NIPG,OBB,IIPG,and SIPG schemes in the context of advection-diffusion equations.展开更多
Solving the advection-diffusion equation in irregular geometries is of great importance for realistic simulations. To this end, we adopt multi-patch parameterizations to describe irregular geometries. Different from t...Solving the advection-diffusion equation in irregular geometries is of great importance for realistic simulations. To this end, we adopt multi-patch parameterizations to describe irregular geometries. Different from the classical multi-patch parameterization method, C 1- continuity is introduced in order to avoid designing interface conditions between adjacent patches. However, singularities of parameterizations can’t always be avoided. Thus, in this paper, a finite volume method is proposed based on smooth multi-patch singular parameterizations. It is called a singular parameterized finite volume method. Firstly, we present a numerical scheme for pure advection equation and pure diffusion equation respectively. Secondly, numerical stability results in L2 norm show that the numerical method is not suffered from the singularities. Thirdly, the numerical method has second order accurate in L2 norm. Finally, three numerical tests in different irregular geometries are presented to show efficiency of this numerical method.展开更多
We present a higher-order in-space characteristic method for the solution of the tran- sient advection diffusion equations in two space dimensions. This method uses biquadratic trial and test functions within the fram...We present a higher-order in-space characteristic method for the solution of the tran- sient advection diffusion equations in two space dimensions. This method uses biquadratic trial and test functions within the framework of the Eulerian-Lagrangian localized Adjoint Methods (ELLAM). It therefore maintains the advantages of previous ELLAM schemes. Namely, it treats general boundary conditions naturally in a systematic manner, conserves mass, and symmetrizes the governing transport equations. Moreover, it generates accurate numerical solutions even if large time steps are used in the simulation. Numerical exper- iments are presented to illustrate the performance of this method and establish its order of convergence numerically.展开更多
文摘The higher-order numerical scheme of nonlinear advection-diffusion equations is studied in this article, where the space fractional derivatives are evaluated by using weighted and shifted Grünwald difference operators and combining the compact technique, in the time direction is discretized by the Crank-Nicolson method. Through the energy method, the stability and convergence of the numerical scheme in the sense of L<sub>2</sub>-norm are proved, and the convergence order is . Some examples are given to show that our numerical scheme is effective.
文摘Atmospheric air pollution turbulent fluxes can be assumed to be proportional to the mean concentration gradient. This assumption, along with the equation of continuity, leads to the advection-diffusion equation. Many models simulating air pollution dispersion are based upon the solution (numerical or analytical) of the advection-diffusion equation as- suming turbulence parameterization for realistic physical scenarios. We present the general steady three-dimensional solution of the advection-diffusion equation considering a vertically inhomogeneous atmospheric boundary layer for arbitrary vertical profiles of wind and eddy-diffusion coefficients. Numerical results and comparison with experimental data are shown.
文摘A comparison between a non-Gaussian puff model and an advanced time-dependent model to simulate the pollutant dispersion in the Planetary Boundary Layer is presented. The puff model is based on a general technique for solving the K-equation, using the truncated Gram-Charlier expansion (type A) of the concentration field and finite set equations for the corresponding moments. The other model (named ADMM: Analytical Dispersion Multilayers Model) is an semi- analytical solution to the time-dependent two-dimensional advection-diffusion equation based on a discretization of the PBL in N sub-layers;in each sub-layers the advection-diffusion equation is solved by the Laplace transform technique, considering an average value for eddy diffusivity and the wind speed. A preliminary performance evaluation is shown in the case of continuous emission from an elevated source in a variable boundary layer. Both models were able to correctly reproduce the concentration field measured and so to be used as operative air pollution models.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project,China(Grant No. S30106)the Innovation Fund for Graduate Student of Shanghai University,China (Grant No. SHUCX120125)
文摘In this paper,an improved complex variable meshless method(ICVMM) for two-dimensional advection-diffusion problems is developed based on improved complex variable moving least-square(ICVMLS) approximation.The equivalent functional of two-dimensional advection-diffusion problems is formed,the variation method is used to obtain the equation system,and the penalty method is employed to impose the essential boundary conditions.The difference method for twopoint boundary value problems is used to obtain the discrete equations.Then the corresponding formulas of the ICVMM for advection-diffusion problems are presented.Two numerical examples with different node distributions are used to validate and investigate the accuracy and efficiency of the new method in this paper.It is shown that ICVMM is very effective for advection-diffusion problems,and has a good convergent character,accuracy,and computational efficiency.
文摘Linear fractional map type (LFMT) nonlinear QCA (NLQCA), one of the simplest reversible NLQCA is studied analytically as well as numerically. Linear advection equation or Time Dependent Schrödinger Equation (TDSE) is obtained from the continuum limit of linear QCA. Similarly it is found that some nonlinear advection-diffusion equations including inviscid Burgers equation and porous-medium equation are obtained from LFMT NLQCA.
基金This project is supported by National Natural Science Foundation of China(No. 19971020) Multidiseipline Scientific Research Foundation of Harbin Institute of Technology, China(No.HIT.MD2001.26).
文摘The element of pesedospectral-multiwavelet-Galerkin method, and how tocombine it with penalty method for treating boundary conditions are given. Multiwavelet bases don'toverlap on the given scale, and possess the same compact set in a group of several functions, sothey can be directly used to the numerical discretion on the finite interval. Numerical tests showthat general boundary conditions can be enforced with the penalty method, and thatpesedospectral-multiwavelet-Galerkin method can well track the solutions' development. This alsoproves that pesedospectral-multiwavelet-Galerkin method is effective.
文摘In this paper, a posteriori error estimates for the generalized Schwartz method with Dirichlet boundary conditions on the interfaces for advection-diffusion equation with second order boundary value problems are proved by using the Euler time scheme combined with Galerkin spatial method. Furthermore, an asymptotic behavior in Sobolev norm is de- duced using Benssoussau-Lions' algorithm. Finally, the results of some numerical experiments are presented to support the theory.
基金supported by the National Natural Science Foundation of China (Grant No. 11171208)the Leading Academic Discipline Project of Shanghai City,China (Grant No. S30106)
文摘The complex variable reproducing kernel particle method (CVRKPM) of solving two-dimensional variable coefficient advection-diffusion problems is presented in this paper. The advantage of the CVRKPM is that the shape function of a two-dimensional problem is formed with a one-dimensional basis function. The Galerkin weak form is employed to obtain the discretized system equation, and the penalty method is used to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional variable coefficient advection-diffusion problems are obtained. Two numerical examples are given to show that the method in this paper has greater accuracy and computational efficiency than the conventional meshless method such as reproducing the kernel particle method (RKPM) and the element- free Galerkin (EFG) method.
文摘In this paper, two modified QUICK schemes, namely Q-QUICK and UQ-QUICK, for improving the preci-sion of convective flux approximation are verified in advection-diffusion equation of pollutants on unstruc-tured grids. The constructed auxiliary nodes for Q-QUICK/UQ-QUICK are composed of two neighboring nodes plus the next upwind node, the later node is generated from intersection of the line of current neighboring nodes and their corresponding interfaces. 2D unsteady advection-diffusion equation of pollut-ants is conducted for their verifications on unstructured grids. The numerical results show that Q-QUICK and UQ-QUICK have similar computational accuracy to the central difference scheme and similar numerical stability to upwind difference scheme after applying the deferred correction method. In addition, their corre-sponding CPU times are approximately equivalent to those of traditional difference schemes and their abili-ties for adapting high grid deformation are robust.
文摘A singularly perturbed advection-diffusion two-point Robin boundary value problem whose solution has a single boundary layer is considered. Based on the piecewise linear polynomial approximation, the finite element method is applied to the problem. Estimation of the error between solution and the finite element approximation is given in energy norm on shishkin-type mesh.
文摘This paper presents discontinuous Legendre wavelet Galerkin (DLWG) approach for solving one-dimensional advection-diffusion equation (ADE). Variational formulation of this type equation and corresponding numerical fluxes are devised by utilizing the advantages of both the Legendre wavelet bases and discontinuous Galerkin (DG) method. The distinctive features of the proposed method are its simple applicability for a variety of boundary conditions and able to effectively approximate the solution of PDEs with less storage space and execution. The results of a numerical experiment are provided to verify the efficiency of the designed new technique.
基金This work was supported by the German Research Foundation(DFG)through the Collaborative Research Center SFB TRR 75 Droplet Dynamics Under Extreme Ambient Conditions
文摘We construct an approximate Riemann solver for scalar advection-diffusion equations with piecewise polynomial initial data.The objective is to handle advection and diffusion simultaneously to reduce the inherent numerical diffusion produced by the usual advection flux calculations.The approximate solution is based on the weak formulation of the Riemann problem and is solved within a space-time discontinuous Galerkin approach with two subregions.The novel generalized Riemann solver produces piecewise polynomial solutions of the Riemann problem.In conjunction with a recovery polynomial,the Riemann solver is then applied to define the numerical flux within a finite volume method.Numerical results for a piecewise linear and a piecewise parabolic approximation are shown.These results indicate a reduction in numerical dissipation compared with the conventional separated flux calculation of advection and diffusion.Also,it is shown that using the proposed solver only in the vicinity of discontinuities gives way to an accurate and efficient finite volume scheme.
文摘By transforming the Caputo tempered fractional advection-diffusion equation into the Riemann–Liouville tempered fractional advection-diffusion equation,and then using the fractional-compact Grünwald–Letnikov tempered difference operator to approximate the Riemann–Liouville tempered fractional partial derivative,the fractional central difference operator to discritize the space Riesz fractional partial derivative,and the classical central difference formula to discretize the advection term,a numerical algorithm is constructed for solving the Caputo tempered fractional advection-diffusion equation.The stability and the convergence analysis of the numerical method are given.Numerical experiments show that the numerical method is effective.
文摘The finite difference method such as alternating group iterative methods is useful in numerical method for evolutionary equations and this is the standard approach taken in this paper. Alternating group explicit (AGE) iterative methods for one-dimensional convection diffusion equations problems are given. The stability and convergence are analyzed by the linear method. Numerical results of the model problem are taken. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show that the behavior of the method with emphasis on treatment of boundary conditions is valuable.
文摘According to Godunov theorem for numerical calculations of advection equations, there exist no high-er-order schemes with constant positive difference coefficients in a family of polynomial schemes with an accuracy exceeding the first-order. In case of advection-diffusion equations, so far there have been not found stable schemes with positive difference coefficients in a family of numerical schemes exceeding the second-order accuracy. We propose a third-order computational scheme for numerical fluxes to guarantee the non-negative difference coefficients of resulting finite difference equations for advection-diffusion equations. The present scheme is optimized so as to minimize truncation errors for the numerical fluxes while fulfilling the positivity condition of the difference coefficients which are variable depending on the local Courant number and diffusion number. The feature of the present optimized scheme consists in keeping the third-order accuracy anywhere without any numerical flux limiter by using the same stencil number as convemtional third-order shemes such as KAWAMURA and UTOPIA schemes. We extend the present method into multi-dimensional equations. Numerical experiments for linear and nonlinear advection-diffusion equations were performed and the present scheme’s applicability to nonlinear Burger’s equation was confirmed.
基金supported by funding from the USDA Forest Service Northern Research Station iTree Spatial Simulation (No. PL-5937)the National Urban and Community Forest Advisory Council iT ree Tool (No. 11-DG-11132544340)The SUNY ESF Department of Environmental Resources Engineering provided computing facilities and logistical support
文摘Flood wave propagation modeling is of critical importance to advancing water resources management and protecting human life and property. In this study, we investigated how the advection-diffusion routing model performed in flood wave propagation on a 16 km long downstream section of the Big Piney River, MO. Model performance was based on gaging station data at the upstream and downstream cross sections. We demonstrated with advection-diffusion theory that for small differences in watershed drainage area between the two river cross sections, inflow along the reach mainly contributes to the downstream hydrograph's rising limb and not to the falling limb. The downstream hydrograph's falling limb is primarily determined by the propagated flood wave originating at the upstream cross section. This research suggests the parameter for the advectiondiffusion routing model can be calibrated by fitting the hydrograph falling limb. Application of the advection diffusion model to the flood wave of January 29, 2013 supports our theoretical finding that the propagated flood wave determines the downstream cross section falling limb, and the model has good performance in our test examples.
文摘We propose and analyze a posteriori energy-norm error estimates for weighted interior penalty discontinuous Galerkin approximations of advection-diffusion-reaction equations with heterogeneous and anisotropic diffusion. The weights, which play a key role in the analysis, depend on the diffusion tensor and are used to formulate the consistency terms in the discontinuous Galerkin method. The error upper bounds, in which all the constants are specified, consist of three terms: a residual estimator which depends only on the elementwise fluctuation of the discrete solution residual, a diffusive flux estimator where the weights used in the method enter explicitly, and a non-conforming estimator which is nonzero because of the use of discontinuous finite element spaces. The three estimators can be bounded locally by the approximation error. A particular attention is given to the dependency on problem parameters of the constants in the local lower error bounds. For moderate advection, it is shown that full robustness with respect to diffusion heterogeneities is achieved owing to the specific design of the weights in the discontinuous Galerkin method, while diffusion anisotropies remain purely local and impact the constants through the square root of the condition number of the diffusion tensor. For dominant advection, it is shown, in the spirit of previous work by Verfiirth on continuous finite elements, that the local lower error bounds can be written with constants involving a cut-off for the ratio of local mesh size to the reciprocal of the square root of the lowest local eignevalue of the diffusion tensor.
文摘We develop a family of characteristic discontinuous Galerkin methods for transient advection-diffusion equations,including the characteristic NIPG,OBB,IIPG,and SIPG schemes.The derived schemes possess combined advantages of EulerianLagrangian methods and discontinuous Galerkin methods.An optimal-order error estimate in the L2 norm and a superconvergence estimate in a weighted energy norm are proved for the characteristic NIPG,IIPG,and SIPG scheme.Numerical experiments are presented to confirm the optimal-order spatial and temporal convergence rates of these schemes as proved in the theorems and to show that these schemes compare favorably to the standard NIPG,OBB,IIPG,and SIPG schemes in the context of advection-diffusion equations.
基金The authors are grateful to the anonymous reviewers for their useful comments and suggestionsThe authors are supported by the NSF of China (No. 11601114, No.11401138)the Anhui Provincial Natural Science Foundation (No. 1608085QA14).
文摘Solving the advection-diffusion equation in irregular geometries is of great importance for realistic simulations. To this end, we adopt multi-patch parameterizations to describe irregular geometries. Different from the classical multi-patch parameterization method, C 1- continuity is introduced in order to avoid designing interface conditions between adjacent patches. However, singularities of parameterizations can’t always be avoided. Thus, in this paper, a finite volume method is proposed based on smooth multi-patch singular parameterizations. It is called a singular parameterized finite volume method. Firstly, we present a numerical scheme for pure advection equation and pure diffusion equation respectively. Secondly, numerical stability results in L2 norm show that the numerical method is not suffered from the singularities. Thirdly, the numerical method has second order accurate in L2 norm. Finally, three numerical tests in different irregular geometries are presented to show efficiency of this numerical method.
文摘We present a higher-order in-space characteristic method for the solution of the tran- sient advection diffusion equations in two space dimensions. This method uses biquadratic trial and test functions within the framework of the Eulerian-Lagrangian localized Adjoint Methods (ELLAM). It therefore maintains the advantages of previous ELLAM schemes. Namely, it treats general boundary conditions naturally in a systematic manner, conserves mass, and symmetrizes the governing transport equations. Moreover, it generates accurate numerical solutions even if large time steps are used in the simulation. Numerical exper- iments are presented to illustrate the performance of this method and establish its order of convergence numerically.