kalast

Thermophysical model for binary systems of asteroids

Demo

In action

Calculations of the temperature at the surface of the secondary body of a binary system of asteroids:

use kalast::{BinarySystem, Properties, Time, World, ASTRONAUMICAL_UNIT, HOUR, MINUTE, V3, YEAR};

fn main() {
    let binary_system = BinarySystem::new(
        "Two spheres",
        V3::new(1.0, 0.0, 0.0) * ASTRONAUMICAL_UNIT * 1.664,
    )
    .with_primary(
        "rsc/obj/sphere_light.obj",
        Properties::new(
            1.0 * YEAR, // rotation period
            2.7 * HOUR, // revolution period
            0.0,        // obliquity
            100.0,      // thermal inertia
            2146.0,     // density
            600.0,      // heat capacity
            0.1,        // albedo
            0.9,        // emissivity
        ),
    )
    .with_secondary(
        "rsc/obj/ellipsoid_light.obj",
        Properties::new(
            11.92 * HOUR, // rotation period
            11.92 * HOUR, // revolution period
            0.0,          // obliquity 162.0 * TAU / 360.
            500.0,        // thermal inertia
            2146.0,       // density
            600.0,        // heat capacity
            0.07,         // albedo
            0.9,          // emissivity
        ),
        V3::new(-1.0, 0.0, 0.0) * 1e3, // relative position
    );

    binary_system.complete_model(true);
    binary_system.primary_mut().fixed(true);
    binary_system.secondary_mut().fixed(true);

    // Time
    let time = Time::new(
        binary_system.secondary().properties.revolution_period() * 50.,
        30.0 * MINUTE,
    );

    // World
    let mut world = World::new(time, binary_system);
    world.start();
    world.save("rsc/data/tmp.txt");

    println!(
        "{} {}",
        world.environment_system.secondary().temperatures().min(),
        world.environment_system.secondary().temperatures().max()
    );
}

Build source

# build as a release
./compile.sh -r

# run
./kalast

TODO

• finit diff, FEM, 2D/3D convection
• occultations & shadowing

License