"""
List of parameters required for 3D MRE as solved according to Daitche's paper
"""
TEST_PARAMETERS = dict(
particle_density = 500,
fluid_density = 972,
particle_radius = 0.00018**0.5,
kinematic_viscosity = 2e-4,
time_scale = 0.1,
char_vel = 0.4,
)
[docs]
class DaitcheParameters:
[docs]
def set_beta(self):
self.beta = self.rho_p /self.rho_f
[docs]
def set_S(self):
self.S = (1.0/3.0)*self.a**2/(self.nu * self.T)
[docs]
def set_R(self):
self.R = 3.0/(1.0 + 2.0*self.beta)
[docs]
def set_gravity(self):
self.g = (self.T/self.U)*9.8
def __init__(self, particle_density, fluid_density, particle_radius,
kinematic_viscosity, time_scale, char_vel):
self.rho_p = particle_density
self.rho_f = fluid_density
self.a = particle_radius
self.nu = kinematic_viscosity
self.T = time_scale
self.U = char_vel
self.set_beta()
self.set_S()
self.set_R()
self.set_gravity()