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

Calculating heat flux

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

Like Tree6Likes
  • 1 Post By hugo17
  • 3 Post By hugo17
  • 2 Post By Nicole

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 22, 2011, 05:17
Default Calculating heat flux
  #1
New Member
 
Join Date: Jul 2011
Posts: 18
Rep Power: 14
hugo17 is on a distinguished road
Hi,

I'm dealing with heat conduction in solids using laplacianFoam and want to calculate the heat flux through walls from Openfoam results. However, unexpected difficulties are arising which I want to explain with the following example:



The analytical result for the heat flux through the walls of the annulus shown in the figure is q"=2460.01 W/m (related to length in z direction).

1. Using wallHeatFluxLaminar (from this thread: http://www.cfd-online.com/Forums/ope...foam1-6-a.html), the following values are calculated:
Outside: q"_a=780.284 W/m^2 --> multiplication by inner circumference gives q"=2451.3 W/m
Inside: q"_i=-1299.32 W/m^2 --> q"=-2449.2 W/m
Compared to the analytical solution, these values are too small.

2. Using foamCalc magGrad T, the following values are calculated:
Outside: magGrad T=78.0284 K/m --> multiplication by k and the inner circumference gives q"=2451.3 W/m (same value as given by wallHeatFluxLaminar)
Inside: magGrad T=129.932 K/m --> q"=-2449.2 W/m
Thus, using magGrad T leads to the same heat flux as wallHeatFluxLaminar for the present example. However, this is not the case for all geometries I dealt with.

3. Calculating the gradient magnitude for myself from mag(grad T) = sqrt(gradTx^2 + gradTy^2), I obtain different heat fluxes:
Outside: sqrt(gradTx^2 + gradTy^2)=78.307 K/m --> q"=2460.09 W/m
Inside: sqrt(gradTx^2 + gradTy^2) --> q"=2460.07 W/m
These heat fluxes are very close to the analytical value.

Now, my questions are:
- Why does the gradient magnitude calculated as sqrt(gradTx^2 + gradTy^2) differ from magGrad T calculated by foamCalc?
- Why is the heat flux value calculated from sqrt(gradTx^2 + gradTy^2) apparantly "better" than this one calculated from WallHeatFluxLaminar or magGrad T?

I really appreciate any comment on this problem.

Regards
Hugo
JuanRodriguez likes this.
hugo17 is offline   Reply With Quote

Old   August 22, 2011, 06:18
Default
  #2
New Member
 
Join Date: Jul 2011
Posts: 18
Rep Power: 14
hugo17 is on a distinguished road
I found the reason for the different results: I had writePrecision set to 6. When I increase it to 8, also wallHeatFluxLaminar and magGrad T give the correct values.

Regards
Hugo
Mojtaba.a, ancipdp and Saleh Mey like this.
hugo17 is offline   Reply With Quote

Old   May 7, 2015, 19:51
Default
  #3
New Member
 
Juan David Rodriguez P
Join Date: Jan 2015
Location: Milano
Posts: 20
Rep Power: 11
JuanRodriguez is on a distinguished road
Hugo, thanks for your contribution.

I want to visualize heat flux lines (orthogonal to isotherms) in paraview. My case is a similar one to yours (heat conduction in 2D from a buried pipe to soil). Do know how to proceed?

Best,
Juan.
JuanRodriguez is offline   Reply With Quote

Old   May 27, 2015, 07:09
Default
  #4
Member
 
Nicole Andrew
Join Date: Sep 2014
Location: Pretoria, South Africa
Posts: 58
Rep Power: 11
Nicole is on a distinguished road
Hi Juan,

Have you had any success yet?

If you are used to modifying the solver, I would recommend just adding a new volScalarField, qFlux or something, calculated from fvc::grad(T) in the solver. I am not sure if this would be the best way.

Apparently there are also calculators in paraview, perhaps one of these would be useful?

Let me know if you want me to elaborate
prandtl Gao and gbope7 like this.
Nicole is offline   Reply With Quote

Old   May 27, 2015, 13:47
Default
  #5
New Member
 
Juan David Rodriguez P
Join Date: Jan 2015
Location: Milano
Posts: 20
Rep Power: 11
JuanRodriguez is on a distinguished road
Hi Nicole,

I finally used funkySetFields to calculate grad(T). I think it is ok and am still learning to modify code in OF, but I will try what you suggest.

Thanks.
Best,
Juan.
JuanRodriguez is offline   Reply With Quote

Reply

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
Sign of Heat Flux at wall Kyung FLUENT 2 February 26, 2016 17:25
Heat Flux at wall in a conjugate heat transfer problem Chander CFX 2 July 9, 2011 23:22
Heat Flux Wall Boundary Confusion. Joee FLUENT 1 August 21, 2010 13:20
Heat flux in ansys cfx juliom OpenFOAM Running, Solving & CFD 2 April 14, 2009 15:30
Calculating Heat flux Andrew Garrard FLUENT 0 July 6, 2004 06:47


All times are GMT -4. The time now is 21:18.