Accurately solving phase interface plays a great role in modeling an immiscible multiphase flow system. In this paper, we propose an accurate interface-capturing lattice Boltzmann method from the perspective of the modified Allen-Cahn equation (ACE). The modified ACE is built based on the commonly used conservative formulation via the relation between the signed-distance function and the order parameter also maintaining the mass-conserved characteristic.
View Article and Find Full Text PDFTwo-dimensional direct numerical simulations of partitioned thermal convection are performed using the thermal lattice Boltzmann method for the Rayleigh number () of 10 and the Prandtl number () of 7.02 (water). The influence of the partition walls on the thermal boundary layer is mainly focused on.
View Article and Find Full Text PDFIn this paper, the lattice Boltzmann method was used to simulate the cardiac flow in children with aseptal defect. The inner wall model of the heart was reconstructed from 210 computed tomography scans. By simulating and comparing the cardiac flow field, the pressure field, the blood oxygen content, and the distribution of entropy generation before and after an operation, the effects of septal defect on pulmonary hypertension(PH), cyanosis, and heart load were analyzed in detail.
View Article and Find Full Text PDFThe effect of rotation on small-scale characteristics and scaling law in the mixing zone of the three-dimensional turbulent Rayleigh-Taylor instability (RTI) is investigated by the lattice Boltzmann method at small Atwood number. The mixing zone width h(t), the root mean square of small scale fluctuation, the spectra, and the structure functions are obtained to analyze the rotating effect. We mainly focus on the process of the development of plumes and discuss the physical mechanism in the mixing zone in rotating and nonrotating systems.
View Article and Find Full Text PDFTime evolution features of kinetic and thermal entropy generation rates in turbulent Rayleigh-Bénard (RB) convection with mixed insulating and conducting boundary conditions at = 10 are numerically investigated using the lattice Boltzmann method. The state of flow gradually develops from laminar flow to full turbulent thermal convection motion, and further evolves from full turbulent thermal convection to dissipation flow in the process of turbulent energy transfer. It was seen that the viscous, thermal, and total entropy generation rates gradually increase in wide range of /τ < 32 with temporal evolution.
View Article and Find Full Text PDFIn this work, the temporal-spatial evolution of kinetic and thermal energy dissipation rates in three-dimensional (3D) turbulent Rayleigh-Taylor (RT) mixing are investigated numerically by the lattice Boltzmann method. The temperature fields, kinetic and thermal energy dissipation rates with temporal-spatial evolution, the probability density functions, the fractal dimension of mixing interface, spatial scaling law of structure function for the kinetic and the thermal energy dissipation rates in 3D space are analysed in detail to provide an improved physical understanding of the temporal-spatial dissipation-rate characteristic in the 3D turbulent Rayleigh-Taylor mixing zone. Our numerical results indicate that the kinetic and thermal energy dissipation rates are concentrated in areas with large gradients of velocity and temperature with temporal evolution, respectively, which is consistent with the theoretical assumption.
View Article and Find Full Text PDFEntropy (Basel)
September 2018
Entropy generation rates in two-dimensional Rayleigh-Taylor (RT) turbulence mixing are investigated by numerical calculation. We mainly focus on the behavior of thermal entropy generation and viscous entropy generation of global quantities with time evolution in Rayleigh-Taylor turbulence mixing. Our results mainly indicate that, with time evolution, the intense viscous entropy generation rate s u and the intense thermal entropy generation rate S θ occur in the large gradient of velocity and interfaces between hot and cold fluids in the RT mixing process.
View Article and Find Full Text PDFA numerical investigation has been carried out to understand the mechanism of the rotation effect on bifurcation and dual solutions in natural convection within a horizontal annulus. A thermal immersed boundary-lattice Boltzmann method was used to resolve the annular flow domain covered by a Cartesian mesh. The Rayleigh number based on the gap width is fixed at 10.
View Article and Find Full Text PDFStatistics of heat transfer in two-dimensional (2D) turbulent Rayleigh-Bénard (RB) convection for Pr=6,20,100 and 106 are investigated using the lattice Boltzmann method (LBM). Our results reveal that the large scale circulation is gradually broken up into small scale structures plumes with the increase of Pr, the large scale circulation disappears with increasing Pr, and a great deal of smaller thermal plumes vertically rise and fall from the bottom to top walls. It is further indicated that vertical motion of various plumes gradually plays main role with increasing Pr.
View Article and Find Full Text PDFEntropy generation in thermal convection with differentially discrete heat boundary conditions at various Rayleigh numbers () are numerically investigated using the lattice Boltzmann method. We mainly focused on the effects of and discrete heat boundary conditions on entropy generation in thermal convection according to the minimal entropy generation principle. The results showed that the presence of the discrete heat source at the bottom boundary promotes the transition to a substantial convection, and the viscous entropy generation rate () generally increases in magnitude at the central region of the channel with increasing .
View Article and Find Full Text PDFThe immersed boundary-lattice Boltzmann method (IB-LBM) was used to examine the motion and deformation of three elastic red blood cells (RBCs) during Poiseuille flow through constricted microchannels. The objective was to determine the effects of the degree of constriction and the Reynolds (Re) number of the flow on the physical characteristics of the RBCs. It was found that, with decreasing constriction ratio, the RBCs experienced greater forced deformation as they squeezed through the constriction area compared to at other parts of the microchannel.
View Article and Find Full Text PDFComput Math Methods Med
March 2017
The flow characteristics in the realistic human upper airway (HUA) with obstruction that resulted from pharyngeal collapse were numerically investigated. The 3D anatomically accurate HUA model was reconstructed from CT-scan images of a Chinese male patient (38 years, BMI 25.7).
View Article and Find Full Text PDF