[moz455]

Hi,
I am looking at the Theta Johnson Jackson boundary condition. I do not see the rho in the code while we have it in the original formula. I have attached the formula and the related part of the code from JohnsonJacksonParticleThetaFvPatchScalarField.C is:
// calculate the reference value and the value fraction
if (restitutionCoefficient_.value() != 1.0)
{
this->refValue() =
(2.0/3.0)
*specularityCoefficient_.value()
*magSqr(U)
/(scalar(1) - sqr(restitutionCoefficient_.value()));

this->refGrad() = 0.0;

scalarField c
(
constant::mathematical:i
*alpha
*gs0
*(scalar(1) - sqr(restitutionCoefficient_.value()))
*sqrt(3*Theta)
/max(4*kappa*alphaMax.value(), SMALL)
);

this->valueFraction() = c/(c + patch().deltaCoeffs());
}

// for a restitution coefficient of 1, the boundary degenerates to a fixed
// gradient condition
else
{
this->refValue() = 0.0;

this->refGrad() =
pos0(alpha - SMALL)
*constant::mathematical:i
*specularityCoefficient_.value()
*alpha
*gs0
*sqrt(3*Theta)
*magSqr(U)
/max(6*kappa*alphaMax.value(), SMALL);

this->valueFraction() = 0;
}

mixedFvPatchScalarField::updateCoeffs();
}