public void PerformIteration(int i)
{
//Основной рассчет
//for (int i = 1; i < modellingParams.Nt; i++)
//{
LayerData result = solver.Solve(true);
//Кладем данные в DataSet
u.Append(result.U.ToArray(), "t");
v.Append(result.V.ToArray(), "t");
w.Append(result.W.ToArray(), "t");
T.Append(result.T.ToArray(), "t");
div.Append(result.Div.ToArray(), "t");
time.Append(new double[1] {
(double)i / modellingParams.Nt
});
dataSet.Commit();
//Переходим на следующий слой
solver = new LayerSolver(prevData, result, modellingParams);
prevData = result;
double temp = 0;
int count = result.Width * result.Height * result.Thickness;
for (int ii = 1; ii < result.Width; ii++)
{
for (int jj = 1; jj < result.Height; jj++)
{
for (int kk = 1; kk < result.Thickness; kk++)
{
temp += result.Div[ii, jj, kk];
}
}
}
temp = temp / count * modellingParams.Dx * modellingParams.Dy * modellingParams.Dz;
Console.WriteLine((double)i / modellingParams.Nt * 100 + "% Error = " + temp);
//}
//dataSet.Commit();
}
public void SolveAll(string ctorString)
{
//dataSet = ProxyDataSet.Open("msds:nc?file=../../../temp.nc");
dataSet = ProxyDataSet.Open("msds:memory");
//Инициализируем DataSet
Variable X = dataSet.AddVariable<double>("X", "x");
Variable Y = dataSet.AddVariable<double>("Y", "y");
Variable Z = dataSet.AddVariable<double>("Z", "z");
Variable time = dataSet.AddVariable<double>("Time", "t");
Variable u = dataSet.AddVariable<double>("U velocity", "x", "y", "z", "t");
Variable v = dataSet.AddVariable<double>("V velocity", "x", "y", "z", "t");
Variable w = dataSet.AddVariable<double>("W velocity", "x", "y", "z", "t");
Variable T = dataSet.AddVariable<double>("Temperature", "x", "y", "z", "t");
Variable div = dataSet.AddVariable<double>("Divergence", "x", "y", "z", "t");
dataSet.Commit();
double[] wArr = new double[modellingParams.Nx];
for (int i = 0; i < modellingParams.Nx; i++)
{
wArr[i] = i * modellingParams.Dx;
}
X.PutData(wArr);
wArr = new double[modellingParams.Ny];
for (int i = 0; i < modellingParams.Ny; i++)
{
wArr[i] = i * modellingParams.Dy;
}
Y.PutData(wArr);
wArr = new double[modellingParams.Nz];
for (int i = 0; i < modellingParams.Nz; i++)
{
wArr[i] = i * modellingParams.Dz;
}
Z.PutData(wArr);
//Инициализируем рассчетный модуль для слоя начальными условиями
LayerSolver solver = new LayerSolver(prevData, modellingParams);
u.Append(prevData.U.ToArray(), "t");
v.Append(prevData.V.ToArray(), "t");
w.Append(prevData.W.ToArray(), "t");
T.Append(prevData.T.ToArray(), "t");
div.Append(prevData.Div.ToArray(), "t");
time.PutData(new double[1] { 0 });
dataSet.Commit();
//Основной рассчет
for (int i = 1; i < modellingParams.Nt; i++)
{
LayerData result = solver.Solve(true);
//Кладем данные в DataSet
u.Append(result.U.ToArray(), "t");
v.Append(result.V.ToArray(), "t");
w.Append(result.W.ToArray(), "t");
T.Append(result.T.ToArray(), "t");
div.Append(result.Div.ToArray(), "t");
time.Append(new double[1] { (double)i / modellingParams.Nt });
dataSet.Commit();
//Переходим на следующий слой
solver = new LayerSolver(prevData, result, modellingParams);
prevData = result;
double temp = 0;
int count = 0;
for (int ii = 1; ii < result.Width; ii++)
{
for (int jj = 1; jj < result.Height; jj++)
{
for (int kk = 1; kk < result.Thickness; kk++)
{
temp += result.Div[ii, jj, kk];
count++;
}
}
}
temp = temp / count * modellingParams.Dx * modellingParams.Dy * modellingParams.Dz;
Console.WriteLine((double)i / modellingParams.Nt * 100 + "% Error = " + temp);
}
dataSet.Commit();
}
public void SolveAll(string ctorString)
{
//dataSet = ProxyDataSet.Open("msds:nc?file=../../../temp.nc");
dataSet = ProxyDataSet.Open("msds:memory");
//Инициализируем DataSet
Variable X = dataSet.AddVariable <double>("X", "x");
Variable Y = dataSet.AddVariable <double>("Y", "y");
Variable Z = dataSet.AddVariable <double>("Z", "z");
Variable time = dataSet.AddVariable <double>("Time", "t");
Variable u = dataSet.AddVariable <double>("U velocity", "x", "y", "z", "t");
Variable v = dataSet.AddVariable <double>("V velocity", "x", "y", "z", "t");
Variable w = dataSet.AddVariable <double>("W velocity", "x", "y", "z", "t");
Variable T = dataSet.AddVariable <double>("Temperature", "x", "y", "z", "t");
Variable div = dataSet.AddVariable <double>("Divergence", "x", "y", "z", "t");
dataSet.Commit();
double[] wArr = new double[modellingParams.Nx];
for (int i = 0; i < modellingParams.Nx; i++)
{
wArr[i] = i * modellingParams.Dx;
}
X.PutData(wArr);
wArr = new double[modellingParams.Ny];
for (int i = 0; i < modellingParams.Ny; i++)
{
wArr[i] = i * modellingParams.Dy;
}
Y.PutData(wArr);
wArr = new double[modellingParams.Nz];
for (int i = 0; i < modellingParams.Nz; i++)
{
wArr[i] = i * modellingParams.Dz;
}
Z.PutData(wArr);
//Инициализируем рассчетный модуль для слоя начальными условиями
LayerSolver solver = new LayerSolver(prevData, modellingParams);
u.Append(prevData.U.ToArray(), "t");
v.Append(prevData.V.ToArray(), "t");
w.Append(prevData.W.ToArray(), "t");
T.Append(prevData.T.ToArray(), "t");
div.Append(prevData.Div.ToArray(), "t");
time.PutData(new double[1] {
0
});
dataSet.Commit();
//Основной рассчет
for (int i = 1; i < modellingParams.Nt; i++)
{
LayerData result = solver.Solve(true);
//Кладем данные в DataSet
u.Append(result.U.ToArray(), "t");
v.Append(result.V.ToArray(), "t");
w.Append(result.W.ToArray(), "t");
T.Append(result.T.ToArray(), "t");
div.Append(result.Div.ToArray(), "t");
time.Append(new double[1] {
(double)i / modellingParams.Nt
});
dataSet.Commit();
//Переходим на следующий слой
solver = new LayerSolver(prevData, result, modellingParams);
prevData = result;
double temp = 0;
int count = 0;
for (int ii = 1; ii < result.Width; ii++)
{
for (int jj = 1; jj < result.Height; jj++)
{
for (int kk = 1; kk < result.Thickness; kk++)
{
temp += result.Div[ii, jj, kk];
count++;
}
}
}
temp = temp / count * modellingParams.Dx * modellingParams.Dy * modellingParams.Dz;
Console.WriteLine((double)i / modellingParams.Nt * 100 + "% Error = " + temp);
}
dataSet.Commit();
}