diff --git a/src/thermophysicalModels/thermophysicalProperties/Make/files b/src/thermophysicalModels/thermophysicalProperties/Make/files
index 78b8fd9b1..ed7e4f9af 100644
--- a/src/thermophysicalModels/thermophysicalProperties/Make/files
+++ b/src/thermophysicalModels/thermophysicalProperties/Make/files
@@ -51,6 +51,7 @@ $(workDir)/liquidProperties/MB/MB.C
$(workDir)/liquidProperties/CH4N2O/CH4N2O.C
$(workDir)/liquidProperties/nC3H8O/nC3H8O.C
$(workDir)/liquidProperties/iC3H8O/iC3H8O.C
+liquidProperties/O2/O2.C
solidProperties/solidProperties/solidProperties.C
diff --git a/src/thermophysicalModels/thermophysicalProperties/liquidProperties/O2/O2.C b/src/thermophysicalModels/thermophysicalProperties/liquidProperties/O2/O2.C
new file mode 100644
index 000000000..0f253bf62
--- /dev/null
+++ b/src/thermophysicalModels/thermophysicalProperties/liquidProperties/O2/O2.C
@@ -0,0 +1,113 @@
+#include "O2.H"
+#include "addToRunTimeSelectionTable.H"
+
+
+namespace Foam
+{
+ defineTypeNameAndDebug(O2, 0);
+ addToRunTimeSelectionTable(liquidProperties, O2,);
+ addToRunTimeSelectionTable(liquidProperties, O2, dictionary);
+}
+
+
+Foam::O2::O2()
+:
+ liquidProperties
+ (
+ 0.031999,
+ 154.5994,
+ 5.0464e+6,
+ 0.07495,
+ 0.294248,
+ 54.361,
+ 146.28,
+ 90.188,
+ -1,
+ 0.0222,
+ 14710.24
+ ),
+
+ pv_(2.68794635e+01, -9.58794777e+02, -1.05052914e+00, -6.90048355e-05, 1.21846002e-02),
+
+ rho_(1.14632554e+02, 2.75655615e-01, 154.5994, 2.77505361e-01),
+
+ hl_(154.5994, 2.89780460e+05, 6.00000000e-01, -7.05482356e-01, 4.75635399e-01, 0.0),
+
+ Cp_(154.5994, 1.09590080e+01, 1.45625217e+03, 2.33843393e+01, -8.99949070e+01),
+
+ h_(-9.68908452e+05, 2.50637084e+04, -5.16097446e+02, 5.55884345e+00, -2.93205446e-02, 6.09448175e-05),
+
+ mu_(-2.11200494e+00, 7.62848791e+01, -1.61677107e+00, -7.70004913e-30, 1.31855757e+01),
+
+ sigma_(154.5994, 3.83537868e-02, 1.21069008e+00, 4.43081971e-02, -5.45942334e-02, 0.0),
+
+ kappa_(-9.19925611e-02, 1.93351101e-02, -4.51801714e-04, 4.79054225e-06, -2.47777620e-08, 5.00382658e-11),
+
+ Cpg_(154.5994, 1.36975300e+01, 3.52746624e+02, -6.80059012e+00, -4.86846476e+01),
+
+ B_(1.20601567e-02, -6.15471876e-04, 1.22024524e-05, -1.17761704e-07, 5.54479099e-10, -1.02125637e-12),
+
+ mug_(-1.76036598e+01, 1.50014787e+01, 1.23492678e+00, 1.17825779e-56, 2.53862523e+01),
+
+ kappag_(-9.52797778e-01, 5.36816017e-02, -1.17434164e-03, 1.25352582e-05, -6.52648492e-08, 1.32897143e-10),
+
+ D_(16.3, 18.5, 31.999, 28.0)
+{}
+
+
+Foam::O2::O2
+(
+ const liquidProperties& l,
+ const NSRDSfunc5& density,
+ const NSRDSfunc1& vapourPressure,
+ const NSRDSfunc6& heatOfVapourisation,
+ const NSRDSfunc14& heatCapacity,
+ const NSRDSfunc0& enthalpy,
+ const NSRDSfunc14& idealGasHeatCapacity,
+ const NSRDSfunc0& secondVirialCoeff,
+ const NSRDSfunc1& dynamicViscosity,
+ const NSRDSfunc1& vapourDynamicViscosity,
+ const NSRDSfunc0& thermalConductivity,
+ const NSRDSfunc0& vapourThermalConductivity,
+ const NSRDSfunc6& surfaceTension,
+ const APIdiffCoefFunc& vapourDiffussivity
+)
+:
+ liquidProperties(l),
+ rho_(density),
+ pv_(vapourPressure),
+ hl_(heatOfVapourisation),
+ Cp_(heatCapacity),
+ h_(enthalpy),
+ Cpg_(idealGasHeatCapacity),
+ B_(secondVirialCoeff),
+ mu_(dynamicViscosity),
+ mug_(vapourDynamicViscosity),
+ kappa_(thermalConductivity),
+ kappag_(vapourThermalConductivity),
+ sigma_(surfaceTension),
+ D_(vapourDiffussivity)
+{}
+
+
+Foam::O2::O2(const dictionary& dict)
+:
+ O2()
+{
+ readIfPresent(*this, dict);
+}
+
+
+
+void Foam::O2::writeData(Ostream& os) const
+{
+ liquidProperties::writeData(*this, os);
+}
+
+
+
+Foam::Ostream& Foam::operator<<(Ostream& os, const O2& l)
+{
+ l.writeData(os);
+ return os;
+}
\ No newline at end of file
diff --git a/src/thermophysicalModels/thermophysicalProperties/liquidProperties/O2/O2.H b/src/thermophysicalModels/thermophysicalProperties/liquidProperties/O2/O2.H
new file mode 100644
index 000000000..8228ced24
--- /dev/null
+++ b/src/thermophysicalModels/thermophysicalProperties/liquidProperties/O2/O2.H
@@ -0,0 +1,151 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration | Website: https://openfoam.org
+ \\ / A nd | Copyright (C) 2011-2019 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::O2
+
+Description
+ water
+
+SourceFiles
+ O2.C
+
+\*---------------------------------------------------------------------------*/
+
+#ifndef O2_H
+#define O2_H
+
+#include "liquidProperties.H"
+#include "NSRDSfunc0.H"
+#include "NSRDSfunc1.H"
+#include "NSRDSfunc5.H"
+#include "NSRDSfunc6.H"
+#include "NSRDSfunc14.H"
+#include "APIdiffCoefFunc.H"
+
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+namespace Foam
+{
+
+
+class O2
+:
+ public liquidProperties
+{
+
+ NSRDSfunc5 rho_;
+ NSRDSfunc1 pv_;
+ NSRDSfunc6 hl_;
+ NSRDSfunc14 Cp_;
+ NSRDSfunc0 h_;
+ NSRDSfunc14 Cpg_;
+ NSRDSfunc0 B_;
+ NSRDSfunc1 mu_;
+ NSRDSfunc1 mug_;
+ NSRDSfunc0 kappa_;
+ NSRDSfunc0 kappag_;
+ NSRDSfunc6 sigma_;
+ APIdiffCoefFunc D_;
+
+
+public:
+
+ friend class liquidProperties;
+
+ TypeName("O2");
+
+
+
+ O2();
+
+ O2
+ (
+ const liquidProperties& l,
+ const NSRDSfunc5& density,
+ const NSRDSfunc1& vapourPressure,
+ const NSRDSfunc6& heatOfVapourisation,
+ const NSRDSfunc14& heatCapacity,
+ const NSRDSfunc0& enthalpy,
+ const NSRDSfunc14& idealGasHeatCapacity,
+ const NSRDSfunc0& secondVirialCoeff,
+ const NSRDSfunc1& dynamicViscosity,
+ const NSRDSfunc1& vapourDynamicViscosity,
+ const NSRDSfunc0& thermalConductivity,
+ const NSRDSfunc0& vapourThermalConductivity,
+ const NSRDSfunc6& surfaceTension,
+ const APIdiffCoefFunc& vapourDiffussivity
+ );
+
+ O2(const dictionary& dict);
+
+ virtual autoPtr clone() const
+ {
+ return autoPtr(new O2(*this));
+ }
+
+
+ inline scalar rho(scalar p, scalar T) const;
+
+ inline scalar pv(scalar p, scalar T) const;
+
+ inline scalar hl(scalar p, scalar T) const;
+
+ inline scalar Cp(scalar p, scalar T) const;
+
+ inline scalar h(scalar p, scalar T) const;
+
+ inline scalar Cpg(scalar p, scalar T) const;
+
+ inline scalar B(scalar p, scalar T) const;
+
+ inline scalar mu(scalar p, scalar T) const;
+
+ inline scalar mug(scalar p, scalar T) const;
+
+ inline scalar kappa(scalar p, scalar T) const;
+
+ inline scalar kappag(scalar p, scalar T) const;
+
+ inline scalar sigma(scalar p, scalar T) const;
+
+ inline scalar D(scalar p, scalar T) const;
+
+ inline scalar D(scalar p, scalar T, scalar Wb) const;
+
+
+ void writeData(Ostream& os) const;
+
+ friend Ostream& operator<<(Ostream& os, const O2& l);
+};
+
+Ostream& operator<<(Ostream& os, const O2& l);
+
+}
+
+#include "O2I.H"
+
+
+#endif
+
+// ************************************************************************* //
\ No newline at end of file
diff --git a/src/thermophysicalModels/thermophysicalProperties/liquidProperties/O2/O2I.H b/src/thermophysicalModels/thermophysicalProperties/liquidProperties/O2/O2I.H
new file mode 100644
index 000000000..a7239d2ca
--- /dev/null
+++ b/src/thermophysicalModels/thermophysicalProperties/liquidProperties/O2/O2I.H
@@ -0,0 +1,110 @@
+/*---------------------------------------------------------------------------*\
+ ========= |
+ \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
+ \\ / O peration | Website: https://openfoam.org
+ \\ / A nd | Copyright (C) 2011-2018 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 .
+
+\*---------------------------------------------------------------------------*/
+
+inline Foam::scalar Foam::O2::rho(scalar p, scalar T) const
+{
+ return rho_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::pv(scalar p, scalar T) const
+{
+ return pv_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::hl(scalar p, scalar T) const
+{
+ return hl_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::Cp(scalar p, scalar T) const
+{
+ return Cp_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::h(scalar p, scalar T) const
+{
+ return h_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::Cpg(scalar p, scalar T) const
+{
+ return Cpg_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::B(scalar p, scalar T) const
+{
+ return B_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::mu(scalar p, scalar T) const
+{
+ return mu_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::mug(scalar p, scalar T) const
+{
+ return mug_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::kappa(scalar p, scalar T) const
+{
+ return kappa_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::kappag(scalar p, scalar T) const
+{
+ return kappag_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::sigma(scalar p, scalar T) const
+{
+ return sigma_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::D(scalar p, scalar T) const
+{
+ return D_.f(p, T);
+}
+
+
+inline Foam::scalar Foam::O2::D(scalar p, scalar T, scalar Wb) const
+{
+ return D_.f(p, T, Wb);
+}
+
+
+// ************************************************************************* //