We present an improved lattice Boltzmann(LB) model for thermal liquid-vapor system.In the new model,the Windowed Fast Fourier Transform(WFFT) and its inverse are used to calculate both the convection term and the exte...We present an improved lattice Boltzmann(LB) model for thermal liquid-vapor system.In the new model,the Windowed Fast Fourier Transform(WFFT) and its inverse are used to calculate both the convection term and the external force term of the LB equation.By adopting the WFFT scheme,Gibbs oscillations can be damped effectively in unsmooth regions while high resolution feature of the spectral method can be retained in smooth regions.As a result,spatial discretization errors are dramatically decreased,conservation of the total energy is much better preserved,and the spurious velocities near the liquid-vapor interface are significantly reduced.The high resolution,together with the low complexity of the WFFT approach,endows the proposed method with considerable potential for studying a wide class of problems in the field of multiphase flows.展开更多
The aims of the present paper are twofold. At first, we further study the Multiple-Relaxation-Time （MRT） Lattice Boltzmann （LB） model proposed in [Europhys. Lett. 90 （2010） 54003]. We discuss the reason why the ...The aims of the present paper are twofold. At first, we further study the Multiple-Relaxation-Time （MRT） Lattice Boltzmann （LB） model proposed in [Europhys. Lett. 90 （2010） 54003]. We discuss the reason why the Gram Schmidt orthogonalization procedure is not needed in the construction of transformation matrix M; point out a reason why the Kataoka-Tsutahara model [Phys. Rev. E 69 （2004） 035701（R）] is only valid in subsonic flows. The yon Neumann stability analysis is performed. Secondly, we carry out a preliminary quantitative study on the Richtmyer- Meshkov instability using the proposed MRT LB model. When a shock wave travels from a light medium to a heavy one, the simulated growth rate is in qualitative agreement with the perturbation model by Zhang-Sohn. It is about half of the predicted value by the impulsive model and is closer to the experimental result. When the shock wave travels from a heavy medium to a light one, our simulation results are also consistent with physical analysis.展开更多
We report the formation of jet-like long spike in the nonlinear evolution of the ablative Rayleigh-Taylor instability （ARTI） experiments by numerical simulations. A preheating model k（T） = KSH[1＋f（T）], where KS...We report the formation of jet-like long spike in the nonlinear evolution of the ablative Rayleigh-Taylor instability （ARTI） experiments by numerical simulations. A preheating model k（T） = KSH[1＋f（T）], where KSH is the Spitzer Harm （SH） electron conductivity and f（T） interprets the preheating tongue effect in the cold plasma ahead of the ablative front [Phys. Rev. E 65 （2002） 57401], is introduced in simulations. The simulation results of the nonlinear evolution of the ARTI are in general agreement with the experiment results. It is found that two factors, i.e., the suppressing of ablative Kelvin Helmholtz instability （AKHI） and the heat flow cone in the spike tips, contribute to the formation of jet-like long spike in the nonlinear evolution of the ARTI.展开更多
We further develop the lattice Boltzmann （LB） model [Physica A 382 （2007） 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook （BGK） collision term in the LB equation, which...We further develop the lattice Boltzmann （LB） model [Physica A 382 （2007） 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook （BGK） collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference （FLFD） scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including （i） The thermal Couette flow; （ii） One- and two-dlmenslonal FLiemann problems. Good agreements are obtained between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems.展开更多
The aims of the present paper are threefold. First, we further study the fast Fourier transform thermal lattice Boltzmann （FFT-TLB） model for van der Waals （VDW） fluids proposed in Phys. Rev. E, 2011, 84（4）： 04...The aims of the present paper are threefold. First, we further study the fast Fourier transform thermal lattice Boltzmann （FFT-TLB） model for van der Waals （VDW） fluids proposed in Phys. Rev. E, 2011, 84（4）： 046715. We analyze the merits of the FFT approach over the traditional finite difference scheme and investigate the effects of smoothing factors on accuracy and stability in detail. Second, we incorporate the VDW equation of state with flexible parameters into the FFT- TLB model. As a result, the revised model may be used to handle multiphase flows with various critical densities and temperatures. Third, we design appropriate boundary conditions for systems with solid walls. These improvements, from the views of numerics and physics, significantly extend the application scope of the model in science and engineering.展开更多
基金Supported by the Science Foundations of Laboratory of Computational Physics and China Academy of Engineering Physics under Grant Nos.2009A0102005 and 2011A0201002National Natural Science Foundation of China under Grant Nos.11075021,91130020,11074300, and 10874242+2 种基金Fundamental research funds for the central university under Grant No.2010YS03Technology Support Program of LangFang under Grant Nos.2010011030,201101118/21/23/24Teaching and Research Foundation of NCIAE under Grant Nos.JY-2011-027-Y and JY-2011-028-Y
文摘We present an improved lattice Boltzmann(LB) model for thermal liquid-vapor system.In the new model,the Windowed Fast Fourier Transform(WFFT) and its inverse are used to calculate both the convection term and the external force term of the LB equation.By adopting the WFFT scheme,Gibbs oscillations can be damped effectively in unsmooth regions while high resolution feature of the spectral method can be retained in smooth regions.As a result,spatial discretization errors are dramatically decreased,conservation of the total energy is much better preserved,and the spurious velocities near the liquid-vapor interface are significantly reduced.The high resolution,together with the low complexity of the WFFT approach,endows the proposed method with considerable potential for studying a wide class of problems in the field of multiphase flows.
基金Support by the Science Foundations of Laboratory of Computational Physics,Science Foundation of China Academy of Engineering Physics under Grant Nos.2009A0102005,2009B0101012National Basic Research Program of China under Grant No.2007CB815105National Natural Science Foundation of China under Grant Nos.11074300,11075021,and 11071024
文摘The aims of the present paper are twofold. At first, we further study the Multiple-Relaxation-Time （MRT） Lattice Boltzmann （LB） model proposed in [Europhys. Lett. 90 （2010） 54003]. We discuss the reason why the Gram Schmidt orthogonalization procedure is not needed in the construction of transformation matrix M; point out a reason why the Kataoka-Tsutahara model [Phys. Rev. E 69 （2004） 035701（R）] is only valid in subsonic flows. The yon Neumann stability analysis is performed. Secondly, we carry out a preliminary quantitative study on the Richtmyer- Meshkov instability using the proposed MRT LB model. When a shock wave travels from a light medium to a heavy one, the simulated growth rate is in qualitative agreement with the perturbation model by Zhang-Sohn. It is about half of the predicted value by the impulsive model and is closer to the experimental result. When the shock wave travels from a heavy medium to a light one, our simulation results are also consistent with physical analysis.
基金Supported by the National Basic Research Program of China under Grant No 2007CB815103, and the National Natural Science Foundation of China under Grant Nos 10935003, 10775020 and 11075024.
文摘We report the formation of jet-like long spike in the nonlinear evolution of the ablative Rayleigh-Taylor instability （ARTI） experiments by numerical simulations. A preheating model k（T） = KSH[1＋f（T）], where KSH is the Spitzer Harm （SH） electron conductivity and f（T） interprets the preheating tongue effect in the cold plasma ahead of the ablative front [Phys. Rev. E 65 （2002） 57401], is introduced in simulations. The simulation results of the nonlinear evolution of the ARTI are in general agreement with the experiment results. It is found that two factors, i.e., the suppressing of ablative Kelvin Helmholtz instability （AKHI） and the heat flow cone in the spike tips, contribute to the formation of jet-like long spike in the nonlinear evolution of the ARTI.
基金Supported by the Science Foundations of LCP and CAEP under Grant Nos. 2009A0102005 and 2009B0101012National Natural Science Foundation of China under Grant Nos. 11075021, 11074300, and 11074303+3 种基金National Basic Research Program (973 Program) under Grant No. 2007CB815105Fundamental Research Funds for the Central University under Grant No. 2010YS03Technology Support Program of LangFang under Grant Nos. 2010011029/30/31Science Foundation of NCIAE under Grant No. 2008-ky-13
文摘We further develop the lattice Boltzmann （LB） model [Physica A 382 （2007） 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook （BGK） collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference （FLFD） scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including （i） The thermal Couette flow; （ii） One- and two-dlmenslonal FLiemann problems. Good agreements are obtained between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems.
文摘The aims of the present paper are threefold. First, we further study the fast Fourier transform thermal lattice Boltzmann （FFT-TLB） model for van der Waals （VDW） fluids proposed in Phys. Rev. E, 2011, 84（4）： 046715. We analyze the merits of the FFT approach over the traditional finite difference scheme and investigate the effects of smoothing factors on accuracy and stability in detail. Second, we incorporate the VDW equation of state with flexible parameters into the FFT- TLB model. As a result, the revised model may be used to handle multiphase flows with various critical densities and temperatures. Third, we design appropriate boundary conditions for systems with solid walls. These improvements, from the views of numerics and physics, significantly extend the application scope of the model in science and engineering.
