Source code for aerosandbox.geometry.nosecone_shapes.haack

import aerosandbox.numpy as np


[docs]def haack_series(
        x_over_L: np.ndarray,
        C=1 / 3
):
    theta = np.arccos(1 - 2 * x_over_L)
    radius = ((
                      theta - np.sin(2 * theta) / 2 + C * np.sin(theta) ** 3
              ) / np.pi) ** 0.5
    return radius


[docs]def karman(
        x_over_L: np.ndarray,
):
    return haack_series(
        x_over_L=x_over_L,
        C=0
    )


[docs]def LV_haack(
        x_over_L: np.ndarray,
):
    return haack_series(
        x_over_L=x_over_L,
        C=1 / 3
    )


[docs]def tangent(
        x_over_L: np.ndarray,
):
    return haack_series(
        x_over_L=x_over_L,
        C=2 / 3
    )


if __name__ == '__main__':
    import matplotlib.pyplot as plt
    import aerosandbox.tools.pretty_plots as p

    fig, ax = plt.subplots()

[docs]    x_over_L = np.cosspace(0, 1, 2000)
    FR = 3

    plt.plot(x_over_L * FR, karman(x_over_L), label="$C=0$ (LD-Haack, Karman)")
    plt.plot(x_over_L * FR, LV_haack(x_over_L), label="$C=1/3$ (LV-Haack)")
    plt.plot(x_over_L * FR, tangent(x_over_L), label="$C=2/3$ (Tangent)")

    p.equal()
    p.show_plot(
        f"Nosecone Haack Series\nFineness Ratio $FR = {FR}$",
        "Length $x$",
        "Radius $r$"
    )