|
[Sponsors] |
Is the mixerVessel2D tutorial from reactingMultiphaseEulerFoam actually working? |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Mazen Draw
Join Date: Sep 2018
Posts: 21
Rep Power: 7 ![]() |
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? |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 931
Rep Power: 12 ![]() |
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. "
__________________
The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki-1/Wiki-2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESI-OpenCFD-Trademark) Bugs/Features: FOAM-Extend (Wikki-FSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror |
|
![]() |
![]() |
![]() |
![]() |
#3 | |
New Member
Mazen Draw
Join Date: Sep 2018
Posts: 21
Rep Power: 7 ![]() |
Quote:
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. |
||
![]() |
![]() |
![]() |
![]() |
#4 |
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 931
Rep Power: 12 ![]() |
no idea, sorry. might consider to ask the Foundation by using one of the links below.
__________________
The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki-1/Wiki-2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESI-OpenCFD-Trademark) Bugs/Features: FOAM-Extend (Wikki-FSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror |
|
![]() |
![]() |
![]() |
![]() |
#5 |
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 931
Rep Power: 12 ![]() |
- Could you please elaborate, if possible, why the results look wrong to you? Is there any benchmark data you do compare with?
__________________
The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki-1/Wiki-2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESI-OpenCFD-Trademark) Bugs/Features: FOAM-Extend (Wikki-FSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror |
|
![]() |
![]() |
![]() |
![]() |
#6 | |
New Member
Mazen Draw
Join Date: Sep 2018
Posts: 21
Rep Power: 7 ![]() |
Quote:
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. |
||
![]() |
![]() |
![]() |
![]() |
#7 | |
New Member
Mazen Draw
Join Date: Sep 2018
Posts: 21
Rep Power: 7 ![]() |
Quote:
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. |
||
![]() |
![]() |
![]() |
Tags |
bug?, mixervessel2d, reactingmultiphaseeulerf |
|
|
![]() |
||||
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 |