Year: 2026
Type: preprint
Abstract: We present a direct simulation Monte Carlo (DSMC) method for the kinetic theory of ordered fluids recently introduced by the authors. The microscopic constituents of these fluids carry an order parameter belonging to a manifold. The method applies Strang splitting to couple a Monte Carlo treatment of the binary collisions with a symplectic integration of the mean-field Vlasov transport on the order-parameter manifold. We use the Nanbu-Babovsky approach for Maxwellian collision kernels, and Bird's acceptance-rejection method for general ones. The mean-field Vlasov force is reconstructed from the particle ensemble at each time step. We also include a Bhatnagar-Gross-Krook (BGK) / Andersen thermostat to control the system temperature. We formulate the method for a generic order-parameter manifold, and specialize it to two-dimensional calamitic (rod-like) molecules. Numerical tests validate the method. We confirm the relaxation to the Maxwellian distribution predicted by the H-theorem, and simulate the mean-field alignment driven by quadratic, Kuramoto, and Onsager interaction potentials. With the Onsager potential we successfully capture the isotropic-nematic phase transition at the correct critical temperature.
Source: arXiv (Cornell University)
Language: English
Subfield: Applied Mathematics
Field: Mathematics
Domain: Physical Sciences
Open Access status: green
Awards: DMS-EPSRC: Stability Analysis for Nonlinear Partial Differential Equations across Multiscale Applications, G2691952350, G3452655474, G4969370889, G5180915604, G5195188938, G567365627, SysGenX: Composable software generation for system-level simulation at Exascale, G6049851770, G6115201310, G6328286017, G7853015121
OpenAlex ID: W7214585950