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

 September 16, 2010, 05:41 Steady state multi species #1 Member   Hrushikesh Khadamkar Join Date: Jul 2010 Location: Mumbai Posts: 66 Rep Power: 9 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. Any suggestions , Please...

 September 16, 2010, 08:27 #2 New Member   Join Date: May 2010 Location: Cologne Posts: 27 Rep Power: 9 as a beginnig you can use this report http://projekter.aau.dk/projekter/fi...784/Report.pdf

 September 27, 2010, 05:06 individual species properties in for loop to calculate rho_mix and nu_mix #3 Member   Hrushikesh Khadamkar Join Date: Jul 2010 Location: Mumbai Posts: 66 Rep Power: 9 hi, 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? for (lable i=0, i

 September 27, 2010, 09:09 #4 Senior Member   Cyprien Join Date: Feb 2010 Location: Stanford University Posts: 262 Rep Power: 11 hi Hrushi! I think you should defined your rho1,rho2,rho3 as a list of scalar : rho=(rho1,rho2,rho3) hence, rho[0]=rho1 ; rho[1|=rho2 ;rho[2]=rho3 and you can use the following loop : Code: ```rho_mix=0.0 for (int i=0; i

 October 24, 2010, 14:54 #6 Senior Member   Cyprien Join Date: Feb 2010 Location: Stanford University Posts: 262 Rep Power: 11 did you try ptrList instead of List ?

 October 27, 2010, 06:24 Error while execution #7 Member   Hrushikesh Khadamkar Join Date: Jul 2010 Location: Mumbai Posts: 66 Rep Power: 9 Hi, I have modified the simpleFoam solver. I have included specie equation. The code compiles giving no compliation error. But when I run my code, it gives following error message: // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 Reading thermophysicalProperties Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting RAS turbulence model kEpsilon kEpsilonCoeffs { Cmu 0.09; C1 1.44; C2 1.92; alphaEps 0.76923; sigmaEps 1.3; } Starting time loop Time = 1 DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 0.000312273, No Iterations 1 DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 0.000245102, No Iterations 1 DILUPBiCG: Solving for Uz, Initial residual = 0.779043, Final residual = 0.000142658, No Iterations 1 DICPCG: Solving for p, Initial residual = 1, Final residual = 0.00865508, No Iterations 87 time step continuity errors : sum local = 0.00383834, global = 3.55093e-06, cumulative = 3.55093e-06 DILUPBiCG: Solving for epsilon, Initial residual = 0.868617, Final residual = 0.0142968, No Iterations 1 DILUPBiCG: Solving for k, Initial residual = 1, Final residual = 0.0286599, No Iterations 1 hanging pointer, cannot dereference#0 Foam::error:rintStack(Foam::Ostream&) in "/export/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #1 Foam::error::abort() in "/export/apps/OpenFOAM/OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so" #2 main in "/home/na16116/run/run/simpleMyFoam/bin" #3 __libc_start_main in "/lib64/libc.so.6" #4 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/na16116/run/run/simpleMyFoam/bin" From function PtrList:perator[] in file /export/apps/OpenFOAM/OpenFOAM-1.6/src/OpenFOAM/lnInclude/PtrListI.H at line 123. FOAM aborting Abort Can anyone explain what is the reason for this error message? Thanks Hrushikesh

 November 2, 2010, 11:09 #9 Senior Member   Cyprien Join Date: Feb 2010 Location: Stanford University Posts: 262 Rep Power: 11 Code: ```createFields.H Info<< nl << "Reading thermophysicalProperties" << endl; PtrList Y(5);``` Is it normal that you have Y(5) instead of Y(3) ??

 November 3, 2010, 00:13 #10 Member   Hrushikesh Khadamkar Join Date: Jul 2010 Location: Mumbai Posts: 66 Rep Power: 9 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? Hrushi

 November 3, 2010, 02:34 #11 Senior Member   Nakul Join Date: Apr 2010 Location: India Posts: 147 Rep Power: 9 Hi Hrushi, I am also trying to implement something similar. Can you tell me what thermodynamics package you are using for your mixture in your thermoPhysicalProperties dictionary? I am using hsCombustionThermo with "multiComponentmixture" but I am not being able to get it right. There is no example as to how multiple species are read by OpenFOAM while using multicomponent mixture. Could you also tell how are you reading in the multiple species of the mixture? Thanx

 November 14, 2010, 06:40 Multi species solver #12 Member   Hrushikesh Khadamkar Join Date: Jul 2010 Location: Mumbai Posts: 66 Rep Power: 9 Hi, Sorry for late reply (was on Diwali vacation). The problem is solved now. Code is compiling and running. I will update you soon with next development. Currently, I am trying to introduce heat transfer effects in above developed solver. Thanks Hrushikesh

 December 1, 2010, 10:05 #13 Member   Antonio Liggieri Join Date: Aug 2010 Posts: 72 Rep Power: 8 Hello Hrushi, 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

December 2, 2010, 00:42
#14
Member

Join Date: Jul 2010
Location: Mumbai
Posts: 66
Rep Power: 9
Quote:
 Originally Posted by alfa_8C Hello Hrushi, 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?

Hrushi

 December 2, 2010, 06:23 #15 Member   Antonio Liggieri Join Date: Aug 2010 Posts: 72 Rep Power: 8 Hello Hrushi, 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

 December 2, 2010, 06:33 hi #16 New Member   Join Date: Dec 2010 Posts: 2 Rep Power: 0 sry...my reply nt regarding thread..i am new to the forum....iwant to learn gambit and fluent.do u knw some institute for training?

December 2, 2010, 07:37
#17
Member

Join Date: Jul 2010
Location: Mumbai
Posts: 66
Rep Power: 9
Quote:
 Originally Posted by alfa_8C Hello Hrushi, 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.

Hrushi

December 2, 2010, 07:38
#18
Member

Join Date: Jul 2010
Location: Mumbai
Posts: 66
Rep Power: 9
Quote:
 Originally Posted by alfa_8C Hello Hrushi, 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.

Hrushi

August 11, 2011, 02:07
Multi Species in OpenFOAM
#19
New Member

prasanth
Join Date: Jul 2010
Location: Chennai, India
Posts: 17
Rep Power: 0
Hello Hrushi,

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

Regards
Prasanth.

Quote:

January 6, 2015, 04:43
#20
New Member

Yu-sen Niu
Join Date: Nov 2014
Posts: 16
Rep Power: 5
Quote:
 Originally Posted by alfa_8C Hello Hrushi, 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.

 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 Mihail CFX 7 September 7, 2014 06:27 Hrushi OpenFOAM 4 June 5, 2014 13:00 Sas CFX 15 July 13, 2010 08:56 saii CFX 2 September 18, 2009 08:07 mrangitschdowcom OpenFOAM Running, Solving & CFD 3 October 10, 2007 05:21

All times are GMT -4. The time now is 23:29.