CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

Wall correction for the temperature gradient

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

LinkBack Thread Tools Display Modes
Old   March 7, 2011, 21:04
Default Wall correction for the temperature gradient
New Member
Hagen Müller
Join Date: Nov 2010
Posts: 28
Rep Power: 7
Hagen is on a distinguished road
Dear all,

I would like to implement a correction for the temperature gradient near the wall in a 2D channel flow. The energy equation is usually formulated as

fvScalarMatrix TEqn
           fvm::div(phi, T)
           - fvm::Sp(fvc::div(phi), T)
           - fvm::laplacian(kappaEff, T)
In order to correct the temperature gradient in the laplacian term it seems to be necessary to rewrite it as laplacian(kappaEff, T) = div(kappaEff * grad(T)) and then manipulate grad(T) at the wall. I found the following formulation to replace the laplacian.

surfaceTensorField kappaEffs = fvc::interpolate(kappaEff);
      surfaceVectorField Sn = mesh.Sf()/mesh.magSf(); 
      surfaceVectorField SfkappaEff = mesh.Sf() & kappaEffs;
      GeometricField<scalar, fvsPatchField, surfaceMesh> SfkappaEffSn = SfkappaEff & Sn;

      volScalarField Dif1
                fvc::div(SfkappaEffSn * fvc::snGrad(T))

       volScalarField Dif2
                fvc::laplacian(kappaEff, T)
Dif 1 and Dif2 are almost identical and grad(T) could be manipulated at the wall. The problem is, that I can not use this formulation for the energy equation because it is in explicit form and an implicit divergence does not exist. Convergence schemes do not seem to work here because they do not accept surfaceFields.

I would be thankful for any suggestion how to express the laplacian in terms of the gradient or any other possibility to correct the temperature gradient at the wall.

Thank you.
Hagen is offline   Reply With Quote


Thread Tools
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 On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with zeroGradient wall BC for temperature - Total temperature loss cboss OpenFOAM 10 March 5, 2015 07:57
udf variable for wall temperature Kiran FLUENT 0 July 31, 2008 08:31
how to calculate the temperature gradient of wall houbaolin FLUENT 0 July 28, 2008 02:51
Non-conformal grid interfaces: wall temperature Evandro FLUENT 0 April 3, 2008 23:43
Variable temperature wall boundary condition Virginie FLUENT 2 October 2, 2007 09:23

All times are GMT -4. The time now is 23:16.