CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Post-Processing

chtMultiRegionFoam wallHeatFlux mismatch

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 7, 2019, 11:47
Default chtMultiRegionFoam wallHeatFlux mismatch
  #1
New Member
 
Join Date: Sep 2019
Posts: 8
Rep Power: 6
foamyva is on a distinguished road
Dear CFD-ehtousiasts,

Currently I am working on a model that should capture the convective heat transfer between a coil, with hot water running through it, to ambient tankwater. With this model I try to capture the local and total heat transfer.

The first setup of the model includes a simple single tube surrounded by tankWater.

I am extracting heat flux values on both sides of my coil and in both directions (just for checking). These value nicely match.
Code:
wallHeatFlux1
{
    type        wallHeatFlux;
    libs        ("libfieldFunctionObjects.so");
    patches     ("coilWall_to_coilWater");
    region      coilWall;
    executeControl  writeTime;
    writeControl    writeTime;
}

wallHeatFlux2
{
    type        wallHeatFlux;
    libs        ("libfieldFunctionObjects.so");
    patches     ("coilWater_to_coilWall");
    region      coilWater;
    executeControl  writeTime;
    writeControl    writeTime;
}

wallHeatFlux3
{
    type        wallHeatFlux;
    libs        ("libfieldFunctionObjects.so");
    patches     ("coilWall_to_tankWater");
    region      coilWall;
    executeControl  writeTime;
    writeControl    writeTime;
}

wallHeatFlux4
{
    type        wallHeatFlux;
    libs        ("libfieldFunctionObjects.so");
    patches     ("tankWater_to_coilWall");
    region      tankWater;
    executeControl  writeTime;
    writeControl    writeTime;
}
Now, the problem:

If I compare the wall heat flux values from the wallHeatFlux utility, to the value computed from the flow rate and temperature at the inlet and outlet as:
Qdot=(mdot*rho*Cp*(inletT-outletT))

I find a significant mismatch between both values (see figure). BTW: I used patchAverage to compute my mean outletT.

Unfortunately, I could not find a solution in other threads on the wallHeatFlux utility.

Any ideas or tips would be much appreciated!
Attached Images
File Type: png wallHeatFlux comparison.png (14.3 KB, 87 views)

Last edited by foamyva; October 8, 2019 at 04:57. Reason: Attached the wrong figure
foamyva is offline   Reply With Quote

Old   October 9, 2019, 03:25
Default Update
  #2
New Member
 
Join Date: Sep 2019
Posts: 8
Rep Power: 6
foamyva is on a distinguished road
I ran the simulation with turbulence turned off. Surprisingly (or not) the wall heat flux and flowrate based heat flux to almost match perfectly (see figure).

Now, I am using an KOmegaSST turbulence model. In order to capture the heat transfer as good as possible I have meshed my system in such a way that y+ is smaller or equal to one everywhere and therefore not using any wall functions.

Any experience on this aspect what could cause the energy imbalance?
Attached Images
File Type: png wallHeatFlux comparison no_turb.png (12.6 KB, 71 views)
foamyva is offline   Reply With Quote

Old   October 9, 2019, 10:38
Default
  #3
Senior Member
 
Join Date: Sep 2013
Posts: 353
Rep Power: 20
Bloerb will become famous soon enough
Your formula is only an approximation instead of temperature you'd need to look at the total Enthalpy between inlet and outlet. And you'd need to integrate instead of average everything.You would also need to include kinetic energy, although this part should probably nearly cancel out or be insignificant. Unless maybe the diameter changed drastically and the flow is quite fast.
Are any of your material or transport properties temperature dependent? Do the wall heat fluxes from the different regions match?
So is the coilWall_to_coilWater heatflux nearly identical to the coilWater_to_ coilWall flux?
Bloerb is offline   Reply With Quote

Old   October 9, 2019, 11:23
Default
  #4
New Member
 
Join Date: Sep 2019
Posts: 8
Rep Power: 6
foamyva is on a distinguished road
Thanks for your help! As mentioned in the original post I checked the heat flux on each side and they do match. Regarding the averaging of the temperature, patchAverage uses areaAverage and therefore does a weighted average calculation. As far as I know this is equivalent to integration, right?

I understand your point about the total enthalpy. However I do not believe that the internal energy/velocity would change so significant that it results in the given energy mismatch. It is still worth checking. Best way would to this would be to check the enthalpy flux over the inlet and outlet? I will have to look into this.

I do use temperature dependent material properties, however I do not believe the problem can be found here, as the no-turbulence model gives good results with the same material properties.

BTW: tube is horizontal and has a uniform cross-section.
foamyva is offline   Reply With Quote

Old   October 9, 2019, 13:22
Default
  #5
Senior Member
 
Join Date: Sep 2013
Posts: 353
Rep Power: 20
Bloerb will become famous soon enough
This is the full energy equation (not a single simplification made):

0 = \nabla\cdot\left(\left(e+\frac{1}{2}\mathbf{v}\cdot\mathbf{v}\right)\rho\mathbf{v}\right)-\nabla\cdot\left(\lambda\nabla T\right)+\nabla\cdot\left(\mathbf{v}p\right)-\underbrace{\rho\mathbf{g}\cdot\mathbf{v}-\rho\dot{q}-\nabla\cdot\left(\boldsymbol{\tau}\cdot\mathbf{v}\right)}_{0}
The last part of it assumed as zero since
  • gravitational acceleration is not active
  • no internal heat source is active inside the fluid
  • you are not adding viscous dissipation via fvOption. This is often negligible, but needs to be accounted for if you are simulating highly viscous substances.
With that said we are left with only these terms:
0 =  \nabla\cdot\left(\left(e+\frac{1}{2}\mathbf{v}\cdot\mathbf{v}\right)\rho\mathbf{v}\right)-\nabla\cdot\left(\lambda\nabla   T\right)+\nabla\cdot\left(\mathbf{v}p\right)

Now we apply the Gauss theorem:

\int \nabla\cdot\mathbf{f}dV = \int \mathbf{n}\cdot\mathbf{f}dA

Which yields:

0 = \int \mathbf{n}\cdot\left(\left(e+\frac{1}{2}\mathbf{v}\cdot\mathbf{v}\right)\rho\mathbf{v}-\lambda\nabla    T+\mathbf{v}p\right)\, dA

This equals:

0 = \int \mathbf{n}\cdot\left(H\rho\mathbf{v}-\lambda\nabla    T\right)\, dA

where

H = h+\frac{1}{2}\mathbf{v}\cdot\mathbf{v}

and the pv term was used to transform e to h. We now need to evaluate that term on the entire boundary of our domain. I'll split it into inlet/outlet and walls. With the difference being the assumption that nothing flows through the wall (n.v=0)


\underbrace{\int_{wall} \mathbf{n}\cdot\left(\lambda\nabla    T\right)\, dA_{wall}}_{\dot{Q}} = \int_{inout} \mathbf{n}\cdot\left(H\rho\mathbf{v}-\lambda\nabla    T\right)\, dA_{inout}


Under the assumption that no temperature gradient is present at inlet and outlet and no change in kinetic energy this equals your formula


\dot{Q}= \int_{inout}  \mathbf{n}\cdot\left(h\rho\mathbf{v}\right)\,  dA_{inout}= \int_{inout}  \mathbf{n}\cdot\left(c_p T \rho\mathbf{v}\right)\,  dA_{inout}


Now you can check if either of these assumptions is not the case for your simulation. And no you need to integrate since T is multiplied and hence weighted by v and cp if neither of these 3 is constant this is inaccurate.
Bloerb is offline   Reply With Quote

Old   October 10, 2019, 06:44
Default Solved!
  #6
New Member
 
Join Date: Sep 2019
Posts: 8
Rep Power: 6
foamyva is on a distinguished road
Thank you very much! I now integrated the inlet and outlet field using patchMassFlowAverage
Code:
 massFlowAverageT
    {
    type patchMassFlowAverage;
    functionObjectLibs ( "libsimpleFunctionObjects.so" );
    fields ( T );
    patches
    (
        inlet
        outlet
    );
    factor 1.0;
    region coilWater;
    outputControlMode outputTime;
    outputInterval 1;
    writeStartTime no;
    verbose true;
    }
In this way you do account for the varying velocity field over the outflow patch together with temperature and indeed the wall heat flux and flow rate based Qdot match.

Naturally, if one wants to be really precise the variation of Cp over the patch should also be taken into account as you showed in the derivation. However, taking into account the minimal variation of Cp (and the hassle to get obtain that variable in post-processing) I left that out for now.

Again thanks a lot for your detailed elaboration!
foamyva is offline   Reply With Quote

Reply

Tags
chtmultiregionfoam, heat exchange, wall heat transfer, wallheatflux

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
OF 6 wallHeatFlux utility not working on chtMultiRegionFoam tutorial troparry OpenFOAM Post-Processing 1 January 10, 2022 07:12
Error in thermophysical properties (chtMultiRegionFoam) mukut OpenFOAM Pre-Processing 28 November 23, 2021 07:34
wallHeatFlux utility with sonicFoam hiuluom OpenFOAM Post-Processing 3 September 18, 2015 08:39
wallHeatFlux utility with chtMultiRegionFoam solver. Usem OpenFOAM Programming & Development 9 January 6, 2014 06:10
wallHeatFlux utility and chtMultiRegionFoam solver Lada OpenFOAM Post-Processing 4 June 7, 2012 10:46


All times are GMT -4. The time now is 00:12.