Simulation de dispersion balistique et gazeuse

z0 altitude explosion (m) · statusCb callback progression UI

Récupère météo. Initialise 180 points sur cercle 1 km autour du point d'explosion avec vitesse Ariane(z0). Pour chaque niveau météo (ordre décroissant) : 180 pts × 3 grains × 3 facteurs vent = jusqu'à 1620 tâches distribuées aux Workers. Coque convexe entre couches. Rend après chaque couche.

Objet res : enveloppes par couche, profil, position explosion, aire sol (m²). Temps calcul, t_first, t_last impacts.

Tableau {sIn, zTarget, Cd0, cdType, m, ventF} · profil météo complet · latitude site (rad)

Découpe en min(hardwareConcurrency−1, 8) lots équilibrés. postMessage par Worker. Promise.all().

Tableau d'états [x, y, z, vx, vy, vz, t] dans l'ordre original des tâches.

État initial [x, y, z, vx, vy, vz, t] · altitude cible · paramètres grain et vent

Boucle RK45 adaptatif. Critère Hairer-Wanner : ATOL=[50m, 50m, 50m, 1m/s, 1m/s, 1m/s, 1s], RTOL=1e-3. DT_MAX=5s. Pas borné par (zLeft/|vz|) pour ne pas dépasser la cible. Interpolation linéaire à l'atterrissage exact sur zTarget. Sécurité 100 000 itérations.

État [x, y, z, vx, vy, vz, t] exactement à zTarget.

État courant s · profil vent · paramètres grain et perturbation vent

ISA → ρ(z), a(z). Vent interpolé × ventF. v_rel = v_vent − v_grain (3D). Ma = |v_rel|/a. Cd(Ma) table selon type. k = ½·ρ·Cd·A/m. Coriolis : Ω décomposé en ENU (Ω_N=Ω·cos(φ), Ω_U=Ω·sin(φ)). Gravité locale g(z) = g₀·(Rₑ/(Rₑ+z))².

[vx, vy, vz, ax, ay, az, 1] — le composant 1 cumule le temps sans coût.

z altitude (m), borné à 0 si négatif

Recherche couche (7 couches OACI 0–86 km). T(z) = T_base + dT·(z−z_base). Pression : isotherme p=p_b·exp(−g·dz/R·T) ou gradient p=p_b·(T/T_b)^(−g/R·dT). Correction humidité : T_v = T·(1+0.378·e/p), HR(z) = HR_SOL·exp(−z/2000).

{rho, a} — densité kg/m³ (air humide) · vitesse son m/s (air sec, convention standard).

Tableau [[x,y],...]

Tri lexicographique O(n log n). Balayage bas (gauche→droite) puis haut (droite→gauche). Élimination concavités par produit vectoriel cross(O,A,B)≤0. Dépile si nécessaire.

Indices des sommets de la coque convexe dans l'ordre. Réduit les N_PTS × 9 combos aux seuls extrêmes pour la couche suivante — complexité bornée.

h altitude (m)

ariane : interpolation linéaire table 11 points (0–50 km). deriveAriane : intégrale trapèzes de tan(pitch)·dh → dérive horizontale dans la direction AZIMUT. tempsAriane : intégrale trapèzes de dh/vz → temps de vol depuis le sol.

{v, pitch} · dérive (m) · temps de vol (s) — utilisés pour position initiale des débris et fraction propergol consommée.

z0 altitude explosion (m) · profil vent open-meteo

Masse : Q = N_boosters × 142t × (1−t_vol/135s) × 18.8% × (1−0.25). Vent moyen u_moy sur couches 0→z0. σ_y(x)=sy_a·x^0.886, σ_z(x)=sz_a·x^sz_b — coefficients issus de PG_CLASSES[classe]. 6 classes A→F disponibles, coefficients Turner (1970). C_sol(t) = 2Q/((2π)^1.5·σ_y²·σ_z)·exp(−z0²/2σ_z²). Seuil opérationnel = AEGL-1/10 = 0.1 mg/m³. Pic par échantillonnage log 200 pts. t_touch et t_end par dichotomie 80 itérations. z0_seuil par dichotomie si noBreach.

Objet puff avec fonctions _Cmax(t) concentration sol · _Ccentre(t) concentration centre nuage (sans atténuation H) · _sy(t) σ_y. Timestamps t_touch, t_peak, t_end. r_aegl1 au pic. Si noBreach : z0_seuil.

Objet météo : u10 vent surface (m/s) · rad_sw radiation solaire (W/m²) · nuage nébulosité (%)

Table Pasquill (1961) / Turner (1970). Nuit si rad_sw < 10 W/m². Ciel couvert si nuage > 80% → classe D. Sinon croisement vent surface × insolation (forte >600 · modérée 300-600 · faible <300 W/m²). Retourne A, B, C, D, E ou F.

Lettre de classe PG. Utilisée si le select UI est sur AUTO. Sinon ignorée — la classe choisie manuellement prime.

z altitude (m) · profil niveaux météo triés par altitude

Recherche linéaire du segment encadrant z. Interpolation linéaire wx, wy. Extrapolation constante hors bornes.

{wx, wy} composantes vent ENU (m/s).

Mêmes paramètres que simuler()

Promise.race entre simuler() et un timeout de 3 minutes. Si le calcul dépasse 3 min (altitude très élevée, machine lente), rejette proprement la promesse et réactive le bouton. N'arrête pas les Workers — ils finissent leur lot courant.

Même retour que simuler(), ou erreur explicite "Calcul interrompu après 3 min".

forceRefresh booléen · coordonnées SITE globales

Court-circuité si PROFIL_VENT_STATIQUE≠null (mode test). Sinon : fetch JSON open-meteo, 16 niveaux hPa + surface_pressure, relative_humidity_2m, wind_speed_10m, shortwave_radiation, cloud_cover. Tranche horaire la plus proche. Conversion direction → ENU. Tri croissant par altitude. Sanity check : vent > 150 km/h → erreur explicite. Cache par nom de site.

{profil, heure, ptMesure, qnh, hr, u10, rad_sw, nuage} — profil de 1 à 16 niveaux + données surface temps réel.

État global : carte.lat/lon/zoom · carte.resultat · carte.puff

Fond. Imagerie Esri World Imagery + étiquettes Esri Boundaries and Places (cache Map borné à 600 tuiles), léger assombrissement. Repère pas de tir. Si resultat : renderEnveloppes() + renderBarbules(). Toujours : renderPuffIsoline().

Canvas 2D complet. Appelée après chaque couche balistique, pan/zoom, resize, et chaque calcul.

carte.resultat.enveloppes[] — polygones ENU par niveau météo

Pour chaque enveloppe : polyEnPix() (metToGeo + geoToPix). Remplissage couleurAlt(z) semi-transparent. Contour plein 2.4px (sol) ou pointillé 1px (couches). Label altitude au sommet (reduce min y). Sol : label "IMPACT SOL" au centroïde. Ligne trajectoire site→explosion. Marqueur explosion vert.

Enveloppes colorées du haut (bleu) vers le sol (rouge). Zone sol pleine avec label.

carte.resultat.profil[] — niveaux vent filtrés ≤ zMax+500m

Repère + croix au point de mesure open-meteo. Colonne verticale 38px/barbule ancrée en dessous. Chaque barbule : convention WMO (fanion △ = 50kt, barbe longue = 10kt, demi = 5kt, cercle = calme). Label altitude (km) + vitesse (m/s) à droite. Couleur = couleurAlt(z).

Colonne lisible, colorée par altitude, lue de haut en bas = bas de l'atmosphère vers le haut.

carte.puff — objet retourné par puffPG()

Si !noBreach : 40 instants entre t_touch et t_end. Pour chaque t : centre dérive (wx_moy·t, wy_moy·t), rayon r=σ_y·√(−2·ln(seuil/C_sol(t))). 24 points par cercle. Coque convexe → haricot rouge pointillé. Croix toutes les 15min de t=0 à t_end : position centre, label "+HhMM · C_centre(t)". Rouge si C_sol(t)>seuil, gris sinon.

Haricot rouge = zone exposée > AEGL-1/10. Croix = dérive du centre du nuage avec concentration maximale dans le nuage. Absent si noBreach.

Déplacement ENU depuis SITE : dxE Est (m) · dyN Nord (m)

d = √(dxE²+dyN²). Azimut = atan2(dxE, dyN). Distance angulaire ang = d/R (R=6371km). Projection directe Vincenty simplifié sphère : lat1 = asin(sin(lat0)·cos(ang)+cos(lat0)·sin(ang)·cos(brg)), lon1 = lon0 + atan2(...).

[lat, lon] degrés. Erreur <0.1% à 200 km (vs approximation plate qui accumule >1% à 100 km).

Contexte 2D · position (px) · direction d'où vient le vent (rad) · vitesse (m/s) · couleur CSS

Conversion m/s → kt, arrondi 5kt. Hampe L=28px dans la direction d'origine. Fanions triangulaires 50kt (remplis), barbes longues 10kt, demi-barbes 5kt. Espacement 6px. Cercle vide si <5kt.

Barbule météorologique standard WMO No.306 dessinée sur canvas.