diff --git a/applications/solvers/combustion/rhoReactingFoam/YEqn.H b/applications/solvers/combustion/rhoReactingFoam/YEqn.H
index 3a371f035c6ad2c6c64e6f09cf99b0d633eb9734..87b25079ff6dc5b31f09a60fbab9345d56aa2bc5 100644
--- a/applications/solvers/combustion/rhoReactingFoam/YEqn.H
+++ b/applications/solvers/combustion/rhoReactingFoam/YEqn.H
@@ -30,6 +30,9 @@ tmp<fv::convectionScheme<scalar> > mvConvection
                 combustion->R(Yi)
             );
 
+            YiEqn.relax();
+            YiEqn.solve(mesh.solver("Yi"));
+
             Yi.max(0.0);
             Yt += Yi;
         }