Sharif Digital Repository

A high-order compact finite-difference lattice Boltzmann method (CFDLBM) is proposed and applied to accurately compute steady and unsteady incompressible flows. Herein, the spatial derivatives in the lattice Boltzmann equation are discretized by using the fourth-order compact FD scheme, and the temporal term is discretized with the fourth-order Runge-Kutta scheme to provide an accurate and efficient incompressible flow solver. A high-order spectral-type low-pass compact filter is used to stabilize the numerical solution. An iterative initialization procedure is presented and applied to generate consistent initial conditions for the simulation of unsteady flows. A sensitivity study is also... 

We extend a hybrid finite-volume-element (FVE) method to treat the laminar reacting flow in cylindrical coordinates considering the collocation of all chosen primitive variables. To approximate the advection fluxes at the cell faces, we use the upwind-biased physical influence scheme PIS and derive a few new extended expressions applicable in the cylindrical frame. These expressions are derived for both the Navier-Stokes and reactive flow governing equations, of which the latter expressions are considered novel in the finite-volume formulation. To validate our derived expressions, the current results are compared with the experimental data and other available numerical solutions. The results... 

In this work, we aim to study the effect of temperature gradients on the combined natural convection-radiation heat transfer problem in participating media. To impose this combined effect, we first solve the radiative transfer equation in an absorbing and emitting media. Then, we suitably add the radiation source terms to the energy equation and solve the fluid flow equations. Literature shows that many incompressible algorithms use the Boussinesq assumption to model the thermobuoyant force; however, the validity of this assumption is limited to cases with low temperature gradient distributions. Evidently, Boussinesq assumption would result in considerable errors in high temperature gradient... 

In this research, numerical optimization of the rear part of a gas turbine, consisting of a single stage axial turbine is carried out. Automated aerodynamic shape optimization is performed by coupling a CFD flow simulation code with the Genetic Algorithm. An effective multi-point optimization method to improve efficiency and/or pressure ratio of the axial turbine is performed. Some variations of optimization parameters such as lean and sweep angels of stator and rotor blades are accomplished. Furthermore, during the optimization process, three-dimensional and turbulent flow field is numerically investigated using a compressible Navier-Stokes solver. The gas turbine experimental... 

In this paper, the ball valve performance is numerically simulated using an unstructured CFD (Computational Fluid Dynamics) code based on the finite volume method. Navier-Stokes equations in addition to a transport equation for the vapor volume fraction were coupled in the RANS solver. Separation is modeled very well with a modification of turbulent viscosity. The results of CFD calculations of flow through a ball valve, based on the concept of experimental data, are described and analyzed. Comparison of the flow pattern at several opening angles is investigated. Pressure drop behind the ball valve and formation of the vortex flow downstream the valve section are also discussed. As the... 

This study describes a multidimensional 3D/lumped parameter (LP) model which contains appropriate inflow/outflow boundary conditions in order to model the entire human arterial trees. A new extensive LP model of the entire arterial network (48 arteries) was developed including the effect of vessel diameter tapering and the parameterization of resistance, conductor and inductor variables. A computer aided-design (CAD) algorithm was proposed to efficiently handle the coupling of two or more 3D models with the LP model, and substantially lessen the coupling processing time. Realistic boundary conditions and Navier-Stokes equations in healthy and stenosed models of carotid artery bifurcation... 

The present work is concerned with the steady incompressible flow through a parallel plate channel with stretching walls under an externally applied magnetic field. The governing continuity and Navier-Stokes equations are reduced to a fourth order nonlinear differential equation by using vorticity definition and similarity solution transformation. The obtained equations are solved by applying the analytical homotopy perturbation method (HPM). The method is called order of magnitude suggested for simplifying series solution to finite expression that is useful in engineering problems. The results are verified by comparing with numerical solutions and demonstrate a good accuracy of the obtained... 

Red Blood Cells (RBCs) are the main cells in human blood, with a main role in the mechanical properties of blood as a fluid. Several methods have been improved to simulate the mechanical behavior of RBC in micro-capillaries. Since, in microscopic scales, using discrete models is more preferred than continuum methods, the Moving Particle Semi-Implicit method (MPS), which is a recent innovative particle based method, can simulate micro-fluidic flows based on NavierStokes equations. Although, by recent developments, the MPS method has turned into a considerable tool for modeling blood flow in micro meter dimensions, some problems, such as a commitment to use small time step sizes, still... 

The present study focuses on the numerical simulation of cavitation around the NACA 0015. The unsteady behaviors of cavitation which have worthwhile applications are investigated. The cavitation patterns, velocity fields and frequency of the cavitating flow around hydrofoil is obtained. For multi phase simulation, single-fluid Navier-Stokes equations, along with a volume fraction transport equation, are employed. The bubble dynamics model is utilized to simulate phase change. SIMPLE algorithm is used for velocity and pressure computations. For discretization of equations the finite-volume approach written in body fitted curvilinear coordinates, on collocated grid, is used. In this study,... 

The present study focuses on the implementation of PISO algorithm to simulate cavitating flows. For simulation of unsteady behaviors of cavitation which have practical applications, the development of unsteady PISO algorithm based on the non-conservative approach is investigated. The effects of mixture compressibility are considered to improve accuracy of simulations. For multi-phase simulation, single-fluid Navier-Stokes equations, along with the volume fraction transport equation, are employed. The bubble dynamics model is utilized to simulate phase change. To prove capabilities of the developed PISO algorithm to simulate cavitating flows, unsteady simulation of cavitation around NACA0015... 

In this paper, the structure of a wall jet deflected by a baffle along with the trajectory of particles has been studied. This baffle is used to produce a stable deflected surface jet, thereby deflecting the high-velocity supercritical stream away from the bed to the surface. An elliptic relaxation turbulence model (ν̄2-f model) has been used to simulate this submerged flow. In recent years, the ν̄2- f turbulence model has become increasingly popular due to its ability to account for near-wall damping without use of damping functions. In addition, it has been proven that the ν̄2- f model is superior to other Reynolds-averaged Navier-Stokes (RANS) methods in many flows where complex flow... 

Landslide-generated waves (LGWs) are among natural hazards that have stimulated attentions and concerns of engineers and researchers during the past decades. At the same period, the application of numerical modeling has been progressively increased to assess, control, and manage the risks of such hazards. This paper represents an overview of numerical studies on LGWs to explore associated recent advances and future challenges. In this review, the main landslide events followed by an LGW hazard are scrutinized. The uncertainty regarding landslide characteristics and the lack of data concerning generated tsunami properties highlights the necessity of probabilistic analysis and numerical... 

Excess pressure drop induced by inertial effects limits the applicability of Darcy's law for modeling of fluid flow through porous media at high velocities. It is expected such additional pressure drop is influenced by pore/morphology of porous media. This work concerns with fundamental understanding of how throat curvature affects intrinsic properties of porous media at non-Darcy flow conditions using network modeling. Conical, parabolic, hyperbolic, and sinusoidal capillary ducts with three types of imposed anisotropy are used to construct the network in a more realistic manner. Solutions of one dimensional Navier-Stokes equation for incompressible fluid flow through converging/diverging... 

In this work, numerical investigation of a gas turbine model combustor (GTMC) was carried out using two different turbulence-chemistry interaction models: the EDC (eddy dissipation concept) and TPDF (transported probability density function). GTMC with good optical access for laser measurements provides a useful database for swirling CH4/Air diffusion flames at atmospheric pressure. Modeling was performed by solving Reynolds-averaged Navier-Stokes (RANS) and Reynolds stress model (RSM) equations for a two-dimensional (2D) axisymmetric computational domain accompanied by swirl and the combustion chamber was investigated for both reacting and nonreacting conditions. A detailed reduced... 

A computational tool is developed to capture the interaction of solid object with two-phase flow. The full two-dimensional Navier–Stokes equations are solved on a regular structured grid to resolve the flow field. The level set and the immersed boundary methods are used to capture the free surface of a fluid and a solid object, respectively. A two-step projection method along with Multi-Processing (OpenMP) is employed to solve the flow equations. The computational tool is verified based on numerical and experimental data with three scenarios: a cylinder falling into a rectangular domain due to gravity, transient vertical oscillation of a cylinder by releasing above its equilibrium position,... 

The first notable megawatt class wind turbine, which was the pioneer of improvement in the blade performance in large wind turbines, appeared in Vermont. Nowadays, modern wind turbines are using blades with multi-airfoils at different sections. In this study, in order to indicate the best airfoil profile for the optimum performance in different sections of a blade, five popular airfoils, including S8xx, FFA and AH series, were studied. On the large-scale profile, shear stress transport K–ω model was applied for the simulation of horizontal axis wind turbines for different wind speeds. The aerodynamic simulation was accomplished using computational fluid dynamic method, which in turn is based... 

The behavior of non-Newtonian fluids is considered as an important subject in micro scale and microfluidic flow researches. Because of the complexity and cost in the numerical works and the experimental set-ups in some instances, the analytical approach can be taken into account as a robust alternative tool to solve the non-Newtonian microfluidic flows in some special cases benefiting from a few simplified assumptions. In this work, we analyze the flow of two non-Newtonian fluids including the power-law and grade-fluid models in microchannels. For the grade-fluid, the stress tensors are defined considering the Rivlin-Ericksen tensor definitions. To avoid the complexities in the entrance... 

Kinetic theory based numerical scheme such as direct simulation Monte Carlo (DSMC) and information preservation (IP) schemes properly solve micro-nano flow problems in transition and free molecular regimes. However, the high computational cost of these methods encourages the researchers toward extending the applicability of the continuumbased equations beyond the slip flow regime. In addition to correct velocity profile, the continuum-based equations should predict accurate mass flow rate magnitude. The secondorder velocity slip models derived from the kinetic theory provide accurate velocity profiles up to Kn=0.5; however, they yield erroneous mass flow rate magnitudes because the basic... 

In this work, we present a level set method to simulate steady and unsteady mass transfer from a single droplet moving in a second phase fluid under buoyant force. We initially use level set to determine the interface between the two phases, where the shape of drop forms. Next, we extend this method to solve a unique mass transfer equation for the entire solution domain without considering the discontinuity appeared at the interface. We use a finite element method incorporated with the characteristic-based split (CBS) algorithm to implement axi-symmetric mass transfer equations on a stationary Eulerian grid. Of course, the convection-diffusion modeling of mass transfer is different from the... 

The current study investigates the conjugate heat transfer characteristics for laminar flow in backward facing step channel. All of the channel walls are insulated except the lower thick wall under a constant temperature. The upper wall includes a insulated obstacle perpendicular to flow direction. The effect of obstacle height and location on the fluid flow and heat transfer are numerically explored for the Reynolds number in the range of 10 ≤ Re ≤ 300. Incompressible Navier-Stokes and thermal energy equations are solved simultaneously in fluid region by the upwind compact finite difference scheme based on flux-difference splitting in conjunction with artificial compressibility method. In...