Hi, I am trying to develop liquid phase multi species solver with/without reaction (no combustion).
Can anybody explain me how to define species, species properties in OpenFOAM?
How to create lookup table for species so that each species is assigned some number?
So that I can solve the transport equation for each species like that done in Yeqn.H code in reactingFoam solver.
Mass fraction given by the above YEqn can then be used to calculate mass weighted average mixture properties like:
rhomix = rho1*Y[1] + rho2*Y[2] + rho3*Y[3]

This average mixture properties then can be used to solve the momentum equation.

Reaction effect can be introduced by adding the source term in Yeqn.H later on.

I would like to use mixture viscosity in the solving momentumm equation. For that I need to calculate mixture density (mass weighted). From that I can calculate mixture kinematic viscosity (nu_mix) . To do this, I need access to species mass fraction, species density, species viscosity.
rho_mix = rho1*Y[1]+rho2*Y[2]+rho3*Y[3]
nu_mix= (rho1*nu1*Y[1]+rho2*nu2*Y[2]+rho3*nu3*Y[3])/rho_mix

My question is how can I call individual species properties in for loop to calculate rho_mix and nu_mix?

I am not able to understand the error message. I tried to look into list.H and list .C but could not get the clue.
Can you (anyone) please help me out in how to define the list of scalars and their assignment?

OpenFOAM is free software; you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by the
Free Software Foundation; either version 2 of the License, or (at your
option) any later version.

OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.

You should have received a copy of the GNU General Public License
along with OpenFOAM; if not, write to the Free Software Foundation,
Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA

Yes, it doesnot make any difference. The problem lies in how to fill up the ptrList? I looked at ptrList code. The set function is used for the same.
I am trying with that now. Do you know the correct way/alternate way to fill up the ptrList?

I followed your thread and it seems you've done exactly the same as I am starting to do right now. Is your solver anyway available as a contribution?

If not, what did you change then from the last post on, to make the solver not only compiling but also running?

Best Regards form Zurich to Mumbai,
Tony

Hi Tony,

No, it is not. The current solver does only the multi species part without reaction and heat transfer. It just transports mixture of species frpom one point to another. Also, the solver is not stable when it is running. The initial residual values do not come below 0.1. They fluctuate between 0.12 to 0.73. The solution is not converged. Do you know any reason for this?

Also now I am trying to include heat transfer in the above solver by introducing enthalpy equation. I am confused with thermophysical model framework available in OF. Can you help me out in this or give suggestions on it?

what is the reason for you to create a solver with multicomponent mixture model on your own? You could easily use reactingFoam. I think this solver would fulfill all your needs, as it provides mixture, heat transfer, reactions on/off, but no radiation by the way. But this is not really a problem since it is very easy to implement a radiation model in OpenFoam. If you need help regarding this just let me know.

what is the reason for you to create a solver with multicomponent mixture model on your own? You could easily use reactingFoam. I think this solver would fulfill all your needs, as it provides mixture, heat transfer, reactions on/off, but no radiation by the way. But this is not really a problem since it is very easy to implement a radiation model in OpenFoam. If you need help regarding this just let me know.

Tony

Hi Tony,

Thanks for giving your time. The reactingFoam solver updates the species properties using gas law. I don't want that. I needed a solver to solve the reaction, heat transfer in liquid phase with steady state approach and which would update physical properties as mass weighted average. I think reactingFoam doesn't provide this option. If it can, please let me know the way for this. The reason behind developing this solver is that I first wanted to introduce YEqn and calculate mixture density and use this in UEqn followed by heat transfer addition and reaction stuff.

Thanks again for sharing your thought on this. I really appreciate that. Will contact you if I need any help.

what is the reason for you to create a solver with multicomponent mixture model on your own? You could easily use reactingFoam. I think this solver would fulfill all your needs, as it provides mixture, heat transfer, reactions on/off, but no radiation by the way. But this is not really a problem since it is very easy to implement a radiation model in OpenFoam. If you need help regarding this just let me know.

Tony

Hi Tony,

Thanks for giving your time. The reactingFoam solver updates the physical and thermophysical properties using gas law. I don't want that. I needed a solver to solve the reaction, heat transfer in liquid phase with steady state approach and which would update physical properties as mass weighted average. I think reactingFoam doesn't provide this option. If it can, please let me know the way for this. The reason behind developing this solver is that I first wanted to introduce YEqn and calculate mixture density and use this in UEqn followed by heat transfer addition and reaction stuff.

Thanks again for sharing your thought on this. I really appreciate that. Will contact you if I need any help.

I want to solve the multSpecies problem using OpenFOAM. I have seen your posting in forum regarding this solver. Have you success on this solver? If So, Could you please tell me what are the changes did you made for this solver

I am not able to understand the error message. I tried to look into list.H and list .C but could not get the clue.
Can you (anyone) please help me out in how to define the list of scalars and their assignment?

what is the reason for you to create a solver with multicomponent mixture model on your own? You could easily use reactingFoam. I think this solver would fulfill all your needs, as it provides mixture, heat transfer, reactions on/off, but no radiation by the way. But this is not really a problem since it is very easy to implement a radiation model in OpenFoam. If you need help regarding this just let me know.

Tony

Hi Tony
I'm a new one in OpenFOAM. Recently I'm looking for a method to simulate jet flow of a nozzle using OpenFOAM. So the mixture of combustion gas and air must be considered, and I don't want to calculate the reaction of them. I have read the example case of reatingFoam, but I don't know how to set species and turn reaction off. Can U give me some suggestions?
Thank you very much.