modelForwardReaction

Chemical Reaction
Diagram of ForwardReaction

Extends from Chemical.Processes.Internal.PartialReactionWithProductsDefinition (Chemical Reaction), Chemical.Interfaces.PartialSolutionSensor, Chemical.Interfaces.ConditionalKinetics (Input of kinetics coefficient vs. parametric kinetics coefficient).

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.StoichiometricNumber[nS]s (from PartialReaction)ones(nS)Stoichiometric coefficients for substrates
Modelica.Units.SI.StoichiometricNumber[nP]p (from PartialReaction)ones(nP)Stoichiometric numbers for products
Chemical.Interfaces.Definitionprocess (from PartialReactionWithProductsDefinition)Chemical.Interfaces.processData(1)Process definition
Realk_forward (from ConditionalKinetics)1Forward rate coefficient (mole-fraction based)
Advanced
Modelica.Units.SI.MolarFlowRaten_flow_reg (from PartialReaction)dropOfCommons.n_flow_regRegularization threshold of mass flow rate
StateSelectn_flowStateSelect (from PartialReaction)StateSelect.defaultState select for n_flow
Modelica.Units.SI.TimeTC (from PartialReaction)dropOfCommons.TCTime constant for electro-chemical potential adaption
Utilities.Units.InertanceL (from PartialReaction)dropOfCommons.LInertance of the flow
Initialization › Molar flow
InitializationMethodsinitN_flow (from PartialReaction)Chemical.Utilities.Types.InitializationMethods.noneInitialization method for n_flow
Modelica.Units.SI.MolarFlowRaten_flow_0 (from PartialReaction)0Initial value for n_flow
Utilities.Units.MolarFlowAccelerationn_acceleration_0 (from PartialReaction)0Initial value for der(n_flow)
General › Ports
IntegernS (from PartialReaction)0Number of substrate types
IntegernP (from PartialReaction)0Number of product types
Products definitions
FirstProductChoicefirstProductFrom (from PartialReactionWithProductsDefinition)FirstProductChoice.ProcessFirst product definition comes from?
Chemical.Interfaces.DefinitionfirstProduct (from PartialReactionWithProductsDefinition)dropOfCommons.DefaultSubstanceFirst product definition as Substance
Chemical.Interfaces.Definition[:]nextProducts (from PartialReactionWithProductsDefinition)fill(((s*ones(nS))/(p*ones(nP)))*dropOfCommons.DefaultSubstance, max(1, nP - 1))Definitions of next products
Chemical solution (of products)
SolutionChoicesolutionFrom (from PartialSolutionSensor)Chemical.Utilities.Types.SolutionChoice.FirstSubstrateChemical solution comes from?
Chemical.Interfaces.SolutionStatesolutionParam (from PartialSolutionSensor)Chemical.Interfaces.SolutionState(phase = Chemical.Interfaces.Phase.Incompressible)Chemical solution state as Parameter
Advanced › Conditional inputs
BooleanuseForwardRateInput (from ConditionalKinetics)falseForward rate coefficient from input?

Connectors

TypeNameDefaultDescription
Chemical.Interfaces.Rear[nS]substrates (from PartialReaction)
Chemical.Interfaces.Fore[nP]products (from PartialReaction)
Chemical.Interfaces.SolutionPortsolution (from PartialSolutionSensor)To connect substance with solution, where is pressented
Modelica.Blocks.Interfaces.RealInputkfInput (from ConditionalKinetics)kf

Components

TypeNameDefaultDescription
Modelica.Units.SI.MolarFlowRaterr (from PartialReaction)Reaction molar flow rate
Chemical.Interfaces.SolutionStatesolutionState (from PartialSolutionSensor)
Realkf (from ConditionalKinetics)Current forward rate coefficient

Revisions

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