Calculadora Inestabilidad Richtmyer-Meshkov Shock

Calculadora inestabilidad Richtmyer-Meshkov RMI shock-driven: γ = k·A*·Δu. ICF NIF, SN 1987A, ENEA Frascati DTT, Avio Vega, Cineca Leonardo.

Parámetros de cálculo
¿Te ha sido útil esta calculadora?

Calculadora inestabilidad Richtmyer-Meshkov shock

La calculadora inestabilidad Richtmyer-Meshkov RMI determina la tasa de crecimiento lineal γ_RM de una perturbación en la interfaz entre dos fluidos de densidad diferente atravesada por una onda de choque, según el modelo lineal de Robert D. Richtmyer 1960 (Communications on Pure and Applied Mathematics, vol. 13, p. 297, "Taylor instability in shock acceleration of compressible fluids", Courant Institute NYU + Los Alamos National Laboratory) confirmado experimentalmente por Evgeny E. Meshkov 1969 (Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, vol. 4, p. 151, VNIIEF Sarov shock tube He/SF6) primera evidencia laboratorio.

Fórmula clásica RMI lineal: γ_RM = k · A* · Δu, donde k = 2π/λ es el número de onda de la perturbación inicial (rad/m), A* = (ρ₂ - ρ₁) / (ρ₂ + ρ₁) es el número de Atwood post-shock (compresibilidad modifica A pre/post), Δu es el salto de velocidad impartido por el shock al contacto (m/s, derivado de relaciones Rankine-Hugoniot para número Mach M_s y razón presiones p₂/p₁).

Diferencia fundamental respecto Rayleigh-Taylor RTI: RMI guiada por impulso instantáneo (shock) no por gravedad continua → crecimiento lineal en tiempo h(t) = h₀ + γ_RM·t (NO exponencial como RTI), y crece incluso con A* < 0 (light-to-heavy o heavy-to-light) porque shock baroclínico genera vorticidad en ambas direcciones.

Régimen no-lineal: cuando amplitud alcanza Δh ~ λ/π → transición a estructuras bubble-spike, acoplamiento modos, finalmente turbulencia Kolmogorov totalmente desarrollada con cascada energética ε ~ Δu³/L.

Física baroclínica: shock genera vorticidad ω = ∇ρ × ∇p / ρ² en interfaz no alineada, depositando instantáneamente energía rotacional que alimenta mezcla turbulenta.

4 categorías aplicaciones:

  • ICF Inertial Confinement Fusión NIF Lawrence Livermore + LMJ francés + ENEA Casaccia 1962 ABC laser, cápsulas DT (deuterio-tritio) implosionan con A=0,7 Δu=10⁵ m/s λ=10⁻⁶ m → γ=7×10⁸ s⁻¹ (mezcla del ablator plástico en el fuel arruina ignition)
  • Astrofísica supernova SN 1987A LMC blast through stellar envelope, A=0,5 Δu=10⁶ m/s λ=10¹⁰ m → γ=10⁻⁴ s⁻¹ (mezcla layer Fe/Si/O/He visible clumpy ejecta INAF Roma 2002 Hubble + JWST 2021)
  • Aeroespacial hipersónico scramjet/ramjet fuel-air mixing Mach 5-15, A=0,3 Δu=300-3000 m/s, combustores turbulentos Avio SpA Colleferro 1908 Vega launcher + Ansaldo Energía Génova 1853 turbinas
  • Sonoluminiscencia shock-driven cavitation bubble collapse + atmospheric meteoritic impact mixing (Tunguska 1908, Cheliábinsk 2013). Las aplicaciones italianas RMI:
  • ENEA Frascati DTT 2028 + FTU 1989 tokamak ELM (Edge Localized Mode) shock-driven plasma instabilities + ENEA Casaccia 1962 ABC laser fusión ICF Italia (cápsula DT 1 J 1,06 μm)
  • INAF Roma 2002 + INAF Capodimonte 1812 SN 1987A LMC clumpy ejecta + SN 1572 Tycho + SN 1604 Kepler shells RMI
  • INFN-LNGS Gran Sasso XENONnT + Cineca Marconi100 30 PFlops + Leonardo HPC pre-exascale 250 PFlops simulaciones multifase RMI 3D
  • Avio SpA Colleferro 1908 Vega-C P120C solid motor + Vega-E Z40 LOX/metano combustores RMI mezcla
  • Ansaldo Energía Génova 1853 GT36 H-class 538 MW combustores turbulentos DLN dry-low-NOx
  • Politecnico Milano 1863 + Torino 1859 + Bologna 1088 + Pisa Scuola Normale 1810 CFD RMI shock tube He/SF6/Air
  • ICRA Sapienza 1985 + ICRANet Pescara 2003 + SISSA Trieste 1978 + ICTP Trieste 1964 Remo Ruffini astrophysical shocks SN remnants. Escalas típicas: laboratorio shock tube He-air (γ=560 s⁻¹), ICF NIF cápsula (γ=7×10⁸ s⁻¹), SN 1987A astrofísica (γ=10⁻⁴ s⁻¹)

Tabla RMI italiana escenarios ENEA INAF Cineca

EscenarioA* AtwoodΔu (m/s)λ (m)γ_RM (s⁻¹)Aplicación italiana
Lab shock tube He-aire0,33001×10⁻²5,6×10²Politecnico Milano CFD
ICF NIF DT cápsula0,71×10⁵1×10⁻⁶4,4×10¹¹ENEA Casaccia ABC laser
ELM tokamak DTT0,91×10⁴1×10⁻³5,7×10⁷ENEA Frascati DTT 2028
Scramjet Mach 8 H₂-aire0,515005×10⁻³9,4×10⁵Avio Vega hipersónico
GT36 turbina combustor0,22005×10⁻³5,0×10⁴Ansaldo Energía Génova
SN 1987A blast envelope0,51×10⁶1×10¹⁰3,1×10⁻⁴INAF Roma 2002 Hubble
SN 1572 Tycho remnant0,43×10⁶1×10¹⁵7,5×10⁻⁹INAF Capodimonte INFN
Sonoluminiscencia H₂O0,9915001×10⁻⁶9,3×10⁹INRIM Torino metrología

(Fuente: Richtmyer 1960 + Meshkov 1969 + ENEA Frascati DTT + INAF Roma 2002 SN 1987A + Cineca Leonardo HPC + Avio Vega Ansaldo Energía).

Ejemplos cálculo italiano RMI shock

  • ICF NIF DT cápsula + ENEA Casaccia 1962 ABC laser ICF Italia: A* = 0,7 (DT/ablator plástico CH ρ₂=1000 ρ₁=180 kg/m³ post-shock), Δu = 10⁵ m/s (impulso láser 1 MJ 192 beam Nd:YAG), λ = 10⁻⁶ m (perturbación superficie ablator micrométrica). k = 2π/λ = 6,28×10⁶ rad/m. γ_RM = k·A*·Δu = 6,28×10⁶ · 0,7 · 10⁵ = 4,4×10¹¹ s⁻¹ (tiempo característico 2,3 ps, mezcla del ablator plástico en fuel DT arruina hot-spot ignition, problema P0 ICF)
  • Supernova SN 1987A INAF Roma 2002 LMC 50 kpc Hubble + JWST 2021: A* = 0,5 (Fe-Ni core/H envelope), Δu = 10⁶ m/s (shock breakout), λ = 10¹⁰ m (interfaz He/H pre-collapse 10 R_sol). γ_RM = 6,28×10⁻¹⁰ · 0,5 · 10⁶ = 3,14×10⁻⁴ s⁻¹ (1 hora). Mezcla layer Fe-Ni/Si/O/He visible clumpy ejecta 56Ni Bullet filaments
  • Lab shock tube He-aire Politecnico Milano CFD OpenFOAM: A* = 0,3, Δu = 300 m/s (Mach 1,2), λ = 1 cm. γ_RM = 628 · 0,3 · 300 = 5,65×10⁴ s⁻¹. Saturación t_NL ≈ 56 μs (validado fast camera schlieren)
  • ENEA Frascati DTT 2028 tokamak ELM shock + ENEA FTU 1989: A* = 0,9 (plasma 100 keV pedestal vs SOL), Δu = 10⁴ m/s (Alfvén ELM crash), λ = 1 mm SOL filaments. γ_RM = 5,7×10⁷ s⁻¹
  • Avio Vega-C scramjet H₂-aire Mach 8 Colleferro 1908: A* = 0,5, Δu = 1500 m/s, λ = 5 mm. γ_RM = 9,4×10⁵ s⁻¹ (mezcla rápida benéfica combustión completa)
  • Ansaldo Energía Génova 1853 GT36 H-class 538 MW DLN combustor: A* = 0,2, Δu = 200 m/s, λ = 5 mm. γ_RM = 5,0×10⁴ s⁻¹
  • Sonoluminiscencia H₂O cavitation + INRIM Torino 1923: A* = 0,99 (vacuum/agua), Δu = 1500 m/s (Mach 4 collapse), λ = 1 μm. γ_RM = 9,3×10⁹ s⁻¹ (T_max 10⁴ K plasma flash 10⁻¹² s)

Investigación italiana RMI Richtmyer-Meshkov

La investigación italiana RMI Richtmyer-Meshkov:

  • ENEA Agenzia Nazionale Frascati 1960 + Casaccia 1962 (4.500 investigadores, 12 centros): DTT Divertor Tokamak Test Frascati 2028 600 mln € PNRR Misión 2 fusión + FTU Frascati Tokamak Upgrade 1989-2020 ELM RMI shock plasma + ABC laser fusión Casaccia 1962 ICF Italia 1 J Nd:YAG 1,06 μm cápsula DT primera ICF Europa + RFX-mod Padua 2004 reverse field pinch CNR-INFN-ENEA
  • INAF Istituto Nazionale Astrofisica Roma 2002 (16 observatorios, 1.700 investigadores): SN 1987A INAF Trieste 1898 + Hubble 1990 + JWST 2021 NIRSpec RMI clumpy ejecta 56Ni Bullet filaments; SN 1572 Tycho remnant + SN 1604 Kepler shells INAF Bologna 1936 radio AGN Centaurus shocks; Crab Nebula M1 INAF Capodimonte 1812 pulsar wind shock RMI fingers; VST KIDS dark matter cluster Capodimonte; TNG La Palma 1998 3,58 m INAF; LBT Mount Graham Arizona 2005 2×8,4 m INAF 25%
  • INFN-LNGS Gran Sasso 1987 XENONnT WIMP + cosmic ray timing shocks
  • Avio SpA Colleferro 1908 Vega-C P120C solid motor 173 ton HTPB + Vega-E Z40 LOX/metano combustores RMI mezcla 25.000 empleados
  • Ansaldo Energía Génova 1853 GT36 H-class 538 MW DLN combustores turbulentos + AE94.3A 340 MW
  • Cineca Casalecchio Bologna 1969 consorcio universitario 70 universidades: Marconi100 30 PFlops 2020 + Leonardo HPC pre-exascale 250 PFlops 2023 EuroHPC (5° TOP500 mundial, 4.992 GPU NVIDIA A100, 13.824 CPU Intel Sapphire Rapids) simulaciones RMI 3D OpenFOAM Athena++ FLASH Code multifase reactiva shock-driven turbulence
  • Politecnico Milano 1863 DAER aeroespacial + Politecnico Torino 1859 DIMEAS + Università Bologna 1088 DIN + Pisa Scuola Normale 1810 + Sapienza Roma 1303 CFD RMI shock tube + scramjet + ICF
  • ICRA Sapienza 1985 + ICRANet Pescara 2003 + SISSA Trieste 1978 + ICTP Trieste 1964 Remo Ruffini astrophysical shocks SN remnants GRB jet-cloud interactions
  • INRIM Torino 1923 metrología time-of-flight shock pressure + AGCOM 1997 + AGCM Comunicaciones regulación
  • Premios Nobel italianos física: Marconi 1909, Fermi 1938, Segrè 1959, Carlo Rubbia 1984, Riccardo Giacconi 2002. Bonus fiscales italianos investigación 2025: Crédito de Impuesto I+D art. 1 c. 200 L. 160/2019 hasta 25%; Patent Box L. 178/2020 50%; Crédito Industria/Transición 5.0 art. 38 DL 39/2024 hasta 35%; PNRR Misión 2 hidrógeno + Misión 4 investigación 6,9 mil mill. €

Actualizado a julio de 2026

Preguntas frecuentes

¿Qué es la inestabilidad Richtmyer-Meshkov RMI?

Inestabilidad en la interfaz entre dos fluidos de densidad diferente atravesada por una onda de choque. Formulada por Robert D. Richtmyer 1960 (Comm. Pure Appl. Math.) Courant Institute NYU + LANL, confirmada experimental por E.E. Meshkov 1969 VNIIEF Sarov shock tube He/SF6. Difiere Rayleigh-Taylor: shock impulsivo no gravedad continua, crecimiento lineal h(t)=h₀+γt no exponencial.

¿Fórmula tasa crecimiento lineal RMI?

γ_RM = k · A* · Δu donde k = 2π/λ número de onda perturbación inicial (rad/m), A* = (ρ₂-ρ₁)/(ρ₂+ρ₁) número Atwood post-shock (compresibilidad modifica pre/post), Δu salto velocidad impartido por shock al contacto (m/s, de Rankine-Hugoniot Mach M_s). Crecimiento lineal h(t) = h₀ + γ_RM·t. Saturación no-lineal cuando h ~ λ/π → turbulencia Kolmogorov.

¿Diferencia RMI vs Rayleigh-Taylor RTI?

RTI guiada gravedad continua g → crecimiento exponencial h(t) = h₀·exp(γt) con γ = sqrt(k·A·g), requiere A>0 (heavy on light); RMI guiada impulso shock instantáneo → crecimiento lineal h(t) = h₀ + γ_RM·t, crece incluso con A<0 (light-to-heavy o heavy-to-light) por vorticidad baroclínica ω = ∇ρ × ∇p / ρ². RMI es el límite impulsivo de RTI con g·t = Δu.

¿Ejemplo RMI ICF NIF cápsula fusión?

ICF NIF Lawrence Livermore + ENEA Casaccia 1962 ABC laser: cápsula DT (deuterio-tritio) + ablator plástico CH, A*=0,7 Δu=10⁵ m/s λ=10⁻⁶ m. γ_RM = 6,28×10⁶ · 0,7 · 10⁵ = 4,4×10¹¹ s⁻¹ (2,3 ps). Mezcla del ablator plástico en fuel DT arruina hot-spot ignition: problema P0 ICF porque destruye confinamiento tiempo-densidad ρR > 0,3 g/cm² Lawson criterion.

¿RMI en supernova SN 1987A INAF?

SN 1987A LMC 50 kpc INAF Roma 2002 + Hubble 1990 + JWST 2021 NIRSpec: shock blast through stellar envelope A*=0,5 Δu=10⁶ m/s λ=10¹⁰ m → γ_RM = 3,14×10⁻⁴ s⁻¹. Mezcla de layer Fe-Ni/Si/O/He visible clumpy ejecta 56Ni Bullet filaments. Mismo mecanismo SN 1572 Tycho INAF Bologna + SN 1604 Kepler + Crab M1 Capodimonte pulsar wind shocks RMI fingers.

¿Aplicaciones italianas RMI ENEA Avio?

Italia RMI: 1) ENEA Frascati DTT 2028 600 mln€ PNRR + FTU 1989 ELM tokamak shocks; 2) ENEA Casaccia ABC laser ICF 1962; 3) INAF Roma SN 1987A clumpy ejecta; 4) Avio SpA Colleferro 1908 Vega-C P120C scramjet H₂-aire Mach 8 mezcla combustores; 5) Ansaldo Energía Génova 1853 GT36 H-class 538 MW DLN turbinas; 6) Cineca Marconi100 30 PFlops + Leonardo HPC 250 PFlops simulaciones 3D OpenFOAM.

¿Bonus fiscales investigación RMI Italia?

Bonus investigación Italia 2025 RMI: Crédito de Impuesto I+D art. 1 c. 200 L. 160/2019 hasta 25%; Patent Box L. 178/2020 50% deducción rentas IP; Crédito de Impuesto Industria/Transición 5.0 art. 38 DL 39/2024 hasta 35% para inversiones energy-saving (combustores Avio Ansaldo); PNRR Misión 2 hidrógeno + Misión 4 Componente 2 investigación 6,9 mil mill. € (DTT Frascati 600 mln €).

Comentarios (0)

Sin comentarios aún. ¡Sé el primero!