Discrete Boltzmann modeling of Rayleigh-Taylor instability in two-component compressible flows

Effects of Inclined Interface Angle on Compressible Rayleigh-Taylor Instability: A Numerical Study Based on the Discrete Boltzmann Method

Rayleigh-Taylor (RT) instability is a basic fluid interface instability that widely exists in nature and in the engineering field. To investigate the impact of the initial inclined interface on compressible RT instability, the two-component discrete Boltzmann method is employed. Both the thermodynamic non-equilibrium (TNE) and hydrodynamic non-equilibrium (HNE) effects are studied. It can be found that the global average density gradient in the horizontal direction, the non-organized energy fluxes, the global average non-equilibrium intensity and the proportion of the non-equilibrium region first increase and then reduce with time. However, the global average density gradient in the vertical direction and the non-organized moment fluxes first descend, then rise, and finally descend. Furthermore, the global average density gradient, the typical TNE intensity and the proportion of non-equilibrium region increase with increasing angle of the initial inclined interface. Physically, there are three competitive mechanisms: (1) As the perturbed interface elongates, the contact area between the two fluids expands, which results in an increasing gradient of macroscopic physical quantities and leads to a strengthening of the TNE effects. (2) Under the influence of viscosity, the perturbation pressure waves on both sides of the material interface decrease with time, which makes the gradient of the macroscopic physical quantity decrease, resulting in a weakening of the TNE strength. (3) Due to dissipation and/or mutual penetration of the two fluids, the gradient of macroscopic physical quantities gradually diminishes, resulting in a decrease in the intensity of the TNE.

Thermodynamic nonequilibrium effects in bubble coalescence: A discrete Boltzmann study

The thermodynamic nonequilibrium (TNE) effects in a coalescence process of two initially static bubbles under thermal conditions are investigated by a discrete Boltzmann model. The spatial distributions of the typical nonequilibrium quantity, i.e., nonorganized momentum fluxes (NOMFs), during evolutions are investigated in detail. The density-weighted statistical method is used to highlight the relationship between the TNE effects and the morphological and kinetics characteristics of bubble coalescence. The results show that the xx component and yy component of NOMFs are antisymmetrical; the xy component changes from an antisymmetric internal and external double quadrupole structure to an outer octupole structure during the coalescence process. Moreover, the evolution of the averaged xx component of NOMFs provides two characteristic instants, which divide the nonequilibrium process into three stages. The first instant, when the averaged xx component of the NOMFs reaches its first local minimum, corresponds to the moment when the mean coalescence speed gets the maximum, and at this time the ratio of minor and major axes is about 1/2. The second instant, when the averaged xx component of the NOMFs gets its second local maximum, corresponds to the moment when the ratio of minor and major axes becomes 1 for the first time. It is interesting to find that the three quantities, TNE intensity, acceleration of coalescence, and the slope of boundary length, show a high degree of correlation and attain their maxima simultaneously. The surface tension and the heat conduction accelerate the process of bubble coalescence, while the viscosity delays it. Both the surface tension and the viscosity enhance the global nonequilibrium intensity, whereas the heat conduction restrains it. These TNE features and findings present some insights into the kinetics of bubble coalescence.

Discrete Boltzmann modeling of Rayleigh-Taylor instability: Effects of interfacial tension, viscosity, and heat conductivity

The two-dimensional Rayleigh-Taylor instability (RTI) in compressible flow with intermolecular interactions is probed via the discrete Boltzmann method. The effects of interfacial tension, viscosity, and heat conduction are investigated. It is found that the influences of interfacial tension on the perturbation amplitude, bubble velocity, and two kinds of entropy production rates all show differences at different stages of RTI evolution. It inhibits the RTI evolution at the bubble acceleration stage, while at the asymptotic velocity stage, it first promotes and then inhibits the RTI evolution. Viscosity and heat conduction inhibit the RTI evolution. Viscosity shows a suppressive effect on the entropy generation rate related to heat flow at the early stage but a first promotive and then suppressive effect on the entropy generation rate related to heat flow at a later stage. Heat conduction shows a promotive effect on the entropy generation rate related to heat flow at an early stage. Still, it offers a first promotive and then suppressive effect on the entropy generation rate related to heat flow at a later stage. By introducing the morphological boundary length, we find that the stage of exponential growth of the interface length with time corresponds to the bubble acceleration stage. The first maximum point of the interface length change rate and the first maximum point of the change rate of the entropy generation rate related to viscous stress can be used as a criterion for RTI to enter the asymptotic velocity stage.

Multiple-relaxation-time finite-difference lattice Boltzmann model for the nonlinear convection-diffusion equation

In this paper, a multiple-relaxation-time finite-difference lattice Boltzmann method (MRT-FDLBM) is developed for the nonlinear convection-diffusion equation (NCDE). Through designing the equilibrium distribution function and the source term properly, the NCDE can be recovered exactly from MRT-FDLBM. We also conduct the von Neumann stability analysis on the present MRT-FDLBM and its special case, i.e., single-relaxation-time finite-difference lattice Boltzmann method (SRT-FDLBM). Then, a simplified version of MRT-FDLBM (SMRT-FDLBM) is also proposed, which can save about 15% computational cost. In addition, a series of real and complex-value NCDEs, including the isotropic convection-diffusion equation, Burgers-Fisher equation, sine-Gordon equation, heat-conduction equation, and Schrödinger equation, are used to test the performance of MRT-FDLBM. The results show that both MRT-FDLBM and SMRT-FDLBM have second-order convergence rates in space and time. Finally, the stability and accuracy of five different models are compared, including the MRT-FDLBM, SMRT-FDLBM, SRT-FDLBM, the previous finite-difference lattice Boltzmann method [H. Wang, B. Shi et al., Appl. Math. Comput. 309, 334 (2017)10.1016/j.amc.2017.04.015], and the lattice Boltzmann method (LBM). The stability tests show that the sequence of stability from high to low is as follows: MRT-FDLBM, SMRT-FDLBM, SRT-FDLBM, the previous finite-difference lattice Boltzmann method, and LBM. In most of the precision test results, it is found that the order from high to low of precision is MRT-FDLBM, SMRT-FDLBM, SRT-FDLBM, and the previous finite-difference lattice Boltzmann method.

Multiple-relaxation-time discrete Boltzmann modeling of multicomponent mixture with nonequilibrium effects

A multiple-relaxation-time discrete Boltzmann model (DBM) is proposed for multicomponent mixtures, where compressible, hydrodynamic, and thermodynamic nonequilibrium effects are taken into account. It allows the specific heat ratio and the Prandtl number to be adjustable, and is suitable for both low and high speed fluid flows. From the physical side, besides being consistent with the multicomponent Navier-Stokes equations, Fick's law, and Stefan-Maxwell diffusion equation in the hydrodynamic limit, the DBM provides more kinetic information about the nonequilibrium effects. The physical capability of DBM to describe the nonequilibrium flows, beyond the Navier-Stokes representation, enables the study of the entropy production mechanism in complex flows, especially in multicomponent mixtures. Moreover, the current kinetic model is employed to investigate nonequilibrium behaviors of the compressible Kelvin-Helmholtz instability (KHI). The entropy of mixing, the mixing area, the mixing width, the kinetic and internal energies, and the maximum and minimum temperatures are investigated during the dynamic KHI process. It is found that the mixing degree and fluid flow are similar in the KHI process for cases with various thermal conductivity and initial temperature configurations, while the maximum and minimum temperatures show different trends in cases with or without initial temperature gradients. Physically, both heat conduction and temperature exert slight influences on the formation and evolution of the KHI morphological structure.

Lattice-Gas-Automaton Modeling of Income Distribution

A simple and effective lattice-gas-automaton (LGA) economic model is proposed for the income distribution. It consists of four stages: random propagation, economic transaction, income tax, and charity. Two types of discrete models are introduced: two-dimensional four-neighbor model (D2N4) and D2N8. For the former, an agent either remains motionless or travels to one of its four neighboring empty sites randomly. For the latter, the agent may travel to one of its nearest four sites or the four diagonal sites. Afterwards, an economic transaction takes place randomly when two agents are located in the nearest (plus the diagonal) neighboring sites for the D2N4 (D2N8). During the exchange, the Matthew effect could be taken into account in the way that the rich own a higher probability of earning money than the poor. Moreover, two kinds of income tax models are incorporated. One is the detailed taxable income brackets and rates, and the other is a simplified tax model based on a fitting power function. Meanwhile, charity is considered with the assumption that a richer agent donates a part of his income to charity with a certain probability. Finally, the LGA economic model is validated by using two kinds of benchmarks. One is the income distributions of individual agents and two-earner families in a free market. The other is the shares of total income in the USA and UK, respectively. Besides, impacts of the Matthew effect, income tax and charity upon the redistribution of income are investigated. It is confirmed that the model has the potential to offer valuable references for formulating financial laws and regulations.

Knudsen Number Effects on Two-Dimensional Rayleigh-Taylor Instability in Compressible Fluid: Based on a Discrete Boltzmann Method

Based on the framework of our previous work [H.L. Lai et al., Phys. Rev. E, 94, 023106 (2016)], we continue to study the effects of Knudsen number on two-dimensional Rayleigh–Taylor (RT) instability in compressible fluid via the discrete Boltzmann method. It is found that the Knudsen number effects strongly inhibit the RT instability but always enormously strengthen both the global hydrodynamic non-equilibrium (HNE) and thermodynamic non-equilibrium (TNE) effects. Moreover, when Knudsen number increases, the Kelvin–Helmholtz instability induced by the development of the RT instability is difficult to sufficiently develop in the later stage. Different from the traditional computational fluid dynamics, the discrete Boltzmann method further presents a wealth of non-equilibrium information. Specifically, the two-dimensional TNE quantities demonstrate that, far from the disturbance interface, the value of TNE strength is basically zero; the TNE effects are mainly concentrated on both sides of the interface, which is closely related to the gradient of macroscopic quantities. The global TNE first decreases then increases with evolution. The relevant physical mechanisms are analyzed and discussed.

Mesoscale modelling of miscible and immiscible multicomponent fluids

A mesoscopic simulation method based on the integration of dissipative particle dynamics (DPD), smoothed particle hydrodynamics (SPH) and computational thermodynamics (CT) has been developed. The kinetic behaviours of miscible and immiscible fluids were investigated. The interaction force between multicomponent mesoscopic particles is derived from the system free energy. The diffusivity of the components in non-ideal solution is determined by the chemical potential. The proposed method provides convincing predictions to the effects of convection, diffusion and microscopic interaction on the non-equilibrium evolution of engineering fluids, and demonstrates a potential to simulate more complicated phenomena in materials processing.

Phase-field method based on discrete unified gas-kinetic scheme for large-density-ratio two-phase flows

In this paper, a phase-field method under the framework of discrete unified gas-kinetic scheme (DUGKS) for incompressible multiphase fluid flows is proposed. Two kinetic models are constructed to solve the conservative Allen-Cahn equation that accounts for the interface behavior and the incompressible hydrodynamic equations that govern the flow field, respectively. With a truncated equilibrium distribution function as well as a temporal derivative added to the source term, the macroscopic governing equations can be exactly recovered from the kinetic models through the Chapman-Enskog analysis. Calculation of source terms involving high-order derivatives existed in the quasi-incompressible model is simplified. A series of benchmark cases including four interface-capturing tests and four binary flow tests are carried out. Results compared to that of the lattice Boltzmann method (LBM) have been obtained. A convergence rate of second order can be guaranteed in the test of interface diagonal translation. The capability of the present method to track the interface that undergoes a severe deformation has been verified. Stationary bubble and spinodal decomposition problems, both with a density ratio as high as 1000, are conducted and reliable solutions have been provided. The layered Poiseuille flow with a large viscosity ratio is simulated and numerical results agree well with the analytical solutions. Variation of positions of the bubble front and spike tip during the evolution of Rayleigh-Taylor instability has been predicted precisely. However, the detailed depiction of complicated interface patterns appearing during the evolution process is failed, which is mainly caused by the relatively large numerical dissipation of DUGKS compared to that of LBM. A high-order DUGKS is needed to overcome this problem.

Mesoscopic Simulation of the (2 + 1)-Dimensional Wave Equation with Nonlinear Damping and Source Terms Using the Lattice Boltzmann BGK Model

In this work, we develop a mesoscopic lattice Boltzmann Bhatnagar-Gross-Krook (BGK) model to solve (2 + 1)-dimensional wave equation with the nonlinear damping and source terms. Through the Chapman-Enskog multiscale expansion, the macroscopic governing evolution equation can be obtained accurately by choosing appropriate local equilibrium distribution functions. We validate the present mesoscopic model by some related issues where the exact solution is known. It turned out that the numerical solution is in very good agreement with exact one, which shows that the present mesoscopic model is pretty valid, and can be used to solve more similar nonlinear wave equations with nonlinear damping and source terms, and predict and enrich the internal mechanism of nonlinearity and complexity in nonlinear dynamic phenomenon.

Discrete Boltzmann trans-scale modeling of high-speed compressible flows

We present a general framework for constructing trans-scale discrete Boltzmann models (DBMs) for high-speed compressible flows ranging from continuum to transition regime. This is achieved by designing a higher-order discrete equilibrium distribution function that satisfies additional nonhydrodynamic kinetic moments. To characterize the thermodynamic nonequilibrium (TNE) effects and estimate the condition under which the DBMs at various levels should be used, two measures are presented: (i) the relative TNE strength, describing the relative strength of the (N+1)th order TNE effects to the Nth order one; (ii) the TNE discrepancy between DBM simulation and relevant theoretical analysis. Whether or not the higher-order TNE effects should be taken into account in the modeling and which level of DBM should be adopted is best described by the relative TNE intensity and/or the discrepancy rather than by the value of the Knudsen number. As a model example, a two-dimensional DBM with 26 discrete velocities at Burnett level is formulated, verified, and validated.