Free-Energy-Based Discrete Unified Gas Kinetic Scheme for van der Waals Fluid

Unlocking the Key to Accelerating Convergence in the Discrete Velocity Method for Flows in the Near Continuous/Continuous Flow Regimes

How to improve the computational efficiency of flow field simulations around irregular objects in near-continuum and continuum flow regimes has always been a challenge in the aerospace re-entry process. The discrete velocity method (DVM) is a commonly used algorithm for the discretized solutions of the Boltzmann-BGK model equation. However, the discretization of both physical and molecular velocity spaces in DVM can result in significant computational costs. This paper focuses on unlocking the key to accelerate the convergence in DVM calculations, thereby reducing the computational burden. Three versions of DVM are investigated: the semi-implicit DVM (DVM-I), fully implicit DVM (DVM-II), and fully implicit DVM with an inner iteration of the macroscopic governing equation (DVM-III). In order to achieve full implicit discretization of the collision term in the Boltzmann-BGK equation, it is necessary to solve the corresponding macroscopic governing equation in DVM-II and DVM-III. In DVM-III, an inner iterative process of the macroscopic governing equation is employed between two adjacent DVM steps, enabling a more accurate prediction of the equilibrium state for the full implicit discretization of the collision term. Fortunately, the computational cost of solving the macroscopic governing equation is significantly lower than that of the Boltzmann-BGK equation. This is primarily due to the smaller number of conservative variables in the macroscopic governing equation compared to the discrete velocity distribution functions in the Boltzmann-BGK equation. Our findings demonstrate that the fully implicit discretization of the collision term in the Boltzmann-BGK equation can accelerate DVM calculations by one order of magnitude in continuum and near-continuum flow regimes. Furthermore, the introduction of the inner iteration of the macroscopic governing equation provides an additional 1-2 orders of magnitude acceleration. Such advancements hold promise in providing a computational approach for simulating flows around irregular objects in near-space environments.

Diffuse interface model for a single-component liquid-vapor system

We elucidate the theoretical relationships among fundamental physical concepts that are involved in the diffuse interface modeling for an isothermal single-component liquid-vapor system, which cover both the equation of state (EOS) and the surface tension force. As an example, a flat surface at equilibrium is discussed both theoretically and numerically by using two different approaches. Particularly, the force structure in the transition region is clearly presented, which demonstrates that the capillary contributions due to the density gradients can suppress the mechanical instability of the thermodynamic pressure and lead to constant hydrodynamic pressure (and chemical potential). Then, by comparing with the van der Waals (vdW) EOS for a flat interface at equilibrium, it is shown that applying the double-well approximation can give qualitative predictions for relatively high density ratio (ρ_{l}/ρ_{g}=7.784) and satisfactory results for relatively low density ratio (ρ_{l}/ρ_{g}=1.774). The main cause for this observation is attributed to the nonlinear variation of the generalized coefficient function in the double-well formulation at different density ratios. In addition, for the latter case, we simulate a droplet impact on a hydrophilic wall by using a recently proposed well-balanced discrete unified gas kinetic scheme (WB-DUGKS), which justifies the applicability of the double-well approximation to complex interfacial dynamics in the low-density-ratio limit. Furthermore, the reason for the inconsistency between the coefficients of the mean-field force expressions in the existing literature is explained.

Kinetic Theory-Based Methods in Fluid Dynamics

Kinetic theory stems from the statistical mechanics established at the mesoscopic scale [...].

Movable and Focus-Tunable Lens Based on Electrically Controllable Liquid: A Lattice Boltzmann Study

Adjusting the focal length by changing the liquid interface of the liquid lens has become a potential method. In this paper, the lattice-Boltzmann-electrodynamic (LB-ED) method is used to numerically investigate the zooming process of a movable and focus-tunable electrowetting-on-dielectrics (EWOD) liquid lens by combining the LBM chemical potential model and the electrodynamic model. The LB method is used to solve the Navier-Stokes equation, and the Poisson-Boltzmann (PB) equation is introduced to solve the electric field distribution. The experimental results are consistent with the theoretical results of the Lippmann-Young equation. Through the simulation of a liquid lens zoom driven by EWOD, it is found that the lens changes from a convex lens to a concave lens with the voltage increases. The focal length change rate in the convex lens stage gradually increases with voltage. In the concave lens stage, the focal length change rate is opposite to that in the convex lens stage. During the zooming process, the low-viscosity liquid exhibits oscillation, and the high-viscosity liquid appears as overdamping. Additionally, methods were proposed to accelerate lens stabilization at low and high viscosities, achieving speed improvements of about 30% and 50%, respectively. Simulations of lens motion at different viscosities demonstrate that higher-viscosity liquids require higher voltages to achieve the same movement speed.