diff --git a/src/TurbulenceModels/turbulenceModels/LES/dynamicLagrangian/dynamicLagrangian.C b/src/TurbulenceModels/turbulenceModels/LES/dynamicLagrangian/dynamicLagrangian.C
index 6a6fb02b93..c0267708f5 100644
--- a/src/TurbulenceModels/turbulenceModels/LES/dynamicLagrangian/dynamicLagrangian.C
+++ b/src/TurbulenceModels/turbulenceModels/LES/dynamicLagrangian/dynamicLagrangian.C
@@ -24,6 +24,7 @@ License
 \*---------------------------------------------------------------------------*/
 
 #include "dynamicLagrangian.H"
+#include "fvOptions.H"
 
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
@@ -42,6 +43,9 @@ void dynamicLagrangian<BasicTurbulenceModel>::correctNut
 {
     this->nut_ = (flm_/fmm_)*sqr(this->delta())*mag(dev(symm(gradU)));
     this->nut_.correctBoundaryConditions();
+    fv::options::New(this->mesh_).correct(this->nut_);
+
+    BasicTurbulenceModel::correctNut();
 }
 
 
@@ -157,6 +161,7 @@ void dynamicLagrangian<BasicTurbulenceModel>::correct()
     // Local references
     const surfaceScalarField& phi = this->phi_;
     const volVectorField& U = this->U_;
+    fv::options& fvOptions(fv::options::New(this->mesh_));
 
     LESeddyViscosity<BasicTurbulenceModel>::correct();
 
@@ -190,11 +195,13 @@ void dynamicLagrangian<BasicTurbulenceModel>::correct()
      ==
         invT*LM
       - fvm::Sp(invT, flm_)
+      + fvOptions(flm_)
     );
 
     flmEqn.relax();
+    fvOptions.constrain(flmEqn);
     flmEqn.solve();
-
+    fvOptions.correct(flm_);
     bound(flm_, flm0_);
 
     volScalarField MM(M && M);
@@ -206,11 +213,13 @@ void dynamicLagrangian<BasicTurbulenceModel>::correct()
      ==
         invT*MM
       - fvm::Sp(invT, fmm_)
+      + fvOptions(fmm_)
     );
 
     fmmEqn.relax();
+    fvOptions.constrain(fmmEqn);
     fmmEqn.solve();
-
+    fvOptions.correct(fmm_);
     bound(fmm_, fmm0_);
 
     correctNut(gradU);