FreeFOAM The Cross-Platform CFD Toolkit
UEqn.H
Go to the documentation of this file.
2  (
3  "muEff",
4  twoPhaseProperties->muf()
6  );
7 
9  (
10  fvm::ddt(rho, U)
11  + fvm::div(rhoPhi, U)
14  - (fvc::grad(U) & fvc::grad(muEff))
15  //- fvc::div(muEff*(fvc::interpolate(dev2(fvc::grad(U))) & mesh.Sf()))
16  );
17 
18  if (finalIter)
19  {
20  UEqn.relax(1);
21  }
22  else
23  {
24  UEqn.relax();
25  }
26 
28  {
29  solve
30  (
31  UEqn
32  ==
34  (
35  (
36  fvc::interpolate(interface.sigmaK())*fvc::snGrad(alpha1)
39  ) * mesh.magSf()
40  ),
41  mesh.solver(U.select(finalIter))
42  );
43  }
44 
45 // ************************ vim: set sw=4 sts=4 et: ************************ //