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