Mercurial > hg > cfcfd3
changeset 2249:cf4b7f61c12b
Eilmer3/4: Fixed bugs in manufacture sources terms for Euler case.
I introduced this bug back in rev 1631 when I changed to using sympy.
I had tested the Navier-Stokes case and was satisfied, but neglected
to run the Euler case. It turns out, I had not set viscosity
and thermal conductivity to zero when generating the source terms
for the Euler equations.
author | Rowan J. Gollan <r.gollan@uq.edu.au> |
---|---|
date | Fri, 05 Jun 2015 15:17:44 -0400 |
parents | 35dd7dc5d4d5 |
children | 314b46321471 |
files | examples/eilmer3/2D/mms/make_source_terms.py examples/eilmer4/2D/manufactured-solution/make_source_terms.py |
diffstat | 2 files changed, 6 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- a/examples/eilmer3/2D/mms/make_source_terms.py Fri Jun 05 15:05:29 2015 -0400 +++ b/examples/eilmer3/2D/mms/make_source_terms.py Fri Jun 05 15:17:44 2015 -0400 @@ -24,6 +24,9 @@ mu, k = symbols('mu k') mu = 10.0 k = Cp*mu/Prandtl +if case == 1 or case == 3: + mu = 0.0 + k = 0.0 # Thermodynamic behvaiour, equation of state and energy equation e, T, et = symbols('e T et')
--- a/examples/eilmer4/2D/manufactured-solution/make_source_terms.py Fri Jun 05 15:05:29 2015 -0400 +++ b/examples/eilmer4/2D/manufactured-solution/make_source_terms.py Fri Jun 05 15:17:44 2015 -0400 @@ -24,6 +24,9 @@ mu, k = symbols('mu k') mu = 10.0 k = Cp*mu/Prandtl +if case == 1 or case == 3: + mu = 0.0 + k = 0.0 # Thermodynamic behvaiour, equation of state and energy equation e, T, et = symbols('e T et')