public static VisualizedData GetVisualizedData(GeometricParameters geometricParameters, ModeParameters modeParameters, EmpiricalCoefficients empiricalCoefficients, MaterialPropertiesParameters materialPropertiesParameters, double z) { double Y = GetShearStrainRate(modeParameters.Vu, geometricParameters.H); double Qa = GetSpecificHeatFluxQa(geometricParameters.W, modeParameters.Tu, empiricalCoefficients); double Qy = GetSpecificHeatFluxQy(geometricParameters, empiricalCoefficients, Y); double QCH = GetVolumetricFlowRate(geometricParameters, modeParameters); double l = geometricParameters.L; List <Point> temperature = new List <Point>(); List <Point> viscosity = new List <Point>(); for (double i = 0; i < l; i += z) { double T = GetMaterialTemperature(empiricalCoefficients, materialPropertiesParameters, Qa, Qy, QCH, geometricParameters.W, i); double Nu = GetMaterialViscosity(empiricalCoefficients, T, Y); temperature.Add(new Point(i, T)); viscosity.Add(new Point(i, Nu)); } double lastT = GetMaterialTemperature(empiricalCoefficients, materialPropertiesParameters, Qa, Qy, QCH, geometricParameters.W, l); double lastNu = GetMaterialViscosity(empiricalCoefficients, lastT, Y); temperature.Add(new Point(l, lastT)); viscosity.Add(new Point(l, lastNu)); return(new VisualizedData(temperature, viscosity)); }
/// <summary> /// Температура материала в канале /// </summary> /// <param name="empiricalCoefficients">Эмперические коэффициенты математической модели</param> /// <param name="materialPropertiesParameters">Параметры свойств материала</param> /// <param name="Qa">Удельный тепловой поток</param> /// <param name="Qy">Удельный тепловой поток</param> /// <param name="QCH">Объемный расход потока материала в канале</param> /// <param name="W">Ширина</param> /// <param name="z">Шаг</param> /// <returns></returns> public static double GetMaterialTemperature(EmpiricalCoefficients empiricalCoefficients, MaterialPropertiesParameters materialPropertiesParameters, double Qa, double Qy, double QCH, double W, double z) { double c = materialPropertiesParameters.c; double p = materialPropertiesParameters.p; double T0 = materialPropertiesParameters.T0; double b = empiricalCoefficients.b; double alphaU = empiricalCoefficients.AlphaU; double Tr = empiricalCoefficients.Tr; double ratio = (Qa * z) / (p * c * QCH); double sublogarithmicExpression = ((b * Qy + W * alphaU) / (b * Qa)) * (1 - Math.Exp(-b * ratio)) + Math.Exp(b * (T0 - Tr - ratio)); return(Tr + (1 / b) * Math.Log(sublogarithmicExpression)); }