基于Fluent的進氣歧管空氣動力學仿真

1、1仿真條件本課程報告使用Fluent對一進任歧管進行空氣動力學仿真。該進氣歧管的入口面積為 5026.55帆/，單個出I I面積為2552.67mm2,如圖所示。根據(jù)四沖程6缸發(fā)動機的工作循環(huán)表，將進氣壓力分成56種工況。曲軸轉(zhuǎn)角（deg ）第一缸第二缸第三缸第四缸第五缸第六缸0-1800-60進氣壓縮作功進氣排氣作功60-120排氣壓縮120-180作功進氣180-360180-240壓縮排氣240-300進氣作功300-360排氣壓縮360-540360-420作功進氣420-480壓縮排氣480-540進氣作功540-720540-600排氣壓縮600-660作功進氣660-720壓縮排

2、氣序號曲軸轉(zhuǎn)角(deg )進氣道末端的壓力值(pa)第一缸第二缸第三缸第四缸第五缸第六缸1206800085000240100000850003608800085000480720005100660006120600007140850006800081608500010000091808500088000102007200011220660001224060000132606800085000142801000008500015300880008500016320720001734066000183606000019380850006800020400850001000002142085000

3、8800022440720002346066000244806000025500680008500026520100000850002754088000850002856072000295806600030600600003162085000680003264085000100000336608500088000346807200035700367206600060000(ed)婦*圖單杠(第一缸)進氣道末端壓力值變化曲線120000100000800006000040000200000050100150200曲軸轉(zhuǎn)角(deg)2仿真步賽基于以上條件對進氣歧管進行空氣動力學仿真。注進氣歧管模型

4、導入到Gambit中。建立實體，進行布爾操作。獲得進氣歧管內(nèi)的空間模型。對模型進行邊角優(yōu)化處理邊角優(yōu)化后的結(jié)果優(yōu)化后對模型劃分網(wǎng)格劃分網(wǎng)格數(shù)165萬設置邊界條件將模型以mesh格式導出o在fluent中導入mesh模型EX檢查網(wǎng)格并定義長度單位21r<TDESKTOP-1461J5E 3d, pors, lam；File Grid Define Snotcriols 9 interface 9 donains, 7onp<, default-i pressure out1 out2 out3 0Ut4 outS out6 wall Fluidparallpl, shell cond

5、utDune.檢查網(wǎng)格Checking cell children.Checking parent-child pointers (cell).Checking parent-child pointers (Face).Checking Checking Checking Checking Chpcking Checking Checking Checking Checking Checking Checking Checking Checking Checking Dune.port Parallel HepFile Grid Define Solve Adapt Surface Displ

6、ay Plot Report Parallel Helpthread consistency across compute nodes, thread layer counts, neighbor counts, exterior faces, intprfacp Facps. interface nodes. Face ids. node ids. corner nodes. interior cell pointers, exterior cell pointers. face count, node count. storage.顯小網(wǎng)格File (Che( Che( Che( Che(

7、 Che( Che( Che( Chp( ChP( Che( Che( Uhe( Che( Che( cne( Che( Chee Done.設置求解器S Parallel FLUENTDESKTOP-1481J5E 3d, pbns, lam- XFile GridDefine Solve Adapt Surface Display Plot Report Paralle- HelpCheckin Checkin Checkin Chpckin Checkin Checkin Checkin Checkin Checkin Checkin Checkin Checkin Checkin Ch

8、eckin Checkin Checkin ChpckinDone.|Models>Solver.Materials.Phases.Operating Conditions.Boundary Conditions.Periodic Conditions.Grid Interfaces.Dynamic Mesh>Mixing Planes.Multiphase.Energy.Viscous-Radiation.Species>Discrete Phase.Solidification & Melting.Acoustics.Turbo Topology.Injectio

9、ns.DTRM Rays.Custom field Functions.Profilec.<>Units.User-Defined>Parallel FLUENT(gFilo Grid DefineCheckingcell cCheckingpa»-cntChpckingparpntCheckingthreadCheckingthreadCheckingneighbCheckingexteriCheckingInterfCheckinginterfCheckingFace iCheckingnndp iCheckingcornerCheckinginteriChec

10、kingexteriCheckingFace cCheckingnode cCheckingstoragVelocity FunnulatiuiiQ Absolute Relative仃 Pressure BasedDensity BasedSpace2DSolver叵 SolverAxisymmetricAxisymmetric SwirlG 3Da Implicit Explicitlime。SteadyUnsteadyFormulationGradient Option( Green-Gauss Cell BasedGreen-Gauss Node Based" Least S

11、quares Cell Basedorouo FormulationDone ( Superficial Velocity Physical VelocityCancel Help使用k-£湍流模型Dune .6fl Pora落Evex邊界條件設置口 Parallel FLUENTDESKTOP-1481 J5Epbns, skeFile Grid Define Solve Adi閏Boundary ConditioiICXJ Helpfluid parallelv shell conduction z Done.656 triangular wal 656 triangular w

12、al 656 triangular wal 656 tridnyuldr wdl1286 triangular pre 652 triangular pre 656 triangular pre 285067 nodes, binary. 285867 node Flags, bidistributing me parts. Faces.ZoneTypedcfault-intcriur|吧out2out3out4out5out6inletvent intake-fan interface mass-tlow-inlet outflow outlet-vent pressure-far-fiel

13、d pressure-inlet八Ipreccure_inlet.1 wallPreccure-outlet巴.symmetry velocity-inlet wallVy IDI9nodes.Set.Copy.CluseHelp岡 Pressure OutletZone Name lout 丁Momentum Thermal Radiation | Species | DPM Multiphase UDSGauge Pressure (pascal)卜33325|constantBackflow DirectionSpecification MethodNormal to BoundaryR

14、adial Equilibrium Pressure DistributionTarget Mass Flow RateTurbulenceSpecification Method |K and Epsilon二jBackflow Turbulent Kinetic Energy (mZ/sZ),p:onstant臼Backflow Turbulent Dissipation Rate (m2/s3) 1constant臼OK | Cancel | Help |求解初始化 Parallel RUENT©DESKTOP 1481J5E 3d. pbns. skoFile GridOUt

15、1UUt4FluishellDone.6566566566561286652656285067285067設置監(jiān)視器口 Print R PlotIterations1000 芻NormalizationNormalise 廠 ScaleoutoutOptionsS Residual MonitorsPlottingWindowIterationsawparall ShellDone.656 tr656 tr 656 tr656 tr 1286 tr652 trAxes. I Curves. I656 tr 285167 no 285167 nodisCheckAbsolute上Residual

16、Monitor Convergence Criteriacontinuity |7|70.001卜-velocity|7“0.001|y-uelocity 0.001|z-uclocity |7|o. 001IkI7I7Fo.00111HConvergence CriterionabsoluteOKPlot Renorm | Cancel | Help求解計算Parallel FLUFNTDF?<TOP-1481 J5F 2d. pbn， sLpFile Grid Define Solve Adopt Surface Display Plot Report Parallel HelpDo

17、ne.656 triangular wall Faces. zone U. binary.h；UDFFIterate I Apply I Close I Help656 triangular656 triangular656 triangular1286 triangular652 triangular656 triangular285067 nudeb, binary285067 nodp Flags, bdistributing n parts, faces.nodes.t：1 ：23iter continuity x-uelocity y-uelocity z-uelocityk eps

18、ilon1 1 .O06Oe*OO 6.2795e*01 3.01 Oe-05 J1.8797e-01 6.782 0c-05 9.5728e*0O求解完成File Grid Defire Solve Adapt Surface Display Plot Report Paralel Helpturbulent uiscosity United to viscosity ratio of1in 10 cells936 1.0030e-632.6We-05 1.1623e-05 1.1193e-051.5259e-05 2.845Oe-056:04reuersed flow in30Faces

19、on prcssurc-outlot 6.reversed flow in29Faces on pressure-outlet ?.turhiilpnt uiscosity931 1.0002e-032.6H169-05 1.16109-05 1.1193»-05reuersed fluu In30reuersed flow in29faces on pressure-outlet 9.United to viscosity ratio of1.06000Oc>0ff； in 9 cells1.519Ze-Ub 2.«62be-UbXlimitpri to uiscn

20、sity ratio offaces un pressure-uuClet 6. |932 bulutiuii is cunueryed932 9.9/62e-U4 2.63/3P-0S 1 .1b91e-Ub 1.119SP-U5i.dQnnnnp*nns in 9 cells1 .52150-05 2.8361ig-05 G:0Ii運 Parallel FLUEN-eDESKTOP-1481 J5E 3d, pbn( eketurbulent viscosityScmedRe&iQjAisJun 29.2017FLUENT 6.3 OU PDflf.lW繪制速度矢量繪制壓力分布圖粉

21、 Paral d FLUEN aDES<TOP-1481 JSE Fluent Inc3 86103270bs 3 32e*03 -6.019*03 1 05e*04 1 41©4Q4 1 "aae-213et04 2 49e*Gfl 2 88Y4-3 2(te 加- 3 5664g- 3 92e-*04- 4 28e心- 4 64e*C4- 5 00e*0 fl- 5心如- 5 72e-0<s- 6 08e*0<- 6 4$e皿-6 80etQ4veiocvvecfots coioreo fcv BtsucPreoeure ipascaDJun 29.

22、FLUERfm 1 pi)n»統(tǒng)計質(zhì)量流量3結(jié)果統(tǒng)計進氣歧管入口流量序號曲軸轉(zhuǎn)角(deg )入口序號曲輪轉(zhuǎn)角(deg )入口質(zhì)量流體積流質(zhì)量流體積流kg/smA3/skg/smA3/s1200.890.7266193800.86450.70572400.44520.3635204000.40490.33063600.69210.565214200.6880.56164800.50410.4115224400.52280.426751000.55310.4516234600.57380.468461200.59820.4883244800.62080.506871400.88440.7

23、219255000.88930.72681600.46170.3769265200.38720.316191800.69490.5672275400.69330.566102000.52080.4252285600.51550.4208112200.57010.4654295800.56610.4621122400.61880.5051306000.61270.5001132600.87550.7147316200.88170.7198142800.3650.298326400.43360.3539153000.68690.5607336600.69810.5699163200.52180.4

24、26346800.51970.4243173400.57370.4683357000.57170.4667183600.62140.5072367200.620.506進氣歧管入口質(zhì)量流量曲線曲軸轉(zhuǎn)角(deg)進氣歧管出口流量，詳細數(shù)據(jù)請看附件序號曲軸轉(zhuǎn)角(deg )第一缸質(zhì)量流體積流kg/smA3/s1200.511450.41752400.05890.0483600.31280.25544800.50410.411551000.55310.451661200.59820.488371400.33850.276381600.36710.299791800.35210.2874第1缸進氣歧管末端質(zhì)量流量曲軸轉(zhuǎn)角(deg)序號曲軸轉(zhuǎn)角(deg )第二缸質(zhì)量流體積流kg/smA3/s255000.52440.428126520-0.0002-0.0001275400.31760.2592285600.51550.4208295800.56610.4621306000.61270.5001316200.34050.278326400.37540.3064336600.36280.2962第2缸進氣歧管末端質(zhì)量流量曲軸轉(zhuǎn)角(deg)序號曲軸轉(zhuǎn)角(deg)第三缸質(zhì)量流體積流kg