Provides methods to calculate thermodynamic properties inlcuding compressibility factors and densities of natural gases. It includes the AGA8 DETAIL and the GERG2008 equations of state described in AGA Report No. 8, Part 1, Third Edition, April 2017.
This crate is a Rust port of NIST's AGA8 code.
To use the AGA8 DETAIL and GERG equiations of state you typically
start by calling the new() and setup() functions to initialize
the calculations. Then you set the gas composition x, the pressure p
and the temperature t. Lastly you call the density() and
properties() functions to calculate the molar density and
the rest of the properies.
All of the calculation results are public fields in the struct that was
created with new().
Note that the gas composition is an array of excactly 21 components that must be in the order shown in the example.
```Rust use aga8::detail::Detail;
let mut aga8_test: Detail = Detail::new();
// Run seup() first to set up internal values aga8test.setup(); // Set the gas composition in mol fraction // The sum of all the components must be 1.0 aga8test.x = [ 0.778240, // Methane 0.020000, // Nitrogen 0.060000, // Carbon dioxide 0.080000, // Ethane 0.030000, // Propane 0.001500, // Isobutane 0.003000, // n-Butane 0.000500, // Isopentane 0.001650, // n-Pentane 0.002150, // Hexane 0.000880, // Heptane 0.000240, // Octane 0.000150, // Nonane 0.000090, // Decane 0.004000, // Hydrogen 0.005000, // Oxygen 0.002000, // Carbon monoxide 0.000100, // Water 0.002500, // Hydrogen sulfide 0.007000, // Helium 0.001000, // Argon ]; // Set pressure in kPa aga8test.p = 50000.0; // Set temperature in K aga8test.t = 400.0; // Run density to calculate the density in mol/l aga8test.density(); // Run properties to calculate all of the // output properties aga8test.properties();
// Molar density assert!((12.807 - aga8test.d).abs() < 1.0e-3); // Compressibility factor assert!((1.173 - aga8test.z).abs() < 1.0e-3); ```