View Issue Details
ID | Project | Category | View Status | Date Submitted | Last Update |
---|---|---|---|---|---|
0003319 | OpenFOAM | Bug | public | 2019-07-26 08:58 | 2020-02-05 20:14 |
Reporter | yfclark | Assigned To | henry | ||
Priority | normal | Severity | minor | Reproducibility | always |
Status | closed | Resolution | no change required | ||
Platform | GNU/Linux | OS | Ubuntu | OS Version | 15.04 |
Summary | 0003319: Bug about surfaceTension implement in reactingmultiphaseEulerFoam | ||||
Description | Hello.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 | ||||
Tags | No tags attached. | ||||
|
|
|
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? |
|
thanks for you thanks for you suggection,I will do some test to make sure which is right! |
|
Orphaned report. |
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 |