CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

Discretisation of the k-epsilon model for neutral athmospheric flows for OpenFOAM

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

LinkBack Thread Tools Display Modes
Old   July 1, 2015, 07:35
Default Discretisation of the k-epsilon model for neutral athmospheric flows for OpenFOAM
New Member
Join Date: Feb 2015
Posts: 11
Rep Power: 2
allett02015 is on a distinguished road
Hello everebody out there,

I'm corrently struggelingwith the k-esplion model applied to neutral athmosperic boundary layers.

In order to check the numerical solution I set up a Couette flow in OpenfFOAM (which matches the zero Pressure gradient assumption) with periodic boundary conditions.

At the walls I applied
nutkAtmRoughWallFunction, the no slip for the velocity, zeroGradient for k
and for espilon a fixeed value of

epsilon[cellI] += w*Cmu75*pow(k[cellI], 1.5)/(kappa_*(y[faceI] + z0_[faceI]));
G[cellI] +=
*(nutw[faceI] + nuw[faceI])
/(kappa_*(y[faceI] + z0_[faceI]));


SigmaEsp I set equal to kappa^2/(Cm^0.5 * (C2 - C1)). I thought that the Log-profile would result for the velocity. But I realized that the profiel was too steep. The reason is that the inclination turbulent length scale was to low.

The problem probably is due to a non neglectable discretization error. See

Richards, P. J., and S. E. Norris. "Appropriate boundary conditions for computational wind engineering models revisited." Journal of Wind Engineering and Industrial Aerodynamics 99.4 (2011): 257-266.


k-epsilon simulations of the neutral ABL: Achieving horizontal homogeneity on practical grid, Sumner and Masson, 2010.

Did anyone solve this issue.

I tried to implement the solution of Richard and Norris (averaging the turbulent shear stress over the cell faces) in OpenFoam. Unfortunatly it didnt work.

Here the relevant piece of the sourcecode:

surfaceVectorField GradU (fvc::snGrad(U_));

dimensionedScalar dims ( "dims", dimensionSet(0,0,-1,0,0,0,0), scalar(1.0));

surfaceVectorField GradUsqr (fvc::snGrad(U_)*dims);

forAll(GradUsqr, i)
GradUsqr[i].component(0) = GradU[i].component(0)*GradU[i].component(0);
GradUsqr[i].component(1) = GradU[i].component(1)*GradU[i].component(1);
GradUsqr[i].component(2) = GradU[i].component(2)*GradU[i].component(2);
// GradUsqr = sqr (GradU);
volVectorField snGradUSqr("snGradUSqr",fvc::average( GradUsqr ) );
volScalarField G(GName(), nut_*( mag ( snGradUSqr ) ) );

Maybe someone has a suggestion how to solve this issue.
allett02015 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
Cavitation model in compressible flows. jinwon park Main CFD Forum 1 May 11, 2015 07:24
Calculation of Epsilon in LES model Yuri Almeida OpenFOAM 2 January 18, 2013 08:43
k and epsilon were hard to converge in multiphase model of Fluent Yanlong Li ANSYS 0 January 2, 2013 06:25
K epsilon model michael OpenFOAM Running, Solving & CFD 1 August 24, 2006 02:57
How about the S-A model for blunt boby flows? maximus Main CFD Forum 0 March 8, 2005 22:47

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