Sharif Digital Repository

Multiphase flows appear in many scientific and engineering applications. On the other hand, there has been a tremendous experimental and numerical investigation trying to understand complex fluid-fluid interfaces better. However, conventional CFD have certain limitations in complicated situations. Fortunately in the past two decades, new approaches, namely SPH and LBM (sometimes called meshless methods) have been developed. Since they are mesoscopic, they have been able to perform much better than conventional CFD. Especially, LBM is becoming more and more popular and has already been divided to many branches  

In this work, the implementation of a high-order compact finite-difference lattice Boltzmann method (CFDLBM) is performed in the generalized curvilinear coordinates to improve the computational efficiency of the solution algorithm to handle curved geometries with non-uniform grids. The incompressible form of the discrete Boltzmann equation with the Bhatnagar-Gross-Krook (BGK) approximation with the pressure as the independent dynamic variable is transformed into the generalized curvilinear coordinates. Herein, the spatial derivatives in the resulting lattice Boltzmann (LB) equation in the computational plane are discretized by using the fourth-order compact finite-difference scheme and the... 

Attempts to simulate compressible flows with moderate Mach number to relatively high ones using Lattice Boltzmann Method (LBM) have been made by numerous researchers in the recent decade. The stability of the LBM is a challenging problem in the simulation of compressible flows with different types of embedded discontinuities. The present study proposes an approach for simulation of inviscid flows by a compressible LB model in order to enhance the robustness using a combination of Essentially NonOscillatory (ENO) scheme and Shock-Detecting Sensor (SDS) procedure. A sensor is introduced with adjustable parameters which is active near the discontinuities and affects less on smooth regions. The... 

Several efforts have been performed to improve LBM defects related to its computational performance. In this work, a new algorithm has been introduced to reduce memory consumption. In the past, most LBM developers have not paid enough attention to retain LBM simplicity in their modified version, while it has been one of the main concerns in developing of the present algorithm. Note, there is also a deficiency in our new algorithm. Besides the memory reduction, because of high memory call back from the main memory, some computational efficiency reduction occurs. To overcome this difficulty, an optimization approach has been introduced, which has recovered this efficiency to the original... 

The search for the development of a reliable mathematical model for understanding bubble dynamics behavior is an ongoing endeavor. A long list of complex phenomena underlies the physics of this problem. In the past decades, the lattice Boltzmann method has emerged as a promising tool to address such complexities. In this regard, we have applied a 121-velocity multiphase lattice Boltzmann model to an asymmetric cluster of bubbles in an acoustic field. A problem as a benchmark is studied to check the consistency and applicability of the model. The problem of interest is to study the deformation and coalescence phenomena in bubble cluster dynamics, as well as the screening effect on an acoustic... 

This paper deals with the investigation of the laminar flow past an elliptical cylinder confined in a channel. In this paper, the Lattice Boltzmann (LB) method is used to simulate flow in two dimensions. The present LB method with the used boundary conditions is validated in simulations of the incompressible flow past a circular cylinder. The simulations are carried out in a range of condition, 0≤θ≤90 (angle of incidence), 5≤Re≤100, (Reynolds number) for AR=0.25,0.5 (aspect ratio). The effects of those parameters on the drag and lift coefficients and other flow characteristics of the cylinder are examined in detail. The results demonstrate that the drag and lift coefficients increase with... 

This paper deals with the investigation of the laminar flow past an elliptical cylinder confined in a channel. In this paper, the Lattice Boltzmann (LB) method is used to simulate flow in two dimensions. The present LB method with the used boundary conditions is validated in simulations of the incompressible flow past a circular cylinder. The simulations are carried out in a range of condition, 0 90 (angle of incidence), 5 Re 100, (Reynolds number) for AR=0.25,0.5 (aspect ratio). The effects of those parameters on the drag and lift coefficients and other flow characteristics of the cylinder are examined in detail. The results demonstrate that the drag and lift coefficients increase with the... 

A multi-component lattice Boltzmann scheme is used to investigate the dynamics of a wettability driven droplet within a microchannel. The driving force for the motion is created by a stepwise change in the wettability of the channel walls. Moreover, an analytical solution is developed for evaluation of the dynamics of the droplet inside the channel. The effects of various parameters such as the height of the channel, the wetting pattern of the channel walls, the viscosity and the density ratio on the dynamics are studied. Also, the effect of grooves of different sizes on the channel surfaces on the dynamics of the droplet is investigated for both hydrophilic and hydrophobic surfaces.... 

In this work, the Chebyshev collocation spectral lattice Boltzmann method is implemented in the generalized curvilinear coordinates to provide an accurate and efficient low-speed LB-based flow solver to be capable of handling curved geometries with non-uniform grids. The low-speed form of the D2Q9 and D3Q19 lattice Boltzmann equations is transformed into the generalized curvilinear coordinates and then the spatial derivatives in the resulting equations are discretized by using the Chebyshev collocation spectral method and the temporal term is discretized with the fourth-order Runge–Kutta scheme to provide an accurate and efficient low-speed flow solver. All boundary conditions are... 

In recent years, application of porous media is highlighted among researchers due to their wide range of usability in micro-scale problems, such as gas reservoirs, micro-filtering, heat exchangers, etc. With this respect, the accurate description of flow behavior using governing equations based on the continuum assumption is not valid since the mean free path is comparable to the characteristics length of the problem. For this purpose, a simple methodology for diffusion reflection boundary condition is developed and validated for two valuable benchmarks, namely micro-channel flow and fractal porous media, where the results were in good agreement with literature. Then, pore-scale simulation... 

Drop formation in cross-junction micro-channels is numerically studied using the lattice Boltzmann method with pseudo-potential model. To verify the simulation, the results are compared to previous numerical and experimental data. Furthermore, the effects of capillary number, flow rate ratio, contact angle, and viscosity ratio on the flow patterns, drop length, and interval between drops are investigated and highlighted. The results show that the drop forming process has different regimes, namely, jetting, drop, and squeezing regimes. Also, it is shown that increasing in the flow rate ratio in the squeezing regime causes increment in drop length and decrement in drops interval distance. On... 

The lattice Boltzmann simulation of natural convection fluid flow and heat transfer is performed. The double-MRT lattice Boltzmann method is employed. The D2Q5 and D2Q9 lattice models are used to simulate the temperature field and flow field, respectively. The configuration of the cavity is U-shaped, and the cavity is filled with CuO-water nanofluid. The KKL (Koo-Kleinstreuer-Li) model is utilized to estimate the thermal conductivity and dynamic viscosity; as well as, the nanoparticle shape effect on the thermo-physical properties of nanofluid is considered. The streamlines, temperature distribution, local entropy generation map, average Nusselt number, magnitude of total entropy generation... 

This paper seeks to make a study on flow control in two-dimensional square cavities having obstacles on their walls. The goal of using these passive controllers is to enhance mixing in an enclosed space. Lattice Boltzmann method is used to simulate the problem. Results are presented for various Reynolds numbers, 400≤Re≤4000 and different arrangements of tiny-obstacles with different heights. To give a perspective on the physics of this problem, time evolution of the flow is studied at Re= 1000. Then, the flow structure is studied for different Reynolds numbers. Findings show that the interaction of the main vortex with the tiny-obstacles inserted on the wall cavity changes the flow pattern... 

This paper seeks to make a study on flow control in two-dimensional square cavities having obstacles on their walls. The goal of using these passive controllers is to enhance mixing in an enclosed space. Lattice Boltzmann method is used to simulate the problem. Results are presented for various Reynolds numbers, 400≤Re≤4000 and different arrangements of tiny-obstacles with different heights. To give a perspective on the physics of this problem, time evolution of the flow is studied at Re= 1000. Then, the flow structure is studied for different Reynolds numbers. Findings show that the interaction of the main vortex with the tiny-obstacles inserted on the wall cavity changes the flow pattern... 

Individual drops are suitable tools to study the liquid-fluid interfacial properties. In this work, force-displacement equation and non-linear oscillations of a pendent drop are numerically investigated. The presented novel force-displacement function allows following the dynamics of a pendent drop and realizing its elastic behavior. The growth and detachment of drop, which is pending due to gravity from a capillary tip, is considered (assuming high density and high viscosity ratios and immiscible two-phase flows). Two-dimensional multi-relaxation time lattice Boltzmann method (MRT-LBM) was used to simulate growth, detachment, and oscillations of the drop using a conservative model for... 

We investigate thin liquid film coating of substrates with trenches via lattice Boltzmann method. The effects of different parameters such as the capillary number, the contact angles and geometric parameters on the results are studied. Our results indicate that for trench depths greater than a critical size coating of the trench is not successful. The critical depth increases by decreasing the capillary number. Moreover we find height, width, capillary number and contact angle under which the coating is successful. The results have been compared with the available results and very close agreement has been achieved  

In this research, the effect of using nanoparticles on the Nusselt number in the natural convection at high Rayleigh numbers has been investigated. A two-component multiple relaxation times lattice Boltzmann code has been developed to investigate the natural convection of nanofluids in a 2D square enclosure. Three main forces, namely, buoyancy, drag, and Brownian forces are included to consider the interactions of the components. The potential force is also added to the model for the interactions among nanoparticles. The effects of Ra number, volume fraction of nanoparticles (ϕ), and size of nanoparticles (λ) have been investigated for two nanofluids (CuO-water and TiO2-water). Five...