solveSolid.H 635 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> hEqn
10
        (
11
12
            fvm::ddt(betav*rho, h)
          - fvm::laplacian(betav*alpha, h, "laplacian(alpha,h)")
13
14
          ==
            sources(rho, h)
15
        );
16

17
        hEqn().relax();
18
19
20

        sources.constrain(hEqn());

21
        hEqn().solve(mesh.solver(h.select(finalIter)));
22
23

        sources.correct(h);
24
25
    }
}
26
27
28

thermo.correct();

29
30
Info<< "Min/max T:" << min(thermo.T()) << ' ' << max(thermo.T()) << endl;

31
32
33
34
if (finalIter)
{
    mesh.data::remove("finalIteration");
}