Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional generalized linear model”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Isotopic Constraints on the Nature of Primary Precipitates in Archean–Early Paleoproterozoic Iron Formations from Determinations of the Iron Phonon Density of States of Greenalite and 2L- and 6L-Ferrihydrite

Iron formations (IFs) are chemical sedimentary rocks that were widely deposited before the Great Oxidation Event (GOE) around 2.4- 2.2 Ga. It is generally thought that IFs precipitated as hydrated Fe 3+ oxides (HFOs) such as ferrihydrite following surface oxidation of Fe 2+ -rich, anoxic deep waters. This model often implicates biological oxidation and underpins reconstructions of marine nutrient concentrations. However, nanoscale petrography indicates that an Fe 2+ silicate, greenalite, is a common primary mineral in well-preserved IFs, motivating an alternative depositional model of anoxic ferrous silicate precipitation. It is unclear, however, if Fe 2+ -rich silicates can produce the Fe isotopic variations in IFs that are well explained by Fe 2+ oxidation. To address this question, we constrain the equilibrium Fe isotopic ( 56 Fe/ 54 Fe) fractionation of greenalite and ferrihydrite by determining the iron phonon densities of states for those minerals. Here, we use ab initio density functional theory (DFT + U) calculations and nuclear resonant inelastic X-ray scattering spectroscopy to show that ferrous greenalite should be isotopically lighter than ferrihydrite by similar to 1-1.2 parts per thousand at equilibrium, and fractionation should scale linearly with increasing Fe 3+ content in greenalite. By anchoring ferrihydrite-greenalite mineral pair fractionations to published experimental Fe isotopic fractionations between HFOs and aqueous Fe 2+ , we show that ferrous greenalite may produce all but the heaviest pre-GOE Fe isotopic compositions and mixed valence greenalites can produce the entire record. Our results suggest that heavy Fe isotope enrichments alone are not diagnostic of primary IF mineralogies, and ferrihydrite and partially oxidized or even purely ferrous greenalite are all viable primary IF mineralogies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Exploring the dependence of gas cooling and heating functions on the incident radiation field with machine learning

ABSTRACT Gas cooling and heating functions play a crucial role in galaxy formation. But, it is computationally expensive to exactly compute these functions in the presence of an incident radiation field. These computations can be greatly sped up by using interpolation tables of pre-computed values, at the expense of making significant and sometimes even unjustified approximations. Here, we explore the capacity of machine learning to approximate cooling and heating functions with a generalized radiation field. Specifically, we use the machine learning algorithm XGBoost to predict cooling and heating functions calculated with the photoionization code cloudy at fixed metallicity, using different combinations of photoionization rates as features. We perform a constrained quadratic fit in metallicity to enable a fair comparison with traditional interpolation methods at arbitrary metallicity. We consider the relative importance of various photoionization rates through both a principal component analysis (PCA) and calculation of SHapley Additive exPlanation (shap) values for our XGBoost models. We use feature importance information to select different subsets of rates to use in model training. Our XGBoost models outperform a traditional interpolation approach at each fixed metallicity, regardless of feature selection. At arbitrary metallicity, we are able to reduce the frequency of the largest cooling and heating function errors compared to an interpolation table. We find that the primary bottleneck to increasing accuracy lies in accurately capturing the metallicity dependence. This study demonstrates the potential of machine learning methods such as XGBoost to capture the non-linear behaviour of cooling and heating functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Deterministic modeling of hybrid nonlinear effects in epsilon-near-zero thin films

In nonlinear optics, significant effort is concentrated on improving the strength and efficiency of interactions; however, experimentally investigating nonlinear materials is a complex, time-consuming, and costly investment. Moreover, it is often challenging to isolate, study, and optimize material parameters in an experiment due to complexities in the growth process. Recently, epsilon-near-zero materials have received a great deal of attention as promising nonlinear optical materials, but like many up-and-coming materials, the ability to explore and optimize their properties has been challenging. Here, we establish a framework to rapidly evaluate the performance of nonlinear epsilon-near-zero materials for both inter- and intraband effects in silico, requiring only an energy-momentum (E-k) diagram, linear optical properties, and experimental conditions. Measured nonlinear reflection and transmission in gallium-doped zinc oxide films are compared to the numerical framework for both intra- and interband excitation to verify accuracy across wavelength and irradiance while two figures of merit (FoMs) are introduced to quickly evaluate the performance of films without a full numerical framework. This capability is used to predict the performance of highly doped gallium nitride, cadmium oxide, zinc oxide, and indium tin oxide films, and efficient intra- and interband operation conditions are identified. Through this numerical framework and the FoMs, the exploration of unstudied epsilon-near-zero materials is enabled without the need for a nonlinear experiment, thereby accelerating the search for more efficient nonlinear materials and excitation conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

3d-3d correspondence for mapping tori

One of the main challenges in 3d-3d correspondence is that no existent approach offers a complete description of 3d N = 2 SCFT T [ M 3 ] — or, rather, a “collection of SCFTs” as we refer to it in the paper — for all types of 3-manifolds that include, for example, a 3-torus, Brieskorn spheres, and hyperbolic surgeries on knots. The goal of this paper is to overcome this challenge by a more systematic study of 3d-3d correspondence that, first of all, does not rely heavily on any geometric structure on M 3 and, secondly, is not limited to a particular supersymmetric partition function of T [ M 3 ]. In particular, we propose to describe such “collection of SCFTs” in terms of 3d N = 2 gauge theories with “non-linear matter” fields valued in complex group manifolds. As a result, we are able to recover familiar 3-manifold invariants, such as Turaev torsion and WRT invariants, from twisted indices and half-indices of T [ M 3 ], and propose new tools to compute more recent q -series invariants Z * ( M 3 ) in the case of manifolds with b 1 > 0. Although we use genus-1 mapping tori as our “case study,” many results and techniques readily apply to more general 3-manifolds, as we illustrate throughout the paper.

Topological Field Theories↗
DOE OSTI2021

Distinct Mechanisms and Hydricities of Cp*Ir-Based CO 2 Hydrogenation Catalysts in Basic Water

Transition metal-catalyzed reversible hydrogenation of CO 2 to formate in aqueous solutions under ambient conditions is an attractive environmentally friendly strategy for the storage and transportation of H 2 in a liquid chemical form. Here, mechanistic details for the CO 2 hydrogenation by a series of pentamethylcyclopentadienyl iridium(III) (Cp*Ir) complexes with picolinamidate ligands are investigated through density functional theory and kinetic isotope effects (KIEs) calculations. Additionally, adopting a speciation approach at a pH of 8.2, CO 2 hydrogenation mechanisms of 37 distinct protonation states and conformers of Cp*Ir type complexes are investigated. This imposes a challenging test for the computational modeling in terms of providing a consistent correlation with experimental kinetics and KIE studies of multiple catalysts instead of a single complex. Overall, H 2 heterolysis to generate an iridium hydride intermediate was demonstrated to be the rate determining step, which proceeds via distinct pathways depending on the stereoelectronic properties of the ligand. A peculiar iridium dihydride route was also uncovered, which could be optimized to accelerate the H 2 heterolysis. We have further computed the thermodynamic and kinetic hydricities of 74 complexes with [Cp*Ir(L)(H)] q ( q = -1, 0, 1 depending on the ligand charge by deprotonation or protonation events) general formula. Hydricity values do not show any noticeable correlation with the thermodynamics (Δ G ) of [Cp*Ir(L)(HCO 2 – )] formation but exhibit linear correlation with the kinetics (Δ G ‡ ) of electrophilic CO 2 attack to iridium hydride species, especially when the charges and local structural effects of the metal hydrides are considered. Insights gained in this study on (i) the factors determining the preferred pathway for the rate-limiting H 2 heterolysis step, (ii) the correlation between thermodynamic hydricity and the kinetics of CO 2 insertion to iridium hydride species, and (iii) mechanistic analysis via combination of free energy profiles and KIE studies for a series of iridium catalysts will provide guidelines for the design of next-generation catalysts for the reversible H 2 storage through CO 2 utilization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Chemical Framework to Design Linear-like Relaxors toward Capacitive Energy Storage

ABO 3 -type perovskite relaxor ferroelectrics (RFEs) have emerged as the preferred option for dielectric capacitive energy storage. However, the compositional design of RFEs with high energy density and efficiency poses significant challenges owing to the vast compositional space and the absence of general rules. Here, we present an atomic-level chemical framework that captures inherent characteristics in terms of radius and ferroelectric activity of ions. By categorizing A/B-site ions as host framework, rattling, ferroelectrically active, and blocking ions and assembling these four types of ions with specific criteria, linear-like relaxors with weak locally correlated and highly extendable unit-cell polarization vectors can be constructed. As example, we demonstrate two new compositions of Bi 0.5 K 0.5 TiO 3- based and BaTiO 3- based relaxors, showing extremely high recoverable energy densities of 17.3 and 12.1 J cm –3 , respectively, both with a high efficiency of about 90%. Further, the role of different types of ions in forming heterogeneous polar structures is identified through element-specific local structure analysis using neutron total scattering combined with reverse Monte Carlo modeling. Further, our work not only opens up new avenues toward rational compositional design of high energy storage performance lead-free RFEs but also sheds light on atomic-level manipulation of functional properties in compositionally complex ferroelectrics.

25 ENERGY STORAGE↗
DOE OSTI2023

Revising the dynamic energy budget theory with a new reserve mobilization rule and three example applications to bacterial growth

Dynamic energy budget (DEB) theory has been applied to model a wide range of organisms, including microbes. In the standard DEB model, biomass is partitioned into reserve and structural compartments, where reserve biomass is mobilized in a pseudolinear manner (while the reserve biomass density, defined as the ratio between reserve and structural biomass, decays linearly) to drive maintenance and the growth of structural biomass (and extracellular enzyme production if it is considered). However, the linear dynamics of the reserve biomass density makes the standard DEB model incapable of explaining the slowdown of microbial growth at high reserve density that is caused by macromolecular crowding effect which reduces biochemical reaction rates (a typical situation occurs when microbes are experiencing severe moisture stress) and is inconsistent with the observation that intracellular enzymatic reactions generally follow non-linear kinetics. By partitioning biomass into reserve, kinetic, and structural compartments, we show here that the Equilibrium Chemistry Approximation (ECA) kinetics can be used to represent enzymatically catalyzed reserve biomass mobilization that can then drive the kinetic and structural biomass synthesis. This revised DEB model better represents the tradeoff in ribosome allocation for structural growth and internal enzyme production, is structurally compatible with metabolic models of cell individuals, and includes the standard DEB model and the popular compromise model as special cases for representing population growth. We then applied the revised DEB model to interpret components of bacterial respiration, their dependence on substrate availability, and emergent microbial carbon use efficiency dynamics for an exponentially growing population. We found that the revised DEB model enables a better understanding of bacterial substrates use (carbon in our examples) than that can be derived from a few other models in the literature. In particular, the revised DEB model explains why carbon use efficiency may first increase, then plateau, and finally decrease with growth rate (and substrate uptake rate), as a function of proteomics. Additionally, the revised DEB model explains why the kinetic biomass compartment needs to be divided to reasonably incorporate proteomic control of microbial growth.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Linking the Dynamic Response of the Carbon Dioxide-Concentrating Mechanism to Carbon Assimilation Behavior in Fremyella diplosiphon

ABSTRACT Cyanobacteria use a carbon dioxide (CO 2 )-concentrating mechanism (CCM) that enhances their carbon fixation efficiency and is regulated by many environmental factors that impact photosynthesis, including carbon availability, light levels, and nutrient access. Efforts to connect the regulation of the CCM by these factors to functional effects on carbon assimilation rates have been complicated by the aqueous nature of cyanobacteria. Here, we describe the use of cyanobacteria in a semiwet state on glass fiber filtration discs—cyanobacterial discs—to establish dynamic carbon assimilation behavior using gas exchange analysis. In combination with quantitative PCR (qPCR) and transmission electron microscopy (TEM) analyses, we linked the regulation of CCM components to corresponding carbon assimilation behavior in the freshwater, filamentous cyanobacterium Fremyella diplosiphon . Inorganic carbon (C i ) levels, light quantity, and light quality have all been shown to influence carbon assimilation behavior in F. diplosiphon . Our results suggest a biphasic model of cyanobacterial carbon fixation. While behavior at low levels of CO 2 is driven mainly by the C i uptake ability of the cyanobacterium, at higher CO 2 levels, carbon assimilation behavior is multifaceted and depends on C i availability, carboxysome morphology, linear electron flow, and cell shape. Carbon response curves (CRCs) generated via gas exchange analysis enable rapid examination of CO 2 assimilation behavior in cyanobacteria and can be used for cells grown under distinct conditions to provide insight into how CO 2 assimilation correlates with the regulation of critical cellular functions, such as the environmental control of the CCM and downstream photosynthetic capacity. IMPORTANCE Environmental regulation of photosynthesis in cyanobacteria enhances organismal fitness, light capture, and associated carbon fixation under dynamic conditions. Concentration of carbon dioxide (CO 2 ) near the carbon-fixing enzyme RubisCO occurs via the CO 2 -concentrating mechanism (CCM). The CCM is also tuned in response to carbon availability, light quality or levels, or nutrient access—cues that also impact photosynthesis. We adapted dynamic gas exchange methods generally used with plants to investigate environmental regulation of the CCM and carbon fixation capacity using glass fiber-filtered cells of the cyanobacterium Fremyella diplosiphon . We describe a breakthrough in measuring real-time carbon uptake and associated assimilation capacity for cells grown in distinct conditions (i.e., light quality, light quantity, or carbon status). These measurements demonstrate that the CCM modulates carbon uptake and assimilation under low-C i conditions and that light-dependent regulation of pigmentation, cell shape, and downstream stages of carbon fixation are critical for tuning carbon uptake and assimilation.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A hybrid calibration approach to Hertz-type contact parameters for discrete element models

This study aims at providing a hybrid calibration framework to estimate Hertz-type contact parameters (particle-scale shear modulus and Poisson ratio) for both two-dimensional and three-dimensional discrete element modelling (DEM). On the basis of statistically isotropic granular packings, a set of analytical formulae between macroscopic material parameters (Young modulus and Poisson ratio) and particle-scale Hertz-type contact parameters for granular systems are derived under small-strain isotropic stress conditions. However, the derived analytical solutions are only estimated values for general models. By viewing each DEM modelling as an implicit mathematical function taking the particle-level parameters as independent variables and employing the derived analytical solutions as the initial input parameters, an automatic iterative scheme is proposed to obtain the calibrated parameters with higher accuracies. Considering highly nonlinear features and discontinuities of the macro-micro relationship in Hertz-based discrete element models, the adaptive moment estimation algorithm is adopted in this study because of its capacity of dealing with noise gradients of cost functions. Here, the proposed method is validated with several numerical cases including randomly distributed monodisperse and polydisperse packings. Noticeable improvements in terms of calibration efficiency and accuracy have been made.

Constitutive law↗
DOE OSTI2023

Noble Gas Transport in the MSRE

This study explores the relationship between the physicochemical parameter known as Henry’s gas constant and gas transport across a two-layer film interface. The investigation utilized the Gibbs free energy, incorporating surface and volume terms to elucidate trends in enthalpy and entropy. Notably, our findings align with experimental data and offer predictive insights into the Henry’s gas constants for helium and krypton, which hold significance for future experiments and theoretical developments. Furthermore, this study enhances the Gibbs free energy theory pertaining to the liquid–gas interface. It underscores the substantial contribution of noble gases in this region to volumetric energy as temperature increases. Additionally, we employed Monte Carlo simulations to analyze the effective thermal neutron multiplication factor, denoted as k eff . Our analysis reveals a linear correlation between graphite density and uniform density as a function of temperature. For the 1D Molten Salt Reactor Experiment (MSRE) system, we employed the Mole code to conduct heat and mass transfer calculations. These computations enable us to ascertain the distribution of fuel temperature based on coefficients and thermal properties. We also studied delayed neutron precursors during fuel cycling, taking into account the drift of cycling fuel through Mole–Griffin coupling. Our model represents k eff and β eff across various volume flow rates and salinity compositions. Finally, this study leveraged xenon-135 for continuous on-line monitoring of fuel salts to investigate the impact of steady-state xenon-135 on the MSRE and to better understand its distribution. These efforts build upon previous research related to removal processes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Multimode two-dimensional vibronic spectroscopy. I. Orientational response and polarization-selectivity

Two-dimensional Electronic–Vibrational (2D EV) spectroscopy and two-dimensional Vibrational–Electronic (2D VE) spectroscopy are among the newest additions to the coherent multidimensional spectroscopy toolbox, and they are directly sensitive to vibronic couplings. In this first of two papers, the complete orientational response functions are developed for a model system consisting of two coupled anharmonic oscillators and two electronic states in order to simulate polarization-selective 2D EV and 2D VE spectra with arbitrary combinations of linearly polarized electric fields. Here, we propose analytical methods to isolate desired signals within complicated spectra and to extract the relative orientation between vibrational and vibronic dipole moments of the model system using combinations of polarization-selective 2D EV and 2D VE spectral features. Time-dependent peak amplitudes of coherence peaks are also discussed as means for isolating desired signals within the time-domain. This paper serves as a field guide for using polarization-selective 2D EV and 2D VE spectroscopies to map coupled vibronic coordinates on the molecular frame.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

MechElastic: A Python library for analysis of mechanical and elastic properties of bulk and 2D materials

We report the MechElastic Python package evaluates the mechanical and elastic properties of bulk and 2D materials using the elastic coefficient matrix ( C ij ) obtained from any ab-initio density-functional theory (DFT) code. The current version of this package reads the output of VASP, ABINIT, and Quantum Espresso codes (but it can be easily generalized to any other DFT code) and performs the appropriate post-processing of elastic constants as per the requirement of the user. This program can also detect the input structure's crystal symmetry and test the mechanical stability of all crystal classes using the Born-Huang criteria. Various useful material-specific properties such as elastic moduli, longitudinal and transverse elastic wave velocities, Debye temperature, elastic anisotropy, 2D layer modulus, hardness, Pugh's ratio, Cauchy's pressure, Kleinman's parameter, and Lame's coefficients, can be estimated using this program. Another existing feature of this program is to employ the ELATE package (2016) [29] and plot the spatial variation of several elastic properties such as Poisson's ratio, linear compressibility, shear modulus, and Young's modulus in three dimensions. Further, the MechElastic package can plot the equation of state (EOS) curves for energy and pressure for a variety of EOS models such as Murnaghan, Birch, Birch-Murnaghan, and Vinet, by reading the inputted energy/pressure versus volume data obtained via numerical calculations or experiments. This package is particularly useful for the high-throughput analysis of elastic and mechanical properties of materials.

2D materials↗
DOE OSTI2023

Probing the rapid formation of black holes and their Galaxy hosts in QSOs

Using the modelling code x-cigale, we reproduced the spectral energy distributions (SEDs) of 1359 SDSS quasi-stellar objects (QSOs) within the redshift range 0 < z < 4, for which we have NIR/MIR fluxes with the highest quality and spectral data characterizing their supermassive black hole (SMBHs). Consistent with a rapid formation of the host galaxies, the star formation histories (SFHs) have small e-folding, at most 750 Myr using an SFH function for spiral or 1000 Myr using one for elliptical. Above z ∼ 1.6, the two solutions are degenerate, the SEDs being dominated by the active galactic nucleus (AGN) continuum and high star formation rates (SFRs), typical of starburst galaxies, while at lower redshifts the starburst nature of the host, independent from its morphology, is better reproduced by a spiral SFH. In general, the SFR increases with the redshift, the mass of the bulge, the AGN luminosity, and Eddington ratio, suggesting there is no evidence of AGN quenching of star formation. Comparing the specific BH accretion rate (BHAR) with specific SFR, all the QSOs at any redshift trace a linear sequence below the Eddington luminosity, in parallel and above the one-to-one relation, implying that QSOs are in a special phase of evolution during which the growth in mass of their SMBH is more rapid than the growth in mass of their galaxy hosts. This particular phase is consistent with a scenario where the galaxy hosts of QSOs in the past grew in mass more rapidly than their SMBHs, suggesting that a high star formation efficiency during their formation was responsible in limiting their masses.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Heterotic instantons for monad and extension bundles

We consider non-perturbative superpotentials from world-sheet instantons wrapped on holomorphic genus zero curves in heterotic string theory. These superpotential contributions feature prominently in moduli stabilization and large field axion inflation, which makes their presence or absence, as well as their functional dependence on moduli, an important issue. We develop geometric methods to compute the instanton superpotentials for heterotic string theory with monad and extension bundles. Using our methods, we find a variety of examples with a non-vanishing superpotential. In view of standard vanishing theorems, we speculate that these results are likely to be attributed to the non-compactness of the instanton moduli space. We test this proposal, for the case of monad bundles, by considering gauged linear sigma models where compactness of the instanton moduli space can be explicitly checked. In all such cases, we find that the geometric results are consistent with the vanishing theorems. Surprisingly, linearly dependent Pfaffians even arise for cases with a non-compact instanton moduli space. This suggests some gauged linear sigma models with a non-compact instanton moduli space may still have a vanishing instanton superpotential.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Multiscale simulations for multi-continuum Richards equations

In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization’s effective coefficients are computed through solving the arising cell problems. Furthermore, to tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Flow dynamics and heat transfer in simplified battery energy storage systems with heated battery modules

Large-scale energy storage systems (ESSs) composed of batteries show promise in addressing current energy challenges, but dissipation of generated heat is important. Here, this paper focuses on buoyant convective flows in simplified ESS battery racks. Natural convection is not generally the primary cooling strategy but can be important in abnormal scenarios where there is module overheat or potentially thermal runaway. We use computational fluid dynamics to investigate the flow dynamics and heat transfer mechanisms in a simplified parameterized rack design. Despite its simplicity, this configuration produces many of the relevant features expected in real ESSs without details of module geometry or hardware, allowing broad conclusions independent of manufacture-specific designs. We start by providing visualizations of the flowfield and measurements of entrainment, heat flux, and pressure. To characterize the dependence on the system parameters, we develop an integral-scale analysis of the average temperature equation to highlight the dominant source terms. We use results from this analysis to derive a steady network model composed of simple algebraic expressions to provide first-order predictions of entrainment through the rack. The network model leads to a linear scaling of the Reynolds number based on convective mass flux with respect to the Grashof number based on the heat source. We deduce empirical relationships that relate the heat exchanged between modules using a surface-averaged Nusselt number as a function of the local Reynolds and Rayleigh numbers. Lastly, we investigate how space between the modules and rack in the spanwise direction creates flow bypass, resulting in different flow pathways.

Battery thermal management↗
DOE OSTI2023

2D k -th nearest neighbour statistics: a highly informative probe of galaxy clustering

ABSTRACT Beyond standard summary statistics are necessary to summarize the rich information on non-linear scales in the era of precision galaxy clustering measurements. For the first time, we introduce the 2D k-th nearest neighbour (kNN) statistics as a summary statistic for discrete galaxy fields. This is a direct generalization of the standard 1D kNN by disentangling the projected galaxy distribution from the redshift-space distortion signature along the line-of-sight. We further introduce two different flavours of 2D kNNs that trace different aspects of the galaxy field: the standard flavour which tabulates the distances between galaxies and random query points, and a ‘DD’ flavour that tabulates the distances between galaxies and galaxies. We showcase the 2D kNNs’ strong constraining power both through theoretical arguments and by testing on realistic galaxy mocks. Theoretically, we show that 2D kNNs are computationally efficient and directly generate other statistics such as the popular two-point correlation function (2PCF), voids probability function, and counts-in-cell statistics. In a more practical test, we apply the 2D kNN statistics to simulated galaxy mocks that fold in a large range of observational realism and recover parameters of the underlying extended halo occupation distribution (HOD) model that includes velocity bias and galaxy assembly bias. We find unbiased and significantly tighter constraints on all aspects of the HOD model with the 2D kNNs, both compared to the standard 1D kNN, and the classical redshift-space 2PCF.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Effect of the Prior and the Experimental Design on the Inference of the Precision Matrix in Gaussian Chain Graph Models

Here, we investigate whether (and how) experimental design could aid in the estimation of the precision matrix in a Gaussian chain graph model, especially the interplay between the design, the effect of the experiment and prior knowledge about the effect. Estimation of the precision matrix is a fundamental task to infer biological graphical structures like microbial networks. We compare the marginal posterior precision of the precision matrix under four priors: flat, conjugate Normal-Wishart, Normal-MGIG and a general independent. Under the flat and conjugate priors, the Laplace-approximated posterior precision is not a function of the design matrix rendering useless any efforts to find an optimal experimental design to infer the precision matrix. In contrast, the Normal-MGIG and general independent priors do allow for the search of optimal experimental designs, yet there is a sharp upper bound on the information that can be extracted from a given experiment. We confirm our theoretical findings via a simulation study comparing (i) the KL divergence between prior and posterior and (ii) the Stein’s loss difference of MAPs between random and no experiment. Our findings provide practical advice for domain scientists conducting experiments to better infer the precision matrix as a representation of a biological network.

54 ENVIRONMENTAL SCIENCES↗