BUG: UListIO: byteSize overflowing on really big faceLists
[OpenFOAM-2.0.x.git] / applications / solvers / combustion / fireFoam / YhsEqn.H
blobaed41a4431359997c531556cab966fd23d2c514e
1 tmp<fv::convectionScheme<scalar> > mvConvection
3     fv::convectionScheme<scalar>::New
4     (
5         mesh,
6         fields,
7         phi,
8         mesh.divScheme("div(phi,Yi_hs)")
9     )
12     radiation->correct();
13     combustion->correct();
14     dQ = combustion->dQ();
15     label inertIndex = -1;
16     volScalarField Yt(0.0*Y[0]);
18     forAll(Y, i)
19     {
20         if (Y[i].name() != inertSpecie)
21         {
22             volScalarField& Yi = Y[i];
23             fvScalarMatrix R(combustion->R(Yi));
25             fvScalarMatrix YiEqn
26             (
27                 fvm::ddt(rho, Yi)
28               + mvConvection->fvmDiv(phi, Yi)
29               - fvm::laplacian(turbulence->alphaEff(), Yi)
30               ==
31                 parcels.SYi(i, Yi)
32               + surfaceFilm.Srho(i)
33               + R
34             );
36             YiEqn.relax();
37             YiEqn.solve(mesh.solver("Yi"));
39             Yi.max(0.0);
40             Yt += Yi;
41         }
42         else
43         {
44             inertIndex = i;
45         }
46     }
48     Y[inertIndex] = scalar(1) - Yt;
49     Y[inertIndex].max(0.0);
51     fvScalarMatrix hsEqn
52     (
53         fvm::ddt(rho, hs)
54       + mvConvection->fvmDiv(phi, hs)
55       - fvm::laplacian(turbulence->alphaEff(), hs)
56      ==
57         DpDt
58       + dQ
59       + radiation->Shs(thermo)
60       + parcels.Sh(hs)
61       + surfaceFilm.Sh()
62     );
64     hsEqn.relax();
65     hsEqn.solve();
67     thermo.correct();
69     Info<< "min/max(T) = " << min(T).value() << ", " << max(T).value() << endl;