public void WBVecCopyTest1()
{
WBVec a = new WBVec(1, 2, 3, 4, 5);
var b = a.Clone();
b[1] = 999;
Assert.AreNotEqual(a[1], b[1]);
}
public void WBVecCopyTest3()
{
WBVec a = new WBVec(1, 2, 3, 4, 5);
var b = a;
b[1] = 999;
Assert.AreEqual(a[1], b[1]);
}
public virtual void SetQ(WBVec q)
{
this.q = q;
ro = q[1];
u = q[2] / ro;
e = q[3] / ro - 0.5 * u * u;
p = GetPressure();
Init_h();
}
public static WBVec F_func(WBVec s)
{
var p = (k - 1) * (s[3] - s[2] * s[2] * 0.5 / s[1]);
return(new WBVec(
s[2],
s[2] * s[2] / s[1] + p,
s[2] * (s[3] + p) / s[1]));
}
public WBVec F(double s1, double s2, double s3)
{
s = new WBVec(s1, s2, s3);
ro = s1;
u = s2 / s1;
p = (k - 1) * (s3 - s2 * s2 * 0.5 / s1);
e = GetE();
InitF();
return(f);
}
public RmLayer StepUp(double tau)
{
InitGrUslov_Wall();
var nxtLayer = (RmLayer)Clone();
nxtLayer.Time = Time + tau;
var s12 = new WBVec[Opt.RealNodesCount + 1];
for (int n = 0; n < Nodes.Count - 1; n++)
{
s12[n] = 0.5 * (Nodes[n].s + Nodes[n + 1].s) - tau / 2 / Opt.H * (Nodes[n + 1].f - Nodes[n].f);
}
for (int n = Opt.LeftNodesCount; n < nxtLayer.Nodes.Count - Opt.RightNodesCount; n++)
{
var s_k1 = Nodes[n].s - tau / Opt.H * (RmNode.F_func(s12[n]) - RmNode.F_func(s12[n - 1])) + RmNode.vyaz * (Nodes[n + 1].s - 2d * Nodes[n].s + Nodes[n - 1].s);
nxtLayer.Nodes[n].F(s_k1);
}
return(nxtLayer);
}
public override void SetQ(WBVec q)
{
this.q = q;
ro = q[1];
u = q[2] / ro;
e = q[3] / ro - 0.5 * u * u;
p = GetPressure();
for (int i = 0; i < mixture.powders.Count; i++)
{
if (q[i + 4] > mixture.powders[i].zk * ro)
{
z[i] = mixture.powders[i].zk;
}
else
{
z[i] = q[i + 4] / ro;
}
}
Init_h();
}
public void SetFlux()
{
flux = AUSMp();
}
public GasCell(GasConstants g, int vecDem)
{
this.g = g;
q = WBVec.Zeros(vecDem);
h = WBVec.Zeros(vecDem);
}
public WBVec F(WBVec S)
{
return(F(S[1], S[2], S[3]));
}
public override WBVec AUSMp()
{
var LeftCell = this.LeftCell as OvCell;
var RightCell = this.RightCell as OvCell;
double r1 = LeftCell.ro;
double u1 = LeftCell.u;
double e1 = LeftCell.e;
double r2 = RightCell.ro;
double u2 = RightCell.u;
double e2 = RightCell.e;
double p1 = LeftCell.p;
double p2 = RightCell.p;
double H1 = LeftCell.H;
double H2 = RightCell.H;
double c1 = LeftCell.CSound;
double c2 = RightCell.CSound;
double[] z1 = new double[LeftCell.mixture.powders.Count];
double[] z2 = new double[LeftCell.mixture.powders.Count];
for (int i = 0; i < LeftCell.mixture.powders.Count; i++)
{
z1[i] = LeftCell.z[i];
z2[i] = RightCell.z[i];
}
double cs = 0.5 * (c1 + c2);
double Mr1 = (u1 - V) / cs;
double Mr2 = (u2 - V) / cs;
double du = u2 - u1;
double M4p = Abs(Mr1) >= 1d
? 0.5 * (Mr1 + Abs(Mr1))
: 0.25 * ((Mr1 + 1.0) * (Mr1 + 1.0)) * (1.0 + 2.0 * 0.25 * (Mr1 - 1.0) * (Mr1 - 1.0));
double P5p = Abs(Mr1) >= 1d
? 0.5 * (Mr1 + Abs(Mr1)) / Mr1
: 0.25 * ((Mr1 + 1.0) * (Mr1 + 1.0)) * ((2.0 - Mr1) + 3.0 * Mr1 * 0.25 * (Mr1 - 1.0) * (Mr1 - 1.0));
double M4m = Abs(Mr2) >= 1d
? 0.5 * (Mr2 - Abs(Mr2))
: -0.25 * ((Mr2 - 1.0) * (Mr2 - 1.0)) * (1.0 + 2.0 * 0.25 * (Mr2 + 1.0) * (Mr2 + 1.0));
double P5m = Abs(Mr2) >= 1d
? 0.5 * (Mr2 - Abs(Mr2)) / Mr2
: 0.25 * ((Mr2 - 1.0) * (Mr2 - 1.0)) * ((2.0 + Mr2) - 3.0 * Mr2 * 0.25 * (Mr2 + 1.0) * (Mr2 + 1.0));
double Mrf = M4p + M4m;
double pf = P5p * p1 + P5m * p2;
double flux1 = 0.5 * (cs * Mrf * (r1 + r2) - cs * Abs(Mrf) * (r2 - r1));
double flux2 = 0.5 * (cs * Mrf * (r1 * u1 + r2 * u2) - cs * Abs(Mrf) * (r2 * u2 - r1 * u1)) + pf;
double flux3 = 0.5 * (cs * Mrf * (r1 * H1 + r2 * H2) - cs * Abs(Mrf) * (r2 * H2 - r1 * H1)) + pf * V;
WBVec vec_flux = WBVec.Zeros(LeftCell.mixture.powders.Count + 3);
vec_flux[1] = flux1;
vec_flux[2] = flux2;
vec_flux[3] = flux3;
double[] fluxz = new double[RightCell.mixture.powders.Count];
for (int i = 0; i < RightCell.mixture.powders.Count; i++)
{
fluxz[i] = 0.5 * (cs * Mrf * (r1 * z1[i] + r2 * z2[i]) - cs * Abs(Mrf) * (r2 * z2[i] - r1 * z1[i]));
vec_flux[4 + i] = fluxz[i];
}
return(vec_flux);
}
