CFD Online Discussion Forums

CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (http://www.cfd-online.com/Forums/openfoam-solving/)
-   -   using Spalart-Allmaras model in interFoam (http://www.cfd-online.com/Forums/openfoam-solving/122856-using-spalart-allmaras-model-interfoam.html)

vishal_s August 29, 2013 11:24

using Spalart-Allmaras model in interFoam
 
I an new to Openfoam , I am trying to use Spalart-Allmaras model for my nozzle flow in interFoam. I incoporated nuTilda and nutildaFinal in the fvsolutions of the system file, but I a getting this error. can any one tell me, how to solve this error.


Reading g
Calculating field g.h

time step continuity errors : sum local = 9.77995e-05, global = -9.77995e-05, cumulative = -9.77995e-05
DICPCG: Solving for pcorr, Initial residual = 1, Final residual = 6.94573e-11, No Iterations 194
time step continuity errors : sum local = 6.7998e-15, global = -4.64687e-17, cumulative = -9.77995e-05
Courant Number mean: 0.00278549 max: 0.285502

Starting time loop

Courant Number mean: 0.00193436 max: 0.198265
Interface Courant Number mean: 0 max: 0
deltaT = 0.000694444
Time = 0.000694444

MULES: Solving for alpha1
Phase-1 volume fraction = 0.0220219 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.0220388 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.0220558 Min(alpha1) = 0 Max(alpha1) = 1
MULES: Solving for alpha1
Phase-1 volume fraction = 0.0220728 Min(alpha1) = 0 Max(alpha1) = 1
DICPCG: Solving for p_rgh, Initial residual = 1, Final residual = 0.0419123, No Iterations 85
time step continuity errors : sum local = 4.60931e-06, global = 9.84271e-09, cumulative = -9.77897e-05
DICPCG: Solving for p_rgh, Initial residual = 0.000734712, Final residual = 3.64035e-05, No Iterations 65
time step continuity errors : sum local = 1.06533e-05, global = 2.14005e-06, cumulative = -9.56496e-05
DICPCG: Solving for p_rgh, Initial residual = 0.00326193, Final residual = 8.82474e-08, No Iterations 123
time step continuity errors : sum local = 4.04865e-09, global = 5.87426e-11, cumulative = -9.56496e-05
smoothSolver: Solving for nuTilda, Initial residual = 0, Final residual = 0, No Iterations 0
--> FOAM Warning :
From function tmp<volScalarField> SpalartAllmaras::k() const
in file SpalartAllmaras/SpalartAllmaras.C at line 243
Turbulence kinetic energy not defined for Spalart-Allmaras model. Returning zero field
--> FOAM Warning :
From function tmp<volScalarField> SpalartAllmaras::k() const
in file SpalartAllmaras/SpalartAllmaras.C at line 243
Turbulence kinetic energy not defined for Spalart-Allmaras model. Returning zero field
ExecutionTime = 0.18 s ClockTime = 1 s

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

This is the what I incorporated in the SYSTEM/FVSOLUTIONS

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.1.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
pcorr
{
solver PCG;
preconditioner DIC;
tolerance 1e-10;
relTol 0;
}

p_rgh
{
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0.05;
}

p_rghFinal
{
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0;
}



nuTildacorr
{
solver PCG;
preconditioner DIC;
tolerance 1e-10;
relTol 0;
}

nuTilda
{
solver PCG;
preconditioner DIC;
tolerance 1e-07;
relTol 0.05;
}

nuTildaFinal
{
solver smoothSolver;
preconditioner DIC;
tolerance 1e-07;
relTol 0;
}

}

PIMPLE
{
momentumPredictor no;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
nAlphaCorr 1;
nAlphaSubCycles 4;
cAlpha 2;
}


// ************************************************** *********************** //


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