Skip to content

Commit

Permalink
delet v/v ratio, it was there too often!
Browse files Browse the repository at this point in the history
  • Loading branch information
joschwar committed Sep 27, 2021
1 parent d49fd8c commit 5b7a542
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion Henry_model_incubation.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -184,7 +184,7 @@ $$\frac{dCO_{2(aq)}}{dt} = k1 \cdot C - kh1 \cdot CO_{2(aq)} + kh2 \cdot CO_{2(g

And the $CO_2$ in the gas phase $\frac{dCO_{2(g)}}{dt}$ outgassing from the water phase $kh1 \cdot CO_{2(aq)}$. The ratio $\frac{V_w}{V_g}$ accounts for the fact, that the volume of the liquid reservoir is different from the gas reservoir, so that the *concentration* change is different when a $CO_2$ molecule passes the phase boundary.

$$\frac{dCO_{2(g)}}{dt} = \frac{V_w}{V_g} \cdot \left( kh1 \cdot CO_{2(aq)} \cdot \frac{V_w}{V_g} - kh2 \cdot CO_{2(g)} \right)$$
$$\frac{dCO_{2(g)}}{dt} = \frac{V_w}{V_g} \cdot \left( kh1 \cdot CO_{2(aq)} - kh2 \cdot CO_{2(g)} \right)$$

These differential equations can then be solved with the `ode()` function from the package `deSolve` in a wrapper function:

Expand Down

0 comments on commit 5b7a542

Please sign in to comment.