Struct PanelMethod 

pub struct PanelMethod { /* private fields */ }

Source panel method for potential flow around a body.

Solves the boundary integral equation: -σ(x)/2 + ∫Γ ∂G/∂n(x,y) σ(y) dΓ(y) = -V∞ · n(x)

where σ is the source strength distribution, V_∞ is the free-stream velocity, and n is the outward normal. After solving, the surface velocity and pressure coefficient can be evaluated.

Implementations

impl PanelMethod

pub fn new(mesh: BoundaryMesh, cfg: PanelMethodConfig) -> Self

Create a new panel method solver.

pub fn solve(&mut self) -> IntegrateResult<()>

Solve for the source strength distribution σ on the boundary.

The no-penetration BC requires V_∞ · n = 0 on the body surface, leading to the RHS: b[i] = -V_∞ · n_i.

pub fn velocity_at(&self, p: [f64; 2]) -> [f64; 2]

Evaluate the induced velocity at a field point p (must be exterior).

u_induced(p) = ∫_Γ σ(y) ∇_p G(p, y) dΓ(y)

pub fn surface_cp(&self) -> IntegrateResult<(Vec<f64>, Vec<f64>)>

Evaluate the surface velocity magnitude and pressure coefficient on each panel’s midpoint.

Returns (velocities, cp_values) where each entry corresponds to a panel.

pub fn source_strengths(&self) -> &[f64]

Return a reference to the computed source strengths.

pub fn potential_at(&self, p: [f64; 2]) -> f64

Compute source panel velocity potential Φ at field point p.