modelActivityCoefficient

Calculate activity coefficient for product[1]
Diagram of ActivityCoefficient

Extends from Modelica.Icons.RectangularSensor (Icon representing a linear measurement device).

Information

s1·S1 + .. + snS·SnS <-> p1·P1 + .. + pnP·PnP

By redefinition of stoichometry as vi = -si, Ai = Si for i=1..nS vi = pi-nS, Ai = Pi-nS for i=nS+1..nS+nP

So the reaction can be written also as 0 = ∑ (vi · Ai)

Equilibrium equation

K = product(a(S).^s) / product( a(P).^s ) = product(a(A).^v) 

dissociation constant

ΔrG = ∑ (vi · ΔfGi) = ΔrH - T·ΔrS = -R·T·log(K)

molar Gibb's energy of the reaction

ΔrH = ∑ (vi · ΔfHi)

molar enthalpy of the reaction

ΔrS = ∑ (vi · ΔfSi) = k·log(Δrω)

molar entropy of the reaction

Notations

Ai

i-th substance

vi

stochiometric coefficients of i-th substance

K

dissociation constant (activity based)

a(Ai)=fi*xi

activity of the substance A

fi

activity coefficient of the substance A

xi

mole fraction of the substance A

ΔfHi

molar enthalpy of formation of i-th substance

ΔfGi

molar Gibbs energy of formation of i-th substance

ΔfSi

molar entropy of formation of i-th substance

Δrω

change of number of microstates of particles by reaction

Parameters

TypeNameDefaultDescription
Modelica.Units.SI.TemperatureTsystem.T_ambientTemperature if not useTemperatureInput
Modelica.Units.SI.AmountOfSubstancen1Amount of all substances in solution per one liter of solution if not useTotalAmountOfSubstancesInput
Modelica.Units.SI.Massm1Mass of solvent per one liter of solution
Modelica.Units.SI.StoichiometricNumber[nS]sones(nS)Stoichiometric reaction coefficient for substrates
Modelica.Units.SI.StoichiometricNumber[nP]pones(nP)Stoichiometric reaction coefficients for products
BooleanMolarityBasedtrueif dissociation coefficient is molarity based
RealDissociationCoefficient_MoleFractionBasedif MolarityBased then DissociationCoefficient_MolarityBased/((n/1)^(p*ones(nP) - s*ones(nS))) else DissociationCoefficient_MolalityBased/((n/m)^(p*ones(nP) - s*ones(nS)))K as ratio of mole fractions
RealDissociationCoefficient_MolalityBased((n/m)^(p*ones(nP) - s*ones(nS)))*DissociationCoefficient_MoleFractionBasedK as ratio of molalities in moles per 1 kg of solvent
RealDissociationCoefficient_MolarityBased((n/1)^(p*ones(nP) - s*ones(nS)))*DissociationCoefficient_MoleFractionBasedK as ratio of molar concentration in moles per liter of solution
Conditional inputs
BooleanuseTemperatureInputfalse=true, if temperature is from input instead of parameter
BooleanuseTotalAmountOfSubstancesInputfalse=true, if total amount of substances in solution is from input instead of parameter
Ports
IntegernS0Number of substrates types
IntegernP0Number of products types

Connectors

TypeNameDefaultDescription
Modelica.Blocks.Interfaces.RealInputtemperature_temperatureTemperature
Modelica.Blocks.Interfaces.RealInputtotalAmountOfSubstances_nTemperature
Interfaces.SubstancePort_b[nP]productsProducts
Interfaces.SubstancePort_b[nS]substratesSubstrates
Modelica.Blocks.Interfaces.RealOutputactivityCoeficientActivity coeficient of one product

Components

TypeNameDefaultDescription
Modelica.Fluid.SystemsystemSystem wide properties
Modelica.Units.SI.MolarEnergyDrGFree Gibbs energy of reaction
RealpK= -log10('mole-fraction based dissociation coefficient')

Revisions

2013-2015 by Marek Matejak, Charles University, Prague, Czech Republic