
[Sponsors] 
EOS and phase change on reactingTwoPhaseEulerFoam 

LinkBack  Thread Tools  Search this Thread  Display Modes 
January 26, 2016, 14:11 
EOS and phase change on reactingTwoPhaseEulerFoam

#1 
New Member
Farhad N.
Join Date: Apr 2010
Posts: 6
Rep Power: 14 
I am struggling with adding a new equation of state to be used with reactingTwoPhaseEulerFoam, since I am dealing with a cryogenic liquid, which is not an ideal liquid, and Peng Robinson is only available for gas. My idea was to replace the Thermo class and write my own class (myThermo) and get rid any possible complications.
I am dealing with phasechange defined in constant/phaseProperties and I chose Antoine equations to consider phasechange at the boiling point according to pressure fluctuations. I don't find, in the code, where phasechange is programmed. Any idea???? Thanks 

May 16, 2016, 13:16 

#2  
New Member
Simone Colucci
Join Date: Mar 2016
Location: Pisa (Italy)
Posts: 23
Rep Power: 8 
Quote:
I'm also struggling with the phasechange in reactingTwoPhaseEulerFoam. Few things I have understood now:  the mass transfer term is somehow similar to the term used in Fluent (Eqn. 17522)  the mass transfer coefficient is defined in reactingEulerFoam/interfacialCompositionModels/massTransferModels/Frossling or .../SphericalMassTransfer.  The equilibrium concentration is in reactingEulerFoam/interfacialCompositionModels/interfaceCompositionModels/Henry or .../Raoult. Anyway, I still do not understand 1) where is the mass transfer term written in the code ? 2) where is the Lewis number used ? 2) why in the Henry law I do notsee any pressure? Any idea?Thanks Simone 

May 24, 2016, 08:24 

#3  
New Member
Simone Colucci
Join Date: Mar 2016
Location: Pisa (Italy)
Posts: 23
Rep Power: 8 
Quote:
 the mass transfer coefficient is calculated by the member function massTransfer in reactingEulerFoam/phaseSystems/phaseSystems/InterfaceCompositionPhaseChangePhaseSystem /InterfaceCompositionPhaseChangePhaseSystem.C This function is called twice from YEqns.H, first before solving the transport equations for components in each phase, Y1iEqn and Y2iEqn, and then at the end to update the mass transfer term. // call function massTransfer autoPtr<phaseSystem::massTransferTable> massTransferPtr(fluid.massTransfer()); ... // call function massTransfer fluid.massTransfer(); // updates interfacial mass flow rates During the first call, the mass transfer term is added to the Y matrix equation as fvm, // Implicit transport through the phase *eqns[name] += phase.rho()*KD*Yf  fvm::Sp(phase.rho()*KD, eqns[name]>psi()); The second call updates the value stored in dmdt adding to dmdtExplicit the part depending on the solution of the Y equation, already solved, dmdt = dmdtSign*phase.rho()*KD*eqns[name]>psi(); Note that dmdt, before entering into the loop, is initialized with the value of dmdtExplicit *this>dmdt_[pair] = *this>dmdtExplicit_[pair]; in this way the explicit part used for calculating dmdt corresponds to the explicit part used in the YEqn.  KD in the phase change term is given by K*D, where K is the mass transfer coefficient calculated in reactingEulerFoam/InterfacialCompositionModels/massTransferModels/SphericalMassTransfer/spericalMassTransfer.C or ../Frossling/Frossling.C and D the mass diffusivity calculated in in reactingEulerFoam/InterfacialCompositionModels/InterfaceCompositionModels/InterfaceCompositionModel.C. Any reference, idea ... about the interface Composition Model? I don't know exactly what is that. Simone 

May 1, 2017, 19:03 

#4 
Senior Member
Join Date: Jan 2013
Posts: 372
Rep Power: 12 
Dear Simone,
I think the interface Composition Model is used to estimate the species mass fraction under saturation at the phase interface which will be used for the mass transfer model. What kind of problems are you solving? gasliquid and gasparticles? 

May 10, 2017, 08:29 

#5 
New Member
Simone Colucci
Join Date: Mar 2016
Location: Pisa (Italy)
Posts: 23
Rep Power: 8 
Hi,
I'm solving for gasliquid only... 

January 2, 2020, 09:13 
Henry's law constant

#6 
Member
Rosario Arnau
Join Date: Feb 2017
Location: Spain
Posts: 57
Rep Power: 7 
Hi Foamers,
I would like to ask you two things about the interface Composition models in this solver: 1. In some tutorials, the interfaceComposition model is just specified for one phase. Example: tutorials\multiphase\reactingTwoPhaseEulerFoam\RAS \bubbleColumnEvaporatingReacting\constant Code:
interfaceComposition ( (gas in liquid) { type Saturated; species ( H2O ); Le 1.0; saturationPressure { type ArdenBuck; } } ); Code:
interfaceComposition ( (gas in liquid) { type Saturated; species ( water ); Le 1.0; saturationPressure { type ArdenBuck; } } (liquid in gas) { type Henry; species (air ); k ( 1.492e2 ); Le 1.0; } ); Code:
(liquid in gas) { type Henry; species (air ); k ( 1.492e2 ); Le 1.0; } Thanks in advance! 

January 27, 2022, 03:00 

#7 
New Member
Sk Hossen Ali
Join Date: Jul 2021
Location: India
Posts: 8
Rep Power: 3 
Dear Simone,
You seem to understood quite a bit about the mass transfer in interfaceCompositionPhaseChangeTwoPhaseSystem (ICPCTPS) and I also do agree with your points. I am struggling for the same questions in case of the thermalPhaseChangeTwoPhaseSystem (TPCTPS). As the description of the solver TPCTPS suggest depending upon the heat transfer between the phases the amount of mass evaporation/condensation happens. When I ran few cases to simulate evaporation due to heat transfer from hot gases to the liquid (at saturation condition), it never happens. After few times steps the temperature of gases near the interface seems to go below saturation temperature and opposite happens (i.e., condensation from hot gases to liquid). I have found no matter how much is the temperature of liquid or gases near the interface always the condensation happens. The phases as liquid (occupied in the lower half of domain) saturated water, gas (occupied in the upper half of domain) mixture of air and water vapor(initially zero amount) I am attaching my whole case which I am running in openFoam2106 version. 1. How is the mass transfer term calculated based on the thermal exchange between the phases? Where are those terms in code like you have shown for the case of ICPCTPS. https://drive.google.com/file/d/1DhB...ew?usp=sharing My main goal is to  simulate evaporation to happen due to heat transfer from hot gases to the liquid interface. Any help would greatly appreciated. Thanks in advance. 

Tags 
equation of state, evaporation, openfoam, phase change, reactingtwophaseeulerfoam 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Calculation of the Governing Equations  Mihail  CFX  7  September 7, 2014 07:27 
Solid EOS  adam  Main CFD Forum  5  March 18, 2008 14:12 