Return to search

An Approximation for the Twenty-One-Moment Maximum-Entropy Model of Rarefied Gas Dynamics

The use of moment-closure methods to predict continuum and moderately rarefied flow offers
many modelling and numerical advantages over traditional methods. The maximum-entropy
family of moment closures offers models described by hyperbolic systems of balance
laws. In particular, the twenty-one moment model of the maximum-entropy hierarchy offers a
hyperbolic treatment of viscous flows exhibiting heat transfer. This twenty-one moment
model has the ability to provide accurate solutions where the Navier-Stokes equations lose
physical validity due to the solution being too far from local equilibrium. Furthermore,
its first-order hyperbolic nature offers the potential for improved numerical accuracy as
well as a decreased sensitivity to mesh quality. Unfortunately, higher-order
maximum-entropy closures cannot be expressed in closed form. The only known affordable
option is to propose approximations. Previous approximations to the fourteen-moment
maximum-entropy model have been proposed [McDonald and Torrilhon,
2014]. Although this fourteen-moment model also predicts viscous flow with heat
transfer, the necessary moments to close the system renders it more difficult to
approximate accurately than the twenty-one moment model. The proposed approximation for
the fourteen-moment model also has realizable states for which hyperbolicity is lost.

Unfortunately, the velocity distribution function associated with the twenty-one moment
model is an exponential of a fourth-order polynomial. Such a function cannot be integrated
in closed form, resulting in closing fluxes that can only be obtained through complex
numerical methods. The goal of this work is to present a new approximation to the closing
fluxes that respect the maximum-entropy philosophy as closely as possible. Preliminary
results show that a proposed approximation is able to provide shock predictions that are
in good agreement with the Boltzmann equation and surpassing the prediction of the
Navier-Stokes equations. Furthermore, Couette flow results as well as lid-driven cavity
flows are computed using a novel approach to Knudsen layer boundary conditions. A
dispersion analysis as well as an investigation of the hyperbolicity of the model is also
shown. The Couette flow results are compared against Navier-Stokes and the free-molecular
analytical solutions for a varying Knudsen number, for which the twenty-one moment model
show good agreement over the domain. The shock-tube problem is also computed for different
Knudsen numbers. The results are compared with the one obtained by directly solving the BGK
equation. Finally, the lid-driven cavity flow computed with the twenty-one moment model
shows good agreement with the direct simulation Monte-Carlo (DSMC) solution.

Identiferoai:union.ndltd.org:uottawa.ca/oai:ruor.uottawa.ca:10393/45657
Date23 November 2023
CreatorsGiroux, Fabien
ContributorsMcDonald, James Gerald
PublisherUniversité d'Ottawa / University of Ottawa
Source SetsUniversité d’Ottawa
LanguageEnglish
Detected LanguageEnglish
TypeThesis
Formatapplication/pdf
RightsAttribution-NonCommercial-NoDerivatives 4.0 International, http://creativecommons.org/licenses/by-nc-nd/4.0/

Page generated in 0.0025 seconds