Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userLattice Boltzmann simulations of melting in a rectangular cavity heated locally from below at high Rayleigh number
83
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
This work presented a block triple-relaxation-time (B-TriRT) lattice Boltzmann model for simulating melting in a rectangular cavity heated from below at high Rayleigh (Ra) number (Ra = 108). The test of benchmark problem shows that present B-TriRT can dramatically reduce the numerical diffusion across the phase interface. In addition, the influences of the location of the heated region are investigated. The results indicate that the location of heated region plays an essential role in melting rate and the full melting occur earliest when the heated region is located in the middle region.
rate research
Read More
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.
The variational multiscale (VMS) formulation is used to develop residual-based VMS large eddy simulation (LES) models for Rayleigh-B{e}nard convection. The resulting model is a mixed model that incorporates the VMS model and an eddy viscosity model. The Wall-Adapting Local Eddy-viscosity (WALE) model is used as the eddy viscosity model in this work. The new LES models were implemented in the finite element code Drekar. Simulations are performed using continuous, piecewise linear finite elements. The simulations ranged from $Ra = 10^6$ to $Ra = 10^{14}$ and were conducted at $Pr = 1$ and $Pr = 7$. Two domains were considered: a two-dimensional domain of aspect ratio 2 with a fluid confined between two parallel plates and a three-dimensional cylinder of aspect ratio $1/4$. The Nusselt number from the VMS results is compared against three dimensional direct numerical simulations and experiments. In all cases, the VMS results are in good agreement with existing literature.
We show that, when a single relaxation time lattice Boltzmann algorithm is used to solve the hydrodynamic equations of a binary fluid for which the two components have different viscosities, strong spurious velocities in the steady state lead to incorrect results for the equilibrium contact angle. We identify the origins of these spurious currents, and demonstrate how the results can be greatly improved by using a lattice Boltzmann method based on a multiple-relaxation-time algorithm. By considering capillary filling we describe the dependence of the advancing contact angle on the interface velocity.
In the recent years the lattice Boltzmann (LB) methodology has been fruitfully extended to include the effects of thermal fluctuations. So far, all studied cases pertain equilibrium fluctuations, i.e. fluctuations with respect to an equilibrium background state. In this paper we take a step further and present results of fluctuating LB simulations of a binary mixture confined between two parallel walls in presence of a constant concentration gradient in the wall-to-wall direction. This is a paradigmatic set-up for the study of non-equilibrium (NE) fluctuations, i.e. fluctuations with respect to a non- equilibrium state. We analyze the dependence of the structure factors for the hydrodynamical fields on the wave vector $boldsymbol{q}$ in both the directions parallel and perpendicular to the walls, as well as the finite-size effects induced by confinement, highlighting the long-range ($sim |boldsymbol{q}|^{-4}$) nature of correlations in the NE framework. Results quantitatively agree with the predictions of fluctuating hydrodynamics. Moreover, in presence of a non-ideal (NI) equation of state of the mixture, we also observe that the (spatially homogeneous) average pressure changes, due to a genuinely new contribution triggered by the long-range nature of NE fluctuations. These NE pressure effects are studied at changing the system size and the concentration gradient. Taken all together, we argue that these findings are instrumental to boost the applicability of the fluctuating LB methodology in the framework of NE fluctuations, possibly in conjunction with experiments.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
