CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

EOS and phase change on reactingTwoPhaseEulerFoam

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

Like Tree6Likes
  • 2 Post By S.Colucci
  • 3 Post By S.Colucci
  • 1 Post By openfoammaofnepo

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 26, 2016, 14:11
Default EOS and phase change on reactingTwoPhaseEulerFoam
  #1
New Member
 
Farhad N.
Join Date: Apr 2010
Posts: 6
Rep Power: 14
farhhad is on a distinguished road
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 phase-change defined in constant/phaseProperties and I chose Antoine equations to consider phase-change at the boiling point according to pressure fluctuations. I don't find, in the code, where phase-change is programmed. Any idea????

Thanks
farhhad is offline   Reply With Quote

Old   May 16, 2016, 13:16
Default
  #2
New Member
 
Simone Colucci
Join Date: Mar 2016
Location: Pisa (Italy)
Posts: 23
Rep Power: 8
S.Colucci is on a distinguished road
Quote:
Originally Posted by farhhad View Post
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 phase-change defined in constant/phaseProperties and I chose Antoine equations to consider phase-change at the boiling point according to pressure fluctuations. I don't find, in the code, where phase-change is programmed. Any idea????

Thanks
Hi Farhad,
I'm also struggling with the phase-change in reactingTwoPhaseEulerFoam. Few things I have understood now:
- the mass transfer term is somehow similar to the term used in Fluent (Eqn. 17-522)
- 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
rarnaunot and HHK_1702 like this.
S.Colucci is offline   Reply With Quote

Old   May 24, 2016, 08:24
Default
  #3
New Member
 
Simone Colucci
Join Date: Mar 2016
Location: Pisa (Italy)
Posts: 23
Rep Power: 8
S.Colucci is on a distinguished road
Quote:
Originally Posted by S.Colucci View Post
Hi Farhad,
I'm also struggling with the phase-change in reactingTwoPhaseEulerFoam. Few things I have understood now:
- the mass transfer term is somehow similar to the term used in Fluent (Eqn. 17-522)
- 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
I can now answer to some of my questions:

- 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
S.Colucci is offline   Reply With Quote

Old   May 1, 2017, 19:03
Default
  #4
Senior Member
 
Join Date: Jan 2013
Posts: 372
Rep Power: 12
openfoammaofnepo is on a distinguished road
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? gas-liquid and gas-particles?
rarnaunot likes this.
openfoammaofnepo is offline   Reply With Quote

Old   May 10, 2017, 08:29
Default
  #5
New Member
 
Simone Colucci
Join Date: Mar 2016
Location: Pisa (Italy)
Posts: 23
Rep Power: 8
S.Colucci is on a distinguished road
Hi,
I'm solving for gas-liquid only...
S.Colucci is offline   Reply With Quote

Old   January 2, 2020, 09:13
Default Henry's law constant
  #6
Member
 
Rosario Arnau
Join Date: Feb 2017
Location: Spain
Posts: 57
Rep Power: 7
rarnaunot is on a distinguished road
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;
        }
    }
);
And in others two models are specified:

Code:
interfaceComposition
(
    (gas in liquid)
    {
        type Saturated;
        species ( water );
        Le 1.0;
        saturationPressure
        {
            type ArdenBuck;
        }
    }

    (liquid in gas)
    {
        type Henry;
        species (air  );
        k ( 
		1.492e-2
		);
        Le 1.0;
    }
);
2. What are the units of the Henry's constant specied in:
Code:
   (liquid in gas)
    {
        type Henry;
        species (air  );
        k ( 
		1.492e-2
		);
        Le 1.0;
    }
Guess that they are dimensionless...

Thanks in advance!
rarnaunot is offline   Reply With Quote

Old   January 27, 2022, 03:00
Default
  #7
New Member
 
Sk Hossen Ali
Join Date: Jul 2021
Location: India
Posts: 8
Rep Power: 3
Sk Hossen Ali is on a distinguished road
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 openFoam-2106 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.
Sk Hossen Ali is offline   Reply With Quote

Reply

Tags
equation of state, evaporation, openfoam, phase change, reactingtwophaseeulerfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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 Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


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


All times are GMT -4. The time now is 07:08.