Merge branch 'FixParamterAllenCahn' into 'master'

Correction for paramter calculation of Allen Cahn model

See merge request pycodegen/lbmpy!178

src/lbmpy/phasefield_allen_cahn/parameter_calculation.py

@@ -366,8 +366,8 @@ def calculate_parameters_taylor_bubble(reference_length=128,
     water_mu_lattice_units = math.sqrt(density_diff * density_heavy * grav_df_cube) / inverse_viscosity_number
     air_mu_lattice_units = water_mu_lattice_units / (water_mu / air_mu)
 
-    dynamic_viscosity_heavy = water_mu_lattice_units / density_heavy
-    dynamic_viscosity_light = air_mu_lattice_units / density_light
+    dynamic_viscosity_heavy = water_mu_lattice_units
+    dynamic_viscosity_light = air_mu_lattice_units
 
     surface_tension_lattice_units = density_diff * gravity_lattice_units * diameter_fluid ** 2 / bond_number
 

tests/phasefield_allen_cahn/test_analytical.py

@@ -44,6 +44,6 @@ def test_analytical():
     assert np.isclose(parameters.density_heavy, 1.0)
     assert np.isclose(parameters.density_light, 0.001207114228456914)
     assert np.isclose(parameters.dynamic_viscosity_heavy, 5.733727652152216e-05)
-    assert np.isclose(parameters.dynamic_viscosity_light, 0.0008630017037694861)
+    assert np.isclose(parameters.dynamic_viscosity_light, 1.0417416358027054e-06)
     assert np.isclose(parameters.gravitational_acceleration, -7.407407407407407e-08)
     assert np.isclose(parameters.surface_tension, 3.149857262258028e-05, rtol=1e-05)