Abstract: In designing a low-energy-consuming buildings, the heat transfer through a large glass or wall becomes critical. Multiple layers of the window glasses and walls are employed for the high insulation. The gravity driven air flow between window glasses or wall layers is a natural heat convection phenomenon being a key of the heat transfer. For the first step of the natural heat transfer analysis, in this study the development and application of a finite volume method for the numerical computation of viscous incompressible flows is presented. It will become a part of the natural convection analysis with high-order scheme, multi-grid method, and dual-time step in the future. A finite volume method based on a fully-implicit second-order is used to discretize and solve the fluid flow on unstructured grids composed of arbitrary-shaped cells. The integrations of the governing equation are discretised in the finite volume manner using a collocated arrangement of variables. The convergence of the SIMPLE segregated algorithm for the solution of the coupled nonlinear algebraic equations is accelerated by using a sparse matrix solver such as BiCGSTAB. The method used in the present study is verified by applying it to some flows for which either the numerical solution is known or the solution can be obtained using another numerical technique available in the other researches. The accuracy of the method is assessed through the grid refinement.
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: At present, vibrations of rotors of gas transmittal unit evade sustainable forecasting. This paper describes elastic oscillation modes in resilient supports and rotor impellers modeled during computational experiments with regard to interference in the system of gas-dynamic flow and compressor rotor. Verification of aeroelastic approach was done on model problem of interaction between supersonic jet in shock tube with deformed plate. ANSYS 15.0 engineering analysis system was used as a modeling tool of numerical simulation in this paper. Finite volume method for gas dynamics and finite elements method for assessment of the strain stress state (SSS) components were used as research methods. Rotation speed and material’s elasticity modulus varied during calculations, and SSS components and gas-dynamic parameters in the dynamic system of gas-dynamic flow and compressor rotor were evaluated. The analysis of time dependence demonstrated that gas-dynamic parameters near the rotor blades oscillate at 200 Hz, and SSS parameters at the upper blade edge oscillate four times higher, i.e. with blade frequency. It has been detected that vibration amplitudes correction in the test points at magnetic bearings by aeroelasticity may correspond up to 50%, and about -π/4 for phases.
Abstract: This paper discusses the applicability of the numerical model for a damage prediction method of the accidental hydrogen explosion occurring in a hydrogen facility. The numerical model was based on an unstructured finite volume method (FVM) code “NuFD/FrontFlowRed”. For simulating unsteady turbulent combustion of leaked hydrogen gas, a combination of Large Eddy Simulation (LES) and a combustion model were used. The combustion model was based on a two scalar flamelet approach, where a G-equation model and a conserved scalar model expressed a propagation of premixed flame surface and a diffusion combustion process, respectively. For validation of this numerical model, we have simulated the previous two types of hydrogen explosion tests. One is open-space explosion test, and the source was a prismatic 5.27 m3 volume with 30% of hydrogen-air mixture. A reinforced concrete wall was set 4 m away from the front surface of the source. The source was ignited at the bottom center by a spark. The other is vented enclosure explosion test, and the chamber was 4.6 m × 4.6 m × 3.0 m with a vent opening on one side. Vent area of 5.4 m2 was used. Test was performed with ignition at the center of the wall opposite the vent. Hydrogen-air mixtures with hydrogen concentrations close to 18% vol. were used in the tests. The results from the numerical simulations are compared with the previous experimental data for the accuracy of the numerical model, and we have verified that the simulated overpressures and flame time-of-arrival data were in good agreement with the results of the previous two explosion tests.
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: A solar receiver is designed for operation under
extremely uneven heat flux distribution, cyclic weather, and cloud
transient cycle conditions, which can include large thermal stress and
even receiver failure. In this study, the effect of different oil velocity
on convection coefficient factor and impact of wind velocity on local
Nusselt number by Finite Volume Method will be analyzed. This
study is organized to give an overview of the numerical modeling
using a MATLAB software, as an accurate, time efficient and
economical way of analyzing the heat transfer trends over stationary
receiver tube for different Reynolds number. The results reveal when
oil velocity is below 0.33m/s, the value of convection coefficient is
negligible at low temperature. The numerical graphs indicate that
when oil velocity increases up to 1.2 m/s, heat convection coefficient
increases significantly. In fact, a reduction in oil velocity causes a
reduction in heat conduction through the glass envelope. In addition,
the different local Nusselt number is reduced when the wind blows
toward the concave side of the collector and it has a significant effect
on heat losses reduction through the glass envelope.
Abstract: In this study, the three-dimensional cavitating
turbulent flow in a complete Francis turbine is simulated using
mixture model for cavity/liquid two-phase flows. Numerical analysis
is carried out using ANSYS CFX software release 12, and standard k-ε
turbulence model is adopted for this analysis. The computational
fluid domain consist of spiral casing, stay vanes, guide vanes, runner
and draft tube. The computational domain is discretized with a threedimensional
mesh system of unstructured tetrahedron mesh. The
finite volume method (FVM) is used to solve the governing equations
of the mixture model. Results of cavitation on the runner’s blades
under three different boundary conditions are presented and
discussed. From the numerical results it has been found that the
numerical method was successfully applied to simulate the cavitating
two-phase turbulent flow through a Francis turbine, and also
cavitation is clearly predicted in the form of water vapor formation
inside the turbine. By comparison the numerical prediction results
with a real runner; it’s shown that the region of higher volume
fraction obtained by simulation is consistent with the region of runner
cavitation damage.
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: This work consists of a numerical simulation of
convective heat transfer in a vertical plane channel filled with a heat
generating porous medium, in the absence of local thermal
equilibrium. The walls are maintained to a constant temperature and
the inlet velocity is uniform. The dynamic range is described by the
Darcy-Brinkman model and the thermal field by two energy
equations model. A dimensionless formulation is developed for
performing a parametric study based on certain dimensionless groups
such as, the Biot interstitial number, the thermal conductivity ratio
and the volumetric heat generation, q '''. The governing equations are
solved using the finite volume method, gave rise to a multitude of
results concerning in particular the thermal field in the porous
channel and the existence or not of the local thermal equilibrium.
Abstract: This paper presents a computational study of steady
state three dimensional very high turbulent flow and heat transfer
characteristics in a constant temperature-surfaced circular duct fitted
with 900 hemispherical inline baffles. The computations are based on
realizable k-ɛ model with standard wall function considering the
finite volume method, and the SIMPLE algorithm has been
implemented. Computational Study are carried out for Reynolds
number, Re ranging from 80000 to 120000, Prandtl Number, Pr of
0.73, Pitch Ratios, PR of 1,2,3,4,5 based on the hydraulic diameter of
the channel, hydrodynamic entry length, thermal entry length and the
test section. Ansys Fluent 15.0 software has been used to solve the
flow field. Study reveals that circular pipe having baffles has a higher
Nusselt number and friction factor compared to the smooth circular
pipe without baffles. Maximum Nusselt number and friction factor
are obtained for the PR=5 and PR=1 respectively. Nusselt number
increases while pitch ratio increases in the range of study; however,
friction factor also decreases up to PR 3 and after which it becomes
almost constant up to PR 5. Thermal enhancement factor increases
with increasing pitch ratio but with slightly decreasing Reynolds
number in the range of study and becomes almost constant at higher
Reynolds number. The computational results reveal that optimum
thermal enhancement factor of 900 inline hemispherical baffle is
about 1.23 for pitch ratio 5 at Reynolds number 120000.It also shows
that the optimum pitch ratio for which the baffles can be installed in
such very high turbulent flows should be 5. Results show that pitch
ratio and Reynolds number play an important role on both fluid flow
and heat transfer characteristics.
Abstract: Hypersonic flows around spatial vehicles during their reentry phase in planetary atmospheres are characterized by intense aerothermodynamics phenomena. The aim of this work is to analyze high temperature flows around an axisymmetric blunt body taking into account chemical and vibrational non-equilibrium for air mixture species and the no slip condition at the wall. For this purpose, the Navier-Stokes equations system is resolved by the finite volume methodology to determine the flow parameters around the axisymmetric blunt body especially at the stagnation point and in the boundary layer along the wall of the blunt body. The code allows the capture of shock wave before a blunt body placed in hypersonic free stream. The numerical technique uses the Flux Vector Splitting method of Van Leer. CFL coefficient and mesh size level are selected to ensure the numerical convergence.
Abstract: Reflux condensation occurs in vertical channels and tubes when there is an upward core flow of vapour (or gas-vapour mixture) and a downward flow of the liquid film. The understanding of this condensation configuration is crucial in the design of reflux condensers, distillation columns, and in loss-of-coolant safety analyses in nuclear power plant steam generators. The unique feature of this flow is the upward flow of the vapour-gas mixture (or pure vapour) that retards the liquid flow via shear at the liquid-mixture interface. The present model solves the full, elliptic governing equations in both the film and the gas-vapour core flow. The computational mesh is non-orthogonal and adapts dynamically the phase interface, thus produces a sharp and accurate interface. Shear forces and heat and mass transfer at the interface are accounted for fundamentally. This modeling is a big step ahead of current capabilities by removing the limitations of previous reflux condensation models which inherently cannot account for the detailed local balances of shear, mass, and heat transfer at the interface. Discretisation has been done based on finite volume method and co-located variable storage scheme. An in-house computer code was developed to implement the numerical solution scheme. Detailed results are presented for laminar reflux condensation from steam-air mixtures flowing in vertical parallel plate channels. The results include velocity and gas mass fraction profiles, as well as axial variations of film thickness.
Abstract: The thermal control in many systems is widely
accomplished applying mixed convection process due to its low cost,
reliability and easy maintenance. Typical applications include the
aircraft electronic equipment, rotating-disc heat exchangers, turbo
machinery, and nuclear reactors, etc. Natural convection in an inclined
square enclosure heated via wall heater has been studied numerically.
Finite volume method is used for solving momentum and energy
equations in the form of stream function–vorticity. The right and left
walls are kept at a constant temperature, while the other parts are
adiabatic. The range of the inclination angle covers a whole revolution.
The method is validated for a vertical cavity. A general power law
dependence of the Nusselt number with respect to the Rayleigh
number with the coefficient and exponent as functions of the
inclination angle is presented. For a fixed Rayleigh number, the
inclination angle increases or decreases is found.
Abstract: Organic Rankine Cycle (ORC) is the most commonly used method for recovering energy from small sources of heat. The investigation of the ORC in supercritical condition is a new research area as it has a potential to generate high power and thermal efficiency in a waste heat recovery system. This paper presents a steady state ORC model in supercritical condition and its simulations with a real engine’s exhaust data. The key component of ORC, evaporator, is modelled using finite volume method, modelling of all other components of the waste heat recovery system such as pump, expander and condenser are also presented. The aim of this paper is to investigate the effects of mass flow rate and evaporator outlet temperature on the efficiency of the waste heat recovery process. Additionally, the necessity of maintaining an optimum evaporator outlet temperature is also investigated. Simulation results show that modification of mass flow rate is the key to changing the operating temperature at the evaporator outlet.
Abstract: We have developed a new computer program in
Fortran 90, in order to obtain numerical solutions of a system
of Relativistic Magnetohydrodynamics partial differential equations
with predetermined gravitation (GRMHD), capable of simulating
the formation of relativistic jets from the accretion disk of matter
up to his ejection. Initially we carried out a study on numerical
methods of unidimensional Finite Volume, namely Lax-Friedrichs,
Lax-Wendroff, Nessyahu-Tadmor method and Godunov methods
dependent on Riemann problems, applied to equations Euler in
order to verify their main features and make comparisons among
those methods. It was then implemented the method of Finite
Volume Centered of Nessyahu-Tadmor, a numerical schemes that
has a formulation free and without dimensional separation of
Riemann problem solvers, even in two or more spatial dimensions,
at this point, already applied in equations GRMHD. Finally, the
Nessyahu-Tadmor method was possible to obtain stable numerical
solutions - without spurious oscillations or excessive dissipation -
from the magnetized accretion disk process in rotation with respect
to a central black hole (BH) Schwarzschild and immersed in a
magnetosphere, for the ejection of matter in the form of jet over a
distance of fourteen times the radius of the BH, a record in terms
of astrophysical simulation of this kind. Also in our simulations,
we managed to get substructures jets. A great advantage obtained
was that, with the our code, we got simulate GRMHD equations in
a simple personal computer.
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: The radiative heat transfer problem is investigated numerically for 2D complex geometry biomass pyrolysis reactor composed of two pyrolysis chambers and a heat recuperator. The fumes are a mixture of carbon dioxide and water vapor charged with absorbing and scattering particles and soot. In order to increase gases residence time and heat transfer, the heat recuperator is provided with many inclined, vertical, horizontal, diffuse and grey baffles of finite thickness and has a complex geometry. The Finite Volume Method (FVM) is applied to study radiative heat transfer. The blocked-off region procedure is used to treat the geometrical irregularities. Eight cases are considered in order to demonstrate the effect of adding baffles on the walls of the heat recuperator and on the walls of the pyrolysis rooms then choose the best case giving the maximum heat flux transferred to the biomass in the pyrolysis chambers. Ray effect due to the presence of baffles is studied and demonstrated to have a crucial effect on radiative heat flux on the walls of the pyrolysis rooms. Shadow effect caused by the presence of the baffles is also studied. The non grey radiative heat transfer is studied for the real existent configuration. The Weighted Sum of The Grey Gases (WSGG) Model of Kim and Song is used as non grey model. The effect of soot volumetric fraction on the non grey radiative heat flux is investigated and discussed.