[romant]

I am not sure I understand. What do you mean with floder and which part of the code. Above, 2 different code snippets have been posted. The code about the volume of the cell is inserted in the createFields.H file, which back then (code 1.7+ to maybe 2.3 or 2.4), existed and basically made field available at runtime.

[Saicharan]

Hi Roman,
I mean in which header file should I add the albanova interface code you specified in your paper? Or should I create a separate header file for adding albanova interface? If I should create one such file where should I add that header file after creating it?
Thanks and regards.

[romant]

Hi Saicharan,


Please have a look at the thesis that describes the code. The code is part of the thesis. http://urn.kb.se/resolve?urn=urn:nbn:se:kth:diva-77980


Best regards,
Roman

[saicharan662000@gmail.com]

Hi roman,
I went through thesis. The albanova file you mentioned in your was said to be extracted from phaseChangeTwoPhaseMixture. Can I know where albanovaInterface.H file is? And also how did u compile the albanovaInterface.C without .H file?
Please help me
Thanks and regards

[Saicharan]

HI roman,
where should I inser the code you gave in interphasechangefoam?
Foam::tmp<Foam::volScalarField>
Foam:haseChangeTwoPhaseMixtures::AlbaNovaInterfa ce::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>("cellVol u");
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)))
);

}
Thanks and regards.
Also thanks for your replies

[romant]

Originally Posted by Saicharan

HI roman,
where should I inser the code you gave in interphasechangefoam?
Foam::tmp<Foam::volScalarField>
Foam:haseChangeTwoPhaseMixtures::AlbaNovaInterfa ce::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>("cellVol u");
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)))
);

}
Thanks and regards.
Also thanks for your replies

Hi,


Please see the attached zip file. It has been a really long time since I wrote this and it was written for OpenFOAM 1.6.x. I cannot make any promises that it works or is correct. Please use this only as an implementation exercise.


It incluse the make files, options files and all header files for the new solver and the models.



Best regards,
Roman

[saicharan662000@gmail.com]

Tank you roman , I will go through them.

[saicharan662000@gmail.com]

Thank you Roman the zip file helped me.

[saicharan662000@gmail.com]

Hi Roman ,
How can you print the interface areas of all the cells using your code. I used it in my own solver myCompressibleInterFoam and simulated sloshing tank test case but I am not getting any output for interface area and interface area is not getting displayed in fields list. can you help me with this.
Thanks in advance.