Prescribed Wall-Heat-Flux Control of Blockage and Impulse in a Rarefied Micro-Nozzle

Prescribed wall heat flux provides an active route for controlling rarefied micro-nozzle flows, but its effect is governed by the coupled wall–bulk thermal response rather than by the imposed flux alone. This work uses direct simulation Monte Carlo (DSMC) simulations to study nitrogen flow in a converging–diverging micro-nozzle with cooling, adiabatic, and heating applied on the diverging wall. The imposed heat flux is scaled by the inlet kinetic-energy flux, $E=0.5{\rho }_i{U}_i^3$, giving $Q_w\mathrm{/}E$ from $-10.5\mathrm{\%}$ to $97.3\mathrm{\%}$; this range spans moderate cooling, weak-to-intermediate heating, and a near-unity thermal-forcing regime. Wall and mass-flux-weighted bulk temperature profiles, film-temperature-based Nusselt and local-viscosity Brinkman-type diagnostics, gradient-length Knudsen indicators, mass-flux thickness, thrust decomposition, and proper orthogonal decomposition (POD) of signed numerical schlieren are analyzed. The results show that heating creates strong wall–bulk stratification: the wall temperature exceeds five times the inlet value, while the bulk temperature responds more gradually. Cooling cases contain locations where $T_w-T_b$ changes sign, making the local Nusselt-type response singular; the raw singular behavior is retained for diagnosis and a validity mask is used only for comparative plotting. Heating contracts the effective mass-carrying core, increasing aerodynamic blockage and reducing mass flow rate. However, strong heating increases the specific impulse from $156$ s to $201$ s because thermal and pressure-thrust augmentation outweigh the mass-flow penalty. The internal compression feature evolves into a finite viscous–thermal compression zone, and its heat-flux-parametric response remains low-dimensional, with the first two POD modes capturing more than $97\mathrm{\%}$ of the fluctuation energy.