shear stress slip BC

include/mm_names.h

@@ -4392,6 +4392,26 @@ struct BC_descriptions BC_Desc[] = {
      {1, 1, 1, 0, 0, 1, 1, 1, 0, 0},
      SINGLE_PHASE,
      DVI_SINGLE_PHASE_DB},
+    // {"SHEAR_STRESS_APPLIED",
+    //  "SHEAR_STRESS_APPLIED_BC",
+    //  STRONG_INT_SURF,
+    //  SHEAR_STRESS_APPLIED_BC,
+    //  R_MOMENTUM1,
+    //  SCALAR,
+    //  R_MOM_TANG1,
+    //  {1, 1, 1, 0, 0, 1, 1, 1, 0, 0},
+    //  SINGLE_PHASE,
+    //  DVI_SINGLE_PHASE_DB},
+    {"SHEAR_STRESS_APPLIED",
+     "SHEAR_STRESS_APPLIED_BC",
+     WEAK_INT_SURF,
+     SHEAR_STRESS_APPLIED_BC,
+     R_MOMENTUM1,
+     VECTOR,
+     NO_ROT,
+     {1, 1, 1, 0, 0, 1, 1, 1, 0, 0},
+     SINGLE_PHASE,
+     DVI_SINGLE_PHASE_DB},
     {"VELO_TANGENT_LS",
      "VELO_TANGENT_LS_BC",
      STRONG_INT_SURF,

include/mm_ns_bc.h

@@ -985,4 +985,12 @@ void dvzdr_zero_deriv_bc(double[DIM], /* func
 void ls_wall_angle_bc(double func[DIM],
                       double d_func[DIM][MAX_VARIABLE_TYPES + MAX_CONC][MDE],
                       const double angle);
+void shear_stress_applied(double func[DIM],
+                          double d_func[MAX_PDIM][MAX_VARIABLE_TYPES + MAX_CONC][MDE],
+                          const int id_side, /* ID of the side of the element             */
+                          double *user,      /* surface tension                           */
+                          struct elem_side_bc_struct *elem_side_bc,
+                          const int iconnect_ptr,
+                          dbl *xi, /* Natural coordinates of the integration point */
+                          const Exo_DB *exo);
 #endif /* GOMA_MM_NS_BC_H */

include/rf_bc_const.h

@@ -819,6 +819,7 @@
 #define CAPILLARY_TABLE_BC         963000007
 #define ELEC_TRACTION_BC           963000001 /* Include Maxwell Stress in CAPILLARY_BC */
 #define ELEC_TRACTION_SOLID_BC     963000003 /* Include Maxwell Stress in solid */
+#define SHEAR_STRESS_APPLIED_BC    963000009
 #define CAP_REPULSE_BC             963100000
 #define CAP_RECOIL_PRESS_BC        963110000
 #define CAP_REPULSE_ROLL_BC        963120000

src/bc_integ.c

@@ -890,6 +890,14 @@ int apply_integrated_bc(double x[],            /* Solution vector for the curren
                         iconnect_ptr, xi, exo);
           break;
 
+        case SHEAR_STRESS_APPLIED_BC:
+          memset(cfunc, 0, MDE * DIM * sizeof(double));
+          memset(d_cfunc, 0, MDE * DIM * (MAX_VARIABLE_TYPES + MAX_CONC) * MDE * sizeof(double));
+
+          shear_stress_applied(func, d_func, elem_side_bc->id_side, bc->BC_Data_Float,
+                               elem_side_bc, iconnect_ptr, xi, exo);
+          break;
+
         case CAPILLARY_BC:
         case CAP_REPULSE_BC:
         case CAP_REPULSE_ROLL_BC:

src/mm_fill_stress.c

@@ -7618,7 +7618,7 @@ int assemble_stress_sqrt_conf(dbl tt, /* parameter to vary time integration from
 
   SUPG_terms supg_terms;
   if (supg != 0.) {
-    supg_tau(&supg_terms, dim, 1e-14, pg_data, dt, TRUE, eqn);
+    supg_tau(&supg_terms, dim, 1e-8, pg_data, dt, TRUE, eqn);
   }
   /* end Petrov-Galerkin addition */
   dbl yzbeta_factor = 0.0;

src/mm_input_bc.c

@@ -2949,6 +2949,10 @@ void rd_bc_specs(FILE *ifp, char *input) {
 
       break;
 
+    case SHEAR_STRESS_APPLIED_BC:
+    // TODO Read Shear Stress Applied BC
+    break;
+
       /*
        *
        * Read in Lagrange multiplier boundary conditions and add