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

Turning interDymFoam into a turbulence solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 29, 2010, 14:06
Default Turning interDymFoam into a turbulence solver
  #1
New Member
 
Saśl Balsa
Join Date: Jun 2010
Posts: 5
Rep Power: 7
Saśl Balsa is on a distinguished road
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
Saśl Balsa is offline   Reply With Quote

Old   June 29, 2010, 14:18
Default
  #2
Member
 
angel
Join Date: May 2009
Location: Spain
Posts: 46
Rep Power: 8
anmartin is on a distinguished road
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
anmartin is offline   Reply With Quote

Old   June 29, 2010, 15:24
Default
  #3
New Member
 
Saśl Balsa
Join Date: Jun 2010
Posts: 5
Rep Power: 7
Saśl Balsa is on a distinguished road
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 is offline   Reply With Quote

Old   June 30, 2010, 16:09
Default
  #4
New Member
 
Saśl Balsa
Join Date: Jun 2010
Posts: 5
Rep Power: 7
Saśl Balsa is on a distinguished road
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
Saśl Balsa 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
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05
Code release: Flow Transition and Turbulence Chaoqun Liu Main CFD Forum 0 September 26, 2008 17:15
Natural convection - Inlet boundary condition max91 CFX 1 July 29, 2008 20:28
Turbulence boundary values lego CFX 9 October 25, 2002 11:55
turbulence modeling questions llowen Main CFD Forum 3 September 11, 1998 04:24


All times are GMT -4. The time now is 20:23.