Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userThree-dimensional discrete Boltzmann models for compressible flows in and out of equilibrium
130
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present a series of three-dimensional discrete Boltzmann (DB) models for compressible flows in and out of equilibrium. The key formulating technique is the construction of discrete equilibrium distribution function through inversely solving the kinetic moment relations that it satisfies. The crucial physical requirement is that all the used kinetic moment relations must be consistent with the non-equilibrium statistical mechanics. The necessity of such a kinetic model is that, with increasing the complexity of flows, the dynamical characterization of non-equilibrium state and the understanding of the constitutive relations need higher order kinetic moments and their evolution. The DB models at the Euler and Navier-Stokes levels proposed by this scheme are validated by several well-known benchmarks, ranging from one-dimension to three-dimension. Particularly, when the local Mach number, temperature ratio, and pressure ratio are as large as $10^2$, $10^4$, and $10^5$, respectively, the simulation results are still in excellent agreement with the Riemann solutions. How to model deeper thermodynamic non-equilibrium flows by DB is indicated. Via the DB method, it convenient to simulate nonequilibrium flows without knowing exact form of the hydrodynamic equations.
rate research
Read More
A two-fluid Discrete Boltzmann Model(DBM) for compressible flows based on Ellipsoidal Statistical Bhatnagar-Gross-Krook(ES-BGK) is presented. The model has flexible Prandtl number or specific heat ratio. Mathematically, the model is composed of two coupled Discrete Boltzmann Equations(DBE). Each DBE describes one component of the fluid. Physically, the model is equivalent to a macroscopic fluid model based on Navier-Stokes(NS) equations, and supplemented by a coarse-grained model for thermodynamic non-equilibrium behaviors. To obtain a flexible Prandtl number, a coefficient is introduced in the ellipsoidal statistical distribution function to control the viscosity. To obtain a flexible specific heat ratio, a parameter is introduced in the energy kinetic moments to control the extra degree of freedom. For binary mixture, the correspondence between the macroscopic fluid model and the DBM may be several-to-one. Five typical benchmark tests are used to verify and validate the model. Some interesting non-equilibrium results, which are not available in the NS model or the single-fluid DBM, are presented.
We present numerical simulations of three-dimensional thermal convective flows in a cubic cell at high Rayleigh number using thermal lattice Boltzmann (LB) method. The thermal LB model is based on double distribution function approach, which consists of a D3Q19 model for the Navier-Stokes equations to simulate fluid flows and a D3Q7 model for the convection-diffusion equation to simulate heat transfer. Relaxation parameters are adjusted to achieve the isotropy of the fourth-order error term in the thermal LB model. Two types of thermal convective flows are considered: one is laminar thermal convection in side-heated convection cell, which is heated from one vertical side and cooled from the other vertical side; while the other is turbulent thermal convection in Rayleigh-Benard convection cell, which is heated from the bottom and cooled from the top. In side-heated convection cell, steady results of hydrodynamic quantities and Nusselt numbers are presented at Rayleigh numbers of $10^6$ and $10^7$, and Prandtl number of 0.71, where the mesh sizes are up to $257^3$; in Rayleigh-Benard convection cell, statistical averaged results of Reynolds and Nusselt numbers, as well as kinetic and thermal energy dissipation rates are presented at Rayleigh numbers of $10^6$, $3times 10^6$, and $10^7$, and Prandtl numbers of 0.7 and 7, where the nodes within thermal boundary layer are around 8. Compared with existing benchmark data obtained by other methods, the present LB model can give consistent results.
In a recent paper, Liu, Zhu and Wu (2015, {it J. Fluid Mech.} {bf 784}: 304) present a force theory for a body in a two-dimensional, viscous, compressible and steady flow. In this companion paper we do the same for three-dimensional flow. Using the fundamental solution of the linearized Navier-Stokes equations, we improve the force formula for incompressible flow originally derived by Goldstein in 1931 and summarized by Milne-Thomson in 1968, both being far from complete, to its perfect final form, which is further proved to be universally true from subsonic to supersonic flows. We call this result the textit{unified force theorem}, which states that the forces are always determined by the vector circulation $pGamma_phi$ of longitudinal velocity and the scalar inflow $Q_psi$ of transverse velocity. Since this theorem is not directly observable either experimentally or computationally, a testable version is also derived, which, however, holds only in the linear far field. We name this version the textit{testable unified force formula}. After that, a general principle to increase the lift-drag ratio is proposed.
Rayleigh-Taylor (RT) instability widely exists in nature and engineering fields. How to better understand the physical mechanism of RT instability is of great theoretical significance and practical value. At present, abundant results of RT instability have been obtained by traditional macroscopic methods. However, research on the thermodynamic non-equilibrium (TNE) effects in the process of system evolution is relatively scarce. In this paper, the discrete Boltzmann method based on non-equilibrium statistical physics is utilized to study the effects of the specific heat ratio on compressible RT instability. The evolution process of the compressible RT system with different specific heat ratios can be analyzed by the temperature gradient and the proportion of the non-equilibrium region. Firstly, as a result of the competition between the macroscopic magnitude gradient and the non-equilibrium region, the average TNE intensity first increases and then reduces, and it increases with the specific heat ratio decreasing; the specific heat ratio has the same effect on the global strength of the viscous stress tensor. Secondly, the moment when the total temperature gradient in y direction deviates from the fixed value can be regarded as a physical criterion for judging the formation of the vortex structure. Thirdly, under the competition between the temperature gradients and the contact area of the two fluids, the average intensity of the non-equilibrium quantity related to the heat flux shows diversity, and the influence of the specific heat ratio is also quite remarkable.
A one-dimensional discrete Boltzmann model for detonation simulation is presented. Instead of numerical solving Navier-Stokes equations, this model obtains the information of flow field through numerical solving specially discretized Boltzmann equation. Several classical benchmarks including Sod shock wave tube, Colella explosion problem, and one-dimensional self-sustainable stable detonation are simulated to validate the new model. Based on the new model, the influence of negative temperature coefficient of reaction rate on detonation is further investigated. It is found that an abnormal detonation with two wave heads periodically appears under negative temperature coefficient condition. The causes of the abnormal detonation are analyzed. One typical cycle of the periodic abnormal detonation and its development process are discussed.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
