Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hierarchy”

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

Multi-fidelity microstructure-induced uncertainty quantification by advanced Monte Carlo methods

Quantifying uncertainty associated with the microstructure variation of a material can be a computationally daunting task, especially when dealing with advanced constitutive models and fine mesh resolutions in the crystal plasticity finite element method (CPFEM). Numerous studies have been conducted regarding the sensitivity of material properties and performance to the mesh resolution and choice of constitutive model. However, a unified approach that accounts for various fidelity parameters, such as mesh resolutions, integration time-steps and constitutive models simultaneously is currently lacking. Further, this paper proposes a novel uncertainty quantification (UQ) approach for computing the properties and performance of homogenized materials using CPFEM, that exploits a hierarchy of approximations with different levels of fidelity. In particular, we illustrate how multi-level sampling methods, such as multi-level Monte Carlo (MLMC) and multi-index Monte Carlo (MIMC), can be applied to assess the impact of variations in the microstructure of polycrystalline materials on the predictions of homogenized materials properties. We show that by adaptively exploiting the fidelity hierarchy, we can significantly reduce the number of microstructures required to reach a certain prescribed accuracy. Finally, we show how our approach can be extended to a multi-fidelity framework, where we allow the underlying constitutive model to be chosen from either a phenomenological plasticity model or a dislocation-density-based model.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Towards performance portability in the Spark astrophysical magnetohydrodynamics solver in the Flash-X simulation framework

Simulations of core-collapse supernovae, and other astrophysical phenomena, are quintessential extreme-scale computing challenges. For core-collapse supernova simulations to be carried out by the ExaStar project under the Exascale Computing Project umbrella, a robust, efficient, and state-of-the-art magnetohydrodynamics solver is a critical requirement. In Flash-X, the primary software instrument for ExaStar, a new magnetohydrodynamics solver has been designed and implemented from the ground up to achieve accuracy and efficiency for simulations of complex astrophysical flows. This new solver, dubbed Spark, uses high-order spatial reconstruction, Runge-Kutta time integration, and an efficient cell-centered approach to satisfying the divergence-free condition for the magnetic fields. Spark was written to be optimized for data locality in cache hierarchy of CPUs. Since data locality optimizations for cache hierarchy are not directly compatible with those of accelerators, we have taken the approach of using program synthesis to avoid massive amounts of code replication that would be necessary if we were to maintain two different versions of the solver. Our program synthesis relies on a simple key-dictionary approach, implemented in python, that enables us to assemble the version of the solver suitable for the target hardware from code fragments identified by specific keys. In this work, we describe the data locality optimizations of the solver for CPUs and accelerators and the program synthesis tools that enable this portability. We also detail the parallel performance of Spark for both CPUs and accelerators.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Review of neutral naturalness

The hierarchy between the mass parameter of the Higgs boson and larger mass scales becomes ever more puzzling as experiments explore higher energies. Neutral naturalness is the umbrella term for symmetry-based explanations for these hierarchies whose quark symmetry partners are not charged under the SU(3) c color gauge group of the Standard Model. Though the first manifestations of this idea predate the physics runs of the Large Hadron Collider, since the Higgs discovery this paradigm has grown and developed to include a wide variety of concrete realizations with connections to intriguing collider signals. Determining the phenomenology of such models often requires the characterization—typically relying on lattice calculations—of a new confining gauge symmetry. This presents additional motivation to further develop our understanding of nonperturbative field theory as well as to pursue specific lattice studies. The wide range of suggested hidden sectors also produces a variety of dark matter candidates, intersections with astrophysics and cosmology, and ties to neutrinos and flavor. In this review, we orient the reader within both this growing collection of specific models and the physical phenomena they produce. We also survey the often less familiar dynamics of hidden-sector glueballs and quirks. In addition to providing a guide to past efforts, we reveal interesting directions for further study.

Batell, Brian [Univ. of Pittsburgh, PA (United Sta↗
DOE OSTI2023

Non-unimodal and non-concave relationships in the network Macroscopic Fundamental Diagram caused by hierarchical streets

Unimodal, concave relationships between average network productivity and accumulation or density aggregated across spatially compact regions of urban networks—so called network Macroscopic Fundamental Diagrams (MFDs)—have recently been shown to exist on homogeneous street networks. When present, MFD relationships facilitate the modeling of traffic congestion at a regional level and have led to the development of various regional traffic control strategies. However, real street networks are not homogeneous—they generally have a hierarchical structure where some streets (e.g., arterials) promote higher mobility than others (e.g., local roads). Here, this paper examines how the presence of hierarchical roadway structures may potentially cause non-unimodal patterns in a network's MFD. These are observed using three types of tools: analytical models of simple network structures, simulations of various idealized roadway networks, and empirical data. The impacts of street hierarchy depend on how vehicles use different roadway types to move within the network; i.e., their routing strategy. The findings suggest that the presence of roadway hierarchies may lead to MFDs that have non-unimodal or non-concave patterns on the free-flow branch when vehicles route themselves according to user equilibrium principles, which is closest to what would be observed in realistic situations. Such patterns are contrary to what is traditionally assumed in most MFD-based modeling frameworks. However, the unimodal and concave MFD should be expected under system optimal routing conditions that maximize network productivity for a given traffic state.

42 ENGINEERING↗
DOE OSTI2021

Scalable ab initio fragmentation methods based on a truncated expansion of the non-orthogonal molecular orbital model

An alternative formulation of the non-orthogonal molecular orbital model of electronic structure theory is developed based on the expansion of the inverse molecular orbital overlap matrix. From this model, a hierarchy of ab initio fragment-based quantum chemistry methods, referred to as the nth-order expanded non-orthogonal molecular orbital methods, are developed using a minimal number of approximations, each of which is frequently employed in intermolecular interaction theory. These novel methods are compared to existing fragment-based quantum chemistry methods, and the implications of those significant differences, where they exist, between the methods developed herein and those already existing methods are examined in detail. Computational complexities and theoretical scaling are also analyzed and discussed. Future extensions for the hierarchy of methods, to account for additional intrafragment and interfragment interactions, are outlined.

Westheimer, Bryce M. (ORCID:0000000267645714)↗
DOE OSTI2022

Dynamics of photosynthetic light harvesting systems interacting with N-photon Fock states

Here, we develop a method to simulate the excitonic dynamics of realistic photosynthetic light harvesting systems, including non-Markovian coupling to phonon degrees of freedom, under excitation by N-photon Fock state pulses. This method combines the input–output and the hierarchical equations of motion formalisms into a double hierarchy of density matrix equations. We show analytically that under weak field excitation relevant to natural photosynthesis conditions, an N-photon Fock state input and a corresponding coherent state input give rise to equal density matrices in the excited manifold. However, an N-photon Fock state input induces no off-diagonal coherence between the ground and excited subspaces, in contrast with the coherences created by a coherent state input. We derive expressions for the probability to absorb a single Fock state photon with or without the influence of phonons. For short pulses (or, equivalently, wide bandwidth pulses), we show that the absorption probability has a universal behavior that depends only upon a system-dependent effective energy spread parameter Δ and an exciton–light coupling constant Γ. This holds for a broad range of chromophore systems and for a variety of pulse shapes. We also analyze the absorption probability in the opposite long pulse (narrow bandwidth) regime. We then derive an expression for the long time emission rate in the presence of phonons and use it to study the difference between collective vs independent emission. Finally, we present a numerical simulation for the LHCII monomer (14-mer) system under single photon excitation that illustrates the use of the double hierarchy equations.

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

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Chemistry↗
DOE OSTI2026

A Fast and Accurate Reduced-Order Model for High-Intensity Transferred Arc Discharges

Arc discharges are widely used in welding, plasma smelting, and other industrial processes, where variations in operating conditions strongly affect arc stability, temperature distribution, and energy transfer. Accurate modeling of these phenomena typically requires computationally expensive high-fidelity simulations. This study presents a hierarchy of three arc discharge models with progressively reduced physical fidelity. The high-fidelity model provides a fully physics-resolved reference, the reduced-order model reproduces the dominant thermal and flow characteristics of the arc with reasonable accuracy, and the Elenbaas-Heller model captures key trends in a simplified, rapid formulation. This hierarchy demonstrates that reduced-order model can effectively balance predictive fidelity and computational efficiency, providing practical tools for arc simulation and parametric studies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Accelerating multilevel Markov Chain Monte Carlo using machine learning models

Here, this work presents an efficient approach for accelerating multilevel Markov Chain Monte Carlo (MCMC) sampling for large-scale problems using low-fidelity machine learning models. While conventional techniques for large-scale Bayesian inference often substitute computationally expensive high-fidelity models with machine learning models, thereby introducing approximation errors, our approach offers a computationally efficient alternative by augmenting high-fidelity models with low-fidelity ones within a hierarchical framework. The multilevel approach utilizes the low-fidelity machine learning model (MLM) for inexpensive evaluation of proposed samples thereby improving the acceptance of samples by the high-fidelity model. The hierarchy in our multilevel algorithm is derived from geometric multigrid hierarchy. We utilize an MLM to accelerate the coarse level sampling. Training machine learning model for the coarsest level significantly reduces the computational cost associated with generating training data and training the model. We present an MCMC algorithm to accelerate the coarsest level sampling using MLM and account for the approximation error introduced. We provide theoretical proofs of detailed balance and demonstrate that our multilevel approach constitutes a consistent MCMC algorithm. Additionally, we derive the expression for cost reduction due to machine learning model to facilitate cost analysis of the hierarchical sampling algorithm. Our technique is demonstrated on a standard benchmark inference problem in groundwater flow, where we estimate the probability density of a quantity of interest using a four-level MCMC algorithm. Our proposed algorithm accelerates multilevel sampling by a factor of two while achieving similar accuracy compared to sampling using the standard multilevel algorithm.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

The DUNE Far Detector Vertical Drift Technology. Technical Design Report

DUNE is an international experiment dedicated to addressing some of the questions at the forefront of particle physics and astrophysics, including the mystifying preponderance of matter over antimatter in the early universe. The dual-site experiment will employ an intense neutrino beam focused on a near and a far detector as it aims to determine the neutrino mass hierarchy and to make high-precision measurements of the PMNS matrix parameters, including the CP-violating phase. It will also stand ready to observe supernova neutrino bursts, and seeks to observe nucleon decay as a signature of a grand unified theory underlying the standard model. The DUNE far detector implements liquid argon time-projection chamber (LArTPC) technology, and combines the many tens-of-kiloton fiducial mass necessary for rare event searches with the sub-centimeter spatial resolution required to image those events with high precision. The addition of a photon detection system enhances physics capabilities for all DUNE physics drivers and opens prospects for further physics explorations. Given its size, the far detector will be implemented as a set of modules, with LArTPC designs that differ from one another as newer technologies arise. In the vertical drift LArTPC design, a horizontal cathode bisects the detector, creating two stacked drift volumes in which ionization charges drift towards anodes at either the top or bottom. The anodes are composed of perforated PCB layers with conductive strips, enabling reconstruction in 3D. Light-trap-style photon detection modules are placed both on the cryostat's side walls and on the central cathode where they are optically powered. This Technical Design Report describes in detail the technical implementations of each subsystem of this LArTPC that, together with the other far detector modules and the near detector, will enable DUNE to achieve its physics goals.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

EC-Bench: A Benchmark for Enzyme Commission Number Prediction

Enzymes are proteins that catalyze specific biochemical reactions in cells. Enzyme Commission (EC) numbers are used to annotate enzymes in a four-level hierarchy that classifies enzymes based on the specific chemical reactions they catalyze. Accurate EC number prediction is essential for understanding enzyme functions. Despite the availability of numerous methods for predicting EC numbers from protein sequences, there is no unified framework for evaluating and studying such methods systematically. This gap limits the ability of the community to identify the most effective approaches for enzyme annotation. We introduce EC-Bench, a benchmark for EC number prediction, consisting of 1) an initial representative set of existing methods (including homology-based, deep learning, contrastive learning, and language model methods), 2) existing and novel accuracy and efficiency performance metrics, and 3) selected datasets to allow for comprehensive comparative study. EC-Bench is open-source and provides a framework for researchers to not only compare among existing methods objectively under uniform conditions, but also to introduce and effectively evaluate performance of new methods in a comparative framework. To demonstrate the utility of EC-Bench, we perform extensive experimentation to compare the existing EC number prediction methods and establish their advantages and disadvantages in a variety of prediction tasks, namely “exact EC number prediction”, “EC number completion” and (partial or additional) “EC number recommendation”. We find wide variation in the performance of different methods, but also subtle but potentially useful differences in the performance of different methods across tasks and for different parts of the EC hierarchy.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Patient-Reported Outcomes in Pediatric Cancer Registration Trials: A US Food and Drug Administration Perspective

Pediatric patient-reported outcome (PRO) data can help inform the US Food and Drug Administration’s (FDA’s) benefit-risk assessment of cancer therapeutics by quantifying symptom and functional outcomes from the patient’s perspective. This study assessed use of PROs in commercial pediatric oncology trials submitted to the FDA for regulatory review. FDA databases were searched to identify pediatric oncology product applications approved between 1997 and 2020. Sponsor-submitted documents were reviewed to determine whether PRO data were collected, which instruments were used, and the quality of collected data (ie, sample size, completion rates, and use of fit-for-purpose instruments). The role of PROs in each trial (endpoint hierarchy) was also recorded in addition to whether any PRO endpoints were included in product labeling. We reviewed 17 pediatric oncology applications, 4 of which included PRO data: denosumab, tisagenlecleucel, larotrectinib, and selumetinib. In these 4 instances, PROs served as exploratory endpoints and were not incorporated in product labeling. Trials that collected PRO data were phase II or phase I/II single-arm studies with sample sizes of 28 to 88 patients. Symptomatic adverse events (AEs) were characterized using clinician-reported Common Terminology Criteria for Adverse Events (CTCAE) without additional patient self-report. PROs were infrequently used in pediatric cancer registration trials. When PROs were used, PRO data were limited by lack of a clear research objective and corresponding prospective statistical analysis plan. Contemporary PRO symptom libraries, such as the National Cancer Institute’s Pediatric PRO-CTCAE, may provide an opportunity to better evaluate the occurrence and impact of symptomatic AEs, from the patient’s perspective, in pediatric oncology trials.

Oncology↗
DOE OSTI2026

Emergence of vorticity and viscous stress in finite-scale quantum hydrodynamics

The Madelung equations offer a hydrodynamic description of quantum systems, from single particles to quantum fluids. In this formulation, the probability density is mapped onto the fluid density and the phase is treated as a scalar potential generating the velocity field. As examples of potential flows, quantum fluids described in this way are inherently irrotational, but quantum vortices may arise at discrete points where the phase is undefined. In this paper, starting from this irrotational description of a quantum fluid, a coarse-graining procedure is applied to arrive at a macroscopic description of the quantum fluid in terms of a hierarchy of moments in which the role of velocity is played by a Favre average of the microscopic velocity field. This hierarchy is truncated using an explicit closure derived from an expansion in a finite length scale. The resulting coarse-grained fields are shown to allow for finite vorticity at any point in the fluid. Additionally, it is shown that this vorticity obeys a similar equation to the vorticity equation in classical hydrodynamics and includes a vortex-stretching term. The particular closure employed here also gives rise to a novel stress term in the fluid equations, which in the appropriate limit appears analogous to an artificial viscous stress from computational fluid dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Entanglement structures in quantum field theories: Negativity cores and bound entanglement in the vacuum

Here, the many-body entanglement between two finite (size-d) disjoint vacuum regions of noninteracting lattice scalar field theory in one spatial dimension, i.e., a (d A × d B ) mixed Gaussian continuous variable system, is locally transformed into a tensor-product core of (1 A × 1 B ) mixed entangled pairs. Accessible entanglement within these core pairs exhibits an exponential hierarchy and as such identifies the structure of dominant region modes from which vacuum entanglement could be extracted into a spatially separated pair of quantum detectors. Beyond the core, the remaining modes of the halo are determined to be AB separable in isolation, as well as separable from the core. However, state preparation protocols that distribute entanglement in the form of (1 A × 1 B ) mixed core pairs are found to require additional entanglement in the halo that is obscured by classical correlations. This inaccessible (bound) halo entanglement is found to mirror the accessible entanglement, but with a step behavior as the continuum is approached. It remains possible that alternate initialization protocols that do not utilize the exponential hierarchy of core-pair entanglement may require less inaccessible entanglement. Entanglement consolidation is expected to persist in higher dimensions and may aid classical and quantum simulations of asymptotically free gauge field theories, such as quantum chromodynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Detectability of collective neutrino oscillation signatures in the supernova explosion of a 8.8 M ⊙ star

In order to investigate the impact of collective neutrino oscillations (CNOs) on the neutrino signal from a nearby supernova, we perform three-flavor neutrino oscillation simulations employing the multiangle effect. The background hydrodynamic model is based on the neutrino hydrodynamic simulation of a 8.8 M ⊙ progenitor star. In this work, we find that CNO commences after some 100 ms post bounce. Before this, CNO is suppressed by matter-induced decoherence. In the inverted mass hierarchy, the spectrum of $\bar{ν}_e$ becomes softer after the onset of CNO. To evaluate the detectability of this modification, we define a hardness ratio between the number of high energy neutrino events and low energy neutrino events adopting a fixed critical energy. We show that Hyper-Kamiokande (HK) can distinguish the effect of CNO for supernova distances out to ~ 10 kpc. On the other hand, for the normal mass hierarchy, the spectrum of ν e becomes softer after the onset of CNO, and we show that DUNE can distinguish this feature for supernova distances out to ~ 10 kpc . More work is necessary to optimize the best value of critical energy for maximum sensitivity. We also show that if the spectrum of $\bar{ν}_e$ in HK becomes softer due to CNO, the spectrum of $ν_e$ in DUNE becomes harder, and vice versa. These synergistic observations in $\bar{ν}_e$ and $ν_e$, by HK and DUNE, respectively, will be an intriguing opportunity to test the occurrence of CNO.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Time dependent signatures of core-collapse supernova neutrinos at HALO

In this work, we calculate the response of a lead-based detector, such as the Helium and Lead Observatory (HALO) or its planned upgrade HALO-1kt to a Galactic core-collapse supernova. We pay particular attention to the time dependence of the reaction rates. All reaction rates decrease as the neutrino luminosity exponentially drops during the cooling period, but the ratio of one-neutron (1n) to two-neutron (2n) event rates in HALO is independent of this overall decrease. Nevertheless, we find that this ratio still changes with time due to the changing character of neutrino flavor transformations with the evolving conditions in the supernova. In the case of inverted hierarchy, this is caused by the fact that the spectral splits become less and less sharp with the decreasing luminosity. In the case of normal hierarchy, it is caused by the passage of the shock wave through the Mikheyev-Smirnov-Wolfenstein resonance region. However, in both cases, we find that the change in the ratio of 1n to 2n event rates is limited to a few percent.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Impact of recent updates to neutrino oscillation parameters on the effective Majorana neutrino mass in 0 ν β β decay

We investigate how recent updates to neutrino oscillation parameters and the sum of neutrino masses influence the sensitivity of neutrinoless double-beta ( 0 ν β β ) decay experiments. Incorporating the latest cosmological constraints on the sum of neutrino masses and laboratory measurements on oscillations, we determine the sum of neutrino masses for both the normal hierarchy (NH) and the inverted hierarchy (IH). Our analysis reveals a narrow range for the sum of neutrino masses, approximately 0.06 eV / c 2 for NH and 0.102 eV / c 2 for IH. Utilizing these constraints, we calculate the effective Majorana masses for both NH and IH scenarios, establishing the corresponding allowed regions. Importantly, we find that the minimum neutrino mass is nonzero, as constrained by the current oscillation parameters. Additionally, we estimate the half-life of 0 ν β β decay using these effective Majorana masses for both NH and IH. Our results suggest that upcoming ton-scale experiments will comprehensively explore the IH scenario, while 100-ton-scale experiments will effectively probe the parameter space for the NH scenario, provided the background index can achieve 1 event/kton-year in the region of interest. Published by the American Physical Society 2024

Astronomy & Astrophysics↗