FreeFOAM The Cross-Platform CFD Toolkit
fuhsEqn.H
Go to the documentation of this file.
1 {
2  // Solve fuel equation
3  // ~~~~~~~~~~~~~~~~~~~
5 
6  {
7  fvScalarMatrix fuEqn
8  (
9  fvm::ddt(rho, fu)
10  + mvConvection->fvmDiv(phi, fu)
11  - fvm::laplacian(turbulence->alphaEff(), fu)
12  ==
13  R
14  );
15 
16  fuEqn.relax();
17  fuEqn.solve();
18  }
19 
20  Info<< "max(fu) = " << max(fu).value() << endl;
21  Info<< "min(fu) = " << min(fu).value() << endl;
22 
23 
24  // Solve sensible enthalpy equation
25  // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
26  radiation->correct();
27  dQ = combustion->dQ(R);
28 
29  {
30  fvScalarMatrix hsEqn
31  (
32  fvm::ddt(rho, hs)
33  + mvConvection->fvmDiv(phi,hs)
34  - fvm::laplacian(turbulence->alphaEff(), hs)
35  ==
36  DpDt
37  + dQ
38  + radiation->Shs(thermo)
39  );
40 
41  hsEqn.relax();
42  hsEqn.solve();
43  }
44 
45  thermo.correct();
46  combustion->correct();
47 }