|
[Sponsors] |
July 5, 2005, 02:49 |
PDF - conserved scalar approach in KIVA
|
#1 |
Guest
Posts: n/a
|
Dear all,
I am trying to simulate turbulent non-premixed flame with equilibrium chemistry. So I`m using pdf-conserved scalar approach (equilibrium flamelet). I have implemented PDF-mixture fraction in KIVA code and obtained all reactive scalars (temperature, enthalpy, density, and species concentration) as functions of the mixture fraction. I already compared the density (obtained using pdf) with the density obtained using cold-flow KIVA simulation for several time steps, and they are in excellent agreement. As far as I understand, the density must be returned back to KIVA to calculate the turbulent flow field, so I can obtain new values of k and epsilon to be used in the transport equation of variance of mixture fraction which is required to construct pdf at the next time step. The current problem is, after I coupled the new density with KIVA (say.. if ro(i4) is density of original kiva, and rbar(i4) is density obtained using pdf, I simply wrote ro(i4)=rbar(i4) ), the code crashed due to temperature overflow. I am trying to trace back the error in the code list, but I was overwhelmed with the complexity of the code. So the questions are: 1. How can I couple the density (obtained using pdf) with kiva? 2. Since original KIVA also calculates the species concentration, internal energy, and temperature, but I already obtained them using pdf, do I have to deactivate all of the KIVA subroutines associated with these scalars? Thank you for your kind help. regards, vega |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 23, 2014 23:48 |
Specifying nonuniform boundary condition | maka | OpenFOAM Running, Solving & CFD | 59 | October 22, 2014 14:52 |
KIVA-conserved scalar | vega | Main CFD Forum | 0 | August 22, 2005 02:17 |
pdf method in kiva | vega | Main CFD Forum | 0 | May 16, 2005 03:45 |
PDF approach (urgent!!!) | anand | FLUENT | 9 | September 27, 2003 00:29 |