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. This work aims at comparing the accuracy and overall performance of a low-Mach CFD solver and a fully-compressible CFD solver for direct numerical simulation (DNS) of nonequilibrium plasma assisted ignition (PAI) using a phenomenological model described in Castela et al. [1]. The phenomenological model describes the impact of nanosecond pulsed plasma discharges by introducing source terms in the reacting flow equations, instead of solving the detailed plasma kinetics at every time step of the discharge. Ultra-fast gas heating and dissociation ofO2 are attributed to the electronic excitation ofN2 and the subsequent quenching to ground state. This process is highly exothermic, and is responsible for dissociation of O2 to form O radicals; both of which promote faster ignition. Another relatively slower process of gas heating associated with vibrational-to-translational relaxation is also accounted for, by solving an additional vibrational energy transport equation. A fully-compressible CFD solver for high Mach (M>0.2) reacting flows, developed by extending the default rhoCentralFoam solver in OpenFOAM, is used to perform DNS of PAI in a 2D domain representing a cross section of a pin-to-pin plasma discharge configuration. The same case is also simulated using a low-Mach, pressure-based CFD solver, built by extending the default reactingFoam solver. The lack of flow or wave dominated transport after the plasma-induced weak shock wave leaves the domain causes inaccurate computation of all the transport variables, with a rather small time step dictated by the CFL condition, with the fully-compressible solver. These issues are not encountered in the low-Mach solver. Finally, the low-Mach solver is used to perform DNS of PAI in lean, premixed, isotropic turbulent mixtures of CH4-air at two different Reynolds numbers of 44 and 395. Local convection of the radicals and vibrational energy from the discharge domain, and straining of the high temperature reaction zones resulted in slower ignition of the case with the higher Re. A cascade effect of temperature reduction in the more turbulent case also resulted in a five - six times smaller value of the vibrational to translational gas heating source term, which further inhibited ignition. Two pulses were sufficient for ignition of the Re = 44 case, whereas three pulses were required for the Re = 395 case; consistent with the results of Ref. [1]. 
    more » « less
  2. 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
  3. The presence of non-local interactions and intermittent signals in the homogeneous isotropic turbulence grant multi-point statistical functions a key role in formulating a new generation of large-eddy simulation (LES) models of higher fidelity. We establish a tempered fractional-order modelling framework for developing non-local LES subgrid-scale models, starting from the kinetic transport. We employ a tempered Lévy-stable distribution to represent the source of turbulent effects at the kinetic level, and we rigorously show that the corresponding turbulence closure term emerges as the tempered fractional Laplacian, $(\varDelta +\lambda )^{\alpha } (\cdot )$ , for $\alpha \in (0,1)$ , $\alpha \neq \frac {1}{2}$ and $\lambda >0$ in the filtered Navier–Stokes equations. Moreover, we prove the frame invariant properties of the proposed model, complying with the subgrid-scale stresses. To characterize the optimum values of model parameters and infer the enhanced efficiency of the tempered fractional subgrid-scale model, we develop a robust algorithm, involving two-point structure functions and conventional correlation coefficients. In an a priori statistical study, we evaluate the capabilities of the developed model in fulfilling the closed essential requirements, obtained for a weaker sense of the ideal LES model (Meneveau, Phys. Fluids , vol. 6, issue 2, 1994, pp. 815–833). Finally, the model undergoes the a posteriori analysis to ensure the numerical stability and pragmatic efficiency of the model. 
    more » « less
  4. Abstract

    We present a theory based on the conventional two-term (i.e. Lorentzian) approximation to the exact solution of the Boltzmann equation in non-magnetized weakly ionized plasma to efficiently obtain the electron rate and transport coefficients in a magnetized plasma for an arbitrary magnitude and direction of applied electric fieldEand magnetic fieldB. The proposed transcendental method does not require the two-term solution of the Boltzmann equation in magnetized plasma, based on which the transport parameters vary as a function of the reduced electric fieldE/N, reduced electron cyclotron frequencyωce/N, and angleE,BbetweenEandBvectors, whereNis the density of neutrals. Comparisons between the coefficients derived from BOLSIG+’s solution (obtained via the two-term expansion whenB0) and coefficients of the presented method are illustrated for air, a mixture of molecular hydrogen (H2) and helium (He) representing the giant gas planets of the Solar System, and pure carbon dioxide (CO2). The new approach may be used in the modeling of magnetized plasma encountered in the context of transient luminous events, e.g. sprite streamers in the atmosphere of Earth and Jupiter, in modeling the propagation of lightning’s electromagnetic pulses in Earth’s ionosphere, and in various laboratory and industrial applications of nonthermal plasmas.

     
    more » « less
  5. Background: When phenotypic characters are described in the literature, they may be constrained or clarified with additional information such as the location or degree of expression, these terms are called “modifiers”. With effort underway to convert narrative character descriptions to computable data, ontologies for such modifiers are needed. Such ontologies can also be used to guide term usage in future publications. Spatial and method modifiers are the subjects of ontologies that already have been developed or are under development. In this work, frequency (e.g., rarely, usually), certainty (e.g., probably, definitely), degree (e.g., slightly, extremely), and coverage modifiers (e.g., sparsely, entirely) are collected, reviewed, and used to create two modifier ontologies with different design considerations. The basic goal is to express the sequential relationships within a type of modifiers, for example, usually is more frequent than rarely, in order to allow data annotated with ontology terms to be classified accordingly. Method: Two designs are proposed for the ontology, both using the list pattern: a closed ordered list (i.e., five-bin design) and an open ordered list design. The five-bin design puts the modifier terms into a set of 5 fixed bins with interval object properties, for example, one_level_more/less_frequently_than, where new terms can only be added as synonyms to existing classes. The open list approach starts with 5 bins, but supports the extensibility of the list via ordinal properties, for example, more/less_frequently_than, allowing new terms to be inserted as a new class anywhere in the list. The consequences of the different design decisions are discussed in the paper. CharaParser was used to extract modifiers from plant, ant, and other taxonomic descriptions. After a manual screening, 130 modifier words were selected as the candidate terms for the modifier ontologies. Four curators/experts (three biologists and one information scientist specialized in biosemantics) reviewed and categorized the terms into 20 bins using the Ontology Term Organizer (OTO) (http://biosemantics.arizona.edu/OTO). Inter-curator variations were reviewed and expressed in the final ontologies. Results: Frequency, certainty, degree, and coverage terms with complete agreement among all curators were used as class labels or exact synonyms. Terms with different interpretations were either excluded or included using “broader synonym” or “not recommended” annotation properties. These annotations explicitly allow for the user to be aware of the semantic ambiguity associated with the terms and whether they should be used with caution or avoided. Expert categorization results showed that 16 out of 20 bins contained terms with full agreements, suggesting differentiating the modifiers into 5 levels/bins balances the need to differentiate modifiers and the need for the ontology to reflect user consensus. Two ontologies, developed using the Protege ontology editor, are made available as OWL files and can be downloaded from https://github.com/biosemantics/ontologies. Contribution: We built the first two modifier ontologies following a consensus-based approach with terms commonly used in taxonomic literature. The five-bin ontology has been used in the Explorer of Taxon Concepts web toolkit to compute the similarity between characters extracted from literature to facilitate taxon concepts alignments. The two ontologies will also be used in an ontology-informed authoring tool for taxonomists to facilitate consistency in modifier term usage. 
    more » « less