Abstract: The efforts to understand the heat transfer behavior of supercritical water in supercritical water cooled reactor (SCWR) are ongoing worldwide to fulfill the future energy demand. The higher thermal efficiency of these reactors compared to a conventional nuclear reactor is one of the driving forces for attracting the attention of nuclear scientists. In this work, a solution procedure has been described for solving supercritical fluid flow problems in complex geometries. The solution procedure is based on non-staggered grid. All governing equations are discretized by finite volume method (FVM) in curvilinear coordinate system. Convective terms are discretized by first-order upwind scheme and central difference approximation has been used to discretize the diffusive parts. k-ε turbulence model with standard wall function has been employed. SIMPLE solution procedure has been implemented for the curvilinear coordinate system. Based on this solution method, 3-D Computational Fluid Dynamics (CFD) code has been developed. In order to demonstrate the capability of this CFD code in supercritical fluid flows, heat transfer to supercritical water in circular tubes has been considered as a test problem. Results obtained by code have been compared with experimental results reported in literature.
Abstract: An essential component of a finite volume method (FVM) is the advection scheme that estimates values on the cell faces based on the calculated values on the nodes or cell centers. The most widely used advection schemes are upwind schemes. These schemes have been developed in FVM on different kinds of structured and unstructured grids. In this research, the physical influence scheme (PIS) is developed for a cell-centered FVM that uses an implicit coupled solver. Results are compared with the exponential differencing scheme (EDS) and the skew upwind differencing scheme (SUDS). Accuracy of these schemes is evaluated for a lid-driven cavity flow at Re = 1000, 3200, and 5000 and a backward-facing step flow at Re = 800. Simulations show considerable differences between the results of EDS scheme with benchmarks, especially for the lid-driven cavity flow at high Reynolds numbers. These differences occur due to false diffusion. Comparing SUDS and PIS schemes shows relatively close results for the backward-facing step flow and different results in lid-driven cavity flow. The poor results of SUDS in the lid-driven cavity flow can be related to its lack of sensitivity to the pressure difference between cell face and upwind points, which is critical for the prediction of such vortex dominant flows.
Abstract: Steady three-dimensional and two free surface waves
generated by moving bodies are presented, the flow problem to be
simulated is rich in complexity and poses many modeling challenges
because of the existence of breaking waves around the ship hull, and
because of the interaction of the two-phase flow with the turbulent
boundary layer. The results of several simulations are reported. The
first study was performed for NACA0012 of hydrofoil with different
meshes, this section is analyzed at h/c= 1, 0345 for 2D. In the second
simulation a mathematically defined Wigley hull form is used to
investigate the application of a commercial CFD code in prediction of
the total resistance and its components from tangential and normal
forces on the hull wetted surface. The computed resistance and wave
profiles are used to estimate the coefficient of the total resistance for
Wigley hull advancing in calm water under steady conditions. The
commercial CFD software FLUENT version 12 is used for the
computations in the present study. The calculated grid is established
using the code computer GAMBIT 2.3.26. The shear stress k-ωSST
model is used for turbulence modeling and the volume of fluid
technique is employed to simulate the free-surface motion. The
second order upwind scheme is used for discretizing the convection
terms in the momentum transport equations, the Modified HRIC
scheme for VOF discretization. The results obtained compare well
with the experimental data.
Abstract: Steady three-dimensional and two free surface waves
generated by moving bodies are presented, the flow problem to be
simulated is rich in complexity and poses many modeling challenges
because of the existence of breaking waves around the ship hull, and
because of the interaction of the two-phase flow with the turbulent
boundary layer. The results of several simulations are reported. The
first study was performed for NACA0012 of hydrofoil with different
meshes, this section is analyzed at h/c= 1, 0345 for 2D. In the second
simulation a mathematically defined Wigley hull form is used to
investigate the application of a commercial CFD code in prediction of
the total resistance and its components from tangential and normal
forces on the hull wetted surface. The computed resistance and wave
profiles are used to estimate the coefficient of the total resistance for
Wigley hull advancing in calm water under steady conditions. The
commercial CFD software FLUENT version 12 is used for the
computations in the present study. The calculated grid is established
using the code computer GAMBIT 2.3.26. The shear stress k-ωSST
model is used for turbulence modeling and the volume of fluid
technique is employed to simulate the free-surface motion. The
second order upwind scheme is used for discretizing the convection
terms in the momentum transport equations, the Modified HRIC
scheme for VOF discretization. The results obtained compare well
with the experimental data.
Abstract: The performance of Advection Upstream Splitting
Method AUSM schemes are evaluated against experimental flow
fields at different Mach numbers and results are compared with
experimental data of subsonic, supersonic and hypersonic flow fields.
The turbulent model used here is SST model by Menter. The
numerical predictions include lift coefficient, drag coefficient and
pitching moment coefficient at different mach numbers and angle of
attacks. This work describes a computational study undertaken to
compute the Aerodynamic characteristics of different air vehicles
configurations using a structured Navier-Stokes computational
technique. The CFD code bases on the idea of upwind scheme for the
convective (convective-moving) fluxes. CFD results for GLC305
airfoil and cone cylinder tail fined missile calculated on above
mentioned turbulence model are compared with the available data.
Wide ranges of Mach number from subsonic to hypersonic speeds are
simulated and results are compared. When the computation is done
by using viscous turbulence model the above mentioned coefficients
have a very good agreement with the experimental values. AUSM
scheme is very efficient in the regions of very high pressure gradients
like shock waves and discontinuities. The AUSM versions simulate
the all types of flows from lower subsonic to hypersonic flow without
oscillations.
Abstract: The spectral action balance equation is an equation that
used to simulate short-crested wind-generated waves in shallow water
areas such as coastal regions and inland waters. This equation consists
of two spatial dimensions, wave direction, and wave frequency which
can be solved by finite difference method. When this equation with
dominating convection term are discretized using central differences,
stability problems occur when the grid spacing is chosen too coarse.
In this paper, we introduce the splitting upwind schemes for avoiding
stability problems and prove that it is consistent to the upwind scheme
with same accuracy. The splitting upwind schemes was adopted
to split the wave spectral action balance equation into four onedimensional
problems, which for each small problem obtains the
independently tridiagonal linear systems. For each smaller system
can be solved by direct or iterative methods at the same time which
is very fast when performed by a multi-processor computer.
Abstract: In this study, a 3D combustion chamber was simulated
using FLUENT 6.32. Aims to obtain accurate information about the
profile of the combustion in the furnace and also check the effect of
oxygen enrichment on the combustion process. Oxygen enrichment is
an effective way to reduce combustion pollutant. The flow rate of air
to fuel ratio is varied as 1.3, 3.2 and 5.1 and the oxygen enriched
flow rates are 28, 54 and 68 lit/min. Combustion simulations
typically involve the solution of the turbulent flows with heat
transfer, species transport and chemical reactions. It is common to
use the Reynolds-averaged form of the governing equation in
conjunction with a suitable turbulence model. The 3D Reynolds
Averaged Navier Stokes (RANS) equations with standard k-ε
turbulence model are solved together by Fluent 6.3 software. First
order upwind scheme is used to model governing equations and the
SIMPLE algorithm is used as pressure velocity coupling. Species
mass fractions at the wall are assumed to have zero normal
gradients.Results show that minimum mole fraction of CO2 happens
when the flow rate ratio of air to fuel is 5.1. Additionally, in a fixed
oxygen enrichment condition, increasing the air to fuel ratio will
increase the temperature peak. As a result, oxygen-enrichment can
reduce the CO2 emission at this kind of furnace in high air to fuel
rates.
Abstract: A new numerical method for solving the twodimensional,
steady, incompressible, viscous flow equations on a
Curvilinear staggered grid is presented in this paper. The proposed
methodology is finite difference based, but essentially takes
advantage of the best features of two well-established numerical
formulations, the finite difference and finite volume methods. Some
weaknesses of the finite difference approach are removed by
exploiting the strengths of the finite volume method. In particular,
the issue of velocity-pressure coupling is dealt with in the proposed
finite difference formulation by developing a pressure correction
equation in a manner similar to the SIMPLE approach commonly
used in finite volume formulations. However, since this is purely a
finite difference formulation, numerical approximation of fluxes is
not required. Results obtained from the present method are based on
the first-order upwind scheme for the convective terms, but the
methodology can easily be modified to accommodate higher order
differencing schemes.
Abstract: A numerical study on the influence of electroosmotic flow on analyte preconcentration by isotachophoresis ( ITP) is made. We consider that the double layer induced electroosmotic flow ( EOF) counterbalance the electrophoretic velocity and a stationary ITP stacked zones results. We solve the Navier-Stokes equations coupled with the Nernst-Planck equations to determine the local convective velocity and the preconcentration dynamics of ions. Our numerical algorithm is based on a finite volume method along with a secondorder upwind scheme. The present numerical algorithm can capture the the sharp boundaries of step-changes ( plateau mode) or zones of steep gradients ( peak mode) accurately. The convection of ions due to EOF reduces the resolution of the ITP transition zones and produces a dispersion in analyte zones. The role of the electrokinetic parameters which induces dispersion is analyzed. A one-dimensional model for the area-averaged concentrations based on the Taylor-Aristype effective diffusivity is found to be in good agreement with the computed solutions.
Abstract: In this study, a 3D combustion chamber was simulated
using FLUENT 6.32. Aim to obtain detailed information on
combustion characteristics and _ nitrogen oxides in the furnace and
the effect of oxygen enrichment in a combustion process. Oxygenenriched
combustion is an effective way to reduce emissions. This
paper analyzes NO emission, including thermal NO and prompt NO.
Flow rate ratio of air to fuel is varied as 1.3, 3.2 and 5.1 and the
oxygen enriched flow rates are 28, 54 and 68 lit/min. The 3D
Reynolds Averaged Navier Stokes (RANS) equations with standard
k-ε turbulence model are solved together by Fluent 6.32 software.
First order upwind scheme is used to model governing equations and
the SIMPLE algorithm is used as pressure velocity coupling. Results
show that for AF=1.3, increase the oxygen flow rate of oxygen
reduction in NO emissions is Lance. Moreover, in a fixed oxygen
enrichment condition, increasing the air to fuel ratio will increase the
temperature peak, but not the NO emission rate. As a result, oxygen
enrichment can reduce the NO emission at this kind of furnace in low
air to fuel rates.
Abstract: In the present work flow past circular cylinder and
cylinder with rectangular and triangular wake splitter is studied to
improve aerodynamic parameters. The Comparison of drag
coefficient is tabulated for bare cylinder, cylinder with rectangular
and triangular wake splitters. Flow past circular cylinder and cylinder
with triangular and rectangular wake splitter is performed at
Reynoldsnumber 5, 20, 40, 50,80, 100.An incompressible PISO finite
volume code employing a non-staggered grid arrangement is used, a
second order upwind scheme is used for convective terms. The time
discretization is implicit and a Second order Crank-Nicholson scheme
is employed. Length of wake splitter in both configurations is taken
to be equal to diameter of cylinder. Wake length is found to be less
with rectangular wake splitter when compared to bare cylinder and
cylinder with triangular wake splitter. Coefficient of drag is found to
be less for triangular wake splitter when compared to bare cylinder &
cylinder with rectangular wake splitter.