跳转到内容

HydroGeoSphere/修改分区属性

来自维基教科书,开放的书籍,开放的世界

有许多指令可以用来修改与分区或分区组关联的属性值。在发出这些指令之前,需要选择合适的介质类型,然后选择要修改的分区。

修改分区属性值

[编辑 | 编辑源代码]

例如,假设您希望为所有分区分配新的多孔介质水力传导率,从而为问题中的所有单元分配新的水力传导率。以下指令集插入到prefix.grok文件中可以完成此操作

    use domain type
    porous media

    clear chosen zones

    choose zones all

    k isotropic
    1.e-5


在这种情况下,我们正在应用指令k 各向同性到多孔介质类型分区,虽然它同样适用于双类型分区。但是,如果您尝试将此指令与断裂或表面类型的分区一起使用,将会向屏幕发出警告消息,并且prefixo.eco文件和执行将停止。

手册的相关部分讨论了在特定情况下有效的指令。例如,第 5.8.2 节描述了可用于定义饱和流动特性的指令。

定义分区属性的另一种方法是使用材料属性文件,该文件应位于与prefix.grok文件相同的目录中。该文件包含介质特定指令列表,这些指令可用于为一个或多个命名材料定义属性。然后,这些材料属性可以分配给当前的已选分区集。要通过使用材料属性文件分配新属性,我们首先需要发出以下指令

属性文件

[编辑 | 编辑源代码]
  1. props_file_name 材料属性文件名称。

将搜索该文件以查找作为输入给出的材料读取属性下面讨论的指令。

• • •

属性文件指令有两个好处:首先,它允许您创建材料属性集并为它们提供有意义的文件名;其次,它允许您只需更改prefix.grok文件中的文件名即可轻松地在材料属性集之间切换。

材料属性文件中任何完全为空白或以感叹号 (!) 开头的行都被视为注释,并将被grok忽略。这允许您在需要时包含注释。

文件中每个不同的材料都由一个唯一的标签标识,并且可能包含要应用于当前分区类型的指令。例如,可以包括在模拟饱和流动时用于定义多孔介质特性的指令(如第 5.8.2.1 节所述)。图 5.9显示了为第 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
图 5.9:.mprops 文件中定义的材料示例

要使用材料属性文件,您需要发出以下指令

读取属性

[编辑 | 编辑源代码]
  1. mat_name 材料名称。

已选分区被分配了在当前材料属性文件中名为mat_name的材料的属性,如属性文件指令中定义的那样。

• • •

因此,例如,以下指令集可以插入到prefix.grok文件中

    use domain type
    porous media

    properties file
    my.mprops

    clear chosen zones
    
    choose zones all

    read properties
    sand


指令读取属性在这种情况下,将搜索多孔介质材料属性文件my.mprops以查找名为sand的材料。如果找到,它将读取定义材料的指令并修改当前已选分区集的多孔介质属性。

示例输出

[编辑 | 编辑源代码]

对每个分区定义的最终数据的摘要列在prefixo.eco文件中,示例见图 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
图 5.10:多孔介质材料的示例输出

在此示例中,由于流动是饱和的,因此无需在材料属性文件中定义可变饱和多孔介质流动属性。此外,使用prefix.grok或材料属性文件中未修改的属性的默认值(例如体积密度、固定区孔隙率等)。

华夏公益教科书