[Marcel_pmd]

Hey everyone,

I have a very basic Problem with the implementation of the thermophysical properties library to the FSI pack of foamExtend 4.0.

In general I want to write a compressible FSI solver to work with some high pressure fluid Models I found in literature. Therefore I implemented these models into the thermophysical Properties library and got them to work correctly in OF 50 with some regular cases of rhopiso and rhosimple Foam.

In the next step I wanted to implement one of the compressible solvers I used ( like rhopiso or rhosimple) to the FSI package in FE40. To do so I need to make the thermophysical model library available to the solver.

As soon as I implement the lines

//- Thermophysical model
autoPtr<basicPsiThermo> pThermo_;

to pisoFluid.H

and

pThermo_
(
basicPsiThermo::New(mesh)
),

to pisoFluid.C in the same way the turbulence model is implemented, the pressure field isn't written correctly any more in the testcase "beamInCrossFlow" for the fsiFoam solver. It only writes the initital 0 Field for the pressure again and again. The velocity field however is written and solved. If I look up the basicThermo.H and .C file there is a definition for pressure etc.
My guess is, that there is a conflict with the namespace definitions.

Does anyone has experience with this kind of problem of compressibility implementation in FSI ?

Every suggestion and idea is welcome.