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

Is the mixerVessel2D tutorial from reactingMultiphaseEulerFoam actually working?

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

Like Tree1Likes
  • 1 Post By MazenDraw

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 7, 2020, 07:23
Default Is the mixerVessel2D tutorial from reactingMultiphaseEulerFoam actually working?
  #1
New Member
 
Mazen Draw
Join Date: Sep 2018
Posts: 21
Rep Power: 7
MazenDraw is on a distinguished road
I tried running the mixerVessel2D tutorial from both OpenFOAM Foundation v 7 and OpenFOAM v 1912.


On OF 7 the tutorial finishes but the results look wrong (see attachments).


On OF 1912 the tutorial crashes at T = 3.8 s, and the results of the older time steps look also wrong.


Can someone confirm this to me please?
Attached Images
File Type: png pressure.png (58.0 KB, 25 views)
File Type: jpg Uairmax.jpg (87.9 KB, 18 views)
File Type: jpg Uwatermax.jpg (96.6 KB, 19 views)
File Type: png alphaAir.png (136.0 KB, 18 views)
MazenDraw is offline   Reply With Quote

Old   May 7, 2020, 12:25
Default
  #2
HPE
Senior Member
 
HPE's Avatar
 
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 12
HPE is on a distinguished road
It blows up for v1912, it was discovered here: https://develop.openfoam.com/Develop...e_requests/347

"
multiphase/interFoam/RAS/DTCHullMoving with hierarchical decomposition throws Floating point exception for unmodified and modified cases.

multiphase/reactingMultiphaseEulerFoam/laminar/mixerVessel2D throws error for unmodified and modified cases.

multiphase/reactingTwoPhaseEulerFoam/RAS/LBend returns different results for separate runs even though I have not modified this tutorial at all.
"
HPE is offline   Reply With Quote

Old   May 7, 2020, 12:42
Default
  #3
New Member
 
Mazen Draw
Join Date: Sep 2018
Posts: 21
Rep Power: 7
MazenDraw is on a distinguished road
Quote:
Originally Posted by HPE View Post
It blows up for v1912, it was discovered here: https://develop.openfoam.com/Develop...e_requests/347

"
multiphase/interFoam/RAS/DTCHullMoving with hierarchical decomposition throws Floating point exception for unmodified and modified cases.

multiphase/reactingMultiphaseEulerFoam/laminar/mixerVessel2D throws error for unmodified and modified cases.

multiphase/reactingTwoPhaseEulerFoam/RAS/LBend returns different results for separate runs even though I have not modified this tutorial at all.
"

Do you also happen to have OF v7? Because there it is not crashing but IMO the results look wrong. Can you also confirm that?
I wanna know if there is an issue with MRF with reactingMultiphaseEulerFoam in general in those releases, because I am trying to create a case setup with MRF in reactingMultiphaseEulerFoam but it's not working. However, the same setup is working with reactingTwoPhaseEulerFoam.
MazenDraw is offline   Reply With Quote

Old   May 7, 2020, 13:17
Default
  #4
HPE
Senior Member
 
HPE's Avatar
 
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 12
HPE is on a distinguished road
no idea, sorry. might consider to ask the Foundation by using one of the links below.
HPE is offline   Reply With Quote

Old   May 8, 2020, 07:16
Default
  #5
HPE
Senior Member
 
HPE's Avatar
 
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 12
HPE is on a distinguished road
- Could you please elaborate, if possible, why the results look wrong to you? Is there any benchmark data you do compare with?
HPE is offline   Reply With Quote

Old   May 8, 2020, 07:27
Default
  #6
New Member
 
Mazen Draw
Join Date: Sep 2018
Posts: 21
Rep Power: 7
MazenDraw is on a distinguished road
Quote:
Originally Posted by HPE View Post
- Could you please elaborate, if possible, why the results look wrong to you? Is there any benchmark data you do compare with?

The pressure field should be low in the rotating region, and high in the stationary region. This pressure gradient is created in stirred tank cases and it causes the flow back to the rotating region. This is not the case here. In this tutorial, you get pressure gradient which is not related to this phenomenon. You can check the results of the mixerVessel2D of reactingTwoPhaseEulerFoam to see what I am talking about.


Also the max velocity magnitude is much higher than it should be. It is 5.6 m/s for air velocity and 1.1 m/s for water velocity. However, the inner region with r = 0.05 m is rotating with speed of 10 rad/s. Therefore, the highest magnitude should be around 0.5 m/s.


Compare the results of mixerVessel2D of rTPEF with the results of mixerVessel2D of rMEF and tell me what you think.
MazenDraw is offline   Reply With Quote

Old   May 8, 2020, 07:45
Default
  #7
New Member
 
Mazen Draw
Join Date: Sep 2018
Posts: 21
Rep Power: 7
MazenDraw is on a distinguished road
Quote:
Originally Posted by MazenDraw View Post
The pressure field should be low in the rotating region, and high in the stationary region. This pressure gradient is created in stirred tank cases and it causes the flow back to the rotating region. This is not the case here. In this tutorial, you get pressure gradient which is not related to this phenomenon. You can check the results of the mixerVessel2D of reactingTwoPhaseEulerFoam to see what I am talking about.
.

I checked the pressure field in rTPEF case, and it looks constant everywhere. Now I am confused!


However, the velocity vector field looks good to me.


I think you are right. I need a validation case to tell whether results are physical or not.
Attached Images
File Type: png pressure.png (28.7 KB, 10 views)
File Type: jpg Uair.jpg (93.9 KB, 7 views)
File Type: jpg Uwater.jpg (94.1 KB, 6 views)
File Type: jpg alphaAir.jpg (43.4 KB, 8 views)
HPE likes this.
MazenDraw is offline   Reply With Quote

Reply

Tags
bug?, mixervessel2d, reactingmultiphaseeulerf

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
Wing external flow (tutorial) D-0679 OpenFOAM Running, Solving & CFD 3 September 16, 2017 15:49
Error during tutorial 2.1.5.3 Mapfields mazamitla OpenFOAM Running, Solving & CFD 3 November 9, 2016 13:12
[Tutorials] New Dakota - OpenFOAM Coupling Tutorial Tobi OpenFOAM Community Contributions 0 October 9, 2016 17:19
DPM parallel is not working but serial is working johnwinter FLUENT 1 March 27, 2012 02:01
Tutorial 17: Unable to run the simulation Mazen CFX 9 April 17, 2009 20:15


All times are GMT -4. The time now is 21:07.