From 3205337e81c181f2deb81b14c0ab5414746bba59 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://cfd.direct>
Date: Tue, 19 Jan 2016 21:20:03 +0000
Subject: [PATCH] scalarTransportFoam: Added support for steady-state solution
and all fvOptions Optional under-relaxation is provided for steady-state
solution. Added missing fvOptions.constrain and fvOptions.correct calls.
---
.../basic/scalarTransportFoam/scalarTransportFoam.C | 11 ++++++++---
1 file changed, 8 insertions(+), 3 deletions(-)
diff --git a/applications/solvers/basic/scalarTransportFoam/scalarTransportFoam.C b/applications/solvers/basic/scalarTransportFoam/scalarTransportFoam.C
index bef59237005..9082964083c 100644
--- a/applications/solvers/basic/scalarTransportFoam/scalarTransportFoam.C
+++ b/applications/solvers/basic/scalarTransportFoam/scalarTransportFoam.C
@@ -2,7 +2,7 @@
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration |
- \\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
+ \\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
@@ -25,7 +25,7 @@ Application
scalarTransportFoam
Description
- Solves a transport equation for a passive scalar
+ Solves the steady or transient transport equation for a passive scalar.
\*---------------------------------------------------------------------------*/
@@ -58,7 +58,7 @@ int main(int argc, char *argv[])
while (simple.correctNonOrthogonal())
{
- solve
+ fvScalarMatrix TEqn
(
fvm::ddt(T)
+ fvm::div(phi, T)
@@ -66,6 +66,11 @@ int main(int argc, char *argv[])
==
fvOptions(T)
);
+
+ TEqn.relax();
+ fvOptions.constrain(TEqn);
+ TEqn.solve();
+ fvOptions.correct(T);
}
runTime.write();
--
GitLab