|
|
|
|
|
|
|
|
t0.png (19,395 bytes) |
|
|
t0.5.png (19,437 bytes) |
|
|
t0.5_patch.png (19,410 bytes) |
|
|
How recent is the OpenFOAM-2.3.x you are using? |
|
|
10th Feb. Commit b9587e0ebc6eeec1e04f1bca4e8c301460905f15. |
|
|
What happens if you run without correctPhi: PIMPLE { momentumPredictor no; nOuterCorrectors 3; nCorrectors 1; nNonOrthogonalCorrectors 0; correctPhi no; } Note that correctPhi is a predictor step, it is not exact and not always beneficial. For these kinds of cases you should run with nOuterCorrectors > 1 to ensure all transport is updated following the changes in the fluxes and commpression effect, particularly for morphing meshes. |
|
|
Running with "correctPhi no" removes the problem, which I guess isn't a surprise. Without correctPhi, fvc::absolute(phi) and (Uf & mesh.Sf()) should be identical. Running with more outer correctors helps, too. I guess this isn't a bug then? Conclusion is that correctPhi shouldn't be used without outer correctors. Thanks for the help. |
|
|
correctPhi may be needed to ensure boundedness of the phase-fraction in the predictor step. Outer-correctors are generally needed for cases with mesh-motion other than solid-body motion. |
|
|
I have studied the changes to divU in your case and believe there is a better way to correctPhi while preserving the coppressibility effect. I will push the change when I have finished testing. |
|
|
Resolved by commit c1417d4c56fbceb4975120fef9f7e74afef22fa1 |
- Anonymous
