interFoam algorithm
Hi :),
For each time step, I found that the order of the computations is a little bit strange. 1) It computes the new surface force term using the alpha field at the time "n" Code:
twoPhaseProperties.correct(); Code:
#include "alphaEqnSubCycle.H" Code:
//  Pressurevelocity PIMPLE corrector loop It seems more logical to update the surface force term (step 1) just after computing the new alpha field. In my mind, the algorithm should be : 1) update properties time "n" 2) compute implicitly velocitypressure field "n+1" 3) compute alpha field "n+1" using the velocitypressure previously computed Am I wrong ? Thank You, Pierre 
I confirm what I said before, this is an error in the "interFoam" solver (and the derivated solvers).
In the H. Rusche thesis, page 162, we can see the following solution procedure : steps 14 : refer to the moving frame step 5 : transport "alpha1" step 6 : update properties (smooth gamma => curvature computation => surface force term) step 7 and 8 : PISOloop to compute pressurevelocity fields. Maybe this error should be reported to the OpenFOAM Foundation to correct it for the following versions. Pierre 
How are step 58 different than what is already implemented?

In the current version of "interFoam", the step 6 is computed before the step 5.
First it computes the surface force term : Code:
twoPhaseProperties.correct(); Code:
#include "alphaEqnSubCycle.H" 
The interface curvature for the surface force term is calculated by calling interface.correct() after the alphaEqnSubCycle.

Ok, I had inverted the "correct" function of twoPhaseMixture for the viscosity with the "correct" function of interfaceProperties for the curvature.
Thank you for your reply 
All times are GMT 4. The time now is 15:40. 