[jayson]

Hello Foamers!
I am fairly newbie.. I want to simulate the combination of two phase flow and conjugate heat transfer. So, I now add solid region to interfoam which has energy equation.

I create the mesh for solid region and the fields for solid.
but I get error.. I think It is fairly correct code..

Foam::Info
<< "Create solid mesh for time ="
<< runTime.timeName() <<Foam::nl <<Foam::endl;

Foam::fvMesh solidMesh
(
Foam::IOobject
(
"wallRegion",
runTime.timeName(),
runTime,
Foam::IOobject::MUST_READ
)
);

Info<< "Reading solid thermophysical properties\n" <<endl;
autoPtr<solidThermo> wWallThermo(solidThermo::New(solidMesh));
solidThermo& wallThermo = wWallThermo();

Info<< "Adding to solidRho\n" <<endl;
volScalarField wallRho
(
IOobject
(
"solidRho",
runTime.timeName(),
solidMesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
wallThermo.rho()
);

Info<< "Adding to solild kappa\n" <<endl;
volScalarField wallK
(
IOobject
(
"solidK",
runTime.timeName(),
solidMesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
wallThermo.Kappa()
);

Info<< "Adding to solidcp\n" <<endl;
volScalarField wallcp
(
IOobject
(
"solidcp",
runTime.timeName(),
solidMesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
wallThermo.Cp()
);


Info<< "solve solid TEqn" << endl;
{
wallRho = wallThermo.rho();
wallK = wallThermo.kappa();
wallcp = wallThermo.Cp();

fvScalarMatrix TEqn
(
fvm::ddt(wallRho*wallcp, T)
-fvm::laplacian(wallK, T)
);

TEqn.relax();

TEqn.solve();

wallThermo.correct();

Info<< "Min/max solid T:" << min(wallThermo.T()).value() << ' '
<< max(wallThermo.T()).value() << endl;

}


/solid/createSolidFields.H: In function ‘int main(int, char**)’:
./solid/createSolidFields.H:31:1: error: no matching function for call to ‘Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject, Foam::tmp<Foam::GeometricField<Foam::Vector<double >, Foam::fvPatchField, Foam::volMesh> >)’
);