diff --git a/applications/solvers/multiphase/compressibleInterFoam/Allwclean b/applications/solvers/multiphase/compressibleInterFoam/Allwclean
index 682ca993f6a6b6db68fca2adbbce62ed93bf6471..aa0f3c01c6ba74c8d0cd77ee40e09214002110f6 100755
--- a/applications/solvers/multiphase/compressibleInterFoam/Allwclean
+++ b/applications/solvers/multiphase/compressibleInterFoam/Allwclean
@@ -2,7 +2,7 @@
 cd ${0%/*} || exit 1    # run from this directory
 set -x
 
-wclean libso twoPhaseThermo
+wclean libso twoPhaseMixtureThermo
 wclean
 wclean compressibleInterDyMFoam
 
diff --git a/applications/solvers/multiphase/compressibleInterFoam/Allwmake b/applications/solvers/multiphase/compressibleInterFoam/Allwmake
index 9ff44117f219c28956a5165cd5d58c2ba5ee1916..a93da26cddb8cbeb7421c7bf086ef7ef86241fc4 100755
--- a/applications/solvers/multiphase/compressibleInterFoam/Allwmake
+++ b/applications/solvers/multiphase/compressibleInterFoam/Allwmake
@@ -2,7 +2,7 @@
 cd ${0%/*} || exit 1    # run from this directory
 set -x
 
-wmake libso twoPhaseThermo
+wmake libso twoPhaseMixtureThermo
 wmake
 wmake compressibleInterDyMFoam
 
diff --git a/applications/solvers/multiphase/compressibleInterFoam/Make/options b/applications/solvers/multiphase/compressibleInterFoam/Make/options
index a5738d03dcaf0209f98a7da780a84d9c853ebe1f..147bf40e5b33b4ad18c248fa683c027547fbc220 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/Make/options
+++ b/applications/solvers/multiphase/compressibleInterFoam/Make/options
@@ -1,5 +1,5 @@
 EXE_INC = \
-    -ItwoPhaseThermo \
+    -ItwoPhaseMixtureThermo \
     -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
     -I$(LIB_SRC)/transportModels/twoPhaseProperties/lnInclude \
     -I$(LIB_SRC)/transportModels/interfaceProperties/lnInclude \
@@ -7,7 +7,7 @@ EXE_INC = \
     -I$(LIB_SRC)/finiteVolume/lnInclude
 
 EXE_LIBS = \
-    -ltwoPhaseThermo \
+    -ltwoPhaseMixtureThermo \
     -lfluidThermophysicalModels \
     -lspecie \
     -ltwoPhaseProperties \
diff --git a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/Make/options b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/Make/options
index 661a30dc1847c3ecb37da7ea7fc70a6e7b20afd3..3ee1e97f92e640bb2a4dc21ba8034d646db71f34 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/Make/options
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/Make/options
@@ -1,6 +1,6 @@
 EXE_INC = \
     -I.. \
-    -I../twoPhaseThermo \
+    -I../twoPhaseMixtureThermo \
     -I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
     -I$(LIB_SRC)/transportModels/twoPhaseProperties/lnInclude \
     -I$(LIB_SRC)/transportModels/interfaceProperties/lnInclude \
@@ -11,7 +11,7 @@ EXE_INC = \
     -I$(LIB_SRC)/finiteVolume/lnInclude
 
 EXE_LIBS = \
-    -ltwoPhaseThermo \
+    -ltwoPhaseMixtureThermo \
     -lfluidThermophysicalModels \
     -lspecie \
     -ltwoPhaseProperties \
diff --git a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
index 3cb787d144ef9369a5ded2416de6469ba43c3333..27d3b87cd1d676f267f5bf811ea297dc3a7cf2e2 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterDyMFoam/compressibleInterDyMFoam.C
@@ -43,7 +43,7 @@ Description
 #include "subCycle.H"
 #include "interfaceProperties.H"
 #include "twoPhaseMixture.H"
-#include "twoPhaseThermo.H"
+#include "twoPhaseMixtureThermo.H"
 #include "turbulenceModel.H"
 #include "pimpleControl.H"
 
diff --git a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
index 85cab52cdb22999435645e3442e4afcaa2744fd8..3148382ad2e23cb2ddc774dd7064b8bc28ce8c38 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
+++ b/applications/solvers/multiphase/compressibleInterFoam/compressibleInterFoam.C
@@ -41,7 +41,7 @@ Description
 #include "rhoThermo.H"
 #include "interfaceProperties.H"
 #include "twoPhaseMixture.H"
-#include "twoPhaseThermo.H"
+#include "twoPhaseMixtureThermo.H"
 #include "turbulenceModel.H"
 #include "pimpleControl.H"
 
diff --git a/applications/solvers/multiphase/compressibleInterFoam/createFields.H b/applications/solvers/multiphase/compressibleInterFoam/createFields.H
index a650b435efb0e418f8543714550ff9b3f89f08ff..102d4d57a99a8bfa111eaca227399c27d2066935 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/createFields.H
+++ b/applications/solvers/multiphase/compressibleInterFoam/createFields.H
@@ -28,8 +28,8 @@
 
     #include "createPhi.H"
 
-    Info<< "Constructing twoPhaseThermo\n" << endl;
-    twoPhaseThermo twoPhaseProperties(mesh);
+    Info<< "Constructing twoPhaseMixtureThermo\n" << endl;
+    twoPhaseMixtureThermo twoPhaseProperties(mesh);
 
     volScalarField& alpha1(twoPhaseProperties.alpha1());
     volScalarField& alpha2(twoPhaseProperties.alpha2());
diff --git a/applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/Make/files b/applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/Make/files
new file mode 100644
index 0000000000000000000000000000000000000000..1cd0a09577c241265ce500c16b0df145a4aa3b2c
--- /dev/null
+++ b/applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/Make/files
@@ -0,0 +1,3 @@
+twoPhaseMixtureThermo.C
+
+LIB = $(FOAM_LIBBIN)/libtwoPhaseMixtureThermo
diff --git a/applications/solvers/multiphase/compressibleInterFoam/twoPhaseThermo/Make/options b/applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/Make/options
similarity index 100%
rename from applications/solvers/multiphase/compressibleInterFoam/twoPhaseThermo/Make/options
rename to applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/Make/options
diff --git a/applications/solvers/multiphase/compressibleInterFoam/twoPhaseThermo/twoPhaseThermo.C b/applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/twoPhaseMixtureThermo.C
similarity index 75%
rename from applications/solvers/multiphase/compressibleInterFoam/twoPhaseThermo/twoPhaseThermo.C
rename to applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/twoPhaseMixtureThermo.C
index e1842b88fd7321947810e9c7b3c98cb805ce7f0f..c48990a721406dc81f55e236e95aec01fce1f7a3 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/twoPhaseThermo/twoPhaseThermo.C
+++ b/applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/twoPhaseMixtureThermo.C
@@ -23,7 +23,7 @@ License
 
 \*---------------------------------------------------------------------------*/
 
-#include "twoPhaseThermo.H"
+#include "twoPhaseMixtureThermo.H"
 #include "gradientEnergyFvPatchScalarField.H"
 #include "mixedEnergyFvPatchScalarField.H"
 
@@ -32,12 +32,12 @@ License
 
 namespace Foam
 {
-    defineTypeNameAndDebug(twoPhaseThermo, 0);
+    defineTypeNameAndDebug(twoPhaseMixtureThermo, 0);
 }
 
 // * * * * * * * * * * * * * * * * Constructors  * * * * * * * * * * * * * * //
 
-Foam::twoPhaseThermo::twoPhaseThermo
+Foam::twoPhaseMixtureThermo::twoPhaseMixtureThermo
 (
     const fvMesh& mesh
 )
@@ -56,13 +56,13 @@ Foam::twoPhaseThermo::twoPhaseThermo
 
 // * * * * * * * * * * * * * * * * Destructor  * * * * * * * * * * * * * * * //
 
-Foam::twoPhaseThermo::~twoPhaseThermo()
+Foam::twoPhaseMixtureThermo::~twoPhaseMixtureThermo()
 {}
 
 
 // * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //
 
-void Foam::twoPhaseThermo::correct()
+void Foam::twoPhaseMixtureThermo::correct()
 {
     thermo1_->he() = thermo1_->he(p_, T_);
     thermo1_->correct();
@@ -76,19 +76,19 @@ void Foam::twoPhaseThermo::correct()
 }
 
 
-bool Foam::twoPhaseThermo::incompressible() const
+bool Foam::twoPhaseMixtureThermo::incompressible() const
 {
     return thermo1_->incompressible() && thermo2_->incompressible();
 }
 
 
-bool Foam::twoPhaseThermo::isochoric() const
+bool Foam::twoPhaseMixtureThermo::isochoric() const
 {
     return thermo1_->isochoric() && thermo2_->isochoric();
 }
 
 
-Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::he
+Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::he
 (
     const volScalarField& p,
     const volScalarField& T
@@ -98,7 +98,7 @@ Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::he
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::he
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::he
 (
     const scalarField& p,
     const scalarField& T,
@@ -111,7 +111,7 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::he
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::he
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::he
 (
     const scalarField& p,
     const scalarField& T,
@@ -124,13 +124,13 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::he
 }
 
 
-Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::hc() const
+Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::hc() const
 {
     return alpha1()*thermo1_->hc() + alpha2()*thermo2_->hc();
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::THE
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::THE
 (
     const scalarField& h,
     const scalarField& p,
@@ -138,12 +138,12 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::THE
     const labelList& cells
 ) const
 {
-    notImplemented("twoPhaseThermo::THE(...)");
+    notImplemented("twoPhaseMixtureThermo::THE(...)");
     return T0;
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::THE
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::THE
 (
     const scalarField& h,
     const scalarField& p,
@@ -151,18 +151,18 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::THE
     const label patchi
 ) const
 {
-    notImplemented("twoPhaseThermo::THE(...)");
+    notImplemented("twoPhaseMixtureThermo::THE(...)");
     return T0;
 }
 
 
-Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::Cp() const
+Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::Cp() const
 {
     return alpha1()*thermo1_->Cp() + alpha2()*thermo2_->Cp();
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::Cp
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::Cp
 (
     const scalarField& p,
     const scalarField& T,
@@ -175,13 +175,13 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::Cp
 }
 
 
-Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::Cv() const
+Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::Cv() const
 {
     return alpha1()*thermo1_->Cv() + alpha2()*thermo2_->Cv();
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::Cv
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::Cv
 (
     const scalarField& p,
     const scalarField& T,
@@ -194,13 +194,13 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::Cv
 }
 
 
-Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::gamma() const
+Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::gamma() const
 {
     return alpha1()*thermo1_->gamma() + alpha2()*thermo2_->gamma();
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::gamma
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::gamma
 (
     const scalarField& p,
     const scalarField& T,
@@ -213,13 +213,13 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::gamma
 }
 
 
-Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::Cpv() const
+Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::Cpv() const
 {
     return alpha1()*thermo1_->Cpv() + alpha2()*thermo2_->Cpv();
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::Cpv
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::Cpv
 (
     const scalarField& p,
     const scalarField& T,
@@ -232,7 +232,7 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::Cpv
 }
 
 
-Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::CpByCpv() const
+Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::CpByCpv() const
 {
     return
         alpha1()*thermo1_->CpByCpv()
@@ -240,7 +240,7 @@ Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::CpByCpv() const
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::CpByCpv
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::CpByCpv
 (
     const scalarField& p,
     const scalarField& T,
@@ -253,13 +253,13 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::CpByCpv
 }
 
 
-Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::kappa() const
+Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::kappa() const
 {
     return alpha1()*thermo1_->kappa() + alpha2()*thermo2_->kappa();
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::kappa
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::kappa
 (
     const label patchi
 ) const
@@ -270,7 +270,7 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::kappa
 }
 
 
-Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::kappaEff
+Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::kappaEff
 (
     const volScalarField& alphat
 ) const
@@ -281,7 +281,7 @@ Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::kappaEff
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::kappaEff
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::kappaEff
 (
     const scalarField& alphat,
     const label patchi
@@ -294,7 +294,7 @@ Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::kappaEff
 }
 
 
-Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::alphaEff
+Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::alphaEff
 (
     const volScalarField& alphat
 ) const
@@ -305,7 +305,7 @@ Foam::tmp<Foam::volScalarField> Foam::twoPhaseThermo::alphaEff
 }
 
 
-Foam::tmp<Foam::scalarField> Foam::twoPhaseThermo::alphaEff
+Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::alphaEff
 (
     const scalarField& alphat,
     const label patchi
diff --git a/applications/solvers/multiphase/compressibleInterFoam/twoPhaseThermo/twoPhaseThermo.H b/applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/twoPhaseMixtureThermo.H
similarity index 93%
rename from applications/solvers/multiphase/compressibleInterFoam/twoPhaseThermo/twoPhaseThermo.H
rename to applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/twoPhaseMixtureThermo.H
index 24966273c65df29d420e55243dd575ea978fefb6..b2e636a4c471ac92b490aa58f1084b72c6dc3888 100644
--- a/applications/solvers/multiphase/compressibleInterFoam/twoPhaseThermo/twoPhaseThermo.H
+++ b/applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/twoPhaseMixtureThermo.H
@@ -22,19 +22,19 @@ License
     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
 Class
-    Foam::twoPhaseThermo
+    Foam::twoPhaseMixtureThermo
 
 Description
 
 SourceFiles
-    twoPhaseThermoI.H
-    twoPhaseThermo.C
-    twoPhaseThermoIO.C
+    twoPhaseMixtureThermoI.H
+    twoPhaseMixtureThermo.C
+    twoPhaseMixtureThermoIO.C
 
 \*---------------------------------------------------------------------------*/
 
-#ifndef twoPhaseThermo_H
-#define twoPhaseThermo_H
+#ifndef twoPhaseMixtureThermo_H
+#define twoPhaseMixtureThermo_H
 
 #include "rhoThermo.H"
 #include "psiThermo.H"
@@ -46,10 +46,10 @@ namespace Foam
 {
 
 /*---------------------------------------------------------------------------*\
-                         Class twoPhaseThermo Declaration
+                         Class twoPhaseMixtureThermo Declaration
 \*---------------------------------------------------------------------------*/
 
-class twoPhaseThermo
+class twoPhaseMixtureThermo
 :
     public psiThermo,
     public twoPhaseMixture
@@ -66,20 +66,20 @@ class twoPhaseThermo
 public:
 
     //- Runtime type information
-    TypeName("twoPhaseThermo");
+    TypeName("twoPhaseMixtureThermo");
 
 
     // Constructors
 
         //- Construct from mesh
-        twoPhaseThermo
+        twoPhaseMixtureThermo
         (
             const fvMesh& mesh
         );
 
 
     //- Destructor
-    virtual ~twoPhaseThermo();
+    virtual ~twoPhaseMixtureThermo();
 
 
     // Member Functions
@@ -112,14 +112,14 @@ public:
             //  Non-const access allowed for transport equations
             virtual volScalarField& he()
             {
-                notImplemented("twoPhaseThermo::he()");
+                notImplemented("twoPhaseMixtureThermo::he()");
                 return thermo1_->he();
             }
 
             //- Enthalpy/Internal energy [J/kg]
             virtual const volScalarField& he() const
             {
-                notImplemented("twoPhaseThermo::he() const");
+                notImplemented("twoPhaseMixtureThermo::he() const");
                 return thermo1_->he();
             }
 
diff --git a/applications/solvers/multiphase/compressibleInterFoam/twoPhaseThermo/Make/files b/applications/solvers/multiphase/compressibleInterFoam/twoPhaseThermo/Make/files
deleted file mode 100644
index 64771e2f6ea0a4eadf36fc495f8e38129e289a9b..0000000000000000000000000000000000000000
--- a/applications/solvers/multiphase/compressibleInterFoam/twoPhaseThermo/Make/files
+++ /dev/null
@@ -1,3 +0,0 @@
-twoPhaseThermo.C
-
-LIB = $(FOAM_LIBBIN)/libtwoPhaseThermo