Esempio n. 1
0
 public My_IsotropicGasMiracleDummy(My_IsotropicGas parent, bool boundarytype = false, double[] cond = null)
 {
     this.parent = parent;
     Vel         = new VelDummy();
     //cond -  плотность давление масса скорость (x,y)
     isboundary = boundarytype;
     if (cond != null)
     {
         Ro    = cond[0];
         P     = cond[1];
         M     = cond[2];
         Vel.X = cond[3];
         Vel.Y = cond[4];
         E     = P / (0.4 * Ro);
     }
 }
Esempio n. 2
0
        private static List <My_IsotropicGas> ConstructBountPart(double boundaryDensity)
        {
            List <My_IsotropicGas> WallParticles = new List <My_IsotropicGas>();

            foreach (var side in My_Sph2D.boundaries)
            {
                Vector2D vector = side.Item2 - side.Item1;
                int      partN  = Convert.ToInt32((Math.Abs(vector.GetLength())) * boundaryDensity);
                double   xstep  = vector.X / partN;
                double   ystep  = vector.Y / partN;
                for (int i = 0; i < partN; ++i)
                {
                    My_IsotropicGas a = new My_IsotropicGas(1, 1, true);
                    a.X = side.Item1.X + i * xstep;
                    a.Y = side.Item1.Y + i * ystep;
                    WallParticles.Add(a);
                }
            }
            return(WallParticles);
        }