Search NASA⌕ Search

SEARCH · Search NASA

Results for “closure relations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A broadly-applicable unified closure relation for Taylor bubble rise velocity in pipes with stagnant liquid

Taylor bubble velocity in slug flow is a closure relation which significantly affects the prediction of liquid holdup (or void fraction) and pressure gradient in mechanistic models of slug flow for oil and gas pipe applications. In this work, we use a validated Computational Fluid Dynamics (CFD) approach to simulate the motion of Taylor bubbles in pipes; the interface is tracked with a Level-Set method implemented in a commercial code. A large numerical database is generated covering the most ample range of fluid properties and pipe inclination angles explored to date (Eo ∊ [10, 700],Mo ∊ [1 x 10 -6 , 5 x 10 3 , and θ ∊ [0°, 90°]). A unified Taylor bubble rise velocity correlation is extracted from that database. The new correlation predicts the numerical database with 8.6% absolute average relative error and a coefficient of determination R 2 = 0.97 and other available experimental data with 13.0% absolute average relative error and R 2 = 0.84 outperforming existing correlations and models.

Mechanics↗
DOE OSTI2023

Uncertainty Quantification for Multiphase Computational Fluid Dynamics Closure Relations with a Physics-Informed Bayesian Approach

Multiphase Computational Fluid Dynamics (MCFD) based on the two-fluid model is considered a promising tool to model complex two-phase flow systems. MCFD simulation can predict local flow features without resolving interfacial information. As a result, the MCFD solver relies on closure relations to describe the interaction between the two phases. Those empirical or semi-mechanistic closure relations constitute a major source of uncertainty for MCFD predictions. In this paper, we leverage a physics-informed uncertainty quantification (UQ) approach to inversely quantify the closure relations’ model form uncertainty in a physically consistent manner. This proposed approach considers the model form uncertainty terms as stochastic fields that are additive to the closure relation outputs. Combining dimensionality reduction and Gaussian processes, the posterior distribution of the stochastic fields can be effectively quantified within the Bayesian framework with the support of experimental measurements. As this UQ approach is fully integrated into the MCFD solving process, the physical constraints of the system can be naturally preserved in the UQ results. Here, in a case study of adiabatic bubbly flow, we demonstrate that this UQ approach can quantify the model form uncertainty of the MCFD interfacial force closure relations, thus effectively improving the simulation results with relatively sparse data support.

42 ENGINEERING↗
NASA NTRS1979

A possible closure relation for heat transport in the solar wind

The objective of the present paper is to search for an empirical closure relation for solar wind heat transport that applies to a microscopic scale. This task is approached by using the quasi-linear wave-particle formalism proposed by Perkins (1973) as a guide to derive an equation relating the relative drift speed between core-electron and proton populations to local bulk flow conditions. The resulting relationship, containing one free parameter, is found to provide a good characterization of Los Alamos Imp electron data measuring during the period from March 1971 through August 1974. An empirical closure relation is implied by this result because of the observed proportionality between heat flux and relative drift speed.

Feldman, W. C.↗
DOE OSTI2022

The Closure Relations in Optical Afterglow of Gamma-Ray Bursts

Gamma-ray bursts (GRBs) are extremely high-energy events that can be observed at very high redshift. In addition to γ rays, they can emit in X-ray, optical, and sometimes radio wavelengths. Here, following the approach in Srinivasaragavan et al.; Dainotti et al.; and Dainotti et al., we consider 82 GRBs from Dainotti et al. that have been observed in optical wavelengths and fitted with a broken power law. We consider the relations between the spectral and temporal indices (closure relations; CRs) according to the synchrotron forward-shock model evolving in the constant-density interstellar medium (ISM; k = 0) and the stellar wind environment (k = 2) in both slow- and fast-cooling regimes, where the density profile is defined as n ∝ r -k . We find the $ν$ > $max${$ν$ c , $ν$ m } regime is most favored, where ν c and ν m are the cooling and characteristic frequencies, respectively. Finally, we test the 2D Dainotti correlation between the rest-frame end time of the plateau and the luminosity at that time on GRBs that fulfill the most-favored CRs. When we compare the intrinsic scatter σ int of those 2D correlations to the scatter presented in Dainotti et al., we see the scatters of our correlations generally agree with the previous values within 1σ, both before and after correction for selection bias. This new information has helped us to pinpoint subsamples of GRBs with features that could drive the GRB emission mechanism, and eventually allow for GRBs to be used as standard candles.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2019

Closure Relations of Gamma-Ray Bursts in High Energy Emission

The synchrotron external shock model predicts the evolution of the spectral (β) and temporal (α) indices during the gamma-ray burst (GRB) afterglow for different environmental density profiles, electron spectral indices, electron cooling regimes, and regions of the spectrum. We study the relationship between α and β, the so-called “closure relations” with GRBs detected by Fermi Large Area Telescope (Fermi-LAT) from 2008 August to 2018 August. The spectral and temporal indices for the >100 MeV emission from the Fermi-LAT as determined in the Second Fermi-LAT Gamma-Ray Burst Catalog (2FLGC) are used in this work. We select GRBs whose spectral and temporal indices are well constrained (58 long-duration GRBs and 1 short-duration GRBs) and classify each GRB into the best-matched relation. As a result, we found that a number of GRBs require a very small fraction of the total energy density contained in the magnetic field (òB <~10−7). The estimated mean and standard deviation of electron spectral index p are 2.40 and 0.44, respectively. The GRBs satisfying a closure relation of the slow cooling tend to have a softer p value compared to those of the fast cooling. Moreover, the Kolmogorov–Smirnov test of the two p distributions from the fast and slow coolings rejects a hypothesis that the two distributions are drawn from the single reference distribution with a significance of 3.2σ. Lastly, the uniform density medium is preferred over the medium that decreases like the inverse of distance squared for long-duration GRBs.

Gamma-ray bursts↗
DOE OSTI2019

Closure Relations of Gamma-Ray Bursts in High Energy Emission

The synchrotron external shock model predicts the evolution of the spectral (β) and temporal (α) indices during the gamma-ray burst (GRB) afterglow for different environmental density profiles, electron spectral indices, electron cooling regimes, and regions of the spectrum. We review the relationship between α and β, the so-called "closure relations" with GRBs detected by Fermi Large Area Telescope (Fermi-LAT) from 2008 August to 2018 August. The spectral and temporal indices for the >100 MeV emission from the Fermi-LAT as determined in the Second Fermi-LAT Gamma-Ray Burst Catalog (2FLGC) are used in this work. We select GRBs whose spectral and temporal indices are well constrained (58 long-duration GRBs and 1 short-duration GRBs) and classify each GRB into the best-matched relation. As a result, we found that a number of GRBs require a very small fraction of the total energy density contained in the magnetic field (ϵ B ≲ 10 –7 ). The estimated mean and standard deviation of electron spectral index p are 2.40 and 0.44, respectively. The GRBs satisfying a closure relation of the slow cooling tend to have a softer p value compared to those of the fast cooling. Furthermore, the Kolmogorov–Smirnov test of the two p distributions from the fast and slow coolings rejects a hypothesis that the two distributions are drawn from the single reference distribution with a significance of 3.2σ. Lastly, the uniform density medium is preferred over the medium that decreases like the inverse of distance squared for long-duration GRBs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

COVID-19–Related School Closures, United States, July 27, 2020–June 30, 2022

As part of a multiyear project that monitored illness-related school closures, we conducted systematic daily online searches during July 27, 2020–June 30, 2022, to identify public announcements of COVID-19–related school closures (COVID-SCs) in the United States lasting ≥1 day. We explored the temporospatial patterns of COVID-SCs and analyzed associations between COVID-SCs and national COVID-19 surveillance data. COVID-SCs reflected national surveillance data: correlation was highest between COVID-SCs and both new PCR test positivity (correlation coefficient [r] = 0.73, 95% CI 0.56–0.84) and new cases (r = 0.72, 95% CI 0.54–0.83) during 2020–21 and with hospitalization rates among all ages (r = 0.81, 95% CI 0.67–0.89) during 2021–22. The numbers of reactive COVID-SCs during 2020–21 and 2021–22 greatly exceeded previously observed numbers of illness-related reactive school closures in the United States, notably being nearly 5-fold greater than reactive closures observed during the 2009 influenza (H1N1) pandemic.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Second Order Closures for the Radiative Transfer Equation: Some Are Unstable

The largest existing simulations of cosmic reionization model radiative transfer with moment methods that require a closure relation. The two most commonly used closure relations are M1 and OTVET; both close the moment hierarchy at the first moment. We explore the properties of a higher, second-order closure. We show that direct generalizations of M1 and OTVET to one higher order are physically unstable - i.e., the closure equations themselves result in unstable solutions, not just their numerical implementation. In fact, a generalization of OTVET to any order higher than the first one is unstable. We are also able to show that any local (i.e., depending only on the local moments of the radiation field, like M1) second-order closure that depends only on the radiation intensity and radiation flux, but does not explicitly depend on the radiation pressure, is physically unstable. This result restricts the choice of possible second-order closure relations.

Gnedin, Nickolay Y. [Fermilab; Chicago U., Astron.↗
DOE OSTI2019

Multiple Components in the Broadband γ -Ray Emission of the Short GRB 160709A

GRB 160709A is one of the few bright short gamma-ray bursts detected by both the Gamma-ray Burst Monitor and the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. The γ-ray prompt emission of GRB 160709A is adequately fitted by combinations of three distinct components: (i) a nonthermal component described by a power law (PL) with a high-energy exponential cutoff, (ii) a thermal component modeled with a Planck function, and (iii) a second nonthermal component shaped by an additional PL crossing the whole γ-ray spectrum. While the thermal component dominates during ~0.12 s of the main emission episode of GRB 160709A with an unusually high temperature of ~340 keV, the nonthermal components dominate in the early and late time. The thermal component is consistent with the photospheric emission resulting in the following parameters: the size of the central engine, ${R}_{0}={3.8}_{-1.8}^{+5.9}$ × 10 8 cm, the size of the photosphere, ${R}_{\mathrm{ph}}={7.4}_{-1.2}^{+0.8}$ × 10 10 cm, and a bulk Lorentz factor, ${\rm{\Gamma }}={728}_{-93}^{+75}$, assuming a redshift of 1. The slope of the additional PL spectrum stays unchanged throughout the burst duration; however, its flux decreases continuously as a function of time. A standard external shock model has been tested for the additional PL component using the relation between the temporal and spectral indices (the closure relation). Each set of spectral and temporal indices from two energy bands (200 keV–40 MeV and 100 MeV–10 GeV) satisfies a distinct closure relation. From the closure relation test we derived the index for the electron spectral distribution, p = 2.5 ± 0.1. The interaction of the jet with the interstellar environment is preferred over the interaction with the wind medium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Influence of phase connectivity on the relationship among capillary pressure, fluid saturation, and interfacial area in two-fluid-phase porous medium systems

Here, multiphase flow in porous medium systems is typically modeled using continuum mechanical representations at the macroscale in terms of averaged quantities. These models require closure relations to produce solvable forms. One of these required closure relations is an expression relating fluid pressures, fluid saturations, and, in some cases, the interfacial area between the fluid phases, and the Euler characteristic. An unresolved question is whether the inclusion of these additional morphological and topological measures can lead to a non-hysteretic closure relation compared to the hysteretic forms that are used in traditional models, which typically do not include interfacial areas, or the Euler characteristic. We develop a lattice-Boltzmann (LB) simulation approach to investigate the equilibrium states of a two-fluid-phase porous medium system, which include disconnected now- wetting phase features. The proposed approach is applied to a synthetic medium consisting of 1,964 spheres arranged in a random, non-overlapping, close-packed manner, yielding a total of 42,908 different equilibrium points. This information is evaluated using a generalized additive modeling approach to determine if a unique function from this family exists, which can explain the data. The variance of various model estimates is computed, and we conclude that, except for the limiting behavior close to a single fluid regime, capillary pressure can be expressed as a deterministic and non-hysteretic function of fluid saturation, interfacial area between the fluid phases, and the Euler characteristic. This work is unique in the methods employed, the size of the data set, the resolution in space and time, the true equilibrium nature of the data, the parameterizations investigated, and the broad set of functions examined. The conclusion of essentially non-hysteretic behavior provides support for an evolving class of two-fluid-phase flow in porous medium systems models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

COVID-19-associated school closures and related efforts to sustain education and subsidized meal programs, United States, February 18-June 30, 2020

Pre-emptive school closures are frontline community mitigation measures recommended by the US Centers for Disease Control and Prevention (CDC) for implementation during severe pandemics. This study describes the spatiotemporal patterns of publicly announced school closures implemented in response to the coronavirus disease 2019 (COVID-19) pandemic and assesses how public K-12 districts adjusted their methods of education delivery and provision of subsidized meals. During February 18–June 30, 2020, we used daily systematic media searches to identify publicly announced COVID-19–related school closures lasting ≥1 day in the United States (US). We also collected statewide school closure policies from state government websites. Data on distance learning and subsidized meal programs were collected from a stratified sample of 600 school districts. The first COVID-19–associated school closure occurred on February 27, 2020 in Washington state. By March 30, 2020, all but one US public school districts were closed, representing the first-ever nearly synchronous nationwide closure of public K-12 schools in the US. Approximately 100,000 public schools were closed for ≥8 weeks because of COVID-19, affecting >50 million K-12 students. Of 600 districts sampled, the vast majority offered distance learning (91.0%) and continued provision of subsidized meal programs (78.8%) during the closures. Despite the sudden and prolonged nature of COVID-19–associated school closures, schools demonstrated flexibility by implementing distance learning and alternate methods to continue subsidized meal programs.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Computational Fluid Dynamics (CFD) Simulations of Taylor Bubbles in Vertical and Inclined Pipes with Upward and Downward Liquid Flow

Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Engineering↗
DOE OSTI2023

Effects of Different Closure Choices in Core-collapse Supernova Simulations

The two-moment method is widely used to approximate the full neutrino transport equation in core-collapse supernova (CCSN) simulations, and different closures lead to subtle differences in the simulation results. In this paper, we compare the effects of closure choices on various physical quantities in 1D and 2D time-dependent CCSN simulations with our multigroup radiation hydrodynamics code FORNAX. We find that choices of the third-order closure relations influence the time-dependent simulations only slightly. Choices of the second-order closure relation have larger consequences than choices of the third-order closure, but these are still small compared to the remaining variations due to ambiguities in some physical inputs such as the nuclear equation of state. We also find that deviations in Eddington factors are not monotonically related to deviations in physical quantities, which means that simply comparing the Eddington factors does not inform one concerning which closure is better.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Quantum maximum entropy closure for small flavor coherence

Quantum angular moment transport schemes are an important avenue toward describing neutrino flavor mixing phenomena in dense astrophysical environments such as supernovae and merging neutron stars. Successful implementation will require new closure relations that go beyond those used in classical transport. In this paper, we derive the first analytic expression for a quantum M1 closure, valid in the limit of small flavor coherence, based on the maximum entropy principle. We verify that the resulting closure relation has the appropriate limits and characteristic speeds in the diffusive and free-streaming regimes. We then use this new closure in a moment linear stability analysis to search for fast flavor instabilities in a binary neutron star merger simulation and find better results as compared with previously designed, ad hoc , semiclassical closures.

astrophysical & cosmological simulations↗
DOE OSTI2023

Machine learning of hidden variables in multiscale fluid simulation

Solving fluid dynamics equations often requires the use of closure relations that account for missing microphysics. For example, when solving equations related to fluid dynamics for systems with a large Reynolds number, sub-grid effects become important and a turbulence closure is required, and in systems with a large Knudsen number, kinetic effects become important and a kinetic closure is required. By adding an equation governing the growth and transport of the quantity requiring the closure relation, it becomes possible to capture microphysics through the introduction of ‘hidden variables’ that are non-local in space and time. The behavior of the ‘hidden variables’ in response to the fluid conditions can be learned from a higher fidelity or ab-initio model that contains all the microphysics. In our study, a partial differential equation simulator that is end-to-end differentiable is used to train judiciously placed neural networks against ground-truth simulations. We show that this method enables an Euler equation based approach to reproduce non-linear, large Knudsen number plasma physics that can otherwise only be modeled using Boltzmann-like equation simulators such as Vlasov or particle-in-cell modeling.

42 ENGINEERING↗
DOE OSTI2025

Do We Know How to Model Reionization?

I compare the power spectra of the radiation fields from two recent sets of fully-coupled simulations that model cosmic reionization: “Cosmic Reionization On Computers” (CROC) and “Thesan”. While both simulations have similar power spectra of the radiation sources, the power spectra of the photoionization rate are significantly different at the same values of cosmic time or the same values of the mean neutral hydrogen fraction. However, the power spectra of the photoionization rate can be matched at large scales for the two simulations when the matching snapshots are allowed to vary independently. I.e., on large scales, the clustering of the radiation field in two simulations evolves similarly, but the exact timing of this evolution is different in different simulations and is not parameterized by an easily interpretable physical quantity like the mean neutral fraction or the mean free path. On small scales, large differences are present and remain partially unexplained. Both CROC and Thesan use the Variable Eddington Tensor approximation for modeling radiative transfer, but adopt different closure relations (optically thin OTVET versus M1). The role of this key difference is tested by using smaller simulations with a new cosmological simulation code that implements both closure relations in a controlled environment (the same hydro, cooling, and gravity solvers and the star formation recipe). In these controlled tests, both the M1 closure and the OTVET ansatz follow the expected behavior from a simple analytical approximation, demonstrating that the differences in the 2-point function of the radiation field induced by the choice of the Eddington tensor are not dominant.

79 ASTRONOMY AND ASTROPHYSICS↗