|July 1, 2015, 07:35||
Discretisation of the k-epsilon model for neutral athmospheric flows for OpenFOAM
Join Date: Feb 2015
Posts: 13Rep Power: 2
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]));
*(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);
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.
|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|