**BrineProp.Brine3salts3gas** is a medium package that, based on Brine.BrineProp.PartialBrineMultiSaltMultiGasTwoPhase, defines a brine property model with 3 salts (NaCl, KCl, CaCl_{2}) and 3 gases (CO_{2}, N_{2}, CH_{4}), which are the main constituents of the geofluid in Gross Schoenebeck, Germany.

It was used for the calculations documented in this PhD thesis.

As it is based on Modelica.Media, the usage differs little from the usage of the two-phase water model:

Create an Instance of the Medium (optionally deactivating range checks, for all options see .PartialFlags ):

package Medium = Brine_Duan_Multi_TwoPhase_ngas_3;

Create an Instance of Medium.Baseproperties:

Medium.BaseProperties props;Use the BaseProperties model to define the actual brine composition(Xi or X), to define the thermodynamic state and calculate the corresponding properties.

props.p = 1e5; props.T = 300; props.Xi = {0.08, 0.004, 0.12, 1-4, 7e-4, 6e-005} "NaCl, KCl, CaCl2, CO2, N2, CH4" d = props.d;

See BrineProp.Examples.BrineProps2phase for more usage examples.

All calculated values are returned in SI units and are mass based.

The enthalpy is calculated as a mass fraction weighted average of the enthalpies of the two phases.

h = x·h_g + (1-x)·h_l

Enthalpy of the gas phase is modelled as the enthalpy of an ideal mixture of ideal gases, i.e. it is calculated as the mass weighted average of the individual gas enthalpies including water.

h_g = sum(h"_{i}·X"_{i})

The individual gas enthalpies are calculated using the functions for ideal gases in Modelica.Media.IdealGases.SingleGases.

h_l = h_{Driesner}+sum(H_{iapp·bi)}

The apparent molar enthalpies of KCl and CaCl2 are calculated from apparent molar heat capacities which are calculated from a 2D fit for data from literature. The contributions of MgCl2 and SrCl2 are neglected, due to their small influence in the GrSbk fluid.

The total density d of the fluid is calculated by combining the densities of both phases (d_g and d_l) according to their volume fractions. The gas phase is assumed to be an Density of the gas phase is assumed to be an ideal mixture of ideal gases.

d = 1/(x/d_g + (1 - x)/d_l)

Helmholtz Centre Potsdam GFZ German Research Centre for Geosciences

Telegrafenberg, D-14473 Potsdam

Germany

francke@gfz-potsdam.de

Name | Description |
---|---|

fugacity_H2O | Calculation of fugacity coefficient according to (Duan 2003) |

setState_pTX | to avoid check error |

setState_phX | to avoid check error |

solubilities_pTX | solubility calculation, returns gas concentration in kg/kg H2O |

density_liq_pTX | |

density_gas_pTX | |

specificEnthalpy_liq_pTX | |

specificEnthalpy_gas_pTX | |

dynamicViscosity_liq | |

dynamicViscosity_gas | |

saturationPressures | |

thermalConductivity | Thermal conductivity of water |

surfaceTension | |

specificHeatCapacityCp_liq | |

specificHeatCapacityCp_gas | calculation of gas specific heat capacity |

Generated at 2024-03-03T19:16:20Z by OpenModelicaOpenModelica 1.22.1 using GenerateDoc.mos