diff --git a/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloud.H b/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloud.H
index 675fc61ce1..19f584337f 100644
--- a/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloud.H
+++ b/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloud.H
@@ -196,6 +196,9 @@ public:
                 inline const BreakupModel<SprayCloud<CloudType> >&
                     breakup() const;
 
+                //- Return non-const-access to the breakup model
+                inline BreakupModel<SprayCloud<CloudType> >& breakup();
+
                 //- Return const-access to the breakup model
                 inline const StochasticCollisionModel<SprayCloud<CloudType> >&
                     stochasticCollision() const;
diff --git a/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloudI.H b/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloudI.H
index f8d6d5f118..c27248aaae 100644
--- a/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloudI.H
+++ b/src/lagrangian/spray/clouds/Templates/SprayCloud/SprayCloudI.H
@@ -49,6 +49,14 @@ Foam::SprayCloud<CloudType>::breakup() const
 }
 
 
+template<class CloudType>
+inline Foam::BreakupModel<Foam::SprayCloud<CloudType> >&
+Foam::SprayCloud<CloudType>::breakup()
+{
+    return breakupModel_();
+}
+
+
 template<class CloudType>
 inline const Foam::StochasticCollisionModel<Foam::SprayCloud<CloudType> >&
 Foam::SprayCloud<CloudType>::stochasticCollision() const
@@ -63,3 +71,5 @@ inline Foam::scalar Foam::SprayCloud<CloudType>::averageParcelMass() const
     return averageParcelMass_;
 }
 
+
+// ************************************************************************* //