CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM

Compute local Reynolds number

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

LinkBack Thread Tools Display Modes
Old   July 12, 2010, 12:44
Default Compute local Reynolds number
Senior Member
Join Date: Mar 2009
Posts: 110
Rep Power: 10
nuovodna is on a distinguished road
Hi, i want to compute local Reynolds number (aka Peclet) using interFoam but i have some problems.
PecletX = rho * DeltaX * Ux / mu
PecletY = rho * DeltaY * Uy / mu

I create two fields : one for Peclet along X axis one along Y axis in this way

volScalarField PecletX
dimensionedScalar("PecletX", dimensionSet(1,0,-2,0,0),1)

volScalarField PecletY
dimensionedScalar("PecletY", dimensionSet(1,0,-2,0,0),1)

Now i have to compute DeltaX and DeltaY for all of cells

forAll(mesh.cells(), celli)
const labelList& cellPts = mesh.cellPoints()[celli];
DeltaX = mesh.points()[cellPts[1]].component(0) - mesh.points()[cellPts[0]].component(0);
DeltaY = mesh.points()[cellPts[4]].component(1) - mesh.points()[cellPts[0]].component(1);
PecletX = (U.component(0) * DeltaX * ;
PecletY = (U.component(1) * DeltaY * ;

The previous way is wrong. I have to make a cycle over whole points of all cells and storage the max and the min value of X (and Y). How can i do this??

nuovodna 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
Obtaining Reynolds Number in FLUENT Emmanuel FLUENT 2 April 7, 2016 01:46
Cell Reynolds Number laliong Main CFD Forum 11 March 19, 2015 12:42
How to write k and epsilon before the abnormal end xiuying OpenFOAM Running, Solving & CFD 8 August 27, 2013 15:33
Turbulence model for low Reynolds number flow? Nokadu Main CFD Forum 3 May 26, 2013 11:42
Reynolds number help Mech FLUENT 0 March 15, 2006 11:28

All times are GMT -4. The time now is 14:15.