Backporting ConstantLewisDiffusivity bugfix

libantioch · Jun 23, 2015 · f425a1f · f425a1f
1 parent 186e5a7
commit f425a1f
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git 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<typename StateType>
     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
@@ -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") );
 }