diff --git a/support/Pipelines/Bbh/ControlId.py b/support/Pipelines/Bbh/ControlId.py
index 1fee4da015f66..f6d46cdf25675 100644
--- a/support/Pipelines/Bbh/ControlId.py
+++ b/support/Pipelines/Bbh/ControlId.py
@@ -298,7 +298,8 @@ def Residual(u):
         F = Residual(u)
 
         # Update the Jacobian using Broyden's method
-        J += np.outer(F, Delta_u) / np.dot(Delta_u, Delta_u)
+        damping_factor = 1.0 - np.exp(-0.5 * iteration)
+        J += np.outer(F, Delta_u) / np.dot(Delta_u, Delta_u) * damping_factor
 
     data_file.close()