def writeTurbulence(self, dh, u=None):
l = 0.07 * dh
if not u:
t = Turbulence.initFieldFoam(self, self.uInf, l)
else:
t = Turbulence.initFieldFoam(self, u, l)
t.write()
def testInitField(self):
"""
Tests the conversion from velocity to Froude number
"""
for input, output in self.uInf_kEpsilonOmega:
result = Turbulence.initField(input['u'], input['l'], I=input['I'])
self.assertAlmostEqual(output['k'], result['k'])
self.assertAlmostEqual(output['epsilon'], result['epsilon'])
self.assertAlmostEqual(output['omega'], result['omega'])
def testInitField(self):
"""
Tests the conversion from velocity to Froude number
"""
for input, output in self.uInf_kEpsilonOmega:
result = Turbulence.initField(
input['u'],
input['l'],
I = input['I']
)
self.assertAlmostEqual(output['k'], result['k'])
self.assertAlmostEqual(output['epsilon'], result['epsilon'])
self.assertAlmostEqual(output['omega'], result['omega'])
