
[Sponsors] 
June 10, 2011, 09:25 
Evaporation model in dieselSpray

#1 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
hii evrybody,
i want to implement a new evaporation model in dieselSpray, the are two models; standartEvaporation model and RutlandFlash boil model wich are both assuming the D2 law. i want to implement a 0D model that calculate the equilibrium state of the droplet after evaporation, i have a c++ code that do this but the implimentation is not to be evident fo me. thank's... Last edited by Tag; June 10, 2011 at 09:57. 

June 10, 2011, 10:20 

#2 
Super Moderator
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 21 
formulate the problem in a way that it can be written as
Code:
d m_d  m_d =   dt tau or relaxationTime in http://foam.sourceforge.net/doc/Doxy...tionModel.html 

June 14, 2011, 03:14 

#3 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
Thank you for your answer Niklas, have you an idea of what the Sevap in the Rho equation contains ?
Sevap += dieselSpray.evaporationSource(i); My 0D model have not time for the moment, it caculate only the properties of droplets at the equilibrium state, it doesn't assume a D2 law or Abramzon model, i use the secant methode to calculate these properties. So if i know what Sevap contains, i think that i can implement it.. thank you 

June 14, 2011, 04:03 

#4 
Super Moderator
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 21 
I suspected that you would say that
actually you can relate your model to a characteristic time and use the existing framework to add your new evaporation model. in FOAM the mass of the droplet is calculated/updated according to this equation Code:
(m^n  m^{n1}) m^n  =   (1) dt tau dt is the integration step and tau is the characteristic evaporation time (evaporation model) if you know that the mass of the drop is updated according to this equation, you can calculate the new mass, m^n, according to your equilibrium conditions and related that to tau. like this: according to (1) you will get Code:
m^n K*dt  = 1/( 1 + dt/tau) = K => tau =  m^{n1} 1  K This way you dont have to calculate the source terms and bother with Sevap as it will be handled automatically. 

June 14, 2011, 05:24 

#5 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
thank you Niklas,
I think it's easier to do that, I will keep you informed on the progress of my implementation. i will haves somme questions about the mass fraction of the liquids phase and gazous phase, temperature, pressure modification and updating ... Thank you again Niklas.... Noureddine. 

June 24, 2011, 10:32 

#6 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
Hii Niklas,
I have successfully implemented my 0D model in OpenFoam as you suggest. for the case of standartEvaporationModel and my 0D evaporation model: i want to know how the temperature and other properties like volume fraction and the diameter of the droplets are calculated and updated, I suspect that the B term in evaporation model is considered constant. Thank you. 

June 25, 2011, 08:00 

#7 
Super Moderator
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 21 
Great,
Im not following, which Bterm? if you want to find out how these properties are updated, just check the updatePropertiesfunction in parcel.C. 

June 29, 2011, 07:55 

#8 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
Dear Niklas,
B is the Spalding mass number. I have to give the initial mass fractions of the liquide, gaz, N2 and O2 (air) and their temperatures in cell to my 0Dmodel that calculate the equilibrium temperature and then the final mass fractions of these species in the same cell, in order to estimate the mass ratio of the liquide : m (n)/m(n1)= Yl(n)/Yl(n1) (Yl mass fraction of liquide). the problème that i don't found yet these propertise in openfoam. Can you give me an idea to how i do that ? Thank's..... 

July 4, 2011, 05:27 

#9 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
Hi,
where can i find the total mass, the gas temperature, the fractions of different species in cell(i), i had to use these properties in one class that I've created in dieselSpray SubModels. Thank U. 

July 4, 2011, 05:59 

#10 
Super Moderator
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 21 
Have you checked parcel.C?
All is in there. if you want the mass of the cell you have to get the density and cell volume Code:
scalar cellV = sDB.mesh().V()[celli]; scalar rho = sDB.rho()[celli]; scalar cellMass = rho*cellV; Code:
scalar Tg0 = sDB.TInterpolator().interpolate(position(), celli, facei); Code:
label j = sDB.liquidToGasIndex()[i]; const volScalarField& Yj = sDB.composition().Y()[j]; scalar Yfg0 = Yj[celli]; 

July 4, 2011, 06:18 

#11 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
Hi Niklas,
thank you for your answer, i suspected that and you confirmed me thses information, what about the mass fraction of vapor ? 

July 4, 2011, 06:20 

#12 
Super Moderator
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 21 
again...
vapor mass fraction of liquid index i (which has index j among the gas species) Code:
label j = sDB.liquidToGasIndex()[i]; const volScalarField& Yj = sDB.composition().Y()[j]; scalar Yfg0 = Yj[celli]; 

July 4, 2011, 06:24 

#13 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
oops ... thank you so much Mr Niklas.


July 4, 2011, 06:25 

#14 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
last question,
know you how the spalding number that used in the evaporation model (ysy/1ys) (somme thing lick that.) is calculated in openFoam? 

July 4, 2011, 08:08 

#15 
Super Moderator
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 21 
standardEvaporationModel.C
B = 1 + Xratio Code:
/* (pressure  partialFuelVaporPressure)/ (pressure  pressureAtSurface) = 1 + Xratio if the pressure @ Surface > pressure this lead to boiling and Xratio > infinity (as it should) ... this is numerically nasty NB! by N. Nordin X_v,s = (p_v,s/p) X_v,d where X_v,d = 1 for single component fuel according to eq (3.136) in D. Clerides Thesis */ scalar Xratio = (Xs  Xf)/max(SMALL, 1.0  Xs); if (Xratio > 0.0) { lgExpr = log(1.0 + Xratio); } 

July 4, 2011, 08:21 

#16 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
thank you Niklas,
i have already seen this in standardEvaporationModel.C, Xs is a variable of the relaxationTime function, but when this function is used we must give the Xs that is, i think is calculated by another function or expression, this is what i want to know. thank you. 

July 5, 2011, 02:27 

#17 
Super Moderator
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 21 
Ah, OK.
Thats done in parcel/setRelaxationTimes.C However, Xs just calls the function which is defined in the liquidMixture class thermophysicalModels/liquidMixture/liquidMixture/liquidMixture.C which just looks like this Code:
Foam::scalarField Foam::liquidMixture::Xs ( const scalar p, const scalar Tg, const scalar Tl, const scalarField& xg, const scalarField& xl ) const { scalarField xs(xl.size(), 0.0); // Raoult's Law forAll(xs, i) { scalar Ti = min(TrMax*properties_[i].Tc(), Tl); xs[i] = properties_[i].pv(p, Ti)*xl[i]/p; } return xs; } 

July 5, 2011, 08:44 

#18 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
Hi Niklas,
To use the molar fraction and gas temperature in my "model0D" class, i have to use an instance of the "spray" class, so what kind of inheritance can do that easily ? Thank's.... 

July 5, 2011, 09:04 

#19 
Super Moderator
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 21 
Im not sure I understand your question, but the access calls to those variables are posted in my previous mails.


July 5, 2011, 09:20 

#20 
New Member
Noureddine TAGUELMIMT
Join Date: Mar 2011
Location: Rouen, France
Posts: 21
Rep Power: 7 
Hi,
This is a part of a header of the class i created called model0D, #ifndef model0D_H #define model0D_H namespace Foam { // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // class spray; // line added to call the spray class class model0D { And, in the model0D.C i added the properties lick gas temperature : spray mysDB; scalar Tg0 = mysDB.TInterpolator().interpolate(position(), celli, facei); but it doesn't work!!, i suspect that i haven't made the necessary operations to make correct calling between the two classes . Thank you.... 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Superlinear speedup in OpenFOAM 13  msrinath80  OpenFOAM Running, Solving & CFD  18  March 3, 2015 06:36 
Water subcooled boiling  Attesz  CFX  7  January 5, 2013 04:32 
multiphaseInterFoam for RAS turbulence model  chiven  OpenFOAM Bugs  7  August 23, 2011 02:52 
multi fluid mixture model issue  rystokes  CFX  3  August 9, 2009 19:13 
Advanced Turbulence Modeling in Fluent, Realizable kepsilon Model  Jonas Larsson  FLUENT  5  March 13, 2000 04:27 