public GasParticleVer3(double h, double hmax, List <BorderSegment> boundaries = null, GasParticleOptions opts = null) : base(hmax)
{
options = opts ?? GasParticleOptions.DiffRoDiffE();
Boundaries = boundaries ?? new List <BorderSegment>();
this.h = h;
Vel = new Position2D();
Vel.Name = "Vel";
AddChild(Vel);
AddDiffVect(Vel);
dVdt = new Position2D();
dVdt.Name = "dV";
AddChild(dVdt);
Vel.AddDiffVect(dVdt, false);
if (options.diffRo)
{
AddDiffPropToParam(pRo, pdRo);
}
if (options.diffE)
{
AddDiffPropToParam(pE, pdE);
}
if (options.diffFullE)
{
AddDiffPropToParam(pFullE, pdFullE);
}
}
/// <summary>
///
/// </summary>
/// <param name="initcond">[0]длина трубы;
/// [1]высота трубы;
/// [2]отностиельная граница по Икс
/// [3]Давление слева
/// [4]Плотность слева
/// [5]Давление справа
/// [6]Плотность справа
/// [7]Интервал между частицами
/// </param>
/// <returns></returns>
public static Tuple <IEnumerable <GasParticleVer3>, IEnumerable <IGasParticleVer3> > TestTrubaParticles(params double[] initcond)
{
double lt = initcond[0]; //длина трубы
double ht = initcond[1]; //высота трубы
double x0t = initcond[2]; //отностиельная граница по Икс
double p1t = initcond[3]; //Давление слева
double ro1t = initcond[4]; //Плотность слева
double p2t = initcond[5]; //Давление справа
double ro2t = initcond[6]; //Плотность справа
double delta0 = initcond[7]; //Плотность справа
options = options ?? GasParticleOptions.DiffOnlyE();
var bounds = new List <BorderSegment>(4);
bounds.Add(new BorderSegment(0, 0, 0, ht));
bounds.Add(new BorderSegment(0, ht, lt, ht));
bounds.Add(new BorderSegment(lt, ht, lt, 0));
bounds.Add(new BorderSegment(lt, 0, 0, 0));
int Nx = (int)Round(lt / delta0);
double dx = lt / Nx;
int Ny = (int)Round(ht / delta0);
double dy = ht / Ny;
double delta = Min(dx, dy);
double m1 = GetMass(delta, h_default, ro1t);
double m2 = GetMass(delta, h_default, ro2t);
var particles = new List <GasParticleVer3>(Nx * Ny + 10);
double x = 0.5 * delta;
do
{
double y = 0.5 * delta;
do
{
var p = new GasParticleVer3(h_default, 2 * h_default, bounds);
p.X = x;
p.Y = y;
p.P = x < lt * x0t ? p1t : p2t;
p.Ro = x < lt * x0t ? ro1t : ro2t;
p.M = x < lt * x0t ? m1 : m2;
p.k = k;
particles.Add(p);
y += delta;
} while(y < ht);
x += delta;
} while(x < lt);
particles.ForEach(p => {
p.E = p.P / ((p.k - 1d) * p.Ro);
});
return(new Tuple <IEnumerable <GasParticleVer3>, IEnumerable <IGasParticleVer3> >(particles, null));
}
public static Sph2D GetTest()
{
//string real = "i1.txt";
//string bound = "b1.txt";
//return new My_Sph2D(real,bound);
SPH2D_Ver3.h_default = 1d / 50;
SPH2D_Ver3.k = 3;
SPH2D_Ver3.options = GasParticleOptions.DiffRoDiffE();
//return SPH2D_Ver3.TestTruba(1,0.2,0.5,1,1,0.1,0.125,SPH2D_Ver3.h_default/2);
//return SPH2D_Ver3.TestGasBall(1,1,0);
//return new Sph2D_improoveIntegr(SPH2D_Ver3.TestTrubaParticles(1,0.2,0.5,3E4,1500,1E4,1200,SPH2D_Ver3.h_default));
return(new Sph2D_improoveIntegr(SPH2D_Ver3.TestGasBallParticles(1, 1, 0)));
}
/// <summary>
///
/// </summary>
/// <param name="initcond"> [0]давление внутри шара;
/// [1]плотность внутри шара;
/// [2]скорость шара (мах)</param>
/// <returns></returns>
public static Tuple <IEnumerable <GasParticleVer3>, IEnumerable <IGasParticleVer3> > TestGasBallParticles(params double[] initcond)
{
options = options ?? GasParticleOptions.DiffOnlyE();
double p1 = initcond[0];
double ro1 = initcond[1];
double mach1 = initcond[2];
//
double e_check = 4.29 * 10E6;
ro1 = 1.63;
p1 = 0.65 * e_check;
double E = p1 / ((k - 1d) * ro1);
double Vx = (k * (k - 1) * E) * mach1;
double ht = 1;
double lt = 1;
int Nx = (int)Round(lt / h_default);
double delta = lt / Nx;
int Ny = Nx;
double m1 = GetMass(delta, h_default, ro1);
var particles = new List <GasParticleVer3>(Nx * Ny + 10);
var center = new Particle2DDummyBase(1)
{
X = lt * 0.5,
Y = lt * 0.5
};
var bounds = new List <BorderSegment>(1);
// bounds.Add(new BorderSegment(lt + 0.5 * delta,-0.25 * lt,lt + 0.5 * delta,1.5 * lt));
double x = 0.5 * delta;
double y;
do
{
y = 0.5 * delta;
do
{
var p = new GasParticleVer3(h_default, 2.01 * h_default, bounds);
p.X = x;
p.Y = y;
p.P = p1;
p.Ro = ro1;
p.M = m1;
p.E = E;
p.Vel.X = Vx;
p.k = k;
if (p.GetDistTo(center) < lt * 0.5)
{
particles.Add(p);
}
y += delta;
} while(y < lt);
x += delta;
} while(x < lt);
var borders = new List <IGasParticleVer3>(400);
borders.Add(new GasParticleVer3Dummy(2 * h_default)
{
X = -0.25 * lt * 5, Y = -0.25 * lt * 5, isboundary = true
});
borders.Add(new GasParticleVer3Dummy(2 * h_default)
{
X = -0.25 * lt * 5, Y = 1.25 * lt * 5, isboundary = true
});
//borders.Add(new GasParticleVer3Dummy(2 * h_default) { X = 1.25 * lt,Y = -0.25 * lt,isboundary = true });
//borders.Add(new GasParticleVer3Dummy(2 * h_default) { X = 1.25 * lt,Y = 1.25 * lt,isboundary = true });
y = -0.25 * lt;
//do {
// borders.Add(new GasParticleVer3Dummy(2 * h_default) { X = bounds[0].p1.X + h_default,Y = y,isboundary = true });
// y += (delta / 5);
//} while(y < 1.25 * lt);
return(new Tuple <IEnumerable <GasParticleVer3>, IEnumerable <IGasParticleVer3> >(particles, borders));
}
