From 22a6bb9ef2ccb5b0b7f31c213fc524823e101fe3 Mon Sep 17 00:00:00 2001 From: Ray Speth Date: Tue, 2 Sep 2014 22:34:16 +0000 Subject: [PATCH] [1D/Python] Fix initial profile generation for diffusion flames Don't assume that the leftmost grid point is z = 0. --- interfaces/cython/cantera/onedim.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/interfaces/cython/cantera/onedim.py b/interfaces/cython/cantera/onedim.py index 698444815..d99f6e685 100644 --- a/interfaces/cython/cantera/onedim.py +++ b/interfaces/cython/cantera/onedim.py @@ -415,7 +415,7 @@ class CounterflowDiffusionFlame(FlameBase): Y = np.zeros((nz, self.gas.n_species)) T = np.zeros(nz) for j in range(nz): - x = zz[j] + x = zz[j] - zz[0] zeta = f * (x - x0) zmix = 0.5 * (1.0 - erf(zeta)) if zmix > zst: @@ -427,7 +427,7 @@ class CounterflowDiffusionFlame(FlameBase): T[0] = T0f T[-1] = T0o - zrel = zz/dz + zrel = (zz - zz[0])/dz self.set_profile('u', [0.0, 1.0], [u0f, -u0o]) self.set_profile('V', [0.0, x0/dz, 1.0], [0.0, a, 0.0])