[oobiols]

Dear OpenFOAMers,

I am using Openfoam-v6 and I want to solve a simple benchmark, let's say, laminar flow around a cylinder (Re~50).
The thing is that I want to give it my own initial condition. In the 0 folder, I would use a customized U field. How can I tell OpenFOAM to compute the phi scalar field from this U vector field? It should be possible because I have U and the mesh is correctly created using blockMesh. For future iterations, phi is computed and shows up in the timesteps folders, but not for the 0 folder. Any ideas?

I know that the SIMPLE algorithm, for instance, computes phi as :

phi = flux(U)

In older versions of OF, it is computed as:

phi = linearInterpolate(U) & mesh.Sf()

Is it possible to add this as a function object or a function in the controlDict? How?

I have tried to compile my own custom function object to get it but I am not being succesful. I am using the OF docker. Thank you very much!

[oobiols]

Actually, I can make the question more general:

- How to compute the phi field from a U field given by me?

Thanks a lot.