diff --git a/src/lagrangian/intermediate/submodels/Thermodynamic/HeatTransferModel/RanzMarshall/RanzMarshall.C b/src/lagrangian/intermediate/submodels/Thermodynamic/HeatTransferModel/RanzMarshall/RanzMarshall.C
index c45e0968fa74bdb0291796e58fd6c253e9fe11b3..9a3015164b67834ff02bcb98fbbae600a5b5c33f 100644
--- a/src/lagrangian/intermediate/submodels/Thermodynamic/HeatTransferModel/RanzMarshall/RanzMarshall.C
+++ b/src/lagrangian/intermediate/submodels/Thermodynamic/HeatTransferModel/RanzMarshall/RanzMarshall.C
@@ -6,6 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2011 OpenFOAM Foundation
+    Copyright (C) 2021 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -36,14 +37,22 @@ Foam::RanzMarshall<CloudType>::RanzMarshall
     CloudType& cloud
 )
 :
-    HeatTransferModel<CloudType>(dict, cloud, typeName)
+    HeatTransferModel<CloudType>(dict, cloud, typeName),
+    a_(this->coeffDict().template getOrDefault<scalar>("a", 2.0)),
+    b_(this->coeffDict().template getOrDefault<scalar>("b", 0.6)),
+    m_(this->coeffDict().template getOrDefault<scalar>("m", 1.0/2.0)),
+    n_(this->coeffDict().template getOrDefault<scalar>("n", 1.0/3.0))
 {}
 
 
 template<class CloudType>
 Foam::RanzMarshall<CloudType>::RanzMarshall(const RanzMarshall<CloudType>& htm)
 :
-    HeatTransferModel<CloudType>(htm)
+    HeatTransferModel<CloudType>(htm),
+    a_(htm.a_),
+    b_(htm.b_),
+    m_(htm.m_),
+    n_(htm.n_)
 {}
 
 
@@ -63,7 +72,7 @@ Foam::scalar Foam::RanzMarshall<CloudType>::Nu
     const scalar Pr
 ) const
 {
-    return 2.0 + 0.6*sqrt(Re)*cbrt(Pr);
+    return a_ + b_*pow(Re, m_)*pow(Pr, n_);
 }
 
 
diff --git a/src/lagrangian/intermediate/submodels/Thermodynamic/HeatTransferModel/RanzMarshall/RanzMarshall.H b/src/lagrangian/intermediate/submodels/Thermodynamic/HeatTransferModel/RanzMarshall/RanzMarshall.H
index 6377a3352702c9d76c534d57f9cc58a84bb2d532..7cdaf14ca9ff2723cc44735ae7d58943fe7a6d83 100644
--- a/src/lagrangian/intermediate/submodels/Thermodynamic/HeatTransferModel/RanzMarshall/RanzMarshall.H
+++ b/src/lagrangian/intermediate/submodels/Thermodynamic/HeatTransferModel/RanzMarshall/RanzMarshall.H
@@ -6,6 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2011-2016 OpenFOAM Foundation
+    Copyright (C) 2021 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -53,6 +54,20 @@ class RanzMarshall
 :
     public HeatTransferModel<CloudType>
 {
+    // Private Data
+
+        //- Correlation coefficient
+        const scalar a_;
+
+        //- Correlation coefficient
+        const scalar b_;
+
+        //- Correlation exponent of particle Reynolds number
+        const scalar m_;
+
+        //- Correlation exponent of Prandtl number
+        const scalar n_;
+
 
 public: