@Test
public void compresibility() {
double pressure = 101325;
double temperature = 300;
double expResult = 0.969775;
vdw.setTemperature(temperature);
vdw.setPressure(pressure);
double result = vdw.calculateCompresibilityFactor();
assertEquals(expResult, result, 1e-3);
}
public MultiComponentTest() {
ethane = new Compound("ethane");
ethane.setAcentricFactor(0.09781);
ethane.setCriticalTemperature(305.43);
ethane.setCriticalPressure(48.1595 * 101325);
ethane.setK_StryjekAndVera(0.02669);
propane = new Compound("propane");
propane.setAcentricFactor(0.15416);
propane.setCriticalTemperature(369.82);
propane.setCriticalPressure(41.9396 * 101325);
propane.setK_StryjekAndVera(0.03136);
nHeptane = new Compound("n-heptane");
nHeptane.setAcentricFactor(0.35);
nHeptane.setCriticalTemperature(540.14);
nHeptane.setCriticalPressure(27.05 * 101325);
nHeptane.setK_StryjekAndVera(0.02325);
Cubic eos = EquationsOfState.pengRobinson();
Alpha alpha = Alphas.getStryjekAndVeraExpression();
InteractionParameter k = new InteractionParameter(true);
k.setValue(ethane, propane, 0.2);
k.setValue(ethane, nHeptane, 0.1);
k.setValue(propane, nHeptane, 0.05);
HashSet<Compound> components = new HashSet();
components.add(ethane);
components.add(propane);
components.add(nHeptane);
MixingRule mr = new VDWMixingRule();
vdw = new Mixture(eos, alpha, components, Phase.VAPOR, mr, k);
vdw.setFraction(ethane, 0.3);
vdw.setFraction(propane, 0.3);
vdw.setFraction(nHeptane, 0.4);
}