CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

KIVA-conserved scalar

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

LinkBack Thread Tools Display Modes
Old   August 22, 2005, 02:17
Default KIVA-conserved scalar
Posts: n/a
Dear all,

I am facing a problem related with coupling of PDF-conserved scalar algorithm with KIVA code. I do hope any of you could help me by giving some hints.

Dealing with diffusion flame, I have succesfully implemented the transport equations of mean and variance of mixture fraction into KIVA. Consequently, a beta-function PDF of mixture fraction can also be obtained. Based on this PDF, I could also obtain the mean reactive scalars (species concentration, density, temperature, and enthalpy). Until this point, I did well.

However, in order to obtain the correct flow field, for the flow equation in KIVA I have to use the mean density obtained by PDF approach. So I specified the mean density in KIVA with the value I obtained using PDF, and I also deactivated subroutine ysolve (to calculate the species mass fraction) since I think this task is also taken over by PDF.

The problem is, when I tried to execute the code, temperature overflow occured. And after checking the code, I found that this temperature overflow is related with subroutine state, where the state equation is solved. And subroutine state is related with subroutine tsolve, where mean temperature is solved based on energy equation (original KIVA scheme), not by PDF integration.

So my questions would be:

1. Which subroutines do I have to disregard when dealing with pdf method in KIVA?

2. In this case, do I still need iteration (known as big iteration in KIVA) to calculate the flow field? As far as I can understand, given that I have the mean density obtained by pdf, I just have to use this mean density in subroutine vsolve (where velocity field is solved) and deactivate tsolve, psolve, and state. In addition, the temperature is defined by the mean temperature I obtained by PDF integration. Is this correct?

3. Do you have any suggestion on this problem?

Thank you very much for your kind help.



  Reply With Quote


Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
dieselFoam problem!! trying to introduce a new heat transfer model vivek070176 OpenFOAM Programming & Development 10 December 24, 2014 00:48
Solving for an additional species CO in coalChemistryFoam N. A. OpenFOAM Programming & Development 2 November 13, 2014 09:43
Specifying nonuniform boundary condition maka OpenFOAM Running, Solving & CFD 59 October 22, 2014 14:52
CFX12 rif errors romance CFX 4 October 26, 2009 14:41
PDF - conserved scalar approach in KIVA vega Main CFD Forum 0 July 5, 2005 02:49

All times are GMT -4. The time now is 18:22.