CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Implementing kEpsilon turbulance mode with the airfoil example.....

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 8, 2014, 01:18
Default Implementing kEpsilon turbulance mode with the airfoil example.....
  #1
Senior Member
 
tareqkh's Avatar
 
Khamlaj
Join Date: Nov 2010
Location: United States
Posts: 137
Rep Power: 6
tareqkh is on a distinguished road
Hey Foamers,

I have implemented kOmega turbulence model in the airfoil2d example. However, when I switched to kepsilon turbulence model following all required steps, simpleFoam was unable to run. In the 0 folder, I put k and epsilon folders, and I deleted the file nuTilda and that is based on what I have read on the CFD forum.

My k is:

Quote:
dimensions [ 0 2 -2 0 0 0 0 ];

internalField uniform 0.01;

boundaryField
{
inlet
{
type fixedValue;
value uniform 0.01;
}

outlet
{
type zeroGradient;
}

walls
{
type kqRWallFunction;
value uniform 0.01;
}

frontAnBack
{
type empty;
}
}
My epsilon file is

Quote:
dimensions [ 0 2 -3 0 0 0 0 ];

internalField uniform 4.85;

boundaryField
{
inlet
{
type fixedValue;
value $internalField;
}

outlet
{
type zeroGradient;
}

walls
{
type epsilonWallFunction;
value $internalField;
}

frontAndBack
{
type empty;
}
}
fscheme is

Quote:
ddtSchemes
{
default steadyState;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss upwind;
div(phi,k) Gauss upwind;
div(phi,epsilon) Gauss upwind;
div(phi,R) Gauss upwind;
div(R) Gauss linear;
div((nuEff*dev(grad(U).T()))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DepsilonEff,epsilon) Gauss linear corrected;
laplacian(DREff,R) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
interpolate(U) linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}
fsolution is

Quote:
solvers
{
p GAMG
{
tolerance 1e-16;
relTol 0.01;

smoother GaussSeidel;

cacheAgglomeration true;

nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;

minIter 0;
maxIter 1000;
nPreSweeps 0;
nPostSweeps 2;
nFinestSweeps 2;
scaleCorrection true;
directSolveCoarsest false;
};

Urel smoothSolver
{
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-16;
relTol 0.01;
maxIter 500;
};

k smoothSolver
{
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-16;
relTol 0.01;
maxIter 500;
};

epsilon smoothSolver
{
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-16;
relTol 0.01;
maxIter 500;
};

omega smoothSolver
{
smoother GaussSeidel;
nSweeps 1;
tolerance 1e-16;
relTol 0.01;
maxIter 500;
};
}

SIMPLE
{
nNonOrthogonalCorrectors 0;
}

relaxationFactors
{
fields
{
p 0.3;
}
equations
{
Urel 0.6;
k 0.5;
epsilon 0.5;
omega 0.5;
R 0.5;
nuTilda 0.5;
}
}
Can somebody explain why it is not working?

Thank you....
tareqkh is offline   Reply With Quote

Reply

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
Problems with the new PyFoam release zfaraday OpenFOAM 13 December 9, 2014 19:58
pyfoam vitorspadetoventurin OpenFOAM Running, Solving & CFD 3 December 2, 2014 08:18
2D FFD Optimization RLangtry SU2 2 August 5, 2014 09:48
SU2 AOA optimization 454514566@qq.com SU2 8 July 7, 2014 07:01
Problem with airfoil shape optimization robyTKD SU2 Shape Design 6 June 13, 2013 19:13


All times are GMT -4. The time now is 02:21.