Resolving Cryogenic and Hypersonic Rarefied Flows via Deep Learning-Accelerated Lennard-Jones DSMC

Abstract

Integrating the physically realistic Lennard--Jones (LJ) potential into Direct Simulation Monte Carlo (DSMC) remains challenging because the long-range potential complicates collision-rate definition and makes repeated scattering-angle evaluation expensive. This study develops an LJ--DSMC framework built around two methodological advances and a transport-level validation of the resulting collision kernel. First, a generalized collision-selection treatment is formulated for Bird's DSMC algorithms (DSMC1, DSMC1S, and DS2V) through a Variable Effective Diameter (VED) model obtained from local Chapman--Enskog viscosity matching. This viscosity-consistent pair-selection model provides a finite DSMC collision-rate closure for the LJ potential and is validated in helium and argon normal shocks, cryogenic supersonic Couette flow, and hypersonic cylinder flows. The results show agreement with VHS in high-temperature repulsive regimes, but reveal clear LJ effects, including reduced shear stress and larger cryogenic wakes, when attractive forces become important. Second, the computational bottleneck of the accepted LJ binary-scattering step is removed by training a Deep Operator Network (DeepONet) to predict the LJ deflection angle from high-fidelity scattering data, replacing the numerical Matsumoto--Koura integral while preserving the standard elastic post-collision update. The surrogate gives a bulk mean wrapped-angle error of \(1.6\times10^{-3}\,\mathrm{rad}\) and a 99th-percentile error of \(9.9\times10^{-3}\,\mathrm{rad}\), accelerates the collision subroutine by 40\%, and reduces total wall time by 36\%. Finally, the same DeepONet--LJ scattering kernel is tested beyond viscosity-controlled flows through diffusion benchmarks.