In this paper we make a close study of the finite analytic method by means of the maximum principles in differential equations and give the proof of the stability and convergence of the finite analytic method.
In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existenc...In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.展开更多
The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain ...The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain decomposition method for two dimensional time dependent convection diffusion equations with variable coefficients. By combining predictor-corrector technique,modified upwind differences with explicitimplicit coupling,the method under consideration provides intrinsic parallelism while maintaining good stability and accuracy.Moreover,for multi-dimensional problems, the method can be readily implemented on a multi-processor system and does not have the limitation on the choice of subdomains required by some other similar predictor-corrector or stabilized schemes.These properties of the method are demonstrated in this work through both rigorous mathematical analysis and numerical experiments.展开更多
Aim To construct a third order upwind scheme for convection equation. Methods Upwind Lagrange interpolation was used. Results and Conclusion The schemes L p stability for p∈ is proved. Numerical exam...Aim To construct a third order upwind scheme for convection equation. Methods Upwind Lagrange interpolation was used. Results and Conclusion The schemes L p stability for p∈ is proved. Numerical examples show that performance of the third order upwind scheme is better than that of most second order schemes.展开更多
A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order...A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.展开更多
This paper studies the time asymptotic behavior of solutions for a nonlinear convection diffusion reaction equation in one dimension.First,the pointwise estimates of solutions are obtained,furthermore,we obtain the op...This paper studies the time asymptotic behavior of solutions for a nonlinear convection diffusion reaction equation in one dimension.First,the pointwise estimates of solutions are obtained,furthermore,we obtain the optimal Lp,1≤ p ≤ +∞,convergence rate of solutions for small initial data.Then we establish the local existence of solutions,the blow up criterion and the sufficient condition to ensure the nonnegativity of solutions for large initial data.Our approach is based on the detailed analysis of the Green function of the linearized equation and some energy estimates.展开更多
We investigate the linearization of systems of n-component nonlinear diffusion equations; such systems have physical applications in soil science, mathematical biology and invariant curve flows. Equivalence transforma...We investigate the linearization of systems of n-component nonlinear diffusion equations; such systems have physical applications in soil science, mathematical biology and invariant curve flows. Equivalence transformations of their auxiliary systems are used to identify the systems that can be linearized. We also provide several examples of systems with two-component equations, and show how to linearize them by nonlocal mappings.展开更多
In this paper,we investigate a stochastic meshfree finite volume element method for an optimal control problem governed by the convection diffusion equations with random coefficients.There are two contributions of thi...In this paper,we investigate a stochastic meshfree finite volume element method for an optimal control problem governed by the convection diffusion equations with random coefficients.There are two contributions of this paper.Firstly,we establish a scheme to approximate the optimality system by using the finite volume element method in the physical space and the meshfree method in the probability space,which is competitive for high-dimensional random inputs.Secondly,the a priori error estimates are derived for the state,the co-state and the control variables.Some numerical tests are carried out to confirm the theoretical results and demonstrate the efficiency of the proposed method.展开更多
In this article, the author studies the Cauchy problem of the damped wave equation with a nonlinear convection term in multi-dimensions. The author shows that a classical solution to the Cauchy problem exists globally...In this article, the author studies the Cauchy problem of the damped wave equation with a nonlinear convection term in multi-dimensions. The author shows that a classical solution to the Cauchy problem exists globally in time under smallness condition on the initial perturbation. Furthermore, the author obtains the L^p (2 ≤ p ≤ ∞) decay estimates of the solution.展开更多
In this paper ,in the space that possesses restoring nucleus, we obtain analyticsolutions in the series form for the steady-state convection diffusion equation The solutions have the following characteristics: (1) the...In this paper ,in the space that possesses restoring nucleus, we obtain analyticsolutions in the series form for the steady-state convection diffusion equation The solutions have the following characteristics: (1) they ave given in the accurate form:(2)they can be calculated in the explicit way, without solving the eguations;(3) the error of the approximate solution will be monotonically decreased under the meaning of the norm of the spaces when a cardinal term is added in the procedure of numerical solution .Finally, we calculated the example in [2] the result shows that our solution is more accurate than that in [2].展开更多
To reduce computational costs, an improved form of the frequency domain boundary element method(BEM) is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitr...To reduce computational costs, an improved form of the frequency domain boundary element method(BEM) is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation(BIE) representation solves the two-dimensional convected Helmholtz equation(CHE) and its fundamental solution, which must satisfy a new Sommerfeld radiation condition(SRC) in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green's kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.展开更多
The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted...The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted with finite element method and time with finite difference scheme was derived, where the numerical solution of velocity, pressure, and temperature can be found together, and a numerical example to simulate the close square cavity is given, which is of practical importance.展开更多
Transport of nonreactive solutes in soils is principally controlled by soil properties, such as particle-size distribution and pore geometry. Surface tension of soil water yields capillary forces that bind the water i...Transport of nonreactive solutes in soils is principally controlled by soil properties, such as particle-size distribution and pore geometry. Surface tension of soil water yields capillary forces that bind the water in the soil pores. Changes in soil water surface tension by contaminants may affect flow of soil water due to decreased capillary forces, caused by lowered soil water surface tension. This study aimed at assessing solute transport in sand columns as affected by effluent surface tension. Miscible displacement (MD) tests were conducted on sand columns repacked with sands sieved from 2.0, 1.0, 0.5 and 0.25 mm screens. The MD tests were conducted with 0.05 M bromide solutions prepared using water with surface tension adjusted to 72.8, 64, 53.5 and 42 dyne/cm2. Obtained breakthrough curves were modeled with the convection-dispersion equation (CDE) model. Coefficient of hydrodynamic dispersion and pore-water velocity responded inconsistently across decreased particle-sizes and water surface tensions and this was attributed to non-uniform effect of lowered effluent surface tension on solute transport in different pore-size distribution.展开更多
In this paper,ETD3-Padéand ETD4-PadéGalerkin finite element methods are proposed and analyzed for nonlinear delayed convection-diffusion-reaction equations with Dirichlet boundary conditions.An ETD-based RK ...In this paper,ETD3-Padéand ETD4-PadéGalerkin finite element methods are proposed and analyzed for nonlinear delayed convection-diffusion-reaction equations with Dirichlet boundary conditions.An ETD-based RK is used for time integration of the corresponding equation.To overcome a well-known difficulty of numerical instability associated with the computation of the exponential operator,the Padéapproach is used for such an exponential operator approximation,which in turn leads to the corresponding ETD-Padéschemes.An unconditional L^(2) numerical stability is proved for the proposed numerical schemes,under a global Lipshitz continuity assumption.In addition,optimal rate error estimates are provided,which gives the convergence order of O(k^(3)+h^(r))(ETD3-Padé)or O(k^(4)+h^(r))(ETD4-Padé)in the L^(2)norm,respectively.Numerical experiments are presented to demonstrate the robustness of the proposed numerical schemes.展开更多
A linear convection equation with discontinuous coefficients arises in wave propagation through interfaces. An interface condition is needed at the interface to select a unique solution. An upwind scheme that builds t...A linear convection equation with discontinuous coefficients arises in wave propagation through interfaces. An interface condition is needed at the interface to select a unique solution. An upwind scheme that builds this interface condition into its numerical flux is called the immersed interface upwind scheme. An ι1-error estimate of such a scheme was first established by Wen et al. (2008). In this paper, we provide a simple analysis on the ι1-error estimate. The main idea is to formulate the solution to the underline initial-value problem into the sum of solutions to two convection equations with constant coefficients, which can then be estimated using classical methods for the initial or boundary value problems.展开更多
In this paper, we study variational discretization for the constrained optimal control problem governed by convection dominated diffusion equations, where the state equation is approximated by the edge stabilization G...In this paper, we study variational discretization for the constrained optimal control problem governed by convection dominated diffusion equations, where the state equation is approximated by the edge stabilization Galerkin method. A priori error estimates are derived for the state, the adjoint state and the control. Moreover, residual type a posteriori error estimates in the L^2-norm are obtained. Finally, two numerical experiments are presented to illustrate the theoretical results.展开更多
In this paper we study the source-type solution for the heat equation with convection: ut = △u + b·▽un for (x,t) ∈ ST→ RN × (0,T] and u(x,0) = δ(x) for x ∈ RN, where δ(x) denotes Dirac meas...In this paper we study the source-type solution for the heat equation with convection: ut = △u + b·▽un for (x,t) ∈ ST→ RN × (0,T] and u(x,0) = δ(x) for x ∈ RN, where δ(x) denotes Dirac measure in = RN,N 2,n 0 and b = (b1,...,bN) ∈ RN is a vector. It is shown that there exists a critical number pc = N+2 such that the source-type solution to the above problem exists and is unique if 0 N n 〈 pc and there exists a unique similarity source-type solution in the case n = N+1 , while such a solution does not exist N if n 〉 pc. Moreover, the asymptotic behavior of the solution near the origin is studied. It is shown that when 0 〈 n 〈 N+1 the convection is too weak and the short time behavior of the source-type solution near the origin N is the same as that for the heat equation without convection.展开更多
Examines the convection diffusion problems using domain decomposition method. Presentation of continuous and discrete convection diffusion equations; Kinds of multiplicative Schwarz algorithms; Optimal order error est...Examines the convection diffusion problems using domain decomposition method. Presentation of continuous and discrete convection diffusion equations; Kinds of multiplicative Schwarz algorithms; Optimal order error estimate results.展开更多
文摘In this paper we make a close study of the finite analytic method by means of the maximum principles in differential equations and give the proof of the stability and convergence of the finite analytic method.
基金supported by the National Basic Research Program under the Grant 2005CB321701the National Natural Science Foundation of China under the Grants 60474027 and 10771211.
文摘In this paper,we investigate a streamline diffusion finite element approxi- mation scheme for the constrained optimal control problem governed by linear con- vection dominated diffusion equations.We prove the existence and uniqueness of the discretized scheme.Then a priori and a posteriori error estimates are derived for the state,the co-state and the control.Three numerical examples are presented to illustrate our theoretical results.
基金the National Natural Science Foundation of China(No.10571017)supported in part by the National Natural Science Foundation of China(No.60533020)supported in part by NSF DMS 0712744
文摘The numerical solution of large scale multi-dimensional convection diffusion equations often requires efficient parallel algorithms.In this work,we consider the extension of a recently proposed non-overlapping domain decomposition method for two dimensional time dependent convection diffusion equations with variable coefficients. By combining predictor-corrector technique,modified upwind differences with explicitimplicit coupling,the method under consideration provides intrinsic parallelism while maintaining good stability and accuracy.Moreover,for multi-dimensional problems, the method can be readily implemented on a multi-processor system and does not have the limitation on the choice of subdomains required by some other similar predictor-corrector or stabilized schemes.These properties of the method are demonstrated in this work through both rigorous mathematical analysis and numerical experiments.
文摘Aim To construct a third order upwind scheme for convection equation. Methods Upwind Lagrange interpolation was used. Results and Conclusion The schemes L p stability for p∈ is proved. Numerical examples show that performance of the third order upwind scheme is better than that of most second order schemes.
基金supported by the National Natural Science Foundation of China (No. 10601022)NSF ofInner Mongolia Autonomous Region of China (No. 200607010106)513 and Science Fund of InnerMongolia University for Distinguished Young Scholars (No. ND0702)
文摘A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.
文摘This paper studies the time asymptotic behavior of solutions for a nonlinear convection diffusion reaction equation in one dimension.First,the pointwise estimates of solutions are obtained,furthermore,we obtain the optimal Lp,1≤ p ≤ +∞,convergence rate of solutions for small initial data.Then we establish the local existence of solutions,the blow up criterion and the sufficient condition to ensure the nonnegativity of solutions for large initial data.Our approach is based on the detailed analysis of the Green function of the linearized equation and some energy estimates.
基金Supported by the National Natural Science Foundation of China under Grant No 10671156, and the Programme for New Century Excellent Talents in University (NCET-04-0968).
文摘We investigate the linearization of systems of n-component nonlinear diffusion equations; such systems have physical applications in soil science, mathematical biology and invariant curve flows. Equivalence transformations of their auxiliary systems are used to identify the systems that can be linearized. We also provide several examples of systems with two-component equations, and show how to linearize them by nonlocal mappings.
基金supported by the National Natural Science Foundation of China(Nos.11701253,11971259,11801216)Natural Science Foundation of Shandong Province(No.ZR2017BA010)。
文摘In this paper,we investigate a stochastic meshfree finite volume element method for an optimal control problem governed by the convection diffusion equations with random coefficients.There are two contributions of this paper.Firstly,we establish a scheme to approximate the optimality system by using the finite volume element method in the physical space and the meshfree method in the probability space,which is competitive for high-dimensional random inputs.Secondly,the a priori error estimates are derived for the state,the co-state and the control variables.Some numerical tests are carried out to confirm the theoretical results and demonstrate the efficiency of the proposed method.
基金supported by Shanghai Municipal Natural Science Foundation 09ZR1413500National Natural Science Foundation of China 11071162
文摘In this article, the author studies the Cauchy problem of the damped wave equation with a nonlinear convection term in multi-dimensions. The author shows that a classical solution to the Cauchy problem exists globally in time under smallness condition on the initial perturbation. Furthermore, the author obtains the L^p (2 ≤ p ≤ ∞) decay estimates of the solution.
文摘In this paper ,in the space that possesses restoring nucleus, we obtain analyticsolutions in the series form for the steady-state convection diffusion equation The solutions have the following characteristics: (1) they ave given in the accurate form:(2)they can be calculated in the explicit way, without solving the eguations;(3) the error of the approximate solution will be monotonically decreased under the meaning of the norm of the spaces when a cardinal term is added in the procedure of numerical solution .Finally, we calculated the example in [2] the result shows that our solution is more accurate than that in [2].
基金supported by National Engineering School of Tunis (No.13039.1)
文摘To reduce computational costs, an improved form of the frequency domain boundary element method(BEM) is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation(BIE) representation solves the two-dimensional convected Helmholtz equation(CHE) and its fundamental solution, which must satisfy a new Sommerfeld radiation condition(SRC) in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green's kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.
文摘The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted with finite element method and time with finite difference scheme was derived, where the numerical solution of velocity, pressure, and temperature can be found together, and a numerical example to simulate the close square cavity is given, which is of practical importance.
文摘Transport of nonreactive solutes in soils is principally controlled by soil properties, such as particle-size distribution and pore geometry. Surface tension of soil water yields capillary forces that bind the water in the soil pores. Changes in soil water surface tension by contaminants may affect flow of soil water due to decreased capillary forces, caused by lowered soil water surface tension. This study aimed at assessing solute transport in sand columns as affected by effluent surface tension. Miscible displacement (MD) tests were conducted on sand columns repacked with sands sieved from 2.0, 1.0, 0.5 and 0.25 mm screens. The MD tests were conducted with 0.05 M bromide solutions prepared using water with surface tension adjusted to 72.8, 64, 53.5 and 42 dyne/cm2. Obtained breakthrough curves were modeled with the convection-dispersion equation (CDE) model. Coefficient of hydrodynamic dispersion and pore-water velocity responded inconsistently across decreased particle-sizes and water surface tensions and this was attributed to non-uniform effect of lowered effluent surface tension on solute transport in different pore-size distribution.
文摘In this paper,ETD3-Padéand ETD4-PadéGalerkin finite element methods are proposed and analyzed for nonlinear delayed convection-diffusion-reaction equations with Dirichlet boundary conditions.An ETD-based RK is used for time integration of the corresponding equation.To overcome a well-known difficulty of numerical instability associated with the computation of the exponential operator,the Padéapproach is used for such an exponential operator approximation,which in turn leads to the corresponding ETD-Padéschemes.An unconditional L^(2) numerical stability is proved for the proposed numerical schemes,under a global Lipshitz continuity assumption.In addition,optimal rate error estimates are provided,which gives the convergence order of O(k^(3)+h^(r))(ETD3-Padé)or O(k^(4)+h^(r))(ETD4-Padé)in the L^(2)norm,respectively.Numerical experiments are presented to demonstrate the robustness of the proposed numerical schemes.
基金supported by National Science Foundation of USA(Grant No.DMS1114546)
文摘A linear convection equation with discontinuous coefficients arises in wave propagation through interfaces. An interface condition is needed at the interface to select a unique solution. An upwind scheme that builds this interface condition into its numerical flux is called the immersed interface upwind scheme. An ι1-error estimate of such a scheme was first established by Wen et al. (2008). In this paper, we provide a simple analysis on the ι1-error estimate. The main idea is to formulate the solution to the underline initial-value problem into the sum of solutions to two convection equations with constant coefficients, which can then be estimated using classical methods for the initial or boundary value problems.
基金support of the Chinese and German Research Foundations through the Sino-German Workshop on Applied Mathematics held in Hangzhou in October 2007support of the German Research Foundation through the grants DFG06-381 and DFG06-382+1 种基金support of the National Basic Research Program under the Grant 2005CB321701the National Natural Science Foundation of China under the Grant 60474027 and 10771211
文摘In this paper, we study variational discretization for the constrained optimal control problem governed by convection dominated diffusion equations, where the state equation is approximated by the edge stabilization Galerkin method. A priori error estimates are derived for the state, the adjoint state and the control. Moreover, residual type a posteriori error estimates in the L^2-norm are obtained. Finally, two numerical experiments are presented to illustrate the theoretical results.
基金supported by National Natural Science Foundation of China (Grant Nos.10671103, 11001142)the Natural Science Foundation of Fujian Province (Grant No. S0650027)
文摘In this paper we study the source-type solution for the heat equation with convection: ut = △u + b·▽un for (x,t) ∈ ST→ RN × (0,T] and u(x,0) = δ(x) for x ∈ RN, where δ(x) denotes Dirac measure in = RN,N 2,n 0 and b = (b1,...,bN) ∈ RN is a vector. It is shown that there exists a critical number pc = N+2 such that the source-type solution to the above problem exists and is unique if 0 N n 〈 pc and there exists a unique similarity source-type solution in the case n = N+1 , while such a solution does not exist N if n 〉 pc. Moreover, the asymptotic behavior of the solution near the origin is studied. It is shown that when 0 〈 n 〈 N+1 the convection is too weak and the short time behavior of the source-type solution near the origin N is the same as that for the heat equation without convection.
基金National Natural Science Foundation of China (No. 10071044) and Major Basic Researches Program of China.
文摘Examines the convection diffusion problems using domain decomposition method. Presentation of continuous and discrete convection diffusion equations; Kinds of multiplicative Schwarz algorithms; Optimal order error estimate results.