HydroGeoSphere/Modifying Zoned Properties

There are a number of instructions which can be used to modify the property values associated with a zone or group of zones. Before these instructions are issued, it is necessary to select the appropriate type of media and then choose the zones which you want to modify.

Modifying zone property values

For example, suppose that you wished to assign a new porous media hydraulic conductivity to all zones, and thus to all elements in the problem. The following set of instructions, inserted in the prefix.grok file would accomplish this:

    use domain type
    porous media

    clear chosen zones

    choose zones all

    k isotropic
    1.e-5

In this case, we are applying the instruction k isotropic to zones of type porous media, although it is equally valid to use it with dual-type zones. However, if you try to use this instruction with zones of type fracture or surface, warning messages will be issued to the screen and prefixo.eco file and execution will halt.

The instructions which are valid in specific situations are discussed in the relevant sections of the manual. For example, instructions which can be used for defining saturated flow properties are described in Section 5.8.2.

Another way to define zone properties is through the use of a material properties file, which should be located in the same directory as the prefix.grok file. This file contains lists of media-specific instructions which can be used to define properties for one or more named materials. These material properties can then be assigned to the current set of chosen zones. To assign new properties through the use of a material properties file, we first need to issue the following instruction:

Properties file

1. props_file_name Name of the material properties file.

This file will be searched for materials given as input to the Read properties instruction discussed below.

• • •

The Properties file instruction has two benefits: first, it allows you to create sets of material properties and give them meaningful file names and second, it allows you to easily switch between material property sets merely by changing the file name given in the prefix.grok file.

Any line in a material properties file which is completely blank or which begins with an exclamation point (!) is treated as a comment and is ignored by grok. This allows you to include comments whenever required.

Each distinct material in the file is identified by a unique label and may contain instructions which are to be applied to the current zone type. For example, instructions which can be used for defining porous media properties when simulating saturated flow (as described in Section 5.8.2.1) may be included. Figure 5.9 shows an example of a material defined for the verification problem discussed in Section 4.5.1:

    !------------------------------------------
    Porous medium

    k isotropic
    500.0

    specific storage
    0.0

    porosity
    1.0

    longitudinal dispersivity
    10.0

    transverse dispersivity
    0.1

    transverse vertical dispersivity
    0.1

    tortuosity
    0.1

    end material

Figure 5.9: Example of a Material Defined in an .mprops File

To make use of the material properties file, you would issue the following instruction:

Read properties

1. mat_name Name of the material.

Chosen zones are assigned properties of the material named mat_name in the current material properties file as defined in a properties file instruction.

• • •

So for example, the following set of instructions could be inserted in the prefix.grok file:

    use domain type
    porous media

    properties file
    my.mprops

    clear chosen zones
    
    choose zones all

    read properties
    sand

The instruction read Properties would, in this case, search the porous media material properties file my.mprops for a material named sand. If found, it would read the instructions defining the material and modify the porous media properties for the current set of chosen zones.

Example output

A summary of the final data which has been defined for each zone is listed in the prefixo.eco file, and an example is shown in Figure 5.10:

    ----------------------------------------------------------
    POROUS MEDIA PROPERTIES
    
    ZONE:  1
    MATERIAL: porous medium
    Consists of          100 elements out of           100
    Kxx:   500.000
    Kyy:   500.000
    Kzz:   500.000
    Specific storage:   0.00000
    Porosity:   1.00000

    Longitudinal dispersivity   10.0000
    Transverse dispersivity  0.100000
    Transverse vertical dispersivity  0.100000
    Tortuosity   0.100000
    Bulk density    2650.00
    Immobile zone porosity    0.00000
    Mass transfer coefficient    0.00000
       100 elements of      100  have been assigned properties

Figure 5.10: Example Output for a Porous Media Material

In this example, because flow is saturated, no variably-saturated porous media flow properties need to be defined in the material properties file. Also, default values for properties (e.g. bulk density, immobile zone porosity etc.) which have not been modified in the prefix.grok or material properties file are used.