[kebsiali]

Hello Dear Foamers

I have a simple question
in interFoam in createFields.H the rho is constructed as

volScalarField rho
(
IOobject
(
"rho",
runTime.timeName(),
mesh,
IOobject::READ_IF_PRESENT
),
alpha1*rho1 + alpha2*rho2,
alpha1.boundaryField().types()
);
rho.oldTime();


my question is what do the last couple of commands do
1. alpha1.boundaryField().types()
2. rho.oldTime();


I'm trying to implement a variable rho1 and I'm comparing interFoam with compressibleInterFoam the later does not have those two commands which is why i was wondering, i have nearly 20 self made constructors which are created randomly (without knowing if i have to add somthing like blabla.oldTime(); or blabla.boundaryFiled().types() or not)

[kebsiali]

please!
any response, suggestions ?

[kebsiali]

Hello everybody

I'm working with interFoam and i want to make a varibal phase density
for this task i started by constructing a volScalarField in creatFields as

volScalarField rhof
(
IOobject
(
"rhof",
runTime.timeName(),
mesh,
IOobject::NO_READ
),
mesh
);

before i put some equation instead of "mesh" but it didnt work because

when i call this rhof into the library incompressibleTwoPhaseMixture.C using
const volScalarField& rhof=U_.mesh().lookupObject<volScalarField>("rhof" );
so that rhof takes place of rho1
now the problem reported by terminal is that rhof is not an object of mesh
this appears when i execute the case as it gives other suggestions that are useable objects and all of them are those for which in the constructor it is as the one shown above with mesh as the last line.

now using this constructor shown above with mesh in the last line so that rhof can be callabe the code searches for a boundary condition for rhof

am i going in the right way or should i just delete everything and give up, nah just joking.

any suggestions are heartfully wellcome

[babakflame]

Originally Posted by kebsiali

Hello Dear Foamers

I have a simple question
in interFoam in createFields.H the rho is constructed as

volScalarField rho
(
IOobject
(
"rho",
runTime.timeName(),
mesh,
IOobject::READ_IF_PRESENT
),
alpha1*rho1 + alpha2*rho2,
alpha1.boundaryField().types()
);
rho.oldTime();


my question is what do the last couple of commands do
1. alpha1.boundaryField().types()
2. rho.oldTime();


I'm trying to implement a variable rho1 and I'm comparing interFoam with compressibleInterFoam the later does not have those two commands which is why i was wondering, i have nearly 20 self made constructors which are created randomly (without knowing if i have to add somthing like blabla.oldTime(); or blabla.boundaryFiled().types() or not)

Hi Fellows

Any progress in this question?

I have the same question.

[santiagomarquezd]

Originally Posted by kebsiali

Hello Dear Foamers

I have a simple question
in interFoam in createFields.H the rho is constructed as

volScalarField rho
(
IOobject
(
"rho",
runTime.timeName(),
mesh,
IOobject::READ_IF_PRESENT
),
alpha1*rho1 + alpha2*rho2,
alpha1.boundaryField().types()
);
rho.oldTime();


my question is what do the last couple of commands do
1. alpha1.boundaryField().types()
2. rho.oldTime();


I'm trying to implement a variable rho1 and I'm comparing interFoam with compressibleInterFoam the later does not have those two commands which is why i was wondering, i have nearly 20 self made constructors which are created randomly (without knowing if i have to add somthing like blabla.oldTime(); or blabla.boundaryFiled().types() or not)

Hi, Kebsiali, regarding these:

   42 // Need to store rho for ddt(rho, U)
   43 volScalarField rho
   44 (
   45     IOobject
   46     (
   47         "rho",
   48         runTime.timeName(),
   49         mesh,
   50         IOobject::READ_IF_PRESENT
   51     ),
   52     alpha1*rho1 + alpha2*rho2,
   53     alpha1.boundaryField().types()
   54 );
   55 rho.oldTime();

in createFields.H alpha1.boundaryField().types() implies copying the boundary conditions from alpha1, and rho.oldTime() forces to call the storeOldTimes() method in order to set the old time field required by ddt(rho, U) as is indicated in 42.

Regards!