A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （200...A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （2004） 056702]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative （NND） scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model is validated by well-known benchmarks, （i） Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; （ii） reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.展开更多
The weighted residuals method was used for obtaining the boundary integral representation of the velocity of the three-dimensional inviscid irrotational flow. It is shown that velocity in an arbitrary point of domain ...The weighted residuals method was used for obtaining the boundary integral representation of the velocity of the three-dimensional inviscid irrotational flow. It is shown that velocity in an arbitrary point of domain can be expressed through its values on the boundary. Boundary integral equations of the second kind for solving boundary-valued problems of the first and second kinds are developed. The result has been also generalised to the case of solenoidal vector fields with potential vorticity. It is shown that the resulting integral equations are Fredholm integral equations of the second kind and allow effective numerical solving of corresponding boundary-valued problems. Examples of numerical solutions for a sphere and an ellipsoid are given for demonstration of efficiency of the offered method.展开更多
Due to the complexity of compressible flows,nonlinear hydrodynamic stability theories in supersonic boundary layers are not sufficient.In order to reveal the nonlinear interaction mechanisms of the rapidly amplified 3...Due to the complexity of compressible flows,nonlinear hydrodynamic stability theories in supersonic boundary layers are not sufficient.In order to reveal the nonlinear interaction mechanisms of the rapidly amplified 3-D disturbances in supersonic boundary layers at high Mach numbers,the nonlinear evolutions of different disturbances in flat-plate boundary layers at Mach number 4.5,6 and 8 are analyzed by numerical simulations.It can be concluded that the 3-D disturbances are amplified rapidly when the amplitude of the 2-D disturbance reaches a certain level.The most rapidly amplified 3-D disturbances are Klebanoff type(K-type)disturbances which have the same frequency as the 2-D disturbance.Among these K-type 3-D disturbances,the disturbances located at the junction of upper branch and lower branch of the neutral curve are amplified higher.Through analyzing the relationship between the amplification rate and the spanwise wavenumber of the 3-D disturbances at different evolution stages,the mechanism of the spanwise wavenumber selectivity of K-type 3-D disturbances in the presence of a finite amplitude 2-D disturbance is explained.展开更多
Based on the finite element method(FEM), some iterative methods related to different Reynolds numbers are designed and analyzed for solving the 2D/3D stationary incompressible magnetohydrodynamics(MHD) numerically. Tw...Based on the finite element method(FEM), some iterative methods related to different Reynolds numbers are designed and analyzed for solving the 2D/3D stationary incompressible magnetohydrodynamics(MHD) numerically. Two-level finite element iterative methods, consisting of the classical m-iteration methods on a coarse grid and corrections on a fine grid, are designed to solve the system at low Reynolds numbers under the strong uniqueness condition. One-level Oseen-type iterative method is investigated on a fine mesh at high Reynolds numbers under the weak uniqueness condition. Furthermore, the uniform stability and convergence of these methods with respect to equation parameters R_e, R_m, S_c, mesh sizes h, H and iterative step m are provided. Finally, the efficiency of the proposed methods is confirmed by numerical investigations.展开更多
基金Supported by the Science Foundations of Laboratory of Computational PhysicalScience Foundation of China Academy of Engineering Physics under Grant Nos. 2009A0102005, 2009B0101012National Natural Science Foundation under Grant Nos. 10775018, 11074300, and 1107521 of China
文摘A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （2004） 056702]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative （NND） scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model is validated by well-known benchmarks, （i） Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; （ii） reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.
文摘The weighted residuals method was used for obtaining the boundary integral representation of the velocity of the three-dimensional inviscid irrotational flow. It is shown that velocity in an arbitrary point of domain can be expressed through its values on the boundary. Boundary integral equations of the second kind for solving boundary-valued problems of the first and second kinds are developed. The result has been also generalised to the case of solenoidal vector fields with potential vorticity. It is shown that the resulting integral equations are Fredholm integral equations of the second kind and allow effective numerical solving of corresponding boundary-valued problems. Examples of numerical solutions for a sphere and an ellipsoid are given for demonstration of efficiency of the offered method.
基金supported by the State Key Program of National Natural Science Foundation of China(Grant No.11332007)
文摘Due to the complexity of compressible flows,nonlinear hydrodynamic stability theories in supersonic boundary layers are not sufficient.In order to reveal the nonlinear interaction mechanisms of the rapidly amplified 3-D disturbances in supersonic boundary layers at high Mach numbers,the nonlinear evolutions of different disturbances in flat-plate boundary layers at Mach number 4.5,6 and 8 are analyzed by numerical simulations.It can be concluded that the 3-D disturbances are amplified rapidly when the amplitude of the 2-D disturbance reaches a certain level.The most rapidly amplified 3-D disturbances are Klebanoff type(K-type)disturbances which have the same frequency as the 2-D disturbance.Among these K-type 3-D disturbances,the disturbances located at the junction of upper branch and lower branch of the neutral curve are amplified higher.Through analyzing the relationship between the amplification rate and the spanwise wavenumber of the 3-D disturbances at different evolution stages,the mechanism of the spanwise wavenumber selectivity of K-type 3-D disturbances in the presence of a finite amplitude 2-D disturbance is explained.
基金National Natural Science Foundation of China (Grant Nos. 11271298 and 11362021)
文摘Based on the finite element method(FEM), some iterative methods related to different Reynolds numbers are designed and analyzed for solving the 2D/3D stationary incompressible magnetohydrodynamics(MHD) numerically. Two-level finite element iterative methods, consisting of the classical m-iteration methods on a coarse grid and corrections on a fine grid, are designed to solve the system at low Reynolds numbers under the strong uniqueness condition. One-level Oseen-type iterative method is investigated on a fine mesh at high Reynolds numbers under the weak uniqueness condition. Furthermore, the uniform stability and convergence of these methods with respect to equation parameters R_e, R_m, S_c, mesh sizes h, H and iterative step m are provided. Finally, the efficiency of the proposed methods is confirmed by numerical investigations.
