From 4f7019edc30b189ae15e4eae71985c802371416c Mon Sep 17 00:00:00 2001
From: sergio <s.ferraris@opencfd.co.uk>
Date: Fri, 22 Apr 2022 10:01:32 -0700
Subject: [PATCH] ENH: Create output for Su VollerPrakash

---
 .../porousModels/VollerPrakash/VollerPrakash.C       | 12 +++++++++++-
 1 file changed, 11 insertions(+), 1 deletion(-)

diff --git a/src/phaseSystemModels/multiphaseInter/phasesSystem/interfaceModels/porousModels/VollerPrakash/VollerPrakash.C b/src/phaseSystemModels/multiphaseInter/phasesSystem/interfaceModels/porousModels/VollerPrakash/VollerPrakash.C
index 721076e61d..2771f3bc0c 100644
--- a/src/phaseSystemModels/multiphaseInter/phasesSystem/interfaceModels/porousModels/VollerPrakash/VollerPrakash.C
+++ b/src/phaseSystemModels/multiphaseInter/phasesSystem/interfaceModels/porousModels/VollerPrakash/VollerPrakash.C
@@ -5,7 +5,7 @@
     \\  /    A nd           | www.openfoam.com
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
-    Copyright (C) 2017 OpenCFD Ltd.
+    Copyright (C) 2017-2020 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -68,6 +68,16 @@ Foam::porousModels::VollerPrakash::S() const
     const volScalarField& solidAlpha =
         mesh_.lookupObject<volScalarField>(solidPhase_);
 
+    if (mesh_.time().outputTime() && debug)
+    {
+        volScalarField Su
+        (
+            "Su" + solidAlpha.name(),
+            Cu_*sqr(solidAlpha)/(pow3(1.0 - solidAlpha) + 1e-3)
+        );
+        Su.write();
+    }
+
     return Cu_*sqr(solidAlpha)/(pow3(1.0 - solidAlpha) + 1e-3);
 }