基于Fluent的进气歧管空气动力学仿真

1、1仿真条件本课程报告使用Fluent对一进任歧管进行空气动力学仿真。该进气歧管的入口面积为 5026.55帆/，单个出I I面积为2552.67mm2,如图所示。根据四冲程6缸发动机的工作循环表，将进气压力分成56种工况。曲轴转角（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、气序号曲轴转角(deg )进气道末端的压力值(pa)第一缸第二缸第三缸第四缸第五缸第六缸1206800085000240100000850003608800085000480720005100660006120600007140850006800081608500010000091808500088000102007200011220660001224060000132606800085000142801000008500015300880008500016320720001734066000183606000019380850006800020400850001000002142085000

3、8800022440720002346066000244806000025500680008500026520100000850002754088000850002856072000295806600030600600003162085000680003264085000100000336608500088000346807200035700367206600060000(ed)妇*图单杠(第一缸)进气道末端压力值变化曲线120000100000800006000040000200000050100150200曲轴转角(deg)2仿真步赛基于以上条件对进气歧管进行空气动力学仿真。注进气歧管模型

4、导入到Gambit中。建立实体，进行布尔操作。获得进气歧管内的空间模型。对模型进行边角优化处理边角优化后的结果优化后对模型划分网格划分网格数165万设置边界条件将模型以mesh格式导出o在fluent中导入mesh模型EX检查网格并定义长度单位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.检查网格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.显小网格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设置监视器口 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»统计质量流量3结果统计进气歧管入口流量序号曲轴转角(deg )入口序号曲轮转角(deg )入口质量流体积流质量流体积流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进气歧管入口质量流量曲线曲轴转角(deg)进气歧管出口流量，详细数据请看附件序号曲轴转角(deg )第一缸质量流体积流kg/smA3/s1200.511450.41752400.05890.0483600.31280.25544800.50410.411551000.55310.451661200.59820.488371400.33850.276381600.36710.299791800.35210.2874第1缸进气歧管末端质量流量曲轴转角(deg)序号曲轴转角(deg )第二缸质量流体积流kg/smA3/s255000.52440.428126520-0.0002-0.0001275400.31760.2592285600.51550.4208295800.56610.4621306000.61270.5001316200.34050.278326400.37540.3064336600.36280.2962第2缸进气歧管末端质量流量曲轴转角(deg)序号曲轴转角(deg)第三缸质量流体积流kg