# Wall correction for the temperature gradient

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 March 7, 2011, 21:04 Wall correction for the temperature gradient #1 Member   Hagen Müller Join Date: Nov 2010 Posts: 33 Rep Power: 8 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 Code: ```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. Code: ```surfaceTensorField kappaEffs = fvc::interpolate(kappaEff); surfaceVectorField Sn = mesh.Sf()/mesh.magSf(); surfaceVectorField SfkappaEff = mesh.Sf() & kappaEffs; GeometricField SfkappaEffSn = SfkappaEff & Sn; volScalarField Dif1 ( IOobject ( "Dif1", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), fvc::div(SfkappaEffSn * fvc::snGrad(T)) ); volScalarField Dif2 ( IOobject ( "Dif2", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE ), 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.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post cboss OpenFOAM 10 March 5, 2015 07:57 Kiran FLUENT 0 July 31, 2008 08:31 houbaolin FLUENT 0 July 28, 2008 02:51 Evandro FLUENT 0 April 3, 2008 23:43 Virginie FLUENT 2 October 2, 2007 09:23

All times are GMT -4. The time now is 22:55.

 Contact Us - CFD Online - Top