[saicharan662000@gmail.com]

Hello guys,
I got interface area code from some resource. I added It in twoPhaseMixtureThermo.C and I declared the Interface area in twoPhaseMixtureThermo.H . These files are in compressible interfoam solver folder . can anyone tell me How to print the interface area and How to make it visible in paraFoam. The solver compiled well but Interface area is not getting displayed.



I am attaching the code here
Foam::tmp<Foam::volScalarField>
Foam::twoPhaseMixtureThermo::interfaceArea() const
{
// return the interfacial area based on model for interfacial area
// returns dimensions Area

// model based on regular volume cells, taking the largest cut area
// as maximum for area, linear increase and decrease with alpha
const volScalarField& cellVolume =
alpha1().db().lookupObject<volScalarField>("cellVo lu");
volScalarField limitedAlpha1 = min(max(alpha1(), scalar(0)), scalar(1));

const dimensionedScalar
areaFactor("areaFactor",dimensionSet(0,2,0,0,0,0,0 ), 0.0);

volScalarField interfaceArea = alpha1() * areaFactor;
volScalarField maxArea = alpha1() * areaFactor;

maxArea = sqrt(3.0)*pow(cellVolume,(2.0/3.0));
return tmp<volScalarField>
(
(neg(limitedAlpha1-0.5)*maxArea*2.0*limitedAlpha1) +
(pos(limitedAlpha1-0.5)*maxArea*(-2.0*( limitedAlpha1 - 1.0)))
);

}
It's urgent. can anyone help me with this.

[saicharan662000@gmail.com]

Problem solved.