@woosh/meep-engine 2.129.4 → 2.131.0

package/src/core/math/physics/mie/ri_air.d.ts

@@ -0,0 +1,7 @@
+/**
+ * ☁️ AIR (Medium)
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_air(result: number[], wavelength_nm: number): void;
+//# sourceMappingURL=ri_air.d.ts.map

package/src/core/math/physics/mie/ri_air.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ri_air.d.ts","sourceRoot":"","sources":["../../../../../../src/core/math/physics/mie/ri_air.js"],"names":[],"mappings":"AA4DA;;;;GAIG;AACH,+BAHW,MAAM,EAAE,iBACR,MAAM,QAKhB"}

package/src/core/math/physics/mie/ri_air.js

@@ -0,0 +1,69 @@
+import {
+    ParameterLookupTable
+} from "../../../../engine/graphics/particles/particular/engine/parameter/ParameterLookupTable.js";
+
+/**
+ * Values for standard air (15°C, 1 atm, 450 ppm CO₂)
+ */
+const AIR_REFRACTION_INDEX_LUT = ParameterLookupTable.from(
+    2, // number of columns in the table [n, k]
+    // RI lookup data (n, k) for standard air across UV to long-wave IR
+    [
+        // UV (short-wave)
+        1.0003280, 0.0, // 200nm
+        1.0003070, 0.0, // 250nm
+        1.0002930, 0.0, // 300nm
+        1.0002850, 0.0, // 350nm
+        // Visible baseline (kept from previous table for continuity)
+        1.0002801, 0.0, // 400nm
+        1.0002785, 0.0, // 450nm
+        1.0002773, 0.0, // 500nm
+        1.0002765, 0.0, // 550nm
+        1.0002759, 0.0, // 600nm
+        1.0002755, 0.0, // 650nm
+        1.0002751, 0.0, // 700nm
+        // Near-IR to long-wave IR
+        1.0002747, 0.0, // 800nm
+        1.0002740, 0.0, // 900nm
+        1.0002734, 0.0, // 1000nm (1µm)
+        1.0002720, 0.0, // 1500nm (1.5µm)
+        1.0002710, 0.0, // 2000nm (2µm)
+        1.0002700, 0.0, // 3000nm (3µm)
+        1.0002690, 0.0, // 5000nm (5µm)
+        1.0002685, 0.0, // 10000nm (10µm)
+        1.0002682, 0.0  // 20000nm (20µm)
+    ],
+    // lookup positions, in nanometers
+    [
+        200,
+        250,
+        300,
+        350,
+        400,
+        450,
+        500,
+        550,
+        600,
+        650,
+        700,
+        800,
+        900,
+        1000,
+        1500,
+        2000,
+        3000,
+        5000,
+        10000,
+        20000
+    ]
+)
+
+/**
+ * ☁️ AIR (Medium)
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_air(result, wavelength_nm) {
+    // can effectively be assumed [1,0] but we're being more precise
+    AIR_REFRACTION_INDEX_LUT.sample(wavelength_nm, result);
+}

package/src/core/math/physics/mie/ri_ammonium_sulfate.d.ts

@@ -0,0 +1,7 @@
+/**
+ * 🌫️ Ammonium sulfate (continental haze surrogate)
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_ammonium_sulfate(result: number[], wavelength_nm: number): void;
+//# sourceMappingURL=ri_ammonium_sulfate.d.ts.map

package/src/core/math/physics/mie/ri_ammonium_sulfate.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ri_ammonium_sulfate.d.ts","sourceRoot":"","sources":["../../../../../../src/core/math/physics/mie/ri_ammonium_sulfate.js"],"names":[],"mappings":"AAoIA;;;;GAIG;AACH,4CAHW,MAAM,EAAE,iBACR,MAAM,QAIhB"}

package/src/core/math/physics/mie/ri_ammonium_sulfate.js

@@ -0,0 +1,140 @@
+import {
+    ParameterLookupTable
+} from "../../../../engine/graphics/particles/particular/engine/parameter/ParameterLookupTable.js";
+
+/**
+ * Ammonium sulfate aerosol ("(NH4)2SO4") complex refractive index m = n + i k, wavelength in nm.
+ * Assumptions: dry crystalline aerosol surrogate for continental haze; ~300 K; morphology not modeled.
+ * Guidance/sources (composite, smoothed for Mie with linear interpolation):
+ * - Toon, Pollack & Khare (1976); d'Almeida et al. (1991); Hess et al. (1998, OPAC) — visible/NIR anchors (n≈1.52, k≈0).
+ * - Palmer & Williams (1975) optical constants compendium; Querry (1985/1991) — UV and MIR/LWIR band behavior.
+ * Notes:
+ * - In the visible, k is extremely small (≈1e-6–1e-5). We use small nonzero values to avoid divide-by-zero corner cases.
+ * - MIR bands associated with sulfate and ammonium groups raise k around ~3 µm (weak) and ~7–10 µm (stronger).
+ * - Coverage: 200–20,000 nm; denser sampling near MIR bands.
+ */
+const AMMONIUM_SULFATE_REFRACTION_INDEX_LUT = ParameterLookupTable.from(
+    2, // [n, k]
+    [
+        // UV (short-wave) — increasing absorption
+        1.560, 2.0e-2, // 200 nm
+        1.550, 8.0e-3, // 220 nm
+        1.545, 3.5e-3, // 240 nm
+        1.540, 1.5e-3, // 260 nm
+        1.535, 8.0e-4, // 280 nm
+        1.532, 3.5e-4, // 300 nm
+        1.530, 2.0e-4, // 320 nm
+        1.528, 1.2e-4, // 340 nm
+        1.5265, 8.0e-5, // 360 nm
+        1.5255, 6.0e-5, // 380 nm
+        // Visible (very low absorption; n≈1.52)
+        1.5250, 5.0e-5, // 400 nm
+        1.5245, 3.0e-5, // 450 nm
+        1.5240, 2.0e-5, // 500 nm
+        1.5235, 2.0e-5, // 550 nm
+        1.5230, 2.0e-5, // 600 nm
+        1.5225, 2.0e-5, // 650 nm
+        1.5220, 2.0e-5, // 700 nm
+        // NIR/SWIR (0.7–2.5 µm): still tiny k; slight dispersion in n
+        1.5218, 2.0e-5, // 750 nm
+        1.5215, 2.0e-5, // 800 nm
+        1.5210, 3.0e-5, // 900 nm
+        1.5205, 4.0e-5, // 1000 nm
+        1.5200, 4.5e-5, // 1064 nm
+        1.5195, 6.0e-5, // 1200 nm
+        1.5190, 8.0e-5, // 1300 nm
+        1.5185, 1.6e-4, // 1450 nm
+        1.5180, 1.3e-4, // 1550 nm
+        1.5178, 1.2e-4, // 1650 nm
+        1.5175, 1.2e-4, // 1750 nm
+        1.5172, 1.8e-4, // 1940 nm
+        1.5170, 1.8e-4, // 2000 nm
+        1.5165, 2.2e-4, // 2200 nm
+        1.5160, 2.8e-4, // 2400 nm
+        1.5158, 3.0e-4, // 2500 nm
+        // Mid-IR (2.7–5 µm): weak-to-moderate absorption near 3 µm
+        1.5155, 6.0e-4, // 2700 nm
+        1.5150, 1.5e-3, // 2900 nm
+        1.5145, 3.5e-3, // 3000 nm
+        1.5140, 5.0e-3, // 3100 nm
+        1.5135, 6.0e-3, // 3200 nm
+        1.5130, 6.5e-3, // 3300 nm
+        1.5125, 6.0e-3, // 3500 nm
+        1.5120, 4.0e-3, // 3700 nm
+        1.5115, 2.5e-3, // 4000 nm
+        1.5110, 1.6e-3, // 4500 nm
+        1.5105, 1.2e-3, // 5000 nm
+        // Long-wave IR (5–20 µm): strong sulfate bands ~7–10 µm
+        1.5100, 4.0e-3, // 6000 nm
+        1.5100, 8.0e-3, // 7000 nm
+        1.5105, 1.8e-2, // 8000 nm
+        1.5110, 3.5e-2, // 9000 nm
+        1.5120, 5.5e-2, // 10000 nm (10 µm)
+        1.5140, 3.5e-2, // 12000 nm
+        1.5160, 1.5e-2, // 15000 nm
+        1.5180, 8.0e-3, // 20000 nm
+    ],
+    [
+        200,
+        220,
+        240,
+        260,
+        280,
+        300,
+        320,
+        340,
+        360,
+        380,
+        400,
+        450,
+        500,
+        550,
+        600,
+        650,
+        700,
+        750,
+        800,
+        900,
+        1000,
+        1064,
+        1200,
+        1300,
+        1450,
+        1550,
+        1650,
+        1750,
+        1940,
+        2000,
+        2200,
+        2400,
+        2500,
+        2700,
+        2900,
+        3000,
+        3100,
+        3200,
+        3300,
+        3500,
+        3700,
+        4000,
+        4500,
+        5000,
+        6000,
+        7000,
+        8000,
+        9000,
+        10000,
+        12000,
+        15000,
+        20000
+    ]
+);
+
+/**
+ * 🌫️ Ammonium sulfate (continental haze surrogate)
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_ammonium_sulfate(result, wavelength_nm) {
+    AMMONIUM_SULFATE_REFRACTION_INDEX_LUT.sample(wavelength_nm, result);
+}

package/src/core/math/physics/mie/ri_brine.d.ts

@@ -0,0 +1,7 @@
+/**
+ * 🌊 Sea-salt (brine) — maritime haze surrogate
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_brine(result: number[], wavelength_nm: number): void;
+//# sourceMappingURL=ri_brine.d.ts.map

package/src/core/math/physics/mie/ri_brine.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ri_brine.d.ts","sourceRoot":"","sources":["../../../../../../src/core/math/physics/mie/ri_brine.js"],"names":[],"mappings":"AAyIA;;;;GAIG;AACH,iCAHW,MAAM,EAAE,iBACR,MAAM,QAIhB"}

package/src/core/math/physics/mie/ri_brine.js

@@ -0,0 +1,145 @@
+import {
+    ParameterLookupTable
+} from "../../../../engine/graphics/particles/particular/engine/parameter/ParameterLookupTable.js";
+
+/**
+ * Sea-salt aerosol (deliquesced brine) complex refractive index m = n + i k, wavelength in nm.
+ * Assumptions: maritime haze surrogate; aqueous brine droplet at ~300 K; morphology not modeled.
+ * Guidance/sources (composite, smoothed):
+ * - Refractive index of seawater/brine in visible n≈1.37–1.40 (dependent on salinity and wavelength).
+ * - Hale & Querry (1973) water bands for k behavior in IR; salt slightly perturbs n, k remains water-like.
+ * - Querry NaCl for dry-salt context; here we target deliquesced brine typical for haze.
+ * Notes:
+ * - Visible k is extremely small; we use tiny nonzero values for numerical stability.
+ * - MIR/LWIR follow water band structure but with slightly different n baseline.
+ * - Coverage: 200–20,000 nm; denser near 3 µm and 6–10 µm regions.
+ */
+const SEASALT_BRINE_REFRACTION_INDEX_LUT = ParameterLookupTable.from(
+    2, // [n, k]
+    [
+        // UV (short-wave)
+        1.4100, 2.0e-4, // 200 nm (strong UV absorption relative to visible)
+        1.4000, 6.0e-5, // 220 nm
+        1.3950, 2.0e-5, // 240 nm
+        1.3900, 8.0e-6, // 260 nm
+        1.3860, 3.0e-6, // 280 nm
+        1.3820, 1.2e-6, // 300 nm
+        1.3795, 8.0e-7, // 320 nm
+        1.3775, 5.0e-7, // 340 nm
+        1.3760, 3.0e-7, // 360 nm
+        1.3750, 2.0e-7, // 380 nm
+        // Visible (400–700 nm): n≈1.38–1.39, very small k
+        1.3745, 1.5e-7, // 400 nm
+        1.3738, 1.2e-7, // 450 nm
+        1.3730, 1.0e-7, // 500 nm
+        1.3725, 1.0e-7, // 550 nm
+        1.3720, 1.2e-7, // 600 nm
+        1.3718, 1.5e-7, // 650 nm
+        1.3715, 2.0e-7, // 700 nm
+        // NIR/SWIR (0.7–2.5 µm): weak absorption until ~970/1450/1940 nm bands
+        1.3712, 3.0e-7, // 750 nm
+        1.3710, 4.0e-7, // 800 nm
+        1.3705, 1.2e-6, // 900 nm
+        1.3700, 3.5e-6, // 950 nm
+        1.3695, 1.5e-5, // 970 nm
+        1.3690, 8.0e-6, // 1000 nm
+        1.3685, 3.0e-6, // 1064 nm
+        1.3678, 4.5e-5, // 1200 nm
+        1.3672, 1.8e-4, // 1300 nm
+        1.3660, 1.2e-3, // 1450 nm
+        1.3650, 2.0e-4, // 1550 nm
+        1.3640, 4.5e-4, // 1650 nm
+        1.3630, 1.8e-3, // 1750 nm
+        1.3610, 1.2e-2, // 1940 nm
+        1.3600, 6.0e-3, // 2000 nm
+        1.3580, 2.0e-3, // 2200 nm
+        1.3560, 9.0e-4, // 2400 nm
+        1.3550, 6.0e-4, // 2500 nm
+        // Mid-IR (2.5–5 µm) — OH stretch complex 2.9–3.2 µm
+        1.3520, 1.0e-2, // 2700 nm
+        1.3480, 8.0e-2, // 2900 nm
+        1.3460, 3.5e-1, // 3000 nm
+        1.3430, 9.0e-1, // 3100 nm
+        1.3400, 1.5e+0, // 3200 nm
+        1.3380, 1.2e+0, // 3300 nm
+        1.3350, 7.0e-1, // 3500 nm
+        1.3320, 3.5e-1, // 3700 nm
+        1.3290, 1.8e-1, // 4000 nm
+        1.3260, 9.0e-2, // 4500 nm
+        1.3240, 6.0e-2, // 5000 nm
+        // Long-wave IR (5–20 µm): elevated k due to vibrational modes
+        1.3200, 1.8e-1, // 6000 nm
+        1.3180, 3.5e-1, // 7000 nm
+        1.3160, 7.0e-1, // 8000 nm
+        1.3140, 9.0e-1, // 9000 nm
+        1.3125, 1.1e+0, // 10000 nm
+        1.3110, 9.0e-1, // 12000 nm
+        1.3095, 7.0e-1, // 15000 nm
+        1.3080, 5.5e-1, // 20000 nm
+    ],
+    [
+        200,
+        220,
+        240,
+        260,
+        280,
+        300,
+        320,
+        340,
+        360,
+        380,
+        400,
+        450,
+        500,
+        550,
+        600,
+        650,
+        700,
+        750,
+        800,
+        900,
+        950,
+        970,
+        1000,
+        1064,
+        1200,
+        1300,
+        1450,
+        1550,
+        1650,
+        1750,
+        1940,
+        2000,
+        2200,
+        2400,
+        2500,
+        2700,
+        2900,
+        3000,
+        3100,
+        3200,
+        3300,
+        3500,
+        3700,
+        4000,
+        4500,
+        5000,
+        6000,
+        7000,
+        8000,
+        9000,
+        10000,
+        12000,
+        15000,
+        20000
+    ]
+);
+
+/**
+ * 🌊 Sea-salt (brine) — maritime haze surrogate
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_brine(result, wavelength_nm) {
+    SEASALT_BRINE_REFRACTION_INDEX_LUT.sample(wavelength_nm, result);
+}

package/src/core/math/physics/mie/ri_dust.d.ts

@@ -0,0 +1,7 @@
+/**
+ *  💨 MINERAL DUST
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_dust(result: number[], wavelength_nm: number): void;
+//# sourceMappingURL=ri_dust.d.ts.map

package/src/core/math/physics/mie/ri_dust.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ri_dust.d.ts","sourceRoot":"","sources":["../../../../../../src/core/math/physics/mie/ri_dust.js"],"names":[],"mappings":"AAiJA;;;;GAIG;AACH,gCAHW,MAAM,EAAE,iBACR,MAAM,QAIhB"}

package/src/core/math/physics/mie/ri_dust.js

@@ -0,0 +1,153 @@
+import {
+    ParameterLookupTable
+} from "../../../../engine/graphics/particles/particular/engine/parameter/ParameterLookupTable.js";
+
+/**
+ * Complex refractive index of mineral dust (silicate-rich, e.g., Saharan): n + i k, wavelength in nm.
+ * Assumptions: dry, ~300 K, bulk effective constants representative of a typical desert mineral dust
+ * (quartz/feldspar/clays mix). Morphology not modeled; values intended for bulk-equivalent Mie.
+ *
+ * Composite sources and notes:
+ * - OPAC mineral dust (Hess et al., 1998) — visible/NIR anchors (m ≈ 1.53 + 0.008i at 0.55 µm).
+ * - Sokolik & Toon (1999), JGR 104, 9423–9444 — spectral behavior and strong Si–O features in 8–12 µm.
+ * - Querry (1985/1991), Optical constants of minerals — silicate trends in MIR/LWIR (quartz/feldspars/clays).
+ * - Patterson, Winchester & Barlow (1977), and Shettle & Fenn (1979) — atmospheric dust guidance.
+ *
+ * Table design:
+ * - 200–20,000 nm coverage; linear interpolation by ParameterLookupTable.
+ * - Preserves your original 400–700 nm values exactly to avoid visible-range regressions in graphics.
+ * - Denser sampling near 3 µm OH-complex and the strong 9–10 µm Si–O reststrahlen band.
+ */
+const DUST_REFRACTION_INDEX_LUT = ParameterLookupTable.from(
+    2, // number of columns in the table [n, k]
+    // Lookup data (n, k) per wavelength below. Wavelength positions are in nm.
+    [
+        // UV (short-wave): increasing absorption towards deep UV
+        1.550, 6.0e-2, // 200 nm
+        1.540, 3.5e-2, // 220 nm
+        1.535, 2.0e-2, // 240 nm
+        1.530, 1.5e-2, // 260 nm
+        1.525, 1.1e-2, // 280 nm
+        1.523, 9.5e-3, // 300 nm
+        1.522, 9.0e-3, // 320 nm
+        1.521, 8.7e-3, // 340 nm
+        1.5205, 8.5e-3, // 360 nm
+        1.5200, 8.4e-3, // 380 nm
+        // Visible (preserve exact legacy values for 400–700 nm)
+        1.5200, 8.2e-3, // 400 nm
+        1.5200, 7.0e-3, // 450 nm
+        1.5200, 6.0e-3, // 500 nm
+        1.5200, 5.0e-3, // 550 nm
+        1.5100, 4.0e-3, // 600 nm
+        1.5100, 3.5e-3, // 650 nm
+        1.5100, 3.0e-3, // 700 nm
+        // NIR/SWIR (0.7–2.5 µm): modest absorption; slight dispersion in n
+        1.5150, 3.2e-3, // 750 nm
+        1.5140, 3.0e-3, // 800 nm
+        1.5130, 3.2e-3, // 900 nm
+        1.5120, 3.5e-3, // 1000 nm
+        1.5120, 3.8e-3, // 1064 nm
+        1.5110, 4.5e-3, // 1200 nm
+        1.5100, 5.5e-3, // 1300 nm
+        1.5090, 8.0e-3, // 1450 nm
+        1.5080, 7.0e-3, // 1550 nm
+        1.5070, 6.5e-3, // 1650 nm
+        1.5060, 6.2e-3, // 1750 nm
+        1.5050, 7.0e-3, // 1940 nm
+        1.5050, 7.5e-3, // 2000 nm
+        1.5040, 9.0e-3, // 2200 nm
+        1.5030, 1.2e-2, // 2400 nm
+        1.5030, 1.4e-2, // 2500 nm
+        // Mid-IR (2.5–5 µm): increasing absorption toward lattice/OH bands
+        1.5030, 2.0e-2, // 2700 nm
+        1.5040, 3.0e-2, // 2900 nm
+        1.5050, 3.5e-2, // 3000 nm
+        1.5060, 4.0e-2, // 3100 nm
+        1.5080, 4.5e-2, // 3200 nm
+        1.5100, 5.0e-2, // 3300 nm
+        1.5120, 5.5e-2, // 3500 nm
+        1.5150, 6.0e-2, // 3700 nm
+        1.5180, 7.0e-2, // 4000 nm
+        1.5220, 8.0e-2, // 4500 nm
+        1.5250, 9.0e-2, // 5000 nm (5 µm)
+        // Long-wave IR (5–20 µm): strong Si–O stretch band near 9–10 µm
+        1.5300, 1.2e-1, // 6000 nm
+        1.5400, 1.8e-1, // 7000 nm
+        1.5500, 3.0e-1, // 8000 nm
+        1.5600, 4.0e-1, // 8500 nm
+        1.5700, 5.5e-1, // 9000 nm
+        1.5800, 7.0e-1, // 9500 nm
+        1.5900, 8.0e-1, // 10000 nm (10 µm)
+        1.5800, 5.0e-1, // 11000 nm
+        1.5600, 3.0e-1, // 12000 nm
+        1.5300, 1.5e-1, // 15000 nm
+        1.5000, 8.0e-2, // 20000 nm (20 µm)
+    ],
+    // Wavelength positions (nm) corresponding to the (n,k) pairs above
+    [
+        200,
+        220,
+        240,
+        260,
+        280,
+        300,
+        320,
+        340,
+        360,
+        380,
+        400,
+        450,
+        500,
+        550,
+        600,
+        650,
+        700,
+        750,
+        800,
+        900,
+        1000,
+        1064,
+        1200,
+        1300,
+        1450,
+        1550,
+        1650,
+        1750,
+        1940,
+        2000,
+        2200,
+        2400,
+        2500,
+        2700,
+        2900,
+        3000,
+        3100,
+        3200,
+        3300,
+        3500,
+        3700,
+        4000,
+        4500,
+        5000,
+        6000,
+        7000,
+        8000,
+        8500,
+        9000,
+        9500,
+        10000,
+        11000,
+        12000,
+        15000,
+        20000
+    ]
+);
+
+/**
+ *  💨 MINERAL DUST
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_dust(result, wavelength_nm) {
+    DUST_REFRACTION_INDEX_LUT.sample(wavelength_nm, result);
+}

package/src/core/math/physics/mie/ri_pollen.d.ts

@@ -0,0 +1,7 @@
+/**
+ * 🌿 POLLEN (Organic)
+ * @param {number[]} result complex number, refraction index
+ * @param {number} wavelength_nm in nanometers
+ */
+export function ri_pollen(result: number[], wavelength_nm: number): void;
+//# sourceMappingURL=ri_pollen.d.ts.map

package/src/core/math/physics/mie/ri_pollen.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ri_pollen.d.ts","sourceRoot":"","sources":["../../../../../../src/core/math/physics/mie/ri_pollen.js"],"names":[],"mappings":"AA0IA;;;;GAIG;AACH,kCAHW,MAAM,EAAE,iBACR,MAAM,QAIhB"}