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.展开更多
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.展开更多
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].展开更多
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.展开更多
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.
Engineering and applied mathematics disciplines that involve differential equations in general,and initial value problems in particular,include classical mechanics,thermodynamics,electromagnetism,and the general theor...Engineering and applied mathematics disciplines that involve differential equations in general,and initial value problems in particular,include classical mechanics,thermodynamics,electromagnetism,and the general theory of relativity.A reliable,stable,efficient,and consistent numerical scheme is frequently required for modelling and simulation of a wide range of real-world problems using differential equations.In this study,the tangent slope is assumed to be the contra-harmonic mean,in which the arithmetic mean is used as a correction instead of Euler’s method to improve the efficiency of the improved Euler’s technique for solving ordinary differential equations with initial conditions.The stability,consistency,and efficiency of the system were evaluated,and the conclusions were supported by the presentation of numerical test applications in engineering.According to the stability analysis,the proposed method has a wider stability region than other well-known methods that are currently used in the literature for solving initial-value problems.To validate the rate convergence of the numerical technique,a few initial value problems of both scalar and vector valued types were examined.The proposed method,modified Euler explicit method,and other methods known in the literature have all been used to calculate the absolute maximum error,absolute error at the last grid point of the integration interval under consideration,and computational time in seconds to test the performance.The Lorentz system was used as an example to illustrate the validity of the solution provided by the newly developed method.The method is determined to be more reliable than the commonly existing methods with the same order of convergence,as mentioned in the literature for numerical calculations and visualization of the results produced by all the methods discussed,Mat Lab-R2011b has been used.展开更多
In this paper,three kinds of discrete formulae for the Caputo fractional derivative are studied,including the modified L1 discretisation forα∈(O,1),and L2 discretisation and L2C discretisation forα∈(1,2).The trunc...In this paper,three kinds of discrete formulae for the Caputo fractional derivative are studied,including the modified L1 discretisation forα∈(O,1),and L2 discretisation and L2C discretisation forα∈(1,2).The truncation error estimates and the properties of the coeffcients of all these discretisations are analysed in more detail.Finally,the theoretical analyses areverifiedby thenumerical examples.展开更多
The present article mainly focuses on the fractional derivatives with an exponential kernel(“exponential fractional derivatives”for brevity).First,several extended integral transforms of the exponential fractional d...The present article mainly focuses on the fractional derivatives with an exponential kernel(“exponential fractional derivatives”for brevity).First,several extended integral transforms of the exponential fractional derivatives are proposed,including the Fourier transform and the Laplace transform.Then,the L2 discretisation for the exponential Caputo derivative with a∈(1,2)is established.The estimation of the truncation error and the properties of the coefficients are discussed.In addition,a numerical example is given to verify the correctness of the derived L2 discrete formula.展开更多
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.展开更多
基金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.
基金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.
文摘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].
文摘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.
文摘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.
文摘Engineering and applied mathematics disciplines that involve differential equations in general,and initial value problems in particular,include classical mechanics,thermodynamics,electromagnetism,and the general theory of relativity.A reliable,stable,efficient,and consistent numerical scheme is frequently required for modelling and simulation of a wide range of real-world problems using differential equations.In this study,the tangent slope is assumed to be the contra-harmonic mean,in which the arithmetic mean is used as a correction instead of Euler’s method to improve the efficiency of the improved Euler’s technique for solving ordinary differential equations with initial conditions.The stability,consistency,and efficiency of the system were evaluated,and the conclusions were supported by the presentation of numerical test applications in engineering.According to the stability analysis,the proposed method has a wider stability region than other well-known methods that are currently used in the literature for solving initial-value problems.To validate the rate convergence of the numerical technique,a few initial value problems of both scalar and vector valued types were examined.The proposed method,modified Euler explicit method,and other methods known in the literature have all been used to calculate the absolute maximum error,absolute error at the last grid point of the integration interval under consideration,and computational time in seconds to test the performance.The Lorentz system was used as an example to illustrate the validity of the solution provided by the newly developed method.The method is determined to be more reliable than the commonly existing methods with the same order of convergence,as mentioned in the literature for numerical calculations and visualization of the results produced by all the methods discussed,Mat Lab-R2011b has been used.
文摘In this paper,three kinds of discrete formulae for the Caputo fractional derivative are studied,including the modified L1 discretisation forα∈(O,1),and L2 discretisation and L2C discretisation forα∈(1,2).The truncation error estimates and the properties of the coeffcients of all these discretisations are analysed in more detail.Finally,the theoretical analyses areverifiedby thenumerical examples.
文摘The present article mainly focuses on the fractional derivatives with an exponential kernel(“exponential fractional derivatives”for brevity).First,several extended integral transforms of the exponential fractional derivatives are proposed,including the Fourier transform and the Laplace transform.Then,the L2 discretisation for the exponential Caputo derivative with a∈(1,2)is established.The estimation of the truncation error and the properties of the coefficients are discussed.In addition,a numerical example is given to verify the correctness of the derived L2 discrete formula.
文摘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.