diff --git a/libraries/lagrangian/spray/parcels/derived/basicSprayParcel/makeBasicSprayParcelSubmodels.C b/libraries/lagrangian/spray/parcels/derived/basicSprayParcel/makeBasicSprayParcelSubmodels.C
index 01a480b..893319a 100644
--- a/libraries/lagrangian/spray/parcels/derived/basicSprayParcel/makeBasicSprayParcelSubmodels.C
+++ b/libraries/lagrangian/spray/parcels/derived/basicSprayParcel/makeBasicSprayParcelSubmodels.C
@@ -42,6 +42,7 @@ License
#include "makeReactingParcelSurfaceFilmModels.H"
// Spray
+#include "DistortedSphereDragForce.H"
#include "makeSprayParcelAtomizationModels.H"
#include "makeSprayParcelBreakupModels.H"
#include "makeSprayParcelCollisionModels.H"
@@ -67,6 +68,7 @@ namespace Foam
makeReactingParcelSurfaceFilmModels(basicSprayCloud);
// Spray sub-models
+ makeParticleForceModelType(DistortedSphereDragForce, basicSprayCloud);
makeSprayParcelAtomizationModels(basicSprayCloud);
makeSprayParcelBreakupModels(basicSprayCloud);
makeSprayParcelCollisionModels(basicSprayCloud);
diff --git a/libraries/lagrangian/spray/submodels/ParticleForces/Drag/DistortedSphereDrag/DistortedSphereDragForce.C b/libraries/lagrangian/spray/submodels/ParticleForces/Drag/DistortedSphereDrag/DistortedSphereDragForce.C
new file mode 100644
index 0000000..5993e7b
--- /dev/null
+++ b/libraries/lagrangian/spray/submodels/ParticleForces/Drag/DistortedSphereDrag/DistortedSphereDragForce.C
@@ -0,0 +1,98 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see .
+
+\*---------------------------------------------------------------------------*/
+
+#include "DistortedSphereDragForce.H"
+
+// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
+
+template
+Foam::scalar Foam::DistortedSphereDragForce::CdRe(const scalar Re) const
+{
+ if (Re > 1000.0)
+ {
+ return 0.424*Re;
+ }
+ else
+ {
+ return 24.0*(1.0 + 1.0/6.0*pow(Re, 2.0/3.0));
+ }
+}
+
+
+// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
+
+template
+Foam::DistortedSphereDragForce::DistortedSphereDragForce
+(
+ CloudType& owner,
+ const fvMesh& mesh,
+ const dictionary& dict
+)
+:
+ ParticleForce(owner, mesh, dict, typeName, false)
+{}
+
+
+template
+Foam::DistortedSphereDragForce::DistortedSphereDragForce
+(
+ const DistortedSphereDragForce& df
+)
+:
+ ParticleForce(df)
+{}
+
+
+// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
+
+template
+Foam::DistortedSphereDragForce::~DistortedSphereDragForce()
+{}
+
+
+// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
+
+template
+Foam::forceSuSp Foam::DistortedSphereDragForce::calcCoupled
+(
+ const typename CloudType::parcelType& p,
+ const scalar dt,
+ const scalar mass,
+ const scalar Re,
+ const scalar muc
+) const
+{
+ forceSuSp value(vector::zero, 0.0);
+
+ // limit the drop distortion to y=0 (sphere) and y=1 (disk)
+ scalar y = min(max(p.y(),0.0),1.0);
+
+ value.Sp() = mass*0.75*muc*CdRe(Re)*(1+2.632*y)/(p.rho()*sqr(p.d()));
+
+ return value;
+}
+
+
+// ************************************************************************* //
diff --git a/libraries/lagrangian/spray/submodels/ParticleForces/Drag/DistortedSphereDrag/DistortedSphereDragForce.H b/libraries/lagrangian/spray/submodels/ParticleForces/Drag/DistortedSphereDrag/DistortedSphereDragForce.H
new file mode 100644
index 0000000..bf7cb47
--- /dev/null
+++ b/libraries/lagrangian/spray/submodels/ParticleForces/Drag/DistortedSphereDrag/DistortedSphereDragForce.H
@@ -0,0 +1,126 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration |
+ \\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
+ \\/ M anipulation |
+-------------------------------------------------------------------------------
+License
+ This file is part of OpenFOAM.
+
+ OpenFOAM is free software: you can redistribute it and/or modify it
+ under the terms of the GNU General Public License as published by
+ the Free Software Foundation, either version 3 of the License, or
+ (at your option) any later version.
+
+ OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
+ ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
+ FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
+ for more details.
+
+ You should have received a copy of the GNU General Public License
+ along with OpenFOAM. If not, see .
+
+Class
+ Foam::DistortedSphereDragForce
+
+Description
+ Drag model based on assumption of distorted spheres according to:
+
+ \verbatim
+ Liu, A.B.; Mather, D.; Reitz, R.D.
+ "Effects of Drop Drag and Breakup on Fuel Sprays"
+ SAE Paper 930072
+ SAE Transactions, Vol. 102, Section 3, Journal of Engines, pp. 63-95, 1993
+ \endverbatim
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef DistortedSphereDragForce_H
+#define DistortedSphereDragForce_H
+
+#include "ParticleForce.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+/*---------------------------------------------------------------------------*\
+ Class DistortedSphereDragForce Declaration
+\*---------------------------------------------------------------------------*/
+
+template
+class DistortedSphereDragForce
+:
+ public ParticleForce
+{
+ // Private Member Functions
+
+ //- Drag coefficient multiplied by Reynolds number
+ scalar CdRe(const scalar Re) const;
+
+
+public:
+
+ //- Runtime type information
+ TypeName("distortedSphereDrag");
+
+
+ // Constructors
+
+ //- Construct from mesh
+ DistortedSphereDragForce
+ (
+ CloudType& owner,
+ const fvMesh& mesh,
+ const dictionary& dict
+ );
+
+ //- Construct copy
+ DistortedSphereDragForce(const DistortedSphereDragForce& df);
+
+ //- Construct and return a clone
+ virtual autoPtr > clone() const
+ {
+ return autoPtr >
+ (
+ new DistortedSphereDragForce(*this)
+ );
+ }
+
+
+ //- Destructor
+ virtual ~DistortedSphereDragForce();
+
+
+ // Member Functions
+
+ // Evaluation
+
+ //- Calculate the coupled force
+ virtual forceSuSp calcCoupled
+ (
+ const typename CloudType::parcelType& p,
+ const scalar dt,
+ const scalar mass,
+ const scalar Re,
+ const scalar muc
+ ) const;
+};
+
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+} // End namespace Foam
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#ifdef NoRepository
+# include "DistortedSphereDragForce.C"
+#endif
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+#endif
+
+// ************************************************************************* //