View Issue Details

IDProjectCategoryView StatusLast Update
0003319OpenFOAMBugpublic2019-07-30 14:25
ReporteryfclarkAssigned To 
PrioritynormalSeverityminorReproducibilityalways
Status newResolutionopen 
PlatformGNU/LinuxOSUbuntuOS Version15.04
Product Version6 
Fixed in Version 
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

View 2 revisions

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!

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 View Revisions