FAIRDO BG example 3.4.cellml

Bond graph example: A simple biochemical reaction

For biochemical systems qⁱ_j is the molar quantity of a chemical species i in solution in compartment j and uⁱ_j is its chemical potential. A biochemical reaction in one compartment ('c'), with forward and reverse affinities A^f and A^r, is shown below along with its bond graph representation. Note the usual representation for energy storage, but the dissipative reaction term (Re_r) now requires explicit dependence on the chemical potential of the reactants and products, rather than depending only on the potential across the dissipation 'resistor'. Therefore, the corresponding Re component must have two ports, in contrast to the standard one-port R component. Flow 1:nodes are therefore needed on either side of the reaction, when multiple species are involved, to provide the two ports.

A biochemical reaction where A^f and A^r are the forward and reverse affinities (a), and its bond graph representation (b). The 0:nodes represent the chemical potentials uⁱ₁ and uⁱ₂ associated with the chemical species qⁱ₁ and qⁱ₂ stored in compartments 1 and 2. The reaction Re_r (the dissipative term for a biochemical system) has associated chemical flux v_r that depends explicitly on the chemical potentials of the reactant (q¹_c) and product (q²_c) since in this simple case u^f_R = A^f = u¹_c and u^r_R = A^r = u²_c.

The Views Available menu to the right provides various options to explore this model here in the Physiome Model Repository. Of particular interest is the Launch with OpenCOR menu item, which will load the simulation experiment shown below directly into the OpenCOR desktop application.

Showing the result of launching the simulation experiment from this exposure in OpenCOR and executing the simulation.