CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > FLUENT

Implementing local chemical equilibrium

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 9, 2020, 10:33
Default Implementing local chemical equilibrium
  #1
New Member
 
Markus Pieber
Join Date: Jul 2015
Posts: 9
Rep Power: 10
mpieber is on a distinguished road
Dear community members,

I am currently trying to implement the local chemical equilibrium approach in a reactive transient Fluent simulation, where each species is represented by a UDS. Target is that both species (UDS) transport equations and chemical equilibrium condition are fulfilled at the end of each time step, i.e. the species should be in local chemical equilibrium at any time.

Does anyone know a way (e.g. using UDFs) to implement local chemical equilibrium in a transient reactive Fluent simulation (when the species are represented by UDSs)?
Many thanks in advance for any help!

Best regards,
Markus
mpieber is offline   Reply With Quote

Old   March 10, 2020, 08:09
Default Chemical Equilibria
  #2
Senior Member
 
vinerm's Avatar
 
Vinerm
Join Date: Jun 2009
Location: Nederland
Posts: 2,946
Blog Entries: 1
Rep Power: 35
vinerm will become famous soon enough
How would you like to implement the equilibrium? In real systems, it is based on Gibbs Free energy. In simpler terms, you have to ensure that at certain conditions, which is defined by minimum Gibbs Free energy, mass transfer between two scalars settle down to be equal.

So, assuming that the said system of scalars has an equilibrium constant, which may or may not be a function of P and T (and the values of scalars themselves, which will make problem really non-linear), you should be able to define mass transfer based on the equilibrium constant.
__________________
Regards,
Vinerm

PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority.
vinerm is offline   Reply With Quote

Old   March 10, 2020, 08:46
Default
  #3
New Member
 
Markus Pieber
Join Date: Jul 2015
Posts: 9
Rep Power: 10
mpieber is on a distinguished road
Dear Vinerm,

Thanks for your answer.

Assuming there is a functional relationship UDS2,eq = f(UDS1,eq) between UDS1 and UDS2 describing (chemical or adsorption) equilibrium, how could the mass transfer (source term) between UDS1 and UDS2 be defined (or how could the equilibrium condition be enforced otherwise)?

Best regards,
Markus
mpieber is offline   Reply With Quote

Old   March 10, 2020, 08:56
Default Exchange in Multiphase
  #4
Senior Member
 
vinerm's Avatar
 
Vinerm
Join Date: Jun 2009
Location: Nederland
Posts: 2,946
Blog Entries: 1
Rep Power: 35
vinerm will become famous soon enough
You can use either DEFINE_MASS_TRANSFER or DEFINE_LINEARIZED_MASS_TRANSFER functions. Example for DEFINE_MASS_TRANSFER are available at

https://www.afs.enea.it/project/nept...udf/node63.htm
__________________
Regards,
Vinerm

PM to be used if and only if you do not want something to be shared publicly. PM is considered to be of the least priority.
vinerm is offline   Reply With Quote

Reply


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
pimpleDyMFoam computation randomly stops babapeti OpenFOAM Running, Solving & CFD 5 January 24, 2018 05:28
Floating point exception error lpz_michele OpenFOAM Running, Solving & CFD 53 October 19, 2015 02:50
Upgraded from Karmic Koala 9.10 to Lucid Lynx10.04.3 bookie56 OpenFOAM Installation 8 August 13, 2011 04:03
[blockMesh] BlockMeshmergePatchPairs hjasak OpenFOAM Meshing & Mesh Conversion 11 August 15, 2008 07:36
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 18:07


All times are GMT -4. The time now is 15:10.