// расчёт вязкости материала
private static double Viscosity(double T, InputDataObject inpObj)
{
return(inpObj.EmpiricalVector.Mu * Math.Exp(-inpObj.EmpiricalVector.b * (T - inpObj.EmpiricalVector.Tgamma))
* Math.Pow((inpObj.ModeVector.Vu / inpObj.GeometryVector.H), inpObj.EmpiricalVector.n - 1));
}
private static double Temperature(double Q, double qa, double qr, double coord, InputDataObject inpObj)
{
return(inpObj.EmpiricalVector.Tgamma + (1 / inpObj.EmpiricalVector.b) * Math.Log(((inpObj.EmpiricalVector.b * qr + inpObj.GeometryVector.W *
inpObj.EmpiricalVector.alphaU) /
(inpObj.EmpiricalVector.b * qa)) *
(1 - Math.Exp(-(inpObj.EmpiricalVector.b * qa * coord) /
(inpObj.MaterialProperties.Ro * inpObj.MaterialProperties.C * Q)))
+
Math.Exp(inpObj.EmpiricalVector.b *
(inpObj.MaterialProperties.T0 - inpObj.EmpiricalVector.Tgamma -
(qa * coord) / (inpObj.MaterialProperties.Ro * inpObj.MaterialProperties.C * Q)))));
}
