13 rhoThermo& thermo1 =
phase1.thermo();
14 rhoThermo& thermo2 =
phase2.thermo();
surfaceScalarField & phi1
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const volScalarField & alpha2
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField