skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


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
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. ; (, Journal of Physics: Conference Series)
    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. ; ; ; ; (, Journal of computational science)
    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
  3. ; ; (, Archive for Rational Mechanics and Analysis)
    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
  4. ; ; ; ; ; ; (, Transport in Porous Media)
    AbstractElectrokinetic in-situ recovery is an alternative to conventional mining, relying on the application of an electric potential to enhance the subsurface flow of ions. Understanding the pore-scale flow and ion transport under electric potential is essential for petrophysical properties estimation and flow behavior characterization. The governing physics of electrokinetic transport is electromigration and electroosmotic flow, which depend on the electric potential gradient, mineral occurrence, domain morphology (tortuosity and porosity, grain size and distribution, etc.), and electrolyte properties (local pH distribution and lixiviant type and concentration, etc.). Herein, mineral occurrence and its associated zeta potential are investigated for EK transport. The new Ek model which is designed to solve the EK flow in complex porous media in a highly parallelizable manner includes three coupled equations: (1) Poisson equation, (2) Nernst–Planck equation, and (3) Navier–Stokes equation. These equations were solved using the lattice Boltzmann method within X-ray computed microtomography images. The proposed model is validated against COMSOL multiphysics in a two-dimensional microchannel in terms of fluid flow behavior when the electrical double layer is both resolvable and unresolvable. A more complex chalcopyrite-silica system is then obtained by micro-CT scanning to evaluate the model performance. The effects of mineral occurrence, zeta potential, and electric potential on the three-dimensional chalcopyrite-silica system were evaluated. Although the positive zeta potential of chalcopyrite can induce a flow of ferric ion counter to the direction of electromigration, the net effect is dependent on the occurrence of chalcopyrite. However, the ion flux induced by electromigration was the dominant transport mechanism, whereas advection induced by electroosmosis made a lower contribution. Overall, a pore-scale EK model is proposed for direct simulation on pore-scale images. The proposed model can be coupled with other geochemical models for full physicochemical transport simulations. Meanwhile, electrokinetic transport shows promise as a human-controllable technique because the electromigration of ions and the applied electric potential can be easily controlled externally. Graphical abstract 
    more » « less
  5. ; (, IMA Journal of Numerical Analysis)
    Abstract In this work we introduce semi-implicit or implicit finite difference schemes for the continuity equation with a gradient flow structure. Examples of such equations include the linear Fokker–Planck equation and the Keller–Segel equations. The two proposed schemes are first-order accurate in time, explicitly solvable, and second-order and fourth-order accurate in space, which are obtained via finite difference implementation of the classical continuous finite element method. The fully discrete schemes are proved to be positivity preserving and energy dissipative: the second-order scheme can achieve so unconditionally while the fourth-order scheme only requires a mild time step and mesh size constraint. In particular, the fourth-order scheme is the first high order spatial discretization that can achieve both positivity and energy decay properties, which is suitable for long time simulation and to obtain accurate steady state solutions. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email