solvePorousSolid.H 343 Bytes
Newer Older
1
2
3
4
5
{
    for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
    {
        tmp<fvScalarMatrix> TEqn
        (
6
          - fvm::laplacian(betav*kappa, T, "laplacian(K,T)")
7
8
9
10
11
12
13
14
15
          + sources(rho, T)
        );

        TEqn().relax();
        TEqn().solve();
    }

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