[hbulus]

Hello everyone,

I want to initialize my solution with potentialFoam for the sonicFoam. However, since potentialFoam solves for incompressible flow, there will be dimension incompability in fluxes. Therefore, i need to convert volumetric flux (phi from potentialFoam) to mass flux(phi for compressible flow). I found a 'convertPhi.C' code(which is attached) which basically modify Phi by multipliying rho, also for other variables you can modify it later. However, when i try to run it , it asks for rhoRef value. I tried everything but it still asks. I tried to add rhoRef in constant/thermophysicalProperties, system/fvOptions, even in the code by basically adding 'scalar rhoRef(1.0)'. Is there anyone who sees what the problem is in here?

Thanks for any suggestions!

[simrego]

Hi!


Just remove it. Then the solver will create the flux field for you from the velocity field.

[hbulus]

Originally Posted by simrego

Hi!


Just remove it. Then the solver will create the flux field for you from the velocity field.

Thanks for reply simrego,

What do you mean by removing ? Could you explain more ?

[simrego]

Sorry. Remove the phi file from the 0 folder. (if you have a parallel case, from every processor*/0 folder.) At least I think you have it in your 0 folder. Of course if for some reason your start time is not zero, then delete it from the proper time folder.


The solver will reconstruct your flux field from the velocity field (which is written by potentialFoam) correctly and you don't have to do any black magic trick.

[hbulus]

Originally Posted by simrego

Sorry. Remove the phi file from the 0 folder. (if you have a parallel case, from every processor*/0 folder.) At least I think you have it in your 0 folder. Of course if for some reason your start time is not zero, then delete it from the proper time folder.


The solver will reconstruct your flux field from the velocity field (which is written by potentialFoam) correctly and you don't have to do any black magic trick.

Thanks for your quick replies,

I tried as you said but the main solver didn't reconstruct phi, again. In this way, there is no meaning of initialization.I think, i have to find a way to convert volumetric flow to mass flow ..

[hbulus]

How do you make initialization for compressible flows, mr simfego?

[simrego]

There is a meaning since potentialFoam write the velocty field, and an approximated pressure field if needed near the flux field.
The velocity field is always written, which is reconstructed from the flux field.
In the sonicFoam solver you include the compressibleCreatePhi.H file which will read the flux field if present, or create it if not (from the velocity field, as: "linearInterpolate(rho*U) & mesh.Sf()".
Thus you have your flux field from the velocity field which is initialized by the potentialFoam solver mr hbulus!

[hbulus]

Oh, finally it makes sense to use in that way, thanks a lot. I dont solve my problem about convertPhi yet, but i am working on it. If i would solve, i am gonna share with you. Thanks for the all information you gave me.

Have a good day