Struct AcousticSolver2D 

pub struct AcousticSolver2D {
    pub p: Grid2D,
    pub vx: Grid2D,
    pub vy: Grid2D,
    /* private fields */
}

Solveur acoustique 2D (FDTD, schéma saute-mouton sur grille décalée).

Disposition : p[x,y] au centre des cellules, vx[x,y] entre p[x,y] et p[x+1,y], vy[x,y] entre p[x,y] et p[x,y+1]. Les frontières sont des murs rigides (vitesse normale nulle), sauf si set_absorbing_boundary est utilisé.

Fields

p: Grid2D

vx: Grid2D

vy: Grid2D

Implementations

impl AcousticSolver2D

pub fn new( nx: usize, ny: usize, dx: f32, dy: f32, default_rho: f32, default_c: f32, ) -> Self

pub fn rho(&self) -> &Grid2D

Densité (lecture seule). Utiliser set_medium pour modifier.

pub fn set_pml_boundary(&mut self, width: usize, max_sigma: f32)

Configure des couches absorbantes PML sur les bords. width : épaisseur en nombre de cellules (ex: 20 à 30). max_sigma : intensité maximale de l’amortissement.

pub fn bulk_modulus(&self) -> &Grid2D

Module de compressibilité K = ρc² (lecture seule).

pub fn set_medium(&mut self, x: usize, y: usize, rho: f32, c: f32)

Définit le milieu d’une cellule : densité rho et vitesse du son c.

pub fn set_medium_rect( &mut self, x0: usize, y0: usize, x1: usize, y1: usize, rho: f32, c: f32, )

Définit le milieu d’un rectangle [x0, x1) × [y0, y1).

pub fn c_max(&self) -> f32

Vitesse du son maximale dans le domaine.

pub fn max_stable_dt(&self) -> f32

Pas de temps maximal stable (condition CFL 2D) pour le milieu actuel.

pub fn stable_dt(&self, safety: f32) -> f32

safety × max_stable_dt(), avec 0 < safety <= 1 (0.9 est un bon choix).

pub fn set_absorbing_boundary(&mut self, width: usize, max_damping: f32)

Couche éponge de width cellules sur les 4 bords.

max_damping (entre 0 et 1) est la fraction d’amplitude retirée à chaque pas sur le bord extérieur ; le profil décroît de façon quadratique vers l’intérieur. Valeurs typiques : width 15–30, max_damping 0.05–0.2. width = 0 désactive l’amortissement (murs rigides).

pub fn add_source( &mut self, x: usize, y: usize, waveform: Waveform, amplitude: f32, )

Ajoute une source ponctuelle pilotée par une forme d’onde.

pub fn inject_source(&mut self, x: usize, y: usize, value: f32)

Injecte instantanément une impulsion de pression.

pub fn add_receiver(&mut self, x: usize, y: usize) -> usize

Ajoute un récepteur ; retourne son identifiant.

pub fn receiver(&self, id: usize) -> &[f32]

Signal de pression enregistré (un échantillon par pas de temps).

pub fn time(&self) -> f64

Temps simulé écoulé (secondes).

pub fn step(&mut self, dt: f32)

Effectue un pas de simulation temporel dt.