mikezero學(xué)習(xí)教程3flow model沙模塊

關(guān)于本文檔（Aboutthis 編寫本文檔的目的 用戶所需的知識背景(AssumedUser MIKEZero模型編輯器程序界面結(jié)構(gòu)(GeneralEditor 系統(tǒng)(Online 引言 概述(General 說明(Getting 示例 概述 環(huán)流示例(Sandtransportincirculation 沙模塊（SandTransport 模型定義（Model 時間參數(shù)（Time 求解格式（Solution 泥沙特性（Sediment 床面阻力 驅(qū)動力參數(shù) 擴散 對流擴散（Vertical 推薦設(shè)置 mended 源匯項 源匯項設(shè)定（Source 備注與提示（Remarksand 初始條件（Initial 泥沙粒徑組分濃度（Fraction 床面活動層厚度（Layer 邊界條件（Boundary 地貌演化模塊（Morphology 輸出 圖形化查看（Geographic 輸出類型（Output 輸出項（Output 地貌演化模塊（Morphology 概述（General 模型定義（Model 時間參數(shù)（Time 河岸侵蝕（Bank 邊界說明（Boundary 關(guān)于本文檔（Aboutthis編寫本文檔的目的用戶所需的知識背景(AssumedUserViewer）以及網(wǎng)格生成工具（MeshGenerator）。上述Mike3基本組件的說明文檔詳見C:\ProgramFiles\DHI\MIKEZero\Manuals\MIKE_ZERO\MzGeneric.pdf。MIKEZero模型編輯器程序界面結(jié)構(gòu)(GeneralEditor系統(tǒng)(Online圖1.1MIKE3FlowModelFM系引言Mike3FlowModelFM沙模塊主要模擬單一波浪、水流或波流共同作用下沙的侵蝕、說明(GettingMIKE3FlowModelFM戶手冊和包括水動力模塊的內(nèi)容翔實的“跟我學(xué)”教學(xué)指南。PDF格式的用戶手冊和教學(xué)DHI軟件的安裝程序中包括兩份內(nèi)容翔實的用戶教學(xué)指南文檔，一份為純水流動力這兩份文檔的PDF格式文件一般于：C:\ProgramFiles\DHI\MIKEZero\Manuals\MIKE_3\示例MIKE3FlowModelFM模型的最好方法之一。因此模型安裝程序中提查看該條目的作用及其變化對模擬結(jié)果的影響。以下為程序提供的所有示例文件的：U型河道示例C:\Program環(huán)流示例C:\Program沙模塊（SandTransport問題是非常必要的（限泥沙補給和高床面抗沖性條件）。

EngelundandHansen(全沙VanRijn(推移質(zhì)+懸移質(zhì)輸沙EngelundandFredsoe(推移質(zhì)+懸移質(zhì)輸沙Meyer-PeterandMuller(推移質(zhì)輸沙當用戶設(shè)定模型中存在推移質(zhì)輸沙時，相應(yīng)輸沙中則包含一推移質(zhì)系 ，該 當用戶設(shè)定模型中存在懸移質(zhì)輸沙時，相應(yīng)輸沙中則包含一懸移質(zhì)系 ，該計算結(jié)果超出的適用范圍。全沙（Total因EngelundandHansen僅計算全沙輸沙率，無法直接得出模型計算所需的推移懸沙濃度上限（umConcentration） 泥沙輸移表（Sedimenttransporttable）“泥沙輸移表”文件是由用戶事先基于Mike21Toolbox中的“GenerationofQ3D給定的“泥沙輸移表”（Sedimenttransporttable）線性插值的結(jié)果，故當用戶所生成的于某些參數(shù)出現(xiàn)異常的情況，程序的計算過程將中斷。此外，對于前200次變量的數(shù)值超過用戶給定最大值的情況，程序?qū)⒂涗浵鄳?yīng)警告信息至log文件中。VS--d--g--泥沙輸移計算的時間步長為時間設(shè)置框中設(shè)定的主時間步長（overalltimestep）乘以時間倍乘系數(shù)（timestepfactor）。及臨界CFL數(shù)的設(shè)置。1?2CFL數(shù)的數(shù)值為近似預(yù)估值，故在這種默認設(shè)置情況下仍然CFL值適當減小必須，當用戶將最小和最大時間步長均設(shè)定為與主時間步長相同時(p23)，模型將小于1。當選擇采用波流共同作用下泥沙輸移模式計算時，用戶還需設(shè)定床沙幾何標準方差g常數(shù)（空間均一d50（mm單位）和幾何標準方差σg2d網(wǎng)格化數(shù)據(jù)文件（dfs2）2d無dfs2文件時，相應(yīng)計算網(wǎng)格上的泥沙特性參 mended言，模型適用的床沙范圍為0.06mm至2mm。通常相對泥沙密度取值為2.65。床面阻力（Bed無床面阻力（光滑壁面系數(shù)（Chezy動床阻力 水動力模型中生成的阻力系數(shù) fromHD 水深變化而變化。系數(shù)（C）或系數(shù)可由下式給出：thetopographyofapointbarcrestto emorerounded(asreportedbyTalmon,1992).系數(shù)C（m1/2/s、系數(shù)（m1/3/s）和動床阻力可基于如下兩種方式設(shè)定常數(shù)（空間均一2d網(wǎng)格化數(shù)據(jù)文件（dfs2）2d無結(jié)構(gòu)數(shù)據(jù)文件（dfsudfsu文件時，當采用dfsu文件時，相應(yīng)計算網(wǎng)格上的泥沙特性參數(shù)基于分段常數(shù)插值方法生成；當采用dfs2文件時，相應(yīng)計算網(wǎng)格上的泥沙特性參數(shù)基于雙線性插值方法生成。 mended無波浪條件（Nowaves）用戶設(shè)定波浪條件（Wavefield）常數(shù)時間上變化，空間上不變（Varyingintime,constantin時間上不變，空間上變化(Constantintime,varyingin時間上、空間上均變化(Varyingintimeand用戶所提供的參數(shù)必須為均或有效。時間上變化，空間上不變（Varyingintime,constantin（dfs0時間步長不一致，則程序?qū)⒆詣舆M行相應(yīng)插值處理（見圖51。時間上不變，空間上變化（Constantintime,varyingin這種情況下，用戶需準備一空間上至少完整包括模擬區(qū)域范圍且包含相關(guān)參數(shù)的2d網(wǎng)（dfsu 這種情況下，用戶需準備一空間上至少完整包括模擬區(qū)域范圍且包含相關(guān)參數(shù)的2d網(wǎng)（dfsu計算的時間步長不一致時，則程序?qū)⒆詣舆M行相應(yīng)插值處理（見圖51。圖5.1（dfsueffectsofnon-resolvedprocessescanbesignificant,inwhichcasethedispersioncoefficientformallyshoulddependonthe3模型中的擴散項主要用來描述因不確定的物理過程引起的輸移問題。區(qū)分平流擴散無擴散(No擴散系數(shù)(Dispersioncoefficient渦粘系數(shù)類比(Scalededdyviscosity常數(shù)（時間、空間上定常）(Constant(inbothtime 空間上變化(Varying mended 點然后單擊“goto…”按鈕進行相關(guān)設(shè)置。源匯項設(shè)定（Sourcespecification） 懸沙濃度是獨立于水體的。源匯項相關(guān)參數(shù)數(shù)值的給定可以采用兩種方式常數(shù)（時間上不變常數(shù)（空間上均一搭建模型時，模型配置編輯器(set-upeditor)首先會檢查網(wǎng)格文件，自動找出網(wǎng)格文件中的邊界，并給每個邊界預(yù)先指定一個名稱，用戶可自行在框中對各個邊界的邊界設(shè)置（BoundarySpecification）在陸地邊界平衡輸沙條件（Euilibrium用戶設(shè)定取值(Specified常數(shù)（時間上、空間上不變線性插值分段立方插值(picewise正常順序逆序(reverse(Softstart在“輸出”框中，用戶可以單擊”NewOutput”命令添加新的輸出條目，或是”Deleteoutput”移除某一輸出條目。對于每個輸出條目，不管該條目是否在本次計算中外用戶亦可以單擊”View”命令基于MIKEZEROVIEWING/EDITING工具查看相關(guān)內(nèi)容。輸出格式（Ouputformat）相應(yīng)輸出結(jié)果文件為dfs0格式。出結(jié)果文件為dfs0格式。5.2ProfileDataDataDataData是否輸出干濕點（Treatmentoffloodand擇使用干濕判斷功能，則相應(yīng)干水深、濕水深和淹沒水深取值均為0。點序列（Pointseries）對于欲輸出數(shù)據(jù)空間點的坐標信息，用戶既可手工從設(shè)置框內(nèi)輸入，亦可從某一文件導(dǎo)入。坐標信息文件為一ascii文本文件，格式為“XYZName”，其中坐標信用戶必須設(shè)定欲輸出數(shù)據(jù)空間點坐標信息對應(yīng)的正確的坐標系信息（如Long/lat，UTM-32等?？臻g點所在網(wǎng)格單元中心點處的原始數(shù)據(jù)，每一網(wǎng)格單元號及其中心點處坐標記錄于log根據(jù)其所處網(wǎng)格單元頂點處數(shù)值進行線性插值的結(jié)果，網(wǎng)格單元各頂點處數(shù)值基于pseudo-Lacian方法（HolmesandConnell，1989）求得。每一網(wǎng)格單元號及其中心點處坐標記錄于log文件中。對于欲輸出線序列的坐標信息，用戶既可手工從設(shè)置框內(nèi)輸入兩端點，亦可從某一文件導(dǎo)入。坐標信息文件為一ascii文本文件，格式為“XYZ Name，信息XYZ列必須采用浮點數(shù)字格式表示且Z坐標僅用于3d輸出項的數(shù)據(jù)輸出，Name用戶必須設(shè)定欲輸出數(shù)據(jù)空間點坐標信息對應(yīng)的正確的坐標系信息（Long/lat，UTM-32等。空間點所在網(wǎng)格單元中心點處的原始數(shù)據(jù)，每一網(wǎng)格單元號及其中心點處坐標記錄于log根據(jù)其所處網(wǎng)格單元頂點處數(shù)值進行線性插值的結(jié)果，網(wǎng)格單元各頂點處數(shù)值基于pseudo-Lacian方法（HolmesandConnell，1989）求得。每一網(wǎng)格單元號及其中心點處坐標記錄于log文件中。用戶可選擇輸出某一封閉多邊形區(qū)域內(nèi)的網(wǎng)格單元上的計算結(jié)果，用戶需設(shè)定該多邊形的各頂點坐標程序默認第一個和最后一個坐標點為封閉所定義多邊形的一條線段的兩端點對于所定義多邊形區(qū)域的坐標信息用戶既可手工從設(shè)置框內(nèi)輸入各頂點坐標亦可從某一文件導(dǎo)入。坐標信息文件為一ai本文件，格式為“XYZ ae，其中坐標信息XYZ列須采用浮點數(shù)字格式表示且Z坐僅用于3d出項的數(shù)據(jù)輸出ae列為相應(yīng)空間點說明信息。用戶必須設(shè)定欲輸出數(shù)據(jù)空間點坐標信息對應(yīng)的正確的坐標系信息（Long/lat，UTM-32等。點坐標，亦可從某一文件導(dǎo)入。坐標信息文件為一ascii文本文件，格式為“XYZ項的數(shù)據(jù)輸出，Name列為相應(yīng)空間點說明信息。此外，用戶亦需設(shè)定所要輸出結(jié)果的起用戶必須設(shè)定欲輸出數(shù)據(jù)空間點坐標信息對應(yīng)的正確的坐標系信息（如Long/lat，UTM-32等。當用戶選擇采用斷面序列形式輸出計算結(jié)果時，需指定所需斷面的起始端點坐標。該坐標信息既可由用戶手工從設(shè)置框內(nèi)輸入兩端點亦可從某一文件導(dǎo)入坐標信息文件為一ai本文件，格式為“XYZ ae，其中坐信息YZ列必須采用浮點數(shù)字格式表示且Z標僅用于3d輸出項的數(shù)據(jù)輸出ae列為相應(yīng)空間點說明信息。此外如果用戶給定坐標信息文件中包含兩個以上的坐標信息則開始兩坐標點信息將被作為所需斷面的兩端點。hecrosconisdenedasaonoflentaeThefaceisincludedinthesectionwhenthelinebetweenthetwoelementcentersofthefacescrossesthelinebetweenthespecifiedfirstandlastpoint.用戶必須設(shè)定欲輸出數(shù)據(jù)空間點坐標信息對應(yīng)的正確的坐標系信息（Long/lat，UTM-32等。區(qū)域序列（標，亦可從某一文件導(dǎo)入。坐標信息文件為一ascii文本文件，格式為“XYZZ坐標，Name列為相應(yīng)空間點說明信息。用戶必須設(shè)定欲輸出數(shù)據(jù)空間點坐標信息對應(yīng)的正確的坐標系信息（如Long/lat，UTM-32等。2d變量(2DfieldSSC-Fraction1Bed ponentBed ponent ponentSuspended ponentTotalload,x-directionTotalload,y-directionRateofbedlevelchangeBedlevelchangeBedlevelBedload,magnitudeBedload,direction[degtrueSuspendedload,magnitudeSuspendeload,direction[degtrueTotalload,magnitudeTotalload,direction[degtrueAcc.bed ponentAcc.bed ponent Acc.total ponentAcc.total ponentCentrifugalaccelerationSurfaceelevationStillwaterdepthTotalwaterdepthDepthaverageU-velocityWaveheightWaveperiodWavedirection[degtrue質(zhì)量收支變量（Massbudget）Totalarea-totalvolume/masswithinWetarea-volume/massintheareawithinpolygonforwhichthedepthislargerthanthedryingRealwetarea-volume/massintheareawithinpolygonforwhichwaterdepthislargerthanthewettingDryarea-volume/massintheareawithinpolygonforwhichthedepthislessthanthedryingTransport-accumulatedvolume/masstransportedoverlallimitsSource-accumulatedvolume/massadded/remo