Micronozzle Acceleration (MNA) Overview

• Micronozzle Acceleration (MNA) is a technique that uses micron-scale nozzle geometries to convert thermal, phase-change, or field energy into directed kinetic energy in fluids or ions.
• It leverages diverse mechanisms—from CFD-optimized cold-gas flow and non-equilibrium molecular dynamics to thermocavitation and laser–plasma interactions—to achieve high thrust, supersonic jets, and GeV-level ion beams.
• Optimization of nozzle parameters, such as throat curvature and expansion ratio, is critical for enhancing performance in applications like microsatellite propulsion, needle-free injection, and compact accelerator designs.

Micronozzle acceleration (MNA) refers to a set of physical mechanisms whereby the geometry of a nozzle at the micron or sub-millimeter scale is tailored to significantly accelerate fluids or charged particles. MNA exploits both continuum and non-equilibrium effects to convert internal (thermal or phase-change) or field (laser, plasma) energy into directed kinetic energy, achieving outcomes such as high thrust in microsatellite cold-gas propulsion, ballistic microjets for needle-free injection, supersonic gas cooling for molecular beams, and multi-hundred MeV–GeV ion acceleration. The diversity of MNA implementations spans classical compressible flow, thermocavitation hydrodynamics, atomistic molecular dynamics, and laser–plasma interaction physics.

1. Continuum Micronozzle Acceleration in Cold Gas Micro-Propulsion

In microsatellite cold-gas propulsion, MNA is realized via sub-millimeter convergent–divergent “Laval” nozzles optimized to maximize thrust as a function of geometry, propellant, and operating conditions. Key geometric parameters include throat radius of curvature ($R_t$), throat and exit widths ($W_t$, $W_e$), expansion ratio ($\varepsilon = A_e/A_t$), and, in advanced designs, placement and dimensions of dual-throat stages.

CFD-driven sensitivity analysis and response-surface optimization demonstrate that introducing a curved throat of $R_t \simeq D_t$ (where $D_t$ is the throat diameter) yields a 25% thrust enhancement (from 113.1 mN to 141 mN) over sharp-cornered geometries, primarily by mitigating adverse pressure gradients and extending boundary-layer development (Niksirat, 2023). Implementation of a dual-throat nozzle, where a secondary convergent section downstream of the first throat is positioned at $L_{\rm conv,2} \approx 374$ μm and a divergence length $L_a \approx 374$ μm (for $D_t\sim1$ mm), induces re-compression and realigns oblique expansion shocks, suppressing flow separation and boosting thrust to 261 mN.

Governing equations for steady compressible flow include conservation of mass, momentum (Navier–Stokes), and energy, subject to turbulence (k–ε model), ideal gas assumptions, and adiabatic walls. The thrust can be computed as

$T = \dot{m} (v_e - v_0) + (p_e - p_a)A_e,$

with isentropic relations

$W_t$0

Statistical DOE and response-surface methods reveal that $W_t$1 (exit width) contributes $W_t$2 to thrust sensitivity, $W_t$3 (throat width) $W_t$4, with remaining influence from $W_t$5, $W_t$6, and inlet width. Optimal expansion ratios $W_t$7 are in the range $W_t$8–$W_t$9, while scaling laws prescribe that $W_e$0 (thrust coefficient) can be enhanced by 15–20% via dual-throat staging. Maximum thrust is

$W_e$1

Boundary-layer control, shock alignment, and careful contouring are critical to maximizing acceleration and minimizing adverse separation, especially at low Reynolds numbers characteristic of micro-propulsion (Niksirat, 2023).

2. Non-Equilibrium and Atomistic Regimes: Molecular Dynamics of Microscopic Laval Nozzles

When nozzle dimensions approach the mean free path (high Knudsen number, $W_e$2–$W_e$3), as in the molecular dynamics study of slit Laval nozzles (Ortmayer et al., 2023), MNA enters a non-equilibrium regime with pronounced molecular-scale effects. Here, a simple monoatomic Lennard–Jones fluid is driven from a GCMC inlet reservoir through a micrometer- to nanometer-scale nozzle with atomically smooth walls; the outgoing flow is sampled as a stationary non-equilibrium ensemble (NEMD).

Key findings include:

• Supersonic acceleration of gas even for throats of $W_e$4–$W_e$5 molecular diameters ($W_e$6), with local Mach number $W_e$7–5 in large nozzles and $W_e$8 in smaller ones.
• The location of the sonic horizon ($W_e$9 where $\varepsilon = A_e/A_t$0) shifts downstream by $\varepsilon = A_e/A_t$1–$\varepsilon = A_e/A_t$2 for the smallest nozzles.
• Temperature anisotropy develops in the supersonic expansion: $\varepsilon = A_e/A_t$3, $\varepsilon = A_e/A_t$4, $\varepsilon = A_e/A_t$5 components differ by 10–20%, with velocity distributions remaining Maxwellian within each direction.
• Spatiotemporal density correlation analysis exposes a well-defined, one-way sonic horizon, whereby upstream density fluctuation correlations vanish downstream.
• Partial ballistic transport (remnants of two-body scattering) is observed, along with a tendency for transient condensation in the diverging section.

The isentropic area–velocity and Mach–area relations,

$\varepsilon = A_e/A_t$6

remain qualitatively valid for larger nozzles ($\varepsilon = A_e/A_t$7), while for smaller ones, rarefaction effects, wall slip, and delayed sonic transition become significant. This suggests micro-nozzle design requires explicit consideration of Knudsen number, thermal anisotropy, and wall characteristics (Ortmayer et al., 2023).

3. Thermocavitation-Driven Micronozzle Acceleration in Microfluidics

In microfluidic jetting, MNA is achieved by laser-induced thermocavitation: tightly focused CW laser irradiation (λ ≈ 450 nm, $\varepsilon = A_e/A_t$8 400–600 mW) into dye-loaded microchambers (depth $\varepsilon = A_e/A_t$9m) rapidly nucleates vapor bubbles, driving phase-change-initiated expulsion of liquid through micro-nozzles (diameter $R_t \simeq D_t$0m) (Galvez et al., 2020). The resulting pressure pulse and meniscus curvature generate high-velocity collimated jets.

Fluid acceleration is dominated by inertial effects (Ohnesorge number $R_t \simeq D_t$1), with typical jet velocities $R_t \simeq D_t$2–$R_t \simeq D_t$3 m/s and accelerations $R_t \simeq D_t$4–$R_t \simeq D_t$5 m/s$R_t \simeq D_t$6. The core acceleration relation is

$R_t \simeq D_t$7

where $R_t \simeq D_t$8 is the pre-expansion fill height. Performance further scales as $R_t \simeq D_t$9 and increases with taper angle, with experiments reporting a gain factor up to 2 for $D_t$0 in the $D_t$1–$D_t$2 range.

Boundary-integral (BI) potential flow models capture the meniscus evolution and jet ejection quantitatively, assuming incompressible, irrotational flow with dynamic boundary conditions set by $D_t$3 (bubble pressure) and interfacial tension. The parameter space is governed by critical dimensionless numbers:

• Ohnesorge $D_t$4
• Weber $D_t$5
• Reynolds $D_t$6

Optimal design for needle-free jet injection and high-fidelity inkjet printing is set by selecting $D_t$7–$D_t$8m, $D_t$9–$L_{\rm conv,2} \approx 374$0, and $L_{\rm conv,2} \approx 374$1m for stable jets with $L_{\rm conv,2} \approx 374$2–$L_{\rm conv,2} \approx 374$3 m/s, $L_{\rm conv,2} \approx 374$4–$L_{\rm conv,2} \approx 374$5, $L_{\rm conv,2} \approx 374$6–$L_{\rm conv,2} \approx 374$7 (Galvez et al., 2020).

4. Laser–Plasma Micronozzle Acceleration for Proton Beams

At relativistic intensities ($L_{\rm conv,2} \approx 374$8 W/cm$L_{\rm conv,2} \approx 374$9), MNA manifests as a three-stage laser–plasma ion acceleration process in structured targets consisting of a micron-scale hydrogen rod embedded within a hollow aluminum nozzle (Murakami et al., 30 Nov 2025).

The target geometry comprises a nozzle with head, neck, and skirt shaped to focus laser energy and guide hot-electron outflow:

• Nozzle head height $L_a \approx 374$0m, neck $L_a \approx 374$1m, skirt $L_a \approx 374$2m
• H-rod diameter $L_a \approx 374$3m, wall thickness $L_a \approx 374$4m

Upon irradiation, the sequence is:

1. Run-up phase ($L_a \approx 374$5 fs): Laser generates relativistic electrons, sheath field initiates proton expansion.
2. Main-drive phase ($L_a \approx 374$6–250 fs): Hot electrons charge the nozzle tail, establishing a longitudinal electrostatic field $L_a \approx 374$7 V/m.
3. Afterburner phase ($L_a \approx 374$8 fs): Comoving field persists, protons absorb further energy from electron thermal expansion, contributing $L_a \approx 374$9300 MeV extra gain.

Particle-in-cell (PIC) simulations reveal that, at $D_t\sim1$0 W/cm$D_t\sim1$1, protons reach cutoff energies $D_t\sim1$2 GeV with characteristic mid-spectrum plateaus ($D_t\sim1$3–$D_t\sim1$4 MeV), angular divergence FWHM $D_t\sim1$5, and conversion efficiency $D_t\sim1$6. Scaling shows $D_t\sim1$7 (plane) and up to $D_t\sim1$8 (Gaussian spot), exceeding typical TNSA behavior (linear scaling). Nozzle shape and rod–nozzle alignment modulate $D_t\sim1$9 by up to 10%. The self-similar plasma expansion (afterburner) is analytically captured as

$T = \dot{m} (v_e - v_0) + (p_e - p_a)A_e,$0

with $T = \dot{m} (v_e - v_0) + (p_e - p_a)A_e,$1 electron temperature, $T = \dot{m} (v_e - v_0) + (p_e - p_a)A_e,$2 sound speed, and $T = \dot{m} (v_e - v_0) + (p_e - p_a)A_e,$3 Debye length. These results indicate that structured MNA targets achieve $T = \dot{m} (v_e - v_0) + (p_e - p_a)A_e,$4–$T = \dot{m} (v_e - v_0) + (p_e - p_a)A_e,$5 higher energies than H-rod or foil alone under identical laser drive (Murakami et al., 30 Nov 2025).

5. Optimization Strategies and Performance Sensitivities

Across physical realizations, MNA efficiency depends critically on geometry, staging, and material/system selection. In continuum micro-nozzles, response surface optimization highlights the dominance of exit width and throat curvature. For atomistic nozzles, large nozzle widths relative to molecular scale ($T = \dot{m} (v_e - v_0) + (p_e - p_a)A_e,$6) are essential to approach isentropic acceleration. In thermocavitation MNA, controlling fill height and taper angle achieve maximal jet speed and collimation, while Ohnesorge and Weber numbers inform stability margins.

In laser–plasma MNA, rod/nozzle geometry, relative positioning (e.g., rod–nozzle gap $T = \dot{m} (v_e - v_0) + (p_e - p_a)A_e,$7), and target material composition can modulate the field amplification and energy transfer, with elliptic rods and tight contact optimizing $T = \dot{m} (v_e - v_0) + (p_e - p_a)A_e,$8 and conversion efficiency (Niksirat, 2023, Ortmayer et al., 2023, Galvez et al., 2020, Murakami et al., 30 Nov 2025).

6. Applications and Implications

MNA underpins key advances in:

• Spacecraft cold-gas micropropulsion: Achieving thrust in the 0.1–0.25 N class for AGN and attitude control thrusters (Niksirat, 2023).
• Miniaturized beam sources and microfluidics: Generating supersonic beams and precision jets for molecular spectroscopy, nano-printing, and needle-free drug delivery (Ortmayer et al., 2023, Galvez et al., 2020).
• Relativistic ion acceleration: Delivering $T = \dot{m} (v_e - v_0) + (p_e - p_a)A_e,$9 GeV proton beams for hadron therapy, compact accelerator injectors, and high-energy-density science at moderate divergence and conversion efficiency (Murakami et al., 30 Nov 2025).

Limitations include requirements for advanced microfabrication, precise flow or field alignment, and addressing pronounced non-equilibrium effects at the molecular scale or in laser–plasma targets.

7. Scaling Laws, Design Guidelines, and Future Directions

The scaling of MNA output metrics as a function of system size, drive energy, and geometric ratios is determinative for both performance and practical realization:

• In gas flows, optimal expansion ratio and throat curvature maximize enthalpy conversion and boundary-layer management.
• In atomistic and non-equilibrium MNA, minimizing Knudsen number and maximizing active area relative to molecular scale restores continuum efficacy.
• For thermocavitation, the inverse dependency on fill height ($W_t$00), and the gain with taper angle, are dominant.
• In laser-plasma MNA, proton energy scales sub-linearly to super-linearly with laser intensity, and geometric optimizations yield up to $W_t$01 enhancements.

A plausible implication is that future MNA research will emphasize integrative optimization over geometry, material, and drive modalities, as well as extension to arrayed and three-dimensional micro-nozzle systems for enhanced output and new functionalities (Niksirat, 2023, Ortmayer et al., 2023, Galvez et al., 2020, Murakami et al., 30 Nov 2025).