# Simulated ASTM Table 54 logic (Product Group B example) # Not official ASTM data – for understanding only def vcf_table54_groupB(temp_c, density_obs_kgm3): """ Simplified volume correction factor (VCF) to 15°C. Based on typical expansion coefficients for middle distillates. """ # Approx cubic expansion coefficient per °C for Group B # Real ASTM uses complex polynomials per density/temp grid alpha = 0.00085 # 1/°C
print(f"VCF (Table 54, Group B): vcf") print(f"Volume at 15°C: volume_at_15C:.2f m³")
# Density at 15°C estimate (simplified) density_15C = density_obs_kgm3 * (1 + alpha * (temp_c - 15))
# Volume Correction Factor vcf = density_15C / density_obs_kgm3 return round(vcf, 5) observed_density = 850.0 # kg/m³ at 30°C observed_temp = 30.0 # °C
vcf = vcf_table54_groupB(observed_temp, observed_density) observed_volume = 1000.0 # m³ volume_at_15C = observed_volume * vcf
