def test_mixture_molar_mass_lgas(): test_mix_molar = pd.read_csv(os.path.join(pp_dir, 'properties', 'lgas', 'molar_mass.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] components_molar_mass = np.concatenate([molar_ch4, molar_n2, molar_co2, molar_c2h6]) components_molar_proportions = np.array([.85, .103, .013, .034]) mix_molar = calculate_mixture_molar_mass( components_molar_mass=components_molar_mass, components_molar_proportions=components_molar_proportions) assert (mix_molar == test_mix_molar) components_mass_proportions = calculate_mass_fraction_from_molar_fraction( components_molar_proportions, components_molar_mass) mix_molar = calculate_mixture_molar_mass( components_molar_mass=components_molar_mass, components_mass_proportions=components_mass_proportions) assert (mix_molar == test_mix_molar)
def test_mixture_density_lgas(): test_mix_density = pd.read_csv(os.path.join(pp_dir, 'properties', 'lgas', 'density.txt'), header=None, sep=' ').values[:, 1] dens_ch4 = pd.read_csv(os.path.join(pp_dir, 'properties', 'methane', 'density.txt'), header=None, sep=' ').values[:, 1] dens_n2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'nitrogen', 'density.txt'), header=None, sep=' ').values[:, 1] dens_co2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'carbondioxid', 'density.txt'), header=None, sep=' ').values[:, 1] dens_c2h6 = pd.read_csv(os.path.join(pp_dir, 'properties', 'ethane', 'density.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_density = np.concatenate([[dens_ch4], [dens_n2], [dens_co2], [dens_c2h6]]) components_molar_proportions = np.array([.85, .103, .013, .034]) components_molar_mass = np.concatenate( [molar_ch4, molar_n2, molar_co2, molar_c2h6]) components_mass_proportions = calculate_mass_fraction_from_molar_fraction( components_molar_proportions, components_molar_mass) mix_dens = calculate_mixture_density( components_density, components_mass_proportions=components_mass_proportions) assert np.allclose(mix_dens, test_mix_density, rtol=0, atol=1e-10) mix_dens = calculate_mixture_density( components_density[:, 0], components_mass_proportions=components_mass_proportions) assert np.allclose(mix_dens, test_mix_density[0], rtol=0, atol=1e-10)
def test_mixture_heat_capacity_hgas(): test_mix_capacity = pd.read_csv(os.path.join(pp_dir, 'properties', 'hgas', 'heat_capacity.txt'), header=None, sep=' ').values[:, 1] cap_ch4 = pd.read_csv(os.path.join(pp_dir, 'properties', 'methane', 'heat_capacity.txt'), header=None, sep=' ').values[:, 1] cap_n2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'nitrogen', 'heat_capacity.txt'), header=None, sep=' ').values[:, 1] cap_co2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'carbondioxid', 'heat_capacity.txt'), header=None, sep=' ').values[:, 1] cap_c2h6 = pd.read_csv(os.path.join(pp_dir, 'properties', 'ethane', 'heat_capacity.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_capacity = np.concatenate([[cap_ch4], [cap_n2], [cap_co2], [cap_c2h6]]) components_molar_proportions = np.array([.964, .005, .005, .026]) components_molar_mass = np.concatenate( [molar_ch4, molar_n2, molar_co2, molar_c2h6]) components_mass_proportions = calculate_mass_fraction_from_molar_fraction( components_molar_proportions, components_molar_mass) mix_cap = calculate_mixture_heat_capacity(components_capacity, components_mass_proportions) assert np.allclose(mix_cap, test_mix_capacity, rtol=0, atol=1e-10)
def test_mixture_density_lgas(): test_mix_density = pd.read_csv(os.path.join(pp_dir, 'properties', 'lgas', 'density.txt'), header=None, sep=' ').values[:, 1] dens_ch4 = pd.read_csv(os.path.join(pp_dir, 'properties', 'methane', 'density.txt'), header=None, sep=' ', comment='#').values dens_n2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'nitrogen', 'density.txt'), header=None, sep=' ', comment='#').values dens_co2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'carbondioxide', 'density.txt'), header=None, sep=' ', comment='#').values dens_c2h6 = pd.read_csv(os.path.join(pp_dir, 'properties', 'ethane', 'density.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] dens_ch4 = interp1d(dens_ch4[:, 0], dens_ch4[:, 1]) dens_ch4 = dens_ch4((values_to_take)) dens_n2 = interp1d(dens_n2[:, 0], dens_n2[:, 1], fill_value='extrapolate') dens_n2 = dens_n2((values_to_take)) dens_co2 = interp1d(dens_co2[:, 0], dens_co2[:, 1]) dens_co2 = dens_co2((values_to_take)) dens_c2h6 = interp1d(dens_c2h6[:, 0], dens_c2h6[:, 1]) dens_c2h6 = dens_c2h6((values_to_take)) components_density = np.concatenate([[dens_ch4], [dens_n2], [dens_co2], [dens_c2h6]]) components_molar_proportions = np.array([.85, .103, .013, .034]) components_molar_mass = np.concatenate([molar_ch4, molar_n2, molar_co2, molar_c2h6]) components_mass_proportions = calculate_mass_fraction_from_molar_fraction( components_molar_proportions, components_molar_mass) mix_dens = calculate_mixture_density( components_density, components_mass_proportions=components_mass_proportions) assert np.allclose(mix_dens, test_mix_density, rtol=0, atol=1e-10) mix_dens = calculate_mixture_density( components_density[:, 0], components_mass_proportions=components_mass_proportions) assert np.allclose(mix_dens, test_mix_density[0], rtol=0, atol=1e-10)
def test_mixture_heat_capacity_hgas(): test_mix_capacity = pd.read_csv(os.path.join(pp_dir, 'properties', 'hgas', 'heat_capacity.txt'), header=None, sep=' ', comment='#').values[:, 1] cap_ch4 = pd.read_csv(os.path.join(pp_dir, 'properties', 'methane', 'heat_capacity.txt'), header=None, sep=' ', comment='#').values cap_n2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'nitrogen', 'heat_capacity.txt'), header=None, sep=' ', comment='#').values cap_co2 = pd.read_csv(os.path.join(pp_dir, 'properties', 'carbondioxide', 'heat_capacity.txt'), header=None, sep=' ', comment='#').values cap_c2h6 = pd.read_csv(os.path.join(pp_dir, 'properties', 'ethane', 'heat_capacity.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] cap_ch4 = interp1d(cap_ch4[:, 0], cap_ch4[:, 1]) cap_ch4 = cap_ch4((values_to_take)) cap_n2 = interp1d(cap_n2[:, 0], cap_n2[:, 1], fill_value='extrapolate') cap_n2 = cap_n2((values_to_take)) cap_co2 = interp1d(cap_co2[:, 0], cap_co2[:, 1]) cap_co2 = cap_co2((values_to_take)) cap_c2h6 = interp1d(cap_c2h6[:, 0], cap_c2h6[:, 1]) cap_c2h6 = cap_c2h6((values_to_take)) components_capacity = np.concatenate([[cap_ch4], [cap_n2], [cap_co2], [cap_c2h6]]) components_molar_proportions = np.array([.964, .005, .005, .026]) components_molar_mass = np.concatenate([molar_ch4, molar_n2, molar_co2, molar_c2h6]) components_mass_proportions = calculate_mass_fraction_from_molar_fraction( components_molar_proportions, components_molar_mass) mix_cap = calculate_mixture_heat_capacity(components_capacity, components_mass_proportions) assert np.allclose(mix_cap, test_mix_capacity, rtol=0, atol=1e-10)