View Issue Details

IDProjectCategoryView StatusLast Update
0002677OpenFOAMBugpublic2017-08-29 14:51
ReporterCarlen Assigned Tohenry  
PriorityhighSeverityminorReproducibilityhave not tried
Status resolvedResolutionfixed 
PlatformGNU/LinuxOSUbuntuOS Version14.04
Summary0002677: Bug in compressibleMultiphaseInterFoam
DescriptionIn pEqn.H just before updating flux phi in the forAllIter loop, I think there is a phasei++ missing. The phasei is initiated as 0 before the forAll loop but is not updated in this loop to update dgdt for each phase. Please see the code below.

            forAllIter
            (
                PtrDictionary<phaseModel>,
                mixture.phases(),
                phase
            )
            {
                phase().dgdt() =
                    pos(phase())
                  *(p_rghEqnComps[phasei] & p_rgh)/phase().thermo().rho();
                //missing phasei++; here
            }

From what I see here only first phase's dgdt is updated at each time step. Maybe I am wrong.
TagsNo tags attached.

Activities

henry

2017-08-29 14:51

manager   ~0008643

Thanks for finding and reporting this bug.

Resolved in OpenFOAM-5.x by commit adac0dbb1d7a816e7976ddc2cd01be23e1fbb50f

Resolved in OpenFOAM-dev by commit 92ef068675f205a00e7367f26d81c5c376575fcf

Issue History

Date Modified Username Field Change
2017-08-25 03:06 Carlen New Issue
2017-08-29 14:51 henry Assigned To => henry
2017-08-29 14:51 henry Status new => resolved
2017-08-29 14:51 henry Resolution open => fixed
2017-08-29 14:51 henry Note Added: 0008643