Sharif Digital Repository

Due to superior accuracy and stability of multiple relaxation time (MRT) collision operator over its single relaxation time (SRT) counterpart, new lifting relations are proposed here to construct single particle distribution functions for MRT-LBM from macroscopic variables. Using these lifting relations, a new hybrid FVM-LB method is presented (called Finite type-LB hybrid method), which is consistent with MRT-LBM. In this new hybrid method, single-particle distribution functions in MRT-LBM sub-domain boundaries are computed, using equilibrium and non-equilibrium moments. These moments are computed in Navier-Stokes/FVM sub-domain boundaries, using macroscopic variables and their derivatives.... 

A new procedure for simulating turbulent flow in three-dimensional arbitrary geometries is presented. Finite volume method using physical covariant velocities on a staggered grid arrangement was used in this investigation. This work is an extension of previous successful work to three-dimensional cases. The ability of the new algorithm was tested using a conventional two-equation turbulence model on a highly separated turbulent flow test case. The low Reynolds number k-ω turbulence model of Wilcox was utilized to evaluate its capability in modeling highly curved flows. Turbulent flow over a three-dimensional hill, which is appropriate in assessment of ability of turbulence models in... 

Film cooling of surfaces appears in many applications. For instance, it is one of the most effective methods to improve the efficiency of gas turbines. As a fundamental study, two different types of film cooling (slot and discrete holes injections) are numerically simulated here. A flat surface is used to model a small portion of a gas turbine blade. Incompressible, stationary, viscous, turbulent flow is assumed using the STAR-CD software with the standard k-ε model and a cell-centered finite volume method on a nonuniform structured grid. The jet flow Reynolds number, based on the jet's hydraulic diameter, is 4.7 × 103. The study of the injection angle and the velocity ratio shows that the... 

Blade film cooling is one of the best methods to improve efficiency of gas turbines. In this work, two different methods of film cooling, namely, slot injection and discrete hole injection have been numerically studied on a flat plate. Incompressible, stationary, viscous, turbulent flow has been simulated using the FLUENT CFD code with the standard k-s model. The study of injection angle and velocity ratio show that the optimum film cooling in both methods, occurs at the jet angle of 30° but with the velocity ratio of 1.5 for slot case and 0.5 for discrete hole case. The study of jet aspect ratio in discrete hole method, shows that stretching the hole in spanwise direction increases the film... 

Calculation of a three-dimensional turbulent flow of square jets injected perpendicularly into a crossflow is performed computationally. The jet-to-crossflow velocity ratios were selected to be 0.5, 1.0, and 1.5, and the jet Reynolds number was 4700 based on the jet diameter. We solved the Reynolds-averaged Navier-Stokes equations in the general form using the SIMPLE finite volume method over a nonuniform staggered grid. Our computational domain included the jet channel flows as well as the flow over a flat plate. For the turbulence modeling, the standard k-ε model with wall functions and the zonal (k-ε)/(k-ω) turbulence model (shear stress transport model) were used. The results of the two... 

Mixing processes are widely used in industries and are being studied by a number of scientists; even though, a fundamental insight of the subject is still lacking. In this paper, some basic thoughts are set on the analogy between fluid mixing phenomena and fractal geometry/ properties/ processes. To this end, an unstable horseshoe pattern, which leads to a continuous extending boundary with maximum stretching and folding, is created and repeated. As a result, self-similarity property along with instability is imposed throughout the region, enhancing the fluid mixing. The procedures are verified computationally using the techniques of capacity and information dimensions, which are measured as... 

In this research, effects of superhydrophobic nano-coating on frictional drag force have been investigated. The result of this study could be considered to be used as a method in applications concerned with fuel consumption reduction, less CO2 emission and environmental problems as well as speed increase; while, its significance can also be of great use in marine applications. A rotating disc apparatus was used as the experimental set-up to compare the frictional drag force on an aluminum disc with TiO2 superhydrophobic nano-coating and a smooth coatless aluminum disc. The superhydrophobic nano-coating was prepared using sol-gel method and was shown to be able to produce a contact angle of... 

The pulsating flow in the exhaust gas of a SI engine causes an unsteady flow at the inlet to the turbocharger turbine. In a four cylinder four stroke engine, the pulse frequency varies between 20 and 200 Hz. Three dimensional pulsating flows in a vane-less turbocharger turbine of a 1.7 liters SI engine are simulated numerically and validated experimentally. Simulations are done for 720 degree engine cycle at three engine speeds. The results are shown the inlet pulsating flow has significant effects on several turbine parameters especially the inlet total pressure, the reduced mass flow rate and the efficiency. The results show a very good agreement between the three-dimensional unsteady... 

A comprehensive study was performed to analyze the unsteady laminar flow characteristics around a porously covered, a fully porous, and a solid squared section cylinder located in the middle of a plane channel. In order to simulate fluid flow inside porous media and porous–fluid interface accurately (minimizing modeling error), the porous region was analyzed in pore scale, using LBM. Additionally, to minimize the LBM-related compressibility error through the porous region, a multi-block multiple relaxation time lattice Boltzmann method (MRT-LBM) was used. Also, to decrease CPU time, a Navier–Stokes flow solver, based on finite volume method and SIMPLE algorithm, was coupled with MRT-LBM to... 

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... 

In this work, two-dimensional laminar flow and heat transfer across a heated square cylinder, covered by a porous layer in a plane channel have been numerically investigated. The flow and thermal fields inside the porous layer were simulated using BrinkmaneForchmeyer extended Darcy model. Simulations were performed in different Reynolds numbers (Re = 60, 120, 160, and 200), porosities (ω = 0.7, 0.87, and 0.96), solid to fluid thermal conductivity ratios (λR = 10, 200, and 2000) and blockage ratios (BR = 0.5, 0.25 and 0.125). The effects of the mentioned parameters on pressure drop and heat transfer rate were investigated in detail. Also, the contribution of each side of the central squared... 

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... 

In this study, an arbitrary Lagrangian-Eulerian (ALE) approach is incorporated with a mixed finite-volume-element (FVE) method to establish a novel moving boundary method for simulating unsteady incompressible flow on non-stationary meshes. The method collects the advantages of both finite-volume and finite-element (FE) methods as well as the ALE approach in a unified algorithm. In this regard, the convection terms are treated at the cell faces using a physical-influence upwinding scheme, while the diffusion terms are treated using bilinear FE shape functions. On the other hand, the performance of ALE approach is improved by using the Laplace method to improve the hybrid grids, involving... 

A series of experiments was performed to investigate the interaction of an under-expanded axisymmetric supersonic jet exhausted from a flat plate with a high subsonic crossflow. The goal was to study the effect of boundary layer thickness (δ) and jet to freestream dynamic pressure ratio (J) on flow field pressure distributions. The resulting measurements upstream of the jet showed that with increasing boundary layer thickness, the magnitude of the pressure coefficient decreases, whereas downstream of the jet, the recovery of the back-pressure moved closer to the nozzle exit. Flow field measurements indicated that with increasing boundary layer thickness, the jet plume dissipation rate... 

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... 

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... 

Three-dimensional bubble dynamics in rotating flow under an accelerating field such as a centrifugal one is studied in this work. We employ the lattice Boltzmann method in two phase flows to simulate bubble dynamics for different Bond and Morton numbers of 0.1, 1, 10, and 100 and 0.001, 0.01, 0.1, 1, 10, and 100, respectively. Another dimensionless number named as dimensionless force, F∗, which is the ratio of buoyancy force to centripetal force is defined to explain the dynamics of the bubbles. In this work, we consider 5×10-7≤F∗≤5. The results show that bubbles in rotating flows have different kinds of motions such as spinning, rotating, and translating. Based on the ratios of the forces... 

Due to the widespread use of rarefied gas flow in micro-porous media in industrial and engineering problems, a pore-scale modeling of rarefied gas flow through two micro-porous media with fractal geometries is presented, using lattice Boltzmann method. For this purpose, square- and circular-based Sierpinski carpets with fractal geometries are selected due to their inherent behavior for real porous media. Diffusive reflection slip model is used and developed for these porous media through this study. With this respect, the planar Poiseuille flow is selected as a benchmark and validated with the literature. The effect of Knudsen number (Kn) on the permeability is investigated and compared in... 

Multiphase problems with high density ratios and complex interfaces deal with numerical instabilities and require accurate considerations for capturing the multiphase interfaces. An Incompressible Smoothed Particle Hydrodynamics (ISPH) scheme is presented to simulate such problems. In order to keep the present scheme simple and stable, well-established formulations are used for discretizing the spatial derivatives and a repulsive force is applied at the multiphase interface between particles of different fluids to maintain the interface sharpness. Special considerations are included to overcome the difficulties to model severe physical discontinuities at the interface and surface tension... 

Hybrid FVM-LBM schemes are developed in the past few years to use capabilities of both Navier-Stokes based finite volume method (FVM) and lattice Boltzmann method (LBM) to solve macro-meso multiscale problems. In this scheme, the major task is to develop some lifting relations that reconstruct distribution functions in LBM sub-domain from macroscopic variables and their derivatives. The macroscopic variables are computed using Navier-Stokes based FVM in macroscale sub-domain, while distribution functions are computed using LBM in mesoscale sub-domain. The pioneer works in this area used the single relaxation time (SRT) version of LBM. However, it is known that the numerical stability and...