Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userValidation and application of the lattice Boltzmann algorithm for a turbulent immiscible Rayleigh-Taylor system
74
0
0.0
(
0
)
Added by
Hugo Saraiva Tavares
Publication date
2021
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
We develop a multicomponent lattice Boltzmann (LB) model for the 2D Rayleigh--Taylor turbulence with a Shan-Chen pseudopotential implemented on GPUs. In the immiscible case this method is able to accurately overcome the inherent numerical complexity caused by the complicated structure of the interface that appears in the fully developed turbulent regime. Accuracy of the LB model is tested both for early and late stages of instability. For the developed turbulent motion we analyze the balance between different terms describing variations of the kinetic and potential energies. Then, we analyze the role of interface in the energy balance, and also the effects of the vorticity induced by the interface in the energy dissipation. Statistical properties are compared for miscible and immiscible flows. Our results can also be considered as a first validation step to extend the application of LB model to 3D immiscible Rayleigh-Taylor turbulence.
rate research
Read More
We studied turbulence induced by the Rayleigh-Taylor (RT) instability for 2D immiscible two-component flows by using a multicomponent lattice Boltzmann method with a Shan-Chen pseudopotential implemented on GPUs. We compare our results with the extension to the 2D case of the phenomenological theory for immiscible 3D RT studied by Chertkov and collaborators ({it Physical Review E 71, 055301, 2005}). Furthermore, we compared the growth of the mixing layer, typical velocity, average density profiles and enstrophy with the equivalent case but for miscible two-component fluid. Both in the miscible and immiscible cases, the expected quadratic growth of the mixing layer and the linear growth of the typical velocity are observed with close long-time asymptotic prefactors but different initial transients. In the immiscible case, the enstrophy shows a tendency to grow like $propto t^{3/2}$, with the highest values of vorticity concentrated close to the interface. In addition, we investigate the evolution of the typical drop size and the behavior of the total length of the interface in the emulsion-like state, showing the existence of a power law behavior compatible with our phenomenological predictions. Our results can also be considered as a first validation step to extend the application of lattice Boltzmann tool to study the 3D immiscible case.
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.
We numerically investigate turbulent Rayleigh-Benard convection within two immiscible fluid layers, aiming to understand how the layer thickness and fluid properties affect the heat transfer (characterized by the Nusselt number $Nu$) in two-layer systems. Both two- and three-dimensional simulations are performed at fixed global Rayleigh number $Ra=10^8$, Prandtl number $Pr=4.38$, and Weber number $We=5$. We vary the relative thickness of the upper layer between $0.01 le alpha le 0.99$ and the thermal conductivity coefficient ratio of the two liquids between $0.1 le lambda_k le 10$. Two flow regimes are observed: In the first regime at $0.04lealphale0.96$, convective flows appear in both layers and $Nu$ is not sensitive to $alpha$. In the second regime at $alphale0.02$ or $alphage0.98$, convective flow only exists in the thicker layer, while the thinner one is dominated by pure conduction. In this regime, $Nu$ is sensitive to $alpha$. To predict $Nu$ in the system in which the two layers are separated by a unique interface, we apply the Grossmann-Lohse theory for both individual layers and impose heat flux conservation at the interface. Without introducing any free parameter, the predictions for $Nu$ and for the temperature at the interface well agree with our numerical results and previous experimental data.
Simulating inhomogeneous flows with different characteristic scales in different coordinate directions using the collide-and-stream based lattice Boltzmann methods (LBM) can be accomplished efficiently using rectangular lattice grids. We develop and investigate a new rectangular central moment LBM based on non-orthogonal moment basis (referred to as RC-LBM). The equilibria to which the central moments relax under collision in this approach are obtained from matching with those corresponding to the continuous Maxwell distribution. A Chapman-Enskog analysis is performed to derive the correction terms to the second order moment equilibria involving the grid aspect ratio and velocity gradients that restores the isotropy of the viscous stress tensor and eliminates the non-Galilean invariant cubic velocity terms of the resulting hydrodynamical equations. A special case of this rectangular formulation involving the raw moments (referred to as the RNR-LBM) is also constructed. The resulting schemes represent a considerable simplification, especially for the transformation matrices and isotropy corrections, and improvement over the existing MRT-LB schemes on rectangular lattice grids that use orthogonal moment basis. Numerical validation study of both the RC-LBM and RNR-LBM for a variety of benchmark flow problems are performed that show good accuracy at various grid aspect ratios. The ability of our proposed schemes to simulate flows using relatively lower grid aspect ratios than considered in prior rectangular LB approaches is demonstrated. Furthermore, simulations reveal the superior stability characteristics of the RC-LBM over RNR-LBM in handling shear flows at lower viscosities and/or higher characteristic velocities. In addition, computational advantages of using our rectangular LB formulation in lieu of that based on the square lattice is shown.
A new lattice Boltzmann model for multicomponent ideal gas mixtures is presented. The model development consists of two parts. First, a new kinetic model for Stefan- Maxwell diffusion amongst the species is proposed and realized as a lattice Boltzmann equation on the standard discrete velocity set. Second, a compressible lattice Boltzmann model for the momentum and energy of the mixture is established. Both parts are consistently coupled through mixture composition, momentum, pressure, energy and enthalpy whereby a passive scalar advection-diffusion coupling is obviated, unlike in previous approaches. The proposed model is realized on the standard three-dimensional lattices and is validated with a set of benchmarks highlighting various physical aspects of compressible mixtures. Stefan-Maxwell diffusion is tested against experiment and theory of uphill diffusion of argon and methane in a ternary mixture with hydrogen. The speed of sound is measured in various binary and ternary compositions. We further validate the Stefan-Maxwell diffusion coupling with hydrodynamics by simulating diffusion in opposed jets and the three-dimensional Kelvin-Helmholtz instability of shear layers in a two-component mixture. Apart from the multicomponent compressible mixture, the proposed lattice Boltzmann model also provides an extension of the lattice Boltzmann equation to the compressible flow regime on the standard three-dimensional lattice.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
