CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

calculation of gradients near wall

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

LinkBack Thread Tools Display Modes
Old   November 5, 2009, 02:51
Default calculation of gradients near wall
New Member
Sumeet Kumar
Join Date: Mar 2009
Posts: 21
Rep Power: 10
crazysumi is on a distinguished road
Hi FOAMers,

I want to calculate gradients at the cell centroids connected to a wall patch.

The quantities I need to evaluate are :-
a) del Utau / del tau
b) del Unormal / del normal
c) del p / del tau

Utau is tangential component of velocity
Unormal is normal component of velocity
tau is a coordinate tangent to local wall face
normal is a coordinate normal to local wall face

Can someone assist me calculation of these quantities at the cell centroids. Hints and suggestions in this regard are greatly appreciated

Thanks & Regards
crazysumi is offline   Reply With Quote

Old   November 5, 2009, 09:54
New Member
Sumeet Kumar
Join Date: Mar 2009
Posts: 21
Rep Power: 10
crazysumi is on a distinguished road
Hi Again

I tried to calculate following quanitity..

1) del p / del tau
vectorField n = patch().nf() / mag( patch().nf() ) ;
volVectorField gradP(fvc::grad(db().lookupObject<volScalarField>( "p")));
const vectorField gradPInternal = gradP.boundaryField()[patchI].patchInternalField();
const scalarField gradPtau = mag(gradPInternal - n * (gradPInternal & n));

The scalarField gradPtau will give me the required derivative

However the above method cannot be reiterated for calculation of following derivatives.-
a) del U(tau) / del tau
b) del U(normal) / del normal

because the U (tau) and U (normal) vectors are only defined at the cell centroids of patchInternalField list. How do I calculate such derivatives. Please someone advice.

crazysumi 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
Fluent3DMeshToFoam simvun OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 48 May 14, 2012 05:20
Errors running allwmake in OpenFOAM141dev with WM_COMPILE_OPTION%3ddebug unoder OpenFOAM Installation 11 January 30, 2008 21:30
Warning 097- AB Siemens 6 November 15, 2004 05:41
Multicomponent fluid Andrea CFX 2 October 11, 2004 05:12
Wall function in adverse pressure gradients stephane baralon Main CFD Forum 11 September 2, 1999 04:05

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