SEARCH · Engineering Papers
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.
This package implements Bayesian Neural Network (BNN) based surrogate models for multi-modal data. These models are fit using a custom Variational Bayes procedure leveraging conjugate posteriors in the last layer of the network, offering more accurate predictions with better-calibrated uncertainty quantification than mean-field Variational Bayes.
The software contains objects that are building blocks for PV modeling and interactive data fitting based on: Optoelectronic models for tandem/multijunction solar cells including resistive and luminescent coupling; simulation of modules composed of 2T, 3T, and 4T tandem solar cells; and energy yield analysis of PV systems composed of tandem solar cells.
Background: Rapid diagnostic tests (RDTs) are widely used for malaria diagnosis of both symptomatic and asymptomatic infections. Although RDTs are a reliable and practical diagnostic tool, the sensitivity of histidine-rich protein 2 (HRP2)-based RDTs can be reduced if pfhrp2 or pfhrp3 (pfhrp2/3) gene deletions exist in the Plasmodium falciparum parasite population. This study evaluated dried blood spot (DBS) samples collected from a national household survey to investigate the presence of pfhrp2/3 deletions and the performance of the RDT used in the cross-sectional survey in a low transmission setting. Methods: The 2015 Ethiopia Malaria Indicator Survey tested household members by RDT and collected DBS samples. DBS (n = 2648) from three regions in northern Ethiopia were tested by multiplex bead-based antigen detection assay after completion of the survey. The multiplex assay detected pan-Plasmodium lactate dehydrogenase (LDH), pAldolase, and HRP2 antigens in samples. Samples suspected for pfhrp2/3 gene deletions (pLDH and/or pAldolase positive but low or absent HRP2) were further investigated by molecular assays for gene deletions. Antigen results were also compared to each individual’s RDT results. Dose–response logistic regression models were fit to estimate RDT level of detection (LOD) antigen concentrations at which 50, 75, 90, and 95% of the RDTs returned a positive result during this survey. Results: Out of 2,648 samples assayed, 29 were positive for pLDH or pAldolase antigens but low or absent for HRP2 signal, and 15 of these samples (51.7%) were successfully genotyped for pfhrp2/3. Of these 15 P. falciparum infections, eight showed single deletions in pfhrp3, one showed a single pfhrp2 deletion, and six were pfhrp2/3 double-deletions. Six pfhrp2 deletions were observed in Tigray and one in Amhara. Twenty-five were positive for HRP2 by the survey RDT while the more sensitive bead assay detected 30 HRP2-positive samples. A lower concentration of HRP2 antigen generated a positive test result by RDT compared to pLDH (95% LOD: 16.9 ng/mL vs. 319.2 ng/mL, respectively). Conclusions: There is evidence of dual pfhrp2/3 gene deletions in the Tigray and Amhara regions of Ethiopia in 2015. As the prevalence of malaria was very low (< 2%), it is difficult to make strong conclusions on RDT performance, but these results challenge the utility of biomarkers in household surveys in very low transmission settings.
This data upload includes a compilation of experiments for quantifying uranium adsorption onto quartz. The provided .csv file has been compiled from the literature in a findable, accessible, interoperable, reusable (FAIR) data format. This was accomplished using the Lawrence Livermore National Laboratory Surface Complexation Database Converter (SCDC) code written in the R programming language (free licensing available at https://ipo.llnl.gov/technologies/software/llnl-surface-complexation-database-converter-scdc). This constitutes all current data compiled on uranium-quartz interactions in the L-SCIE (LLNL Surface Complexation/Ion Exchange) database (as of 02/28/2022). This data was used to develop surface complexation models that fit the global community dataset (https://doi.org/10.1021/acs.est.1c07109). The FAIR-formatted dataset also enables the implementation of alternative machine-learning approaches that can be explored in the future.
The development of insulation materials with low effective thermal conductivity is essential for energy savings in various applications, including buildings, food services, pipe insulation, and refrigeration. Such materials can be developed by using micro- or nanoporous structures, as well as low-thermal conductivity gases and/or reduced pressure. A variety of effective medium approximation (EMA) models have been developed to study thermal transport through solid in porous structures. However, in many cases, the impacts of porosity on solid conduction are not well predicted by EMA models because of the assumptions made for simplification. Furthermore, the results can vary by a factor of up to 1.5, depending on the morphology of the pores. Hence, proper guidance is needed to choose the appropriate EMA model for a given morphology. This work presents a finite element method study using COMSOL Multiphysics software for various geometries, including hollow cubes, overlapping hollow spheres (normally stacked and tightly packed), and Voronoi structures, and compares the results to those obtained from various EMA models. Simulation results for the hollow cubes and Voronoi structures matched well with the Maxwell-Eucken and Russell model. The results for hollow spheres closely matched the Bauer model and the Glicksman model with fitting parameters. This work provides guidance on analyzing and designing insulation materials for energy savings in the future.
A new copper equation of state is developed utilizing the available experimental data in addition to recent theoretical calculations. Semi-empirical models are fit to the data and the results are tabulated in the SNL SESAME format. Comparison to other copper EOS tables are given, along with recommendations of which tables provide the best accuracy.
The carrier dynamics in InGaN/GaN light-emitting diodes (LEDs) are directly tied to their efficiency and maximum modulation speed, which are important metrics for solid-state lighting, displays, and optical communication. In this work, we measure the carrier dynamics of a variety of InGaN/GaN LEDs using small-signal electroluminescence methods [1]. A rate equation approach and associated small-signal circuit are used to model carrier injection, recombination in the active region, recombination in the cladding regions, and carrier escape. The model is fit to the measured optical frequency response (S21) and input impedance (S11) of the LEDs to extract the various carrier lifetimes, the carrier density, and the radiative and non-radiative recombination rates. We specifically study planar nonpolar and semipolar LEDs, which show record-high modulation speeds for III-nitride LEDs and present the modulation characteristics of core-shell nanowire-based LEDs. The planar nonpolar m-plane ($101\bar{0}$) micro-LEDs achieve a record-high -3dB modulation bandwidth for a III-nitride LED of 1.5 GHz [2]. The -3dB response of an electrically injected nanowire-based micro-LED with nonpolar facets is also reported, showing a -3dB bandwidth of 1.2 GHz [3]. The high speed is attributed to the shorter carrier lifetime associated with the nonpolar orientation. We also study the carrier dynamics in semipolar ($20\bar{2}\bar{1}$) LEDs for various temperatures [4]. Finally, we present carrier dynamics measurements on commercial-grade c-plane epitaxy for various active region designs, including a wavelength series and a growth quality series. The wavelength series offers insight into the contributions of the quantum confined Stark effect (QCSE) and InGaN material quality on the green gap [5]. The growth quality series investigates the role of non-radiative centers on the LED performance. Extraction of the carrier dynamics using small-signal electroluminescence offers insight into the factors limiting the efficiency and high-speed performance of III-nitride emitters and can be leveraged to ultimately improve the devices.
As a result of a multi-year NA-22 project, CGMF was integrated into MCNP6.2 and publicly released. CGMF was open-sourced and publicly released and MCNP6.3 was updated to include the latest version and is in the process of being publicly released. Current and future plans include global optimization and uncertainty quantification within CGMF, model parameter fitting such that CGMF may be used in ENDF/B evaluations, and improving both standalone and MCNP-integrated CGM (non-fission) simulations.
The carrier dynamics in InGaN/GaN LEDs are directly tied to their efficiency and maximum modulation speed, which are important metrics for solid-state lighting, displays, and optical communication. We measure the carrier dynamics of nonpolar, semipolar, and commercial c-plane InGaN/GaN LEDs using small-signal electroluminescence methods. Rate equations and a small-signal circuit are used to model the carrier dynamics. The model is fit to the optical frequency response and input impedance of the LEDs to extract the carrier lifetimes and the recombination rates. The results offer insight into the underlying causes of efficiency droop and the green gap, and inform device design strategies.
Research Element 1: Denitrification is a difficult flux to measure when fluxes are low and when the ecosystem is difficult to sample. These attributes describe Nyack floodplain, a large gravel-bedded aquifer adjacent to the Middle Fork Flathead River. Given that denitrification is simply respiration with nitrate as a terminal electron acceptor, we reasoned that N 2 production should negatively covary with oxygen respiration so long that oxygen is available. We sampled 16 wells 8 times throughout a year and analyzed these water samples on a membrane-inlet mass spectrometer for N 2 :Ar and O 2 . We fit multilevel models to N 2 :Ar as a function of dissolved oxygen. Within wells, there was no relationship between O 2 and N 2 :Ar. Among wells N 2 :Ar increased as O 2 declined. Based on the slope of this relationship, ~ -0.07, and estimated oxygen respiration in Nyack derived from prior studies, a possible range for denitrification is 0-8 mmol m -2 d -1 corresponding to a turnover rate for nitrate of 0-0.2/d. This estimated turnover rate is in the lower 20% of rates measured in the LINX 2 study, yet shows that even in a low nitrate (0.005 mmol/L) aquifer denitrification can be a substantial sink for dissolved N. Research Element 3 and 4: We developed a new theoretical approach to describe the relationship between water inflow, water outflow, and water age in porous media control volumes under dynamic flow conditions. Our approach had three distinct steps. First we described the application of well-established equations borrowed from Chemical Engineering to describe water age distributions in the hyporheic zone under steady state conditions. We published the associated results in the scientific journal PLOS ONE. Second, we developed a new theory ("flow spans") to describe how predicted steady-state water age distributions will respond to time varying imbalances in hyporheic discharge and recharge. Finally, we developed a proof-of-concept simulation model based on the new theory and used the model to examine the effects of flow regulation on a simulated river with characteristics representative of western U.S. montane rivers with expansive hyporheic zones.
We present galaxy-galaxy lensing measurements of a sample of low surface brightness galaxies (LSBGs) drawn from the Dark Energy Survey Year 3 (Y3) data. LSBGs are diffuse galaxies with a surface brightness dimmer than the ambient night sky. Given their faintness, the use of standard observational techniques proves challenging. Weak gravitational lensing probes both the baryonic and dark matter content of galaxies, rendering it a powerful technique to estimate LSBG masses. The LSBG lens sample consists of 23, 790 total extended galaxies separated into red and blue color types at \(g-i\ge 0.60\) and \(g-i\le 0.60\), respectively. We use the Y3 \sc{metacalibration} shape catalog as the source sample, with a number density of \(5.59 \mathrm{gal}/\mathrm{arcmin}^2\). We measure the tangential shear around the lens galaxies across angular scales of \(0.25-400\) arcmin and find a signal-to-noise of 6.67 for red galaxies, 2.17 for blue galaxies, and 5.30 for the total sample. We fit a model built from two NFW profiles corresponding to the LSBG dark matter subhalo and host halo to the red LSBG shear measurements with an MCMC. We estimate the host halo mass at \(7.3 ^{+2.0}_{-1.6}\times 10^{12} M_\mathrm{\odot}\). We place a 95\% upper bound on the subhalo mass at \(3\times 10^{11} M_\mathrm{\odot}\). We utilize the lens sample photometry to obtain an estimate of the red LSBG stellar mass distribution. We compare the ratio between the stellar mass and the subhalo mass to the parameterized, satellite-specific SHMR. This work represents the first example of an attempted constraint on the masses of LSBGs using weak gravitational lensing.
Cryogenic Process Design Kits (PDKs) are an indispensable tool in the design of complex integrated circuits across a wide spectrum of applications, from noble element detectors to Quantum Information Science, Superconducting Nanowire Single Photon Detectors (SNSPDs), and precision atomic clocks. The development of PDK-compatible SPICE models is a complex endeavor requiring test structures, measurements, model extraction and fitting. We will present the cryogenic modeling and development of a cryo-PDK for a 22nm FDSOI CMOS process for operation at 3.8 Kelvin.
Background: In case-cohort studies with binary outcomes, ordinary logistic regression analyses have been widely used because of their computational simplicity. However, the resultant odds ratio estimates cannot be interpreted as relative risk measures unless the event rate is low. The risk ratio and risk difference are more favorable outcome measures that are directly interpreted as effect measures without the rare disease assumption. Methods: We provide pseudo-Poisson and pseudo-normal linear regression methods for estimating risk ratios and risk differences in analyses of case-cohort studies. These multivariate regression models are fitted by weighting the inverses of sampling probabilities. Also, the precisions of the risk ratio and risk difference estimators can be improved using auxiliary variable information, specifically by adapting the calibrated or estimated weights, which are readily measured on all samples from the whole cohort. Finally, we provide computational code in R (R Foundation for Statistical Computing, Vienna, Austria) that can easily perform these methods. Results: Through numerical analyses of artificially simulated data and the National Wilms Tumor Study data, accurate risk ratio and risk difference estimates were obtained using the pseudo-Poisson and pseudo-normal linear regression methods. Also, using the auxiliary variable information from the whole cohort, precisions of these estimators were markedly improved. Conclusion: The ordinary logistic regression analyses may provide uninterpretable effect measure estimates, and the risk ratio and risk difference estimation methods are effective alternative approaches for case-cohort studies. These methods are especially recommended under situations in which the event rate is not low.
We obtain the light-front wavefunctions for the nucleon in the valence quark Fock space from an effective Hamiltonian, which includes the transverse and longitudinal confinement and the one gluon exchange interaction with fixed coupling. The wavefunctions are generated by solving the eigenvalue equation in a basis light front quantization. Fitting the model parameters, the wavefunctions lead to good simultaneous description of electromagnetic form factors, radii, and parton distribution functions for the proton.
The Sandia Matrices are a free alternative to the Raven’s Progressive Matrices (RPMs). This study offers a psychometric review of Sandia Matrices items focused on two of the most commonly investigated issues regarding the RPMs: (a) dimensionality and (b) sex differences. Model-data fit of three alternative factor structures are compared using confirmatory multidimensional item response theory (IRT) analyses, and measurement equivalence analyses are conducted to evaluate potential sex bias. Although results are somewhat inconclusive regarding factor structure, results do not show evidence of bias or mean differences by sex. Finally, although the Sandia Matrices software can generate infinite items, editing and validating items may be infeasible for many researchers. Further, to aide implementation of the Sandia Matrices, we provide scoring materials for two brief static tests and a computer adaptive test. Implications and suggestions for future research using the Sandia Matrices are discussed.
Electrochemical impedance spectroscopy (EIS) is a valuable technique to detect the health status and aging phenomena in lithium-ion batteries (LiB). Equivalent circuit modeling (ECM) is conventionally used when interpreting EIS data and gaining physical insights into the aging mechanisms. However, performing ECM is resource intensive and expert-level of knowledge is usually required to select suitable models and fitting parameters. This article presents a quick and user-friendly data analysis algorithm as an alternative to ECM by mathematically fitting geometric features in Nyquist plots and obtaining the growth trends of the features. The evolving trends in the Nyquist plots, such as chord lengths of the arcs and interception points, are consistent with the growth of resistance components obtained using ECM with R 2 values from 0.67 to 0.99, and therefore can be used as indicators of battery aging. Our results show that the quick-fitting approach is suitable for analyzing a series of EIS data acquired during battery cycling and identifying the underlying aging mechanisms.
One of the most common types of models that helps us to understand neuron behavior is based on the Hodgkin–Huxley ion channel formulation (HH model). A major challenge with inferring parameters in HH models is non-uniqueness: many different sets of ion channel parameter values produce similar outputs for the same input stimulus. Such phenomena result in an objective function that exhibits multiple modes (i.e., multiple local minima). This non-uniqueness of local optimality poses challenges for parameter estimation with many algorithmic optimization techniques. HH models additionally have severe non-linearities resulting in further challenges for inferring parameters in an algorithmic fashion. To address these challenges with a tractable method in high-dimensional parameter spaces, we propose using a particular Markov chain Monte Carlo (MCMC) algorithm, which has the advantage of inferring parameters in a Bayesian framework. The Bayesian approach is designed to be suitable for multimodal solutions to inverse problems. We introduce and demonstrate the method using a three-channel HH model. We then focus on the inference of nine parameters in an eight-channel HH model, which we analyze in detail. We explore how the MCMC algorithm can uncover complex relationships between inferred parameters using five injected current levels. The MCMC method provides as a result a nine-dimensional posterior distribution, which we analyze visually with solution maps or landscapes of the possible parameter sets. The visualized solution maps show new complex structures of the multimodal posteriors, and they allow for selection of locally and globally optimal value sets, and they visually expose parameter sensitivities and regions of higher model robustness. We envision these solution maps as enabling experimentalists to improve the design of future experiments, increase scientific productivity and improve on model structure and ideation when the MCMC algorithm is applied to experimental data.
Recent alluvial sediments in riverbeds play a significant role in controlling hydrologic exchange flows (HEFs) in river systems. The alluvial layer is usually associated with strong heterogeneity in physical properties (e.g., permeability and hydraulic conductivity), which affects local HEFs and therefore biogeochemical processes. The spatial distribution of these physical properties needs to be determined to inform the numerical models used to reveal the realistic hydro-biogeochemical behaviors. Such information can be obtained based on the intrinsic link between sediment grain-size distribution and hydraulic properties where sediment texture information is available. However, grain-size measurements are usually spatially sparse and do not have adequate coverage and resolution, particularly for a relatively large domain such as the Hanford Reach of the Columbia River. In this paper, we adopted machine learning (ML) approaches for categorizing and mapping the spatial distributions of riverbed substrate grain size and filling in missing areas of substrate data using the ML models along the reach. Such ML models for substrate size mapping were trained at 13,372 locations using measured substrate sizes along with observed and simulated attributes, including bathymetric attributes (e.g., elevation, slope, and aspect ratio) from LIDAR and bathymetric surveys, and hydrodynamic properties (e.g., water depth, velocity, shear stress, and their statistical moments). An ensemble bagging-based ML technique, Random Forest, was adopted to identify the most influential factors as predictors to develop the predictive models with over-fitting issues addressed. The models were evaluated with respect to each individual substrate size class and the lumped group, and then used to generate the final substrate size maps covering all the grid cells in the numerical modeling domain.