def test_mixture_viscosity_lgas():
test_mix_viscos = pd.read_csv(os.path.join(pp_dir, 'properties', 'lgas',
'viscosity.txt'),
header=None,
sep=' ').values[:, 1]
viscos_ch4 = pd.read_csv(os.path.join(pp_dir, 'properties', 'methane',
'viscosity.txt'),
header=None,
sep=' ').values[:, 1]
viscos_n2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'nitrogen',
'viscosity.txt'),
header=None,
sep=' ').values[:, 1]
viscos_co2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'carbondioxid',
'viscosity.txt'),
header=None,
sep=' ').values[:, 1]
visco_c2h6 = pd.read_csv(os.path.join(pp_dir, 'properties', 'ethane',
'viscosity.txt'),
header=None,
sep=' ').values[:, 1]
molar_ch4 = pd.read_csv(os.path.join(pp_dir, 'properties', 'methane',
'molar_mass.txt'),
header=None,
sep=' ').values[0]
molar_n2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'nitrogen',
'molar_mass.txt'),
header=None,
sep=' ').values[0]
molar_co2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'carbondioxid',
'molar_mass.txt'),
header=None,
sep=' ').values[0]
molar_c2h6 = pd.read_csv(os.path.join(pp_dir, 'properties', 'ethane',
'molar_mass.txt'),
header=None,
sep=' ').values[0]
components_viscosities = np.concatenate([[viscos_ch4], [viscos_n2],
[viscos_co2], [visco_c2h6]])
components_molar_proportions = np.array([.85, .103, .013, .034])
components_molar_mass = np.concatenate(
[molar_ch4, molar_n2, molar_co2, molar_c2h6])
mix_viscos = calculate_mixture_viscosity(components_viscosities,
components_molar_proportions,
components_molar_mass)
assert np.allclose(mix_viscos, test_mix_viscos, rtol=0, atol=1e-8)
mix_viscos = calculate_mixture_viscosity(components_viscosities[:, 0],
components_molar_proportions,
components_molar_mass)
assert np.allclose(mix_viscos, test_mix_viscos, rtol=0, atol=2e-5)
def test_mixture_viscosity_lgas():
test_mix_viscos = pd.read_csv(os.path.join(pp_dir, 'properties', 'lgas', 'viscosity.txt'),
header=None, sep=' ', comment='#').values[:, 1]
viscos_ch4 = pd.read_csv(os.path.join(pp_dir, 'properties', 'methane', 'viscosity.txt'),
header=None, sep=' ', comment='#').values
viscos_n2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'nitrogen', 'viscosity.txt'),
header=None, sep=' ', comment='#').values
viscos_co2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'carbondioxide', 'viscosity.txt'),
header=None, sep=' ', comment='#').values
viscos_c2h6 = pd.read_csv(os.path.join(pp_dir, 'properties', 'ethane', 'viscosity.txt'),
header=None, sep=' ', comment='#').values
molar_ch4 = pd.read_csv(os.path.join(pp_dir, 'properties', 'methane', 'molar_mass.txt'),
header=None, sep=' ', comment='#').values[0]
molar_n2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'nitrogen', 'molar_mass.txt'),
header=None, sep=' ', comment='#').values[0]
molar_co2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'carbondioxide', 'molar_mass.txt'),
header=None, sep=' ', comment='#').values[0]
molar_c2h6 = pd.read_csv(os.path.join(pp_dir, 'properties', 'ethane', 'molar_mass.txt'),
header=None, sep=' ', comment='#').values[0]
values_to_take = [263, 273, 283, 293, 303, 313, 323, 333]
viscos_ch4 = interp1d(viscos_ch4[:, 0], viscos_ch4[:, 1])
viscos_ch4 = viscos_ch4((values_to_take))
viscos_n2 = interp1d(viscos_n2[:, 0], viscos_n2[:, 1], fill_value='extrapolate')
viscos_n2 = viscos_n2((values_to_take))
viscos_co2 = interp1d(viscos_co2[:, 0], viscos_co2[:, 1])
viscos_co2 = viscos_co2((values_to_take))
viscos_c2h6 = interp1d(viscos_c2h6[:, 0], viscos_c2h6[:, 1])
viscos_c2h6 = viscos_c2h6((values_to_take))
components_viscosities = np.concatenate([[viscos_ch4], [viscos_n2], [viscos_co2], [viscos_c2h6]])
components_molar_proportions = np.array([.85, .103, .013, .034])
components_molar_mass = np.concatenate([molar_ch4, molar_n2, molar_co2, molar_c2h6])
mix_viscos = calculate_mixture_viscosity(
components_viscosities, components_molar_proportions, components_molar_mass)
assert np.allclose(mix_viscos, test_mix_viscos, rtol=0, atol=1e-8)
mix_viscos = calculate_mixture_viscosity(
components_viscosities[:, 0], components_molar_proportions, components_molar_mass)
assert np.allclose(mix_viscos, test_mix_viscos, rtol=0, atol=2e-5)