This paper deals with numerical computation and analysis for the initial boundary problems of two dimensional(2D)Sobolev equations with piecewise continuous argument.Firstly,a two-level high-order compact difference m...This paper deals with numerical computation and analysis for the initial boundary problems of two dimensional(2D)Sobolev equations with piecewise continuous argument.Firstly,a two-level high-order compact difference method(HOCDM)with computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4))is suggested,whereτ,h_(x),h_(y) denote the temporal and spatial stepsizes of the method,respectively.In order to improve the temporal computational accuracy of this method,the Richardson extrapolation technique is used and thus a new two-level HOCDMis derived,which is proved to be convergent of order four both in time and space.Although the new two-level HOCDM has the higher computational accuracy in time than the previous one,it will bring a larger computational cost.To overcome this deficiency,a three-level HOCDM with computational accuracy O(τ^(4)+h_(x)^(4)+h_(y)^(4))is constructed.Finally,with a series of numerical experiments,the theoretical accuracy and computational efficiency of the above methods are further verified.展开更多
This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called...This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4)),where τ,h_(x )and h_(y) are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.展开更多
To improve the spectral characteristics of the high-order weighted compact nonlinear scheme(WCNS),optimized flux difference schemes are proposed.The disadvantages in previous optimization routines,i.e.,reducing formal...To improve the spectral characteristics of the high-order weighted compact nonlinear scheme(WCNS),optimized flux difference schemes are proposed.The disadvantages in previous optimization routines,i.e.,reducing formal orders,or extending stencil widths,are avoided in the new optimized schemes by utilizing fluxes from both cell-edges and cell-nodes.Optimizations are implemented with Fourier analysis for linear schemes and the approximate dispersion relation(ADR)for nonlinear schemes.Classical difference schemes are restored near discontinuities to suppress numerical oscillations with use of a shock sensor based on smoothness indicators.The results of several benchmark numerical tests indicate that the new optimized difference schemes outperform the classical schemes,in terms of accuracy and resolution for smooth wave and vortex,especially for long-time simulations.Using optimized schemes increases the total CPU time by less than 4%.展开更多
In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered gri...In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition （ABC） to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.展开更多
The Spalart-Allmaras （S-A） turbulence model, the shear-stress transport （SST） turbulence model and their compressibility corrections are revaluated for hypersonic compression comer flows by using high-order differ...The Spalart-Allmaras （S-A） turbulence model, the shear-stress transport （SST） turbulence model and their compressibility corrections are revaluated for hypersonic compression comer flows by using high-order difference schemes. The compressibility effect of density gradient, pressure dilatation and turbulent Mach number is accounted. In order to reduce confusions between model uncertainties and discretization errors, the formally fifth-order explicit weighted compact nonlinear scheme （WCNS-E-5） is adopted for convection terms, and a fourth-order staggered central difference scheme is applied for viscous terms. The 15° and 34° compression comers at Mach number 9.22 are investigated. Numerical results show that the original SST model is superior to the original S-A model in the resolution of separated regions and predictions of wall pressures and wall heat-flux rates. The capability of the S-A model can be largely improved by blending Catris＇ and Shur＇s compressibility corrections. Among the three corrections of the SST model listed in the present paper, Catris＇ modification brings the best results. However, the dissipation and pressure dilatation corrections result in much larger separated regions than that of the experiment, and are much worse than the original SST model as well as the other two corrections. The correction of turbulent Mach number makes the separated region slightly smaller than that of the original SST model. Some results of low-order schemes are also presented. When compared to the results of the high-order schemes, the separated regions are smaller, and the peak wall pressures and peak heat-flux rates are lower in the region of the reattachment points.展开更多
Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with ...Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with increased order of accuracy. Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note that there exist some very small coefficients. With the order increasing, the number of these small coefficients increases, however, the values decrease sharply. An error analysis demonstrates that omitting these small coefficients not only maintain approximately the same level of accuracy of finite difference but also reduce computational cost significantly. Moreover, it is easier to truncate for the high-order finite-difference formulas than for the pseudospectral for- mulas. Thus this study proposes a truncated high-order finite-difference method, and then demonstrates the efficiency and applicability of the method with some numerical examples.展开更多
A high-order accuracy explicit difference scheme for solving 4-dimensional heatconduction equation is constructed. The stability condition is r = △t/△x^2 = △t/△y^2 = △t/△z^2 = △t/△w^2 〈 3/8, and the truncatio...A high-order accuracy explicit difference scheme for solving 4-dimensional heatconduction equation is constructed. The stability condition is r = △t/△x^2 = △t/△y^2 = △t/△z^2 = △t/△w^2 〈 3/8, and the truncation error is O（△t^2 ＋ △x^4）.展开更多
The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass cons...The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass conservation, this model is not satisfactory during the simulation computations. In this paper, a mass correction is introduced into the ZSC model to make up the mass leakage, while a high-order difference is used to calculate the gradient of the order parameter to improve the accuracy. To verify the improved model, several three-dimensional multiphase flow simulations are carried out,including a bubble in a stationary flow, the merging of two bubbles, and the bubble rising under buoyancy. The numerical simulations show that the results from the present model are in good agreement with those from previous experiments and simulations. The present model not only retains the good properties of the original ZSC model, but also achieves the mass conservation and higher accuracy.展开更多
Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, th...Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, the properties of high-order finite difference schemes are studied based on idealized numerical testing, for the purpose of their application in the Global/Regional Assimilation and Prediction System(GRAPES) model. It is found that the pros and cons due to grid staggering choices diminish with higher-order schemes based on linearized analysis of the one-dimensional gravity wave equation. The improvement of higher-order difference schemes is still obvious for the mesh with smooth varied grid distance. The results of discontinuous square wave testing also exhibits the superiority of high-order schemes. For a model grid with severe non-uniformity and non-orthogonality, the advantage of high-order difference schemes is inapparent, as shown by the results of two-dimensional idealized advection tests under a terrain-following coordinate. In addition, the increase in computational expense caused by high-order schemes can be avoided by the precondition technique used in the GRAPES model. In general, a high-order finite difference scheme is a preferable choice for the tropical regional GRAPES model with a quasi-uniform and quasi-orthogonal grid mesh.展开更多
A family of high-order accuracy explict difference schemes for solving 3-dimension parabolic P. D. E. is constructed. The stability condition is r = Deltat/Deltax(2) Deltat/Deltay(2) = Deltat/Deltaz(2) < 1/2 ,and t...A family of high-order accuracy explict difference schemes for solving 3-dimension parabolic P. D. E. is constructed. The stability condition is r = Deltat/Deltax(2) Deltat/Deltay(2) = Deltat/Deltaz(2) < 1/2 ,and the truncation error is 0(<Delta>t(2) + Deltax(4)).展开更多
In this paper, a new three-level explicit difference scheme with high-order accuracy is proposed for solving three-dimensional parabolic equations. The stability condition is r = Delta t/Delta x(2) = Delta t/Delta gam...In this paper, a new three-level explicit difference scheme with high-order accuracy is proposed for solving three-dimensional parabolic equations. The stability condition is r = Delta t/Delta x(2) = Delta t/Delta gamma(2) = Delta t/Delta z(2) less than or equal to 1/4, and the truncation error is O(Delta t(2) + Delta x(4)).展开更多
We find a lower bound for the essential norm of the difference of two composition operators acting on H 2（BN ） or As2 （BN ） （s 1）. This result plays an important role in proving a necessary and sufficient condit...We find a lower bound for the essential norm of the difference of two composition operators acting on H 2（BN ） or As2 （BN ） （s 1）. This result plays an important role in proving a necessary and sufficient condition for the difference of linear fractional composition operators to be compact, which answers a question posed by MacCluer and Weir in 2005.展开更多
A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite differen...A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.展开更多
In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws...In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.展开更多
In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstructio...In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4 est based on the forest of octrees is adopted. Several two-and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.展开更多
In the past, several authors studied spaces of m-th order difference sequences, among them, H.Polat and F.Basar （[17]） defined the Euler spaces of m-th order difference sequences e r 0 （△ （ m ） ）, e r c （△ （...In the past, several authors studied spaces of m-th order difference sequences, among them, H.Polat and F.Basar （[17]） defined the Euler spaces of m-th order difference sequences e r 0 （△ （ m ） ）, e r c （△ （ m ） ） and e r ∞ （△ （ m ） ） and characterized some classes of matrix transformations on them. In our paper, we add a new supplementary aspect to their research by characterizing classes of compact operators on those spaces. For that purpose, the spaces are treated as the matrix domains of a triangle in the classical sequence spaces c 0 , c and ∞ . The main tool for our characterizations is the Hausdorff measure of noncompactness.展开更多
A new method was proposed for constructing total variation diminishing （TVD） upwind schemes in conservation forms. Two limiters were used to prevent nonphysical oscillations across discontinuity. Both limiters can e...A new method was proposed for constructing total variation diminishing （TVD） upwind schemes in conservation forms. Two limiters were used to prevent nonphysical oscillations across discontinuity. Both limiters can ensure the nonlinear compact schemes TVD property. Two compact TVD （CTVD） schemes were tested, one is thirdorder accuracy, and the other is fifth-order. The performance of the numerical algorithms was assessed by one-dimensional complex waves and Riemann problems, as well as a twodimensional shock-vortex interaction and a shock-boundary flow interaction. Numerical results show their high-order accuracy and high resolution, and low oscillations across discontinuities.展开更多
Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to ...Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to eliminate the errors, a method of rotated staggered-grid finite-difference scheme（RSG） is proposed. However, the RSG brings serious numerical dispersion. The compact staggered-grid finite-difference scheme（CSG） is an implicit difference scheme, which use fewer grid points to suppress dispersion more effectively than the SSG. This paper combines the CSG with the RSG to derive a rotated staggered-grid compact finite-difference scheme（RSGC）. The numerical experiments indicate that the RSGC has weaker numerical dispersion and better accuracy than the RSG.展开更多
A global seventh-order dissipative compact finite-difference scheme is optimized in terms of time stability. The dissipative parameters appearing in the boundary closures are assumed to be different, resulting in an o...A global seventh-order dissipative compact finite-difference scheme is optimized in terms of time stability. The dissipative parameters appearing in the boundary closures are assumed to be different, resulting in an optimization problem with several parameters determined by applying a generic algorithm. The optimized schemes are analyzed carefully from the aspects of the eigenvalue distribution, the ε-pseudospectra, the short time behavior, and the Fourier analysis. Numerical experiments for the Euler equations are used to show the effectiveness of the final recommended scheme.展开更多
The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechan...The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechanics （English Edition）, 2007, 28（7）, 943-953, has the same performance as the conventional finite difference schemes. It is just another expression of the conventional finite difference schemes. The proposed expression does not have the advantages of a compact difference scheme. Nonetheless, we can more easily obtain and implement compared with the conventional expression in which the coefficients can only be obtained by solving equations, especially for higher accurate schemes.展开更多
文摘This paper deals with numerical computation and analysis for the initial boundary problems of two dimensional(2D)Sobolev equations with piecewise continuous argument.Firstly,a two-level high-order compact difference method(HOCDM)with computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4))is suggested,whereτ,h_(x),h_(y) denote the temporal and spatial stepsizes of the method,respectively.In order to improve the temporal computational accuracy of this method,the Richardson extrapolation technique is used and thus a new two-level HOCDMis derived,which is proved to be convergent of order four both in time and space.Although the new two-level HOCDM has the higher computational accuracy in time than the previous one,it will bring a larger computational cost.To overcome this deficiency,a three-level HOCDM with computational accuracy O(τ^(4)+h_(x)^(4)+h_(y)^(4))is constructed.Finally,with a series of numerical experiments,the theoretical accuracy and computational efficiency of the above methods are further verified.
文摘This paper deals with the numerical solutions of two-dimensional(2D)semi-linear reaction-diffusion equations(SLRDEs)with piecewise continuous argument(PCA)in reaction term.A high-order compact difference method called Ⅰ-type basic scheme is developed for solving the equations and it is proved under the suitable conditions that this method has the computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4)),where τ,h_(x )and h_(y) are the calculation stepsizes of the method in t-,x-and y-direction,respectively.With the above method and Newton linearized technique,a Ⅱ-type basic scheme is also suggested.Based on the both basic schemes,the corresponding Ⅰ-and Ⅱ-type alternating direction implicit(ADI)schemes are derived.Finally,with a series of numerical experiments,the computational accuracy and efficiency of the four numerical schemes are further illustrated.
基金Project supported by the National Key Project(No.GJXM92579)the Defense Industrial Technology Development Program(No.C1520110002)the State Administration of Science,Technology and Industry for National Defence,China。
文摘To improve the spectral characteristics of the high-order weighted compact nonlinear scheme(WCNS),optimized flux difference schemes are proposed.The disadvantages in previous optimization routines,i.e.,reducing formal orders,or extending stencil widths,are avoided in the new optimized schemes by utilizing fluxes from both cell-edges and cell-nodes.Optimizations are implemented with Fourier analysis for linear schemes and the approximate dispersion relation(ADR)for nonlinear schemes.Classical difference schemes are restored near discontinuities to suppress numerical oscillations with use of a shock sensor based on smoothness indicators.The results of several benchmark numerical tests indicate that the new optimized difference schemes outperform the classical schemes,in terms of accuracy and resolution for smooth wave and vortex,especially for long-time simulations.Using optimized schemes increases the total CPU time by less than 4%.
基金supported by the National Natural Science Foundation of China(NSFC)(Grant No. 41074100)the Program for New Century Excellent Talents in University of Ministry of Education of China(Grant No. NCET-10-0812)
文摘In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition （ABC） to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.
基金Foundation items： National Basic Research Program of China （2009CB723801） National Natural Science Foundation of China （11072259）
文摘The Spalart-Allmaras （S-A） turbulence model, the shear-stress transport （SST） turbulence model and their compressibility corrections are revaluated for hypersonic compression comer flows by using high-order difference schemes. The compressibility effect of density gradient, pressure dilatation and turbulent Mach number is accounted. In order to reduce confusions between model uncertainties and discretization errors, the formally fifth-order explicit weighted compact nonlinear scheme （WCNS-E-5） is adopted for convection terms, and a fourth-order staggered central difference scheme is applied for viscous terms. The 15° and 34° compression comers at Mach number 9.22 are investigated. Numerical results show that the original SST model is superior to the original S-A model in the resolution of separated regions and predictions of wall pressures and wall heat-flux rates. The capability of the S-A model can be largely improved by blending Catris＇ and Shur＇s compressibility corrections. Among the three corrections of the SST model listed in the present paper, Catris＇ modification brings the best results. However, the dissipation and pressure dilatation corrections result in much larger separated regions than that of the experiment, and are much worse than the original SST model as well as the other two corrections. The correction of turbulent Mach number makes the separated region slightly smaller than that of the original SST model. Some results of low-order schemes are also presented. When compared to the results of the high-order schemes, the separated regions are smaller, and the peak wall pressures and peak heat-flux rates are lower in the region of the reattachment points.
基金supported by China Scholarship Council and partially by the National "863" Program of China under contract No. 2007AA06Z218.
文摘Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with increased order of accuracy. Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note that there exist some very small coefficients. With the order increasing, the number of these small coefficients increases, however, the values decrease sharply. An error analysis demonstrates that omitting these small coefficients not only maintain approximately the same level of accuracy of finite difference but also reduce computational cost significantly. Moreover, it is easier to truncate for the high-order finite-difference formulas than for the pseudospectral for- mulas. Thus this study proposes a truncated high-order finite-difference method, and then demonstrates the efficiency and applicability of the method with some numerical examples.
基金NSF of the Education Department of Henan Province(20031100010)
文摘A high-order accuracy explicit difference scheme for solving 4-dimensional heatconduction equation is constructed. The stability condition is r = △t/△x^2 = △t/△y^2 = △t/△z^2 = △t/△w^2 〈 3/8, and the truncation error is O（△t^2 ＋ △x^4）.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11862003 and 81860635)the Key Project of the Natural Science Foundation of Guangxi Zhuang Autonomous Region,China(Grant No.2017GXNSFDA198038)+3 种基金the Project of Natural Science Foundation of Guangxi Zhuang Autonomous Region,China(Grant No.2018GXNSFAA281302)the Project for Promotion of Young and Middle-aged Teachers’Basic Scientific Research Ability in Guangxi Universities,China(Grant No.2019KY0084)the“Bagui Scholar”Teams for Innovation and Research Project of Guangxi Zhuang Autonomous Region,Chinathe Graduate Innovation Program of Guangxi Normal University,China(Grant No.JXYJSKT-2019-007)。
文摘The Z–S–C multiphase lattice Boltzmann model [Zheng, Shu, and Chew(ZSC), J. Comput. Phys. 218, 353(2006)]is favored due to its good stability, high efficiency, and large density ratio. However, in terms of mass conservation, this model is not satisfactory during the simulation computations. In this paper, a mass correction is introduced into the ZSC model to make up the mass leakage, while a high-order difference is used to calculate the gradient of the order parameter to improve the accuracy. To verify the improved model, several three-dimensional multiphase flow simulations are carried out,including a bubble in a stationary flow, the merging of two bubbles, and the bubble rising under buoyancy. The numerical simulations show that the results from the present model are in good agreement with those from previous experiments and simulations. The present model not only retains the good properties of the original ZSC model, but also achieves the mass conservation and higher accuracy.
基金supported by the National Natural Science Foundation of China (Grant No. U1811464)。
文摘Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, the properties of high-order finite difference schemes are studied based on idealized numerical testing, for the purpose of their application in the Global/Regional Assimilation and Prediction System(GRAPES) model. It is found that the pros and cons due to grid staggering choices diminish with higher-order schemes based on linearized analysis of the one-dimensional gravity wave equation. The improvement of higher-order difference schemes is still obvious for the mesh with smooth varied grid distance. The results of discontinuous square wave testing also exhibits the superiority of high-order schemes. For a model grid with severe non-uniformity and non-orthogonality, the advantage of high-order difference schemes is inapparent, as shown by the results of two-dimensional idealized advection tests under a terrain-following coordinate. In addition, the increase in computational expense caused by high-order schemes can be avoided by the precondition technique used in the GRAPES model. In general, a high-order finite difference scheme is a preferable choice for the tropical regional GRAPES model with a quasi-uniform and quasi-orthogonal grid mesh.
文摘A family of high-order accuracy explict difference schemes for solving 3-dimension parabolic P. D. E. is constructed. The stability condition is r = Deltat/Deltax(2) Deltat/Deltay(2) = Deltat/Deltaz(2) < 1/2 ,and the truncation error is 0(<Delta>t(2) + Deltax(4)).
文摘In this paper, a new three-level explicit difference scheme with high-order accuracy is proposed for solving three-dimensional parabolic equations. The stability condition is r = Delta t/Delta x(2) = Delta t/Delta gamma(2) = Delta t/Delta z(2) less than or equal to 1/4, and the truncation error is O(Delta t(2) + Delta x(4)).
基金Supported by the National Natural Science Foundation of China (10971219)Shanghai Education Research and Innovation Project (10YZ185)Shanghai University Research Special Foundation for Outstanding Young Teachers (sjr09015)
文摘We find a lower bound for the essential norm of the difference of two composition operators acting on H 2（BN ） or As2 （BN ） （s 1）. This result plays an important role in proving a necessary and sufficient condition for the difference of linear fractional composition operators to be compact, which answers a question posed by MacCluer and Weir in 2005.
基金the National Natural Science Foundation of China
文摘A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.
基金Project supported by the National Natural Science Foundation of China(No.11571366)the Basic Research Foundation of National Numerical Wind Tunnel Project(No.NNW2018-ZT4A08)
文摘In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.
基金supported by the National Natural Science Foundation of China(Nos.91752114 and 11672160)
文摘In the present paper, high-order finite volume schemes on unstructured grids developed in our previous papers are extended to solve three-dimensional inviscid and viscous flows. The highorder variational reconstruction technique in terms of compact stencil is improved to reduce local condition numbers. To further improve the efficiency of computation, the adaptive mesh refinement technique is implemented in the framework of high-order finite volume methods. Mesh refinement and coarsening criteria are chosen to be the indicators for certain flow structures. One important challenge of the adaptive mesh refinement technique on unstructured grids is the dynamic load balancing in parallel computation. To solve this problem, the open-source library p4 est based on the forest of octrees is adopted. Several two-and three-dimensional test cases are computed to verify the accuracy and robustness of the proposed numerical schemes.
基金supported by the research project#144003 of the Serbian Ministry of Science, Technology and Development
文摘In the past, several authors studied spaces of m-th order difference sequences, among them, H.Polat and F.Basar （[17]） defined the Euler spaces of m-th order difference sequences e r 0 （△ （ m ） ）, e r c （△ （ m ） ） and e r ∞ （△ （ m ） ） and characterized some classes of matrix transformations on them. In our paper, we add a new supplementary aspect to their research by characterizing classes of compact operators on those spaces. For that purpose, the spaces are treated as the matrix domains of a triangle in the classical sequence spaces c 0 , c and ∞ . The main tool for our characterizations is the Hausdorff measure of noncompactness.
基金Project supported by the National Natural Science Foundation of China (Nos. 10172015 and 90205010)
文摘A new method was proposed for constructing total variation diminishing （TVD） upwind schemes in conservation forms. Two limiters were used to prevent nonphysical oscillations across discontinuity. Both limiters can ensure the nonlinear compact schemes TVD property. Two compact TVD （CTVD） schemes were tested, one is thirdorder accuracy, and the other is fifth-order. The performance of the numerical algorithms was assessed by one-dimensional complex waves and Riemann problems, as well as a twodimensional shock-vortex interaction and a shock-boundary flow interaction. Numerical results show their high-order accuracy and high resolution, and low oscillations across discontinuities.
文摘Numerical simulation in transverse isotropic media with tilted symmetry axis（TTI） using the standard staggered-grid finite-difference scheme（SSG）results in errors caused by averaging or interpolation. In order to eliminate the errors, a method of rotated staggered-grid finite-difference scheme（RSG） is proposed. However, the RSG brings serious numerical dispersion. The compact staggered-grid finite-difference scheme（CSG） is an implicit difference scheme, which use fewer grid points to suppress dispersion more effectively than the SSG. This paper combines the CSG with the RSG to derive a rotated staggered-grid compact finite-difference scheme（RSGC）. The numerical experiments indicate that the RSGC has weaker numerical dispersion and better accuracy than the RSG.
基金Project supported by the National Natural Science Foundation of China(Nos.11601517,11502296,61772542,and 61561146395)the Basic Research Foundation of National University of Defense Technology(No.ZDYYJ-CYJ20140101)
文摘A global seventh-order dissipative compact finite-difference scheme is optimized in terms of time stability. The dissipative parameters appearing in the boundary closures are assumed to be different, resulting in an optimization problem with several parameters determined by applying a generic algorithm. The optimized schemes are analyzed carefully from the aspects of the eigenvalue distribution, the ε-pseudospectra, the short time behavior, and the Fourier analysis. Numerical experiments for the Euler equations are used to show the effectiveness of the final recommended scheme.
基金Supported by the National Natural Science Foundation of China (Nos.50876114 and 10602043)the Program for New Century Excellent Talents in University,and the Scientific Research Key Project Fund of Ministry of Education (No.106142)
文摘The explicit compact difference scheme, proposed in Three-point explicit compact difference scheme with arbitrary order of accuracy and its application in CFD by Lin et al., published in Applied Mathematics and Mechanics （English Edition）, 2007, 28（7）, 943-953, has the same performance as the conventional finite difference schemes. It is just another expression of the conventional finite difference schemes. The proposed expression does not have the advantages of a compact difference scheme. Nonetheless, we can more easily obtain and implement compared with the conventional expression in which the coefficients can only be obtained by solving equations, especially for higher accurate schemes.
