solveSolid.H 482 Bytes
Newer Older
1
2
3
4
5
if (finalIter)
{
    mesh.data::add("finalIteration", true);
}

6
7
8
{
    for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
    {
9
        tmp<fvScalarMatrix> TEqn
10
        (
11
            fvm::ddt(rho*cp, T)
12
          - fvm::laplacian(kappa, T)
13
        );
14
        TEqn().relax();
15
        TEqn().solve(mesh.solver(T.select(finalIter)));
16
    }
mattijs's avatar
new bc    
mattijs committed
17

Andrew Heather's avatar
Andrew Heather committed
18
    Info<< "Min/max T:" << min(T) << ' ' << max(T) << endl;
19
}
20
21
22
23
24
25
26

thermo.correct();

if (finalIter)
{
    mesh.data::remove("finalIteration");
}