diff --git a/src/lbmpy/phasefield_allen_cahn/parameter_calculation.py b/src/lbmpy/phasefield_allen_cahn/parameter_calculation.py index ead4c4b95eec2bd194d7be78e6bea71f681c7ef7..d4077a99030f1dbf66196e91ddfbc4b73564f15b 100644 --- a/src/lbmpy/phasefield_allen_cahn/parameter_calculation.py +++ b/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 diff --git a/tests/phasefield_allen_cahn/test_analytical.py b/tests/phasefield_allen_cahn/test_analytical.py index a9f4d2ee1af9e35a8d344f6f075181bc0df18b53..ec187936f6952aff61a29353fc7996dc3d14bae0 100644 --- a/tests/phasefield_allen_cahn/test_analytical.py +++ b/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)