Within the projection schemes for the incompressible Navier-Stokes equations(namely"pressure-correction"method),we consider the simplest method(of order one in time)which takes into account the pressure in b...Within the projection schemes for the incompressible Navier-Stokes equations(namely"pressure-correction"method),we consider the simplest method(of order one in time)which takes into account the pressure in both steps of the splitting scheme.For this scheme,we construct,analyze and implement a new high order compact spatial approximation on nonstaggered grids.This approach yields a fourth order accuracy in space with an optimal treatment of the boundary conditions(without error on the velocity)which could be extended to more general splitting.We prove the unconditional stability of the associated Cauchy problem via von Neumann analysis.Then we carry out a normal mode analysis so as to obtain more precise results about the behavior of the numerical solutions.Finally we present detailed numerical tests for the Stokes and the Navier-Stokes equations(including the driven cavity benchmark)to illustrate the theoretical results.展开更多
In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids....In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids. Firstly, the original equation is transformed from the physical domain (with a nonuniform mesh) to the computational domain (with a uniform mesh) by using a coordinate transformation. Then, a fourth order compact difference scheme is proposed to solve the transformed elliptic equation on uniform girds. After that, a multigrid method is employed to solve the linear algebraic system arising from the difference equation. At last, the numerical experiments on some elliptic problems with interior/boundary layers are conducted to show high accuracy and high efficiency of the present method.展开更多
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.展开更多
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.展开更多
To ensure time stability of a seventh-order dissipative compact finite difference scheme, fourth-order boundary closures are used near domain boundaries previously. However, this would reduce the global convergence ra...To ensure time stability of a seventh-order dissipative compact finite difference scheme, fourth-order boundary closures are used near domain boundaries previously. However, this would reduce the global convergence rate to fifth-order only. In this paper, we elevate the boundary closures to sixth-order to achieve seventh-order global accuracy. To keep the improved scheme time stable, the simultaneous approximation terms (SATs) are used to impose boundary conditions weakly. Eigenvalue analysis shows that the improved scheme is time stable. Numerical experiments for linear advection equations and one-dimensional Euler equations are implemented to validate the new scheme.展开更多
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.展开更多
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.展开更多
We have successfully attempted to solve the equations of full-MHD model within the framework of Ψ- ωformulation with an objective to evaluate the performance of a new higher order scheme to predict b...We have successfully attempted to solve the equations of full-MHD model within the framework of Ψ- ωformulation with an objective to evaluate the performance of a new higher order scheme to predict better values of control parameters of the flow. In particular for MHD flows, magnetic field and electrical conductivity are the control parameters. In this work, the results from our efficient high order accurate scheme are compared with the results of second order method and significant discrepancies are noted in separation length, drag coefficient and mean Nusselt number. The governing Navier-Stokes equation is fully nonlinear due to its coupling with Maxwell’s equations. The momentum equation has several highly nonlinear body-force terms due to full-MHD model in cylindrical polar system. Our high accuracy results predict that a relatively lower magnetic field is sufficient to achieve full suppression of boundary layer and this is a favorable result for practical applications. The present computational scheme predicts that a drag-coefficient minimum can be achieved when β=0.4 which is much lower when compared to the value β=1 as given by second order method. For a special value of β=0.65, it is found that the heat transfer rate is independent of electrical conductivity of the fluid. From the numerical values of physical quantities, we establish that the order of accuracy of the computed numerical results is fourth order accurate by using the method of divided differences.展开更多
An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in ...An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.展开更多
A new class of finite difference schemes--the weighted compact schemes are proposed. According to the idea of the WENO schemes, the weighted compact scheme is constructed by a combination of the approximations of deri...A new class of finite difference schemes--the weighted compact schemes are proposed. According to the idea of the WENO schemes, the weighted compact scheme is constructed by a combination of the approximations of derivatives on candidate stencils with properly assigned weights so that the non oscillatory property is achieved when discontinuities appear. The primitive function reconstruction method of ENO schemes is applied to obtain the conservative form of the weighted compact scheme. This new scheme not only preserves the characteristic of standard compact schemes and achieves high order accuracy and high resolution using a compact stencil, but also can accurately capture shock waves and discontinuities without oscillation. Numerical examples show that the new scheme is very promising and successful.展开更多
By using Richardson extrapolation and fourth-order compact finite difference scheme on different scale grids, a sixth-order solution is computed on the coarse grid. Other three techniques are applied to obtain a sixth...By using Richardson extrapolation and fourth-order compact finite difference scheme on different scale grids, a sixth-order solution is computed on the coarse grid. Other three techniques are applied to obtain a sixth-order solution on the fine grid, and thus give out three kinds of Richardson extrapolation-based sixth order compact computation methods. By carefully analyzing the truncation errors respectively on 2D Poisson equation, we compare the accuracy of these three sixth order methods theoretically. Numerical results for two test problems are discussed.展开更多
In this paper, we introduce high-order finite volume methods for the multi-term time fractional sub-diffusion equation. The time fractional derivatives are described in Caputo’s sense. By using some operators, we obt...In this paper, we introduce high-order finite volume methods for the multi-term time fractional sub-diffusion equation. The time fractional derivatives are described in Caputo’s sense. By using some operators, we obtain the compact finite volume scheme have high order accuracy. We use a compact operator to deal with spatial direction;then we can get the compact finite volume scheme. It is proved that the finite volume scheme is unconditionally stable and convergent in L<sub>∞</sub>-norm. The convergence order is O(τ<sup>2-α</sup> + h<sup>4</sup>). Finally, two numerical examples are given to confirm the theoretical results. Some tables listed also can explain the stability and convergence of the scheme.展开更多
In this paper, a novel multisymplectic scheme is proposed for the coupled nonlinear Schrodinger-KdV （CNLS-KdV） equations. The CNLS-KdV equations are rewritten into the multisymplectic Hamiltonian form by introducing...In this paper, a novel multisymplectic scheme is proposed for the coupled nonlinear Schrodinger-KdV （CNLS-KdV） equations. The CNLS-KdV equations are rewritten into the multisymplectic Hamiltonian form by introducing some canonical momenta. To simulate the problem efficiently, the CNLS-KdV equations are approximated by a high order compact method in space which preserves N semi-discrete multisymplectic conservation laws. We then discretize the semi-discrete system by using a symplectic midpoint scheme in time. Thus, a full-discrete multisymplectic scheme is obtained for the CNLS-KdV equations. The conservation laws of the full-discrete scheme are analyzed. Some numerical experiments are presented to further verify the convergence and conservation laws of the new scheme.展开更多
A high-order finite difference Pade scheme also called compact scheme for solving Korteweg-de Vries （KdV） equations, which preserve energy and mass conservations, was developed in this paper. This structure-preservi...A high-order finite difference Pade scheme also called compact scheme for solving Korteweg-de Vries （KdV） equations, which preserve energy and mass conservations, was developed in this paper. This structure-preserving algorithm has been widely applied in these years for its advantage of maintaining the inherited properties. For spatial discretization, the authors obtained an implicit compact scheme by which spatial derivative terms may be approximated through combining a few knots. By some numerical examples including propagation of single soliton and interaction of two solitons, the scheme is proved to be effective.展开更多
文摘Within the projection schemes for the incompressible Navier-Stokes equations(namely"pressure-correction"method),we consider the simplest method(of order one in time)which takes into account the pressure in both steps of the splitting scheme.For this scheme,we construct,analyze and implement a new high order compact spatial approximation on nonstaggered grids.This approach yields a fourth order accuracy in space with an optimal treatment of the boundary conditions(without error on the velocity)which could be extended to more general splitting.We prove the unconditional stability of the associated Cauchy problem via von Neumann analysis.Then we carry out a normal mode analysis so as to obtain more precise results about the behavior of the numerical solutions.Finally we present detailed numerical tests for the Stokes and the Navier-Stokes equations(including the driven cavity benchmark)to illustrate the theoretical results.
文摘In this paper, a high order compact difference scheme and a multigrid method are proposed for solving two-dimensional (2D) elliptic problems with variable coefficients and interior/boundary layers on nonuniform grids. Firstly, the original equation is transformed from the physical domain (with a nonuniform mesh) to the computational domain (with a uniform mesh) by using a coordinate transformation. Then, a fourth order compact difference scheme is proposed to solve the transformed elliptic equation on uniform girds. After that, a multigrid method is employed to solve the linear algebraic system arising from the difference equation. At last, the numerical experiments on some elliptic problems with interior/boundary layers are conducted to show high accuracy and high efficiency of the present method.
基金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 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 National Natural Science Foundation of China(No.11601517)the Basic Research Foundation of National University of Defense Technology(No.ZDYYJ-CYJ20140101)
文摘To ensure time stability of a seventh-order dissipative compact finite difference scheme, fourth-order boundary closures are used near domain boundaries previously. However, this would reduce the global convergence rate to fifth-order only. In this paper, we elevate the boundary closures to sixth-order to achieve seventh-order global accuracy. To keep the improved scheme time stable, the simultaneous approximation terms (SATs) are used to impose boundary conditions weakly. Eigenvalue analysis shows that the improved scheme is time stable. Numerical experiments for linear advection equations and one-dimensional Euler equations are implemented to validate the new scheme.
基金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.
基金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.
文摘We have successfully attempted to solve the equations of full-MHD model within the framework of Ψ- ωformulation with an objective to evaluate the performance of a new higher order scheme to predict better values of control parameters of the flow. In particular for MHD flows, magnetic field and electrical conductivity are the control parameters. In this work, the results from our efficient high order accurate scheme are compared with the results of second order method and significant discrepancies are noted in separation length, drag coefficient and mean Nusselt number. The governing Navier-Stokes equation is fully nonlinear due to its coupling with Maxwell’s equations. The momentum equation has several highly nonlinear body-force terms due to full-MHD model in cylindrical polar system. Our high accuracy results predict that a relatively lower magnetic field is sufficient to achieve full suppression of boundary layer and this is a favorable result for practical applications. The present computational scheme predicts that a drag-coefficient minimum can be achieved when β=0.4 which is much lower when compared to the value β=1 as given by second order method. For a special value of β=0.65, it is found that the heat transfer rate is independent of electrical conductivity of the fluid. From the numerical values of physical quantities, we establish that the order of accuracy of the computed numerical results is fourth order accurate by using the method of divided differences.
基金supported by the National Natural Science Foundation of China(No.51174236)the National Basic Research Program of China(973 Program)(No.2011CB606306)the Opening Project of State Key Laboratory of Porous Metal Materials(No.PMM-SKL-4-2012)
文摘An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.
文摘A new class of finite difference schemes--the weighted compact schemes are proposed. According to the idea of the WENO schemes, the weighted compact scheme is constructed by a combination of the approximations of derivatives on candidate stencils with properly assigned weights so that the non oscillatory property is achieved when discontinuities appear. The primitive function reconstruction method of ENO schemes is applied to obtain the conservative form of the weighted compact scheme. This new scheme not only preserves the characteristic of standard compact schemes and achieves high order accuracy and high resolution using a compact stencil, but also can accurately capture shock waves and discontinuities without oscillation. Numerical examples show that the new scheme is very promising and successful.
文摘By using Richardson extrapolation and fourth-order compact finite difference scheme on different scale grids, a sixth-order solution is computed on the coarse grid. Other three techniques are applied to obtain a sixth-order solution on the fine grid, and thus give out three kinds of Richardson extrapolation-based sixth order compact computation methods. By carefully analyzing the truncation errors respectively on 2D Poisson equation, we compare the accuracy of these three sixth order methods theoretically. Numerical results for two test problems are discussed.
文摘In this paper, we introduce high-order finite volume methods for the multi-term time fractional sub-diffusion equation. The time fractional derivatives are described in Caputo’s sense. By using some operators, we obtain the compact finite volume scheme have high order accuracy. We use a compact operator to deal with spatial direction;then we can get the compact finite volume scheme. It is proved that the finite volume scheme is unconditionally stable and convergent in L<sub>∞</sub>-norm. The convergence order is O(τ<sup>2-α</sup> + h<sup>4</sup>). Finally, two numerical examples are given to confirm the theoretical results. Some tables listed also can explain the stability and convergence of the scheme.
基金This work is supported by the NNSFC （Nos. 11771213, 41504078, 11301234, 11271171）, the National Key Research and Development Project of China （No. 2016YFC0600310）, the Major Projects of Natural Sciences of University in Jiangsu Province of China （No. 15KJA110002） and the Priority Academic Program Development of Jiangsu Higher Education Institutions, the Provincial Natural Science Foundation of Jiangxi （Nos. 20161ACB20006, 20142BCB23009, 20151BAB 201012）.
文摘In this paper, a novel multisymplectic scheme is proposed for the coupled nonlinear Schrodinger-KdV （CNLS-KdV） equations. The CNLS-KdV equations are rewritten into the multisymplectic Hamiltonian form by introducing some canonical momenta. To simulate the problem efficiently, the CNLS-KdV equations are approximated by a high order compact method in space which preserves N semi-discrete multisymplectic conservation laws. We then discretize the semi-discrete system by using a symplectic midpoint scheme in time. Thus, a full-discrete multisymplectic scheme is obtained for the CNLS-KdV equations. The conservation laws of the full-discrete scheme are analyzed. Some numerical experiments are presented to further verify the convergence and conservation laws of the new scheme.
文摘A high-order finite difference Pade scheme also called compact scheme for solving Korteweg-de Vries （KdV） equations, which preserve energy and mass conservations, was developed in this paper. This structure-preserving algorithm has been widely applied in these years for its advantage of maintaining the inherited properties. For spatial discretization, the authors obtained an implicit compact scheme by which spatial derivative terms may be approximated through combining a few knots. By some numerical examples including propagation of single soliton and interaction of two solitons, the scheme is proved to be effective.
