From f425a1f1a884fc66e1e87cef100b725499ae8082 Mon Sep 17 00:00:00 2001 From: "Paul T. Bauman" Date: Tue, 23 Jun 2015 15:47:46 -0500 Subject: [PATCH] Backporting ConstantLewisDiffusivity bugfix --- src/transport/include/antioch/constant_lewis_diffusivity.h | 2 +- test/constant_lewis_unit.C | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/src/transport/include/antioch/constant_lewis_diffusivity.h b/src/transport/include/antioch/constant_lewis_diffusivity.h index 2b0241ae..577ea756 100644 --- a/src/transport/include/antioch/constant_lewis_diffusivity.h +++ b/src/transport/include/antioch/constant_lewis_diffusivity.h @@ -52,7 +52,7 @@ namespace Antioch template ANTIOCH_AUTO(StateType) D( const StateType& rho, const StateType& cp, const StateType& k ) const - ANTIOCH_AUTOFUNC(StateType, _Le*k/(rho*cp)) + ANTIOCH_AUTOFUNC(StateType, k/(_Le*rho*cp)) }; diff --git a/test/constant_lewis_unit.C b/test/constant_lewis_unit.C index 4d950b0c..75d7a399 100644 --- a/test/constant_lewis_unit.C +++ b/test/constant_lewis_unit.C @@ -75,7 +75,7 @@ int tester() const Scalar D = diff.D( rho, cp, k ); - const Scalar D_exact = Le*k/(rho*cp); + const Scalar D_exact = k/(Le*rho*cp); return test_val( D, D_exact, tol, std::string("D") ); }