Sharif Digital Repository

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

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, a Chebyshev collocation spectral lattice Boltzmann method (CCSLBM) is developed and assessed for the computation of low speed flows. Both steady and unsteady flows are considered here. The discrete Boltzmann equation (DBE) with the Bhatnagar-Gross-Krook (BGK) approximation based on the pressure distribution function is considered and the space discretization is performed by the Chebyshev collocation spectral method to achieve a highly accurate flow solver. To provide accurate unsteady solutions, the time integration of the temporal term in the LB equation is made by the fourth-order Runge-Kuta scheme. To achieve numerical stability and accuracy, the physical boundary... 

In this study,thecut-cells method is developed for simulating two-dimensional, inviscid, compressible flows with immersed boundaries. A finite volume method based on the second-order accurate central-difference scheme and the Runge-Kutta time stepping scheme is used. The key aspects that require to be considered in thedevelopment of such a solver are implementation of boundary conditions on the immersed boundaries and correct discretizing the governing equations in those cells cut by the boundaries. An appropriate interpolation procedure is applied to preserve the second-order spatial accuracy of the solver. The solution procedure is validated vs. well documented test problems for a wide... 

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

In the present study, the simulation of incompressible Electro-Magneto-hydrodynamic flows is performed using a finite volume lattice Boltzmann method (FVLBM). The Boltzmann transport equation is solved using a cell-centered finite volume method on structured meshes. A central difference scheme is used to discretize the spatial derivatives and the fourth-order numerical dissipation term is added to stabilize the solution. To discretize the temporal derivative, the fourth-order Runge-Kutta time stepping scheme is applied. The standard collision-streaming lattice Boltzmann method has been used to simulate EMHD flows in the literature, however, it has several deficiencies such as the... 

In this study, the numerical simulation of the two-dimensional incompressible flow around an oscillating airfoil is performed. For this aim, the incompressible Navier-Stokes equations based on the artificial compressibility approach written in the arbitrary Lagrangian-Eulerian form are considered. Then, the vorticity confinement method is incorporated in the formulation and the resulting system of equations is solved by a second-order central-difference finite volume method with the controllable dissipation terms. For the time integration, the implicit dual-time stepping scheme is implemented. At first, the numerical solution of the incompressible flow over the oscillating NACA0012 airfoil... 

In the present work, the numerical solution of 2D inviscid compressible Magneto-hydrodynamic flow is performed by using the spectral difference (SD) method on quadrilateral grids. In this numerical method, similar to the discontinuous Galerkin (DG) and spectral volume (SV) methods, the concept of the discontinuous and high-order local representations is used to achieve conservation property and high-order accuracy. In the SD method, the test function or the surface integral is not involved and thus it has a simpler formulation than the DG and SV methods. In this numerical method, two sets of structured points, namely unknown points and flux points, are defined in each cell to support the... 

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, two finite-difference lattice Boltzmann methods (FDLBM) are applied and assessed for the simulation of two-phase liquid-vapor flows with high density ratios. For this aim, the He-Shan-Doolen type lattice Boltzmann multiphase model is used and the spatial derivatives in the resulting system of equations are discretized by using the second-order central difference and modified Lax-Wendroff schemes. Suitable numerical dissipation terms and filters are applied to regularize the numerical solution and remove spurious waves generated by flow nonlinearities in smooth regions and at the same time to remove the numerical oscillations in the interface region of the two phases.Three... 

In the present study, an incompressible smoothed particle hydrodynamics method based on vorticity-stream function (VSF-SPH) formulation is developed and assessed for simulating steady and unsteady incompressible flows. The vorticity-stream function formulation in the Eulerian reference frame is written in a Lagrangian reference frame to provide an appropriate incompressible SPH algorithm. The advantage of the proposed smoothed particle hydrodynamics method based on the vorticity-stream function (VSF-SPH) formulation over the weakly compressible SPH (WCSPH) is that the VSF-SPH method is a truly incompressible SPH algorithm and it does not involve any approximate enforcement of the... 

In the present study, an incompressible smoothed particle hydrodynamics based on the artificial compressibility method is applied for simulating the incompressible turbulent flows. The Reynolds-averaged incompressible Navier–Stokes equations using the artificial compressibility method in the Eulerian reference frame are written in the Lagrangian reference frame to provide an appropriate incompressible SPH algorithm for the turbulent flow computations. Here, the k-L_m turbulence model, which is a simplified k-ϵ turbulence model, is used and formulated in the Lagrangian reference frame. The SPH formulation implemented here is based on an implicit dual-time stepping scheme to be capable of... 

In the present work, the deformation and collapse of a single cavitating bubble near solid boundaries is simulated by solving the preconditioned, homogenous, multiphaseNavier-Stokes equations. Up to now, all studies in the literature performed by the volume of fluid (VOF)approach to capture the bubble surface have been based on the pressure-based category in which the flow variables are calculated through solving the Poisson equation. Here, the density-based category is applied and the solution methodology is based on the artificial compressibility approach. The compressible form of the Navier-Stokes equations is applied inside the bubble and the liquid phase is assumed to be incompressible.... 

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

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 the current study, the numerical simulation of the turbulent cavitating flows is performed by solving the preconditioned, homogenous, multiphase Navier-Stokes equations. For the turbulence modeling, the standard two-equation k-ϵ model is used. The baseline differential equations system is comprised of the mixture volume, mixture momentum and constituent volume fraction equations together with two equations for the turbulence kinetic energy k and the turbulence energy dissipation rate ϵ. For the calculation of the eddy viscosity near the wall boundary, appropriate turbulence damping functions are applied to modify the source terms of the ϵ equation. The system of governing equations is...