FreeFOAM The Cross-Platform CFD Toolkit
correctViscosity.H
Go to the documentation of this file.
1 {
2  mul = muc +
4  (
5  plasticViscosityCoeff,
6  plasticViscosityExponent,
7  alpha
8  );
9 
10  if (BinghamPlastic)
11  {
13  (
14  yieldStressCoeff,
15  yieldStressExponent,
16  yieldStressOffset,
17  alpha
18  );
19 
20  mul =
21  tauy/
22  (
23  mag(fvc::grad(U))
24  + 1.0e-4*
25  (
26  tauy
28  (
29  "deltaTauy",
30  tauy.dimensions(),
31  1.0e-15
32  )
33  )/mul
34  )
35  + mul;
36  }
37 }
38 
39 // ************************ vim: set sw=4 sts=4 et: ************************ //