Modeling a Coil Heat Exchanger

August 21, 2014

Coil heat exchangers are simple and easy to manufacture. Here, we consider an axially wound coil of copper carrying hot water that heats air inside of a circular duct. Since the geometry is almost invariant about the centerline, the model is solved in the 2D axisymmetric plane. Additional expressions are added to compute the temperature drop between turns of the coil, which greatly simplifies the modeling.

A Coil Heat Exchanger Model

The coil heat exchanger we’ll consider is shown in the figure below.

Model of coil heat exchanger.
A copper coil carries hot water through a duct carrying cold air.

Copper tubing is helically wound so that it can be inserted along the axis of a circular air duct. Cold air is moving through the duct, and hot water is pumped through the tubing. The air flow pattern and the temperature of the air and copper pipes will be computed using the Conjugate Heat Transfer interface. Since the geometry is almost axisymmetric, we can simplify our modeling by assuming that the geometry and the air flow are entirely axisymmetric. Thus, we can use the 2D axisymmetric Conjugate Heat Transfer interface. Since the airspeed is high, a turbulent flow model is used; in this case, it is the k-epsilon model.

We can assume that the water flowing inside of the pipe is a fully developed flow. We can also assume that the temperature variation of the water is small enough that the density does not change, hence the average velocity will be constant. Therefore, we do not need to model the flow of the water at all and can instead model the heat transfer between the fluid and the pipe walls via a forced convective heat transfer correlation.

The Convective Heat Flux boundary condition uses a Nusselt number correlation for forced internal convection to compute the heat transfer between the water and copper tubing. This boundary condition is applied at all inside boundaries of the copper piping. As inputs, it takes the pipe dimensions, fluid type, fluid velocity, and fluid temperature. With the exception of the fluid temperature, all of these quantities remain constant between the turns of the tubing.

Modeling the Temperature Drop Between Turns

As the hot water is being pumped through the copper coils, it cools down. However, since the model is axisymmetric, each turn of the coil is independent of the others, unless we explicitly pass information between them. That is, we must apply a separate Convective Heat Flux boundary condition at the inside boundaries of each coil turn.

This raises the question: How do we compute the temperature drop between each turn and incorporate this information into our model?

Consider the water passing through one turn of the copper coil. The heat lost by the water equals the heat transfer into the copper pipes. Under the assumption of constant material properties, and neglecting viscous losses, the temperature drop of the water passing through one turn of the pipe is:

\Delta T = \frac{Q}{\dot m C_p} = \frac{\int q'' dA}{\dot m C_p}

where \dot m is the mass flow rate, C_p is the specific heat of water, and Q is the total heat lost by the water, which is equal to the integral of the heat flux into the copper, integrated over the inside boundaries of the coil. This integral can be evaluated via the Integration Component Coupling, defined over the inside coil boundaries.

Integration Component Coupling in COMSOL Multiphysics.
The Integration Component Coupling defined over a boundary. Note: The integral is computed in the revolved geometry.

Using these coupling operators, we can define a set of user-defined variables for the temperature drop:

DT1 = intop1(-nitf.nteflux/mdot0/Cp0)

This evaluates the temperature drop along the first turn of the pipe. We can define a different temperature drop variable for each turn of the pipe and use them sequentially for each turn.

Defining different temperature drop variables in software.
The water temperature in the sixth turn considers the temperature drop in the first five turns.

Plots show relationship between flow field and temperature.
Flow field and temperature plot (left) and the temperature along a line through the center of the coils (right).

Since this is a 2D axisymmetric model, it will solve very quickly. We can examine the temperature and the flow fields and plot the temperature drop along a line down the center of the coils. We can observe that the water cools down between each turn of the coil, and the air heats up.

This can be considered a parallel-flow heat exchanger, since the hot and cold fluids flow in the same overall direction. If we wanted to change this model to the counter-flow configuration, we could simply switch the air inlet and outlet conditions so that the fluids travel in opposite directions.

What other kinds of heat exchanger configurations do you think this technique can be applied to?


Comments (13)

Leave a Comment
Log In | Registration
Loading...
Simon Stubkier
Simon Stubkier
September 8, 2014

Very interesting blog – would it be possible to have a look at the model file?

I am currently working on a very similar problem – but in 3D. Do you have any experiance in meshing 3D geometries similar to this?

Could you provide the cylindrical 3D plot of the temperatures in the coil ?

Br

Walter Frei
Walter Frei
September 8, 2014 COMSOL Employee

Dear Simon,
Please contact your local COMSOL Support Staff, and they will be happy to address your questions.
Best,
Walter

amin aram
amin aram
November 29, 2017

Dear Bro!
I need the Comsol file if this project. is it possible?
amin.hadidi@yahoo.com
@amin_aram telegram

Schubert Su
Schubert Su
February 7, 2019

Dear Walter Frei,
I also need the COMSOL file for this project. Is it possible?
m91720029@email.ntou.edu.tw

thank you

Long Chen
Long Chen
May 27, 2019

Dear Bro!
I also need the Comsol file of this project. is it possible?
18209291445@163.com
thank you

SME Coils
SME Coils
September 23, 2020

Your topic is very nice and helpful to us … Thank you for the information you wrote.

YGS JAGAJO
YGS JAGAJO
June 10, 2022

This far we could conclude that pressure drop should be calculated too? since for a reason, and with such angle with the coil bend, pressure drop might even affected by the temperature we leave by, isn’t it?

Walter Frei
Walter Frei
June 10, 2022 COMSOL Employee

We do assume here that the fluid inside the pipe, water, has constant density (is incompressible) so pressure should not enter the problem, for this case.

Dang Hai
Dang Hai
September 4, 2022

Dear Walter Frei
In the figure above, the coil structure is not axisymmetric, but the assume it is axisymmetric. Do you have any reference taking how we can assume it is axisymmentric?

Walter Frei
Walter Frei
September 6, 2022 COMSOL Employee

Dear Dang Hai, The nice thing about assumptions is that we can make them without having to find references to justify them! Now, people are free to disagree with your assumptions, of course. One approach to settle the matter is to build both a 2D and 3D model, and compare.

Adnan Shaikh
Adnan Shaikh
December 4, 2022

Dear Walter,
Kindly send me comsol file of this problem on
arkp.adnan@gmail.com

Thank You

Amber Zaheer
Amber Zaheer
October 30, 2023

Hi,

Can I please have the COMSOL file for this model sent at amber.zaheer@manchester.ac.uk? Compatible with COMSOL 6.0.

Many thanks,
Amber.

Ruben Giovani Tambunan
Ruben Giovani Tambunan
September 7, 2024

Dear Walter Frei,
I also need the COMSOL file for this project. Is it possible?
rubengvnt01@gmail.com

thank you

EXPLORE COMSOL BLOG