[quarkz]

Hi,

I am trying to convert my original laminar 3D double cone case to a RAS case, with the k Epsilon model. The laminar case can run but it diverges after a short while. FYI, I'm using rhoCentralFoam with compressible high speed flow.

I tried to view some youtube tutorials on implementing a RAS case. However I got this error which I can't solve. So my simulation can't start.

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 6960 average: 6960
bounding epsilon, min: 0 max: 183 average: 183
kEpsilonCoeffs
{
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 -0.33;
sigmak 1;
sigmaEps 1.3;
}

No finite volume options present



--> FOAM FATAL ERROR:
Different dimensions for =
dimensions : [0 2 -1 0 0 0 0] = [1 -1 -1 0 0 0 0]


From function bool Foam::dimensionSet:perator=(const Foam::dimensionSet&) const
in file dimensionSet/dimensionSet.C at line 171.

FOAM aborting

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::error::abort() at ??:?
#2 Foam::dimensionSet:perator=(Foam::dimensionSet const&) const at ??:?
#3 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>:perator=(Foam::tmp<Foam::Geometri cField<double, Foam::fvPatchField, Foam::volMesh> > const&) in "/home/tsltaywb/OpenFOAM/OpenFOAM-4.1/platforms/linux64GccDPInt32Opt/bin/rhoCentralFoam"
#4 Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foa m::CompressibleTurbulenceModel<Foam::fluidThermo> > >::correctNut() at ??:?
#5 Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Fo am::ThermalDiffusivity<Foam::CompressibleTurbulenc eModel<Foam::fluidThermo> > > >::correctNut() at ??:?
#6 ? in "/home/tsltaywb/OpenFOAM/OpenFOAM-4.1/platforms/linux64GccDPInt32Opt/bin/rhoCentralFoam"
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? in "/home/tsltaywb/OpenFOAM/OpenFOAM-4.1/platforms/linux64GccDPInt32Opt/bin/rhoCentralFoam"
Aborted (core dumped)


So can someone tell me where the error can't from? So that I can look at the right place. Thanks!

[BSengupta]

Hi Can any one please help me to overcome it.. Than you very much.



Overall domain bounding box (0 0 0) (1 0.4 0.1)
p { margin-bottom: 0.1in; line-height: 115%; }




Error:
Time = 10


GAMG: Solving for Ux, Initial residual = 0.424802, Final residual = 0.017996, No Iterations 1
GAMG: Solving for Uy, Initial residual = 0.332272, Final residual = 0.0140346, No Iterations 1
GAMG: Solving for Uz, Initial residual = 0.559406, Final residual = 0.0285128, No Iterations 1
GAMG: Solving for e, Initial residual = 0.574963, Final residual = 0.0496275, No Iterations 1
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#4 Foam::hePsiThermo<Foam:siThermo, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#5 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7 ? in "/opt/openfoam6/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
Floating point exception (core dumped)


/*---------------------------------------------------------------------------------*/
Model: rhoSimpleFoam
/*---------------------------------------------------------------------------------*/
Boundary Condition:
1. Pressure
frontandback
{
type symmetry;


}
inlet
{
type zeroGradient;
//type mixed;
//refValue uniform 100000;
//refGradient uniform 0;
//valueFraction uniform 0.3;
}
outlet
{
type fixedValue;
value uniform 100000;
}
wall
{
type zeroGradient;
}


2. Velocity
frontandback
{
type symmetry;


}
inlet
{
type fixedValue;
value uniform (5 0 0);
rhoInlet 1.5;
}
outlet
{
type zeroGradient;
}
wall
{
type noSlip;
}


3. Temperature
frontandback
{
type symmetry;


}
inlet
{
type fixedValue;
value uniform 300;
}
outlet
{
type zeroGradient;
}
wall
{
type zeroGradient;
}


4. k
frontandback
{
type symmetry;


}
inlet
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.14;
value uniform 0.1;
}
outlet
{
type inletOutlet;
inletValue uniform 1;
value uniform 1;
}
wall
{
type kqRWallFunction;
value uniform 1;
}


5. epsilon
frontandback
{
type symmetry;


}
inlet
{
type turbulentMixingLengthDissipationRateInlet;
mixingLength 0.14;
value uniform 0.006;
}
outlet
{
type inletOutlet;
inletValue uniform 0.006;
value uniform 0.006;
}
wall
{
type epsilonWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0.006;
}


6. nut


frontandback
{
type symmetry;


}
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
wall
{
type nutkWallFunction;
Cmu 0.09;
kappa 0.41;
E 9.8;
value uniform 0;
}


7. alphat
frontandback
{
type symmetry;


}
inlet
{
type calculated;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
wall
{
type compressible::alphatWallFunction;
Prt 0.85;
value uniform 0;
}


/*------------------------------------------------------------------------------*/
fvSolution
ddtSchemes
{
default steadyState;
}


gradSchemes
{
default Gauss linear;
}


divSchemes
{
default none;


div(phi,U) bounded Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
div(phi,e) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;


div(phid,p) Gauss upwind;
div(phi,Ekp) bounded Gauss upwind;
div((phi|interpolate(rho)),p) Gauss upwind;
}


laplacianSchemes
{
default Gauss linear corrected;
}


interpolationSchemes
{
default linear;
}


snGradSchemes
{
default corrected;
}


/*----------------------------------------------------------------------*/
fvSolution
solvers
{
p
{
solver GAMG;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
nCellsInCoarsestLevel 20;
}


"(U|e|k|epsilon)"
{
solver GAMG;
tolerance 1e-08;
relTol 0.1;
smoother GaussSeidel;
nCellsInCoarsestLevel 20;
}
}


SIMPLE
{
nNonOrthogonalCorrectors 0;
rhoMin 0.1;
rhoMax 1.5;
//pMaxFactor 2;
//pMinFactor 0.1;
transonic yes;
consistent yes;


residualControl
{
p 1e-3;
U 1e-4;
e 1e-3;


// possibly check turbulence fields
"(k|epsilon|omega)" 1e-3;
}
}


relaxationFactors
{
fields
{
p 0.3;
rho 0.7;
}
equations
{
p 0.7;
U 0.7;
e 0.7;
k 0.7;
epsilon 0.7;
}
}






p { margin-bottom: 0.1in; line-height: 115%; }