waveSuperposition: Fixed bug in gas pressure

2018-12-06 18:18:36 +00:00
parent 2e1e8837c0
commit 4796ce32bd
1 changed files with 15 additions and 1 deletions
--- a/src/waves/waveSuperposition/waveSuperposition.C
+++ b/src/waves/waveSuperposition/waveSuperposition.C
@ -347,7 +347,21 @@ Foam::tmp<Foam::scalarField> Foam::waveSuperposition::pGas
    const vectorField& p
 ) const
 {
-    return - pLiquid(t, p);
+    tensor axes;
+    scalar u;
+    vectorField xyz(p.size());
+    transformation(p, axes, u, xyz);
+
+    axes = tensor(- axes.x(), - axes.y(), axes.z());
+
+    if (heightAboveWave_)
+    {
+        xyz.replace(2, height(t, p));
+    }
+
+    xyz.replace(2, - xyz.component(2));
+
+    return pressure(t, xyz);
 }