Sharif Digital Repository

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

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

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 study, a high-order finite-difference lattice Boltzmann solver is applied for simulating steady and unsteady three-dimensional incompressible flows. To achieve an accurate and robust flow solver, the incompressible form of the lattice Boltzmann equation in the three-dimensional generalized curvilinear coordinates is discretized spatially based on the fifth-order weighted essentially non-oscillatory (WENO) finite-difference scheme. To ensure the stability and temporal accuracy of the flow solver, the fourth-order Runge-Kutta method is used for the time integration. To examine the accuracy and performance of the flow solver, different three-dimensional incompressible flow... 

In the present study, the numerical simulation of the panel flutter in compressible inviscid flow is performed by the compact finite difference method. For this purpose, the 2D compressible Euler equations written in the arbitrary Lagrange-Eulerian form are considered and the resulting system of equations in the generalized curvilinear coordinates is solved by the fourth-order compact finite-difference method. An appropriate nonlinear filter is applied for the shock capturing and for the solution to be stable. The governing equation for the panel is also numerically solved by using the fourth-order compact finite difference method. The time integration in the flow domain is made by the... 

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

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, the numerical simulation of 2D inviscid compressible cavitation flow is performed by using the compact finite-difference method. The problem formulation is based on the multiphase compressible Euler equations with the assumption of the homogeneous equilibrium model and the system of baseline differential equations is comprised of the continuity, momentum and energy equations for the vapor-liquid mixture. To complete the system of governing equations, the ideal gas relation is used for the vapor phase and the Tait relation is applied for the liquid phase, and therefore, the compressibility effects are considered for both the vapor and liquid phases. To analyze the flow... 

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

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

In this study, the Shokov-BGK model of the Boltzmann equation is reformulated and generalized to consider the real gas effects. At first, the formulation is performed to consider an arbitrary specific heats ratio and the correct Prandtl number for polyatomic gases. Here, the resulting equations of the present formulation are numerically solved by applying the high-order finite-difference weighted essentially non-oscillatory (WENO) scheme. The present solution method is tested by computing the one-dimension Reiman problem with different specific heats ratios for a wide range of the Knudsen numbers. The results are compared with the available gas-kinetic results which show good agreement. It... 

In the present study, the computation of the viscous compressible flow over two-dimensional geometries is performed by using the immersed boundary method and applying a second-order finite volume scheme. For the solution of the governing equations, a uniform Cartesian grid that is not coincident with the body surface is used and the boundary conditions on the wall are satisfied by the ghost-cell immersed boundary method. The spatial discretization of the fluid equations is carried out using the second-order central difference finite volume scheme and the time integration is performed by applying the fourth-order Runge-Kutta method. To stabilize the solution algorithm and reduce unwanted... 

In this study, 2-D compressible viscous and inviscid flows are simulated by using a finite volume Lattice Boltzmann method. Two different models, namely, the Qu model and Watari model are employed for compressible flows simulations. The first model includes 13 discrete velocity vectors and 2 energy levels in which the Maxwellian function is replaced with a simple function for describing the distribution function that is suitable for inviscid flow simulations. The second model is a thermal multi-velocity model with isotropic tensors up to seventh rank that is suitable for compressible viscous and inviscid flow simulations with arbitrary specific heats ratio. In both the models, lattice... 

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

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

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