[cig-commits] commit: Do not apply damping force to the bottom

[Mercurial]

changeset:   452:cb0cc087a6c9
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Jun 20 13:24:16 2012 -0700
files:       Utils/src/BuoyancyDampingTerm.cxx
description:
Do not apply damping force to the bottom


diff -r 4853368a2fda -r cb0cc087a6c9 Utils/src/BuoyancyDampingTerm.cxx
--- a/Utils/src/BuoyancyDampingTerm.cxx	Wed Jun 13 15:08:55 2012 -0700
+++ b/Utils/src/BuoyancyDampingTerm.cxx	Wed Jun 20 13:24:16 2012 -0700
@@ -217,6 +217,18 @@ void _BuoyancyDampingTerm_AssembleElemen
         {
           const unsigned int face((local_norm==0 ? 1 :
                                    (local_norm==1 ? 0 : 2))*2+(sgn+1)/2);
+          IJK ijk;
+          RegularMeshUtils_Node_1DTo3D
+            (variable1->feMesh,
+             Mesh_DomainToGlobal(variable1->feMesh,MT_VERTEX,
+                                 elementType->faceNodes[face][0]),ijk);
+
+          /* Do not apply a damping force to the bottom boundary.
+             Assume that it is counteracted by whatever is below
+             it. */
+          if(local_norm==1 && sgn==-1 && ijk[1]==0)
+            continue;
+          
           for(int i=0; i<3; ++i)
             for(int j=0;j<(dim==2 ? 1 : 3); ++j)
               {