Highly accurate spatial discretization is essentially required to perform numerical climate and weather prediction. The difference between the differential and the finite-difference operator is however a primitive err...Highly accurate spatial discretization is essentially required to perform numerical climate and weather prediction. The difference between the differential and the finite-difference operator is however a primitive error source in the numerics. This paper presents an optimization of centered finite-difference operator based on the principle of constrained cost function, which can reduce the truncation error to minimum. In the optimization point of view, such optimal operator is in fact an attempt to minimize spatial truncation er-rors in atmospheric modeling, in a simple way and indeed a quite innovative way to implement Variational Continuous Assimilation (VCA) technique. Furthermore, the optimizing difference operator is consciously designed to be meshing-independent, so that it can be used for most Arakawa-mesh configurations, such as un-staggered (Arakawa-A) or com-monly staggered (Arakawa-B, Arakawa-C, Arakawa-D) mesh. But for the calibration purpose, the pro-posed operator is implemented on an un-staggered mesh in which the truncation oscillation is mostly ex-cited, and it thus makes a severe and indeed a benchmark test for the proposed optimal scheme. Both theo-retical investigation and practical modeling indicate that the aforementioned numerical noise can be significantly eliminated.展开更多
The truncation error associated with a given sampling representation is defined as the difference between the signal and on approximating sumutilizing a finite number of terms. In this paper we give uniform bound for ...The truncation error associated with a given sampling representation is defined as the difference between the signal and on approximating sumutilizing a finite number of terms. In this paper we give uniform bound for truncation error of bandlimited functions in the n dimensional Lebesgue space Lp(Rn) associated with multidimensional Shannon sampling representation.展开更多
A new so called truncation error reduction method (TERM) is developed in this work. This is an iterative process which uses a coarse grid (2 h ) to estimate the truncation error and then reduces the error on the or...A new so called truncation error reduction method (TERM) is developed in this work. This is an iterative process which uses a coarse grid (2 h ) to estimate the truncation error and then reduces the error on the original grid ( h ). The purpose is to use coarse grids to get more accurate results and to develop a new method which could do coarse grid direct numerical simulation (DNS) for more accurate and acceptable DNS solutions.展开更多
This paper investigates double sampling series derivatives for bivariate functions defined on R2 that are in the Bernstein space. For this sampling series, we estimate some of the pointwise and uniform bounds when the...This paper investigates double sampling series derivatives for bivariate functions defined on R2 that are in the Bernstein space. For this sampling series, we estimate some of the pointwise and uniform bounds when the function satisfies some decay conditions. The truncated series of this formula allow us to approximate any order of partial derivatives for function from Bernstein space using only a finite number of samples from the function itself. This sampling formula will be useful in the approximation theory and its applications, especially after having the truncation error well-established. Examples with tables and figures are given at the end of the paper to illustrate the advantages of this formula.展开更多
The damage identification is made by the numerical simulation analysis of a five-storey-and-two-span RC frame structure, using improved and unimproved direct analytical method respectively; and the fundamental equatio...The damage identification is made by the numerical simulation analysis of a five-storey-and-two-span RC frame structure, using improved and unimproved direct analytical method respectively; and the fundamental equations were solved by the minimal least square method （viz. general inverse method）. It demonstrates that the feasibility and the accuracy of the present approach were impoved significantly, compared with the result of unimproved damage identification.展开更多
Let B^pΩ, 1 ≤ p 〈 ∞, be the space of all bounded functions from Lp（R） which can be extended to entire functions of exponential type Ω. The uniform error bounds for truncated Whittaker-Kotelnikov-Shannon series ...Let B^pΩ, 1 ≤ p 〈 ∞, be the space of all bounded functions from Lp（R） which can be extended to entire functions of exponential type Ω. The uniform error bounds for truncated Whittaker-Kotelnikov-Shannon series based on local sampling are derived for functions f ∈ B^pΩ without decay assumption at infinity. Then the optimal bounds of the aliasing error and truncation error of Whittaker-Kotelnikov-Shannon expansion for non-bandlimited functions from Sobolev classes L/（Wp（R）） are determined up to a logarithmic factor.展开更多
Let B^pΩ,1≤Р≤∞,be the set of all bounded functions in L^p（R）which can be extended to entire unctions of exponential typeΩ. The unitbrm bounds for truncation error of Shannon sampling expansion fromlocal averag...Let B^pΩ,1≤Р≤∞,be the set of all bounded functions in L^p（R）which can be extended to entire unctions of exponential typeΩ. The unitbrm bounds for truncation error of Shannon sampling expansion fromlocal averages are obtained for functions f∈BpΩwith the decay condition f（t）≤A/t^δ,t≠0,where A and δare positive constants. Furthermore we also establish similar results for non-bandlimit functions in Besov classes with the same decay condition as above.展开更多
This paper establishes relations between the stability and the high-order truncated corrections for modeling of the mass conservation equation with the tworelaxation-times(TRT)collision operator.First we propose a sim...This paper establishes relations between the stability and the high-order truncated corrections for modeling of the mass conservation equation with the tworelaxation-times(TRT)collision operator.First we propose a simple method to derive the truncation errors from the exact,central-difference type,recurrence equations of the TRT scheme.They also supply its equivalent three-time-level discretization form.Two different relationships of the two relaxation rates nullify the third(advection)and fourth(pure diffusion)truncation errors,for any linear equilibrium and any velocity set.However,the two relaxation times alone cannot remove the leading-order advection-diffusion error,because of the intrinsic fourth-order numerical diffusion.The truncation analysis is carefully verified for the evolution of concentration waves with the anisotropic diffusion tensors.The anisotropic equilibrium functions are presented in a simple but general form,suitable for the minimal velocity sets and the d2Q9,d3Q13,d3Q15 and d3Q19 velocity sets.All anisotropic schemes are complemented by their exact necessary von Neumann stability conditions and equivalent finite-difference stencils.The sufficient stability conditions are proposed for the most stable(OTRT)family,which enables modeling at any Peclet numbers with the same velocity amplitude.The heuristic stability analysis of the fourth-order truncated corrections extends the optimal stability to larger relationships of the two relaxation rates,in agreementwith the exact(one-dimensional)and numerical(multi-dimensional)stability analysis.A special attention is put on the choice of the equilibrium weights.By combining accuracy and stability predictions,several strategies for selecting the relaxation and free-tunable equilibrium parameters are suggested and applied to the evolution of the Gaussian hill.展开更多
A family of symmetric (hybrid) two step sixth P-stable methods for the accurate numerical integration of second order periodic initial value problems have been considered in this paper. These methods, which require on...A family of symmetric (hybrid) two step sixth P-stable methods for the accurate numerical integration of second order periodic initial value problems have been considered in this paper. These methods, which require only three (new) function evaluation per iteration and per step integration. These methods have minimal local truncation error (LTE) and smaller phase-lag of sixth order than some sixth orders P-stable methods in [1-3,10-11]. The theoretical and numerical results show that these methods in this paper are more accurate and efficient than some methods proposed in [1-3,10].展开更多
This paper presents truncation errors among Corrector Formula for left Rectangular rule and Corrector Formula for middle Rectangular rule respectively. It also displays an analysis on convergence order of compound cor...This paper presents truncation errors among Corrector Formula for left Rectangular rule and Corrector Formula for middle Rectangular rule respectively. It also displays an analysis on convergence order of compound corrector formulas for rectangular rule. Examples of numerical calculation have validated theoretical analysis.展开更多
The truncation error of improved Cotes formula is presented in this paper. It also displays an analysis on convergence order of improved Cotes formula. Examples of numerical calculation is given in the end.
In this paper,we establish a new multivariate Hermite sampling series involving samples from the function itself and its mixed and non-mixed partial derivatives of arbitrary order.This multivariate form of Hermite sam...In this paper,we establish a new multivariate Hermite sampling series involving samples from the function itself and its mixed and non-mixed partial derivatives of arbitrary order.This multivariate form of Hermite sampling will be valid for some classes of multivariate entire functions,satisfying certain growth conditions.We will show that many known results included in Commun Korean Math Soc,2002,17:731-740,Turk J Math,2017,41:387-403 and Filomat,2020,34:3339-3347 are special cases of our results.Moreover,we estimate the truncation error of this sampling based on localized sampling without decay assumption.Illustrative examples are also presented.展开更多
The solution of time-dependent hyperbolic conservation laws on cut cell meshes causes the small cell problem:standard schemes are not stable on the arbitrarily small cut cells if an explicit time stepping scheme is us...The solution of time-dependent hyperbolic conservation laws on cut cell meshes causes the small cell problem:standard schemes are not stable on the arbitrarily small cut cells if an explicit time stepping scheme is used and the time step size is chosen based on the size of the background cells.In May and Berger(J Sci Comput 71:919–943,2017),the mixed explicit-implicit approach in general and MUSCL-Trap(monotonic upwind scheme for conservation laws and trapezoidal scheme)in particular have been introduced to solve this problem by using implicit time stepping on the cut cells.Theoretical and numerical results have indicated that this might lead to a loss in accuracy when switching between the explicit and implicit time stepping.In this contribution,we examine this in more detail and will prove in one dimension that the specific combination MUSCL-Trap of an explicit second-order and an implicit second-order scheme results in a fully second-order mixed scheme.As this result is unlikely to hold in two dimensions,we also introduce two new versions of mixed explicit-implicit schemes based on exchanging the explicit scheme.We present numerical tests in two dimensions where we compare the new versions with the original MUSCL-Trap scheme.展开更多
The classic state methods for trajectory estimation in boost phase with multi-range-rate system include method of point-by-point manner and that of spline-model-based manner. Both are deficient in terms of model-appro...The classic state methods for trajectory estimation in boost phase with multi-range-rate system include method of point-by-point manner and that of spline-model-based manner. Both are deficient in terms of model-approximation accuracy and systematic error determination thus resulting in the estimation errors well beyond the requirements, especially, concerning the maneuvering trajectory. This article proposes a new high-precision estimation approach based on the residual error analysis. The residual error comprises three components, i. e. systematic error, model truncation error and random error. The approach realizes self-adaptive estimation of systematic errors in measurements following the theory of sparse representation of signals to minimize the low-frequency components of residual errors. By taking median- and high-frequency components as indexes, the spline model-approximation is improved by optimizing node sequence of the spline function and the weight selection for data fusion through iteration. Simulation has validated the performances of the proposed method.展开更多
In this study, numerical prediction of surges associated with a storm was made through the method of lines（MOL） in coordination with the newly proposed RKARMS（4, 4） method for the meghna estuarine region, along th...In this study, numerical prediction of surges associated with a storm was made through the method of lines（MOL） in coordination with the newly proposed RKARMS（4, 4） method for the meghna estuarine region, along the coast of Bangladesh. For this purpose, the vertically integrated shallow water equations（SWEs） in Cartesian coordinates were firstly transformed into ordinary differential equations（ODEs） of initial valued, which were then soloved using the new RKARMS（4, 4） method. Nested grid technique was employed for resolving the complexities of the region of interest with minimum cost. Fresh water discharge through the lower Meghna River was taken into account along the north east corner of the innermost child scheme. Numerical experiments were performed with the severe cyclone on April 1991 that crossed the coast over the study area. Simulated results by the study were found to be in good agreement with some reported data and were found to compare well with the results obtained by the MOL in addition with the classical 4th order Runge-Kutta（RK（4, 4）） method and the standard finite difference method（FDM）.展开更多
A new Runge-Kutta （PK） fourth order with four stages embedded method with error control is presentea m this paper for raster simulation in cellular neural network （CNN） environment. Through versatile algorithm, si...A new Runge-Kutta （PK） fourth order with four stages embedded method with error control is presentea m this paper for raster simulation in cellular neural network （CNN） environment. Through versatile algorithm, single layer/raster CNN array is implemented by incorporating the proposed technique. Simulation results have been obtained, and comparison has also been carried out to show the efficiency of the proposed numerical integration algorithm. The analytic expressions for local truncation error and global truncation error are derived. It is seen that the RK-embedded root mean square outperforms the RK-embedded Heronian mean and RK-embedded harmonic mean.展开更多
A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh rat...A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.展开更多
In this paper, a class of explicit difference schemes with parameters for solving five-dimensional heat-conduction equation are constructed and studied.the truncation error reaches O（τ^2＋ h%4）, and the stability c...In this paper, a class of explicit difference schemes with parameters for solving five-dimensional heat-conduction equation are constructed and studied.the truncation error reaches O（τ^2＋ h%4）, and the stability condition is given. Finally, the numerical examples and numerical results are presented to show the advantage of the schemes and the correctness of theoretical analysis.展开更多
In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We the...In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.展开更多
In this paper, an explicit three_level symmetrical differencing scheme with parameters for solving parabolic partial differential equation of three_dimension will be considered. The stability condition and local trunc...In this paper, an explicit three_level symmetrical differencing scheme with parameters for solving parabolic partial differential equation of three_dimension will be considered. The stability condition and local truncation error for the scheme are r<1/2 and O( Δ t 2+ Δ x 4+ Δ y 4+ Δ z 4) ,respectively.展开更多
基金We acknowledge the anonymous reviewers for their helpful comments and criticism on an earlier manuscript.The authors are indebted to the supports from the National Natural Science Foundation of China under Grant Nos.40175025and 40028504the State key Bas
文摘Highly accurate spatial discretization is essentially required to perform numerical climate and weather prediction. The difference between the differential and the finite-difference operator is however a primitive error source in the numerics. This paper presents an optimization of centered finite-difference operator based on the principle of constrained cost function, which can reduce the truncation error to minimum. In the optimization point of view, such optimal operator is in fact an attempt to minimize spatial truncation er-rors in atmospheric modeling, in a simple way and indeed a quite innovative way to implement Variational Continuous Assimilation (VCA) technique. Furthermore, the optimizing difference operator is consciously designed to be meshing-independent, so that it can be used for most Arakawa-mesh configurations, such as un-staggered (Arakawa-A) or com-monly staggered (Arakawa-B, Arakawa-C, Arakawa-D) mesh. But for the calibration purpose, the pro-posed operator is implemented on an un-staggered mesh in which the truncation oscillation is mostly ex-cited, and it thus makes a severe and indeed a benchmark test for the proposed optimal scheme. Both theo-retical investigation and practical modeling indicate that the aforementioned numerical noise can be significantly eliminated.
基金Projcct supported by the Natural Science Foundation of China (Grant No. 10371009 ) of Beijing Educational Committee (No. 2002KJ112).
文摘The truncation error associated with a given sampling representation is defined as the difference between the signal and on approximating sumutilizing a finite number of terms. In this paper we give uniform bound for truncation error of bandlimited functions in the n dimensional Lebesgue space Lp(Rn) associated with multidimensional Shannon sampling representation.
文摘A new so called truncation error reduction method (TERM) is developed in this work. This is an iterative process which uses a coarse grid (2 h ) to estimate the truncation error and then reduces the error on the original grid ( h ). The purpose is to use coarse grids to get more accurate results and to develop a new method which could do coarse grid direct numerical simulation (DNS) for more accurate and acceptable DNS solutions.
文摘This paper investigates double sampling series derivatives for bivariate functions defined on R2 that are in the Bernstein space. For this sampling series, we estimate some of the pointwise and uniform bounds when the function satisfies some decay conditions. The truncated series of this formula allow us to approximate any order of partial derivatives for function from Bernstein space using only a finite number of samples from the function itself. This sampling formula will be useful in the approximation theory and its applications, especially after having the truncation error well-established. Examples with tables and figures are given at the end of the paper to illustrate the advantages of this formula.
文摘The damage identification is made by the numerical simulation analysis of a five-storey-and-two-span RC frame structure, using improved and unimproved direct analytical method respectively; and the fundamental equations were solved by the minimal least square method （viz. general inverse method）. It demonstrates that the feasibility and the accuracy of the present approach were impoved significantly, compared with the result of unimproved damage identification.
基金Supported by the National Natural Science Foundation of China (10971251, 11101220 and 11271199)the Program for new century excellent talents in University of China (NCET-10-0513)
文摘Let B^pΩ, 1 ≤ p 〈 ∞, be the space of all bounded functions from Lp（R） which can be extended to entire functions of exponential type Ω. The uniform error bounds for truncated Whittaker-Kotelnikov-Shannon series based on local sampling are derived for functions f ∈ B^pΩ without decay assumption at infinity. Then the optimal bounds of the aliasing error and truncation error of Whittaker-Kotelnikov-Shannon expansion for non-bandlimited functions from Sobolev classes L/（Wp（R）） are determined up to a logarithmic factor.
基金Supported by the National Natural Science Foundation of China(Nos.61379014 and 11271199)
文摘Let B^pΩ,1≤Р≤∞,be the set of all bounded functions in L^p（R）which can be extended to entire unctions of exponential typeΩ. The unitbrm bounds for truncation error of Shannon sampling expansion fromlocal averages are obtained for functions f∈BpΩwith the decay condition f（t）≤A/t^δ,t≠0,where A and δare positive constants. Furthermore we also establish similar results for non-bandlimit functions in Besov classes with the same decay condition as above.
基金The author is thankful to D.d’Humi`eres for his parallel work on the Fourier analysis of the TRT AADE model and to anonymous referee for constructive suggestions.
文摘This paper establishes relations between the stability and the high-order truncated corrections for modeling of the mass conservation equation with the tworelaxation-times(TRT)collision operator.First we propose a simple method to derive the truncation errors from the exact,central-difference type,recurrence equations of the TRT scheme.They also supply its equivalent three-time-level discretization form.Two different relationships of the two relaxation rates nullify the third(advection)and fourth(pure diffusion)truncation errors,for any linear equilibrium and any velocity set.However,the two relaxation times alone cannot remove the leading-order advection-diffusion error,because of the intrinsic fourth-order numerical diffusion.The truncation analysis is carefully verified for the evolution of concentration waves with the anisotropic diffusion tensors.The anisotropic equilibrium functions are presented in a simple but general form,suitable for the minimal velocity sets and the d2Q9,d3Q13,d3Q15 and d3Q19 velocity sets.All anisotropic schemes are complemented by their exact necessary von Neumann stability conditions and equivalent finite-difference stencils.The sufficient stability conditions are proposed for the most stable(OTRT)family,which enables modeling at any Peclet numbers with the same velocity amplitude.The heuristic stability analysis of the fourth-order truncated corrections extends the optimal stability to larger relationships of the two relaxation rates,in agreementwith the exact(one-dimensional)and numerical(multi-dimensional)stability analysis.A special attention is put on the choice of the equilibrium weights.By combining accuracy and stability predictions,several strategies for selecting the relaxation and free-tunable equilibrium parameters are suggested and applied to the evolution of the Gaussian hill.
基金State Key Laboratory of Oil/Gas Reservoir Geology and Exploitation (PLN0115).
文摘A family of symmetric (hybrid) two step sixth P-stable methods for the accurate numerical integration of second order periodic initial value problems have been considered in this paper. These methods, which require only three (new) function evaluation per iteration and per step integration. These methods have minimal local truncation error (LTE) and smaller phase-lag of sixth order than some sixth orders P-stable methods in [1-3,10-11]. The theoretical and numerical results show that these methods in this paper are more accurate and efficient than some methods proposed in [1-3,10].
文摘This paper presents truncation errors among Corrector Formula for left Rectangular rule and Corrector Formula for middle Rectangular rule respectively. It also displays an analysis on convergence order of compound corrector formulas for rectangular rule. Examples of numerical calculation have validated theoretical analysis.
文摘The truncation error of improved Cotes formula is presented in this paper. It also displays an analysis on convergence order of improved Cotes formula. Examples of numerical calculation is given in the end.
文摘In this paper,we establish a new multivariate Hermite sampling series involving samples from the function itself and its mixed and non-mixed partial derivatives of arbitrary order.This multivariate form of Hermite sampling will be valid for some classes of multivariate entire functions,satisfying certain growth conditions.We will show that many known results included in Commun Korean Math Soc,2002,17:731-740,Turk J Math,2017,41:387-403 and Filomat,2020,34:3339-3347 are special cases of our results.Moreover,we estimate the truncation error of this sampling based on localized sampling without decay assumption.Illustrative examples are also presented.
文摘The solution of time-dependent hyperbolic conservation laws on cut cell meshes causes the small cell problem:standard schemes are not stable on the arbitrarily small cut cells if an explicit time stepping scheme is used and the time step size is chosen based on the size of the background cells.In May and Berger(J Sci Comput 71:919–943,2017),the mixed explicit-implicit approach in general and MUSCL-Trap(monotonic upwind scheme for conservation laws and trapezoidal scheme)in particular have been introduced to solve this problem by using implicit time stepping on the cut cells.Theoretical and numerical results have indicated that this might lead to a loss in accuracy when switching between the explicit and implicit time stepping.In this contribution,we examine this in more detail and will prove in one dimension that the specific combination MUSCL-Trap of an explicit second-order and an implicit second-order scheme results in a fully second-order mixed scheme.As this result is unlikely to hold in two dimensions,we also introduce two new versions of mixed explicit-implicit schemes based on exchanging the explicit scheme.We present numerical tests in two dimensions where we compare the new versions with the original MUSCL-Trap scheme.
基金National Natural Science Foundation of China(60604020)
文摘The classic state methods for trajectory estimation in boost phase with multi-range-rate system include method of point-by-point manner and that of spline-model-based manner. Both are deficient in terms of model-approximation accuracy and systematic error determination thus resulting in the estimation errors well beyond the requirements, especially, concerning the maneuvering trajectory. This article proposes a new high-precision estimation approach based on the residual error analysis. The residual error comprises three components, i. e. systematic error, model truncation error and random error. The approach realizes self-adaptive estimation of systematic errors in measurements following the theory of sparse representation of signals to minimize the low-frequency components of residual errors. By taking median- and high-frequency components as indexes, the spline model-approximation is improved by optimizing node sequence of the spline function and the weight selection for data fusion through iteration. Simulation has validated the performances of the proposed method.
文摘In this study, numerical prediction of surges associated with a storm was made through the method of lines（MOL） in coordination with the newly proposed RKARMS（4, 4） method for the meghna estuarine region, along the coast of Bangladesh. For this purpose, the vertically integrated shallow water equations（SWEs） in Cartesian coordinates were firstly transformed into ordinary differential equations（ODEs） of initial valued, which were then soloved using the new RKARMS（4, 4） method. Nested grid technique was employed for resolving the complexities of the region of interest with minimum cost. Fresh water discharge through the lower Meghna River was taken into account along the north east corner of the innermost child scheme. Numerical experiments were performed with the severe cyclone on April 1991 that crossed the coast over the study area. Simulated results by the study were found to be in good agreement with some reported data and were found to compare well with the results obtained by the MOL in addition with the classical 4th order Runge-Kutta（RK（4, 4）） method and the standard finite difference method（FDM）.
基金supported as a part of Technical Quality Improvement Programme (TEQIP)
文摘A new Runge-Kutta （PK） fourth order with four stages embedded method with error control is presentea m this paper for raster simulation in cellular neural network （CNN） environment. Through versatile algorithm, single layer/raster CNN array is implemented by incorporating the proposed technique. Simulation results have been obtained, and comparison has also been carried out to show the efficiency of the proposed numerical integration algorithm. The analytic expressions for local truncation error and global truncation error are derived. It is seen that the RK-embedded root mean square outperforms the RK-embedded Heronian mean and RK-embedded harmonic mean.
文摘A class of two-level explicit difference schemes are presented for solving three-dimensional heat conduction equation. When the order of truncation error is 0(Deltat + (Deltax)(2)), the stability condition is mesh ratio r = Deltat/(Deltax)(2) = Deltat/(Deltay)(2) = Deltat/(Deltaz)(2) less than or equal to 1/2, which is better than that of all the other explicit difference schemes. And when the order of truncation error is 0((Deltat)(2) + (Deltax)(4)), the stability condition is r less than or equal to 1/6, which contains the known results.
基金Supported by NSF of the Education Department of Henan Province(20031100010)
文摘In this paper, a class of explicit difference schemes with parameters for solving five-dimensional heat-conduction equation are constructed and studied.the truncation error reaches O（τ^2＋ h%4）, and the stability condition is given. Finally, the numerical examples and numerical results are presented to show the advantage of the schemes and the correctness of theoretical analysis.
基金The project supported by the China NKBRSF(2001CB409604)
文摘In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.
文摘In this paper, an explicit three_level symmetrical differencing scheme with parameters for solving parabolic partial differential equation of three_dimension will be considered. The stability condition and local truncation error for the scheme are r<1/2 and O( Δ t 2+ Δ x 4+ Δ y 4+ Δ z 4) ,respectively.
