FLUENT算例 (15)
- 格式:doc
- 大小:161.00 KB
- 文档页数:5
fluent经典算例Tutorial:Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel CellIntroductionThe purpose of this tutorial is to provide guidelines and recommendations for using the PEM Fuel Cell add-on module available in FLUENT.This tutorial demonstrates how to do the following:Build and generate a mesh for a single-channel PEM fuel cell.Assign zone names and types that are required in the FLUENT PEM fuel cell add-onmodule.Import the mesh into FLUENT,set up the case,and run the calculations.Postprocess the results.PrerequisitesThis tutorial assumes you are familiar with both the GAMBIT and FLUENT interfaces and that you have a good understanding of the basic setup and solution procedures.Details of most steps will not be shown explicitly.If you have not used FLUENT before,it would be helpful to?rst refer to the FLUENT6.3 User’s Guide and theFLUENT6.3Tutorial Guide.You can also refer to the GAMBIT User’s Guide,GAMBIT Tutorial Guide,and GAMBIT Modeling Guide.Problem DescriptionA schematic description of a single,straight-channel,counter-?ow PEM fuel cell is shown inFigure1.The dimensions of the fuel cell are2.4×2.88×125mm in the x,y,and z directions, respectively.This model represents a repeating channel of a larger counter-?ow PEM fuel cell.The cross-sectional area of the membrane electrode assemblyis2.4×125=300mm2.All the zones as marked in Figure1(with exception of the optional cooling channel)must be de?ned as separate zones within GAMBIT.Note that there can be multiples of each zone.For example,there can be multiple cathode gas di?usion layers,each with di?erent material properties.In FLUENT6.3,all continuum zones must be of the?uid type except for the current collectors which may be either solid or? uid.For simplicity,it is recommended that you make the current collectors solid continuum zones.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel CellFigure1:PEM Fuel Cell ZonesFor details about the FLUENT PEMFC model,refer to o f FLUENT6.3Fuel Cell ModulesManual.Note:If you would like to bypass all GAMBIT steps and proceed to the FLUENT portion of the tutorial,skip to Step. Preparation in GAMBIT1.Create a folder called pem-single-channel to contain the?les generated in thistutorial.2.Start GAMBIT and specify this folder as the working folder.To simplify the geometry construction and meshing steps,a journal?le,pem-single-channel.jou is provided.This?le contains GAMBIT instructions which will create the geometry and gen-erate the mesh.It is recommended that you step through the journal?le to understand eachof the steps and recognize the assigned zone boundary types pertinent to fuel cell problems.Geometry Creation and Mesh GenerationStep1:Geometry Creation1.Create a rectangular face primitive in the xy plane.Operation?→Geometry?→Face?→Create Real Rectangular Face(a)Enter2.4for Width and1.2for Height.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel Cell(b)Select+X+Y for Direction,(c)Click Apply and close the Create Real Rectangular Face panel.2.Create another face by sweeping the uppermost edge by0.21units in the y-direction.Operation?→Geometry?→Face?→Sweep Edges3.Create another face by sweeping the uppermost edge by0.012units in the y-direction.4.Create another face by sweeping the uppermost edge by0.036units in the y-direction.5.Create another face by sweeping the uppermost edge by0.012units in the y-direction.6.Create another face by sweeping the uppermost edge by0.21units in the y-direction.7.Create another face by sweeping the uppermost edge by1.2units in the y-direction.8.Create a rectangular face with dimensions(x,y)=(0.8,0.6).Use+x+y for the di-rection(so that the global coordinate system origin is in the lower left corner of theface).9.Move the newly created face using a vector of(0.8,0.6,0).Operation?→Geometry?→Face?→Move/Copy Faces10.Copy the face by1.08units in the y-direction.11.Split the lowermost face with the face you just created.Operation?→Geometry?→Face?→Split Face12.Similarly,split the uppermost face with its internal face.The resulting mesh should appear as shown in Figure2.Step2:Mesh Generation(Manual)1.Mesh the edges as shown in Figure2.The numbers of cells along each edge are indicated below.The geometry and mesh are symmetric about its horizontal midplane.Operation?→Mesh?→Edge?→Mesh Edges2.Mesh the nine faces using the Quad Submap scheme.Operation?→Mesh?→Face?→Mesh Faces3.Create an edge by sweeping any one of the vertices by125units in the positive z direction.Operation?→Geometry?→Edge?→Sweep VerticesModeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel CellFigure2:Edge and Face Mesh(a)Select Vector for Path and click the De?ne button.(b)Enable Magnitude and enter125.(c)Select Positive for Z from the Direction list.(d)Click Apply in the Vector De?nition form.(e)Click Apply and close the Sweep Vertices form.4.Mesh the newly created edge using a double-sided graded edge mesh that consists of 60elements.(a)Enable Double sided.(b)Enter1.1for Ratio1and Ratio2.(c)Click Apply and close the Mesh Edges form.5.Create volumes by sweeping the nine faces along the newly created edge. Operation?→Geometry?→Volume?→Sweep Faces(a)Enable With Mesh.(b)Click Apply and close the Sweep Faces form.The volume mesh is as shown in Figure3.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel CellFigure3:Volume MeshStep3:Zone Assignments and Mesh ExportFLUENT’s PEM fuel cell add-on module requires that boundary and continuum zones berigorously de?ned.Care should be taken in order to have a logical system of naming eachzone to represent each of the regions shown in Figure1.At a minimum,the boundary zones that are required include the following:Inlet and outlet zones for the anode gas channel.Inlet and outlet zones for the cathode gas channel.Surfaces representing anode and cathode terminals.Optional boundary zones that could be de?ned include any voltage jump surfaces,interiorow surfaces,or non-conformal interfaces that are required.The following continuum zones are also required:Flow channels for anode-and cathode-sideuids.Anode and cathode current collectors.Anode and cathode gas diusion layers.Anode and cathode catalyst layers.Electrolyte membrane.The inlets should all be assigned the boundary zone type MASS FLOW INLET and outletsshould be assigned the PRESSURE OUTLET type.The terminals are the regions where thevoltage(or current?ux density)is known.Normally,the anode is grounded(V=0)andModeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel Cell the cathode terminal is at a?xed potential that is less than the open-circuit potential.Bothof the terminals should be assigned the WALL boundary type.Voltage jump zones can optionally be placed between the various components(such asbetween the gas di?usion layer and the current collector).Faces which represent?uid/solidinterfaces must be of type WALL.Additional interior zones may be de?ned for purposes of post-processing.If such interiorzones are de?ned,they should have no e?ect on the solution.FLUENT’s PEM add-on module supports the use of non-conformal grid interfaces.In such cases,it is recommended that the INTERFACE boundary type be used halfway between the membrane continuum zone,and that the mesh on opposite sides of the interface have similar size,aspect ratio,and orientation.In such cases,the membrane will consist of two uid zones.1.Assign boundary zones according to the de?nitions listed in Table1.Operation?→Zones?→Specify Boundary TypesTable1:Boundary Zone AssignmentsAnode-side inlet(z=125,upper)inlet-a MASS FLOW INLETCathode-side inlet(z=0,lower)inlet-c MASS FLOW INLETAnode-side outlet(z=0,upper)outlet-a PRESSURE OUTLETCathode-side outlet(z=125,lower)outlet-c PRESSURE OUTLETAnode terminal(y=2.88)wall-terminal-a WALLCathode terminal(y=0)wall-terminal-c WALLAnode-side?ow channel walls wall-ch-a WALLCathode-side?ow channel walls wall-ch-c WALLFuel cell ends wall-ends WALLAnode-side di?usion layer walls wall-gdl-a WALLCathode-side di?usion layer boundaries wall-gdl-c WALLLateral boundaries of the fuel cell wall-sides WALL2.Assign continuum zones according to the de?nitions listed in Table2.Refer to Figure1 when selecting which volumes to assign to each zone.Operation?→Zones?→Specify Continuum TypesModeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel Cell Table2:Continuum Zone AssignmentsAnode-side catalyst layer catalyst-a FLUIDCathode-side catalyst layer catalyst-c FLUIDAnode-side?ow channel channel-a FLUIDCathode-side?ow channel channel-c FLUIDAnode-side gas di?usion layer gdl-a FLUIDCathode-side gas di?usion layer gdl-c FLUIDElectrolyte membrane membrane FLUIDAnode current collector current-a SOLIDCathode current collector current-c SOLID3.Export the mesh?le as pem-single-channel.msh.Setup and Solution in FLUENTFLUENT’s PEM Fuel Cell(PEMFC)model is provided as an add-on module with thestandard FLUENT licensed software.A special license is required in order to use thismodel.The module is installed as part of the standard FLUENT installation in the folder/addons/fuelcells2.2within the FLUENT installation folder.The PEMFC model con-sists of a user-de?ned function(UDF)library and a compiled Scheme library which can beloaded using a text user interface(TUI)command.trtitlePreparation1.Copy the mesh?le pem-single-channel.msh.gz to your working folder.If you worked through the GAMBIT portion of this tutorial,an uncompressed versionof this?le will already be in place.2.Start the3DDP(3ddp)version of FLUENT.!If you wish to solve this case in parallel,you will need to set up and save the case?le in serial mode?rst.Once this is done,you can start a parallel FLUENTsession and proceed with the calculations.Step1:Grid1.Read the mesh?le,pem-single-channel.msh.gz.File?→Read?→Case...FLUENT will perform various checks on the mesh and will report the progress in theconsole window.Make sure that the reported minimum volume is a positive number.2.Check the grid.Grid?→CheckA grid check should always be performed in order to verify the integrity of the meshle.Specically,you should verify that the minimum cell volume is a positive value.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel Cell3.Since the mesh was created in units of millimeters,it must be scaled.Grid?→Scale...(a)Select mm from the Grid Was Created In drop-down list.(b)Click the Change Length Units button.(c)Click Scale.(d)Verify that the value of Zmax is125mm and close the Scale Grid panel.Step2:Models1.To load the PEMFC model using the text user interface,enter the following commandin the console:/define/models/addon-module3This command will load a Scheme library which contains the PEM model GUI anda UDF library.Upon successful execution,the following message will be displayed inthe console:Addon Module:fuelcells2.2...loaded!2.Calculate the surface area of the membrane for post-processing.In this case,the membrane area is equal to the surface area of the cathode terminal. This surface is named wall-terminal-c.Report?→Projected Areas...(a)Select Y from the Projection Direction group box.(b)Select wall-terminal-c from the Surfaces list.(c)Click Compute.The projected area reported is0.0003m2.(d)Close the Projected Surface Areas panel.3.Change the Solution Zones for user-de?ned scalars2and3to all zones.De?ne?→User-De?ned?→Scalars...This is an optional step which allows UDS-2and UDS-3to be postprocessed on bothuid and solid zones.4.Con?gure the PEM model.De?ne?→Models?→PEMFC...(a)Click the Anode tab.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel Cell i.Select Current Collector in the Anode Zone Type group box and select current-a from the Zone(s)selection list.ii.Select Flow Channel in the Anode Zone Type group box and select channel-afrom the Zone(s)selection list.iii.Select Di?usion Layer in the Anode Zone Type group box and select gdl-a fromthe Zone(s)selection list.iv.Select Catalyst Layer in the Anode Zone Type group box and select catalyst-afrom the Zone(s)selection list.(b)Click the Membrane tab.i.Select membrane from the Membrane Zone(s)selection list.(c)Click the Cathode tab.i.Select Current Collector in the Cathode Zone Type group box and selectcurrent-c from the Zone(s)selection list.ii.Select Flow Channel in the Cathode Zone Type and select channel-c from theZone(s)selection list.iii.Select Di?usion Layer in the Cathode Zone Type and select gdl-c from theZone(s)selection list.iv.Select Catalyst Layer in the Cathode Zone Type and select catalyst-c from theZone(s)selection list.(d)Click the Reports tab.i.Specify the value of Membrane-Electrode-Assembly Projected Area to0.0003m2. Recall that this value was obtained earlier in the tutorial.ii.Select wall-terminal-a from the Anode selection list and select wall-terminal-cfrom the Cathode selection list.(e)Click OK to close the PEM panel.FLUENT reports in the console that the energy equation has been enabled auto-matically.For convenience,the PEMFC model also automatically enables species transport and creates default materials.Step3:MaterialsRetain the default settings for the materials.Step4:Operating ConditionsDe?ne?→Operating Conditions...1.Set the Operating Pressure to200000pascal.2.Click OK to close the Operating Conditions panel.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel Cell Step5:Boundary ConditionsDe?ne?→Boundary Conditions...There are several zones which must be speci?ed in the boundary conditions panel.These are the anode and cathode voltage terminals,as well as the inlets and outlets.1.Set boundary conditions for the anode voltage terminal,wall-terminal-a.At this surface,the voltage is grounded and the temperature is constant.(a)Click the Thermal tab and enter353K for Temperature.(b)Click the UDS tab.i.Select Speci?ed Value from the Electric Potential drop-down list User-De?nedScalar Boundary Condition group box.ii.Enter0for Electric Potential in the User-De?ned Scalar Boundary Value groupbox.This boundary condition represents a grounded terminal.(c)Click OK to close the Wall panel.2.Set boundary conditions for the cathode voltage terminal,wall-terminal-c.At this surface,the voltage is maintained at a constant,positive value.(a)Click the Thermal tab and enter353K for Temperature.(b)Click the UDS tab.i.Select Speci?ed Value from the Electric Potential drop-down list User-De?nedScalar Boundary Condition group box.ii.Enter0.75for Electric Potential in the User-De?ned Scalar Boundary Valuegroup box.This boundary condition represents a terminal operating at0.75Volts.iii.Click OK to close the Wall panel.To calculate an IV polarization curve,you should vary the Electric Potential forthe cathode,starting from a voltage near the open circuit voltage and gradually decreasing it,converging the solution each time you change the value.3.Set boundary conditions for the anode gas?ow inlet,inlet-a.At this inlet,a humidi?ed hydrogen stream enters the fuel cell.No liquid enters the channel.(a)Enter6.0e-7kg/s for Mass Flow Rate and0for Supersonic/Initial Gauge Pressure.(b)Click the Thermal tab and enter353K for Temperature.(c)Click the Species tab and set the mass fractions of h2,o2and h2o to0.8,0.0,and0.2,respectively.(d)Click the UDS tab and select Speci?ed Value from the Water Saturation drop-down list in the User-De?ned Scalar Boundary Condition group box.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel Cell(e)Enter0for Water Saturation in the User-De?ned Scalar Boundary Value group box.(f)Click OK to close the Mass-Flow Inlet panel.4.Set boundary conditions for the cathode gas?ow inlet,inlet-c.At this inlet,a humidi?ed air stream enters the fuel cell.No liquid enters the channel.(a)Enter5.0e-6kg/s for Mass Flow Rate.(b)Click the Thermal tab and enter353K for Total Temperature.(c)Click the Species tab and set the mass fractions of h2,o2and h2o to0.0,0.2,and0.1,respectively.(d)Click the UDS tab and select Speci?ed Value from the Water Saturation drop-down list in the User De?ned Scalar Boundary Condition group box.(e)Enter0for Water Saturation in the User-De?ned Scalar Boundary Value group box.(f)Click OK to close the Mass-Flow Inlet panel.5.Set boundary conditions for the anode gas?ow outlet,outlet-a.(a)Click the Thermal tab and enter353K for Back?ow Total Temperature.(b)Click OK to close the Pressure Outlet panel.6.To set boundary conditions for the cathode gas?ow outlet,copy the boundary con-ditions from outlet-a to outlet-c7.Close the Boundary Conditions panel.Step6:Solution ControlsThe default solver settings are not su?cient to obtain a converged solution.Therefore,the following modi?cations must be made.1.Set the under-relaxation factor for Pressure to0.7,Momentum to0.3,Protonic Po-tential to0.95,and Water Content to0.95.Solve?→Controls?→Solution...2.Modify the multigrid settings.Solve?→Controls?→Multigrid...(a)Select F-Cycle from the Cycle Type drop-down lists for all equations.You will need to scroll down to set all equations.(b)Enter0.001for Termination Restriction for h2,o2,h2o,and Water Saturation.(c)Select BCGSTAB from the Stabilization Method drop-down list for h2,o2,h2o,Water Saturation,Electric Potential and Protonic Potential.(d)Enter0.0001for Termination Restriction for Electric Potential and Protonic Po-tential.(e)Increase the value of Max Cycles to50in the Algebraic Multigrid Controls group box.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel Cell(f)Click OK to close the Multigrid Controls panel.3.Enable the plotting of residuals.Solve?→Monitors?→Residual...Note:The PEMFC model automatically disables convergence checking for all equa-tions.4.Initialize the solution.Solve?→Initialize?→Initialize...(a)Set Temperature to353K.(b)Click Apply.(c)Click Init and close the Solution Initialization panel.5.Save the case and data?les as pem-single-channel.cas.gz andpem-single-channel.dat.gz.File?→Write?→Case&Data...!If you want to run the calculations in parallel,exit FLUENT and start a parallel session at this point.Open the case and data? les you saved in the previousstep and proceed.6.Request200iterations.The solution residuals will drop to acceptable values.Solve?→Iterate...The solution residual plot should resemble that shown in Figure4.The average currentdensity is displayed in the console at the end of each iteration.At the end of the calculations,the current density is reported as approximately0.324A/cm2.Step7:Postprocessing1.Create surfaces for postprocessing.Surface?→Iso-Surface...(a)Select Grid...and Z-Coordinate from the Surface of Constant drop-down lists.(b)Click Compute.(c)Enter62.5for Iso-Values(mm).(d)Enter plane-xy for New Surface Name.(e)Click Create.(f)Similarly,create another surface along the length of the fuel cell.This surfaceshould be a surface of constant X-Coordinate,with a value /doc/d88704d349649b6648d74700.html this surface plane-yz.(g)Close the Iso-Surface panel.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel CellFigure4:Residual Plot2.Create custom vectors for display.Display?→Vectors...(a)Click the Custom Vectors...button to open the Custom Vectors panel.i.Enter current-?ux-density for Vector Name.ii.Select User De?ned Memory...and X Current Flux Density from the X Com-ponent drop-down lists.iii.Select User De?ned Memory...and Y Current Flux Density from the Y Com-ponent drop-down lists.iv.Select User De?ned Memory...and Z Current Flux Density from the Z Com-ponent drop-down lists.v.Click De?ne and close the Custom Vectors panel.(b)Select current-?ux-density from the Vectors of drop-down list.(c)Select?lled-arrow from the Style drop-down list.(d)Click the Vector Options...button to open the Vector Options panel.i.Enter0.5for Scale Head.ii.Click Apply and close the Vector Options panel.(e)Select User-De?ned Memory and Current Flux Density Magnitude from the Colorby drop-down lists.(f)Enable Draw Grid from the Options list to open the Grid Display panel.i.Deselect all surfaces except plane-xy from the Surfaces selection list.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel CellFigure5:Current Flux at a Cross-Section Midway Down the Length of the PEM Channel ii.Ensure that Edges is enabled from the Options list and Feature is selectedfrom the Edge Type list.iii.Click Display and close the Grid Display panel.(g)Select plane-xy in the Surfaces selection list.(h)Click Display and close the Vectors panel.3.Auto-?t the display to the graphics window by pressing Ctrl-A on the keyboard./doc/d88704d349649b6648d74700.html pare your results with those shown in Figure5.Note:The maximum current density occurs in the regions between the channels and also that the plot units are A/m2.5.Plot contours of hydrogen mass fraction along the surface plane-yzDisplay?→Contours...(a)Enable Filled and Draw Grid from the Options list.(b)Select plane-yz from the Surfaces selection list in the Grid Display panel.(c)Click Display and close the Grid Display panel.(d)Restore the right view.Display?→View...i.Select right from the Views selection list.ii.Click Apply and close the Views panel.The display updates to the right-hand view.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel CellFigure6:Contours of Hydrogen Mass Fraction Along the Channel Length(e)Select plane-yz from the Surfaces selection list.(f)Select Species...and Mass fraction of h2from the Contours of drop-down lists.(g)Click Display.The resulting display is di?cult to visualize since the aspect ratio of the channelis large.You can change the way FLUENT displays data using the following steps:Display?→Scene...i.Select all entries from the Names selection list.ii.Click the Transform...button to open the Transformations panel.iii.Set Z to0.1in the Scale group box.iv.Click Apply and close the Transformations panel.The graphics display will be scaled accordingly.(h)Auto-?t the image to the window by pressing Ctrl-A in the graphics window.(i)Compare your result with that shown in Figure6.The?ow in the anode(upper)channel is from left to right.Note that the hydrogenmass fraction decreases in the direction of?ow.This is due to water being pulled through the membrane along with hydrogen as it is consumed in the fuel cell.6.Plot contours of oxygen mass fraction along the surface plane-yzDisplay?→Contours...(a)Follow the same procedure to generate the oxygen mass fraction contour plotthat you used to generate the hydrogen mass fraction contours.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel CellFigure7:Contours of Oxygen Mass Fraction Along the Channel Length(b)Compare your result with that shown in Figure7.The?ow in the cathode(lower) channel is from right to left.As expected,the oxygen mass fraction decreases inthe direction of?ow.7.Verify that global conservation of mass is observed.This will be done using a few text user interface(TUI)commands and basic electrochemistry concepts.(a)Compute the net oxygen consumption.Enter the TUI command:/report/species-mass-flowIf you have read a data?le instead of performing iterations,you must performat least one iteration in order to populate this data from the solver.The output is as follows:zone22(inlet-a):(4.8e-0701.2e-07)zone21(inlet-c):(01e-065e-07)zone20(outlet-a):(-4.6958153e-07-8.9805489e-11-2.8157297e-07)zone19(outlet-c):(-2.6251551e-10-9.1930928e-07-4.2904232e-07)zone53(wall-ch-a-shadow):(000)zone15(wall-ch-c):(000)zone29(wall-ends:029):(000)zone30(wall-ends:030):(000)zone31(wall-ends:031):(000)zone32(wall-ends:032):(000)zone33(wall-ends:033):(000)Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel Cell zone51(wall-gdl-a-shadow):(000) zone11(wall-gdl-c):(000)zone1(wall-sides:001):(000)zone23(wall-sides:023):(000)zone25(wall-sides:025):(000)zone26(wall-sides:026):(000)zone27(wall-sides:027):(000)net species-mass-flow:(1.0155954e-088.0600915e-08-9.061529e-08)To interpret the output,each line can be read aszone num(zone-name):(˙m1˙m2...˙m n)where the subscripts1,2,...,nrefer to each species being calculated.Here,we are considering three species,namely h2,o2and h2o.Therefore,the second value in each list is the calculatedoxygen mass?ow rate in kilograms per second.In addition,a negative numberindicates?ow out of the domain from that boundary.From the last line,the netoxygen consumed is8.06×10?8kg/s.The molecular weight of oxygen is31.9988kg/kmol.Also,since the valence of a diatomic oxygen molecule is4,there are4kmol of electrons released per kmol of oxygen.Finally Faraday’s constant is9.6485×107C/kmol-electrons.Thus,the total release of electrons(which isequivalent to the current in Amperes),isI=˙mv FM=(8.06×10?8)(4.0)(9.6485×107)31.9988=0.972A(1)The total current is obtained by integrating the current density over the surface of the terminal.This integral value can be calculated in several ways.One way is to multiply the membrane area by the reported current density.This givesI=(0.0003)(0.3241)(100)2=0.972A(2) Alternatively,you can integrate the user memory Y Current Flux Density(the y-component of current density)over the terminal surface.This integration yields an accurate result since the y-direction is normal to the terminal.To do this, use the Surface Integrals panel.Report?→Surface Integrals...i.Select Report Type from the Integral drop-down list.ii.Select User De?ned Memory...and Y Current Flux Density from the Field Variable drop-down lists.iii.Select wall-terminal-a from the Surfaces selection list and click Compute.The absolute value of the number reported is approximately0.972A.We have electrochemical balance in the calculations. iv.Close the Surface Integrals panel.Modeling a Single-Channel,Counter-Flow Polymer Electrolyte Membrane(PEM)Fuel CellSummaryIn this tutorial,you learned how to set up and model a single-channel PEM fuel cell.Themodel provides detailed information on the distribution of current and voltage on all theelectrically conducting regions,along with species and current?ux density distributionthroughout the fuel cell.。
1.GAMBIT建模(1)操作1;双击桌面“GAMBIT快捷方式”进入操作空间,如下图图1几何模型创建遵循原则:按组成几何模型的几何拓扑结构,由低层向高层创建,即按由点到面,再到体的顺序创建(这样减的好处是便于网格的划分和Fluent求解时边界条件的设定,缺点是步骤繁琐)。
对于较简单的模型,可省略低拓扑结构,直接建立最终模型。
(2)创建节点操作2:鼠标左键依次单击Operation Geometry进入下图图2操作3:左键单击上图Apply,创建第一个点(坐标(0、0、0)),如下图,可以发现坐标原点显示白色图3操作4:在上图2中Global下输入x:200,y:0,z: 0,左键单击Apply, 如下图,操作5:左键单击(作用:窗口显示),如下图重复操作4和操作5,以此建立点(0,-0.0025,0)、(0.05,-0.0025,0)、(0.05、-0.0125)(0.15,0.0025,0),(0.15,0.0125,0),(0.2,0.0025),如下图,(3)节点成线操作1:左键依次单击,然后Shift+鼠标左键依次(为沿围成图形各点顺序)单击所创建的点,如下图,左键单击Apply,如下图,(4)连线成面操作1:左键单击,Shift+左键依次单击图中各线段,如下图,左键单击Apply,如下图,(5)网格划分网格划分遵循原则与模型创建类似操作1:左键依次单击,如下图,操作2:Shift+左键选中模型中一条边操作3:左键单击,弹出菜单中选择interval count,在左侧输入框中输入;200(为该条边上网格节点数),如下图重复操作2和操作3,将沿y方向的短边和长边节点数分别设定为20和80.注意:对于相互平行两边,节点设置可只在一边进行;如果两边均设定,切记两边节点数要一致。
操作4:左键单击,Shift+左键单击图中任一边,所有边显示红色,然后左键单击Apply,如下图,(6)设置边界类型操作1: 左键单击,如下图,操作2:Shift+左键单击模型最左侧沿y方向的短边(显红色),单击,弹出菜单中选择PRESSURE-INLET,左键单击Apply,如下图,操作3:Shift+左键单击模型最右侧沿y方向的短边(显红色),单击,弹出菜单中选择PRESSURE-OUTLET,左键单击Apply,如下图,操作4:其它边不设置,默认为壁面条件。
fluent流固耦合传热算例摘要:I.引言- 介绍fluent 软件和流固耦合传热算例II.流固耦合传热的基本概念- 解释流固耦合传热- 说明流固耦合传热在工程领域的重要性III.fluent 软件介绍- 介绍fluent 软件的背景和功能- 说明fluent 软件在流固耦合传热计算方面的应用IV.流固耦合传热算例- 介绍一个具体的流固耦合传热算例- 详细描述算例的步骤和结果V.结论- 总结流固耦合传热算例的重要性- 提出进一步研究的建议正文:I.引言fluent 软件是一款专业的流体动力学模拟软件,广泛应用于航空航天、汽车制造、能源等行业。
在fluent 中,流固耦合传热是一个重要的计算功能。
本文将介绍fluent 软件和流固耦合传热算例,并通过一个具体的算例详细说明流固耦合传热在工程领域中的应用。
II.流固耦合传热的基本概念流固耦合传热是指在流体流动过程中,由于流体和固体壁面之间的温度差而产生的热传递现象。
在实际工程中,流体和固体之间的热传递过程往往是非常复杂的,需要通过数值模拟来进行分析。
fluent 软件提供了一种流固耦合传热计算的功能,可以帮助工程师更好地理解和优化工程过程中的热传递现象。
III.fluent 软件介绍fluent 软件由美国ANSYS 公司开发,是一款功能强大的流体动力学模拟软件。
fluent 软件可以模拟多种流体流动和传热现象,包括稳态和瞬态模拟、层流和紊流模拟、等温、绝热和热传导模拟等。
在fluent 中,用户可以自定义模型和求解器,以满足不同工程需求。
在流固耦合传热方面,fluent 软件提供了一种耦合求解器,可以将流体流动和固体传热两个问题同时求解。
这种耦合求解器可以大大提高计算效率,并更好地模拟实际工程中的热传递过程。
IV.流固耦合传热算例下面我们通过一个具体的算例来说明fluent 软件在流固耦合传热计算方面的应用。
算例描述:一个矩形通道中,流体流动与固体壁面的热传递过程。
经过痛苦的一段经历,终究将局部成绩本相大白,为了使保位同仁不再经过我之痛苦,此刻将本人多孔介质经验公布如下,但愿各位能加精:之杨若古兰创作1.Gambit中划分网格以后,定义须要做为多孔介质的区域为fluid,与缺省的fluid分别开来,再定义其名称,我习气将名称定义为porous;2.在fluent中定义鸿沟条件define-boundary condition-porous(刚定义的名称),将其设置鸿沟条件为fluid,点击set 按钮即弹出与fluid鸿沟条件一样的对话框,选中porous zone与laminar复选框,再点击porous zone标签即出现一个带有滚动条的界面;3.porous zone设置方法:1)定义矢量:二维定义一个矢量,第二个矢量方向不必定义,是与第一个矢量方向正交的;三维定义二个矢量,第三个矢量方向不必定义,是与第一、二个矢量方向正交的;(如何晓得矢量的方向:打开grid图,看看X,Y,Z的方向,如果是X向,矢量为1,0,0,同理Y向为0,1,0,Z向为0,0,1,如果所须要的方向与坐标轴正向相反,则定义矢量为负)圆锥坐标与球坐标请参考fluent帮忙.2)定义粘性阻力1/a与内部阻力C2:请参看本人上一篇博文“终究搞清fluent中多孔粘性阻力与内部阻力的计算方法”,此处不赘述;3)如果了定义粘性阻力1/a与内部阻力C2,就不必定义C1与C0,由于这是两种分歧的定义方法,C1与C0只在幂率模型中出现,该处坚持默认就行了;4)定义孔隙率porousity,默认值1暗示全开放,此值按实验测值填写即可.完了,其他设置与普通k-e或RSM不异.总结一下,与君共享!Tutorial 7. Modeling Flow Through Porous Media IntroductionMany industrial applications involve the modeling of flow through porous media, suchas filters, catalyst beds, and packing. This tutorial illustrates how to set up and solve aproblem involving gas flow through porous media.The industrial problem solved here involves gas flow through acatalytic converter. Catalyticconverters are commonly used to purify emissions from gasoline and diesel enginesby converting environmentally hazardous exhaust emissions to acceptable substances.Examples of such emissions include carbon monoxide (CO), nitrogen oxides (NOx), andunburned hydrocarbon fuels. These exhaust gas emissions are forced through a substrate,which is a ceramic structure coated with a metal catalyst such as platinum or palladium.The nature of the exhaust gas flow is a very important factor in determining the performanceof the catalytic converter. Of particular importance is the pressure gradientand velocity distribution through the substrate. Hence CFD analysis is used to designefficient catalytic converters: by modeling the exhaust gas flow, the pressure drop andthe uniformity of flow through the substrate can be determined. In this tutorial, FLUENTis used to model the flow of nitrogen gas through a catalytic converter geometry, so thatthe flow field structure may be analyzed.This tutorial demonstrates how to do the following:_ Set up a porous zone for the substrate with appropriate resistances._ Calculate a solution for gas flow through the catalytic converterusing the pressurebasedsolver._ Plot pressure and velocity distribution on specified planes of the geometry._ Determine the pressure drop through the substrate and the degree of non-uniformityof flow through cross sections of the geometry using X-Y plots and numerical reports.Problem DescriptionThe catalytic converter modeled here is shown in Figure 7.1. The nitrogen flows inthrough the inlet with a uniform velocity of 22.6 m/s, passes through a ceramic monolithsubstrate with square shaped channels, and then exits through the outlet.While the flow in the inlet and outlet sections is turbulent, the flow through the substrateis laminar and is characterized by inertial and viscous loss coefficients in the flow (X)direction. The substrate is impermeable in other directions, which is modeled using losscoefficients whose values are three orders of magnitude higher than in the X direction.Setup and SolutionStep 1: Grid1. Read the mesh file (catalytic converter.msh).File /Read /Case...2. Check the grid.Grid /CheckFLUENT will perform various checks on the mesh and report the progress in theconsole. Make sure that the minimum volume reported is a positive number.3. Scale the grid.Grid!Scale...(a) Select mm from the Grid Was Created In drop-down list.(b) Click the Change Length Units button.All dimensions will now be shown in millimeters.(c) Click Scale and close the Scale Grid panel.4. Display the mesh.Display /Grid...(a) Make sure that inlet, outlet, substrate-wall, and wall are selected in the Surfacesselection list.(b) Click Display.(c) Rotate the view and zoom in to get the display shown in Figure7.2.(d) Close the Grid Display panel.The hex mesh on the geometry contains a total of 34,580 cells.Step 2: Modelsfine /Models /Solver...2. Select the standard k-εfine/ Models /Viscous...Step 3: Materials1. Add nitrogen to the list of flfine /Materials...(a) Click the Fluent Database... button to open the Fluent Database Materialspanel.i. Select nitrogen (n2) from the list of Fluent Fluid Materials.ii. Click Copy to copy the information for nitrogen to your list of fluid materials.iii. Close the Fluent Database Materials panel.(b) Close the Materials panel.Step 4: Boundary Conditions.Define /Boundary Conditions...1. Set the boundary conditions for the fluid (fluid).(a) Select nitrogen from the Material Name drop-down list.(b) Click OK to close the Fluid panel.2. Set the boundary conditions for the substrate (substrate).(a) Select nitrogen from the Material Name drop-down list.(b) Enable the Porous Zone option to activate the porous zone model.(c) Enable the Laminar Zone option to solve the flow in the porous zone withoutturbulence.(d) Click the Porous Zone tab.i. Make sure that the principal direction vectors are set as shown in e the scroll bar to access the fields that are not initially visible in thepanel.ii. Enter the values in Table 7.2 for the Viscous Resistance and Inertial Resistance.Scroll down to access the fields that are not initially visible in the panel.(e) Click OK to close the Fluid panel.3. Set the velocity and turbulence boundary conditions at the inlet (inlet).(a) Enter 22.6 m/s for the Velocity Magnitude.(b) Select Intensity and Hydraulic Diameter from the Specification Method dropdownlist in the Turbulence group box.(c) Retain the default value of 10% for the Turbulent Intensity.(d) Enter 42 mm for the Hydraulic Diameter.(e) Click OK to close the Velocity Inlet panel.4. Set the boundary conditions at the outlet (outlet).(a) Retain the default setting of 0 for Gauge Pressure.(b) Select Intensity and Hydraulic Diameter from the Specification Method dropdownlist in the Turbulence group box.(c) Enter 5% for the Backflow Turbulent Intensity.(d) Enter 42 mm for the Backflow Hydraulic Diameter.(e) Click OK to close the Pressure Outlet panel.5. Retain the default boundary conditions for the walls (substrate-wall and wall) andclose the Boundary Conditions panel.Step 5: Solution1. Set the solution parameters.Solve /Controls /Solution...(a) Retain the default settings for Under-Relaxation Factors.(b) Select Second Order Upwind from the Momentum drop-down list in the Discretizationgroup box.(c) Click OK to close the Solution Controls panel./Monitors /Residual...(a) Enable Plot in the Options group box.(b) Click OK to close the Residual Monitors panel.3. Enable the plotting of the mass flow rate at the outlet.Solve / Monitors /Surface...(a) Set the Surface Monitors to 1.(b) Enable the Plot and Write options for monitor-1, and click the Define... buttonto open the Define Surface Monitor panel.i. Select Mass Flow Rate from the Report Type drop-down list. ii. Select outlet from the Surfaces selection list.iii. Click OK to close the Define Surface Monitors panel.(c) Click OK to close the Surface Monitors panel.4. Initialize the solution from the inlet.Solve /Initialize /Initialize...(a) Select inlet from the Compute From drop-down list.(b) Click Init and close the Solution Initialization panel.5. Save the case file (catalytic converter.cas).File /Write /Case...6. Run the calculation by requesting 100 iterations.Solve /Iterate...(a) Enter 100 for the Number of Iterations.(b) Click Iterate.The FLUENT calculation will converge in approximately 70 iterations. By thispoint the mass flow rate monitor has attended out, as seen in Figure 7.3.(c) Close the Iterate panel.7. Save the case and data files (catalytic converter.cas and catalytic converter.dat).File /Write /Case & Data...Note: If you choose a file name that already exists in the current folder, FLUENTwill prompt you for confirmation to overwrite the file.Step 6: Post-processing1. Create a surface passing through the centerline for post-processing purposes.Surface/Iso-Surface...(a) Select Grid... and Y-Coordinate from the Surface of Constant drop-down lists.(b) Click Compute to calculate the Min and Max values.(c) Retain the default value of 0 for the Iso-Values.(d) Enter y=0 for the New Surface Name.(e) Click Create.2. Create cross-sectional surfaces at locations on either side of the substrate, as wellas at its center.Surface /Iso-Surface...(a) Select Grid... and X-Coordinate from the Surface of Constant drop-down lists.(b) Click Compute to calculate the Min and Max values.(c) Enter 95 for Iso-Values.(d) Enter x=95 for the New Surface Name.(e) Click Create.(f) In a similar manner, create surfaces named x=130 and x=165 with Iso-Valuesof 130 and 165, respectively. Close the Iso-Surface panel after all the surfaceshave been created.3. Create a line surface for the centerline of the porous media. Surface /Line/Rake...(a) Enter the coordinates of the line under End Points, using the starting coordinateof (95, 0, 0) and an ending coordinate of (165, 0,0), as shown.(b) Enter porous-cl for the New Surface Name.(c) Click Create to create the surface.(d) Close the Line/Rake Surface panel.4. Display the two wall zones (substrate-wall and wall).Display/Grid...(a) Disable the Edges option.(b) Enable the Faces option.(c) Deselect inlet and outlet in the list under Surfaces, and make sure that onlysubstrate-wall and wall are selected.(d) Click Display and close the Grid Display panel.(e) Rotate the view and zoom so that the display is similar to Figure 7.2.5. Set the lighting for the display.Display /Options...(a) Enable the Lights On option in the Lighting Attributes group box.(b) Retain the default selection of Gourand in the Lighting drop-down list.(c) Click Apply and close the Display Options panel.6. Set the transparency parameter for the wall zones (substrate-wall and wall).Display/Scene...(a) Select substrate-wall and wall in the Names selection list.(b) Click the Display... button under Geometry Attributes to openthe DisplayProperties panel.i. Set the Transparency slider to 70.ii. Click Apply and close the Display Properties panel.(c) Click Apply and then close the Scene Description panel.7. Display velocity vectors on the y=0 surface.Display /Vectors...(a) Enable the Draw Grid option.The Grid Display panel will open.i. Make sure that substrate-wall and wall are selected in the list under Surfaces.ii. Click Display and close the Display Grid panel.(b) Enter 5 for the Scale.(c) Set Skip to 1.(d) Select y=0 from the Surfaces selection list.(e) Click Display and close the Vectors panel.The flow pattern shows that the flow enters the catalytic converter as a jet, withrecirculation on either side of the jet. As it passes through the porous substrate, itdecelerates and straightens out, and exhibits a more uniform velocity distribution.This allows the metal catalyst present in the substrate to be more effective.Figure 7.4: Velocity Vectors on the y=0 Plane8. Display filled contours of static pressure on the y=0 plane. Display /Contours...(a) Enable the Filled option.(b) Enable the Draw Grid option to open the Display Grid panel.i. Make sure that substrate-wall and wall are selected in the list under Surfaces.ii. Click Display and close the Display Grid panel.(c) Make sure that Pressure... and Static Pressure are selected from the Contoursof drop-down lists.(d) Select y=0 from the Surfaces selection list.(e) Click Display and close the Contours panel.Figure 7.5: Contours of the Static Pressure on the y=0 planeThe pressure changes rapidly in the middle section, where the fluid velocity changesas it passes through the porous substrate. The pressure drop can be high, due to theinertial and viscous resistance of the porous media. Determining this pressure dropis a goal of CFD analysis. In the next step, you will learn how to plot the pressuredrop along the centerline of the substrate.9. Plot the static pressure across the line surface porous-cl.Plot /XY Plot...(a) Make sure that the Pressure... and Static Pressure are selectedfrom the Y AxisFunction drop-down lists.(b) Select porous-cl from the Surfaces selection list.(c) Click Plot and close the Solution XY Plot panel.Figure 7.6: Plot of the Static Pressure on the porous-cl Line SurfaceIn Figure 7.6, the pressure drop across the porous substrate can be seen to beroughly 300 Pa.10. Display filled contours of the velocity in the X direction on the x=95, x=130 andx=165 surfaces.Display /Contours...(a) Disable the Global Range option.(b) Select Velocity... and X Velocity from the Contours of drop-down lists.(c) Select x=130, x=165, and x=95 from the Surfaces selection list, and deselecty=0.(d) Click Display and close the Contours panel.The velocity profile becomes more uniform as the fluid passes through the porousmedia. The velocity is very high at the center (the area in red) just before thenitrogen enters the substrate and then decreases as it passes through and exits thesubstrate. The area in green, which corresponds to a moderate velocity, increasesin extent.Figure 7.7: Contours of the X Velocity on the x=95,x=130, and x=165 Surfaces11. Use numerical reports to determine the average, minimum, and maximum of thevelocity distribution before and after the porous substrate.Report /Surface Integrals...(a) Select Mass-Weighted Average from the Report Type drop-down list.(b) Select Velocity and X Velocity from the Field Variable drop-down lists.(c) Select x=165 and x=95 from the Surfaces selection list.(d) Click Compute.(e) Select Facet Minimum from the Report Type drop-down list and click Computeagain.(f) Select Facet Maximum from the Report Type drop-down list and click Computeagain.(g) Close the Surface Integrals panel.The numerical report of average, maximum and minimum velocity can be seen inthe main FLUENT console, as shown in the following example:The spread between the average, maximum, and minimum values for X velocitygives the degree to which the velocity distribution isnon-uniform. You can also usethese numbers to calculate the velocity ratio (i.e., the maximum velocity divided bythe mean velocity) and the space velocity (i.e., the product of the mean velocity andthe substrate length).Custom field functions and UDFs can be also used to calculate more complex measuresof non-uniformity, such as the standard deviation and the gamma uniformityindex.SummaryIn this tutorial, you learned how to set up and solve a problem involving gas flow throughporous media in FLUENT. You also learned how to perform appropriate post-processingto investigate the flow field, determine the pressure drop across the porous media andnon-uniformity of the velocity distribution as the fluid goes through the porous media.Further ImprovementsThis tutorial guides you through the steps to reach an initial solution. You may be ableto obtain a more accurate solution by using an appropriate higher-order discretizationscheme and by adapting the grid. Grid adaption can also ensure that the solution isindependent of the grid. These steps are demonstrated in Tutorial 1.。
fluent算例集Fluent算例集是一个强大的工具,可以帮助人们更好地理解和运用流畅的英语。
在这个算例集中,我们将探讨一些常见的英语用法和表达,以及如何在日常交流中更加流利地运用它们。
第一部分:常用口语表达1. 意思是什么?在日常对话中,我们经常会遇到一些不熟悉的词汇或表达。
这时,我们可以用这个句子来询问对方的意思。
比如,如果对方说了一个生词,我们可以问:“这个词的意思是什么?”2. 你是什么意思?有时候,在交流中我们可能会遇到一些误解或不清楚对方的意图。
这个句子可以用来询问对方的具体意思,以便更好地理解对方的观点。
3. 我没听懂。
当我们在听别人说话时,有时候会遇到一些听不懂的地方。
这个句子可以用来表达自己的困惑,并希望对方再解释一遍。
4. 你能再说一遍吗?类似于上一个句子,这个句子也可以用来请求对方再重复一遍刚才的话。
常用在电话或者嘈杂的环境中。
5. 这是什么意思?当我们看到一些不熟悉的符号、图标或者标志时,可以用这个句子来询问它们的具体含义。
比如,在一家陌生的餐厅,我们可以指着菜单上的某个词汇问服务员:“这是什么意思?”第二部分:语法用法示例1. 一般现在时一般现在时是英语中最基础的时态之一。
它用来表达经常性的动作、习惯、真理、现在的状态等。
比如,我们可以说:“我每天早上六点起床。
”或者“太阳从东方升起。
”2. 现在进行时3. 过去时过去时用来表达已经发生的动作或者状态。
比如,我们可以说:“昨天我去了一趟图书馆。
”或者“他在上个月搬家了。
”4. 将来时将来时用来表达将来将要发生的动作或者状态。
比如,我们可以说:“明天我会去参加一个会议。
”或者“他下个月要去巴黎旅行。
”第三部分:写作技巧1. 引入段在写作中,引入段是非常重要的一部分。
它用来引入文章的主题,并吸引读者的注意力。
可以使用一些有趣的事实、引用或者问题来引起读者的兴趣。
2. 主体段落主体段落是文章的核心部分,用来展开和论证文章的主题。
前言为了使学生尽快熟悉计算流体软件FLUENT以及更好的掌握计算流体力学的计算模型,本书编制了几个简单的模型,包括了组分燃烧、管内流动、换热和房间温度场四个方面的内容。
其中概括了二维和三维的模型,描述详细,可根据步骤建模、划分网格和计算以及后处理。
本书不可能面面具到并进行详细讲解,但相信读者通过本书的学习,一定能领会其中的技巧。
目录前言﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍1 燃烧器内甲烷和空气的燃烧﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍3 管内层流流动数值计算﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍ 38 蒸汽喷射器内的传热模拟﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍ 52 组分传输与气体燃烧算例﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍ 75 空调房间温度场的模拟﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍﹍102燃烧器内甲烷和空气的燃烧问题描述这个问题在图1中以图解的形式表示出来。
此几何体包括一个简化的向燃烧腔加料的燃料喷嘴,由于几何结构对称可以仅做出燃烧室几何体的1/4模型。
喷嘴包括两个同心管,其直径分别是4个单位和10个单位,燃烧室的边缘与喷嘴下的壁面融合在一起。
图1:问题图示一、利用GAMBIT建立计算模型启动GAMBIT。
第一步:选择一个解算器选择用于进行CFD计算的求解器。
操作:Solver -> FLUENT5/6第二步:生成两个圆柱体1、生成一个柱体以形成燃烧室操作:GEOMETRY-> VOLUME-> CREATE VOLUMER打开Create Real Cylinder 窗口,如图2所示a) 在柱体的Height 中键入值1.2。
b) 在柱体的Radius 1 中键入值0.4。
Radius 2的文本键入框可留为空白,GAMBIT 将默认设定为Radius 1值相等。
c) 选择Positive Z (默认)作为Axis Location 。
d) 点击Apply 按钮。
2、按照上述步骤以生成一个Height =2,Radius 1 =1并以positive z 为轴的柱体。
AnsysFluent基础详细⼊门教程(附简单算例)Ansys Fluent基础详细⼊门教程(附简单算例)当你决定使FLUENT解决某⼀问题时,⾸先要考虑如下⼏点问题:定义模型⽬标:从CFD模型中需要得到什么样的结果?从模型中需要得到什么样的精度;选择计算模型:你将如何隔绝所需要模拟的物理系统,计算区域的起点和终点是什么?在模型的边界处使⽤什么样的边界条件?⼆维问题还是三维问题?什么样的⽹格拓扑结构适合解决问题?物理模型的选取:⽆粘,层流还湍流?定常还是⾮定常?可压流还是不可压流?是否需要应⽤其它的物理模型?确定解的程序:问题可否简化?是否使⽤缺省的解的格式与参数值?采⽤哪种解格式可以加速收敛?使⽤多重⽹格计算机的内存是否够⽤?得到收敛解需要多久的时间?在使⽤CFD分析之前详细考虑这些问题,对你的模拟来说是很有意义的。
第01章fluent介绍及简单算例 (2)第02章fluent⽤户界⾯22 (3)第03章fluent⽂件的读写 (5)第04章fluent单位系统 (8)第05章fluent⽹格 (10)第06章fluent边界条件 (36)第07章fluent流体物性 (55)第08章fluent基本物理模型 (63)第11章传热模型 (75)第22章fluent 解算器的使⽤ (82)第01章fluent介绍及简单算例FLUENT是⽤于模拟具有复杂外形的流体流动以及热传导的计算机程序。
对于⼤梯度区域,如⾃由剪切层和边界层,为了⾮常准确的预测流动,⾃适应⽹格是⾮常有⽤的。
FLUENT解算器有如下模拟能⼒:●⽤⾮结构⾃适应⽹格模拟2D或者3D流场,它所使⽤的⾮结构⽹格主要有三⾓形/五边形、四边形/五边形,或者混合⽹格,其中混合⽹格有棱柱形和⾦字塔形。
(⼀致⽹格和悬挂节点⽹格都可以)●不可压或可压流动●定常状态或者过渡分析●⽆粘,层流和湍流●⽜顿流或者⾮⽜顿流●对流热传导,包括⾃然对流和强迫对流●耦合热传导和对流●辐射热传导模型●惯性(静⽌)坐标系⾮惯性(旋转)坐标系模型●多重运动参考框架,包括滑动⽹格界⾯和rotor/stator interaction modeling的混合界⾯●化学组分混合和反应,包括燃烧⼦模型和表⾯沉积反应模型●热,质量,动量,湍流和化学组分的控制体源●粒⼦,液滴和⽓泡的离散相的拉格朗⽇轨迹的计算,包括了和连续相的耦合●多孔流动●⼀维风扇/热交换模型●两相流,包括⽓⽳现象●复杂外形的⾃由表⾯流动上述各功能使得FLUENT具有⼴泛的应⽤,主要有以下⼏个⽅⾯●Process and process equipment applications●油/⽓能量的产⽣和环境应⽤●航天和涡轮机械的应⽤●汽车⼯业的应⽤●热交换应⽤●电⼦/HV AC/应⽤●材料处理应⽤●建筑设计和⽕灾研究总⽽⾔之,对于模拟复杂流场结构的不可压缩/可压缩流动来说,FLUENT是很理想的软件。
计算流体力学作业FLUENT 模拟燃烧
一、模拟对象描述
圆柱型火焰燃烧器的结构图1所示。
火焰是湍流扩散火焰,在进口的中心处有一个小喷嘴。
甲烷以80m/s的速度从小喷嘴中射入,周围空气以0.5m/s 的速度流入燃烧器,过量空气系数为1.28。
在甲烷和空气之间用一层外墙隔开。
甲烷流动的雷诺数为5700.甲烷与空气的反应采用最常见的单步总包反应,而且认为反应是扩散控制的,因此使用涡耗散模型对其进行模拟。
图1 二维湍流扩散燃烧器中甲烷空气燃烧
二、实例操作步骤
1.利用GAMBIT建立计算区域和指定边界条件类型。
2.利用FLUENT求解器求解
步骤1:网格的相关操作
启动二维FLUENT,在菜单中点击File-Read-Case…,在相应目录中,找到自己生成的gascomb.msh。
点击Grid-Check,检查网格。
点击Grid-Scale…,设定网格尺寸,将网格改为按毫米生成。
点击Grid-Check,检查一下计算域是否正确:X的最大值是1.8,Y的最大值是0.225.然后关闭对话框。
点击Display-Grid…显示网格
步骤2:模型的设定
步骤3:材料属性设定
步骤4边界条件的设定
步骤5:设定初始条件和其他求解控制参数设置
残差随迭代逐渐收敛情况步骤6:结果显示
温度等值线云图。