WEBVTT 00:00:00.000 --> 00:00:05.000 The letters "USGS" and the words "Science for a Changing World" appear in green then fade to black. 00:00:05.000 --> 00:00:11.666 This video describes how permeability is defined in SUTRA and specified in ModelMuse 00:00:11.666 --> 00:00:18.766 This sphere represents the permeability in three dimensions of the element in SUTRA. 00:00:18.766 --> 00:00:24.032 The permeability can have different values in different directions. 00:00:24.033 --> 00:00:28.066 In the direction of maximum permeability I'm going to set this at 3. 00:00:28.066 --> 00:00:33.632 In the direction of middle permeability, I will set it to 1.5. 00:00:33.633 --> 00:00:40.766 Right now the direction of middle permeability is back into the screen. 00:00:40.766 --> 00:00:46.699 The direction of minimum permeability I'll set at 0.5. 00:00:46.700 --> 00:00:49.166 That's up and down right now. 00:00:49.166 --> 00:00:56.999 we can also specify the directions of these three components of the permeability 00:00:57.000 --> 00:01:04.233 Angle 1 represents the direction of maximum permeability in the horizontal plane. 00:01:04.233 --> 00:01:09.133 I'm going to set that to 70. 00:01:09.133 --> 00:01:14.066 Now you see it is rotated around the horizontal plane back to an angle of 70°. 00:01:14.066 --> 00:01:21.766 Angle 2. represents the direction of the maximum permeability above the horizontal plane. 00:01:21.766 --> 00:01:28.532 Negative angles are used if the direction of the maximum permeability is downward below the horizontal plane. 00:01:28.533 --> 00:01:32.266 I'm going to set Angle 2 to 20. 00:01:32.266 --> 00:01:40.666 See that the ellipsoid has been rotated upward above the horizontal plaine in this direction. 00:01:40.666 --> 00:01:50.666 Finally, Angle 3 represents rotation of the permeability ellipsoid around the axis of maximum permeability. 00:01:50.666 --> 00:01:56.732 I'm going to set Angle 3 to 40. 00:01:56.733 --> 00:02:03.766 For 2D models, there is no middle permeability and only Angle 1 is needed. 00:02:03.766 --> 00:02:08.266 Let's animate this whole process of setting these angles. 00:02:08.266 --> 00:02:15.666 First, it rotates around this way around the horizontal plane. 00:02:15.666 --> 00:02:17.399 Then up 00:02:17.400 --> 00:02:21.566 Then finally around the long axis. 00:02:21.566 --> 00:02:26.166 That's how the permeability of this particular element is defined. 00:02:26.166 --> 00:02:31.032 The permeability can be sent separately for each element 00:02:31.033 --> 00:02:38.633 ModelMuse provides a way to use objects to help specify an Angle 1 of the permeability orientation. 00:02:38.633 --> 00:02:45.866 Here's how. First this capability is only needed if the permeability is anisotropic. 00:02:45.866 --> 00:03:01.666 For example, we'll set the value of the minimum hydraulic conductivity to be 1/3 the value of the maximum hydraulic conductivity. 00:03:01.666 --> 00:03:06.032 We'll do this for the middle hydraulic conductivity too. 00:03:06.033 --> 00:03:11.366 Divide by 3. 00:03:11.366 --> 00:03:15.999 Divide by 3. 00:03:16.000 --> 00:03:18.500 To see what the permeability looks like, 00:03:18.500 --> 00:03:27.100 we'll use Data Visualization dialog box to show the directions maximum and minimum hydraulic conductivity. 00:03:27.100 --> 00:03:32.766 So we go to "Data|Data Visualization". 00:03:32.766 --> 00:03:37.066 Select "Vectors" 00:03:37.066 --> 00:03:39.366 Hydraulic Conductivity 00:03:39.366 --> 00:03:43.866 Apply 00:03:43.866 --> 00:03:53.166 OK so these lines illustrate the directions of maximum, middle, and minimum hydraulic conductivity. 00:03:53.166 --> 00:03:57.332 They are black at the beginning and extend... 00:03:57.333 --> 00:04:03.266 Half way through they change to a different color. Red is for the maximum permeability 00:04:03.266 --> 00:04:07.932 Green is for the middle and blue is for the minimum hydraulic conductivity. 00:04:07.933 --> 00:04:15.199 Next we'll create a dataset for the maximum permeability direction and choose an interpolation method for it. 00:04:15.200 --> 00:04:18.166 We'll call it MaxPermeabilityDirection 00:04:18.166 --> 00:04:24.232 "Data|Edit Data Sets" "Add" 00:04:24.233 --> 00:04:38.966 Call it MaxPermeabilityDirection 00:04:38.966 --> 00:04:45.732 It is going to be evaluated at elements because the permeability is evaluated at elements. 00:04:45.733 --> 00:04:51.066 It will be a 2D dataset so we can choose an interpolation method. 00:04:51.066 --> 00:04:55.266 I'm going to choose Inverse Distance Squared for this. 00:04:55.266 --> 00:05:04.666 Then I will use this for the AngleHorizontal. 00:05:04.666 --> 00:05:08.066 This corresponds to Angle1 of the SUTRA input file. 00:05:08.066 --> 00:05:17.232 Edit the formula.. Data Sets|User Defined|MaxPermeabilityDirection 00:05:17.233 --> 00:05:18.833 Apply 00:05:18.833 --> 00:05:20.466 Close 00:05:20.466 --> 00:05:25.832 That didn't change anything because the default formula for MaxPermeabilityDirection 00:05:25.833 --> 00:05:31.233 was zero which is the same thing as the other formula was. 00:05:31.233 --> 00:05:41.999 One thing I forgot to do, though, was to set the units of MaxPermeabilityDirection to degrees 00:05:42.000 --> 00:05:50.500 and that's important. Interpolation won't be done properly unless the units are set to degrees or radians 00:05:50.500 --> 00:05:52.100 Click Apply 00:05:52.100 --> 00:05:53.766 Close 00:05:53.766 --> 00:05:59.132 And now we'll draw line objects. 00:05:59.133 --> 00:06:10.066 I'll uncheck "Set values of intersected elements." 00:06:10.066 --> 00:06:13.832 I'll set the number of Z formulas to zero. 00:06:13.833 --> 00:06:17.266 Set values by interpolation. 00:06:17.266 --> 00:06:21.332 Data Sets, User Defined, MaxPermeabilityDirection 00:06:21.333 --> 00:06:25.033 Here's a little bit tricky part. 00:06:25.033 --> 00:06:39.899 You have to set the formula to ObjectCurrentSegmentAngleDegrees 00:06:39.900 --> 00:06:41.800 Click OK 00:06:41.800 --> 00:06:45.100 OK. Look what's happened to the angles. 00:06:45.100 --> 00:06:50.666 They are now pointing in the same direction as the nearest segment on this object. 00:06:50.666 --> 00:06:56.599 Let's try another one. 00:06:56.600 --> 00:07:09.800 Do the same thing again. 00:07:09.800 --> 00:07:13.966 Now these are pointing in the direction we wanted them to go 00:07:13.966 --> 00:07:25.866 and there is a smooth variation in the permeability direction based on how far away from the object it is 00:07:25.866 --> 00:07:30.932 Now watch what happens if I set the units of MaxPermeabilityDirection incorrectly 00:07:30.933 --> 00:07:34.233 I'll go to "Data|Edit Data Sets..." 00:07:34.233 --> 00:07:42.933 I'll select the data set here and just wipe out the units. Click Apply 00:07:42.933 --> 00:07:45.033 Close 00:07:45.033 --> 00:07:48.399 See all these angles are missed up. 00:07:48.400 --> 00:07:51.566 Well, that's because the units aren't set correctly 00:07:51.566 --> 00:08:01.799 It somewhat different method for interpolating when you are using degrees or radians instead of using some other units., 00:08:01.800 --> 00:08:05.866 So I'm going to undo that change. Here we are back again. 00:08:05.866 --> 00:08:13.066 For more information, read the topic "Specifying Angles for SUTRA" in the ModelMuse help. 00:08:13.066 --> 00:08:18.066 The letters "USGS" and the words "Science for a Changing World" appear in green.