[shyam]

Hi Foamers,

I need to use interFoam to solve a problem, but it blows up within few time steps, which I believe is because of poor initialization . Hence I used simpleFoam to obtain the velocity and pressure fields. Now I want use the pressure field of simpleFoam (with units m^2/S^2) to obtain the pressure field of interFoam (kg/m/s^2) which is nothing but the product of pressure and density. How do I obtain this in openfoam?

[nimasam]

some hints you need to investigate yourself
1) Are you familiar with funkySetFields ? its great tool to make non uniform field
2) as a little I know ,you can easily change the dimension in p file ,so you can apply this pressure to interFoam
3)if you want to use pSimple*rho as new in interFoam, you should modified your solver!
something like that
volScalarField pSimple
(
IOobject
(
"pSimple",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::NO_WRITE
),
mesh
);

p = pSimple*rho;

[Bernhard]

Originally Posted by nimasam

2) as a little I know ,you can easily change the dimension in p file ,so you can apply this pressure to interFoam

But then your off by a factor of the density is suppose?

Originally Posted by nimasam

3)if you want to use pSimple*rho as new in interFoam, you should modified your solver!
something like that

If you do this one, I assume the foamCalc utility can do the trick.

[shyam]

Thanks nima and Bernhard for your quick thoughts. I will try option 2 of nima and foamCalc and post my experience.