關(guān)于本文檔（Aboutthis 編寫本文檔的目的 用戶所需的知識(shí)背景(AssumedUser MIKEZero模型編輯器程序界面結(jié)構(gòu)(GeneralEditor 系統(tǒng)(Online 引言 概述(General 說明(Getting 示例 概述 U型河道示例(U-shaped 環(huán)流示例(Sandtransportincirculation 沙模塊（SandTransport 模型定義（Model 概述（General 波流共同作用工況(Combinedcurrentand 時(shí)間參數(shù)（Time 求解格式（Solution 泥沙特性（Sediment mended 床面阻力（Bed 概述（General mended 驅(qū)動(dòng)力參數(shù) 擴(kuò)散 對(duì)流擴(kuò)散（Vertical mended 源匯項(xiàng) 源匯項(xiàng)設(shè)定（Source 備注與提示（Remarksand 初始條件（Initial 泥沙粒徑組分濃度（Fraction 床面活動(dòng)層厚度（Layer 邊界條件（Boundary 邊界設(shè)置（Boundary 地貌演化模塊（Morphology 輸出 圖形化查看（Geographic 輸出類型（Output 輸出項(xiàng)（Output 地貌演化模塊（Morphology 概述（General 模型定義（Model 時(shí)間參數(shù)（Time 河岸侵蝕（Bank 邊界條件（Boundary 邊界說明（Boundary 關(guān)于本文檔（Aboutthis編寫本文檔的目的MIKE3FlowModelFM（非粘性泥沙輸移的計(jì)算）用戶所需的知識(shí)背景(AssumedUser文件及其編輯器、繪圖工具（PlotComposer）、MikeZeroToolbox、數(shù)據(jù)查看工具（DataViewer）以及網(wǎng)格生成工具（MeshGenerator）Mike3基本組件的說明文C:\ProgramFiles\DHI\MIKEZero\Manuals\MIKE_ZERO\MzGeneric.pdf。MIKEZero模型編輯器程序界面結(jié)構(gòu)(GeneralEditor編輯窗體(Editor供了一個(gè)文本菜單項(xiàng)供用戶選擇是否顯示圖形化的地形、網(wǎng)格、GIS背景圖層或圖例；從系統(tǒng)(Online調(diào)出系統(tǒng)的方法有多種，用戶可按照個(gè)人需要選擇任意法使用在幫助頁面中手工打開系在幫助菜單中，選中“helptopics”菜單項(xiàng)即可調(diào)出系統(tǒng)圖1.1MIKE3FlowModelFM系引言Mike3FlowModelFM沙模塊主要模擬單一波浪、水流或波流共同作用下沙的侵蝕、沙模塊主要用于研究非粘性泥沙輸移的工程問題說明(GettingMIKE3FlowModelFM模型的水動(dòng)力模塊計(jì)用戶手冊(cè)和包括水動(dòng)力模塊的內(nèi)容翔實(shí)的“跟我學(xué)”教學(xué)指南。PDF格式的用戶手冊(cè)和教學(xué)指南文檔一般于：C:\ProgramFiles\DHI\MIKEZero\Manuals\MIKE_3\C:\ProgramFiles\DHI\MIKEZero\Manuals\MIKE_3\FlowModel_FM\DHI型。這兩份文檔的PDF格式文件一般于：C:\ProgramFiles\DHI\MIKEZero\Manuals\MIKE_3\C:\ProgramFiles\DHI\MIKEZero\Manuals\MIKE_3\FlowModel_FM\示例MIKE3FlowModelFM模型的最好方法之一。因此模型安裝程序中：U型河道示例C:\Program環(huán)流示例C:\Program環(huán)流示例(Sandtransportincirculation沙模塊（SandTransport移問題是非常必要的（限泥沙補(bǔ)給和高床面抗沖性條件）。 正比于二者比值的平方：5.1

3d泥沙組分定義（FractionEngelundandHansen(全沙VanRijn(推移質(zhì)+懸移質(zhì)輸沙EngelundandFredsoe(推移質(zhì)+懸移質(zhì)輸沙Meyer-PeterandMuller(推移質(zhì)輸沙推移質(zhì)(Bed中則包含一推移質(zhì)系 ，該中則包含一推移質(zhì)系 ，該10.52.0計(jì)算結(jié)果超計(jì)算結(jié)果超 的適用范圍懸移質(zhì)（Suspended當(dāng)用戶設(shè)定模型中存在懸移質(zhì)輸沙時(shí)，相應(yīng)輸沙中則包含一懸移質(zhì)系 ，該1，0.52.0。注意該參數(shù)取值超出推薦取值范圍時(shí)可能會(huì)造成計(jì)算結(jié)果超出的適用范圍。全沙（Total因EngelundandHansen僅計(jì)算全沙輸沙率，無法直接得出模型計(jì)算所需的推移 分別代表懸移質(zhì)、推移質(zhì)在全沙中所占的比例，如設(shè) 分別代表懸移質(zhì)、推移質(zhì)在全沙中所占的比例，如設(shè) umConcentration）所給定的“泥沙輸移表”（Sedimenttransporttable）線性插值的結(jié)果。 綴相同，后綴分別為“.lon”和“.crs”的ASCII碼格式文件。泥沙輸移表（SedimenttransportMike21Toolbox中的“GenerationofQ3DSedimentTables”組件生成的。用戶在使用該組件生成“泥沙輸移表”文件時(shí)可以綜合考備注與提示（Remarksand戶所給定的“泥沙輸移表”（Sedimenttransporttable）線性插值的結(jié)果，故當(dāng)用戶所生log文件中。5.1VS--d--g--泥沙輸移計(jì)算的開始模型時(shí)間為用戶在時(shí)間設(shè)置框中設(shè)定的開始時(shí)間泥沙輸移計(jì)算的時(shí)間步長(zhǎng)為時(shí)間設(shè)置框中設(shè)定的主時(shí)間步長(zhǎng)（overalltimestep）乘以時(shí)間倍乘系數(shù)（timestepfactor）。CFL數(shù)(27頁)1，為保證所有網(wǎng)格點(diǎn)CFL數(shù)均滿足該限制條件，模型中時(shí)間步長(zhǎng)的取值采用一浮動(dòng)范圍的方式，因此模型中用用戶可在HydrodynamicModel的求解格式框中進(jìn)行最小和最大時(shí)間步長(zhǎng)范圍以CFL數(shù)的設(shè)置。備注與提示（Remarksand1?23-4倍。一般來說采用高階CFL1CFL1時(shí)，模型即可保持計(jì)算CFL數(shù)的數(shù)值為近似預(yù)估值，故在這種默認(rèn)設(shè)置情況下仍CFL值適當(dāng)減?。?1之間），此外用戶亦可適當(dāng)減小所設(shè)定的最大時(shí)間步長(zhǎng)。必須，當(dāng)用戶將最小和最大時(shí)間步長(zhǎng)均設(shè)定為與主時(shí)間步長(zhǎng)相同時(shí)(p23)，模型將CFL1。泥沙輸移計(jì)算需要用戶事先設(shè)定模擬區(qū)域的泥沙特性d50當(dāng)選擇采用波流共同作用下泥沙輸移模式計(jì)算時(shí)，用戶還需設(shè)定床沙幾何標(biāo)準(zhǔn)方差常數(shù)（空間均一d50（mm單位）和幾何標(biāo)準(zhǔn)方差σg2d網(wǎng)格化數(shù)據(jù)文件（dfs2）2d（dfsu mended沙模塊主要基于當(dāng)前相對(duì)較為成非粘性泥沙輸移理論進(jìn)行泥沙輸移計(jì)算。一般而言，0.06mm2mm。0.30.7之間。2.65。床面阻力（Bed無床面阻力（光滑壁面系數(shù)（Chezy系數(shù)(Manning動(dòng)床阻力 水動(dòng)力模型中生成的阻力系數(shù) fromHD當(dāng)選擇采用“動(dòng)床阻力”形式設(shè)定床面阻力時(shí)，用戶需基于系數(shù)（m1/3/s）的形動(dòng)床阻力（Alluvial 隨水深變化而變化。系數(shù)（C）或系數(shù)可由下式給出：式中，Ha、b增大而相對(duì)偏大。一般而言，thisreducesovershooteffectswhenbedscourcommences,andcausesthetopographyofapointbarcrestto emorerounded(asreportedbyTalmon,當(dāng)采用動(dòng)床阻力的形式設(shè)定床面阻力時(shí)，用戶需設(shè)定系數(shù)取值的上下限。當(dāng)基于式5.2計(jì)算所得系數(shù)超過用戶規(guī)定的上下限時(shí)，計(jì)算中相應(yīng)取值實(shí)際上為所對(duì)應(yīng)如果水動(dòng)力模型的率定結(jié)果表明，系數(shù)與局地水深表現(xiàn)為一明顯的對(duì)應(yīng)關(guān)C（m1/2/s常數(shù)（空間均一（dfsudfs2文件時(shí)，相應(yīng)計(jì)算網(wǎng)格上的泥沙特性參數(shù)基于雙線性插值方法生成。 mended為避免因引入動(dòng)床阻力而引起模型計(jì)算的失穩(wěn)問題，最大系數(shù)取值應(yīng)盡可能不大無波浪條件（Nowaves）(Hrms)min(Tp)min插值計(jì)算求得，相應(yīng)輸沙設(shè)定為0度。用戶設(shè)定波浪條件（Wavefield）常數(shù)時(shí)間上變化，空間上不變（Varyingintime,constant 時(shí)間上不變，空間上變化(Constantintime,varying 時(shí)間上、空間上均變化(Varyingintime 時(shí)間上變化，空間上不變（Varyingintime,constant （dfs05.1時(shí)間上不變，空間上變化（Constantintime,varying （dfsu時(shí)間上、空間均變化（Varyingintime 這種情況下，用戶需準(zhǔn)備一空間上至少完整包括模擬區(qū)域范圍且包含相關(guān)參數(shù)的2d網(wǎng)格化數(shù)據(jù)文件（）或2d無結(jié)構(gòu)數(shù)據(jù)文件（文件時(shí)，當(dāng)采用文件時(shí)，相應(yīng)計(jì)算網(wǎng)格上的泥沙特性參數(shù)基于分段常數(shù)插值方法生成；當(dāng)采用2文件時(shí)，相應(yīng)計(jì)算網(wǎng)格上的泥沙特性參數(shù)基于雙線性插值方法生成。當(dāng)或u型計(jì)算的時(shí)間步長(zhǎng)不一致時(shí)，則程序?qū)⒆詣?dòng)進(jìn)行相應(yīng)插值處理（見圖5.1圖5.1用戶需準(zhǔn)備一空間上至少完整包括模擬區(qū)域范圍且包含相關(guān)參數(shù)的2d（dfsu尤為顯著。因此，擴(kuò)散系數(shù)在形式上應(yīng)遵從這種關(guān)系。Especiallyinthehorizontaldirectionstheeffectsofnon-resolvedprocessescanbesignificant,inwhichcasethedispersioncoefficientformallyshoulddependonthe平流擴(kuò)散(Horizontal無擴(kuò)散(No擴(kuò)散系數(shù)(Dispersioncoefficient渦粘系數(shù)類比(Scalededdyviscosity當(dāng)采用擴(kuò)散系數(shù)時(shí)，用戶需給定擴(kuò)散系數(shù)的數(shù)值（單位m2/s）當(dāng)采用渦粘系數(shù)類比時(shí)，水流方程求解時(shí)所需的擴(kuò)散系數(shù)通過將渦粘系數(shù)乘常數(shù)（時(shí)間、空間上定常）(Constant(inbothtime 空間上變化(Varying 比例系數(shù)(Scalingfactor)對(duì)流擴(kuò)散問題的相關(guān)設(shè)置與平流擴(kuò)散問題相關(guān)設(shè)置完全相同，可參考5.7.1 mended當(dāng)基于渦粘系數(shù)類比（如Smagorinsky或k-ε模型）的方法進(jìn)行擴(kuò)散系數(shù)設(shè)置時(shí)，用這里提到的比例因子可以由數(shù)的倒數(shù) 得出。程序中數(shù)的準(zhǔn)確的直接預(yù)估出該參數(shù)的合理取值是較為的。但由雷諾相似準(zhǔn)則可知，擴(kuò)散系數(shù)可（Rodi，1980）匯點(diǎn)然后單擊“goto…”按鈕進(jìn)行相關(guān)設(shè)置。源匯項(xiàng)設(shè)定（Sourcespecification）量在水動(dòng)力模塊的源匯 框中設(shè)置量在水動(dòng)力模塊的源匯 框中設(shè)置流入水體）時(shí)，沙模塊中源匯項(xiàng)的懸沙濃度取值為用戶設(shè)定值；當(dāng)用戶設(shè)定的源強(qiáng)通量為負(fù)（水流流出水體）時(shí)，沙模塊中源匯項(xiàng)的懸沙濃度取值為源匯點(diǎn)處的水體點(diǎn)源中的懸沙濃度是獨(dú)立于水體的。強(qiáng)通量為負(fù)（水流流出水體）時(shí)，源匯項(xiàng)的懸沙濃度取值等于指定增加量加上點(diǎn)源處的水體懸沙濃度值，這種情況下的點(diǎn)源實(shí)際上是一個(gè)孤立的入流點(diǎn)；當(dāng)用戶設(shè)定的源強(qiáng)通量為負(fù)（水流流出水體）時(shí)，沙模塊中源匯項(xiàng)的懸沙濃度取值為源匯點(diǎn)處的源匯項(xiàng)相關(guān)參數(shù)數(shù)值的給定可以采用兩種方式常數(shù)（時(shí)間上不變備注與提示（Remarksand常數(shù)（空間上均一邊界條件（Boundary搭建模型時(shí)，模型配置編輯器(set-upeditor)首先會(huì)檢查網(wǎng)格文件，自動(dòng)找出網(wǎng)格文件中的邊界，并給每個(gè)邊界預(yù)先指定一個(gè)名稱，用戶可自行在框中對(duì)各個(gè)邊界的然后單擊“goto…”邊界設(shè)置（BoundarySpecification）陸地邊界平衡輸沙條件（Euilibrium用戶設(shè)定取值(Specified零梯度（Zerogradient）常數(shù)（時(shí)間上、空間上不變時(shí)間上變化、空間上（邊界上）線性插值分段立方插值(picewise正常順序逆序(reverse(Softstart常為海量數(shù)據(jù)，故不加選擇地整個(gè)模擬范圍內(nèi)所有時(shí)間段上的輸出結(jié)果是不現(xiàn)實(shí)的。以單擊”View”MIKEZEROVIEWING/EDITING工具查看相關(guān)內(nèi)容。結(jié)果類型（Field3d2d3d形式輸出。輸出格式（Ouput2d3ddfs0格式。dfs0格式。5.2TimeProfileDataDataDataData輸出文件（Output是否輸出干濕點(diǎn)（Treatmentoffloodand對(duì)于干點(diǎn)和濕點(diǎn)位置的2d和3d計(jì)算結(jié)果，用戶有三種輸出0。點(diǎn)序列（Pointseries）件導(dǎo)入。坐標(biāo)信息文件為一ascii文本文件，格式為“XYZ Name，Y、Z列必須采用浮點(diǎn)數(shù)字格式表示，且Z坐標(biāo)僅用于3d輸出項(xiàng)的數(shù)據(jù)輸出，Name列為相32等。據(jù)其所處網(wǎng)格單元頂點(diǎn)處數(shù)值進(jìn)行線性插值的結(jié)果，網(wǎng)格單元各頂點(diǎn)處數(shù)值基于pseudo-log文件中。線序列（Line文件導(dǎo)入。坐標(biāo)信息文件為一ascii文本文件，格式為“XYZName，其中坐標(biāo)信息32等。據(jù)其所處網(wǎng)格單元頂點(diǎn)處數(shù)值進(jìn)行線性插值的結(jié)果，網(wǎng)格單元各頂點(diǎn)處數(shù)值基于pseudo-log文件中。面序列（AreaName32等。體積序列（Volume頂點(diǎn)坐標(biāo)，亦可從某一文件導(dǎo)入。坐標(biāo)信息文件為一ascii文本文件，格式為“XYZName32等。斷面序列（CrosssectionName，斷面的兩端點(diǎn)。Thecrosssectionisdefinedasasectionofelementfaces.Thefaceisincludedinthesectionwhenthelinebetweenthetwoelementcentersofthefacescrossesthelinebetweenthespecifiedfirstandlastpoint.32等。 坐標(biāo)亦可從某一文件導(dǎo)入坐標(biāo)信息文件為一ascii文本文件格式“XYZ Name其中坐標(biāo)信息X、Y、Z列必須采用浮點(diǎn)數(shù)字格式表示，但本輸出類型中并不使用Z坐標(biāo)，Name列為相應(yīng)空間點(diǎn)說明信息。32等。輸出項(xiàng)（Output2d變量(2DfieldSSC-Fraction1Bed ponentBed ponentSuspended ponentSuspended ponentTotalload,y-directionRateofbedlevelchangeBedlevelchangeBedlevel5個(gè)基本變量輸出。Bedload,magnitudeBedload,direction[degtrueSuspendedload,magnitudeSuspendeload,direction[degtrueTotalload,magnitudeTotalload,direction[degtrueAcc.bed ponentAcc.bed ponentAcc.suspended ponentAcc.suspended ponentAcc.total ponentAcc.total ponentDeviationCentrifugalaccelerationSurfaceelevationStillwaterdepthDepthaverageU-velocityDepthaverageV-velocityWaveheightWaveperiodWavedirection[degtrue0度，為波浪傳來的方向。質(zhì)量收支變量（Massbudget）Totalarea-totalvolume/masswithinWetarea-volume/massintheareawithinpolygonforwhichthedepthislargerthanthedryingRealwetarea-volume/massintheareawithinpolygonforwhichwaterdepthislargerthanthewettingDryarea-volume/massintheareawithinpolygonforwhichthedepthislessthanthedryingTransport-accumulatedvolume/masstransportedover llimitsSource-accumulatedvolume/massadded/remove