View Issue Details

IDProjectCategoryView StatusLast Update
0001971OpenFOAMBugpublic2016-01-17 11:04
Reporterrlange Assigned Tohenry  
PrioritynormalSeveritycrashReproducibilityalways
Status resolvedResolutionfixed 
PlatformGNU/LinuxOSUbuntuOS Version14.04
Summary0001971: Possible bug in the use of transonic in twoPhaseEulerFoam
DescriptionWhen transonic is enabled in twoPhaseEulerFoam,
It showed the error message for all cases tested and tutorials.

--> FOAM FATAL ERROR:
cannot scale a matrix containing a faceFluxCorrection
      
    From function fvMatrix<Type>::operator*=(const DimensionedField<scalar, volMesh>&)
          
in file src/finiteVolume/lnInclude/fvMatrix.C, at line 1220.
FOAM aborting.
I don't know, but this looks like a bug to me.
TagsNo tags attached.

Activities

henry

2016-01-14 10:28

manager   ~0005830

Which version of OpenFOAM are you running and have you also tested reactingTwoPhaseEulerFoam which will replace twoPhaseEulerFoam after it is fully validated.

rlange

2016-01-14 10:33

reporter   ~0005831

I tested versions 2.3.1 and 3.0.1.
I didn't test reactingTwoPhaseEulerFoam. I'll test!

henry

2016-01-14 11:30

manager   ~0005832

I have some thoughts on this. What we need is a suitable test-case (small-ish, 2D) which is transonic with a modest max Mach number; do you have something like this setup? What twoPhaseEulerFoam cases are you running which need transonic support in the pressure equation?

rlange

2016-01-14 17:19

reporter   ~0005833

Last edited: 2016-01-14 17:24

I am simulating a tube rupture of a industrial heat exchanger, the tube pressure is 88 bar and shell side has a pressure of 1 bar. The tube has natural gas and shell has thermal oil.

About the test-case 2D...
I guess so, before, i had tested a case 2D with VOF.
Mesh made by blockMesh, 90000 cells, setup with the above described values, basically.
Within a certain time of simulation, the Mach number was less than 0.5.

henry

2016-01-14 17:32

manager   ~0005834

If you can provide a suitable 2D test-case, ideally smaller than 90000 cells I will investigate this further.

henry

2016-01-17 11:04

manager   ~0005838

Resolved in OpenFOAM-dev by commit 7f67b3e46e6b2114db8ba7c5b863ae509797b430

Issue History

Date Modified Username Field Change
2016-01-14 10:16 rlange New Issue
2016-01-14 10:28 henry Note Added: 0005830
2016-01-14 10:33 rlange Note Added: 0005831
2016-01-14 11:30 henry Note Added: 0005832
2016-01-14 17:19 rlange Note Added: 0005833
2016-01-14 17:24 rlange Note Edited: 0005833
2016-01-14 17:32 henry Note Added: 0005834
2016-01-17 11:04 henry Note Added: 0005838
2016-01-17 11:04 henry Status new => resolved
2016-01-17 11:04 henry Resolution open => fixed
2016-01-17 11:04 henry Assigned To => henry