(PDF) ANSYS Fluent CFD Modeling of Solar Air-Heater Thermoaerodynamics
(PDF) ANSYS Fluent CFD Modeling of Solar Air-Heater Thermoaerodynamics
Dmitry Pashchenko
Guangdong Technion - Israel Institute of Technology
Abstract
Numerical modeling of thermal and aerodynamic processes taking place in a solar air
heater with light-absorbing L-shaped fins is performed. The study is carried out by means Discover the world's
research
of an ANSYS Fluent Solver with an integrated solar calculator for the city of Samara,
Russia. The influence of the design (the step between the fins) and technological variables
25+ million
(the Reynolds number) on heat-exchange processes and flow aerodynamics was members
established during the CFD analysis. Pictorial contours are determined for the velocity
distribution, pressure, and other characteristics of the airflow in the solar air-heater box. 160+ million
Thermoaerodynamic characteristics and effective values of technological parameters for publication
the effective use of a solar air heater are studied. The maximum thermoaerodynamic pages
characteristic for the solar air heater with a finned lightabsorbing surface is 1.91 for Re = 2.3+ billion Join for free
1500, and the minimum test fin pitch is 30 mm. citations
Public Full-text 1
ISSN 0003-701X, Applied Solar Energy, 2018, Vol. 54, No. 1, pp. 32–39. © Allerton Press, Inc., 2018.
Original Russian Text © D.I. Pashchenko, 2018, published in Geliotekhnika, 2018, No. 1, pp. 37–47.
Abstract—Numerical modeling of thermal and aerodynamic processes taking place in a solar air heater with
light-absorbing L-shaped f ins is performed. The study is carried out by means of an ANSYS Fluent Solver
with an integrated solar calculator for the city of Samara, Russia. The inf luence of the design (the step
between the fins) and technological variables (the Reynolds number) on heat-exchange processes and f low
aerodynamics was established during the CFD analysis. Pictorial contours are determined for the velocity
distribution, pressure, and other characteristics of the airf low in the solar air-heater box. Thermoaerody-
namic characteristics and effective values of technological parameters for the effective use of a solar air heater
are studied. The maximum thermoaerodynamic characteristic for the solar air heater with a finned light-
absorbing surface is 1.91 for Re = 1500, and the minimum test fin pitch is 30 mm.
Download full-text PDF Download citation Copy link
DOI: 10.3103/S0003701X18010103
32
Inlet Outlet
40
20
40
I II III
200
Fig. 1. Computational domain: (I) air-inlet zone, (II) heat-exchange zone, and (III) outlet zone.
34 PASHCHENKO
(a)
(b)
Fig. 2. Structure of the CFD grid: (a) appearance of the CFD grid of the air channel; (b) boundar y layer near the fin surface.
(a) (b)
60
Roughness coecient 10–3,
40 10
30
8
20
10 6
5 10 15 20 5 10 15 20
Reynolds number 104, Re Reynolds number 104, Re
Fig. 4. Comparison of experimental data and the results of CFD modeling for the smooth light-absorbing surface (unfinned).
36 PASHCHENKO
Velocity contour 1
8.555e–1
8.020e–1
7.485e–1
6.951e–1
6.416e–1
5.881e–1
5.347e–1
4.812e–1
4.277e–1
3.743e–1
3.208e–1
2.673e–1
2.139e–1
1.604e–1
1.069e–1
5.347e–2
0.000e+0
[m s–1]
Fig. 5. Airf low-velocity contour for the surface with a 40-mm fin pitch for Reynolds number Re = 1200.
Fig. 6. Vector field of the airf low velocity for the fir st five fins with a pitch of 40 mm along the airf low, Re = 1200.
on the Reynolds number for light-absorbing surfaces Figure 9 shows the dependence of the air tempera
with different fin pitches obtained by the developed ture drop at the inlet and outlet of the solar heater fo
model for various Reynolds numbers 400 > Re> 2000. various shapes of the light-absorbing surface durin
the day. As initial data for the analysis of the inf luenc
It can be seen from Fig. 8 that the presence of fins of these parameters on the air heating value, the fol
significantly increases the heat exchange between the lowing values were used: date March 20, 2017; Samar
light-absorbing surface and the airf low in comparison (GMT + 4); absolutely clear sky (sunshine factor = 1)
with a smooth (unf inned) surface. Thus, for a 40-mm Re = 800; level angle of the solar air heater to th
ground surface α = 60°.
fin pitch, the Nusselt number increases by more than
three times compared to a smooth surface. Figure 8 From Fig. 9 it can be seen that the maximum ai
heating is observed in the interval between 1 and 2 pm
also shows that the average Nusselt number increases This is because the angle cosin e between the sunbeam
with the increase of the f in pitch. This is because the and the SAH surface at this time has the maximum
increase in the distance between the fins leads to a value for the middle of the spring in Samara. The max
decrease in the kinetic energy of the turbulent flow. imum air heating value of 25.4°C is achieved for th
Fig. 7. Contour of the air-f low-turbulence kinetic energy for the surface with a fin pitch of 40 mm, Re = 1200.
(a) (b)
50 Unfinned 50 Re = 2000
30 mm Re = 1600
40 mm Re = 1200
40 40
Nusselt number, Nu
Nusselt number, Nu
50 mm Re = 800
60 mm Re = 400
30 30
20 20
10 10
0 0
5 10 15 20 30 40 50 60
Reynolds number ×103, Re Pitch, mm
Fig. 8. Variation of the Nusselt number with the Reynolds number for various fin pitches.
minimum inlet air velocity and the minimum 30-mm where ΔР is pressure losses in SAH, Pa; L is a charac
test pitch. teristic length, m; deqv is an equivalent diameter, m
ρ is air density, kg/m3; and u is air velocity, m/s.
According to Fig. 10, th e L-shaped f in setting on
ROUGHNESS COEFFICIENT light-absorbing surface of a solar heater creates a
additional roughness that destroys the laminar wal
One of the main parameters taking into account layer. This leads to the appearance of extra aerody
the efficiency of using the f inned light-absorbing sur- namic resistance and, as a result, to a signif ican
face is a roughness coeff icient. Figure 10 shows a increase in the roughness coeff icient in comp arison
graph of variance of the friction coefficient for a with a smooth surface. So for Re = 8000, the rough
finned light-ab sorbing surface having a different pitch ness coefficient is 3.8–4.2 times higher for differen
and for a smooth surface. The friction coefficient in fin pitches than for a smooth (unf inned) surface. A
increase in the Reynolds number causes suppression
the present CFD experiment is determined from the of the laminar sublayer and, consequently, the value o
Darcy–Weisbach formula [19, 20]: the roughness coeff icient decreases with an increase in
the Reynolds number, as shown in Fig. 10b. As dis
ΔР
λ= , (6) cussed above, the increase in the Reynolds numbe
L u2 increases the turbulence intensity and, as a result, th
ρ heat exchange grows with a smaller value of the fric
d eqv 2
tion coefficient.
A pressure contour obtained with the CFD analysi
30 for Re = 1200 and a 40-mm fin pitch is shown in Fig. 11
30 mm Figure 11 shows that the maximum aerodynami
Temperature variation, C
40 mm
resistance along the airf low is offered by a fin. Mini
50 mm
60 mm mization of this resistance is possible du e to the use o
20 a fir st fin of another shape, for example, an L-shap e
fin, inverted in the airf low direction.
10 THERMOAERODYNAMIC CHARACTERISTIC
It was found in the CFD study that was carried ou
that the artificial roughness in the form of additiona
0 fins on a light-abs orbing surface greatly enhances th
10 11 12 13 14 15 16 heat exchange between th e wall and the heated ai
Time, h Nevertheless, the additional roughness leads to highe
aerodynamic resistance and, as a consequence, to a
Fig. 9. Air-temperature variation during the day for differ- increase in the energy consumption for air pumpin
Downloadent light-absorbing
full-text PDF surfaces. Download citation Copythrough
link the SAH air channel. Therefore, for an objec
APPLIED SOLAR ENERGY Vol. 54 No. 1 2018
38 PASHCHENKO
(a) (b)
50 Unfinned 40
Roughness coecient 10–2,
Pressure contour 1
2.334e+1
1.833e+1
1.331e+1
8.298e+0
3.284e+0
–1.731e+0
–6.745e+0
–1.176e+1
–1.677e+1
–2.179e+1
–2.680e+1
–3.182e+1
–3.683e+1
[Pa]
Fig. 11. Pressure contour in the air channel of the solar air heater for Re = 1200 and a 40-mm f in pitch.
tive analysis, it is necessary to evaluate a generalized increase in the Reynolds number is usually associate
thermoaerodynamic characteristic of the simulated air with an increase in airf low velocity, which has a signif
heater, and not only thermal and aerodynamic ones. icant effect on the roughness coefficient in accordanc
A characteristic of the system based on the interre- with relationship (7).
lation between thermal and aerodynamic parameters
is called a “thermoaerodynamic” characteristic, for 2.0
which quantitative evaluation of which the following 60 mm
expression is used [21]: 50 mm
Thermoaerodynamic characteristic
40 mm
Nu Nu uf
. η= (7) 1.9
30 mm
f uf ) (f
1/3
... In recent years, models that are significantly superior to traditional ones both in terms of computational efficiency and
in terms of accuracy of solving a wide range of forced convection problems have gained widespread popularity. These
include, in particular, a one-parameter model based on the use of the transport equation for turbulent viscosity, the SA
model [23] , and the two-parameter Menter model [24], which is a combination of the standard model and famous
model. The wide capabilities of these models and their significant superiority over standard models in relation to solving
a number of complex problems of forced convection have been demonstrated in many works. ...
... where are the velocity components, and in the momentum equation P is the pressure, μ and μ t are the dynamic and
turbulent viscosity of air, respectively, and in the energy equation T, and are the temperature, molecular thermal
diffusivity, and turbulent thermal diffusivity of air, respectively and i, j = 1, 2, and 3. Turbulence is modeled using the SA
model. The oneequation model is specified as follows [22, 23, 30]: (4) and turbulent eddy viscosity is calculated by
formula: ...
... (5) Other values are presented in article [23] . Continuity, momentum, and energy equations are given to analyze the
passage of air flow through a solar air heater, since they allow us to describe the main physical processes occurring in
the system. ...
Simulation of Turbulent Natural Convection in Photovoltaic Solar Panels Based on the Spalart–Allmares (SA) Turbulence
Model
View
... Solar energy system is widely investigated by researchers with many tools, ANSYS-FLUENT was used to investigate
the thermodynamic properties and efficiency of the solar system. For example, the thermal and aerodynamic processes
occurring in a solar air heater with L-shaped fins that absorb light are modeled numerically [30] . For the Russian city of
Samara, the analysis is conducted using an ANSYS-FLUENT Solver with an inbuilt solar calculator. ...
... The MATLAB/SIMULINK is widely used to model and validation of Photovoltaic solar cell to produce electricity,
pumping system, and many electric devices like Figure 13. CFD model and grid [30] . ...
... Velocity contour and velocity vector at Re = 1200 [30] . electric vehicles. ...
... However, with the continual evolution of computational fluid dynamics (CFD) technology, CFD software platforms
such as ANSYS Fluent [43] , OpenFOAM [44], and Phoenics [45] have emerged as indispensable analytical tools in
aerospace engineering. These advanced computational tools enable engineers to simulate and analyze complex
aerodynamic phenomena with unprecedented accuracy and efficiency, facilitating the optimization of aerospace vehicle
designs and enhancing their performance characteristics. ...
Article
Apr 2024
Yining Zhang · Yifan Xie · Guanyixuan Zhao · Yang Yang
... With industrialisation and the advancement of modern technology, the global energy protocol has been facing
unprecedented challenges in coping with the current energy needs, particularly due to the overconsumption of natural
resources and also with the related adverse impacts to the nature. Phase change materials (PCMs), which can absorb
and release thermal energy during the heat transfer and conversion process via phase transitions, such as paraffin wax
and polyethylene glycol (PEG), are potential smart energy options that can be utilised in energy storage applications,
heat transfer systems (by diffusion of thermal energy) [1][2][3][4], and solar-thermal energy conversion applications
[5, 6]; and also they can be selected based on different melting temperatures to suit multiple applications [7,8]. For
example, at low temperatures (< 20 • C), they can be used in industrial refrigeration, dissipation for spacecraft systems,
as well as heat and cold storage units in buildings [9]. ...
Enhancing the thermal performance of polyethylene glycol phase change material with carbon-based fillers
... In this investigation, quality check criteria were set to monitored with skewness correlation and aspect ratio. The
polyhedral grid type has an optimum performing skewness factor of 0.89 and an aspect ratio of 0.15 which is
recommended by ANSYS fluent manual [37, 38]. Furthermore, the selected grid produces a lower cell count which
drastically reduces computational time as shown in Table 2. ...
Battery thermal runaway propagation time delay strategy using phase change material integrated with pyro block lining: Dual
functionality battery thermal design
Article
Apr 2023
Virendra Talele · Mahesh Suresh Patil · Satyam Panchal · Roydon Fraser
... Pashchenko conducted a CFD-based analysis for SAE using ANSYS-FLUENTand concluded that the optimal efficiency
occurs at an inclination angle of 60 degrees between the absorption surface and the Earth [19]. Moreover, he discovered
that using an L-shaped fin on SAE with a pitch of 30 mm and Re of 1500 improves its thermodynamic performance
[20] . Using MATLAB tools, the finned and roughened SAE designs are examined, and the design configurations are
optimized based on energy and exergy efficiency [21,22]. ...
Genetic Algorithm Integrated Fuzzy AHP-VIKOR Approach for the Investigation of W-Cut Insert Heat Exchanger for Cooling
of Dielectric Fluid Used in Ultra-High Voltage Transformer
A highly extensible smoothed particle hydrodynamics with meshless fluid-end general interpolation method based on
precise code interaction coupling environment for fluid–structure interaction surface coupling
Article
Sep 2024
Sifan Long · Xiao-Wei Guo · Xiaokang Fan · Canqun Yang
3D numerical model of a Concentrated photovoltaic thermal (CPV/T) system for thermal and electrical performance
optimization
Article
Jul 2024
Jhanbux Manek Variava · Jayesh K. Ratnadhariya · Mohd. Irfanul Haque Siddiqui · Kishor Kumar Sadasivuni
View
Experimental and Numerical Analysis of the Solar Collector with Stainless Steel Scourers Added to the Absorber Surface
Mar 2024
Ömer FARUK Can · Nevin Celik · Filiz Ozgen · Ali Taskiran
Experimental investigation of double-glazed double-pass solar airheater (DG-DPSAH) with multi-v ribs having trapezoidal
roughness geometry
Article
May 2022
Swati Mor · Niraj Kumar · Gulshan Sachdeva
Show more
Download full-text PDF Download citation Copy link
Experimental and numerical studies on flow from axisymmetric nozzle flow with sudden expansion for m...
January 2016 · International Journal of Energy, Environment and Economics
Researches over late in the field of aerodynamic vehicles have been concerned with the problem of flow separation over its base
which has led to the establishment of low pressure circulation region. This pressure is found to be noticeably lower in comparison with
the free stream atmospheric pressure. Base drag due to such pressure differences can be up to two-thirds of the total drag which is ...
[Show full abstract]
View full-text
Experimental and Numerical Studies On Flow From Axisymmetric Nozzles Flow With Sudden Expansion For...
Researches over late in the field of aerodynamic vehicles have been concerned with the problem of flow separation over its base
which has lead to the establishment of low pressure circulation region. This pressure is found to be noticeably lower in comparison
with the free stream atmospheric pressure. Base drag due to such pressure differences can be up to two-thirds of the total drag which
is ... [Show full abstract]
View full-text
Experimental and Numerical Studies on flow from Axisymmetric Nozzle Flow With Sudden Expansion for M...
February 2016
Jaimon Dennis Quadros · Sher Afghan Khan · Jolene S Vas
Researches over the recent past in the field of aerodynamic vehicles have been concerned with the problem of flow separation over its
base which has lead to the formation of low pressure circulation region. This pressure is found to be noticeably lower than the free
stream atmospheric pressure. Base drag due to such pressure differences can be up to two-thirds of the total drag and is principally ...
[Show full abstract]
View full-text
Article
Dmitry Pashchenko
Rapid development of the technologies of use of renewable energy sources in the recent decade has opened up prospects for
satisfying our society′s ever growing needs for fuel and power resources. In the present work, the author has shown results of
investigation into the thermal and aerodynamic characteristics of a solar air heater. A computer model whose adequacy has been
confirmed using ... [Show full abstract]
Read more
Article
Experimental and CFD investigation of convection heat transfer in solar air heater with reverse L-sh...
A solar air heater is a thermal system which uses artificial roughness in the form of repeated ribs on the absorber plate to enhance the
heat transfer rate. Forced convection heat transfer of air in a solar air heater with reverse L-shaped ribs has been carried out
experimentally and numerically. Thermal performance of solar air heater is studied with design variables such as relative roughness ...
[Show full abstract]
[Show full abstract]
Download full-text PDF Download citation Copy link
Read more
Article
Thermal performance evaluation of solar air heater using combined square and equilateral triangular...
September 2018 · Australian Journal of Mechanical Engineering
Heat transfer and fluid flow characteristics in roughened solar air heater duct with combined square and equilateral triangular ribs as
roughness is investigated using experimental and computational fluid dynamics (CFD) approaches. Relative roughness pitch (7.14 ≤
P/e ≤ 17.86), roughness height (e/D = 0.042), Reynolds number (3800≤ Re ≤ 18,000) and heat flux (1000 W/m²) are used as operating
... [Show full abstract]
Read more
Article
Experimental and CFD-based thermal performance prediction of solar air heater provided with chamfere...
August 2015 · Journal of the Brazilian Society of Mechanical Sciences and Engineering
An experimental and two-dimensional computational fluid dynamics (CFD) analysis of a solar air heater has been carried out using
right-angle triangular ribs as artificial roughness on the absorber plate as a convenient method for enhancement of thermal
performance of solar air heater. The relative roughness pitch (P/e = 7.14–35.71), Reynolds number (Re = 3800–18000), and relative
roughness height ... [Show full abstract]
Read more
© 2008-2024 ResearchGate GmbH. All rights reserved. Terms · Privacy · Copyright · Imprint · Consent preferences