skip to main content


Title: Inverse design of mesoscopic models for compressible flow using the Chapman-Enskog analysis
Abstract

In this paper, based on simplified Boltzmann equation, we explore the inverse-design of mesoscopic models for compressible flow using the Chapman-Enskog analysis. Starting from the single-relaxation-time Boltzmann equation with an additional source term, two model Boltzmann equations for two reduced distribution functions are obtained, each then also having an additional undetermined source term. Under this general framework and using Navier-Stokes-Fourier (NSF) equations as constraints, the structures of the distribution functions are obtained by the leading-order Chapman-Enskog analysis. Next, five basic constraints for the design of the two source terms are obtained in order to recover the NSF system in the continuum limit. These constraints allow for adjustable bulk-to-shear viscosity ratio, Prandtl number as well as a thermal energy source. The specific forms of the two source terms can be determined through proper physical considerations and numerical implementation requirements. By employing the truncated Hermite expansion, one design for the two source terms is proposed. Moreover, three well-known mesoscopic models in the literature are shown to be compatible with these five constraints. In addition, the consistent implementation of boundary conditions is also explored by using the Chapman-Enskog expansion at the NSF order. Finally, based on the higher-order Chapman-Enskog expansion of the distribution functions, we derive the complete analytical expressions for the viscous stress tensor and the heat flux. Some underlying physics can be further explored using the DNS simulation data based on the proposed model.

 
more » « less
Award ID(s):
1706130
NSF-PAR ID:
10212405
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
Advances in Aerodynamics
Volume:
3
Issue:
1
ISSN:
2524-6992
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Todorov, M D (Ed.)
    Christov functions are a complete orthonormal set of functions on L^2(-∞,∞) that allow us to expand derivatives, nonlinear products, and nonlocal (integro-differential) terms back into the same basis. These properties are beneficial when solving nonlinear evolution equations using Galerkin spectral methods. In this work, we demonstrate such a “Christov expansion method” for the Benjamin–Ono (BO) equation. In the BO equation, the dispersion term is nonlocal, given by the Hilbert transform of the second spatial derivative of the unknown function. The Hilbert transform of the Christov functions can be computed using complex integration and Cauchy’s residue theorem to obtain simple relations. Then, a Galerkin spectral expansion can be used to the solve the BO equation. Time integration is performed using a Crank–Nicolson-type scheme. Importantly, the Christov expansion method yields a banded matrix for the spatial discretization, even though the spatial terms are nonlocal. To demonstrate the approach and its implementation, we perform numerical experiments showing the steady propagation of single and the overtaking interaction of multiple BO solitary waves. 
    more » « less
  2. For the quantum kinetic system modelling the Bose-Einstein Condensate that accounts for interactions between condensate and excited atoms, we use the Chapman-Enskog expansion to derive its hydrodynamic approximations, include both Euler and Navier-Stokes approximations. The hydrodynamic approximations describe not only the macroscopic behavior of the BEC but also its coupling with the non-condensates, which agrees with the Landau two-fluid theory. 
    more » « less
  3. Finch, a domain specific language and code generation framework for partial differential equations (PDEs), is demonstrated here to solve two classical problems: steady-state advection diffusion equation (single PDE) and the phonon Boltzmann transport equation (coupled PDEs). Both finite volume and finite element methods are explored. In addition to work presented at the 2022 International Conference on Computational Science (Heisler et al., 2022), we include recent developments for solving nonlinear equations using both automatic and symbolic differentiation, and demonstrate the capability for the Bratu (nonlinear Poisson) equation. 
    more » « less
  4. Boundary effects play an important role in the study of hydrodynamic limits in the Boltzmann theory. Based on a systematic study of the viscous layer equations and the L2 to L∞ framework, we establish the validity of the Hilbert expansion for the Boltzmann equation with specular reflection boundary conditions, which leads to derivations of compressible Euler equations and acoustic equations in half-space. 
    more » « less
  5. In this paper, we study ultra-weak discontinuous Galerkin methods with generalized numerical fluxes for multi-dimensional high order partial differential equations on both unstructured simplex and Cartesian meshes. The equations we consider as examples are the nonlinear convection-diffusion equation and the biharmonic equation. Optimal error estimates are obtained for both equations under certain conditions, and the key step is to carefully design global projections to eliminate numerical errors on the cell interface terms of ultra-weak schemes on general dimensions. The well-posedness and approximation capability of these global projections are obtained for arbitrary order polynomial space based on a wide class of generalized numerical fluxes on regular meshes. These projections can serve as general analytical tools to be naturally applied to a wide class of high order equations. Numerical experiments are conducted to demonstrate these theoretical results.

     
    more » « less