exaFOAM/foam-extend4.1-coherent-io
Archived
1
1
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Bugfix: clean-up

Browse source
This commit is contained in:
Hrvoje Jasak 2016-04-20 12:55:08 +01:00
parent 89dd971095
commit 25ac1973ca
4 changed files with 9 additions and 10 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
applications/solvers/coupled/conjugateHeatSimpleFoam/conjugateHeatSimpleFoam.C
View file

@ -81,17 +81,17 @@ int main(int argc, char *argv[])
// Update thermal conductivity in the solid // Update thermal conductivity in the solid
solidThermo.correct(); solidThermo.correct();
ksolid = solidThermo.k(); kSolid = solidThermo.k();
// Coupled patches // Coupled patches
# include "attachPatches.H" # include "attachPatches.H"
kappaEff.correctBoundaryConditions(); kappaEff.correctBoundaryConditions();
ksolid.correctBoundaryConditions(); kSolid.correctBoundaryConditions();
// Interpolate to the faces and add thermal resistance // Interpolate to the faces and add thermal resistance
surfaceScalarField ksolidf = fvc::interpolate(ksolid); surfaceScalarField kSolidf = fvc::interpolate(kSolid);
solidThermo.modifyResistance(ksolidf); solidThermo.modifyResistance(kSolidf);
# include "solveEnergy.H" # include "solveEnergy.H"

4
applications/solvers/coupled/conjugateHeatSimpleFoam/createFields.H
View file

@ -100,7 +100,9 @@
( (
"rhok", "rhok",
runTime.timeName(), runTime.timeName(),
mesh mesh,
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
), ),
1.0 - beta*(T - TRef) 1.0 - beta*(T - TRef)
); );

2
applications/solvers/coupled/conjugateHeatSimpleFoam/createSolidFields.H
View file

@ -16,7 +16,7 @@
thermalModel solidThermo(Tsolid); thermalModel solidThermo(Tsolid);
Info<< "Reading solid diffusivity k\n" << endl; Info<< "Reading solid diffusivity k\n" << endl;
volScalarField ksolid volScalarField kSolid
( (
IOobject IOobject
( (

5
applications/solvers/coupled/conjugateHeatSimpleFoam/solveEnergy.H
View file

@ -20,19 +20,16 @@
+ 3.0*radiation->Rp()*pow4(T) + 3.0*radiation->Rp()*pow4(T)
); );
//Called automatically in 1.6.1-ext
//TFluidEqn->boundaryManipulate(T.boundaryField());
TFluidEqn->relax(); TFluidEqn->relax();
fvScalarMatrix* TSolidEqn = new fvScalarMatrix fvScalarMatrix* TSolidEqn = new fvScalarMatrix
( (
- fvm::laplacian(ksolidf, Tsolid, "laplacian(k,T)") - fvm::laplacian(kSolidf, Tsolid, "laplacian(k,T)")
+ fvm::SuSp(-solidThermo.S()/Tsolid, Tsolid) + fvm::SuSp(-solidThermo.S()/Tsolid, Tsolid)
); );
TSolidEqn->relax(); TSolidEqn->relax();
// Add fluid equation // Add fluid equation
TEqns.set(0, TFluidEqn); TEqns.set(0, TFluidEqn);