diff --git a/src/lagrangian/parcel/phaseProperties/phaseProperties/phaseProperties.C b/src/lagrangian/parcel/phaseProperties/phaseProperties/phaseProperties.C
index 146595dbf4..b2fd9da731 100644
--- a/src/lagrangian/parcel/phaseProperties/phaseProperties/phaseProperties.C
+++ b/src/lagrangian/parcel/phaseProperties/phaseProperties/phaseProperties.C
@@ -71,6 +71,7 @@ void Foam::phaseProperties::reorder(const wordList& specieNames)
     {
         Y_.setSize(names_.size());
         Y_ = 0;
+        carrierIds_.setSize(names_.size(), -1);
     }
 
     // Set the mass-fraction for each specie in the list to the corresponding
diff --git a/src/lagrangian/parcel/submodels/Momentum/ParticleForces/Scaled/ScaledForce.C b/src/lagrangian/parcel/submodels/Momentum/ParticleForces/Scaled/ScaledForce.C
index abf329fb75..3090c7c142 100644
--- a/src/lagrangian/parcel/submodels/Momentum/ParticleForces/Scaled/ScaledForce.C
+++ b/src/lagrangian/parcel/submodels/Momentum/ParticleForces/Scaled/ScaledForce.C
@@ -72,7 +72,7 @@ Foam::ScaledForce<CloudType>::ScaledForce
 )
 :
     ParticleForce<CloudType>(df),
-    model_(nullptr),
+    model_(df.model_().clone().ptr()),
     factor_(1)
 {}
 
diff --git a/src/thermophysicalModels/thermophysicalProperties/liquidProperties/liquid/liquid.C b/src/thermophysicalModels/thermophysicalProperties/liquidProperties/liquid/liquid.C
index c1f6efd12d..b23c16e396 100644
--- a/src/thermophysicalModels/thermophysicalProperties/liquidProperties/liquid/liquid.C
+++ b/src/thermophysicalModels/thermophysicalProperties/liquidProperties/liquid/liquid.C
@@ -83,6 +83,27 @@ Foam::liquid::liquid(const dictionary& dict)
 {}
 
 
+Foam::liquid::liquid(const liquid& lm)
+:
+    liquidProperties(lm),
+    rho_(lm.rho_().clone().ptr()),
+    pv_(lm.pv_().clone().ptr()),
+    hl_(lm.hl_().clone().ptr()),
+    Cp_(lm.Cp_().clone().ptr()),
+    h_(lm.h_().clone().ptr()),
+    Cpg_(lm.Cpg_().clone().ptr()),
+    B_(lm.B_().clone().ptr()),
+    mu_(lm.mu_().clone().ptr()),
+    mug_(lm.mug_().clone().ptr()),
+    kappa_(lm.kappa_().clone().ptr()),
+    kappag_(lm.kappag_().clone().ptr()),
+    sigma_(lm.sigma_().clone().ptr()),
+    D_(lm.D_().clone().ptr()),
+    Hf_(lm.Hf_)
+{}
+
+
+
 // * * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * //
 
 void Foam::liquid::write(Ostream& os) const
diff --git a/src/thermophysicalModels/thermophysicalProperties/liquidProperties/liquid/liquid.H b/src/thermophysicalModels/thermophysicalProperties/liquidProperties/liquid/liquid.H
index eb089a480f..cfdc4887bc 100644
--- a/src/thermophysicalModels/thermophysicalProperties/liquidProperties/liquid/liquid.H
+++ b/src/thermophysicalModels/thermophysicalProperties/liquidProperties/liquid/liquid.H
@@ -94,6 +94,9 @@ public:
         //- Construct from dictionary
         liquid(const dictionary& dict);
 
+        //- Copy constructor
+        liquid(const liquid& lm);
+
         //- Construct and return clone
         virtual autoPtr<liquidProperties> clone() const
         {