example_fuel_qar_data.ipynb

Kernel: .venv

import pandas as pd
import openap
import matplotlib.pyplot as plt

pd.options.display.max_columns = 100

flight = pd.read_csv("data/flight_a320_qar.csv").query("ALTI_STD_FT>500")
flight.head()

fuelflow = openap.FuelFlow(ac="a320", engine="CFM56-5B4/P")

acc = (
    flight.TRUE_AIR_SPD_KT.diff().bfill()
    * openap.aero.kts
    / flight.FLIGHT_TIME.diff().bfill()
)

ff_estimate = fuelflow.enroute(
    mass=flight.MASS_KG,
    tas=flight.TRUE_AIR_SPD_KT,
    alt=flight.ALTI_STD_FT,
    vs=flight.VERT_SPD_FTMN,
)

plt.plot(
    flight.FLIGHT_TIME,
    flight.FUEL_FLOW_KGH * 2 / 3600,
    lw=1,
    label="actual (kg/s",
)
plt.plot(
    flight.FLIGHT_TIME,
    ff_estimate,
    alpha=0.8,
    lw=1,
    label="estimate (kg/s)",
)
plt.legend()

<matplotlib.legend.Legend at 0x70397d796720>