該算例來自comsol自帶算例，
C:\Program Files\COMSOL\Multiphysics\applications\Subsurface_Flow_Module\Fluid_Flow\ discrete_fracture.mph
模型建立流程
1. 設置求解環境
（1）啓動Comsol Multiphysics，在“新建”任務窗格中單擊“模型導向”。
（2）在“選擇空間維度”任務窗格中單擊“三維”選項，在打開的“選擇物理場”任務窗格中單擊“流體流動”|“多孔介質和地下水流”|“Darcy定律（dl）”模塊，單擊“增加”按鈕，然後單擊“研究”按鈕。
（3）在“選擇研究”任務窗格中單擊 “預置研究”|“瞬態”選項，然後單擊“完成”按鈕，進入“模型開發器”界面。
2.模型創建
（1）在“模型開發器”任務窗格的“組件1”選項中，右擊“定義”，在彈出的快捷菜單選擇“變量”。選擇”變量”，在“設定”任務窗格中設定“長度單位”爲“m”。
（2）在“模型開發器”任務窗格的“組件1”選項中，選擇“幾何1”，在彈出的快捷菜單中選擇“工作平面”命令。
（3）選擇“工作平面1”|“平面幾何”，右擊“平面幾何”，在彈出的快捷菜單中選擇“貝氏多線段”命令。選擇“貝氏多線段”，在右側“設定”任務窗格中設定貝氏多線段參數，點擊“全部構建”。同理增加貝氏多線段2，設定參數。
（4）選擇“工作平面1”|“平面幾何”，右擊“平面幾何”在彈出的快捷菜單中選擇“長方體”命令。
（5）在“模型開發器”任務窗格的“組件1”選項中，選擇“幾何1” |“工作平面1”，右擊“平面幾何”，在彈出的快捷菜單中選擇“布爾運算與分割”|“差集”。然後在”設定”任務窗格中設定“差集”輸入對象。
（6）在“模型開發器”任務窗格的“組件1”選項中，選擇“幾何1”|“形成聯合體”，在“設定”任務窗格中點擊“全部構建”。
在“模型開發器”任務窗格的“組件1”選項中，選擇“幾何1”，右擊“工作平面”，在彈出的快捷菜單中選擇“拉伸”。然後在“設定”任務窗格中設定“拉伸”參數，點擊“構建所有對象”。
4. Darcy定律設置
（1）在“模型開發器”任務窗格的“組件1”|“Darcy定律”選項中，選擇“初始值1”，在“設定”任務窗格中設定參數。
（2）在“主菜單”中選擇“物理場”|“域”|“流體和基本屬性”，單擊“流體和基本屬性”增加到組件1。同理選擇選擇“物理場”|“域”|“儲水模型”增加到組件1。
（3）在“模型開發器”任務窗格的“組件1”選項中，選擇“Darcy定律（dl）”|“儲水模型1”，在右側“設定”任務窗格中設置參數。
（4）在“主菜單”中選擇“物理場”|“邊界”|“裂隙流”，單擊“裂隙流”增加到組件1。
在“模型開發器”任務窗格的“組件1”選項中，選擇“Darcy定律（dl）”|“裂隙流”，在右側“設定”任務窗格中輸入對象（中間裂隙）。
（5）在“模型開發器”任務窗格的“組件1”選項中，選擇“Darcy定律（dl）”，右擊“裂隙流1”，在彈出的快捷菜單中選擇“儲水模型”。選擇“儲水模型”，在右側“設定”任務窗格中設定參數。
（6）在“模型開發器”任務窗格的“組件1”選項中，選擇“Darcy定律（dl）”，右擊“裂隙流1”，在彈出的快捷菜單中選擇“壓力”。選擇“壓力1”，在右側“設定”任務窗格中設定參數。同理增加壓力2，選擇“壓力2”，在右側“設定”任務窗格中設定參數。
5. 網格參數設置
在“模型開發器”任務窗格的“組件1”選項中，右擊”網格1”，在彈出的快捷菜單中選擇”自由剖分四面體網格”。選擇“網格1”|“尺寸”， 在“設定”任務窗格中設定“單元尺寸”爲“粗化”。選擇“網格1”|“自由剖分四面體網格1”，在“設定”任務窗格中選擇“全部構建”。
7. 求解設置
在“模型開發器”任務窗格的“研究1”選項中，選擇“研究1”，在”設定”任務窗格中，單擊“計算”。
8. 結果分析
（1）在“模型開發器”任務窗格的“組件1”|“結果”選項中，選擇“壓力（dl），在“設定”任務窗格中，選擇時間1000s，選擇“切片1”，在“設定”任務窗格中設定參數，點擊“繪製”，即可得到壓力切片圖。
（2）在“模型開發器”任務窗格的“組件1”|“結果”選項中，選擇“壓力（dl）1，在“設定”任務窗格中，選擇時間1000s，在“設定”任務窗格中設定參數，點擊“繪製”，即可得到壓力錶面圖。
（3）在“模型開發器”任務窗格的“組件1”選項中，右擊“結果”，在彈出的快捷菜單中選擇“三維繪圖組”命令，選擇“三維繪圖組”在“設定”任務窗格中設定參數。右擊“裂隙壓力等值面圖”， 在彈出的快捷菜單中選擇“等值面”命令，選擇“等值面”，在“設定”任務窗格中設定參數。右擊“裂隙壓力等值面圖”， 在彈出的快捷菜單中選擇“面箭頭”命令。點擊“繪製”，得到裂隙壓力等值面圖。
作者：ZhaoH