Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical 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 1,045 records · Page 58

The Development of the Pooled Rideshare Acceptance Model (PRAM)

Due to the advancements in real-time information communication technologies and sharing economies, rideshare services have gained significant momentum by offering dynamic and/or on-demand services. Rideshare service companies evolved from personal rideshare, where riders traveled solo or with known individuals, into pooled rideshare (PR), where riders can travel with one to multiple unknown riders. Similar to other shared economy services, pooled rideshare is beneficial as it efficiently utilizes resources, resulting in reduced energy usage, as well as reduced costs for the riders. However, previous research has demonstrated that riders have concerns about using pooled rideshare, especially regarding personal safety. A U.S. national survey with 5385 participants was used to understand human factor-related barriers and user preferences to develop a novel Pooled Rideshare Acceptance Model (PRAM). This model used a covariance-based structural equation model (CB-SEM) to identify the relationships between willingness to consider PR factors (time/cost, privacy, safety, service experience, and traffic/environment) and optimizing one’s experience of PR factors (vehicle technology/accessibility, convenience, comfort/ease of use, and passenger safety), resulting in the higher-order factor trust service. We examined the factors’ relative contribution to one’s willingness/attitude towards PR and user acceptance of PR. Privacy, safety, trust service, and convenience were statistically significant factors in the model, as were the comfort/ease of use factor and the service experience, traffic/environment, and passenger safety factors. The only two non-significant factors in the model were time/cost and vehicle technology/accessibility; it is only when a rider feels safe that individuals then consider the additional non-significant variables of time, cost, technology, and accessibility. Privacy, safety, and service experience were factors that discouraged the use of PR, whereas the convenience factor greatly encouraged the acceptance of PR. Despite the time/cost factor’s lack of significance, individual items related to time and cost were crucial when viewed within the context of convenience. This highlights that while user perceptions of privacy and safety are paramount to their attitude towards PR, once safety concerns are addressed, and services are deemed convenient, time and cost elements significantly enhance their trust in pooled rideshare services. This study provides a comprehensive understanding of user acceptance of PR services and offers actionable insights for policymakers and rideshare companies to improve their services and increase user adoption.

dynamic rideshare system↗

Low-energy enhancement of the magnetic dipole radiation in odd-mass lanthanides

We compute the magnetic dipole (M1) $\gamma$-ray strength functions ($\gamma$SF) for the odd-mass lanthanides $^{\textrm{143-151}}$Nd and $^{\textrm{147-153}}$Sm using the shell-model Monte Carlo method in combination with the static-path approximation and the maximum-entropy method. In particular, we quantify the statistical uncertainties in the calculated M1 $\gamma$SFs and show that they are under control for the excitation energies relevant to the experiments despite a Monte Carlo sign problem that originates in the projection onto an odd number of neutrons. We identify a low-energy enhancement (LEE) in the M1 $\gamma$SFs of these odd-mass lanthanides, which was recently observed experimentally in some of them. We also find a scissors mode resonance (SR) in the strongly deformed isotopes. We observe that the decrease in the LEE strength with neutron number along an isotopic chain is compensated for by an increase in the SR strength in the deformed nuclei. Furthermore, we compare our results with recent experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Demographics, Stellar Populations, and Star Formation Histories of Fast Radio Burst Host Galaxies: Implications for the Progenitors

We present a comprehensive catalog of observations and stellar population properties for 23 highly secure host galaxies of fast radio bursts (FRBs). Our sample comprises 6 repeating FRBs and 17 apparent nonrepeaters. We present 82 new photometric and 8 new spectroscopic observations of these hosts. Using stellar population synthesis modeling and employing nonparametric star formation histories (SFHs), we find that FRB hosts have a median stellar mass of ≈10 9.9 M ⊙ , mass-weighted age ≈5.1 Gyr, and ongoing star formation rate ≈1.3 M ⊙ yr -1 but span wide ranges in all properties. Classifying the hosts by degree of star formation, we find that 87% (20 of 23 hosts) are star-forming, two are transitioning, and one is quiescent. The majority trace the star-forming main sequence of galaxies, but at least three FRBs in our sample originate in less-active environments (two nonrepeaters and one repeater). Across all modeled properties, we find no statistically significant distinction between the hosts of repeaters and nonrepeaters. However, the hosts of repeating FRBs generally extend to lower stellar masses, and the hosts of nonrepeaters arise in more optically luminous galaxies. While four of the galaxies with the clearest and most prolonged rises in their SFHs all host repeating FRBs, demonstrating heightened star formation activity in the last ≲100 Myr, one nonrepeating host shows this SFH as well. Our results support progenitor models with short delay channels (i.e., magnetars formed via core-collapse supernova) for most FRBs, but the presence of some FRBs in less-active environments suggests a fraction form through more delayed channels.

79 ASTRONOMY AND ASTROPHYSICS↗

Accelerating phase field simulations through a hybrid adaptive Fourier neural operator with U-net backbone

Prolonged contact between a corrosive liquid and metal alloys can cause progressive dealloying. For one such process as liquid-metal dealloying (LMD), phase field models have been developed to understand the mechanisms leading to complex morphologies. However, the LMD governing equations in these models often involve coupled non-linear partial differential equations (PDE), which are challenging to solve numerically. In particular, numerical stiffness in the PDEs requires an extremely refined time step size (on the order of 10 -12 s or smaller). This computational bottleneck is especially problematic when running LMD simulation until a late time horizon is required. This motivates the development of surrogate models capable of leaping forward in time, by skipping several consecutive time steps at-once. In this paper, we propose a U-shaped adaptive Fourier neural operator (U-AFNO), a machine learning (ML) based model inspired by recent advances in neural operator learning. U-AFNO employs U-Nets for extracting and reconstructing local features within the physical fields, and passes the latent space through a vision transformer (ViT) implemented in the Fourier space (AFNO). We use U-AFNOs to learn the dynamics of mapping the field at a current time step into a later time step. We also identify global quantities of interest (QoI) describing the corrosion process (e.g., the deformation of the liquid-metal interface, lost metal, etc.) and show that our proposed U-AFNO model is able to accurately predict the field dynamics, in spite of the chaotic nature of LMD. Most notably, our model reproduces the key microstructure statistics and QoIs with a level of accuracy on par with the high-fidelity numerical solver, while achieving a significant 11, 200 × speed-up on a high-resolution grid when comparing the computational expense per time step. Finally, we also investigate the opportunity of using hybrid simulations, in which we alternate forward leaps in time using the U-AFNO with high-fidelity time stepping. We demonstrate that while advantageous for some surrogate model design choices, our proposed U-AFNO model in fully auto-regressive settings consistently outperforms hybrid schemes.

36 MATERIALS SCIENCE↗

Mapping Incidence and Prevalence Peak Data for SIR Modeling Applications

Infectious disease modeling and forecasting have played a key role in helping assess and respond to epidemics and pandemics. Recent work has leveraged data on disease peak infection and peak hospital incidence to fit compartmental models for the purpose of forecasting and describing the dynamics of a disease outbreak. Incorporating these data can greatly stabilize a compartmental model fit on early observations, where slight perturbations in the data may lead to model fits that forecast wildly unrealistic peak infection. We introduce a new method for incorporating historic data on the value and time of peak incidence of hospitalization into the fit for a Susceptible-Infectious-Recovered (SIR) model by formulating the relationship between an SIR model’s starting parameters and peak incidence as a system of two equations that can be solved computationally. We demonstrate how to calculate SIR parameter estimates – which describe disease dynamics such as transmission and recovery rates – using this method, and determine that there is a noticeable loss in accuracy whenever prevalence data is misspecified as incidence data. To exhibit the modeling potential, we update the Dirichlet-Beta State Space modeling framework to use hospital incidence data, as this framework was previously formulated to incorporate only data on total infections. This approach is assessed for practicality in terms of accuracy and speed of computation via simulation.

97 MATHEMATICS AND COMPUTING↗

Materials Characterization, Prediction and Control Project: Summary Report on Data Analytics Framework

This report summarizes the activities performed under the data analytics Vertex in the Materials Characterization, Prediction and Control Project funded under laboratory directed research and development at Pacific Northwest National Laboratory. The data analytics Vertex developed models for associating global or local process parameters, microstructural features, and performance properties of friction-stir-processed 316L stainless steel plates. Statistical, machine learning, and deep learning models, as well as generative artificial intelligence approaches, were used to develop the associations between the process-structure-property data streams. These associations formed the basis for predicting global properties of parts manufactured under different process envelopes, providing a basis for predicting performance using data driven as well as physics-informed and physics-constrained approaches. Additionally, the associations were used to predict local process parameters and microstructural features of the product, predictive relationships that have the potential to form the basis of a control framework that could eventually modulate a friction-stir process to maintain product quality.

316L stainless steel↗

Development of QSAR models to predict blood-brain barrier permeability

Assessing drug permeability across the blood-brain barrier (BBB) is important when evaluating the abuse potential of new pharmaceuticals as well as developing novel therapeutics that target central nervous system disorders. One of the gold-standard in vivo methods for determining BBB permeability is rodent log BB; however, like most in vivo methods, it is time-consuming and expensive. In the present study, two statistical-based quantitative structure-activity relationship (QSAR) models were developed to predict BBB permeability of drugs based on their chemical structure. The in vivo BBB permeability data were harvested for 921 compounds from publicly available literature, non-proprietary drug approval packages, and University of Washington’s Drug Interaction Database. The cross-validation performance statistics for the BBB models ranged from 82 to 85% in sensitivity and 80–83% in negative predictivity. Additionally, the performance of newly developed models was assessed using an external validation set comprised of 83 chemicals. Overall, performance of individual models ranged from 70 to 75% in sensitivity, 70–72% in negative predictivity, and 78–86% in coverage. The predictive performance was further improved to 93% in coverage by combining predictions across the two software programs. These new models can be rapidly deployed to predict blood brain barrier permeability of pharmaceutical candidates and reduce the use of experimental animals.

59 BASIC BIOLOGICAL SCIENCES↗

Map-based cosmology inference with lognormal cosmic shear maps

ABSTRACT Most cosmic shear analyses to date have relied on summary statistics (e.g. ξ+ and ξ−). These types of analyses are necessarily suboptimal, as the use of summary statistics is lossy. In this paper, we forward-model the convergence field of the Universe as a lognormal random field conditioned on the observed shear data. This new map-based inference framework enables us to recover the joint posterior of the cosmological parameters and the convergence field of the Universe. Our analysis properly accounts for the covariance in the mass maps across tomographic bins, which significantly improves the fidelity of the maps relative to single-bin reconstructions. We verify that applying our inference pipeline to Gaussian random fields recovers posteriors that are in excellent agreement with their analytical counterparts. At the resolution of our maps – and to the extent that the convergence field can be described by the lognormal model – our map posteriors allow us to reconstruct all summary statistics (including non-Gaussian statistics). We forecast that a map-based inference analysis of LSST-Y10 data can improve cosmological constraints in the σ8–Ωm plane by $\approx\!{30}{{\ \rm per\ cent}}$ relative to the currently standard cosmic shear analysis. This improvement happens almost entirely along the $S_8=\sigma _8\Omega _{\rm m}^{1/2}$ directions, meaning map-based inference fails to significantly improve constraints on S8.

79 ASTRONOMY AND ASTROPHYSICS↗

Dark and luminous mass components of Omega Centauri from stellar kinematics

Here we combine proper motion data from Gaia EDR 3 and HST with line-of-sight velocity data to study the stellar kinematics of the ω Cen globular cluster. Using a steady-state, axisymmetric dynamical model, we measure the distribution of both the dark and luminous mass components. Assuming both Gaussian and Navarro–Frenk–White mass profiles, depending on the data set, we measure an integrated mass of ≲10 6 M ⊙ within the ω Cen half-light radius for a dark component that is distinct from the luminous stellar component. For the HST and radial velocity data, models with a non-luminous mass component are strongly statistically preferred relative to a stellar mass-only model with a constant mass-to-light ratio. While a compact core of stellar remnants may account for a dynamical mass up to ~5 x 10 5 M ⊙ , they likely cannot explain the higher end of the range. This leaves open the possibility that this non-luminous dynamical mass component comprises non-baryonic dark matter. In comparison to the dark matter distributions around dwarf spheroidal galaxies, the ω Cen dark mass component is much more centrally concentrated. Interpreting the non-luminous mass distribution as particle dark matter, we use these results to obtain the J-factor, which sets the sensitivity to the annihilation cross-section. For the data sets considered, the range of median J-factors is ~10 22 -10 24 GeV 2 cm -5 , which is larger than that obtained for any dwarf spheroidal galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

Prospects for Neutron Reactions on Excited States in High-Density Plasmas

With the reactions of high flux neutrons, such as in a DT plasma, there is a prospect of seeing new kinds of neutron-nucleus reactions for the first time. If neutrons excite a heavy nucleus, for example, there is possibility of a second neutron reacting on excited states of the residual nucleus before that nucleus has de-excited to its ground state. The possibility of such reactions on excited states has rarely been considered. The cross section for neutron induced fission on the isomeric state of 235 U has been measured (D’Eer et al., Phys. Rev. C, 1988, 38: 1270–1276) and calculated (Younes et al., 2003, Maslov, 2007), and reactions on rotationally excited nuclear states has been calculated (Kawano et al., Phys. Rev. C, 2009, 80: 024611). In high flux plasmas, however, a much wider range of reactions is possible. We therefore need to consider excited states at much higher-energies than previously modeled, and then estimate whether second neutrons are likely to rescatter on those excited states. To determine the likelihood of such rescattering events, we first need to know the probable time series of nuclear decays of those excited states. The lifetimes of many low-lying states have been measured experimentally, but now we need to know the lifetimes of the many higher excited states that could be produced from incident 14 MeV neutrons. These are too numerous to be measured and also too numerous to be calculated individually, so statistical Hauser-Feshbach decay models are used. I show some lifetime calculations for 89 Y, 169 Tm, and 197 Au targets, and predictions for the number of rescattering events in plausible plasma scenarios.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement-based quantum simulation of Abelian lattice gauge theories

Numerical simulation of lattice gauge theories is an indispensable tool in high energy physics, and their quantum simulation is expected to become a major application of quantum computers in the future. In this work, for an Abelian lattice gauge theory in d d spacetime dimensions, we define an entangled resource state (generalized cluster state) that reflects the spacetime structure of the gauge theory. We show that sequential single-qubit measurements with the bases adapted according to the former measurement outcomes induce a deterministic Hamiltonian quantum simulation of the gauge theory on the boundary. Our construction includes the (2+1) ( 2 + 1 ) -dimensional Abelian lattice gauge theory simulated on three-dimensional cluster state as an example, and generalizes to the simulation of Wegner’s lattice models M_{(d,n)} M ( d , n ) that involve higher-form Abelian gauge fields. We demonstrate that the generalized cluster state has a symmetry-protected topological order with respect to generalized global symmetries that are related to the symmetries of the simulated gauge theories on the boundary.Our procedure can be generalized to the simulation of Kitaev’s Majorana chain on a fermionic resource state. We also study the imaginary-time quantum simulation with two-qubit measurements and post-selections, and a classical-quantum correspondence, where the statistical partition function of the model M_{(d,n)} M ( d , n ) is written as the overlap between the product of two-qubit measurement bases and the wave function of the generalized cluster state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

OmicsMLMentor: A Web Application for Guided Machine Learning Analysis of Omics Data

Expression-based omics technologies (e.g. proteomics, metabolomics, transcriptomics, etc.) increasingly rely on supervised and unsupervised machine learning (ML) models to find key biomolecules distinguishing conditions, identify natural groupings in biological data, or generate predictions for outcomes of interest. Fitting ML models to omics data presents several challenges, including handling missing data, selecting a normalization method, choosing a valid model, and optimizing hyperparameters, all requiring statistical programming skills to address these challenges. Thus, the open-source web application SLOPE was designed to lower the barrier to ML modeling for omics data. SLOPE supports the fitting of 15 ML models (10 supervised and 5 unsupervised) tailored to omics datasets, such as proteomics, metabolomics, lipidomics, and transcriptomics. SLOPE offers several omics-specific features, including methods for handling missingness (imputation, conversion, removal), normalization tests, ranking of models based on the structure of a user’s data and user input, and optimal hyperparameter selections using cross-validation splits. By streamlining ML workflows for omics analysis, SLOPE address critical gaps in existing online web tools, facilitating a broader adoption of these models for omics research. Here, SLOPE is applied to data from a lignin exposure study to highlight the workflow for fitting both supervised and unsupervised models to data.

lipidomics↗

Mapping the kinetic evolution of metastable grain boundaries under non-equilibrium processing

The kinetic evolution of a multiplicity of metastable grain boundaries (GBs) under fast driving conditions are studied by atomistic modeling. Assisted with an enhanced statistical analysis, the energetic evolution of GBs over a broad metastability-temperature space is mapped out, wherein two distinct regimes—an ageing regime and a rejuvenating regime—are retrieved with high fidelity. By comparing the results under various conditions (e.g. random perturbations, isothermal annealing, and fast heating/cooling), it is shown that such ageing/rejuvenating mechanism map is universal, irrespective of the actual stimuli used to elicit the metastable GBs. Here, the ageing/rejuvenating phenomena are demonstrated to stem from the energy imbalance of uphill climbing and downhill dropping during sequential transitions in the system's potential energy landscape. Without the necessity of introducing free parameters, such model can reconcile experimentally measured hardness variation of nanocrystalline metals subjected to femto-second laser irradiation, and it therefore provides a novel perspective on achieving a plurality of interfacial states and facilitating previously inaccessible property regimes.

36 MATERIALS SCIENCE↗

An Improved Risk Priority Measure for the Management of Project Risk

Monte Carlo-based project risk models are essential tools in project risk management, particularly for large, complex projects in which it is difficult to glean intuitively the importance of individual risks to the overall project. Conventional correlation-based methods for measuring the importance of contributing risks in such models are argued here to be both statistically flawed and unintuitive in the importance metrics they produce. Drawing on and adapting methods more traditionally associated with engineering risk models, Risk Reduction Worth (RRW) is introduced as an importance measure. RRW is shown to be statistically robust and provide an easily interpretable, quantitative measure of the importance of a risk to the overall performance of a project; specifically, measuring the extent to which removal of a risk would improve confidence in project outcomes. To ensure practicality of the method, an approach is described in which RRWs can be produced for all contributing risks based on a single Monte Carlo sample. This is demonstrated through application to a large, complex project risk model from which intuitive risk importance insights provide a robust, transparent basis for the allocation of risk management resources.

Risk↗

Microscopic constraints for the equation of state and structure of neutron stars: A Bayesian model mixing framework

Bayesian model mixing (BMM) is a statistical technique that can combine constraints from different regions of an input space in a principled way. Here we extend our BMM framework for the equation of state (EOS) of strongly interacting matter from symmetric nuclear matter to asymmetric matter, specifically focusing on zero-temperature, charge-neutral, 𝛽-equilibrated matter. We use Gaussian processes (GPs) to infer constraints on the neutron-star matter EOS at intermediate densities from two different microscopic theories: chiral effective-field theory (𝜒⁢EFT) at baryon densities around nuclear saturation, 𝑛 𝐵 ∼ 𝑛 0 , and perturbative QCD at asymptotically high baryon densities, 𝑛 𝐵 ⩾ 20⁢𝑛 0 . The uncertainties of the 𝜒⁢EFT and pQCD EOSs are obtained using the BUQEYE truncation error model. We demonstrate the flexibility of our framework through the use of two categories of GP kernels: conventional stationary kernels and a nonstationary changepoint kernel. We use the latter to explore potential constraints on the dense matter EOS by including exogenous data representing theory predictions and heavy-ion collision measurements at densities ⩾ 2⁢𝑛 0 . We also use our EOSs to obtain neutron-star mass-radius relations and their uncertainties. Finally, our framework, whose implementation will be available through a GitHub repository, provides a prior distribution for the EOS that can be used in large-scale neutron-star inference frameworks.

Bayesian methods↗

Evaluating Tropical Precipitation Relations in CMIP6 Models with ARM Data

Abstract A set of diagnostics based on simple, statistical relationships between precipitation and the thermodynamic environment in observations is implemented to assess phase 6 of the Coupled Model Intercomparison Project (CMIP6) model behavior with respect to precipitation. Observational data from the Atmospheric Radiation Measurement (ARM) permanent field observational sites are augmented with satellite observations of precipitation and temperature as an observational baseline. A robust relationship across observational datasets between column water vapor (CWV) and precipitation, in which conditionally averaged precipitation exhibits a sharp pickup at some critical CWV value, provides a useful convective onset diagnostic for climate model comparison. While a few models reproduce an appropriate precipitation pickup, most models begin their pickup at too low CWV and the increase in precipitation with increasing CWV is too weak. Convective transition statistics compiled in column relative humidity (CRH) partially compensate for model temperature biases—although imperfectly since the temperature dependence is more complex than that of column saturation. Significant errors remain in individual models and weak pickups are generally not improved. The conditional-average precipitation as a function of CRH can be decomposed into the product of the probability of raining and mean precipitation during raining times (conditional intensity). The pickup behavior is primarily dependent on the probability of raining near the transition and on the conditional intensity at higher CRH. Most models roughly capture the CRH dependence of these two factors. However, compensating biases often occur: model conditional intensity that is too low at a given CRH is compensated in part by excessive probability of precipitation.

54 ENVIRONMENTAL SCIENCES↗

Indirect Tool Condition Monitoring Using Ensemble Machine Learning Techniques

Abstract Tool condition monitoring (TCM) has become a research area of interest due to its potential to significantly reduce manufacturing costs while increasing process visibility and efficiency. Machine learning (ML) is one analysis technique which has demonstrated advantages for TCM applications. However, the commonly studied individual ML models lack generalizability to new machining and environmental conditions, as well as robustness to the unbalanced datasets which are common in TCM. Ensemble ML models have demonstrated superior performance in other fields, but have only begun to be evaluated for TCM. As a result, it is not well understood how their TCM performance compares to that of individual models, or how homogeneous and heterogeneous ensemble models’ performances compare to one another. To fill in these research gaps, milling experiments were conducted using various cutting conditions, and the model groups were compared across several performance metrics. Statistical t-tests were also used to evaluate the significance of model performance differences. Through the analysis of four individual ML models and five ensemble models, all based on the processes’ sound, spindle power, and axial load signals, it was found that on average, the ensemble models performed better than the individual models, and that the homogeneous ensembles outperformed the heterogeneous ensembles.

Engineering↗