View Issue Details

IDProjectCategoryView StatusLast Update
0003319OpenFOAMBugpublic2020-02-05 20:14
Reporteryfclark Assigned Tohenry  
PrioritynormalSeverityminorReproducibilityalways
Status closedResolutionno change required 
PlatformGNU/LinuxOSUbuntuOS Version15.04
Summary0003319: Bug about surfaceTension implement in reactingmultiphaseEulerFoam
DescriptionHello.I am studying the reactingmultiphaseEulerFoam,
the code for surfaceTension implement is:
https://github.com/OpenFOAM/OpenFOAM-5.x/blob/7f7d351b741bf6406366a043cac98de56d2d44dd/applications/solvers/multiphase/reactingEulerFoam/reactingMultiphaseEulerFoam/pU/pEqn.H#L137

forAll(phases, phasei)
{
phaseModel& phase = phases[phasei];

    phigs.set
    (
        phasei,
        (
            alpharAUfs[phasei]
           *(
               ghSnGradRho
             - (fvc::interpolate(phase.rho() - rho))*(g & mesh.Sf())
             - fluid.surfaceTension(phase)*mesh.magSf()
            )
        ).ptr()
    );
}
finnaly,the velocity is updated by:

           phase.U() =
                HbyAs[phasei]
              + fvc::reconstruct
                (
                    alpharAUfs[phasei]*mSfGradp
                  - phiFs[phasei]
                  - phigs[phasei]
                );
It seems that the code updated velocity as:
https://user-images.githubusercontent.com/17508779/61933300-29b3c000-afb8-11e9-9f0c-f60759476f3c.png
My derivation mathematical formula is:
https://user-images.githubusercontent.com/17508779/61933381-536ce700-afb8-11e9-9255-1931029aca27.png

TagsNo tags attached.

Activities

yfclark

2019-07-26 08:58

reporter  

openfoam5.0.gif (1,931 bytes)   
openfoam5.0.gif (1,931 bytes)   
my.gif (1,892 bytes)   
my.gif (1,892 bytes)   

henry

2019-07-26 09:07

manager   ~0010648

Last edited: 2019-07-26 09:21

The formulation you have looks inconsistent with that used in VoF mixture velocity, if you combine two phase velocity expressions to create the corresponding mixture velocity equation what is the resulting net surface tension force?

Do you have any test cases which demonstrate that your formulation is better than the current?

yfclark

2019-07-26 09:17

reporter   ~0010649

thanks for you thanks for you suggection,I will do some test to make sure which is right!

henry

2020-02-05 20:14

manager   ~0011120

Orphaned report.

Issue History

Date Modified Username Field Change
2019-07-26 08:58 yfclark New Issue
2019-07-26 08:58 yfclark File Added: openfoam5.0.gif
2019-07-26 08:58 yfclark File Added: my.gif
2019-07-26 09:07 henry Note Added: 0010648
2019-07-26 09:17 yfclark Note Added: 0010649
2019-07-26 09:21 henry Note Edited: 0010648
2020-02-05 20:14 henry Assigned To => henry
2020-02-05 20:14 henry Status new => closed
2020-02-05 20:14 henry Resolution open => no change required
2020-02-05 20:14 henry Note Added: 0011120