[NJG]

Hello Community,

I am presently in the middle of trying to alter the bouyantBoussinesqPimpleFoam solver (I will call this the "reference solver" below) to handle not only convective heat transfer but changes in air composition as well. The idea is to treat things as incompressible in the typical way the Boussinesq approximation does with density only affecting buoyant body forces but have this buoyant force be the result of density change both from temperature change and CO2 mole fraction. I want to model a single occupant within a closed room breathing out warm relatively CO2 rich air.

Things so far are more or less coming along. I basically just copied the base solver and started making my own additions. For a first attempt, I am just duplicating the relationships for thermal diffusivity as molecular diffusivity follows pretty much identical forms. As such I have added x, Dx, and Dxt variables to match the T, kappa, and kappat for temperature. I also created a new wall function just like the kappatJayatillekeWallFunction but for Dxt (based on the Schmidt number, turbulent and laminar). I have added a molecular diffusion equation. Etc. etc.

So, I have a bunch of questions, but I will just start in the middle with what I am dealing with now.

Looking in the reference solver tutorial files, I have a T.org file in the 0 folder as well as a T file. The T file has a long list of temperature values, one for every finite volume if I understand correctly. My first question is a) how do I populate this T file with initial values without actually going through and adding them all by hand (my final run case will have ~10^6 volumes, so that would suck), and b) what is T.org and is it related to this?

A follow up question is that I plan to have an x (CO2 mole fraction) file just like I have a T file. Would I also need an x.org file then, too?

Thanks for the help! If I get this running I will make it available to anyone interested in playing with it.

-NG