[calandrea]

Hello,

I am trying to implement the SIMPLE scheme for multiphase flows.
I have started combining the already available solvers rasInterFoam and transientSimpleFoam, unfortunately without success.

The error message I get is:

temporary deallocated
From function T& tmp<t>::operator()()
in file /home/califano/OpenFOAM/OpenFOAM-1.5/src/OpenFOAM/lnInclude/tmpI.H at line 164

I think this is due to the way UEqn is solved for multiphase flows:

UEqn == fvc::reconstruct( (...)*mesh.magSf())

In fact, I don't get any error substituting the UEqn for incompressible solvers:

UEqn() == -fvc::grad(p)

It seems that the error is related to the way the tmp class is defined, but my skills in OpenFOAM and C++ are too poor to let me go deeper into the problem.

Anyone can help me with this problem and/or give me a hint to where I can look for a solution?
Thank in advance,
Andrea