stormbird 0.9.0

A library for modeling modern wind propulsion devices
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

// Copyright (C) 2024, NTNU 
// Author: Jarle Vinje Kramer <jarlekramer@gmail.com; jarle.a.kramer@ntnu.no>
// License: GPL v3.0 (see separate file LICENSE or https://www.gnu.org/licenses/gpl-3.0.html)


use stormath::spatial_vector::{
    SpatialVector,
    transformations::RotationType,
    geometry_functions,
};

use stormath::type_aliases::Float;

use super::super::panel::Panel;

#[test]
/// Tests that checks that the near and far field calculations are similar for a panel with unit 
/// strength.
fn compare_near_and_far_field() {
    // Set up a panel in a local coordinate system
    let panel_points_0 = [
        SpatialVector::from([0.0, 0.0, 0.0]),
        SpatialVector::from([1.0, 0.0, 0.0]),
        SpatialVector::from([1.0, 1.0, 0.0]),
        SpatialVector::from([0.0, 1.0, 0.0]),
    ];

    // Set a control point in the local coordinate system
    let ctrl_point_0 = SpatialVector::new(5.0, 2.0, 1.0);

    // Add some rotation to test a "arbitrary" oriented panel (the angles should not matter)
    let rotation = SpatialVector::new(
        Float::from(45.0).to_radians(),
        Float::from(-10.0).to_radians(),
        0.0,
    );

    let rotation_type = RotationType::XYZ;

    let panel_points = [
        panel_points_0[0].rotate(rotation, rotation_type),
        panel_points_0[1].rotate(rotation, rotation_type),
        panel_points_0[2].rotate(rotation, rotation_type),
        panel_points_0[3].rotate(rotation, rotation_type),
    ];

    let ctrl_point = ctrl_point_0.rotate(rotation, rotation_type);

    let panel_area = geometry_functions::area_of_quadrilateral(&panel_points);

    let panel = Panel::new(
        panel_points,
        3.0,
        0.1,
    );

    assert!(panel_area - 1.0 < 1e-10);

    // Compare velocity calculations
    let u_i_near = panel.induced_velocity_as_vortex_lines_with_unit_strength(ctrl_point);
    let u_i_far  = panel.induced_velocity_as_point_doublet_with_unit_strength(ctrl_point);

    let difference = (u_i_near - u_i_far).length();

    dbg!(u_i_near);
    dbg!(u_i_far);
    dbg!(panel_area);

    let relative_error = difference / u_i_near.length();

    dbg!(relative_error);

    assert!(relative_error < 0.02);
}