From dfc10f67a45fb4b82da42c6bb631744b63c113e1 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://openfoam.org>
Date: Thu, 14 Feb 2019 10:11:00 +0000
Subject: [PATCH] equationOfState::perfectFluid: Added Cp departure function

---
 .../specie/equationOfState/perfectFluid/perfectFluidI.H     | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/src/thermophysicalModels/specie/equationOfState/perfectFluid/perfectFluidI.H b/src/thermophysicalModels/specie/equationOfState/perfectFluid/perfectFluidI.H
index 770ed5a2ef..9204872f40 100644
--- a/src/thermophysicalModels/specie/equationOfState/perfectFluid/perfectFluidI.H
+++ b/src/thermophysicalModels/specie/equationOfState/perfectFluid/perfectFluidI.H
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     | Website:  https://openfoam.org
-    \\  /    A nd           | Copyright (C) 2012-2018 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2012-2019 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -101,7 +101,9 @@ inline Foam::scalar Foam::perfectFluid<Specie>::H(scalar p, scalar T) const
 template<class Specie>
 inline Foam::scalar Foam::perfectFluid<Specie>::Cp(scalar p, scalar T) const
 {
-    return 0;
+    return
+        sqr(p/(rho(p, T)*T))/this->R()
+      - sqr(Pstd/(rho(Pstd, T)*T))/this->R();
 }