def __init__(self, time_divisor=0.5):
self.c = Constants(drop_species='water', gas_species='air', T_G=298,
rho_G=1.184, C_p_G=1007, Re_d=0)
self.c.drop_properties()
self.c.gas_properties()
self.c.get_reference_conditions()
self.c.add_drop_properties()
self.c.add_gas_properties()
self.c.add_properties()
self.p1 = particle(self.c, [0, 0, 0], velocity=[0, 0, 0],
D=np.sqrt(1.1)/1000, T_d=282,
ODE_solver=0, coupled=0)
self.p2 = particle(self.c, [0, 0, 0], velocity=[0, 0, 0],
D=np.sqrt(1.1)/1000, T_d=282,
ODE_solver=1, coupled=0)
self.p3 = particle(self.c, [0, 0, 0], velocity=[0, 0, 0],
D=np.sqrt(1.1)/1000, T_d=282,
ODE_solver=2, coupled=0)
self.p4 = particle(self.c, [0, 0, 0], velocity=[0, 0, 0],
D=np.sqrt(1.1)/1000, T_d=282,
ODE_solver=3, coupled=0)
self.p5 = particle(self.c, [0, 0, 0], velocity=[0, 0, 0],
D=np.sqrt(1.1)/1000, T_d=282,
ODE_solver=4, coupled=0)
self.time_divisor = time_divisor
self.div = self.p2.get_tau()*(self.time_divisor)
self.N = 10000
self.p2_error = []
self.p3_error = []
self.p4_error = []
self.p5_error = []
def __init__(self, drop_species='decane', gas_species='air', T_G=1000,
rho_G=0.3529, C_p_G=1135, Re_d=17, T_d=315, D=0.004):
self.drop_species = drop_species
self.gas_species = gas_species
self.T_G = T_G
self.rho_G = rho_G
self.C_p_G = C_p_G
self.Re_d = Re_d
self.T_d = T_d
self.D = D
self.c = Constants(self.drop_species, self.gas_species, self.T_G,
self.rho_G, self.C_p_G, self.Re_d)
self.c.drop_properties()
self.c.gas_properties()
self.c.get_reference_conditions()
self.c.add_drop_properties()
self.c.add_gas_properties()
self.c.add_properties()
self.p = particle(self.c, [0, 0, 0], velocity=[0, 0, 0],
D=self.D, T_d=self.T_d,
ODE_solver=2, coupled=2)
self.div = self.p.get_tau()/32
self.N = 10000
def __init__(self, time_divisor=0.0625):
self.c = Constants()
self.c.drop_properties()
self.c.gas_properties()
self.c.get_reference_conditions()
self.c.add_drop_properties()
self.c.add_gas_properties()
self.c.add_properties()
self.p2 = particle(self.c, [0, 0, 0], velocity=[0, 0, 0],
D=np.sqrt(1.1)/1000, T_d=282,
ODE_solver=1, coupled=2)
self.p3 = particle(self.c, [0, 0, 0], velocity=[0, 0, 0],
D=np.sqrt(1.1)/1000, T_d=282,
ODE_solver=2, coupled=2)
self.p4 = particle(self.c, [0, 0, 0], velocity=[0, 0, 0],
D=np.sqrt(1.1)/1000, T_d=282,
ODE_solver=3, coupled=2)
self.p5 = particle(self.c, [0, 0, 0], velocity=[0, 0, 0],
D=np.sqrt(1.1)/1000, T_d=282,
ODE_solver=4, coupled=2)
self.time_divisor = time_divisor
self.div = self.p2.get_tau()*self.time_divisor
self.N = 10000