[Mohadali]

Hi all,


I am trying to read a mesh from another case in order to interpolate between both cases (fine and coarse)


Foam::fvMesh meshSource
(
Foam::IOobject
(
Foam::fvMesh::defaultRegion,
path, // I am giving the exact path for the source caseDirectory
mesh,
Foam::IOobject::MUST_READ
)
);


I tried a lot of possibilities without succeed.


Can someone help please?

[ashishmagar600]

Hi

If you just want to interpolate between the same geometries but different mesh sizes, you can use mapFields utility.

Eg: you want to interpolate values for fine grid case, use:

mapFields /path/to/coarse/grid -consistent

The "-consistent" arg is used if both the meshes with patches and BCs are identical.

Try

mapfields -help

[Mohadali]

Thank you Ashish,


I am trying to interpolate at each time step, so the mapFielads utility is not useful here.


Any other suggestion.
Thx

[ashishmagar600]

It could be used if you just want to interpolate between two meshes.

I did it once for simulation of 15s (save time 0.001; i.e, 15000 timesteps) for
separate coarse and fine mesh.

You could write a bash script for creating a new time folder & use mapFields -sourceTime ../pathToCase in a loop. The mapFields utility maps -sourceTime specified to the latest Time folder specified in current dir.

Hope this helps.

[Mohadali]

Thank you Ashish for your reply,


I am trying to do interpolation inside the time loop to avoid the bash scripts.

I saw that people are doing it, but my problem is how to define a new mesh (as foam object) inside the main solver.

[ashishmagar600]

I think using a script to run mapFields is easier.

Anyway, you could take a look at mapFields source code, you might get some idea how do they do it and implement it in your code.

[Mohadali]

Thank you Ashish for your help