Commit f719c0df authored by sergio's avatar sergio
Browse files

BUG: porousSolid T Ew modifications and lookup of T instead of the

thermo
parent 011b804b
...@@ -9,7 +9,7 @@ ...@@ -9,7 +9,7 @@
tmp<volScalarField> tcp = thermo.Cp(); tmp<volScalarField> tcp = thermo.Cp();
const volScalarField& cp = tcp(); const volScalarField& cp = tcp();
tmp<volScalarField> tkappa = betav*thermo.K(); tmp<volScalarField> tkappa = thermo.K();
//tmp<volSymmTensorField> tkappa = thermo.directionalK()*betav; //tmp<volSymmTensorField> tkappa = thermo.directionalK()*betav;
const volScalarField& kappa = tkappa(); const volScalarField& kappa = tkappa();
......
...@@ -3,7 +3,7 @@ ...@@ -3,7 +3,7 @@
{ {
tmp<fvScalarMatrix> TEqn tmp<fvScalarMatrix> TEqn
( (
- fvm::laplacian(kappa, T, "laplacian(K,T)") - fvm::laplacian(betav*kappa, T, "laplacian(K,T)")
+ sources(rho, T) + sources(rho, T)
); );
......
...@@ -9,7 +9,7 @@ ...@@ -9,7 +9,7 @@
tmp<volScalarField> tcp = thermo.Cp(); tmp<volScalarField> tcp = thermo.Cp();
const volScalarField& cp = tcp(); const volScalarField& cp = tcp();
tmp<volScalarField> tkappa = betav*thermo.K(); tmp<volScalarField> tkappa = thermo.K();
//tmp<volSymmTensorField> tkappa = thermo.directionalK()*betav; //tmp<volSymmTensorField> tkappa = thermo.directionalK()*betav;
const volScalarField& kappa = tkappa(); const volScalarField& kappa = tkappa();
......
...@@ -9,7 +9,7 @@ if (finalIter) ...@@ -9,7 +9,7 @@ if (finalIter)
tmp<fvScalarMatrix> TEqn tmp<fvScalarMatrix> TEqn
( (
fvm::ddt(betav*rho*cp, T) fvm::ddt(betav*rho*cp, T)
- fvm::laplacian(kappa, T, "laplacian(K,T)") - fvm::laplacian(betav*kappa, T, "laplacian(K,T)")
+ sources(rho, T) + sources(rho, T)
); );
......
...@@ -167,16 +167,13 @@ void Foam::interRegionHeatTransferModel::addSup ...@@ -167,16 +167,13 @@ void Foam::interRegionHeatTransferModel::addSup
const fvMesh& secondaryMesh = const fvMesh& secondaryMesh =
mesh_.time().lookupObject<fvMesh>(mapRegionName_); mesh_.time().lookupObject<fvMesh>(mapRegionName_);
basicThermo& secondaryThermo = const volScalarField& Tsecondary =
const_cast<basicThermo&> secondaryMesh.lookupObject<volScalarField>("T");
(
secondaryMesh.lookupObject<basicThermo>("thermophysicalProperties")
);
secondaryToPrimaryInterpPtr_->interpolateInternalField secondaryToPrimaryInterpPtr_->interpolateInternalField
( (
Tmapped, Tmapped,
secondaryThermo.T(), Tsecondary,
meshToMesh::MAP, meshToMesh::MAP,
eqOp<scalar>() eqOp<scalar>()
); );
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment