[yapalparvi]

Hi,

I have been able to add energy equation to the interFoam solver and do some test cases with that and it works OK. In this case the surface tension was a constant. Now I would like to make it dependent on temperature. Any help?

Thanks
Ramesh

[isabel]

If you read the files "interfaceProperties.C" and "interfaceProperties.H" you will see how to set the surface tension, sigma, as temperature dependent.

[AlGates]

Think, I have the same problem here. interfaceProperties reads the sigma-values in the constructor, but this happens only at the beginning of the simulation and the constructor is never called again.

Should there be added a new function to be able to change the values of sigma?

[isabel]

I have implemented this code:


volVectorField gradgamma = fvc::grad(gamma);
surfaceVectorField gradgammaf = fvc::interpolate(gradgamma);
surfaceVectorField nHatfv = gradgammaf/(mag(gradgammaf) + SMALL); // normal vector
surfaceScalarField nHatf_ = nHatfv & mesh.Sf();
volScalarField K_ = -fvc::div(nHatf_); //curvature


The momentum equation would change this line:


solve
(
UEqn
==
fvc::reconstruct
(
(
fvc::interpolate(sigma*K_)*fvc::snGrad(gamma)
- ghf*fvc::snGrad(rho)
- fvc::snGrad(pd)
) * mesh.magSf()
)
);


And the pressure equation would change this:


phi = phiU +
(
fvc::interpolate(sigma*K_)*fvc::snGrad(gamma)
- ghf*fvc::snGrad(rho)
)*rUAf*mesh.magSf();

[vahid.najafi]

Hi Dear Isabel.

I want to add surface tension(sigma) in a code,for this added :
#include ''fvCFD.H''
#include ''interfaceproperties.H''
fvc::interpolate(interface.sigma())

in this code:
Foam::tmp<Foam::volScalarField>
Foam:haseChangeTwoPhaseMixtures::SchnerrSauer: Coeff
(
const volScalarField& p
) const
{
volScalarField limitedAlpha1(min(max(alpha1_, scalar(0)), scalar(1)));
volScalarField rho
(
limitedAlpha1*rho1() + (scalar(1) - limitedAlpha1)*rho2()
);
return

//......I want to change it( <<sigma>> surface tension multiple in it):
(3*rho1()*rho2())*sqrt(2/(3*rho1()))* (fvc::interpolate(interface.sigma()))
*rRb(limitedAlpha1)/(rho*sqrt(mag(p - pSat()) + 0.01*pSat()));
//.................................................. ......
}
dont successful wmake, and seen(was not declared ):
phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C:113: error: 'interface' was not declared in this scope
make: *** [Make/linux64GccDPOpt/SchnerrSauer.o] Error 1
please help me,How to correct this problem???

[isabel]

Dear vahid.najafi,

I am not sure which your the problem. Your solver should compile Ok (I am working with OpenFOAM 1.7 version). Try these ideas:

- Write "interfaceProperties.H" instead of "interfaceproperties.H"
- Add the line #include "twoPhaseMiture.H"

If these changes do not work do not hesitate to answer again

[vahid.najafi]

Hi Dear Isabel.
Thanks for your answer.
I worked what ,you said to me :

I Changed solver, As follows:

#include "fvCFD.H"
#include "twoPhaseMixture.H"
#include "interfaceProperties.H"


// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
Foam::tmp<Foam::volScalarField>
Foam:haseChangeTwoPhaseMixtures::SchnerrSauer::r Rb
(
const volScalarField& limitedAlpha1
) const
{
return pow
(
((4*constant::mathematical:i*n_)/3)*fvc::interpolate(interface.sigma())
*limitedAlpha1/(1.0 + alphaNuc() - limitedAlpha1),
1.0/3.0
);
}


but i seen again:
geTwoPhaseMixtures::SchnerrSauer::rRb(const Foam::volScalarField&) const’:
phaseChangeTwoPhaseMixtures/SchnerrSauer/SchnerrSauer.C:83: error: 'interface' was not declared in this scope
make: *** [Make/linux64GccDPOpt/SchnerrSauer.o] Error 1

an other question :

twoPhaseMixture.H in directury in src files:
transportModels/incompressible/incompressibleTwoPhaseMixture \

but sigma is not Definition in this code!!!

The directory above is correct???That is your purpose(for twoPhaseMiture.H)?

regards
Vahid Najafi

[block]

thanks found my problem