Sharif Digital Repository

A high-order compact finite-difference lattice Boltzmann method (CFDLBM) is proposed and applied to accurately compute steady and unsteady incompressible flows. Herein, the spatial derivatives in the lattice Boltzmann equation are discretized by using the fourth-order compact FD scheme, and the temporal term is discretized with the fourth-order Runge-Kutta scheme to provide an accurate and efficient incompressible flow solver. A high-order spectral-type low-pass compact filter is used to stabilize the numerical solution. An iterative initialization procedure is presented and applied to generate consistent initial conditions for the simulation of unsteady flows. A sensitivity study is also... 

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

In this work, we aim to study the effect of temperature gradients on the combined natural convection-radiation heat transfer problem in participating media. To impose this combined effect, we first solve the radiative transfer equation in an absorbing and emitting media. Then, we suitably add the radiation source terms to the energy equation and solve the fluid flow equations. Literature shows that many incompressible algorithms use the Boussinesq assumption to model the thermobuoyant force; however, the validity of this assumption is limited to cases with low temperature gradient distributions. Evidently, Boussinesq assumption would result in considerable errors in high temperature gradient... 

In this research, numerical optimization of the rear part of a gas turbine, consisting of a single stage axial turbine is carried out. Automated aerodynamic shape optimization is performed by coupling a CFD flow simulation code with the Genetic Algorithm. An effective multi-point optimization method to improve efficiency and/or pressure ratio of the axial turbine is performed. Some variations of optimization parameters such as lean and sweep angels of stator and rotor blades are accomplished. Furthermore, during the optimization process, three-dimensional and turbulent flow field is numerically investigated using a compressible Navier-Stokes solver. The gas turbine experimental... 

In this paper, the ball valve performance is numerically simulated using an unstructured CFD (Computational Fluid Dynamics) code based on the finite volume method. Navier-Stokes equations in addition to a transport equation for the vapor volume fraction were coupled in the RANS solver. Separation is modeled very well with a modification of turbulent viscosity. The results of CFD calculations of flow through a ball valve, based on the concept of experimental data, are described and analyzed. Comparison of the flow pattern at several opening angles is investigated. Pressure drop behind the ball valve and formation of the vortex flow downstream the valve section are also discussed. As the... 

This study describes a multidimensional 3D/lumped parameter (LP) model which contains appropriate inflow/outflow boundary conditions in order to model the entire human arterial trees. A new extensive LP model of the entire arterial network (48 arteries) was developed including the effect of vessel diameter tapering and the parameterization of resistance, conductor and inductor variables. A computer aided-design (CAD) algorithm was proposed to efficiently handle the coupling of two or more 3D models with the LP model, and substantially lessen the coupling processing time. Realistic boundary conditions and Navier-Stokes equations in healthy and stenosed models of carotid artery bifurcation... 

The present work is concerned with the steady incompressible flow through a parallel plate channel with stretching walls under an externally applied magnetic field. The governing continuity and Navier-Stokes equations are reduced to a fourth order nonlinear differential equation by using vorticity definition and similarity solution transformation. The obtained equations are solved by applying the analytical homotopy perturbation method (HPM). The method is called order of magnitude suggested for simplifying series solution to finite expression that is useful in engineering problems. The results are verified by comparing with numerical solutions and demonstrate a good accuracy of the obtained... 

Flow hydrodynamics of laminar falling film of aqueous Li Br solution (Li Br - H2O) on a horizontal elliptical tube has been investigated in this research. The film velocity distribution and film thickness, namely, the flow characteristics are determined by solving analytically simultaneous simplified Navier - Stokes equations and continuity equation in polar and Cartesian coordinates. The analysis is based on steady state laminar flow of falling liquid film of Li Br - H2O on a horizontal elliptical tube in polar model and Cartesian model (CM), for cases in which traction on the film surface is considered negligible. Models are compared with each other in three cases of aspect ratios (Ar),... 

Effects of sprue base size and design on flow pattern during aluminum gravity casting have been investigated by employing different sprue base sizes and using computational fluid dynamics (CFD). Calculations was carried out using SUTCAST simulation software based on solving Navier-Stokes equation and tracing the free surface using SOLA-VOF algorithm. Flow pattern was analyzed with focusing on streamlines and velocity distribution in sprue base, runner and in-gate. Increasing well size was produced a vortex flow at the bottom of sprue base which increased the surface velocity of liquid metal in runner. Using a rather big sprue well could eliminate vena contracta, but in-gate velocity was... 

In this study, the instability of spinning cylindrical shells partially filled with viscous liquid is investigated. Based on the Navier-Stokes equations for the incompressible flow, a 2D model is developed for liquid motion at each section of the cylinder. The governing equations of the cylinder vibrations are obtained based on the first-order shear deformable shell theory. The nonpenetration and no-slip boundary conditions of the flow on the wetted surface of the cylinder relate the liquid motion to the shell vibrations. Also the liquid pressure exerted on the cylinder wall combines the vibrations of the rotary cylinder to the liquid motion. By using the obtained coupled liquid-structure... 

An analytical approach is employed to investigate the transient and steady-state stresses in an isotropic, homogeneous half-space subjected to moving concentrated loads with subsonic speeds. Applying the Stokes-Helmholtz resolution to the Navier's equation of motion for the half-space results in a system of wavetype partial differential equations. Based on the new moving coordinate system, a modified system of partial differential equations is obtained. Applying a concurrent two-sided and one-sided Laplace transformation, this system is modified to a system of ordinary differential equations, the solutions of which are obtained with respect to boundary conditions. The transformed transient... 

Pressure transient analysis in vertically heterogeneous reservoirs is examined. The inclusion of a separate model for the free fluid flow in the wellbore is essential to allow for hydraulic communication and mixing of the fluid issuing from different reservoir layers. A two-dimensional model coupling Darcy flow in the reservoir with Navier-Stokes flow in the wellbore is developed and solved by the finite element technique. The coupled wellbore-reservoir flow model is used to analyze a layered reservoir with an abrupt change in permeability and a thick formation showing a gradual change in permeability with depth. Contrary to conventional reservoir models, this new model is able to capture... 

In this article, smoothed particle hydrodynamics method is applied in order to study the free surface flow generated by two-dimensional planing flat plate. For this purpose, a two-dimensional smoothed particle hydrodynamics code is developed and validated by the well-known dam breaking problem. Four trim angles and three different velocities are considered to perform a parametric study to examine their physical effects. The obtained results from smoothed particle hydrodynamics are compared against the corresponding Reynolds-averaged Navier Stokes solutions. It is observed that at lower velocities, there exists a good agreement between the smoothed particle hydrodynamics and Reynolds-averaged... 

This paper is concerned with a quasi-3D design method for the radial and axial diffusers of a centrifugal compressor on the meridional plane. The method integrates a novel inverse design algorithm, called ball-spine algorithm (BSA), and a quasi-3D analysis code. The Euler equation is solved on the meridional plane for a numerical domain, of which unknown boundaries (hub and shroud) are iteratively modified under the BSA until a prescribed pressure distribution is reached. In BSA, unknown walls are composed of a set of virtual balls that move freely along specified directions called spines. The difference between target and current pressure distributions causes the flexible boundary to deform... 

A computational fluid dynamics code was developed to compute the flow inside and around a supersonic external compression axisymmetric intake. The code solves the Reynolds-averaged Navier-Stokes equations using an explicit finite volume method in a structured grid and uses the Baldwin-Lomax algebraic model to compute the turbulent viscosity coefficient. Experiments were performed to validate the predicted results and good agreements are achieved. In the next part of the research, a parametric study was undertaken using the designed base case at a constant Mach number of 2 and at 0° angle of attack. The effects of various important parameters such as free stream Mach number, spike deflection... 

The hydraulic efficiency of sedimentation basins is reduced by short-circuiting, circulation zones and bottom particleladen jets. Baffles are used to improve the sediment tank performance. In this study, laboratory experiments were used to examine the hydrodynamics of several baffle configurations. An accompanying numerical analysis was performed based on the 2-D Reynolds-averaged Navier-Stokes equations along with the k-ε turbulence closure model. The numerical model was supplemented with the volume-of-fluid technique, and the advection-diffusion equation to simulate the dynamics of particle-laden flow. Model predictions compared well with the experimental data. An empirical function was... 

In this article, the effects of volute cross section shape and centroid profile of a radial ow compressor volute were investigated. The performance characteristics of a turbocharger compressor were obtained experimentally by measuring rotor speed and ow parameters at the inlet and outlet of the compressor. The three-dimensional ow field model of the compressor was obtained numerically solving Navier-Stokes equations with SST turbulence model. The compressor characteristic curves were plotted. For model verification, the results were compared with experimental data, showing good agreement. Modification of a volute was performed by introducing a shape factor for volute cross section geometry.... 

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

In this article, the effects of volute cross section shape and centroid profile of a centrifugal compressor volute were investigated. The performance characteristics of a turbocharger compressor were obtained experimentally by measuring rotor speed and flow parameters at the inlet and outlet of the centrifugal compressor. The three dimensional flow field model of the compressor was obtained numerically solving Navier- Stokes equations with SST turbulence model. The compressor characteristic curves were plotted. For model verification, the results were compared with experimental data, showing good agreement. Modification of a volute was performed by introducing a shape factor for volute cross... 

In this study, a moving mesh finite-volume-based finite-element (FVBFE) method is suitably extended to simulate the effect of supercooled liquid water droplet content on ice formation and growth on wing sections. The method benefits from the advantages of both finite-volume and finiteelement methods, which promote achieving a more accurate solution and a higher efficient procedure in ice accretion calculations. The method solves the time-dependent Navier-Stokes (NS) equations on unstructured hybrid grid distributions. In this method, the convection terms are approximated at the cell faces using a physical influence upwinding scheme. We also use linear spring approach to move the hybrid mesh....