Abstract: The purpose of this work is to simulate the flow at the exit of Vulcan 1 engine of European launcher Ariane 5. The geometry of the propellant nozzle is already determined using the characteristics method. The pressure in the outlet section of the nozzle is less than atmospheric pressure on the ground, causing the existence of oblique and normal shock waves at the exit. During the rise of the launcher, the atmospheric pressure decreases and the shock wave disappears. The code allows the capture of shock wave at exit of nozzle. The numerical technique uses the Flux Vector Splitting method of Van Leer to ensure convergence and avoid the calculation instabilities. The Courant, Friedrichs and Lewy coefficient (CFL) and mesh size level are selected to ensure the numerical convergence. The nonlinear partial derivative equations system which governs this flow is solved by an explicit unsteady numerical scheme by the finite volume method. The accuracy of the solution depends on the size of the mesh and also the step of time used in the discretized equations. We have chosen in this study the mesh that gives us a stationary solution with good accuracy.
Abstract: In this study, we have analyzed the transport of analytes
under a two dimensional steady incompressible flow of power-law
fluids through rectangular nanochannel. A mathematical model
based on the Cauchy momentum-Nernst-Planck-Poisson equations is
considered to study the combined effect of mixed electroosmotic
(EO) and pressure driven (PD) flow. The coupled governing
equations are solved numerically by finite volume method. We
have studied extensively the effect of key parameters, e.g., flow
behavior index, concentration of the electrolyte, surface potential,
imposed pressure gradient and imposed electric field strength on
the net average flow across the channel. In addition to study
the effect of mixed EOF and PD on the analyte distribution
across the channel, we consider a nonlinear model based on
general convective-diffusion-electromigration equation. We have also
presented the retention factor for various values of electrolyte
concentration and flow behavior index.
Abstract: Natural circulation loops (NCLs) are buoyancy driven flow systems without any moving components. NCLs have vast applications in geothermal, solar and nuclear power industry where reliability and safety are of foremost concern. Due to certain favorable thermophysical properties, especially near supercritical regions, carbon dioxide can be considered as an ideal loop fluid in many applications. In the present work, a high temperature NCL that uses supercritical carbon dioxide as loop fluid is analysed. The effects of relevant design and operating variables on loop performance are studied. The system operating under steady state is modelled taking into account the axial conduction through loop fluid and loop wall, and heat transfer with surroundings. The heat source is considered to be a heater with controlled heat flux and heat sink is modelled as an end heat exchanger with water as the external cold fluid. The governing equations for mass, momentum and energy conservation are normalized and are solved numerically using finite volume method. Results are obtained for a loop pressure of 90 bar with the power input varying from 0.5 kW to 6.0 kW. The numerical results are validated against the experimental results reported in the literature in terms of the modified Grashof number (Grm) and Reynolds number (Re). Based on the results, buoyancy and friction dominated regions are identified for a given loop. Parametric analysis has been done to show the effect of loop diameter, loop height, ambient temperature and insulation. The results show that for the high temperature loop, heat loss to surroundings affects the loop performance significantly. Hence this conjugate heat transfer between the loop and surroundings has to be considered in the analysis of high temperature NCLs.
Abstract: The objective of this research work is to investigate
for one dimensional transient radiative transfer equations with
conduction using finite volume method. Within the infrastructure of
finite-volume, we obtain the conservative discretization of the terms
in order to preserve the overall conservative property of finitevolume
schemes. Coupling of conductive and radiative equation
resulting in fluxes is governed by the magnitude of emissivity,
extinction coefficient, and temperature of the medium as well as
geometry of the problem.
The problem under consideration has been solved, for a slab
dominating radiation coupled with transient conduction based on
finite volume method. The boundary conditions are also chosen so as
to give a good model of the discretized form of radiation transfer
equation. The important feature of the present method is flexibility in
specifying the control angles in the FVM, while keeping the
simplicity in the solution procedure.
Effects of various model parameters are examined on the
distributions of temperature, radiative and conductive heat fluxes and
incident radiation energy etc. The finite volume method is considered
to effectively evaluate the propagation of radiation intensity through
a participating medium.
Abstract: Unsteady flow and heat transfer from a circular
cylinder in cross-flow is studied numerically. The governing
equations are solved by using finite volume method. Reynolds
number varies in range of 50 to 200; in this range flow is considered
to be laminar and unsteady. Al2O3 nanoparticle with volume fraction
in range of 5% to 20% is added to pure water. Effects of adding
nanoparticle to pure water on lift and drag coefficient and Nusselt
number is presented. Addition of Al2O3 has inconsiderable effect on
the value of drags and lift coefficient. However, it has significant
effect on heat transfer; results show that heat transfer of Al2O3
nanofluid is about 9% to 36% higher than pure water.
Abstract: This work deals with the problem of MHD mixed
convection in a completely porous and differentially heated vertical
channel. The model of Darcy-Brinkman-Forchheimer with the
Boussinesq approximation is adopted and the governing equations are
solved by the finite volume method. The effects of magnetic field and
buoyancy force intensities are given by the Hartmann and Richardson
numbers respectively, as well as the Joule heating represented by
Eckert number on the velocity and temperature fields, are examined.
The main results show an augmentation of heat transfer rate with the
decrease of Darcy number and the increase of Ri and Ha when Joule
heating is neglected.
Abstract: Radiative heat transfer in participating medium was
carried out using the finite volume method. The radiative transfer
equations are formulated for absorbing and anisotropically scattering
and emitting medium. The solution strategy is discussed and the
conditions for computational stability are conferred. The equations
have been solved for transient radiative medium and transient
radiation incorporated with transient conduction. Results have been
obtained for irradiation and corresponding heat fluxes for both the
cases. The solutions can be used to conclude incident energy and
surface heat flux. Transient solutions were obtained for a slab of heat
conducting in slab and by thermal radiation. The effect of heat
conduction during the transient phase is to partially equalize the
internal temperature distribution. The solution procedure provides
accurate temperature distributions in these regions. A finite volume
procedure with variable space and time increments is used to solve
the transient radiation equation. The medium in the enclosure
absorbs, emits, and anisotropically scatters radiative energy. The
incident radiations and the radiative heat fluxes are presented in
graphical forms. The phase function anisotropy plays a significant
role in the radiation heat transfer when the boundary condition is
non-symmetric.
Abstract: Heat transfer from flat tube is studied numerically. Reynolds number is defined base on equivalent circular tube which is varied in range of 100 to 300. In these range of Reynolds number flow is considered to be laminar, unsteady, and incompressible. Equations are solved by using finite volume method. Results show that increasing l/D from 1 to 2 has insignificant effect on heat transfer and Nusselt number of flat tube is slightly lower than circular tube. However, thermal-hydraulic performance of flat tube is up to 2.7 times greater than circular tube.
Abstract: This paper investigates the natural convection heat transfer performance in a complex-wavy-wall cavity filled with power-law fluid. In performing the simulations, the continuity, Cauchy momentum and energy equations are solved subject to the Boussinesq approximation using a finite volume method. The simulations focus specifically on the effects of the flow behavior index in the power-law model and the Rayleigh number on the flow streamlines, isothermal contours and mean Nusselt number within the cavity. The results show that pseudoplastic fluids have a better heat transfer performance than Newtonian or dilatant fluids. Moreover, it is shown that for Rayleigh numbers greater than Ra=103, the mean Nusselt number has a significantly increase as the flow behavior index is decreased.
Abstract: The present study is an analysis of the forced convection heat transfer in porous channel with an oriented jet at the inlet with uniform velocity and temperature distributions. The upper wall is insulated when the bottom one is kept at constant temperature higher than that of the fluid at the entrance. The dynamic field is analysed by the Brinkman-Forchheimer extended Darcy model and the thermal field is traduced by the energy one equation model. The numerical solution of the governing equations is obtained by using the finite volume method. The results mainly concern the effect of Reynolds number, jet angle and thermal conductivity ratio on the flow structure and local and average Nusselt numbers evolutions.
Abstract: Pressure loss in ductworks is an important factor to be considered in design of engineering systems such as power-plants, refineries, HVAC systems to reduce energy costs. Ductwork can be composed by straight ducts and different types of fittings (elbows, transitions, converging and diverging tees and wyes). Duct fittings are significant sources of pressure loss in fluid distribution systems. Fitting losses can be even more significant than equipment components such as coils, filters, and dampers. At the present work, a conventional 90o round elbow under turbulent incompressible airflow is studied. Mass, momentum, and k-e turbulence model equations are solved employing the finite volume method. The SIMPLE algorithm is used for the pressure-velocity coupling. In order to validate the numerical tool, the elbow pressure loss coefficient is determined using the same conditions to compare with ASHRAE database. Furthermore, the effect of Reynolds number variation on the elbow pressure loss coefficient is investigated. These results can be useful to perform better preliminary design of air distribution ductworks in air conditioning systems.
Abstract: One of the major parts of a jet engine is air intake,
which provides proper and required amount of air for the engine to
operate. There are several aerodynamic parameters which should be
considered in design, such as distortion, pressure recovery, etc. In
this research, the effects of lip ice accretion on pitot intake
performance are investigated. For ice accretion phenomenon, two
supervised multilayer neural networks (ANN) are designed, one for
ice shape prediction and another one for ice roughness estimation
based on experimental data. The Fourier coefficients of transformed
ice shape and parameters include velocity, liquid water content
(LWC), median volumetric diameter (MVD), spray time and
temperature are used in neural network training. Then, the subsonic
intake flow field is simulated numerically using 2D Navier-Stokes
equations and Finite Volume approach with Hybrid mesh includes
structured and unstructured meshes. The results are obtained in
different angles of attack and the variations of intake aerodynamic
parameters due to icing phenomenon are discussed. The results show
noticeable effects of ice accretion phenomenon on intake behavior.
Abstract: The present work is a numerical simulation of
nanofluids flow in a double pipe heat exchanger provided with
porous baffles. The hot nanofluid flows in the inner cylinder, whereas
the cold nanofluid circulates in the annular gap. The Darcy-
Brinkman-Forchheimer model is adopted to describe the flow in the
porous regions, and the governing equations with the appropriate
boundary conditions are solved by the finite volume method. The
results reveal that the addition of metallic nanoparticles enhances the
rate of heat transfer in comparison to conventional fluids but this
augmentation is accompanied by an increase in pressure drop. The
highest heat exchanger performances are obtained when
nanoparticles are added only to the cold fluid.
Abstract: In this paper 2D Simulation of catalytic Fixed Bed Reactor in Fischer-Tropsch Synthesis of GTL technology has been performed utilizing computational fluid dynamics (CFD). Synthesis gas (a mixture of carbon monoxide and hydrogen) has been used as feedstock. The reactor was modeled and the model equations were solved employing finite volume method. The model was validated against the experimental data reported in literature. The comparison showed a good agreement between simulation results and the experimental data. In addition, the model was applied to predict the concentration contours of the reactants and products along the length of reactor.
Abstract: We have studied the migration of a charged permeable aggregate in electrolyte under the influence of an axial electric field and pressure gradient. The migration of the positively charged aggregate leads to a deformation of the anionic cloud around it. The hydrodynamics of the aggregate is governed by the interaction of electroosmotic flow in and around the particle, hydrodynamic friction and electric force experienced by the aggregate. We have computed the non-linear Nernest-Planck equations coupled with the Dracy- Brinkman extended Navier-Stokes equations and Poisson equation for electric field through a finite volume method. The permeability of the aggregate enable the counterion penetration. The penetration of counterions depends on the volume charge density of the aggregate and ionic concentration of electrolytes at a fixed field strength. The retardation effect due to the double layer polarization increases the drag force compared to an uncharged aggregate. Increase in migration sped from the electrophretic velocity of the aggregate produces further asymmetry in charge cloud and reduces the electric body force exerted on the particle. The permeability of the particle have relatively little influence on the electric body force when Double layer is relatively thin. The impact of the key parameters of electrokinetics on the hydrodynamics of the aggregate is analyzed.
Abstract: In this work, axisymetric CFD simulation of fixed bed
GTL reactor has been conducted, using computational fluid dynamics
(CFD). In fixed bed CFD modeling, when N (tube-to-particle
diameter ratio) has a large value, it is common to consider the packed
bed as a porous media. Synthesis gas (a mixture of predominantly
carbon monoxide and hydrogen) was fed to the reactor. The reactor
length was 20 cm, divided to three sections. The porous zone was in
the middle section of the reactor. The model equations were solved
employing finite volume method. The effects of particle diameter,
bed voidage, fluid velocity and bed length on pressure drop have
been investigated. Simulation results showed these parameters could
have remarkable impacts on the reactor pressure drop.
Abstract: Solid fuel transient burning behavior under oxidizer
gas flow is numerically investigated. It is done using analysis of the
regression rate responses to the imposed sudden and oscillatory
variation at inflow properties. The conjugate problem is considered
by simultaneous solution of flow and solid phase governing
equations to compute the fuel regression rate. The advection
upstream splitting method is used as flow computational scheme in
finite volume method. The ignition phase is completely simulated to
obtain the exact initial condition for response analysis. The results
show that the transient burning effects which lead to the combustion
instabilities and intermittent extinctions could be observed in solid
fuels as the solid propellants.
Abstract: Double-diffusive steady convection in a partially
porous cavity with partially permeable walls and under the combined
buoyancy effects of thermal and mass diffusion was analysed
numerically using finite volume method.
The top wall is well insulated and impermeable while the bottom
surface is partially well insulated and impermeable and partially
submitted to constant temperature T1 and concentration C1. Constant
equal temperature T2 and concentration C2 are imposed along the
vertical surfaces of the enclosure. Mass suction/injection and
injection/suction are respectively considered at the bottom of the
porous centred partition and at one of the vertical walls.
Heat and mass transfer characteristics as streamlines and average
Nusselt numbers and Sherwood numbers were discussed for different
values of buoyancy ratio, Rayleigh number, and injection/suction
coefficient.
It is especially noted that increasing the injection factor
disadvantages the exchanges in the case of the injection while the
transfer is augmented in case of suction. On the other hand, a critical
value of the buoyancy ratio was highlighted for which heat and mass
transfers are minimized.
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: In this research the separation efficiency of deoiling hydrocyclone is evaluated using three-dimensional simulation of multiphase flow based on Eulerian-Eulerian finite volume method. The mixture approach of Reynolds Stress Model is also employed to capture the features of turbulent multiphase swirling flow. The obtained separation efficiency of Colman's design is compared with available experimental data and showed that the separation curve of deoiling hydrocyclones can be predicted using numerical simulation.