Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin（DG） discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta（IRK） scheme...Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin（DG） discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta（IRK） scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.展开更多
In this paper, based on the idea of El-Mistikawy and Werle[1] we construct a difference scheme for a singularly perturbed self-adjoint ordinary differential equation in conservation form. We prove that it is a uniform...In this paper, based on the idea of El-Mistikawy and Werle[1] we construct a difference scheme for a singularly perturbed self-adjoint ordinary differential equation in conservation form. We prove that it is a uniformly convergent second order scheme.展开更多
In this paper a singularly perturbed linear second order hyperbolic problem with zeroth order reduced equation is discussed. Firstly, an energy inequality of the solution and an estimate of the remainder term of the a...In this paper a singularly perturbed linear second order hyperbolic problem with zeroth order reduced equation is discussed. Firstly, an energy inequality of the solution and an estimate of the remainder term of the asymptotic solution are given. Then an exponentially fitted difference scheme is developed in an equidistant mesh. Finally, uniform convergence in small parameter is proved in the sense of discrete energy norm.展开更多
Computational aeroacoustics （CAA） is an interdiscipline of aeroacoustics and computational fluid dynamics （CFD） for the investigation of sound generation and propagation from various aeroacoustics problems. In thi...Computational aeroacoustics （CAA） is an interdiscipline of aeroacoustics and computational fluid dynamics （CFD） for the investigation of sound generation and propagation from various aeroacoustics problems. In this review, the foundation and research scope of CAA are introduced firstly. A review of the early advances and applications of CAA is then briefly surveyed, focusing on two key issues, namely, high order finite difference scheme and non-reflecting boundary condition. Furthermore, the advances of CAA during the past five years are highlighted. Finally, the future prospective of CAA is briefly discussed.展开更多
This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspecti...This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternative time-clocking. With the time-clocking perspective, unsteady flow solutions from different passages of the same blade row are mapped to flow solutions of the same passage at different time instants or phase angles. Accordingly, the time derivative of the unsteady flow equation is discretized in time directly, which is more natural than transforming the time derivative to a spatial one as with the original STG method. To improve the solution accuracy, a ninth order difference scheme has been investigated for discretizing the time derivative. To achieve a stable solution for the high order scheme, the implicit solution method of Lower-Upper Symmetric GaussSeidel/Gauss-Seidel(LU-SGS/GS) has been employed. The NASA Stage 35 and its blade-countreduced variant are used to demonstrate the validity of the time-clocking based passage reordering and the advantages of the high order difference scheme for the STG method. Results from an existing harmonic balance flow solver are also provided to contrast the two methods in terms of solution stability and computational cost.展开更多
A new shock-capturing method is proposed which is based on upwind schemes and flux-vector splittings. Firstly, original upwind schemes are projected along characteristic directions. Secondly, the amplitudes of the cha...A new shock-capturing method is proposed which is based on upwind schemes and flux-vector splittings. Firstly, original upwind schemes are projected along characteristic directions. Secondly, the amplitudes of the characteristic decompositions are carefully controlled by limiters to prevent non-physical oscillations. Lastly, the schemes are converted into conservative forms, and the oscillation-free shock-capturing schemes are acquired. Two explicit upwind schemes （2nd-order and 3rd-order） and three compact upwind schemes （3rd-order, 5th-order and 7th-order） are modified by the method for hyperbolic systems and the modified schemes are checked on several one-dimensional and two-dimensional test cases. Some numerical solutions of the schemes are compared with those of a WENO scheme and a MP scheme as well as a compact-WENO scheme. The results show that the method with high order accuracy and high resolutions can capture shock waves smoothly.展开更多
The flow field with a high order scheme is usually calculated so as to solve complex flow problems and describe the flow structure accurately. However, there are two problems, i.e., the reduced-order boundary is inevi...The flow field with a high order scheme is usually calculated so as to solve complex flow problems and describe the flow structure accurately. However, there are two problems, i.e., the reduced-order boundary is inevitable and the order of the scheme at the discontinuous shock wave contained in the flow field as the supersonic flow field is low. It is questionable whether the reduced-order boundary and the low-order scheme at the shock wave have an effect on the numerical solution and accuracy of the flow field inside. In this paper, according to the actual situation of the direct numerical simulation of the flow field, two model equations with the exact solutions are solved, which are steady and unsteady, respectively, to study the question with a high order scheme at the interior of the domain and the reduced-order method at the boundary and center of the domain. Comparing with the exact solutions, it is found that the effect of reduced-order exists and cannot be ignored. In addition, the other two model equations with the exact solutions, which are often used in fluid mechanics, are also studied with the same process for the reduced-order problem.展开更多
Two types of implicit algorithms have been im- proved for high order discontinuous Galerkin （DG） method to solve compressible Navier-Stokes （NS） equations on tri- angular grids. A block lower-upper symmetric Gauss...Two types of implicit algorithms have been im- proved for high order discontinuous Galerkin （DG） method to solve compressible Navier-Stokes （NS） equations on tri- angular grids. A block lower-upper symmetric Gauss- Seidel （BLU-SGS） approach is implemented as a nonlin- ear iterative scheme. And a modified LU-SGS （LLU-SGS） approach is suggested to reduce the memory requirements while retain the good convergence performance of the origi- nal LU-SGS approach. Both implicit schemes have the sig- nificant advantage that only the diagonal block matrix is stored. The resulting implicit high-order DG methods are applied, in combination with Hermite weighted essentially non-oscillatory （HWENO） limiters, to solve viscous flow problems. Numerical results demonstrate that the present implicit methods are able to achieve significant efficiency improvements over explicit counterparts and for viscous flows with shocks, and the HWENO limiters can be used to achieve the desired essentially non-oscillatory shock tran- sition and the designed high-order accuracy simultaneously.展开更多
Numerical techniques play an important role in CFD. Some of them are reviewed in this paper. The necessity of using high order difference scheme is demonstrated for the study of high Reynolds number viscous how. Physi...Numerical techniques play an important role in CFD. Some of them are reviewed in this paper. The necessity of using high order difference scheme is demonstrated for the study of high Reynolds number viscous how. Physical guide lines are provided for the construction of these high order schemes. To avoid unduly ad hoc treatment in the boundary region the use of compact scheme is recommended because it has a small stencil size compared with the traditional finite difference scheme. Besides preliminary Fourier analysis shows the compact scheme can also yield better space resolution which makes it more suitable to study flow with multiscales e.g. turbulence. Other approaches such as perturbation method and finite spectral method are also emphasized. Typical numerical simulations were carried out. The first deals with Euler equations to show its capabilities to capture flow discontinuity. The second deals with Navier-Stokes Equations studying the evolution of a mixing layer, the pertinent structures at different times are shown. Asymmetric break down occurs and also the appearance of small vortices.展开更多
The convergence results of block iterative schemes from the EG （Explicit Group） family have been shown to be one of efficient iterative methods in solving any linear systems generated from approximation equations. A...The convergence results of block iterative schemes from the EG （Explicit Group） family have been shown to be one of efficient iterative methods in solving any linear systems generated from approximation equations. Apart from block iterative methods, the formulation of the MSOR （Modified Successive Over-Relaxation） method known as SOR method with red-black ordering strategy by using two accelerated parameters, ω and ω′, has also improved the convergence rate of the standard SOR method. Due to the effectiveness of these iterative methods, the primary goal of this paper is to examine the performance of the EG family without or with accelerated parameters in solving second order two-point nonlinear boundary value problems. In this work, the second order two-point nonlinear boundary value problems need to be discretized by using the second order central difference scheme in constructing a nonlinear finite difference approximation equation. Then this approximation equation leads to a nonlinear system. As well known that to linearize nonlinear systems, the Newton method has been proposed to transform the original system into the form of linear system. In addition to that, the basic formulation and implementation of 2 and 4-point EG iterative methods based on GS （Gauss-Seidel）, SOR and MSOR approaches, namely EGGS, EGSOR and EGMSOR respectively are also presented. Then, combinations between the EG family and Newton scheme are indicated as EGGS-Newton, EGSOR-Newton and EGMSOR-Newton methods respectively. For comparison purpose, several numerical experiments of three problems are conducted in examining the effectiveness of tested methods. Finally, it can be concluded that the 4-point EGMSOR-Newton method is more superior in accelerating the convergence rate compared with the tested methods.展开更多
In this paper, a two-dimensional nonlinear coupled Gray Scott system is simulated with a finite difference scheme and a finite volume technique. Pre and post-processing lead to a new solution called GSmFoam by underst...In this paper, a two-dimensional nonlinear coupled Gray Scott system is simulated with a finite difference scheme and a finite volume technique. Pre and post-processing lead to a new solution called GSmFoam by understandin<span>g geometry settings and mesh information. The concentration profile chan</span>ges over time, as does the intensity of the contour patterns. The OpenFoam solver gives you the confidence to compare the pattern result with efficient numerical algorithms on the Gray Scott model.展开更多
In this paper,we consider numerical solutions of the fractional diffusion equation with theαorder time fractional derivative defined in the Caputo-Hadamard sense.A high order time-stepping scheme is constructed,analy...In this paper,we consider numerical solutions of the fractional diffusion equation with theαorder time fractional derivative defined in the Caputo-Hadamard sense.A high order time-stepping scheme is constructed,analyzed,and numerically validated.The contribution of the paper is twofold:1)regularity of the solution to the underlying equation is investigated,2)a rigorous stability and convergence analysis for the proposed scheme is performed,which shows that the proposed scheme is 3+αorder accurate.Several numerical examples are provided to verify the theoretical statement.展开更多
With the undetermined coefficient method, a new high order scheme for convection equation named as HAUC was worke d out. The accuracy of the new scheme was analyzed by comparing the computed res ults with the exact s...With the undetermined coefficient method, a new high order scheme for convection equation named as HAUC was worke d out. The accuracy of the new scheme was analyzed by comparing the computed res ults with the exact solutions of 1-D pure convection equation and nonlinear con vection equation. The effectiveness of the HAUC was also examined by comparing w ith the results obtained by other schemes. Finally, it was applied to simulate 1 -D dam break flow with differenct ratios of initial upstream water depth to do wnstream one. The results show that the scheme has the ability to simulate both undular bores and moving hydraulic jumps.展开更多
The computation of compressible flows at all Mach numbers is a very challenging problem.An efficient numerical method for solving this problem needs to have shock-capturing capability in the high Mach number regime,wh...The computation of compressible flows at all Mach numbers is a very challenging problem.An efficient numerical method for solving this problem needs to have shock-capturing capability in the high Mach number regime,while it can deal with stiffness and accuracy in the low Mach number regime.This paper designs a high order semi-implicit weighted compact nonlinear scheme(WCNS)for the all-Mach isentropic Euler system of compressible gas dynamics.To avoid severe Courant-Friedrichs-Levy(CFL)restrictions for low Mach flows,the nonlinear fluxes in the Euler equations are split into stiff and non-stiff components.A third-order implicit-explicit(IMEX)method is used for the time discretization of the split components and a fifth-order WCNS is used for the spatial discretization of flux derivatives.The high order IMEX method is asymptotic preserving and asymptotically accurate in the zero Mach number limit.One-and two-dimensional numerical examples in both compressible and incompressible regimes are given to demonstrate the advantages of the designed IMEX WCNS.展开更多
In this paper,we introduce the discrete Maxwellian equilibrium distribution function for incompressible flow and force term into the two-stage third-order Discrete Unified Gas-Kinetic Scheme(DUGKS)for simulating low-s...In this paper,we introduce the discrete Maxwellian equilibrium distribution function for incompressible flow and force term into the two-stage third-order Discrete Unified Gas-Kinetic Scheme(DUGKS)for simulating low-speed turbulent flows.The Wall-Adapting Local Eddy-viscosity(WALE)and Vreman sub-grid models for Large-Eddy Simulations(LES)of turbulent flows are coupled within the present framework.Meanwhile,the implicit LES are also presented to verify the effect of LES models.A parallel implementation strategy for the present framework is developed,and three canonical wall-bounded turbulent flow cases are investigated,including the fully developed turbulent channel flow at a friction Reynolds number(Re)about 180,the turbulent plane Couette flow at a friction Re number about 93 and lid-driven cubical cavity flow at a Re number of 12000.The turbulence statistics,including mean velocity,the r.m.s.fluctuations velocity,Reynolds stress,etc.are computed by the present approach.Their predictions match precisely with each other,and they are both in reasonable agreement with the benchmark data of DNS.Especially,the predicted flow physics of three-dimensional lid-driven cavity flow are consistent with the description from abundant literature.The present numerical results verify that the present two-stage third-order DUGKS-based LES method is capable for simulating inhomogeneous wall-bounded turbulent flows and getting reliable results with relatively coarse grids.展开更多
A higher-order compact scheme on the nine point 2-D stencil is developed for the steady stream-function vorticity form of the incompressible Navier-Stokes(NS)equations in spherical polar coordinates,which was used ear...A higher-order compact scheme on the nine point 2-D stencil is developed for the steady stream-function vorticity form of the incompressible Navier-Stokes(NS)equations in spherical polar coordinates,which was used earlier only for the cartesian and cylindrical geometries.The steady,incompressible,viscous and axially symmetric flow past a sphere is used as a model problem.The non-linearity in the N-S equations is handled in a comprehensive manner avoiding complications in calculations.The scheme is combined with the multigrid method to enhance the convergence rate.The solutions are obtained over a non-uniform grid generated using the transformation r=ex while maintaining a uniform grid in the computational plane.The superiority of the higher order compact scheme is clearly illustrated in comparison with upwind scheme and defect correction technique at high Reynolds numbers by taking a large domain.This is a pioneering effort,because for the first time,the fourth order accurate solutions for the problem of viscous flow past a sphere are presented here.The drag coefficient and surface pressures are calculated and compared with available experimental and theoretical results.It is observed that these values simulated over coarser grids using the present scheme aremore accuratewhen compared to other conventional schemes.It has also been observed that the flow separation initially occurred at Re=21.展开更多
In this paper, a difference scheme with nonuniform meshes is proposed for the initial-boundary problem of the nonlinear parabolic system. It is proved that the difference scheme is second order convergent in both spac...In this paper, a difference scheme with nonuniform meshes is proposed for the initial-boundary problem of the nonlinear parabolic system. It is proved that the difference scheme is second order convergent in both space and time.展开更多
We present a new high ordermethod in space and time for solving the wave equation,based on a newinterpretation of the“Modified Equation”technique.Indeed,contrary to most of the works,we consider the time discretizat...We present a new high ordermethod in space and time for solving the wave equation,based on a newinterpretation of the“Modified Equation”technique.Indeed,contrary to most of the works,we consider the time discretization before the space discretization.After the time discretization,an additional biharmonic operator appears,which can not be discretized by classical finite elements.We propose a new Discontinuous Galerkinmethod for the discretization of this operator,andwe provide numerical experiments proving that the new method is more accurate than the classicalModified Equation technique with a lower computational burden.展开更多
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.展开更多
Various compact difference schemes （both old and new, explicit and implicit, one-level and two-level）, which approximate the diffusion equation and SchrSdinger equation with periodical boundary conditions are constr...Various compact difference schemes （both old and new, explicit and implicit, one-level and two-level）, which approximate the diffusion equation and SchrSdinger equation with periodical boundary conditions are constructed by means of the general approach. The results of numerical experiments for various initial data and right hand side are presented. We evaluate the real order of their convergence, as well as their stability, effectiveness, and various kinds of monotony. The optimal Courant number depends on the number of grid knots and on the smoothness of solutions. The competition of various schemes should be organized for the fixed number of arithmetic operations, which are necessary for numerical integration of a given Cauchy problem. This approach to the construction of compact schemes can be developed for numerical solution of various problems of mathematical physics.展开更多
文摘Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin（DG） discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta（IRK） scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.
文摘In this paper, based on the idea of El-Mistikawy and Werle[1] we construct a difference scheme for a singularly perturbed self-adjoint ordinary differential equation in conservation form. We prove that it is a uniformly convergent second order scheme.
文摘In this paper a singularly perturbed linear second order hyperbolic problem with zeroth order reduced equation is discussed. Firstly, an energy inequality of the solution and an estimate of the remainder term of the asymptotic solution are given. Then an exponentially fitted difference scheme is developed in an equidistant mesh. Finally, uniform convergence in small parameter is proved in the sense of discrete energy norm.
基金Project supported by the National Basic Research Program of China(No.2012CB720202)the National Natural Science Foundation of China(No.51476005)the 111 Project of China(No.B07009)
文摘Computational aeroacoustics （CAA） is an interdiscipline of aeroacoustics and computational fluid dynamics （CFD） for the investigation of sound generation and propagation from various aeroacoustics problems. In this review, the foundation and research scope of CAA are introduced firstly. A review of the early advances and applications of CAA is then briefly surveyed, focusing on two key issues, namely, high order finite difference scheme and non-reflecting boundary condition. Furthermore, the advances of CAA during the past five years are highlighted. Finally, the future prospective of CAA is briefly discussed.
基金co-supported by the National Natural Science Foundation of China(No.51976172)the National Science and Technology Major Project of China(No.2017-Ⅱ-0009-0023)。
文摘This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternative time-clocking. With the time-clocking perspective, unsteady flow solutions from different passages of the same blade row are mapped to flow solutions of the same passage at different time instants or phase angles. Accordingly, the time derivative of the unsteady flow equation is discretized in time directly, which is more natural than transforming the time derivative to a spatial one as with the original STG method. To improve the solution accuracy, a ninth order difference scheme has been investigated for discretizing the time derivative. To achieve a stable solution for the high order scheme, the implicit solution method of Lower-Upper Symmetric GaussSeidel/Gauss-Seidel(LU-SGS/GS) has been employed. The NASA Stage 35 and its blade-countreduced variant are used to demonstrate the validity of the time-clocking based passage reordering and the advantages of the high order difference scheme for the STG method. Results from an existing harmonic balance flow solver are also provided to contrast the two methods in terms of solution stability and computational cost.
基金Project supported by the National Natural Science Foundation of China (Nos.10321002 and 10672012)
文摘A new shock-capturing method is proposed which is based on upwind schemes and flux-vector splittings. Firstly, original upwind schemes are projected along characteristic directions. Secondly, the amplitudes of the characteristic decompositions are carefully controlled by limiters to prevent non-physical oscillations. Lastly, the schemes are converted into conservative forms, and the oscillation-free shock-capturing schemes are acquired. Two explicit upwind schemes （2nd-order and 3rd-order） and three compact upwind schemes （3rd-order, 5th-order and 7th-order） are modified by the method for hyperbolic systems and the modified schemes are checked on several one-dimensional and two-dimensional test cases. Some numerical solutions of the schemes are compared with those of a WENO scheme and a MP scheme as well as a compact-WENO scheme. The results show that the method with high order accuracy and high resolutions can capture shock waves smoothly.
基金Project supported by the National Key Research and Development Project of China(No.2016YFA0401200)the National Natural Science Foundation of China(Nos.11672205 and11332007)
文摘The flow field with a high order scheme is usually calculated so as to solve complex flow problems and describe the flow structure accurately. However, there are two problems, i.e., the reduced-order boundary is inevitable and the order of the scheme at the discontinuous shock wave contained in the flow field as the supersonic flow field is low. It is questionable whether the reduced-order boundary and the low-order scheme at the shock wave have an effect on the numerical solution and accuracy of the flow field inside. In this paper, according to the actual situation of the direct numerical simulation of the flow field, two model equations with the exact solutions are solved, which are steady and unsteady, respectively, to study the question with a high order scheme at the interior of the domain and the reduced-order method at the boundary and center of the domain. Comparing with the exact solutions, it is found that the effect of reduced-order exists and cannot be ignored. In addition, the other two model equations with the exact solutions, which are often used in fluid mechanics, are also studied with the same process for the reduced-order problem.
基金supported by the National Basic Research Program of China(2009CB724104)
文摘Two types of implicit algorithms have been im- proved for high order discontinuous Galerkin （DG） method to solve compressible Navier-Stokes （NS） equations on tri- angular grids. A block lower-upper symmetric Gauss- Seidel （BLU-SGS） approach is implemented as a nonlin- ear iterative scheme. And a modified LU-SGS （LLU-SGS） approach is suggested to reduce the memory requirements while retain the good convergence performance of the origi- nal LU-SGS approach. Both implicit schemes have the sig- nificant advantage that only the diagonal block matrix is stored. The resulting implicit high-order DG methods are applied, in combination with Hermite weighted essentially non-oscillatory （HWENO） limiters, to solve viscous flow problems. Numerical results demonstrate that the present implicit methods are able to achieve significant efficiency improvements over explicit counterparts and for viscous flows with shocks, and the HWENO limiters can be used to achieve the desired essentially non-oscillatory shock tran- sition and the designed high-order accuracy simultaneously.
基金The project supported by the National Natural Science Foundation of China(No.19393100)
文摘Numerical techniques play an important role in CFD. Some of them are reviewed in this paper. The necessity of using high order difference scheme is demonstrated for the study of high Reynolds number viscous how. Physical guide lines are provided for the construction of these high order schemes. To avoid unduly ad hoc treatment in the boundary region the use of compact scheme is recommended because it has a small stencil size compared with the traditional finite difference scheme. Besides preliminary Fourier analysis shows the compact scheme can also yield better space resolution which makes it more suitable to study flow with multiscales e.g. turbulence. Other approaches such as perturbation method and finite spectral method are also emphasized. Typical numerical simulations were carried out. The first deals with Euler equations to show its capabilities to capture flow discontinuity. The second deals with Navier-Stokes Equations studying the evolution of a mixing layer, the pertinent structures at different times are shown. Asymmetric break down occurs and also the appearance of small vortices.
文摘The convergence results of block iterative schemes from the EG （Explicit Group） family have been shown to be one of efficient iterative methods in solving any linear systems generated from approximation equations. Apart from block iterative methods, the formulation of the MSOR （Modified Successive Over-Relaxation） method known as SOR method with red-black ordering strategy by using two accelerated parameters, ω and ω′, has also improved the convergence rate of the standard SOR method. Due to the effectiveness of these iterative methods, the primary goal of this paper is to examine the performance of the EG family without or with accelerated parameters in solving second order two-point nonlinear boundary value problems. In this work, the second order two-point nonlinear boundary value problems need to be discretized by using the second order central difference scheme in constructing a nonlinear finite difference approximation equation. Then this approximation equation leads to a nonlinear system. As well known that to linearize nonlinear systems, the Newton method has been proposed to transform the original system into the form of linear system. In addition to that, the basic formulation and implementation of 2 and 4-point EG iterative methods based on GS （Gauss-Seidel）, SOR and MSOR approaches, namely EGGS, EGSOR and EGMSOR respectively are also presented. Then, combinations between the EG family and Newton scheme are indicated as EGGS-Newton, EGSOR-Newton and EGMSOR-Newton methods respectively. For comparison purpose, several numerical experiments of three problems are conducted in examining the effectiveness of tested methods. Finally, it can be concluded that the 4-point EGMSOR-Newton method is more superior in accelerating the convergence rate compared with the tested methods.
文摘In this paper, a two-dimensional nonlinear coupled Gray Scott system is simulated with a finite difference scheme and a finite volume technique. Pre and post-processing lead to a new solution called GSmFoam by understandin<span>g geometry settings and mesh information. The concentration profile chan</span>ges over time, as does the intensity of the contour patterns. The OpenFoam solver gives you the confidence to compare the pattern result with efficient numerical algorithms on the Gray Scott model.
基金supported by the Fujian Alliance of Mathematics(Grant No.2024SXLMMS03)by the Natural Science Foundation of Fujian Province of China(Grant No.2022J01338)+4 种基金supported by the NSFC(Grant Nos.12361083,62341115)by the Foundation of Guizhou Science and Technology Department(Grant No.QHKJC-ZK[2024]YB497)by the Natural Science Research Project of Department of Education of Guizhou Province(Grant No.QJJ2023012)by the Science Research Fund Support Project of the Guizhou Minzu University(Grant No.GZMUZK[2023]CXTD05)supported by the NSFC(Grant No.12371408).
文摘In this paper,we consider numerical solutions of the fractional diffusion equation with theαorder time fractional derivative defined in the Caputo-Hadamard sense.A high order time-stepping scheme is constructed,analyzed,and numerically validated.The contribution of the paper is twofold:1)regularity of the solution to the underlying equation is investigated,2)a rigorous stability and convergence analysis for the proposed scheme is performed,which shows that the proposed scheme is 3+αorder accurate.Several numerical examples are provided to verify the theoretical statement.
文摘With the undetermined coefficient method, a new high order scheme for convection equation named as HAUC was worke d out. The accuracy of the new scheme was analyzed by comparing the computed res ults with the exact solutions of 1-D pure convection equation and nonlinear con vection equation. The effectiveness of the HAUC was also examined by comparing w ith the results obtained by other schemes. Finally, it was applied to simulate 1 -D dam break flow with differenct ratios of initial upstream water depth to do wnstream one. The results show that the scheme has the ability to simulate both undular bores and moving hydraulic jumps.
基金the National Numerical Wind Tunnel Project(No.NNW2018-ZT4A08)the National Natural Science Foundation of China(Nos.11872323 and 11971025)the Natural Science Foundation of Fujian Province(No.2019J06002)。
文摘The computation of compressible flows at all Mach numbers is a very challenging problem.An efficient numerical method for solving this problem needs to have shock-capturing capability in the high Mach number regime,while it can deal with stiffness and accuracy in the low Mach number regime.This paper designs a high order semi-implicit weighted compact nonlinear scheme(WCNS)for the all-Mach isentropic Euler system of compressible gas dynamics.To avoid severe Courant-Friedrichs-Levy(CFL)restrictions for low Mach flows,the nonlinear fluxes in the Euler equations are split into stiff and non-stiff components.A third-order implicit-explicit(IMEX)method is used for the time discretization of the split components and a fifth-order WCNS is used for the spatial discretization of flux derivatives.The high order IMEX method is asymptotic preserving and asymptotically accurate in the zero Mach number limit.One-and two-dimensional numerical examples in both compressible and incompressible regimes are given to demonstrate the advantages of the designed IMEX WCNS.
基金The National Numerical Wind Tunnel Project,the National Natural Science Foundation of China(Grant No.11902264,11902266,12072283)the 111 Project of China(B17037)as well as the ATCFD Project(2015-F-016).
文摘In this paper,we introduce the discrete Maxwellian equilibrium distribution function for incompressible flow and force term into the two-stage third-order Discrete Unified Gas-Kinetic Scheme(DUGKS)for simulating low-speed turbulent flows.The Wall-Adapting Local Eddy-viscosity(WALE)and Vreman sub-grid models for Large-Eddy Simulations(LES)of turbulent flows are coupled within the present framework.Meanwhile,the implicit LES are also presented to verify the effect of LES models.A parallel implementation strategy for the present framework is developed,and three canonical wall-bounded turbulent flow cases are investigated,including the fully developed turbulent channel flow at a friction Reynolds number(Re)about 180,the turbulent plane Couette flow at a friction Re number about 93 and lid-driven cubical cavity flow at a Re number of 12000.The turbulence statistics,including mean velocity,the r.m.s.fluctuations velocity,Reynolds stress,etc.are computed by the present approach.Their predictions match precisely with each other,and they are both in reasonable agreement with the benchmark data of DNS.Especially,the predicted flow physics of three-dimensional lid-driven cavity flow are consistent with the description from abundant literature.The present numerical results verify that the present two-stage third-order DUGKS-based LES method is capable for simulating inhomogeneous wall-bounded turbulent flows and getting reliable results with relatively coarse grids.
文摘A higher-order compact scheme on the nine point 2-D stencil is developed for the steady stream-function vorticity form of the incompressible Navier-Stokes(NS)equations in spherical polar coordinates,which was used earlier only for the cartesian and cylindrical geometries.The steady,incompressible,viscous and axially symmetric flow past a sphere is used as a model problem.The non-linearity in the N-S equations is handled in a comprehensive manner avoiding complications in calculations.The scheme is combined with the multigrid method to enhance the convergence rate.The solutions are obtained over a non-uniform grid generated using the transformation r=ex while maintaining a uniform grid in the computational plane.The superiority of the higher order compact scheme is clearly illustrated in comparison with upwind scheme and defect correction technique at high Reynolds numbers by taking a large domain.This is a pioneering effort,because for the first time,the fourth order accurate solutions for the problem of viscous flow past a sphere are presented here.The drag coefficient and surface pressures are calculated and compared with available experimental and theoretical results.It is observed that these values simulated over coarser grids using the present scheme aremore accuratewhen compared to other conventional schemes.It has also been observed that the flow separation initially occurred at Re=21.
基金Supported by the National Natural Science Foundation of China(No.10671060,No.10871061)the Youth Foundation of Hunan Education Bureau(No.06B037)the Construct Program of the Key Discipline in Hunan Province
文摘In this paper, a difference scheme with nonuniform meshes is proposed for the initial-boundary problem of the nonlinear parabolic system. It is proved that the difference scheme is second order convergent in both space and time.
基金supported by the Conseil General des Pyrenees Atlantiques.
文摘We present a new high ordermethod in space and time for solving the wave equation,based on a newinterpretation of the“Modified Equation”technique.Indeed,contrary to most of the works,we consider the time discretization before the space discretization.After the time discretization,an additional biharmonic operator appears,which can not be discretized by classical finite elements.We propose a new Discontinuous Galerkinmethod for the discretization of this operator,andwe provide numerical experiments proving that the new method is more accurate than the classicalModified Equation technique with a lower computational burden.
文摘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.
文摘Various compact difference schemes （both old and new, explicit and implicit, one-level and two-level）, which approximate the diffusion equation and SchrSdinger equation with periodical boundary conditions are constructed by means of the general approach. The results of numerical experiments for various initial data and right hand side are presented. We evaluate the real order of their convergence, as well as their stability, effectiveness, and various kinds of monotony. The optimal Courant number depends on the number of grid knots and on the smoothness of solutions. The competition of various schemes should be organized for the fixed number of arithmetic operations, which are necessary for numerical integration of a given Cauchy problem. This approach to the construction of compact schemes can be developed for numerical solution of various problems of mathematical physics.
