From 81f08283f4d9bf15e5347ba0889d8f5bbeba62b2 Mon Sep 17 00:00:00 2001
From: henry <Henry Weller h.weller@opencfd.co.uk>
Date: Tue, 7 Oct 2008 10:07:43 +0100
Subject: [PATCH] Added momentum equation relaxation to allow forcing
diagonal-dominance.
---
applications/solvers/incompressible/oodles/oodles.C | 3 +++
1 file changed, 3 insertions(+)
diff --git a/applications/solvers/incompressible/oodles/oodles.C b/applications/solvers/incompressible/oodles/oodles.C
index eb5fa65dcf6..8c51065db1d 100644
--- a/applications/solvers/incompressible/oodles/oodles.C
+++ b/applications/solvers/incompressible/oodles/oodles.C
@@ -66,6 +66,9 @@ int main(int argc, char *argv[])
+ sgsModel->divDevBeff(U)
);
+ // Optionally ensure diagonal-dominance of the momentum matrix
+ UEqn.relax();
+
if (momentumPredictor)
{
solve(UEqn == -fvc::grad(p));
--
GitLab