Sharif Digital Repository

In the present study, the preconditioned incompressible Navier-Stokes equations with the artificial compressibility method formulated in the generalized curvilinear coordinates are numerically solved by using a high-order compact finite-difference scheme for accurately and efficiently computing the incompressible flows in a wide range of Reynolds numbers. A fourth-order compact finite-difference scheme is utilized to accurately discretize the spatial derivative terms of the governing equations, and the time integration is carried out based on the dual time-stepping method. The capability of the proposed solution methodology for the computations of the steady and unsteady incompressible... 

The Eulerian methods are susceptible to generate the nonphysical spurious currents in the multiphase flow simulations near the interfaces. This paper presents a new Eulerian method to accurately simulate the velocity fields, especially near the multiphase flow interfaces and prevents the numerical results from generating the nonphysical currents. A Eulerian central difference finite-volume scheme equipped with the suitable numerical dissipation terms is used to simulate incompressible multiphase flows. The interface is captured by Flux Corrected Transport-Volume of Fluid method (FCT-VOF). Increasing the accuracy near the sharp gradients, such as interface, the conservative form of... 

This study investigates GH2/LOX coaxial jet flame at trans- and supercritical conditions using the Reynolds averaged Navier-Stokes approach. Four two-equation-turbulence models, three real equation of states, two chemical mechanisms, and three different chamber pressures are examined. Predictions show good agreement with measurements qualitatively and quantitatively. Based on the results, the predictions of the Soave-Redlich-Kwong equation of state (EOS) are closer to the experiment, while the Aungier-Redlich-Kwong EOS has more deviation than the others. Moreover, the k-ω shear stress transport model has better performance than the other turbulence models. It is also found that the flow... 

Spontaneous imbibition is well-known to be one of the most effective processes of oil recovery in fractured reservoirs. However, the detailed pore-scale mechanisms of the counter-current imbibition process and the effects of different fluid/rock parameters on this phenomenon have not yet been deeply addressed. Thiswork presents the results of a newpore-level numerical study of counter-current spontaneous imbibition, using coupled Cahn-Hilliard phase field and Navier-Stokes equations, solved by a finite element method. A 2D fractured medium was constructed consisting of a nonhomogeneous porous matrix, in which the grains were represented by an equilateral triangular array of circles with... 

Numerical simulations of the turbulent flow around single and side-by-side piles at different spacing ratios (centre-to-centre distance to the pile diameter) with flow Reynolds number of 105 on the fixed flat-bed are presented. The calculations are performed using the computational fluid dynamics model, FLOW-3D, which solves the Navier–Stokes equations in three dimensions with a finite-volume method. The numerical results of time-averaged flow patterns around single and side-by-side piles are validated using the available experimental measurements. At the downstream of the single pile, dimensionless vortex shedding frequency (Strouhal number) is estimated as 0.22. The maximum values of bed... 

The present study numerically investigates the feasibility of using multiple dielectric barrier discharge (multi-DBD) plasma actuators as a novel approach for active flow control over a large horizontal axis wind turbine rotor. The National Renewable Energy Laboratory 5 MW offshore wind turbine is used as the baseline case. This turbine uses pitch control system to adjust the generated power above its rated wind speeds, but at lower speeds, this system remains inactive. In this paper, the operational condition speed is considered lower than the rated wind speed. The mathematical electro-static model is implemented to simulate the effects of plasma actuator on the external flow and the... 

Numerical analysis has been conducted to simulate and capture Buzz phenomenon in a supersonic mixed compression air intake. Buzz is an unsteady self-sustained phenomenon occurred in supersonic intakes, especially when operating its subcritical condition, during which the system of compression and shock waves oscillate and move upstream and downstream along the intake. An axisymmetric and unsteady numerical simulation was used to solve Navier-Stokes equations in combination with URANS SST k-ω turbulence model The simulations were performed at M=2.0 and at a specific subcritical point of the intake operation where buzz was detected experimentally. Results are compared with experimental... 

In this paper, numerical simulation is performed to investigate the effects of micro air jets on mixing of the micro hydrogen jet in a transverse supersonic flow. The fundamental flow feature of the interaction between an array of fuel and air jets is investigated in a Mach 4.0 crossflow with a fuel global equivalence ratio of 0.5. Parametric studies were conducted on the various air jet conditions by using the Reynolds-averaged Navier–Stokes equations with Menter's Shear Stress Transport (SST) turbulence model. Numerical study of eight streamwise transverse sonic fuel and air jets in a fully turbulent supersonic flow revealed an extremely complex feature of fuel and air jet interaction. The... 

Purpose - The purpose of this paper is to present a fast, economical and practical method for mathematical modeling of aerodynamic characteristics of rectangular wing in ground (WIG) effect. Design/methodology/approach - Reynolds averaged Navier-Stokes (RANS) equations were converted to Bernoulli equation by reasonable assumptions. Also, Helmbold's equation has been developed for calculation of the slope of wing lift coefficient in ground effect by defining equivalent aspect ratio (ARe). Comparison of present work results against the experimental results has shown good agreement. Findings - A practical mathematical modeling with lower computational time and higher accuracy was presented for... 

Numerical study of pollutant emissions (NO and CO) in a Jet Stirred Reactor (JSR) combustor for methane oxidation under Elevated Pressure Lean Premixed (EPLP) conditions is presented. A Detailed Flow-field Simplified Chemistry (DFSC) method, a low computational cost method, is employed for predicting NO and CO concentrations. Reynolds Averaged Navier Stokes (RANS) equations with species transport equations are solved. Improved-coefficient five-step global mechanisms derived from a new evolutionary-based approach were taken as combustion kinetics. For modeling turbulent flow field, Reynolds Stress Model (RSM), and for turbulence chemistry interactions, finite rate-Eddy dissipation model are... 

In this study, we simulate the flow and heat transfer during hot-wire anemometry and investigate its thermal behavior and physics using the Computational Fluid Dynamics (CFD) tool. In this regard, we use the finite-volume method and solve the compressible Navier-Stokes equations numerically in slightly non-continuum flow fields. We do not use any slip flow model to include the transitional flow physics in our simulations. Using the CFD method, we simulate the flow over hot-wire and evaluate the uncertainty of CFD in thermal simulation of hot-wire in low transitional flow regimes. The domain sizes and the mesh distributions are carefully chosen to avoid boundary condition error appearances.... 

In this research, the problems of adiabatic film-cooling the flat, convex, and concave surfaces are investigated numerically. Two different radii of curvature and one row of vertical injection holes are considered. The Navier-Stokes equations are solved using a fine nonuniform multiblock staggered curvilinear grid and the SIMPLE-based finite volume method. The blowing rates are 0.5 and 1.0 and the mainstream Reynolds number is 10,000. The obtained results indicated that at a low blowing ratio, the cooling effectiveness enhances over the convex surface and reduces over the concave surface compared to the flat surface case. In comparison with the low blowing ratio, the curvature effects at a... 

Recently, the mechanisms of natural undulatory locomotion of aquatic animal swimming have become one of the most significant issues for the researchers, swimmers and engineers. This study aims to elucidate and compare the propulsive vortical signature and performance of backward (negative undulation) and forward (positive undulation) travelling waves through a typical fishlike propulsor by a systematic numerical study. The numerical approach uses a pressure-based finite volume method solver to solve Navier–Stokes equations in an arbitrary Lagrangian–Eulerian framework domain containing a two-dimensional NACA 0012 foil moving with prescribed kinematics. Some of the important findings are: (1)... 

A novel reduced order formulation is proposed for the vibration analysis of conical shells containing stationary fluid. Hamiltonian approach is followed to obtain the governing equations of motion for the structure. Utilizing the Navier-Stokes equations and simplifying for irrotational, compressible and inviscid assumptions, the final fluid equation is obtained. A general solution based on the Galerkin method is proposed for the conical shell in vacuum. Several boundary conditions are investigated to show the capability of the proposed solution. A novel reduced order formulation based on the finite element method is developed for solution of the fluid equation. Static condensation technique... 

Nature presents a variety of propulsion, maneuvering, and stabilization mechanisms which can be inspired to design and construction of manmade vehicles and the devices involved in them, such as stabilizers or control surfaces. This study aims to elucidate and compare the propulsive vortical signature and performance of a foil in two important natural mechanisms: flapping and undulation. Navier–Stokes equations are solved in an ALE framework domain containing a 2D NACA 0012 foil moving with prescribed kinematics. All simulations are carried out using a pressure-based finite volume method solver. The results of time-averaged inline force versus Strouhal number (St) show that in a given... 

With the rise in oil price and population growth, renewable energies are assumed to be the main source of energies for the next generation. Wind as a natural, eco-friendly and renewable source of energy has been at the center of concentration for decades. Recently there has been some ideas regarding the self-regulating urban wind turbines. Researchers have shown that a proper enclosure increases the wind velocity and therefore more torques can be exerted on the rotors, therefore more power can be generated. These enclosures are light and cheap, therefore they are applicable and effective. In the present study, an enclosure was modified to increase the exerted torque with implementing a step... 

Understanding and investigation of a high-speed craft dynamics, their longitudinal dynamic instabilities in calm water and behavior in waves are of a great importance. Determination of motion equation coefficients will help to analyze the dynamics of these kinds of vessels and the factors affecting their dynamic stabilities. Therefore, it could be useful in controlling the vessel instabilities. The main purpose of this research is to determine the coefficients of longitudinal motions of a planing catamaran using computational fluid dynamics (CFD) and evaluating the interceptor effect on hydrodynamic coefficients of that, which is widely utilized in controlling the motions of a high-speed... 

Recently, the mechanisms of natural undulatory locomotion of aquatic animal swimming have become one of the most significant issues for the researchers, swimmers and engineers. This study aims to elucidate and compare the propulsive vortical signature and performance of backward (negative undulation) and forward (positive undulation) travelling waves through a typical fishlike propulsor by a systematic numerical study. The numerical approach uses a pressure-based finite volume method solver to solve Navier–Stokes equations in an arbitrary Lagrangian–Eulerian framework domain containing a two-dimensional NACA 0012 foil moving with prescribed kinematics. Some of the important findings are: (1)... 

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 axial turbine is one of the most challenging components of gas turbines for industrial and aerospace applications. With the ever-increasing requirement for high-aerodynamic performance blades, three-dimensional aerodynamic shape optimization is of great importance. In this research, the rear part of a gas turbine consisting of a one-stage axial turbine is optimized numerically. A useful optimization algorithm is presented to improve the efficiency and/or pressure ratio of the axial turbine with two different objective functions. The three-dimensional blade-shape optimization is employed to study the effects of the turbine stator and rotor lean and sweep angles on the turbine performance....