def test_Cfd(self):
''' Test for Cfd '''
gravity_model = meta_class.GravityModel()
meta_obj = meta_class.Cfd(gravity_model=gravity_model)
# Test setting the attribute on init
self.assertEqual(meta_obj.gravity_model, gravity_model)
self.assertEqual(
meta_obj.definition,
'Computational fluid dynamics general (set of ) equations for momentum, mass and energy'
) # noqa
# Test the default values
self.assertIsInstance(meta_obj.compressibility_model,
meta_class.IncompressibleFluidModel)
self.assertIsInstance(meta_obj.thermal_model,
meta_class.IsothermalModel)
self.assertIsInstance(meta_obj.turbulence_model,
meta_class.LaminarFlowModel)
self.assertIsInstance(meta_obj.multiphase_model,
meta_class.SinglePhaseModel)
self.assertIsInstance(meta_obj.rheology_model,
meta_class.NewtonianFluidModel)
self.assertIsInstance(meta_obj.electrostatic_model,
meta_class.ConstantElectrostaticFieldModel)
self.check_cuds_item(meta_obj)
self.check_cuds_component(meta_obj)
self.check_model_equation(meta_obj)
def test_assign_vector(self):
'''Test for assigning value to a CUBA with vector+float type'''
# GravityModel.acceleration is a vector of float
gravity_model = meta_class.GravityModel()
with self.assertRaises(ValueError):
# shape should be (3)
gravity_model.acceleration = (1.0, 2.0)
# Make sure values assigned can be obtained within
# numerical precision
expected = (1.e-10, 1.e-10, 1.e-10)
gravity_model.acceleration = expected
actual = gravity_model.acceleration
numpy.testing.assert_allclose(actual, expected)
bp.material = cuds.get_by_name('sludge').uid
# water is Newtonian so no need to add that to rheology models
cfd.rheology_model = bp
cuds.add([cfd])
sm = api.StandardStressModel(name='standard_stress_model')
cuds.add([sm])
rvm = api.SimpleRelativeVelocityModel(name='simple_rel_vel_model')
rvm.diffusion_velocity = (0, -0.002, 0)
rvm.linear_constant = 285.0
cuds.add([rvm])
gm = api.GravityModel(name='gravitation')
gm.acceleration = (0, -9.81, 0)
cuds.add([gm])
# boundary condtitions
vel_inlet = api.ConstantVelocityCondition((0.0191, 0, 0), water,
name='vel_inlet')
pres_inlet = api.ZeroGradientPressureCondition(0.0, water, name='pres_inlet')
vf_inlet = api.ConstantVolumeFractionCondition(0.001, water, name='vf_inlet')
vel_outlet = api.InletOutletVelocity((0, 0, 0), water, name='vel_outlet')
pres_outlet = api.ConstantPressureCondition(0.0, water, name='pres_outlet')
vf_outlet = api.InletOutletVolumeFraction(0.001, water, name='vf_outlet')
vel_walls = api.ConstantVelocityCondition((0, 0, 0), water, name='vel_walls')