From 7f413189494bb410730d72533ad6ef277a56403a Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Fri, 12 Jul 2013 13:12:57 +0100
Subject: [PATCH] ENH: InjectionModel - stabilised volume fraction calculation

---
 .../Kinematic/InjectionModel/InjectionModel/InjectionModel.C  | 4 +++-
 1 file changed, 3 insertions(+), 1 deletion(-)

diff --git a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C
index 94eb19932e..f49e790867 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C
+++ b/src/lagrangian/intermediate/submodels/Kinematic/InjectionModel/InjectionModel/InjectionModel.C
@@ -72,7 +72,9 @@ bool Foam::InjectionModel<CloudType>::prepareForNextTimeStep
     newParcels = this->parcelsToInject(t0, t1);
 
     // Volume of parcels to inject
-    newVolumeFraction = this->volumeToInject(t0, t1)/volumeTotal_;
+    newVolumeFraction =
+        this->volumeToInject(t0, t1)
+       /(volumeTotal_ + ROOTVSMALL);
 
     if (newVolumeFraction > 0)
     {