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/)
-   -   Turning interDymFoam into a turbulence solver (http://www.cfd-online.com/Forums/openfoam-solving/77644-turning-interdymfoam-into-turbulence-solver.html)

Saśl Balsa June 29, 2010 14:06

Turning interDymFoam into a turbulence solver
 
I am working with the interDymSolver, working with the sloshing2Dtank tutorial suited to my specific problem. I already performed a simulation in laminar regime and everything worked OK. Now I need to perform the same simulation with a kEpsilon turbulence model. I found some threads on the forum providing same guidance on it, namely, H.Weller stated the next,

"interDyMFoam is setup for laminar flow. If you want to use it for turbulent flow simply transfer the turbulence model parts from either rasInterFoam for RAS or lesInterFoam for LES. "

I copied all the files on the rasInterFoam model (dambreak) and adapted them all to my case. When I started running the case the next error message came up:


MacBook-Pro-de-Saul:BMark_Roll_Keps saul$ interDyMFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.5 |
| \\ / A nd | Web: http://www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Exec : interDyMFoam
Date : Jun 29 2010
Time : 19:46:53
Host : MacBook-Pro-de-Saul.local
PID : 663
Case : /Users/saul/OpenFOAM/saul-1.5/run/tutorials/interDyMFoam/BMark_Roll_Keps
nProcs : 1

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Selecting dynamicFvMesh solidBodyMotionFvMesh
Selecting solid-body motion function SDA

Reading environmentalProperties
Reading field pd

Reading field gamma

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting RAS turbulence model kEpsilon


keyword pRefProbe is undefined in dictionary "/Users/saul/OpenFOAM/saul-1.5/run/tutorials/interDyMFoam/BMark_Roll_Keps/system/fvSolution::PISO"

file: /Users/saul/OpenFOAM/saul-1.5/run/tutorials/interDyMFoam/BMark_Roll_Keps/system/fvSolution::PISO from line 70 to line 75.

From function dictionary::subDict(const word& keyword) const
in file db/dictionary/dictionary.C at line 271.

FOAM exiting

I have lack of experience with the OpenFOAM programming language so after checking the fvSolution I had no many other tricks to use ... Can anyone tell me where the mistake could be?

Ths!
Saśl

anmartin June 29, 2010 14:18

Hello,

Probably you need to define it, something like

PISO
{
momentumPredictor no;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
nAlphaCorr 1;
nAlphaSubCycles 4;
cAlpha 2;
correctPhi no;

pRefPoint (0 0.02 0.251);
pRefValue 1e5;
}

best

Saśl Balsa June 29, 2010 15:24

Thx Angel,

What I've done so far on this regard is the next,


PISO
{
momentumPredictor no;
nCorrectors 3;
nNonOrthogonalCorrectors 0;
nGammaCorr 1;
nGammaSubCycles 4;
cGamma 2;
}

...this is what is written on the fvSolution dictionary. I did not add your last two sentences.

pRefPoint (0 0.02 0.251);
pRefValue 1e5;

...before adding them onto the dictionary, could you please give me an idea of their purpose?

Thx.
Best
Saśl

Saśl Balsa June 30, 2010 16:09

I found that those lines are for the PISO algorithm to set up a reference pressure for calculations ... so , bearing that in mind, I wrote those two lines adapting to my case but it still does not work.

The OpenFOAM complains about the same error ... not pRefProbe defined on controlDict dictionary ...

Any idea on where is mistake? What I did do far is to copy all the files from the rasInterFoam to the interDymFoam, but suited to my case. Then I switched the model to turbulence (kEpsilon) on the RASProperties dict.

Dont know wether is too much left to do. Any further guidance I'd need to know???

Thx
Saśl


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