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); }