Feature/oduran/test fem stress consistency

inputFiles/lagrangianContactMechanics/robustness_cases/cases/rc_hex_c1_lagr_contact.xml

@@ -0,0 +1,71 @@
+<?xml version="1.0" ?>
+<Problem>
+  <Mesh>
+    <VTKMesh name="mesh1" file="../meshes/fractured_mesh_hex_DFN_1.vtu" nodesetNames="{ f1_node_set }"/>
+  </Mesh>
+  <Geometry>
+    <Box name="xnegFace" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  0.01,  1.01,  1.01 }"/>
+    <Box name="xposFace" xMin="{  0.99, -0.01, -0.01 }" xMax="{  1.01,  1.01,  1.01 }"/>
+    <Box name="ynegFace" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  1.01,  0.01,  1.01 }"/>
+    <Box name="yposFace" xMin="{ -0.01,  0.99, -0.01 }" xMax="{  1.01,  1.01,  1.01 }"/>
+    <Box name="znegFace" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  1.01,  1.01,  0.01 }"/>
+    <Box name="zposFace" xMin="{ -0.01, -0.01,  0.99 }" xMax="{  1.01,  1.01,  1.01 }"/>
+    <Box name="fracture" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  1.01,  1.01,  1.01 }"/>
+  </Geometry>
+  <Solvers gravityVector="{0.0, 0.0, 0.0}">
+    <SolidMechanicsAugmentedLagrangianContact 
+      name="mechSolver"
+      simultaneous="1"
+      symmetric="1"
+      iterPenaltyN="1.0e1"
+      iterPenaltyT="1.0e-1"
+      tolNormalTrac="1.e-08"
+      tolTauLimit="1.e-08"
+      tolJumpN="1.e-8"
+      tolJumpT="1.e-8" 
+      discretization="FE1" 
+      targetRegions="{ Region, Fracture }" 
+      logLevel="1">
+      <NonlinearSolverParameters newtonTol="1.0e-5" newtonMaxIter="20" logLevel="1"/>
+      <LinearSolverParameters solverType="gmres" preconditionerType="amg" krylovTol="1.0e-10" logLevel="1"/>
+    </SolidMechanicsAugmentedLagrangianContact>
+    <SurfaceGenerator name="SurfaceGen" targetRegions="{ Region, Fracture }" fractureRegion="Fracture" initialRockToughness="10.0e9" logLevel="1"/>
+  </Solvers>
+  <NumericalMethods>
+    <FiniteElements>
+      <FiniteElementSpace name="FE1" order="1"/>
+    </FiniteElements>
+  </NumericalMethods>
+  <ElementRegions>
+    <CellElementRegion name="Region" cellBlocks="{ * }" materialList="{ rock }"/>
+    <SurfaceElementRegion name="Fracture" defaultAperture="1.0e-4" faceBlock="faceElementSubRegion" materialList="{ fractureContact }"/>
+  </ElementRegions>
+  <Constitutive>
+    <ElasticIsotropic name="rock" defaultDensity="2500" defaultYoungModulus="10e9" defaultPoissonRatio="0.2"/>
+    <Coulomb name="fractureContact" cohesion="0.0" frictionCoefficient="0.0"/>
+  </Constitutive>
+  <FieldSpecifications>
+    <FieldSpecification name="separableFace" fieldName="isFaceSeparable" initialCondition="1" setNames="{ f1_node_set }" objectPath="faceManager" scale="1" />
+    <FieldSpecification name="frac" initialCondition="1" setNames="{ f1_node_set }" objectPath="faceManager" fieldName="ruptureState" scale="1" />
+    <FieldSpecification name="xneg_disp" component="0" setNames="{ xnegFace }" objectPath="nodeManager" fieldName="totalDisplacement" scale="0.0"/>
+    <FieldSpecification name="yneg_disp" component="1" setNames="{ ynegFace }" objectPath="nodeManager" fieldName="totalDisplacement" scale="0.0"/>
+    <FieldSpecification name="zneg_disp" component="2" setNames="{ znegFace }" objectPath="nodeManager" fieldName="totalDisplacement" scale="0.0"/>
+    <Traction name="xpos_traction" setNames="{ xposFace }" objectPath="faceManager" tractionType="normal" scale="-10.0e6"/>
+    <Traction name="ypos_traction" setNames="{ yposFace }" objectPath="faceManager" tractionType="normal" scale="-10.0e6"/>
+    <Traction name="zpos_traction" setNames="{ zposFace }" objectPath="faceManager" tractionType="normal" scale="-10.0e6"/>
+  </FieldSpecifications>
+  <Tasks>
+    <DynamicFieldSpecification name="apply_fracture_updates" fieldSpecificationNames="{frac}"/>
+    <SolidMechanicsAugmentedLagrangianContactInitialization name="ELASTICITY.PRE.INIT.STEP" solidSolverName="mechSolver" logLevel="1"/>
+  </Tasks>  
+  <Outputs>
+    <VTK name="vtkOutput" />
+  </Outputs>
+  <Events minTime="-1.0e11" maxTime="1.0">
+    <SoloEvent name="preFracture" target="/Solvers/SurfaceGen" targetTime="-1.0e11" beginTime="-1.0e11"/>
+    <SoloEvent name="TriggerFractureUpdate" target="/Tasks/apply_fracture_updates" targetTime="-1.0e11" beginTime="-1.0e11"/>
+    <SoloEvent name="ELASTICITY.PRE.INIT.STEP" targetTime="-1.0e11" beginTime="-1.0e11" target="/Tasks/ELASTICITY.PRE.INIT.STEP"/>
+    <PeriodicEvent name="outputs" timeFrequency="1.0" target="/Outputs/vtkOutput"/>
+    <PeriodicEvent name="solverApplications" target="/Solvers/mechSolver" forceDt="1.0"/>
+  </Events>
+</Problem>

inputFiles/lagrangianContactMechanics/robustness_cases/cases/rc_hex_c1_lagr_contact_rs.xml

@@ -0,0 +1,80 @@
+<?xml version="1.0" ?>
+<Problem>
+  <Mesh>
+    <InternalMesh
+      name="mesh1"
+      elementTypes="{ C3D8 }"
+      xCoords="{ 0.0, 1.0 }"
+      yCoords="{ 0.0, 1.0 }"
+      zCoords="{ 0.0, 1.0 }"
+      nx="{ 4 }"
+      ny="{ 4 }"
+      nz="{ 4 }"
+      cellBlockNames="{ cb1 }"/>
+  </Mesh>
+  <Geometry>
+    <Box name="xnegFace" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  0.01,  1.01,  1.01 }"/>
+    <Box name="xposFace" xMin="{  0.99, -0.01, -0.01 }" xMax="{  1.01,  1.01,  1.01 }"/>
+    <Box name="ynegFace" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  1.01,  0.01,  1.01 }"/>
+    <Box name="yposFace" xMin="{ -0.01,  0.99, -0.01 }" xMax="{  1.01,  1.01,  1.01 }"/>
+    <Box name="znegFace" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  1.01,  1.01,  0.01 }"/>
+    <Box name="zposFace" xMin="{ -0.01, -0.01,  0.99 }" xMax="{  1.01,  1.01,  1.01 }"/>
+    <Box name="fracture" xMin="{ +0.49, -0.01, -0.01 }" xMax="{  0.51,  1.01,  1.01 }"/>
+  </Geometry>
+  <Solvers gravityVector="{0.0, 0.0, 0.0}">
+    <SolidMechanicsAugmentedLagrangianContact 
+      name="mechSolver"
+      simultaneous="1"
+      symmetric="1"
+      iterPenaltyN="1.0e1"
+      iterPenaltyT="1.0e-1"
+      tolNormalTrac="1.e-08"
+      tolTauLimit="1.e-08"
+      tolJumpN="1.e-8"
+      tolJumpT="1.e-8" 
+      discretization="FE1" 
+      targetRegions="{ Region, Fracture }" 
+      logLevel="1">
+      <NonlinearSolverParameters newtonTol="1.0e-5" newtonMaxIter="20" logLevel="1"/>
+      <LinearSolverParameters solverType="gmres" preconditionerType="amg" krylovTol="1.0e-10" logLevel="1"/>
+    </SolidMechanicsAugmentedLagrangianContact>
+    <SurfaceGenerator name="SurfaceGen" targetRegions="{ Region, Fracture }" fractureRegion="Fracture" initialRockToughness="10.0e9" logLevel="1"/>
+  </Solvers>
+  <NumericalMethods>
+    <FiniteElements>
+      <FiniteElementSpace name="FE1" order="1"/>
+    </FiniteElements>
+  </NumericalMethods>
+  <ElementRegions>
+    <CellElementRegion name="Region" cellBlocks="{ * }" materialList="{ rock }"/>
+    <SurfaceElementRegion name="Fracture" defaultAperture="1.0e-4" faceBlock="faceElementSubRegion" materialList="{ fractureContact }"/>
+  </ElementRegions>
+  <Constitutive>
+    <ElasticIsotropic name="rock" defaultDensity="2500" defaultYoungModulus="10e9" defaultPoissonRatio="0.2"/>
+    <Coulomb name="fractureContact" cohesion="0.0" frictionCoefficient="0.0"/>
+  </Constitutive>
+  <FieldSpecifications>
+    <FieldSpecification name="separableFace" fieldName="isFaceSeparable" initialCondition="1" setNames="{ fracture }" objectPath="faceManager" scale="1" />
+    <FieldSpecification name="frac" initialCondition="1" setNames="{ fracture }" objectPath="faceManager" fieldName="ruptureState" scale="1" />
+    <FieldSpecification name="xneg_disp" component="0" setNames="{ xnegFace }" objectPath="nodeManager" fieldName="totalDisplacement" scale="0.0"/>
+    <FieldSpecification name="yneg_disp" component="1" setNames="{ ynegFace }" objectPath="nodeManager" fieldName="totalDisplacement" scale="0.0"/>
+    <FieldSpecification name="zneg_disp" component="2" setNames="{ znegFace }" objectPath="nodeManager" fieldName="totalDisplacement" scale="0.0"/>
+    <Traction name="xpos_traction" setNames="{ xposFace }" objectPath="faceManager" tractionType="normal" scale="-10.0e6"/>
+    <Traction name="ypos_traction" setNames="{ yposFace }" objectPath="faceManager" tractionType="normal" scale="-10.0e6"/>
+    <Traction name="zpos_traction" setNames="{ zposFace }" objectPath="faceManager" tractionType="normal" scale="-10.0e6"/>
+  </FieldSpecifications>
+  <Tasks>
+    <SolidMechanicsAugmentedLagrangianContactInitialization name="ELASTICITY.PRE.INIT.STEP" solidSolverName="mechSolver" reset_displacements="0" logLevel="1"/>
+    <SolidMechanicsAugmentedLagrangianContactInitialization name="ELASTICITY.PRE.INIT.STEP.II" solidSolverName="mechSolver" reset_displacements="0" logLevel="1"/>    
+  </Tasks>  
+  <Outputs>
+    <VTK name="vtkOutput" />
+  </Outputs>
+  <Events minTime="0.0" maxTime="1.0">
+    <SoloEvent name="ELASTICITY.PRE.INIT.STEP" targetTime="0.0" beginTime="0.0" target="/Tasks/ELASTICITY.PRE.INIT.STEP"/>    
+    <SoloEvent name="preFracture" target="/Solvers/SurfaceGen" targetTime="0.0" beginTime="0.0" />
+    <SoloEvent name="ELASTICITY.PRE.INIT.STEP.II" targetTime="0.0" beginTime="0.0" target="/Tasks/ELASTICITY.PRE.INIT.STEP.II"/>
+    <PeriodicEvent name="outputs" timeFrequency="1.0" target="/Outputs/vtkOutput"/>
+    <PeriodicEvent name="solverApplications" target="/Solvers/mechSolver" forceDt="1.0"/>
+  </Events>
+</Problem>

inputFiles/lagrangianContactMechanics/robustness_cases/cases/rc_hex_c1_lagr_contact_testFEM.xml

@@ -0,0 +1,75 @@
+<?xml version="1.0" ?>
+<Problem>
+  <Mesh>
+    <VTKMesh name="mesh1" file="../meshes/fractured_wavy_mesh_hex_DFN_123.vtu" nodesetNames="{ f1_node_set, f2_node_set, f3_node_set }" partitionMethod="parmetis"/>
+  </Mesh>
+  <Geometry>
+    <Box name="xnegFace" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  0.01,  1.01,  1.01 }"/>
+    <Box name="xposFace" xMin="{  0.99, -0.01, -0.01 }" xMax="{  1.01,  1.01,  1.01 }"/>
+    <Box name="ynegFace" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  1.01,  0.01,  1.01 }"/>
+    <Box name="yposFace" xMin="{ -0.01,  0.99, -0.01 }" xMax="{  1.01,  1.01,  1.01 }"/>
+    <Box name="znegFace" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  1.01,  1.01,  0.01 }"/>
+    <Box name="zposFace" xMin="{ -0.01, -0.01,  0.99 }" xMax="{  1.01,  1.01,  1.01 }"/>
+    <Box name="fracture" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  1.01,  1.01,  1.01 }"/>
+  </Geometry>
+  <Solvers gravityVector="{0.0, 0.0, 0.0}">
+    <SolidMechanicsAugmentedLagrangianContact 
+      name="mechSolver"
+      simultaneous="1"
+      symmetric="1"
+      iterPenaltyN="1.0e1"
+      iterPenaltyT="1.0e-1"
+      tolNormalTrac="1.e-03"
+      tolTauLimit="1.e-03"
+      tolJumpN="1.e-3"
+      tolJumpT="1.e-3" 
+      discretization="FE1" 
+      targetRegions="{ Region, Fracture }" 
+      logLevel="1">
+      <NonlinearSolverParameters newtonTol="1.0e-5" newtonMaxIter="20" logLevel="1"/>
+      <LinearSolverParameters directParallel="0"/>
+    </SolidMechanicsAugmentedLagrangianContact>
+    <SurfaceGenerator name="SurfaceGen" targetRegions="{ Region, Fracture }" fractureRegion="Fracture" initialRockToughness="10.0e9"/>
+  </Solvers>
+  <NumericalMethods>
+    <FiniteElements>
+      <FiniteElementSpace name="FE1" order="1" useHighOrderQuadratureRule="1" />
+    </FiniteElements>
+  </NumericalMethods>
+  <ElementRegions>
+    <CellElementRegion name="Region" cellBlocks="{ * }" materialList="{ rock }"/>
+    <SurfaceElementRegion name="Fracture" defaultAperture="1.0e-4" faceBlock="faceElementSubRegion" materialList="{ fractureContact }"/>
+  </ElementRegions>
+  <Constitutive>
+    <ElasticIsotropic name="rock" defaultDensity="2500" defaultYoungModulus="10e9" defaultPoissonRatio="0.2"/>
+    <Coulomb name="fractureContact" cohesion="1.0e10" frictionCoefficient="0.5"/>
+  </Constitutive>
+  <FieldSpecifications>
+    <FieldSpecification name="separableFace" fieldName="isFaceSeparable" initialCondition="1" setNames="{ f1_node_set, f2_node_set, f3_node_set }" objectPath="faceManager" scale="1" />
+    <FieldSpecification name="frac" initialCondition="1" setNames="{ f1_node_set, f2_node_set, f3_node_set }" objectPath="faceManager" fieldName="ruptureState" scale="1" />
+
+    <FieldSpecification name="xneg_disp" component="0" setNames="{ xnegFace }" objectPath="nodeManager" fieldName="totalDisplacement" scale="0.0"/>
+    <FieldSpecification name="yneg_disp" component="1" setNames="{ ynegFace }" objectPath="nodeManager" fieldName="totalDisplacement" scale="0.0"/>
+    <FieldSpecification name="zneg_disp" component="2" setNames="{ znegFace }" objectPath="nodeManager" fieldName="totalDisplacement" scale="0.0"/>
+
+    <Traction name="xpos_traction" setNames="{ xposFace }" objectPath="faceManager" tractionType="normal" scale="-1.0e7"/>
+    <Traction name="ypos_traction" setNames="{ yposFace }" objectPath="faceManager" tractionType="normal" scale="-0.5e7"/>
+    <Traction name="zpos_traction" setNames="{ zposFace }" objectPath="faceManager" tractionType="normal" scale="-2.0e7"/>
+  </FieldSpecifications>
+  <Tasks>
+    <SolidMechanicsAugmentedLagrangianContactInitialization name="ELASTICITY.PRE.SPLIT.STEP" solidSolverName="mechSolver" logLevel="1" reset_displacements = "0"/>
+    <SolidMechanicsAugmentedLagrangianContactInitialization name="ELASTICITY.POST.SPLIT.STEP" solidSolverName="mechSolver" logLevel="1" reset_displacements = "0"/>
+  </Tasks>  
+  <Outputs>
+    <VTK name="vtkOutputM" outputRegionType="cell" />
+    <VTK name="vtkOutputF" outputRegionType="surface"/>
+  </Outputs>
+  <Events minTime="0.0" maxTime="1.0">
+    <SoloEvent name="ELASTICITY.PRE.SPLIT.STEP" targetTime="0.0" beginTime="0.0" target="/Tasks/ELASTICITY.PRE.SPLIT.STEP"/>
+    <SoloEvent name="preFracture" target="/Solvers/SurfaceGen" targetTime="0.0" beginTime="0.0"/> 
+    <SoloEvent name="ELASTICITY.POST.SPLIT.STEP" targetTime="0.0" beginTime="0.0" target="/Tasks/ELASTICITY.POST.SPLIT.STEP"/>
+    <PeriodicEvent name="solverApplications" target="/Solvers/mechSolver" forceDt="1.0"/>
+    <PeriodicEvent name="outputsM" timeFrequency="1.0" target="/Outputs/vtkOutputM"/>
+    <PeriodicEvent name="outputsF" timeFrequency="1.0" target="/Outputs/vtkOutputF"/>
+  </Events>
+</Problem>

inputFiles/lagrangianContactMechanics/robustness_cases/cases/rc_hex_c1_testFlow.xml

@@ -0,0 +1,101 @@
+<?xml version="1.0" ?>
+<Problem>
+  <Mesh>
+    <VTKMesh name="mesh1" file="../meshes/fractured_mesh_hex_DFN_123.vtu" nodesetNames="{ f1_node_set, f2_node_set, f3_node_set }" />
+  </Mesh>
+
+  <Geometry>
+    <Box name="xnegFace" xMin="{ -0.01, -0.01, -0.01 }" xMax="{  0.01,  1.01,  1.01 }"/>
+    <Box name="xposFace" xMin="{  0.99, -0.01, -0.01 }" xMax="{  1.01,  1.01,  1.01 }"/>
+  </Geometry>
+
+  <Solvers gravityVector="{0.0, 0.0, 0.0}">
+    <SinglePhaseFVM
+      name="flowSolver"
+      targetRegions="{ Region, Fracture }"
+      discretization="tpfa"
+      logLevel="1">
+      <NonlinearSolverParameters newtonTol="1.0e-6" newtonMaxIter="20"/>
+      <LinearSolverParameters directParallel="0"/>
+    </SinglePhaseFVM>
+    <SurfaceGenerator name="SurfaceGen" targetRegions="{ Region, Fracture }" fractureRegion="Fracture" initialRockToughness="10.0e9" logLevel="1"/>
+  </Solvers>
+
+  <NumericalMethods>
+    <FiniteVolume>
+      <TwoPointFluxApproximation name="tpfa"/>
+    </FiniteVolume>
+  </NumericalMethods>
+
+  <ElementRegions>
+    <CellElementRegion name="Region" cellBlocks="{ * }" materialList="{ rockMatrix, water }"/>
+    <SurfaceElementRegion name="Fracture" faceBlock="faceElementSubRegion" materialList="{ fractureMaterial, water }" defaultAperture="1.0e-7"/>
+  </ElementRegions>
+
+  <Constitutive>
+    <CompressibleSinglePhaseFluid
+      name="water"
+      defaultDensity="1.0"
+      defaultViscosity="1.0"
+      referenceDensity="1.0"
+      referenceViscosity="1.0"
+      referencePressure="0.0"
+      compressibility="0.0"
+      viscosibility="0.0"/>
+
+    <CompressibleSolidConstantPermeability
+      name="rockMatrix"
+      solidModelName="nullSolid"
+      porosityModelName="rockPorosity"
+      permeabilityModelName="rockPerm"/>
+
+    <CompressibleSolidConstantPermeability
+      name="fractureMaterial"
+      solidModelName="nullSolid"
+      porosityModelName="fracPorosity"
+      permeabilityModelName="fracPerm"/>
+
+    <ConstantPermeability name="rockPerm" permeabilityComponents="{ 1.0, 1.0, 1.0 }"/>
+    <ConstantPermeability name="fracPerm" permeabilityComponents="{ 1.0, 1.0, 1.0 }"/>
+
+    <NullModel name="nullSolid"/>
+
+    <PressurePorosity
+      name="rockPorosity"
+      defaultReferencePorosity="0.1"
+      referencePressure="0.0"
+      compressibility="0.0"/>
+
+    <PressurePorosity
+      name="fracPorosity"
+      defaultReferencePorosity="0.1"
+      referencePressure="0.0"
+      compressibility="0.0"/>
+  </Constitutive>
+
+  <FieldSpecifications>
+    <FieldSpecification name="separableFace" fieldName="isFaceSeparable" initialCondition="1" setNames="{ f1_node_set, f2_node_set, f3_node_set }" objectPath="faceManager" scale="1" />
+    <FieldSpecification name="frac" initialCondition="1" setNames="{ f1_node_set, f2_node_set, f3_node_set }" objectPath="faceManager" fieldName="ruptureState" scale="1" />
+    <FieldSpecification name="initialP" fieldName="pressure" initialCondition="1" setNames="{ all }" objectPath="ElementRegions" scale="1.5"/>
+    <FieldSpecification name="initialPf" fieldName="pressure" initialCondition="1" setNames="{ all }" objectPath="ElementRegions/Fracture/faceElementSubRegion" scale="2.0"/>
+    <FieldSpecification name="inletP" fieldName="pressure" setNames="{ xnegFace }" objectPath="faceManager" scale="2.0"/>
+    <FieldSpecification name="outletP" fieldName="pressure" setNames="{ xposFace }" objectPath="faceManager" scale="1.0"/>
+  </FieldSpecifications>
+
+  <Tasks>
+    <DynamicFieldSpecification name="apply_fracture_updates" fieldSpecificationNames="{initialPf}"/>
+  </Tasks>
+
+  <Outputs>
+    <VTK name="vtkOutputM" fieldNames="{ pressure }" outputRegionType="cell" onlyPlotSpecifiedFieldNames="1"/>
+    <VTK name="vtkOutputF" fieldNames="{ pressure }" outputRegionType="surface" onlyPlotSpecifiedFieldNames="1"/>
+  </Outputs>
+
+  <Events minTime="-1.0" maxTime="1.0">
+    <SoloEvent name="preFracture" target="/Solvers/SurfaceGen" targetTime="-1.0" beginTime="-1.0" />
+    <SoloEvent name="TriggerFractureUpdate" target="/Tasks/apply_fracture_updates" targetTime="-1.0" beginTime="-1.0"/>
+    <PeriodicEvent name="outputsM" timeFrequency="1.0" target="/Outputs/vtkOutputM"/>
+    <PeriodicEvent name="outputsF" timeFrequency="1.0" target="/Outputs/vtkOutputF"/>
+    <PeriodicEvent name="solverApplications" target="/Solvers/flowSolver" forceDt="1.0"/>
+  </Events>
+</Problem>