From 4b7f92935e1fcda5dd10a34aeeb6cd51a02a73e9 Mon Sep 17 00:00:00 2001
From: mattijs <mattijs>
Date: Thu, 13 Jan 2022 10:16:31 +0000
Subject: [PATCH] BUG: fixedJump: calculate jump in he. Fixes #2327

---
 .../energyJump/energyJump/energyJumpFvPatchScalarField.C   | 7 ++++++-
 .../energyJumpAMI/energyJumpAMIFvPatchScalarField.C        | 5 ++++-
 2 files changed, 10 insertions(+), 2 deletions(-)

diff --git a/src/thermophysicalModels/basic/derivedFvPatchFields/energyJump/energyJump/energyJumpFvPatchScalarField.C b/src/thermophysicalModels/basic/derivedFvPatchFields/energyJump/energyJump/energyJumpFvPatchScalarField.C
index 5fd72a0f7f..50952eb0e3 100644
--- a/src/thermophysicalModels/basic/derivedFvPatchFields/energyJump/energyJump/energyJumpFvPatchScalarField.C
+++ b/src/thermophysicalModels/basic/derivedFvPatchFields/energyJump/energyJump/energyJumpFvPatchScalarField.C
@@ -6,6 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2012-2017 OpenFOAM Foundation
+    Copyright (C) 2022 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -125,7 +126,11 @@ void Foam::energyJumpFvPatchScalarField::updateCoeffs()
 
         const labelUList& faceCells = this->patch().faceCells();
 
-        setJump(thermo.he(pp, Tbp.jump(), faceCells));
+        setJump
+        (
+            thermo.he(pp, Tbp+Tbp.jump(), faceCells)
+          - thermo.he(pp, Tbp, faceCells)
+        );
     }
 
     fixedJumpFvPatchField<scalar>::updateCoeffs();
diff --git a/src/thermophysicalModels/basic/derivedFvPatchFields/energyJump/energyJumpAMI/energyJumpAMIFvPatchScalarField.C b/src/thermophysicalModels/basic/derivedFvPatchFields/energyJump/energyJumpAMI/energyJumpAMIFvPatchScalarField.C
index ba45de7efd..0ff842683a 100644
--- a/src/thermophysicalModels/basic/derivedFvPatchFields/energyJump/energyJumpAMI/energyJumpAMIFvPatchScalarField.C
+++ b/src/thermophysicalModels/basic/derivedFvPatchFields/energyJump/energyJumpAMI/energyJumpAMIFvPatchScalarField.C
@@ -6,6 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2012-2017 OpenFOAM Foundation
+    Copyright (C) 2022 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -125,7 +126,9 @@ void Foam::energyJumpAMIFvPatchScalarField::updateCoeffs()
 
         const labelUList& faceCells = this->patch().faceCells();
 
-        jump_ = thermo.he(pp, Tbp.jump(), faceCells);
+        jump_ =
+            thermo.he(pp, Tbp+Tbp.jump(), faceCells)
+          - thermo.he(pp, Tbp, faceCells);
     }
 
     fixedJumpAMIFvPatchField<scalar>::updateCoeffs();