# Reacting Flows using PDF Methods

 Register Blogs Members List Search Today's Posts Mark Forums Read

 November 21, 2005, 22:53 Reacting Flows using PDF Methods #1 Gede Guest   Posts: n/a I'm trying to develop numerical method in reacting flow using PDF method. I am confused with the SIMPLE algorithm. From the mixture fraction PDF I will obtain density, temperature, and reacting scalars. Do I still need to execute SIMPLE algorithm to obtain velocity field or just solve the momentum equation , since I can express the pressure field through state equation?

 November 22, 2005, 06:58 Re: Reacting Flows using PDF Methods #2 Salvador Guest   Posts: n/a Remember compressible versus incompressible. In most turbulent combustion applications you use low Mach number assumption, which allows to un-coupled the momentum from the energy equation and therefore use SIMPLE type algorithms (see the Wiki for tons of examples). The pressure you solve with SIMPLE in the momentum equation is NOT the thermodynamic pressure. If you are dealing with a pure compressible flow (Ma>0.7) then you have to change the numerical method to something like Riemann solvers or so. In this case methods based on scalar PDFs are questionable (they may assume constant pressure internally) and you must check thet the density obtained satisfies continuity.

 November 22, 2005, 09:53 Re: Reacting Flows using PDF Methods #3 Confused. Guest   Posts: n/a Hello Salvador, Could you explain further what you mean by: 'The pressure you solve with SIMPLE in the momentum equation is NOT the thermodynamic pressure'? What is the difference between these two 'pressures' you seem to be referring to? Furthermore, the say that the density obtained using a PDF transport equation in high Mach No. flows does not necessarily satisify continuity: Does this mean that the PDF equation is not conservative? How can we make the density field continuity satisfying? I am primarily interested in an answer to the first of my questions.

 November 22, 2005, 12:55 Re: Reacting Flows using PDF Methods #4 Salvador Guest   Posts: n/a 'reduced' pressure, which contains a 1/3 div(v) inside. This term appears from the stress modelling and dissapears if no variations of density are used. The pressure solved in the momentum equation is the pressure created/needed by moving/to move the fluid. The absolute value of pressure is of no significance in incompressible flow only differences. Pressure is only used to satisfy continuity. In compressible flow you need the continuity equation and and an energy equation as well. Thre are many PDF type methods, if you are using a flamelet solution Y(Z). This flamlet profiels are obtained a priori assuming certain Pressure (thermodynamic) and thereore giving certain profiles rho(Z). If the pressure used a priori and the actual are not the same the density field will not satisfy continuity. In many open flames pressure is assumed constant so low Mach number assumptions are valid.

 November 23, 2005, 23:01 Re: Reacting Flows using PDF Methods #5 Gede Guest   Posts: n/a So, is it valid when I calculate the pressure field from the state equation using mixture fraction transport equation (Z transport equation)? Then I calculate the velocity field from the momentum equation using the pressure obtained before. Thank you for the response.

 November 24, 2005, 07:05 Re: Reacting Flows using PDF Methods #6 Salvador Guest   Posts: n/a Well the velocit field obtained will not satisfy continuity. Best (low Mach) --------- 1)Compute u,v,w 2) Then P (or pressure correction depends method) 3) Correct velocities u.v,w (now they wil satisfy continuity) 4) Advance Z (and Z") using new u,v,w 5) Compute T(Z),rho(Z) using flamelets Go to (1) and repeat until steady state or until you want to advance the solution one time step, in the latter yopur new rho will enter in equation 2) torugh the term drho/dt. If your code is unsteady and have large denity changes ypu may have isntability problems

 November 29, 2005, 22:13 Re: Reacting Flows using PDF Methods #7 Gede Guest   Posts: n/a thanks for the details

 Thread Tools Display Modes Linear Mode

 Posting Rules You may not post new threads You may not post replies You may not post attachments You may not edit your posts BB code is On Smilies are On [IMG] code is On HTML code is OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post chandramurthy OpenFOAM Running, Solving & CFD 0 September 10, 2009 03:34 Hanxsi FLUENT 0 February 26, 2006 21:51 kiran FLUENT 0 March 12, 2005 05:24 Karthick FLUENT 2 September 10, 2004 09:44 Mounir Main CFD Forum 1 January 21, 1999 11:08

All times are GMT -4. The time now is 04:20.