|
[Sponsors] |
March 9, 2020, 10:33 |
Implementing local chemical equilibrium
|
#1 |
New Member
Markus Pieber
Join Date: Jul 2015
Posts: 9
Rep Power: 10 |
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! MarkusBest regards, |
|
March 10, 2020, 08:09 |
Chemical Equilibria
|
#2 |
Senior Member
|
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. |
|
March 10, 2020, 08:46 |
|
#3 |
New Member
Markus Pieber
Join Date: Jul 2015
Posts: 9
Rep Power: 10 |
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 |
|
March 10, 2020, 08:56 |
Exchange in Multiphase
|
#4 |
Senior Member
|
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. |
|
|
|
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 |