Improvement of calculation in interface curvature in interFoam and InterIsoFoam
Functionality to add/problem to solve
This implementation reduces the calculation error of interface curvature in interFoam and interIsoFoam. This implementation was reported by Hoang et al. Comput. Fluids 86 (2013) 28-36. I implemented this method. This implementation reduces the spurious current drastically.
This problem is also reported in issue #2479 (closed).
Especially in the case of interIsoFoam, this filter drastically reduces the velocity of spurious current. (Yamamoto and Komarov, Investigation for optimal use of VOF method in OpenFOAM, Proc. of OpenCAE Symposium 2021 (2021) The OpenCAE Society of Japan)
Target audience
To who want to reduce the velocity of spurious current.
Proposal
I attach the repair patch for $FOAM_SRC/transportModels/interfaceProperties/interfaceProperties.C and $FOAM_SRC/transportModels/interfaceProperties/interfaceProperties.H interfacePropertiesC.patch
What does success look like, and how can we measure that?
I simulated the static bubble in water. In this case, the spurious current is generated due to the imbalance of surface tension force. We attached the both case with and without smoother.
With smoother
Without smoother
How to control the number of Laplacian filter
In system/fvSolution, one needs to add the following sentence (nSmooth) in addition to other parameters such as nAlphaCorr, nAlphaSubCycles.
nAlphaCorr 5; nAlphaSubCycles 1; nSmooth 3;
The larger number of nSmooth smooths the original distribution of alpha. Therefore, too large nSmooth has possibility to worsen the calculation accuracy.
Available solvers
When one uses this patch, the following solvers can use this smoother as an option.
- interFoam
- interIsoFoam
- overInterDyMFoam
- compressibleInterFoam
- compressibleInterDyMFoam
- compressibleInterIsoFoam
- overCompressibleInterDyMFoam
- compressibleInterFilmFoam
- interPhaseChangeFoam