[CIG-SHORT] problem definition for Landers-Hector Mine benchmark

[Carl W. Gable]

      > 
      > Hi all:
      > 
      > I would like to put the problem description together for the 
      > Landers-Hector Mine benchmark.
      > 
      > I have the following tidbits in my notes from the Golden workshop:
      > 
      > * Use mesh provided by Carl

What I'm planning to provide is:

1) Mesh with only the Landers Fault, low resolution
2) Mesh with only the Landers Fault, high resolution
3) Mesh with Landers and Hector Mine Faults, high resolution.

      > * Use Landers slip model provided by Yuri
      > * Use layered material model provided by Brad
      >   - viscosity is 10^18 Pa-s in lower-crust and mantle
      > * Boundaries (lateral sides and bottom) are fixed displacements
      > * Solution is computed at 0, 5, 10, 15 years
      > * no gravity
      > 
      > Issues that need to be resolved (please comment):
      > 
      > 1. How is the slip model going to be defined? Are people getting the 
      > slip model info from attributes embedded in the LaGriT output or is 
      > Yuri providing a separate data file?

Yuri provided a slip model. Jiangning Lu (MIT) converted the files
Yuri sent (strike,dip,rake) into AVS format files. I have then used the slip
defined on those surfaces that do not exactly match the
CFM fault surface, and interpolated the displacements onto
the CFM surfaces. I then use these surfaces to define
slip components (x,y,z displacement) for fault nodes
in the meshes I built.

So, slip model info will be provided as node attributes embedded in the LaGriT output.


      > 
      > 2. What depth range are we using for the lower-crust and mantle? I am 
      > guessing that we will use a viscosity of 10^25 Pa-s in the upper 
      > crust as we did for the reverse and strike-slip benchmarks, so at 
      > what depth do we jump to a viscosity of 10^18 Pa-s?
      > 
      > 3. What comprises the "solution"? I am guessing that we want 
      > displacements and velocities at all of the nodes as we have done for 
      > the other benchmarks. In this case, we have additional faults (i.e., 
      > Hector Mine faults), so we will want stress changes on those 
      > surfaces, right? I am guessing we want the traction vector computed 
      > at the quadrature points of the faces making up the Hector Mine 
      > faults.
      > 
      > Brad
      > _______________________________________________
      > CIG-SHORT mailing list
      > CIG-SHORT at geodynamics.org
      > http://geodynamics.org/cgi-bin/mailman/listinfo/cig-short