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

I can't add turbulence model in multiPhaseEulerFoam

Register Blogs Community New Posts Updated Threads Search

Like Tree1Likes
  • 1 Post By geth03

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 12, 2022, 11:07
Post I can't add turbulence model in multiPhaseEulerFoam
  #1
New Member
 
saeed sangchooly
Join Date: Feb 2022
Posts: 17
Rep Power: 4
saeed sangchooly is on a distinguished road
hi everyone

I'm trying to simulate a 3phase case with multiphaseEulerFoam solver. water, oil and air. with interfaceCompression for (air and water) and (air and oil).

I can't add komegaSSTSato or any other turbulence model. in the first time step when it comes to solving transport equation for turbulence parameters (k,omega,..) the solver give me the Floating point exception error.

I've tried using the results of a laminar case as initial condition for my turbulence case.
I've tried many hexahedral meshes with really good quality up to 1m cells.
I've tried many turbulence parameters estimations for initial and boundary condition for each phase.

I really dont know what I'm missing here.

it would be great if anyone can give me a clue .
saeed sangchooly is offline   Reply With Quote

Old   April 14, 2022, 21:34
Post
  #2
New Member
 
saeed sangchooly
Join Date: Feb 2022
Posts: 17
Rep Power: 4
saeed sangchooly is on a distinguished road
Quote:
Originally Posted by saeed sangchooly View Post
hi everyone

I'm trying to simulate a 3phase case with multiphaseEulerFoam solver. water, oil and air. with interfaceCompression for (air and water) and (air and oil).

I can't add komegaSSTSato or any other turbulence model. in the first time step when it comes to solving transport equation for turbulence parameters (k,omega,..) the solver give me the Floating point exception error.

I've tried using the results of a laminar case as initial condition for my turbulence case.
I've tried many hexahedral meshes with really good quality up to 1m cells.
I've tried many turbulence parameters estimations for initial and boundary condition for each phase.

I really dont know what I'm missing here.

it would be great if anyone can give me a clue .
can anyone help me?
saeed sangchooly is offline   Reply With Quote

Old   April 20, 2022, 05:18
Default
  #3
Senior Member
 
Join Date: Dec 2019
Location: Cologne, Germany
Posts: 369
Rep Power: 8
geth03 is on a distinguished road
can you tell which OF version you use?

can you use the banana-trick and show us the ouput?
should look like this:
--> FOAM FATAL ERROR:
Unknown RASModel type banana

Valid RASModel types:

2
(
kEpsilon
kOmegaSST
)

here you can see that multiphaseEulerFoam for my custom solver can only use 2 tubulence models. is komegaSSTSato a valid option in your case?
geth03 is offline   Reply With Quote

Old   April 24, 2022, 13:51
Post
  #4
New Member
 
saeed sangchooly
Join Date: Feb 2022
Posts: 17
Rep Power: 4
saeed sangchooly is on a distinguished road
Quote:
Originally Posted by geth03 View Post
can you tell which OF version you use?

can you use the banana-trick and show us the ouput?
should look like this:
--> FOAM FATAL ERROR:
Unknown RASModel type banana

Valid RASModel types:

2
(
kEpsilon
kOmegaSST
)

here you can see that multiphaseEulerFoam for my custom solver can only use 2 tubulence models. is komegaSSTSato a valid option in your case?
dear geth03
currently Im trying to use komegaSSTSato
and this is the output of banana-trick:

--> FOAM FATAL ERROR:
[4] Unknown RASModel type kOmegaSSSato

Valid RASModel types:

8
(
LaheyKEpsilon
continuousGasKEpsilon
kEpsilon
kOmegaSST
kOmegaSSTSato
kineticTheory
mixtureKEpsilon
phasePressure
)
saeed sangchooly is offline   Reply With Quote

Old   April 25, 2022, 08:27
Default
  #5
Senior Member
 
Join Date: Dec 2019
Location: Cologne, Germany
Posts: 369
Rep Power: 8
geth03 is on a distinguished road
ok, so kOmegaSSTSato is a valid option.

the floating point exception error ist most likely the result of a division by zero. when you initialize your turbulence values such as k and omega, do not set the values in the domain to zero.

for an initial test run do not use wall functions for turbulence properties, use zero gradient instead and see if it is running.

try these suggestions out and let me know if it works.
Ahyar likes this.
geth03 is offline   Reply With Quote

Old   May 15, 2022, 14:58
Post
  #6
New Member
 
saeed sangchooly
Join Date: Feb 2022
Posts: 17
Rep Power: 4
saeed sangchooly is on a distinguished road
Im sorry it took me a while to answer.

I've tried it (zerogradient) , it didn't work.
I've no idea...
saeed sangchooly is offline   Reply With Quote

Old   May 9, 2024, 14:38
Default
  #7
New Member
 
Raphael Santos
Join Date: Oct 2013
Posts: 20
Rep Power: 13
Raphael_Santos is on a distinguished road
Hi,

check if your setup has no "0" value for omega or k for each phase.

At least, use something as "value uniform 1e-14;", but avoid "0".
Raphael_Santos is offline   Reply With Quote

Reply

Tags
komegasstsato, multiphaseeulerfaom, openfoam, openfoam 1806, turbulence model


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
gamma-ReTheta turbulence model for predicting transitional flows FelixL OpenFOAM Programming & Development 123 August 30, 2022 12:50
NEW turbulence TRANSITIONAL model giammy92 OpenFOAM 3 June 30, 2016 10:47
Turbulence Model and limitation to Reynolds number qascapri FLUENT 0 January 24, 2011 11:48
Discussion: Reason of Turbulence!! Wen Long Main CFD Forum 3 May 15, 2009 10:52
Why Turbulence models are not universal. Senthil Main CFD Forum 4 July 5, 2000 05:34


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