Hostname: page-component-7c8c6479df-r7xzm Total loading time: 0 Render date: 2024-03-28T19:23:46.270Z Has data issue: false hasContentIssue false

An approach to entropy consistency in second-order hydrodynamic equations

Published online by Cambridge University Press:  01 March 2004

RAMESH BALAKRISHNAN
Affiliation:
Center for Simulation of Advanced Rockets, University of Illinois at Urbana–Champaign, Illinois 61801, USA Present address: National Center for Supercomputing Applications, 152 Computing Applications Building, 605 West Springfield Avenue, Champaign, Illinois 61820, USA.

Abstract

As the flow becomes rarefied it has been seen that predictions of continuum formulations, such as the Navier–Stokes equations, become inaccurate. These inaccuracies stem from the linear approximations to the stress and heat flux in the viscous flux terms in the Navier–Stokes equations. Hence, it has long been conjectured that the inclusion of higher-order terms in the constitutive relations for the stress and heat flux may improve the predictive capabilities of such continuum formulations. Following this approach, second-order systems of hydrodynamic equations, such as the Burnett and Woods equations, were applied to the shock structure problem. While it was observed that these equations afforded a better description of the shock structure on coarse grids, they were prone to small wavelength instabilities when the grids were refined. The cause of this instability was subsequently traced to the fact that these equations can potentially violate the second law of thermodynamics when the local Knudsen number exceeds a critical limit. This leads to the fundamental question: is entropy consistency achievable in a system of second-order hydrodynamic equations? To answer this question, a novel set of equations, known as the BGK-Burnett equations, is constructed by taking moments of the Boltzmann equation for the second-order distribution function.

The formulation of second-order hydrodynamic equations by moment methods is beset by three hurdles: (i) the highly nonlinear collision integral in the Boltzmann equation needs to be evaluated, (ii) the second-order distribution function does not satisfy the moment closure criterion and (iii) identification of the approximations to the material derivatives in the second-order distribution function that will correctly account for the difference in time scales between the first- and second-order fluxes. The first three terms of the Chapman–Enskog expansion, that defines the second-order distribution function, and the Bhatnagar–Gross–Krook model of the collision integral, form the basis of the BGK-Burnett equations. The entropy-consistent behaviour of the equations depend on the moment closure coefficients and the approximations to the material derivatives in the second-order fluxes. The requirement of moment closure alone, however, results in non-unique closure coefficients and a family of BGK-Burnett equations, from which an entropy-consistent set must be identified. From this family, two sets of BGK-Burnett equations have been considered, and this paper presents the details of the formulation of these two sets of equations, the identification of entropy-consistent approximations to the material derivatives by a novel entropy consistent relaxation technique, and shock structure computations in a monatomic hard sphere gas for a range of Mach numbers.

Type
Papers
Copyright
© 2004 Cambridge University Press

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)