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檔案相同的目錄中。該檔案包含介質特定指令列表,這些指令可用於為一個或多個命名材料定義屬性。然後,這些材料屬性可以分配給當前的已選分割槽集。要透過使用材料屬性檔案分配新屬性,我們首先需要發出以下指令
- 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 檔案中定義的材料示例
要使用材料屬性檔案,您需要發出以下指令
- 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或材料屬性檔案中未修改的屬性的預設值(例如體積密度、固定區孔隙率等)。