ACMHT Lalot 1998
ACMHT Lalot 1998
ACMHT Lalot 1998
net/publication/260712323
CITATIONS READS
0 625
2 authors, including:
Jean-Michel Bergheau
École nationale d'ingénieurs de Saint-Étienne
496 PUBLICATIONS 4,127 CITATIONS
SEE PROFILE
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Jean-Michel Bergheau on 12 March 2014.
Abstract
This paper presents a new finite element used to model heat exchangers and
implemented in SYSTUS+. The fluid is represented by a set of one-dimensional
elements, and the wall of the pipe, in which the fluid flows, by 2D elements. One
element in the fluid is connected with several wall elements. It is well known
that the usual Galerkin method for the finite element modeling of
diffusion-convection problems leads to spatial oscillations for high Peclet
numbers, typically, with first order elements, when the Peclet number is higher
than 2. To avoid these oscillations, the formulation in the fluid uses a
Petrov-Galerkin method involving discontinuous test functions. The heat transfer
computation is based on the average fluid temperature (as the usual correlations
do) but on the local wall temperature. All these elements are fully coupled with
the other elements of the studied structure for conduction, standard convection
and radiation. The computation of the heat transfer coefficient uses various
models. All these models represent the generalization of well known
correlations. It is possible to simulate free convection, forced convection, mixed
convection and change of phase. The efficiency of this approach is discussed on
numerical applications concerning : a condenser, a counterflow heat exchanger, a
1-2 heat exchanger, natural convection over a radiating surface. Very accurate
results are obtained in all the cases.
Fluid element
Macro-element
II Theoretical background
II-1 General assumptions
(1)
(2)
(3)
A
[K f ] = fluid [K ] represents the classical thermal conductivity matrix,
e
f
elements
(4)
(5)
(6)
From assumption A3 and equation (1), the heat balance inside the kth
macro-element gives :
(7)
The non linear coupled problem constituted by equations (3) and (8) can
be solved using a Newton-Raphson procedure. First, equations (3) and (8)
are re-written as :
We give here the main models that are available. They represent :
⎢ ⎜ ⎝ Pr ⎠ ⎟ ⎥
⎣ ⎝ ⎠ ⎦
- forced convection : Nu = b1 Pr b2 Re b3 ; Nu = b1 + b2 Pr b3 Gr b4 + b5 Pr b6 Gr b7
[
- mixed convection : Nu = b1 b2 Re b3 + Gr b4 ]
b5
Pr b6
Other models are available for change of phase, critical heat flux
calculation, natural convection, forced convection, user defined model
(via a FORTRAN subroutine).
III Validation
III-1 Study of a condenser
Firstly, we study the case where all the fluid characteristics are constant.
In this case we have (INCROPERA and DE WITT [3]) :
This shows that, in this case, SYSTUS+ gives very accurate results, even
if the Peclet number is about 2 106.
We know that if the oil has to be at 60°C at the outlet of the exchanger,
the length of the tube is 66.5 meters. The water outlet temperature is then
40.2°C, the convection coefficients are 2250 W/m².K (for the water) and
38.4 W/m².K (for the oil). In SYSTUS+, the water is supported by 50 1D
elements, the oil is supported by 50 1D elements, the wall is represented
by two sets of 2000 2D elements. One set of 2D elements is connected to
the water, the other set is connected to the oil, and each node of one set is
linked to the corresponding node of the other set of elements to get a
unique temperature at each point of the surface. We only have to give the
inlet temperatures, and SYSTUS+ computes the following data :
Water Oil
Outlet Convection Outlet Convection
temperature coefficient temperature coefficient
40.18 °C 2249 W/m².K 60.06 °C 38.36 W/m².K
The 1-2 heat exchanger we study here, is 10 meters long, has one pass for
the oil (represented by 20 1D elements) and 2 passes (5 tubes per pass)
for the water (represented by 20 1D elements per tube). Each tube is
represented by 20 two sets of 200 2D elements : one set for the
convection with the water and one set for the convection with the oil. The
We study here a flat plate. This plate (0.5m x 0.5 m) is heated by one
side, and cooled by the upper side by natural convection and radiation (ε
= 0.7). The temperature of the infinite medium is T∞ = 27 °C . We know
that in this case the surface temperature is the solution of :
IV Conclusion
We have shown that we obtain very accurate results with the new
element. Even if this element has been developed for forced convection
in heat exchangers, we have shown that we still obtain good results for
free convection. We have also tested the element for electrical heaters in
transient states and we have obtained very accurate results (LALOT and
KALUBA [5]). Future studies will concern the transient behavior of heat
exchangers.
Acknowledgments
The authors wish to thank VULCANIC company for its technical and
financial support.
References
[1] SYSTUS+ V1.1, User’s manual, SYSTUS International, 1997.
[2] Hughes T.J.R. and Brooks A., A theoretical framework for Petrov-
Galerkin methods with discontinuous weighting functions :
application to the streamline-upwind procedure, Finite Element in
Fluids, Vol. 4, Chap. 3, 1982, R.H. Gallagher Ed.