[eivind]

Hi,
I am currently running a 3D simpleFoam case of a rough channel.
In order to find the appropriate inlet conditions, I have also done a cyclic 2D simulation of a slice of my geometry.

I now want to map the solution at the inlet patch of the cyclic 2D simulation, onto the inlet of the 3D case and use as boundary condition here.
(The 3D simulation is not cyclic due to inhomogenities in spanwise direction.)

Does anynoe have any ideas on how to achieve this?

I have been thinking maybe mapFields may be a part of the solution, but I have never used this utility before...

What do you think?

Any help would be really much appreciated!

Cheers,
Eivind

[piu58]

An easy way to achieve this is approximatng the velocity field by a polynomial and using the uniformValue-polynomial condition for the inlet.

[eivind]

Originally Posted by piu58

An easy way to achieve this is approximatng the velocity field by a polynomial and using the uniformValue-polynomial condition for the inlet.

Thank you for your reply!

I have already tried this method, but it was difficult to fit the turbulence fields (especially epsilon) to polynomials. Then the eddy viscosity may get unrealistic values. It works - but it is not the ideal/most elegant solution.

Anyone familiar with mapFields? Would it do the trick?

[hxaxtma]

My workaround for cases like this is the following:

1) sample your 2D inlet profiles (U, p, etc..)
2) write out the coordinates of your 3D inlet with writecellCentres
3)Load your 2D profile into python
4)Interpolate your profile on the 3D inlet coordinates in Python,
5)Write out a nonuniform vector list and paste this into your 3D inlet U or p file

That's it!

You can also extend this to your internalField and so on!

mapFields cannot do this kind of interpolation, maybe swak4foam is able to, but I do not use this.