CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Error running RAS kEpsilon

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 13, 2017, 22:17
Default Error running RAS kEpsilon
  #1
Senior Member
 
TWB
Join Date: Mar 2009
Posts: 399
Rep Power: 17
quarkz is on a distinguished road
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!
quarkz is offline   Reply With Quote

Old   December 31, 2018, 04:10
Default
  #2
Member
 
Bidesh Sengupta
Join Date: Sep 2018
Location: Sngapore
Posts: 76
Rep Power: 7
BSengupta is on a distinguished road
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%; }
BSengupta is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Error in LRR (RAS) model in OpenFOAM v4.1 Inquisitor101 OpenFOAM 4 November 29, 2023 15:11
multiphaseInterFoam for RAS turbulence model chiven OpenFOAM Bugs 8 December 6, 2017 03:08
rhoSimpleFoam searches for RAS model hamsteri15 OpenFOAM Running, Solving & CFD 3 February 10, 2016 12:00
Explanation for terms like RAS turbulence, kEpsilon,... tH3f0rC3 OpenFOAM Pre-Processing 2 March 16, 2011 10:16
Turbulence dampening due to magnetic field in LES and RAS eelcovv OpenFOAM 0 June 8, 2010 12:35


All times are GMT -4. The time now is 10:55.