[cig-commits] commit: Fix upwind parameter for Q2 elements

[Mercurial]

changeset:   836:8cd0adf15d97
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Mon Jul 02 17:03:03 2012 -0700
files:       SLE/ProvidedSystems/AdvectionDiffusion/src/UpwindParameter.cxx
description:
Fix upwind parameter for Q2 elements


diff -r 823ca31cccd8 -r 8cd0adf15d97 SLE/ProvidedSystems/AdvectionDiffusion/src/UpwindParameter.cxx
--- a/SLE/ProvidedSystems/AdvectionDiffusion/src/UpwindParameter.cxx	Wed Jun 20 13:23:34 2012 -0700
+++ b/SLE/ProvidedSystems/AdvectionDiffusion/src/UpwindParameter.cxx	Mon Jul 02 17:03:03 2012 -0700
@@ -152,6 +152,10 @@ double AdvDiffResidualForceTerm_UpwindDi
       upwindDiffusivity+=xiUpwind*velocityCentre[dim_I]*lengthScale;
     }
   upwindDiffusivity*=ISQRT15;         /* See Eq. 3.3.11 */
+
+  upwindDiffusivity/=2; /* This factor is because we are using Q2
+                           instead of Q1 elements, so the spacing
+                           needs to be adjusted. */
 	
   return upwindDiffusivity;
 }