openphys 0.1.3

OpenPhys is a simple physics engine library written in Rust.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98

use openphys::modules::heat_transfer::simulate_heat_transfer;
use openphys::utils::object::Object;

#[test]
fn test_heat_transfer_simulation() {
    let mut obj1 = Object {
        mass: 1.0,
        temperature: 373.15,
        specific_heat_capacity: 4186.0,
        ..Object::default()
    };
    let mut obj2 = Object {
        mass: 1.0,
        temperature: 273.15,
        specific_heat_capacity: 4186.0,
        ..Object::default()
    };
    let time_step = 0.1; // 0.1 seconds
    let equilibrium_threshold = 1e-6; // The equilibrium threshold is 10^-6 because
                                      // trying to check if two f64 values are equal
                                      // would cause an infinite (or close to
                                      // infinite) loop

    let result = simulate_heat_transfer(&mut obj1, &mut obj2, time_step, equilibrium_threshold);

    println!("Time steps: {}", result.time_steps);
    println!("Total time elapsed: {:.2}s", result.total_time);
    println!(
        "Total heat transferred: {:.2}K",
        result.total_heat_transferred
    );
    println!("Final temp of obj1: {:.2}K", result.final_temp_obj1);
    println!("Final temp of obj2: {:.2}K", result.final_temp_obj2);
}

#[test]
fn test2_heat_transfer_simulation() {
    let mut obj1 = Object {
        mass: 50.0,
        temperature: 480.15,
        specific_heat_capacity: 12049.0,
        ..Object::default()
    };
    let mut obj2 = Object {
        mass: 82.0,
        temperature: 582.62,
        specific_heat_capacity: 12049.0,
        ..Object::default()
    };
    let time_step = 0.1; // 0.1 seconds
    let equilibrium_threshold = 1e-6; // The equilibrium threshold is 10^-6 because
                                      // trying to check if two f64 values are equal
                                      // would cause an infinite (or close to
                                      // infinite) loop

    let result = simulate_heat_transfer(&mut obj1, &mut obj2, time_step, equilibrium_threshold);

    println!("Time steps: {}", result.time_steps);
    println!("Total time elapsed: {:.2}s", result.total_time);
    println!(
        "Total heat transferred: {:.2}K",
        result.total_heat_transferred
    );
    println!("Final temp of obj1: {:.2}K", result.final_temp_obj1);
    println!("Final temp of obj2: {:.2}K", result.final_temp_obj2);
}

#[test]
fn test3_heat_transfer_simulation_high_accuracy() {
    let mut obj1 = Object {
        mass: 50.0,
        temperature: 480.15,
        specific_heat_capacity: 12049.0,
        ..Object::default()
    };
    let mut obj2 = Object {
        mass: 82.0,
        temperature: 582.62,
        specific_heat_capacity: 12049.0,
        ..Object::default()
    };
    let time_step = 0.1; // 0.1 seconds
    let equilibrium_threshold = 1e-12; // The equilibrium threshold is 10^-12 because
                                       // trying to check if two f64 values are equal
                                       // would cause an infinite (or close to
                                       // infinite) loop

    let result = simulate_heat_transfer(&mut obj1, &mut obj2, time_step, equilibrium_threshold);

    println!("Time steps: {}", result.time_steps);
    println!("Total time elapsed: {:.2}s", result.total_time);
    println!(
        "Total heat transferred: {:.2}K",
        result.total_heat_transferred
    );
    println!("Final temp of obj1: {:.2}K", result.final_temp_obj1);
    println!("Final temp of obj2: {:.2}K", result.final_temp_obj2);
}