《水流对沙滩冲刷过程的数值模拟.docx》由会员分享，可在线阅读，更多相关《水流对沙滩冲刷过程的数值模拟.docx（20页珍藏版）》请在三一文库上搜索。
1、5.3水流对沙滩冲刷过程的数值模拟5.3.1 案例简介本案例是对水流冲刷沙滩过程的气固液三相流进行数值模拟,如图5-3-1所示,这是一个简化的二维模型,区域总长度为2000mm总高度为500mm下半部为一倾斜的沙子区域.水流从左上角的100mnW的进口流入,进入区域冲刷沙子,然后从右侧500mm高的出口流出.通过模拟,可清楚的看到水流对沙滩的冲刷过程,已经气固液三相的分布情况.图5-3-1水流冲刷模型5.3.2 Fluent求解计算设置1 .启动 Fluent-2D(1)双击桌面Fluent14.0 图标,进入启动界面.(2)选中Dimension 2D单项选择按钮,取消对 Display O
2、ptions下的三个复选按钮的选择.(3)其它保持默认设置即可,单击OK按钮进入Fluent 14.0 主界面窗口.2 .读入并检查网格(1)执行菜单栏中的File fReacRMesh命令,在弓t出的 Select File对话框中读入convection.msh 二维网格文件.(2)执行菜单栏中的MeshInfo 一Size命令,得到如图 5-3-3所示的模型网格信息：共有14342个节点,28411个网格面,14070个网格单元.(2)执行菜单栏中的MeshCheck命令.反应信息如图 5-3-4所示,可以看到计算域三维坐标的上下限,检查最小体积和最小面积是否为负数.!L puiplCe
3、llsFacFflKirlrsPar t i t i ons? nrll.91H070N8TH1* 取IN1zomrs, fi fie?-图5-3-3网格数量信息Donkin Lxtents:x-conrilinril p： nin ( m = flanidx (n =九 口皿口ft新中触y-CQCirdinate; niri (fi) = 0.000000+00B nw (n = 5*GOOBOte-fii Uolune statistic*ninimini uulunn (m3) : 3；,333241b-D5Hdxliuuki uuluiM? (3): 9 -037896(*-05 to
4、tal yoiurw心兀 Fir? statistics:ninlraudi Face arem(nZ): 3.33333 0-03nxlmuR Face area (m2); 1 pUt)UljtKde-02 匚heizki ng me!&h - - - - - - -lur*图5-3-4 Fluent 网格信息3 .求解器参数设置(1)单击选择左边workspace 中 Problem SetupfGeneral 命令,在出现的General面板中进行求解器的设置.(2) General面板中,开启重力加速度.单击选中Gravity 复选按钮,Y(m/s2)文本框输入-9.8 , Time
5、下选中Transient单项选择框,其它求解参数保持默认设置,如图 5-3-6所示.图5-3-6求解参数设置(4)单击选择Problem Setup Model命令,对求解模型进行设置,如图 5-3-7所示.(5)双击 Models fViscous-Laminar 选项(或选中 Viscous-Laminar ,然后点击 Edit ), 翻开 Viscous Model 窗口.(6)在弹出的窗口中, Model下选择k-epsilon(2 eqn)单项选择按钮,其他保持默认,如图5-3-8所示,单击OK按钮,启动k-e湍流方程.M 口 dels5-3-7计算模型选择图5-3-8湍流模型选择(
6、7)再次回到 Models 面板,双击 Multiphase-off 选项,翻开 Multiphase Model 窗口. 单击选中 Model 下 Eulerian 单项选择按钮, Number of Eulerian Phases 设置为 3,如图 5-3-9 所示,单击OK完成设置.图5-3-9多相流模型选择4 .定义材料物性(1)单击选择 Problem SetupMaterials命令,在出现的 Materials 面板中对所需材料进行设置,如图 5-3-10所示.(2)双击面板中的 Materials fFluid选项,弹出材料物性参数设置对话框,如图5-3-11所示.(3)材料物
7、性参数设置对话框中,单击 FLUENTDatabase按钮,翻开 FLUENTDatabase Materials 窗口,FLUENTFluid Materials 下选择 water-liquid(h2o) 选项,如图 5-3-12 所示,单击Copy,复制水的物性参数.(4)单击Change/Create按钮,保存对水物性参数的更改,单击 Close关闭窗口.(5) Name文本框输入 sand, Chemical Formula 清空,Density(kg/m3) 设置为 2500, Viscosity(kg/m-s) 设置为10,如图5-3-13所示,单击 Change/Create按
8、钮,完成设置.Mate-rielsMiLgns| .-d. ,t. F, l图5-3-10材料选择面板图5-3-11空气物性参数|HuH,ChmcairSmiua图5-3-12水物性参数设置图5-3-13沙子物性设置-luw Fluid Maw i北ui.#曲yid百厘勺；比屈h炳旭-hicHh fhne clih?rh) 3nylidBnEHbnd (ch2【二 2、料 z t F千nr wDDd-volatiles选项,翻开 Primary Phase 窗口,Name water-liquid 选项,如图 5-3-16 所示.选项,翻开 Primary Phase 窗口,Name sand
9、选项,选中 Granular 复选按钮,5.两相属性设置(1)单击选择Problem Setup Phases命令,在出现的 Phases面板中对所需材料进行设置,如图5-3-14所示.(2)双击 Phases 下 phase-1-Primary Phase 选项,翻开 Primary Phase 窗口,Name 中输入air , Phase Material 右侧下拉菜单项选择择 air选项,如图5-3-15所示.(3)双击 Phases 下 phase-2-Primary Phase 中输入water , Phase Material右侧下拉菜单项选择择(4)双击 Phases 下 pha
10、se-3-Primary Phase 中输入 sand , Phase Material右侧下拉菜单项选择择Properties 下 Diameter(m)设置为 0.000111 , Granular Viscosity(kg/m-s) 选择 gidaspow 选项,Granular Bulk Viscosity(kg/m-s)选择 lun-et-al 选项,Packing Limit 设置为 0.6 ,如图5-3-17所示,单击 OK(5)单击 Phases 面板中的 Interaction 按钮,弹出 Phase Interaction 窗口,Drag 选项卡下,相之间白拽力都设置为 s
11、chiller-naumann 选项,如图5-3-18所示,单击 OK完 成设置.Pfl-a-SEEEd J (iri；占白性身：p图5-3-15气相设置图5-3-16液相设置图5-3-17沙子设置图5-3-14气固液相设置面板图5-3-18三相间拖拽力设置6.边界条件设置(1)单击选择Problem Setup f Boundary ConditionsConditions 面板中对边界条件进行设置,如图 5-3-19所示.命令,在才T开的BoundaryBoundAiy CondhinkBdufsiutFELnrdeF-iut r：iM北u卜“丁 ElOutl Utl：CESflUftl0
12、图5-3-19边界选择(2)双击面板中的 Zone一in选项,弹出 Velocity inlet 都设置为0.01 ,如图5-3-20所示,单击 OK完成设置.窗口,湍动能和湍流耗散率图5-3-20进口边界条件设置(4) Boundary Conditions 面板中,Phase下方的下拉菜单项选择择 water选项,单击Edit 按钮,在弹出的窗口中,速度值设置为5m/s,如图5-3-21所示,单击选择 Multiphase 选项卡,设置进口水的体积分数为 1,如图5-3-22所示.图5-3-21进口速度设置图5-3-22进口水体积分数设置5.3.3求解计算1.求解限制参数(1)单击选择 S
13、olution Solution Methods 命令,在弹出的 Solution Methods 面板 中对求解限制参数进行设置.(2) Scheme下选择Coupled算法,其他保持默认设置,如图 5-3-23所示.2.求解松弛因子设置(1)单击选择 Solution Solution Controls 命令,在弓t出的Solution Controls面板中对求解松弛因子进行设置.(2)面板中 Flow Courant Number 设置为 40, Momentum口 Pressure 都设置为 0.5 , Volume Fraction 设置为 0.4 ,如图 5-3-24 所示.Sn
14、lLitiort Melhod-s；PrBUbV.M 厉 CmEidwmpCnie9 3Cb-?0CnritaLn-ol -blirmt fncinn |ar (Titw=5.9 = JDp-(11iWMMFUUEMT lJlO 111Lnniaj-s -zilVc t-me- frjrtnn arl Jime=i= JQDJp-Qil iST&fiLUFNT MG ?FiMMlI r-Mefll,：IMICT4M1 tOOMJi - fll1 HMIi tDU I MMil j：UU)l：04)lIIL色如Ml,T WE TMMI1B利河1 frJClt-L1 方加附 国一2AJe4i i硼
15、Ml1工四i I . * iCqFue qI ?Wul v lr: dm wihi：, Tirre=Q JDO. tjCLr:二口 UF 由用由口n产*学 ieH?班号TS FLU&JT IdOCh州目市FIXEMTHQ 口riMJoo图5-3-40初始液相云图图5-3-41 0.299s 时刻液相云图,髭*OJ vieMn a V*JJ *3+00 iO3XCjufr-ai I -301! 43ZUi 01*01H 2B Dda - i- H - I I3XM9 DCfrOl 色mn 5-OXH 丁3015-0fl 5DTdCnnicws i VcLms fr-j-zinn 力创七口口图5-
16、3-42 0.599s 时刻液相云图 图5-3-43 0.900s 时刻液相云图mnEZ1 MW 丁 5azi1ESQZH-5国1, 5-mrfi3iQZ1 m + Wfrfl14 I1SIMN1陋皿1WH 川M用即N15731?ia4i 1 m闻 1 san1 JDb-41 &HP43- 历H?CiMTfcajr . afVciijme hdi 归hdj iTep=i OCKK ?+flClWfYBFUJENT ill 11Cpnfcxrs alAti ume liacaon -and- 11ieh 工品011?1,MMBrSFUJEMTllQl图5-3-44初始固相云图图5-3-45 0.299s时刻固相云图I JOv-O It M* Ra iC-：17ntel4 B9HII* 5m* 第1XWI1 i5n4i典!WJI 5BhOI q*3XTI -l1HM1I?A /WKMObAOtauCn vSMH4Nirf9HJJIKT HQQ4.L：pntccrs nfz-mne fiaibon ；5-snd；iJCnnr?-r3 dYz-MT t Ir-iztDr sa-hdl 4Tlmc=9 JJDCh-011图5-3-46 0.599s时刻固相云图图5-3-47 0.900s时刻固相云图