The paper introduces a mathematical model that describes the cavitation process occurring during the passage of a water steam flow in various geometric configurations of a hydrodynamic device. The flow experiences a localized constriction (convergent nozzle) followed by expansion (divergent nozzle), exemplified by a Venturi tube or a Laval nozzle. A narrow flow channel connecting the convergent and divergent sections is equipped with a narrow-section nozzle for injecting methane molecules into the high-speed steam flow. As the steam-gas mixture passes through this zone, it is irradiated with an electron beam and sprayed into a cylindrical chamber at atmospheric pressure, where the distribution of methane molecules in water vapor forms an aerosol. Key geometric parameters of the constriction and expansion zones of the hydraulic system (cavitation-jet chamber) are determined to ensure the uniform distribution of dispersed-phase particles (methane) in the dispersion medium (water vapor). Velocity and pressure distributions of the mixed steam-gas flow are calculated using a turbulent mathematical model, specifically the k-ω model, while the motion of methane particles is simulated using a particle tracing method. The uniformity of methane molecule distribution in water vapor is assessed using Ripley’s K-function. The best performance of the hydrogen-producing chamber was observed when the cavitation-inducing nozzle’s convergence angle exceeded 50 degrees. The divergence angle of the nozzle within the range of 30–40 degrees provided the best distribution in terms of uniformity of the methane particles in the chamber.