Sharif Digital Repository

Multiple square cross section jets into a cross flow at three different velocity ratios, namely 0.5, 1.0 and 1.5, have been computationally simulated, using the Large Eddy Simulation (LES) approach. The finite volume method is applied in the computational methodologies, using an unsteady SIMPLE algorithm and employing a non-uniform staggered grid. All spatial and temporal terms in the Navier-Stokes equations have been discretized using the Power-Law and Crank-Nicolson schemes, respectively. Mean velocity profiles at different X-locations are compared with the existing experimental and Reynolds Averaged Navier-Stokes (RANS) computational results. Although the RANS computations require much... 

The numerical solution of the parabolized Navier-Stokes (PNS) and globally iterated PNS (IPNS) equations for accurate computation of hypersonic axisymmetric flowfields is obtained by using the fourth-order compact finite-difference method. The PNS and IPNS 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 both the compact PNS and IPNS schemes to obtain accurate solutions in the vicinity of the shock. The main advantage of the present formulation is that the basic flow variables and their first and second derivatives are... 

The three-dimensional supersonic turbulent flows over wrap-around fin missiles have been computed using the Thin Layer Navier-Stokes (TLNS) equations to reduce the computational efforts compared to those of the Full Navier-Stokes (FNS) equations. In this research, the missile configuration is divided into multi regions to enable fluid flow simulation using Personal Computers (PC). It also makes it possible to use a different number of nodes and distribution of grids in each region to enhance the accuracy. The Thin Layer Navier-Stokes equations in the generalized coordinate system were solved using an efficient, implicit, finite-difference factored algorithm of the Beam and Warming. For the... 

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

We extend a hybrid DSMC/Navier-Stokes (NS) approach to unify the DSMC and the NS simulators in one framework capable of solving the mixed non-equilibrium and near-equilibrium flow regions efficiently. Furthermore, we use a one-way state-based coupling (Dirichlet-Dirichlet boundary-condition coupling) to transfer the required information from the continuum region to the rarefied one. The current hybrid DSMC-NS frame is applied to the hypersonic flows over nanoflat plate and microcylinder cases. The achieved solutions are compared with the pure DSMC and NS solutions. The results show that the current hybrid approach predicts the surface heat transfer rate and shear stress magnitudes very... 

The main purpose of this work is to improve the efficiency and performance of a primitive finite volume element method which provides superior capability on a single platform. This method is suitably extended in order to use the advantages of parallel computing on multiprocessors or multicomputers. The method is fully implicit which renders huge sparse linear algebraic kernels. Nevertheless, the attention is focused on solving the sparse system rather than constructing it. The current method is a cell-centered scheme. Since the grid is unstructured, each non-boundary node engages with nodes on three or more surrounding elements around that node. Depending on global node numbering, the... 

A numerical upwind-biased procedure which respects the essence of upwinding is suitably extended in order to reduce the false diffusion induced by a first-order approximation. In this regard, some arbitrarily first and second order gradient terms are added to the primary upwind approximation. The additional terms are then discretized using second-order schemes which essentially produce dispersive errors. The suitable choices for the weights of the new added terms result in lowering the dissipative role of the original upwind scheme. Additionally, the implicit appearance of the third-order terms, which are the consequences of second-order discretizations, helps to reduce the dissipative... 

Analytical analysis of fluid flow in cylindrical microchannels subjected to uniform wall injection at various Reynolds numbers is presented. The classical Navier-Stokes equations are used in the present study. Mathematically, using an appropriate change of variable, Navier-Stokes equations are transformed to a set of nonlinear ordinary differential equations. The governing equations are solved analytically using series solution method. The presented analytical results can be used for the prediction of velocity profiles and pressure drops in the cylindrical micro channels. The results are validated against available data in the literature and have shown good agreement. Copyright © 2004 by... 

In this study, we apply a hybrid direct simulation Monte Carlo (DSMC)/Navier-Stokes (NS) frame to simulate the effects of catalyst or splitter plates in propulsive efficiency of micro/nanopropulsion systems. Our hybrid frame uses the local Knudsen number based on the gradient of the flow properties (KnGLL) to distinct the continuum and molecular regions. This frame also uses the state-based coupling (Dirichlet-Dirichlet boundary-condition coupling) to transfer the information between the two regions. We simulate typical micro/nanopropulsion systems consisting of channels, catalyst or splitter plates, and convergent-divergent nozzles. According to the Kn GLL, we apply the NS solver to the... 

Green iron oxide pellets are indurated with thermal treatment in pelletizing plant to achieve sufficient mechanical properties and to be used in iron-steel industry. In the first stage, the pellets are dried with hot air coming from the firing unit. High heat transfer makes the evaporation rate exceeding the outgoing steam rate from the pellets and as a result, increasing the inner pressure of the pellets and fragment. On the other hand, low heat transfer causes low production rate. Using the related equations and considering the pellets' moisture and porosity and applying the pellets' condition on parameters such as specific heat capacity and heat conductivity, an acceptable standard is... 

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

The S-shaped air intakes are very common shapes due to their ease in the engine-body integration or Radar Cross Section, RCS, specifications especially for fighter aircrafts. The numerical shape optimization of an S-shaped air intake using adjoint method is conducted. The flow of a specified air intake that uses S-duct M2129 is simulated using three dimensional (3D) numerical solution of Reynolds-Averaged Navier-Stokes equation along with k-ω SST turbulence model. The main purpose of this optimization scheme is to maximize the total pressure recovery (TPR). Further, the scheme is developed in such a way that would be applicable in industry thru satisfying specified constraint requirements.... 

Moderate or intense low-oxygen dilution (MILD) combustion of liquid fuels has attracted attention because of its advantages in industrial burners and gas turbine applications. Here, numerical investigation has been conducted on an experimental MILD turbulent spray burner. The H∥ flame of Delft spray in a hot co-flow burner is selected, and the Reynolds averaged Navier-Stokes/eddy dissipation concept framework with 40 species/180 reversible reactions through a skeletal chemical mechanism is used in addition to unsteady Lagrangian tracking of spray droplets to investigate the flame structure and chemical kinetic of reacting flow field. At first, current numerical results were compared with... 

Linear stability analysis is used to characterize the dynamics of a plane jet by incorporating non-modal stability analysis besides classical global temporal stability analysis. It is explained that similar shapes of different global modes are the result of non-normal characteristics of linearized Navier Stokes equations. Optimal initial disturbances and their eigenfunctions together with transient energy growth are obtained for different time horizons and Reynolds numbers of the jet in the linear unstable configuration. These structures are localized at the upstream of the jet nozzle at the boundary layer. The transient growth of the inlet perturbation in limited time bounds is found in the... 

Moderate or Intense Low-oxygen Dilution (MILD) combustion is a technology with important characteristics such as significant low emission and high-efficiency combustion. The hydrogen enrichment of conventional fuels is also of interest due to its favorable characteristics, such as low carbon-containing pollutants, high reaction intensity, high flammability, and thus fuel usage flexibility. In this study, the effects of adding hydrogen to methane and syngas fuels have been investigated under conditions of MILD combustion through numerical simulation of a well-set-up MILD burner. The Reynolds-Averaged Navier-Stokes (RANS) approach is adopted along the Eddy Dissipation Concept (EDC) combustion... 

Because reversible fans must be able to provide similar performance in either rotation direction, they need to have completely symmetrical blades. In the present study, an elliptic profile with itsmaximum thickness of 8%chord length at themiddle of the chord line is studied. For verification, theNACA0012 airfoilwas once analyzed two-dimensionally in the isolatedformationandonce in cascade formation with0.55 and 0.83 solidities, andthe obtainedresultswere comparedwith numerical and experimental results. Three well-known Reynolds-averaged Navier-Stokes turbulence models including Spalart-Allmaras, realizable κ-ϵ, and shear-stress transport κ-ω were used to find the suitable turbulence method.... 

Because reversible fans must be able to provide similar performance in either rotation direction, they need to have completely symmetrical blades. In the present study, an elliptic profile with itsmaximum thickness of 8%chord length at themiddle of the chord line is studied. For verification, theNACA0012 airfoilwas once analyzed two-dimensionally in the isolatedformationandonce in cascade formation with0.55 and 0.83 solidities, andthe obtainedresultswere comparedwith numerical and experimental results. Three well-known Reynolds-averaged Navier-Stokes turbulence models including Spalart-Allmaras, realizable κ-ϵ, and shear-stress transport κ-ω were used to find the suitable turbulence method.... 

The accuracy of Reynolds averaged Navier-Stokes (RANS) turbulence models to predict the behavior of two-dimensional (2-D) density current has been examined. In this work, a steady density current is simulated by the κ -ε , κ -ε RNG, two-layer κ -ε and modified v 2̄ - f models. All models are compared with available experimental data. Density current with uniform velocity and concentration enters a channel via a sluice gate into a lighter ambient fluid and moves forward down-slope. The eddy-viscosity concept cannot accurately simulate this flow because of two stress production structures in it. Results show that all isotropic models have a weak outcome for this current, but with improving the... 

This study seeks to explore the ability of unsteady Reynolds-averaged Navier-Stokes (URANS) simulation approach for resolving two-dimensional (2D) gravity currents on fine computational meshes. A 2D URANS equations closed by a buoyancy-modified k-ε turbulence model are integrated in time with a second-order fractional step approach coupled with a direct implicit method and discretized in space on a staggered mesh using a second-order accurate finite volume approach incorporating a high resolution semi-Lagrangian technique for the convective terms. A series of 2D simulations are carried out for gravity currents from both discontinuous and continuous sources. Comparisons with experimental...