scirs2-integrate 0.4.0

Numerical integration module for SciRS2 (scirs2-integrate)
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

//! Boundary condition types for fluid dynamics simulations.

/// Boundary condition types for fluid dynamics simulations
#[derive(Debug, Clone, Copy)]
pub enum FluidBoundaryCondition {
    /// No-slip (velocity = 0 at boundary)
    /// Used for solid wall boundaries where fluid adheres to the surface
    NoSlip,

    /// Free-slip (normal velocity = 0, tangential stress = 0)
    /// Used for slip walls or symmetry boundaries
    FreeSlip,

    /// Periodic boundary
    /// Used when the flow field repeats across the boundary
    Periodic,

    /// Inflow with specified velocity
    /// Takes (u, v) velocity components for 2D or (u, v) for the first two components in 3D
    Inflow(f64, f64),

    /// Outflow (zero gradient)
    /// Used for outlets where flow naturally exits the domain
    Outflow,
}

impl FluidBoundaryCondition {
    /// Check if the boundary condition is a wall type (NoSlip or FreeSlip)
    pub fn is_wall(&self) -> bool {
        matches!(
            self,
            FluidBoundaryCondition::NoSlip | FluidBoundaryCondition::FreeSlip
        )
    }

    /// Check if the boundary condition requires velocity specification
    pub fn requires_velocity(&self) -> bool {
        matches!(self, FluidBoundaryCondition::Inflow(_, _))
    }

    /// Get the inflow velocity if this is an inflow boundary condition
    pub fn get_inflow_velocity(&self) -> Option<(f64, f64)> {
        match self {
            FluidBoundaryCondition::Inflow(u, v) => Some((*u, *v)),
            _ => None,
        }
    }
}

impl Default for FluidBoundaryCondition {
    /// Default boundary condition is no-slip
    fn default() -> Self {
        FluidBoundaryCondition::NoSlip
    }
}