public double[,] Calculate(TurbulenceModel model, double[,] Ux, double[,] Uy)
{
switch (model)
{
case TurbulenceModel.Secundova:
Secundova(Ux, Uy);
break;
case TurbulenceModel.KE:
KE(Ux, Uy);
break;
}
return(nuT);
}
///<summary>Расчет поля скоростей</summary>
/// <param name="hcm">схема расчета уравнения Гельмгольца</param>
/// <param name="tm">модель турбулентности (turbulenceModel = 0 если турбулентность не расчитывается)</param>
/// <param name="tmax">время расчета</param>
public void Calculation(HelmholtzCalcMethod hcm, TurbulenceModel tm, double tmax)
{
double t = 0;
do
{
if (tm != 0) //turbulenceModel == 0 если турбулентность не расчитывается
{
nuT = turb.Calculate(tm, Ux, Uy);
}
Temp = temp.CalcTemp();
Vortex();
CurrentFunction();
Speeds();
t += tau;
} while (t <= tmax);
}
///<summary>Расчет поля скоростей</summary>
/// <param name="pcm">метод расчета поля давления</param>
/// <param name="nscm">схема расчета уравнения Навье-Стокса</param>
/// <param name="tm">модель турбулентности (turbulenceModel = 0 если турбулентность не расчитывается)</param>
/// <param name="tmax">время расчета</param>
public void Calculation(PressureCalcMethod pcm, NavierStokesCalcMethod nscm, TurbulenceModel tm, double tmax)
{
double t = 0;
do
{
if (tm != 0) //turbulenceModel = 0 если турбулентность не расчитывается
{
nuT = turb.Calculate(tm, Ux, Uy);
}
switch (pcm)
{
case PressureCalcMethod.Poisson:
Poisson();
break;
case PressureCalcMethod.WeakСompressibility:
WeakСompressibility();
break;
}
Temp = temp.CalcTemp();
switch (nscm)
{
case NavierStokesCalcMethod.ExplicitScheme:
ExplicitSpeeds();
break;
case NavierStokesCalcMethod.ImplicitScheme:
ImplicitSpeeds();
break;
}
t += tau;
} while (t <= tmax);
}
public TurbulenceProperties(TurbulenceModel model)
{
this.model = model;
}
