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.展开更多
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.展开更多
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.展开更多
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 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.展开更多
Based on the improved mode superposition method proposed by Z. D. Ma and I.Hagiwara, a precisely compensated efficient mode synthesis method is developed. The calculationprocedure is discussed in detail and the trunca...Based on the improved mode superposition method proposed by Z. D. Ma and I.Hagiwara, a precisely compensated efficient mode synthesis method is developed. The calculationprocedure is discussed in detail and the truncation error is also analyzed. By comparison, it isshown that this method has a higher accuracy and a less calculation time than the general used ones.展开更多
In previous decades, many of the practical problems arising in scientific fields such as physics, engineering, and mathematics have been related to nonlinear fractional partial differential equations. One of these non...In previous decades, many of the practical problems arising in scientific fields such as physics, engineering, and mathematics have been related to nonlinear fractional partial differential equations. One of these nonlinear partial differential equations, the dissipative wave equation, has been found to have a plethora of useful applications in different fields. A special class of solutions has been studied for the dissipative wave equation including exact solutions and approximate solutions. The aim of this article is to compare the non-polynomial spline method and the cubic B-spline method with the solution of a nonlinear dissipative wave equation. We will conduct a comparison of the stability of the two methods using the Von Neumann stability analysis. In addition, a numerical example will be presented to illustrate the accuracy of these methods.展开更多
In this paper,a stochastic linear two-step scheme has been presented to approximate backward stochastic differential equations(BSDEs).A necessary and sufficient condition is given to judge the L 2-stability of our num...In this paper,a stochastic linear two-step scheme has been presented to approximate backward stochastic differential equations(BSDEs).A necessary and sufficient condition is given to judge the L 2-stability of our numerical schemes.This stochastic linear two-step method possesses a family of 3-order convergence schemes in the sense of strong stability.The coefficients in the numerical methods are inferred based on the constraints of strong stability and n-order accuracy(n∈N^(+)).Numerical experiments illustrate that the scheme is an efficient probabilistic numerical method.展开更多
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.展开更多
The Sobolev space H^(■)(R^(d)),where■>d/2,is an important function space that has many applications in various areas of research.Attributed to the inertia of a measurement instrument,it is desirable in sampling t...The Sobolev space H^(■)(R^(d)),where■>d/2,is an important function space that has many applications in various areas of research.Attributed to the inertia of a measurement instrument,it is desirable in sampling theory to recover a function by its nonuniform sampling.In the present paper,based on dual framelet systems for the Sobolev space pair(H^(s)(R^(d)),H^(-s)(R^(d))),where d/2<s<■,we investigate the problem of constructing the approximations to all the functions in H^(■)(R^(d))by nonuniform sampling.We first establish the convergence rate of the framelet series in(H^(s)(R^(d)),H^(-s)(R^(d))),and then construct the framelet approximation operator that acts on the entire space H^(■)(R^(d)).We examine the stability property for the framelet approximation operator with respect to the perturbations of shift parameters,and obtain an estimate bound for the perturbation error.Our result shows that under the condition d/2<s<■,the approximation operator is robust to shift perturbations.Motivated by Hamm(2015)’s work on nonuniform sampling and approximation in the Sobolev space,we do not require the perturbation sequence to be in■^(α)(Z^(d)).Our results allow us to establish the approximation for every function in H^(■)(R^(d))by nonuniform sampling.In particular,the approximation error is robust to the jittering of the samples.展开更多
基金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.
文摘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.
基金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.
文摘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.
基金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.
文摘Based on the improved mode superposition method proposed by Z. D. Ma and I.Hagiwara, a precisely compensated efficient mode synthesis method is developed. The calculationprocedure is discussed in detail and the truncation error is also analyzed. By comparison, it isshown that this method has a higher accuracy and a less calculation time than the general used ones.
文摘In previous decades, many of the practical problems arising in scientific fields such as physics, engineering, and mathematics have been related to nonlinear fractional partial differential equations. One of these nonlinear partial differential equations, the dissipative wave equation, has been found to have a plethora of useful applications in different fields. A special class of solutions has been studied for the dissipative wave equation including exact solutions and approximate solutions. The aim of this article is to compare the non-polynomial spline method and the cubic B-spline method with the solution of a nonlinear dissipative wave equation. We will conduct a comparison of the stability of the two methods using the Von Neumann stability analysis. In addition, a numerical example will be presented to illustrate the accuracy of these methods.
基金supported by the National Natural Science of China No.11971263,11871458Shandong Provincial Natural Science Foundation No.ZR2019ZD41National Key R&D Program of China No.2018YFA0703900。
文摘In this paper,a stochastic linear two-step scheme has been presented to approximate backward stochastic differential equations(BSDEs).A necessary and sufficient condition is given to judge the L 2-stability of our numerical schemes.This stochastic linear two-step method possesses a family of 3-order convergence schemes in the sense of strong stability.The coefficients in the numerical methods are inferred based on the constraints of strong stability and n-order accuracy(n∈N^(+)).Numerical experiments illustrate that the scheme is an efficient probabilistic numerical method.
基金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.
基金supported by National Natural Science Foundation of China(Grant Nos.61961003,61561006 and 11501132)Natural Science Foundation of Guangxi(Grant Nos.2018JJA110110 and 2016GXNSFAA380049)+1 种基金the talent project of Education Department of Guangxi Government for Young-Middle-Aged Backbone Teacherssupported by National Science Foundation of USA(Grant No.DMS-1712602)。
文摘The Sobolev space H^(■)(R^(d)),where■>d/2,is an important function space that has many applications in various areas of research.Attributed to the inertia of a measurement instrument,it is desirable in sampling theory to recover a function by its nonuniform sampling.In the present paper,based on dual framelet systems for the Sobolev space pair(H^(s)(R^(d)),H^(-s)(R^(d))),where d/2<s<■,we investigate the problem of constructing the approximations to all the functions in H^(■)(R^(d))by nonuniform sampling.We first establish the convergence rate of the framelet series in(H^(s)(R^(d)),H^(-s)(R^(d))),and then construct the framelet approximation operator that acts on the entire space H^(■)(R^(d)).We examine the stability property for the framelet approximation operator with respect to the perturbations of shift parameters,and obtain an estimate bound for the perturbation error.Our result shows that under the condition d/2<s<■,the approximation operator is robust to shift perturbations.Motivated by Hamm(2015)’s work on nonuniform sampling and approximation in the Sobolev space,we do not require the perturbation sequence to be in■^(α)(Z^(d)).Our results allow us to establish the approximation for every function in H^(■)(R^(d))by nonuniform sampling.In particular,the approximation error is robust to the jittering of the samples.