37 frictionalStressModel,
38 SchaefferFrictionalStress,
48 const dictionary& dict
51 frictionalStressModel(dict)
76 *
pow(
Foam::max(alpha - alphaMinFriction, scalar(0)), 10.0);
93 *
pow(
Foam::max(alpha - alphaMinFriction, scalar(0)), 9.0);
106 const scalar I2Dsmall = 1.0e-15;
110 tmp<volScalarField> tmuf
117 alpha.mesh().time().timeName(),
129 if (alpha[celli] > alphaMax.value() - 5
e-2)
132 0.5*pf[celli]*
sin(phi.value())
134 sqrt(1.0/6.0*(
sqr(D[celli].
xx() - D[celli].yy())
135 +
sqr(D[celli].yy() - D[celli].zz())
136 +
sqr(D[celli].zz() - D[celli].
xx()))
137 +
sqr(D[celli].xy()) +
sqr(D[celli].xz())
138 +
sqr(D[celli].yz())) + I2Dsmall
143 muff.correctBoundaryConditions();