Sharif Digital Repository

In recent years, requirement to increase projectile performance, leads to a great interest in high angle of attack Aerodynamics. Projectile Maneuverability is a concept which is defined by having the capability to perform at high angles of attack while maintaining attached flow over the tail, rather than a capability of performing a mission in minimum possible duration. Interaction between body, wing and tail vortices can delay vortex breakdown over the wing and as a result may enhance projectile maneuverability. Studies show that wing aspect ratio and tail location have a remarkable influence on the vortices interaction and thus projectile maneuverability. In the present study, a series of... 

Potential flow solvers simplify the mathematical formulation and achieve efficient solutions. The prediction of aerodynamic of dual rotor systems using computational fluid dynamic methods is difficult task due to the interference effects between the wakes shed from the rotors. In the present work, a free wake vortex lattice method is used to predict the vertical wake and blade loading of dual rotors in hover. In this approach the blades are modeled as flat plates with zero thickness and ring vortices are distributed on the surface of each blade. When the blades rotate, vortices are shed into the wake and freely move with a local velocity induced by the effects of the vortices on the blades... 

The numerical solution of the parabolized Navier-Stokes (PNS) equations for accurate computation of hypersonic axisymmetric flowfield with real gas effects is obtained by using the fourth-order compact finite-difference method. The PNS equations in the general curvilinear coordinates are solved by using the implicit finite-difference algorithm of Beam and Warming type with a high-order compact accuracy. A shock fitting procedure is utilized in the compact PNS scheme to obtain accurate solutions in the vicinity of the shock. To stabilize the numerical solution, numerical dissipation term and filters are used. The main advantage of the present formulation is that the basic flow variables... 

A collocation spectral scheme is used to compute inviscid supersonic high temperature flow in diverging nozzles. In order to include the real gas effects of high temperature air, chemical equilibrium is implemented via Tannehill et al. correlations. The proposed method is used for numerical simulation of internal inviscid flow of air in quasi one-dimensional and two-dimensional forms. Chebyshev spectral method is applied to the governing equations in order to discritize spatial differentiation terms and an explicit four stage Runge-Kutta method is chosen for time integration. The equations in primitive formulation are modified to include real gas effects. The effects of equilibrium or... 

Nowadays, one of the most important desires of the human being is to reduce his living environmental pollution. Using the diluted combustion systems in new gas turbines in order to produce the minimum amount of has been done to satisfy this desire. It should be noted that reducing this amount and using the lower flame temperature will result in some consequences. The most important problem occurred in industrial and aerial gas turbines are the instability of the combustion due to dilution of the fuel to air mixture which it results in heat release fluctuations. If the heat release fluctuations and acoustic pressure are in the same phases, the amplitude of the fluctuations will increase which... 

In the present study, the numerical solution of 2D inviscid compressible ideal magnetohydrodynamic (MHD) flows by using the spectral difference (SD) method on unstructured meshes is performed. The SD method combines the most desirable features of structured and unstructured grid methods to have computational efficiency and geometric flexibility to accurately compute flow over complex geometries. In the SD method, two sets of structured points, namely “unknown points” and “flux points”, are defined in each cell to support the reconstruction of given order of accuracy. The differential form of the conservation laws is satisfied at nodal unknown points while the flux derivatives expressed in... 

In the present study, a high-order accurate numerical solution of steady incompressible slip flow and heat transfer in 2D microchannels is presented. The numerical method used is an alternating direction implicit operator scheme which is efficiently implemented to solve the incompressible Navier-Stokes equations in the primitive variables formulation using the artificial compressibility method. To stabilize the numerical solution, numerical filters are used. The present methodology considers the solution of the Navier-Stokes equations with¬ employing different slip boundary condition¬¬ (Maxwell,¬ ¬¬Hyperbolic tangent function of Knudsen number¬ and Beskok slip models)¬ ¬¬on the wall to model... 

Cavitation can occur in many fluid systems such as pumps, nozzles, hydrofoils and submarine vehicles and therefore, numerical modeling of this phenomenon has a significant importance. In this study, the numerical simulation of the cavitating flows through the Euler/Navier-Stokes equations employing the interface capturing method associated with a barotropic state law is performed. The system of governing equations is discretized using a cell-centered finite-volume algorithm and the fluxes are evaluated using a central-difference scheme. To account for density jumps across the cavity interface, the numerical dissipation terms with suitable density and pressure sensors are used. Since... 

In this study, fluid-solid interaction (FSI) is simulated computationally by using a high-order accurate numerical method. The two-dimensional incompressible viscous flows are considered in the fluid domain. The primary problem with solutions of the incompressible Navier–Stokes equations is the difficulty of coupling changes in the velocity field with changes in the pressure field while satisfying the continuity equation. Herein, the artificial compressibility method is used to overcome this difficulty. Preconditioning is implemented to reduce the stiffness of the system of equations to increase the convergence rate of the solution. Using preconditioning, physical solutions even at low... 

In this study, the numerical simulation of natural and ventilated cavitating flows is performed. The algorithm employs the homogenous, multiphase Navier-Stokes equations with appropriate mass transfer terms.The base line differential equations system is comprised of the mixture volume, mixture momentum and constituent volume fraction equations. A three species differential formulation is considered for constituent volume fraction transport/generation of liquid, condensable vapor and non-condensable gas fields.The system of governing equations is discretized using a cell-centered finite volume Roe’s upwind scheme. Both laminar and turbulent cavitating flows are considered in this study. For... 

In this study, the simulation of two-dimensional supersonic flows through microchannels in slip flow regime is performed using a lattice Boltzmann model (LBM). Traditional LB models have been used to simulate incompressible fluid flows and there are not suitable for modeling compressible or thermo-fluid flows. Herein, a recently developed LB model, namely, the finite difference lattice Boltzmann method (FDLBM), is employed to simulate compressible flows with embedded shocks. In this model, one can select particle velocities independently from the lattice configuration, and therefore, a correct and numerically stable multispeed thermal model by adopting more isotropic particle velocities can... 

The study of evaporation of fuel droplet and determination of the rate of vaporization are important in designing combustion chambers. For achieving high performance of a combustor, the evaporation of fuel droplets takes place within a high pressure environment. At these conditions, the use of low-pressure models is not appropriate and many effects that are assumed negligible at low ambient pressures become very important. For example, the solubility of the ambient gas into the liquid phase is increased by increasing the ambient pressure. In addition, the ideal gas assumption is not valid for these conditions and one should use an appropriate equation of state (EOS) that can predict the... 

In the present work, the numerical simulation of 2D inviscid compressible flows by using the spectral difference (SD) method on quadrilateral meshes is performed. The SD method combines the most desirable features of structured and unstructured grid methods to attain computational efficiency and geometric flexibility. Similar to the discontinuous Galerkin (DG) and spectral volume (SV) methods, the SD scheme utilizes the concept of discontinuous and high-order local representations to achieve conservation and high accuracy. The SD method is based on the finite-difference formulation and thus its formulation is simpler than the DG and SV methods ... 

In this study, the numerical simulation of cavitating flows for cryogens fluids is performed. The algorithm employs the homogenous, multiphase Euler/Navier-Stokes equations with the interface capturing method. The thermodynamic and thermal effects substantially impact the cavitation dynamics of cryogenic fluids and therefore these effects should be considered by solving the energy equation in conjunction with the mass and momentum conservation, and updating the fluid physical properties. Here, two cavitation modeling strategies, namely, the barotropic cavitation model and the transport equation-based model are used. Both laminar and turbulent cavitating flows are studied in this work. For... 

In the study, the simulation of two-dimensional cavitating flows is performed by applying a high-order accurate numerical method to the preconditioned, homogenous, multiphase Navier-Stokes equations. The baseline differential equations system is comprised of the mixture volume, mixture momentum and constituent volume fraction equations. A coordinate transformation is applied and the resulting system of governing equations in curvilinear coordinates is discretized using a fourth-order compact finite-difference scheme. The high-order accurate numerical scheme employing the suitable linear and nonlinear filters to account for density jumps across the cavity interface is shown to yield an... 

Two-dimensional incompressible flow analysis is one the most important applied issues in engineering and applied science field. Numerical solution of governing equations of flow requires exact computational grid generation.In complex geometries, generation of the grid which is coincident to the body is very difficult and time consuming. Immersed boundary method is an appropriate superseded method of body conformal grid generation in flow field numerical solution. In this method a grid which is not coincidentto bodyis generated and flow field properties are modified on points adjacent to the boundary of the object (Ghost Cell Method) to satisfy boundary conditions.
The purpose of this... 

In this study, the solution of compressible flows is performed using finite volume lattice Boltzmann method (FVLBM). A model associated with 13 discrete velocity vectors and 2 energy levels is used and the Boltzmann transport equation is solved using a cell-centered finite volume on structured meshes. The values of distribution functions on each cell faceare determined by averaging from their values at the two control points located on the center of two neighboring cells. The fourth-order Runge-Kutta time-stepping scheme is applied to discretize temporal derivative term. The second- and fourth-order numerical dissipation termsareadded to the algorithm to stabilize the solution when solving... 

In the present study, the numerical simulation of the compressible inviscid flow around helicopter rotor is performed using the solution of the preconditioned Euler equations. Three preconditioners proposed by Eriksson, Choi and Merkel, and Turkel are implemented in two- and three-dimensional upwind Euler flow solvers on unstructured meshes. The mathematical formulations of these preconditioning schemes for different sets of primitive variables are drawn and their eigenvalues and eigenvectors are compared with each others. For this aim, these preconditioning schemes are expressed in a unified formulation. A cell-centered finite volume Roe's upwind method is used for the discretization of the... 

In this study, the numerical simulation of cavitating flows with compressibility effects is performed. The algorithm employs the multiphase Euler equations with homogeneous equilibrium model. The baseline differential equations system is similar to the one-phase system of equations and comprised of the mixture density, mixture momentums and mixture energy equations. Thephases considered for cavitating flows is liquid-vapor and liquid-gas fields. The system of governing equations is discretized using a cell-centered finite volume AUSM’s upwind scheme. The computations are presented for steady noncavitating/cavitating flows around 1D/2Dproblems for different conditions. A sensitivity study is... 

In the present study, the numerical simulation of incompressible laminar and turbulent flows using a high-order finite difference lattice Boltzmann method is presented. To handle curved geometries with non uniform grids, the incompressible form of lattice Boltzmann equation is transformed into the generalized curvilinear coordinates and the spatial derivatives of the resulting equation are discretized using the fifth-order WENO scheme. The advantage of using the WENO-LBM developed is that it needs less number of grid points and remains stable even at high Reynolds number flows. For the temporal term, the fourth-order explicit Rung-Kutta scheme is adopted for laminar flow calculations and...