Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “space computing”

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 253 records · Page 14

Evidence for WW/WZ vector boson scattering in the decay channel ℓνqq produced in association with two jets in proton-proton collisions at s = 13 TeV

Evidence is reported for electroweak (EW) vector boson scattering in the decay channel $\ell\nu$qq of two weak vector bosons WV (V = W or Z), produced in association with two parton jets. The search uses a data set of proton-proton collisions at 13 TeV collected with the CMS detector during 2016-2018 with an integrated luminosity of 138 fb$^{-1}$. Events are selected requiring one lepton (electron or muon), moderate missing transverse momentum, two jets with a large pseudorapidity separation and a large dijet invariant mass, and a signature consistent with the hadronic decay of a W/Z boson. The cross section is computed in a fiducial phase space defined at parton level requiring all parton transverse momenta $p_\mathrm{T}$$\gt$ 10 GeV and at least one pair of outgoing partons with invariant mass $m_\mathrm{qq}$$\gt$ 100 GeV. The measured and expected EW WV production cross sections are 1.90 $^{+0.53}_{-0.46}$ pb and 2.23 $^{+0.08}_{-0.11}$ (scale) $\pm$ 0.05 (PDF) pb, respectively, where PDF is the parton distribution function. The observed EW signal strength is $m_\mathrm{EW}$ = 0.85 $\pm$ 0.12 (stat) $^{+0.19}_{-0.17}$ (syst), corresponding to a signal significance of 4.4 standard deviations with 5.1 expected, and it is measured keeping the quantum chromodynamics (QCD) associated diboson production fixed to the standard model prediction. This is the first evidence of vector boson scattering in the $\ell\nu$qq decay channel at LHC. The simultaneous measurement of the EW and QCD associated diboson production agrees with the standard model prediction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A systems engineering vision for floating offshore wind cost optimization

The U.S. offshore wind resource potential is immense, and often in close proximity to densely-populated coastal load centers. In many U.S. coastal areas, water depths favor the deployment of floating over fixed-bottom offshore wind technology. Floating offshore wind plants have the potential to be cost-competitive with fixed-bottom installations, but because the technology has not yet been deployed at commercial scale, it is not clear when and with what configurations this potential cost parity can be achieved. This article first reviews the state of floating offshore wind technology and deployments to identify key gaps that must be addressed to bring down the overall cost of energy produced. The article then puts forth a long-term vision for a research program and design methodology that may be able to push floating wind plants toward a lower levelized cost of energy than fixed-bottom offshore wind. The method involves a fully integrated systems-engineering and techno-economic design approach to capture the complex interactions between the physics, manufacturing, installation, and operation of floating wind turbines to achieve transformational cost reductions. The approach also envisions multifidelity and uncertainty management strategies to examine the most robust and viable concepts in the design trade-space. To better focus the computational resources, engineering lessons learned from existing offshore wind systems and concept studies are used to develop a set of criteria that can be applied to prefilter candidate technology building blocks that have the greatest cost reduction potential.

17 WIND ENERGY↗

The “Fresnel Equations” for Diffuse radiation on Inclined photovoltaic Surfaces (FEDIS)

Here the well-known Fresnel equations solve for the reflection and transmission of light for precise incident angles. The transmission of diffuse radiation incident on a planar or domed surface is often needed for real-world applications. Due to the complexity of the Fresnel equations, the analytical solution of the integration has hitherto been unobtainable over the last centuries. Therefore, this problem was numerically solved by integrating the angular transmittances in space often leading to substantial computing burden and bias in the results. To efficiently estimate the solar energy resource for a glass-covered photovoltaic (PV) module, we derive an analytical solution of diffuse transmission based on the rigorous integration of an alternate form of the Fresnel equations. The approach leads to a simple yet accurate relative transmittance model that reconciles the solar energy sensed by pyranometers and PV panels. With limited and clearly stated approximations, the complex mathematical derivation resulted in an elegant solution. An experiment using 1-year of data at the National Renewable Energy Laboratory's (NREL's) Solar Radiation Research Laboratory (SRRL) shows that the new model dramatically decreases the disparity between the solar radiation measurements by a Kipp and Zonen CM Pyranometer 22 (CMP22) and an IMT reference cell on a 1-axis tracking system. The solution in this paper can be widely used in scientific and engineering research, development, and applications wherever the Fresnel equations are used.

14 SOLAR ENERGY↗

Orthorhombic Structures as Inorganic Halide Perovskite Models for High-Throughput DFT Investigations

Although high-symmetry Pm-3m space group cubic models are computationally efficient for high-throughput density functional theory (DFT) calculations of inorganic ternary (ABX 3 ) halide perovskites (HPs), they frequently predict band gaps (E g ) that disagree with experiment. Conversely, while low-symmetry cubic polymorphous networks (PN) comprised of 160 to 320 atoms incur significantly greater computational cost, they predict E g ’s that are more closely aligned with experiment. Here, in this study, we compare the DFT total energies and E g ’s predicted by four high-symmetry structure models (Pnma orthorhombic, R3m trigonal, P4/mbm tetragonal, and Pm-3m cubic) to cubic PNs for 5 experimentally characterized ternary HPs and find that the orthorhombic model computes Eg’s with the smallest MAD of 0.23 eV relative to the PNs. Pair distribution functions and DFT-computed total energies show that octahedral tilting, which is present in the 20-atom orthorhombic and 160-atom cubic PN models but not in the 5-atom cubic models, stabilizes all 5 compositions in our benchmarking set. We also find that imposing PN constraints when generating and optimizing these orthorhombic structures by fixing the unit cell lattice vectors and displacing the atoms prior to ionic relaxation with DFT, which we call the orthorhombic surrogate model (OSM), lowers the MAD of E g predictions to 0.09 eV. Our OSM predicts the PN band gaps of an additional 95 theoretical inorganic ternary HPs with MAD of 0.08 eV, supporting its usage in high-throughput DFT investigations to closely estimate PN band gaps with much less computational expense.

14 SOLAR ENERGY↗

Efficient generation of grids and traversal graphs in compositional spaces towards exploration and path planning

Abstract Diverse disciplines across science and engineering deal with problems related to compositions, which exist in non-Euclidean simplex spaces, rendering many standard tools inaccurate or inefficient. This work explores such spaces conceptually in the context of materials discovery, quantifies their computational feasibility, and implements several essential methods specific to simplex spaces through a new high-performance open-source library . Most significantly, we derive and implement an algorithm for constructing a novel n-dimensional simplex graph data structure, containing all discretized compositions and possible neighbor-to-neighbor transitions. Critically, no distance or neighborhood calculations are performed, instead leveraging pure combinatorics and order in procedurally generated simplex grids, keeping the algorithm $${\mathcal{O}}(N)$$ O ( N ) , with minimal memory, enabling rapid construction of graphs with billions of transitions in seconds. Additionally, we demonstrate how such graph representations can be combined to homogeneously express complex path-planning problems, while facilitating efficient deployment of existing high-performance gradient descent, graph traversal, and other optimization algorithms.

Krajewski, Adam M. (ORCID:0000000222660099)↗

Iterative computational design and crystallographic screening identifies potent inhibitors targeting the Nsp3 macrodomain of SARS-CoV-2

The nonstructural protein 3 (NSP3) of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) contains a conserved macrodomain enzyme (Mac1) that is critical for pathogenesis and lethality. While small-molecule inhibitors of Mac1 have great therapeutic potential, at the outset of the COVID-19 pandemic, there were no well-validated inhibitors for this protein nor, indeed, the macrodomain enzyme family, making this target a pharmacological orphan. Here, we report the structure-based discovery and development of several different chemical scaffolds exhibiting low- to sub-micromolar affinity for Mac1 through iterations of computer-aided design, structural characterization by ultra-high-resolution protein crystallography, and binding evaluation. Potent scaffolds were designed with in silico fragment linkage and by ultra-large library docking of over 450 million molecules. Both techniques leverage the computational exploration of tangible chemical space and are applicable to other pharmacological orphans. Overall, 160 ligands in 119 different scaffolds were discovered, and 153 Mac1-ligand complex crystal structures were determined, typically to 1 Å resolution or better. Our analyses discovered selective and cell-permeable molecules, unexpected ligand-mediated conformational changes within the active site, and key inhibitor motifs that will template future drug development against Mac1.

60 APPLIED LIFE SCIENCES↗

Feedforward equilibrium trajectory optimization with GSPulse

One of the common tasks required for designing new plasma scenarios or evaluating capabilities of a tokamak is to design the desired equilibria using a Grad-Shafranov (GS) equilibrium solver. However, most standard equilibrium solvers are time-independent and do not include dynamic effects such as plasma current flux consumption, induced vessel currents, or voltage constraints. Another class of tools, plasma equilibrium evolution simulators, do include time-dependent effects. These are generally structured to solve the forward problem of evolving the plasma equilibrium given feedback-controlled voltages. In this work, we introduce GSPulse, a novel algorithm for equilibrium trajectory optimization, that is more akin to a pulse planner than a pulse simulator. GSPulse includes time-dependent effects and solves the inverse problem: given a user-specified set of target equilibrium shapes, as well as limits on the coil currents and voltages, the optimizer returns trajectories of the voltages, currents, and achievable equilibria. This task is useful for scoping performance of a tokamak and exploring the space of achievable pulses. The computed equilibria satisfy both Grad-Shafranov force balance and axisymmetric circuit dynamics. The optimization is performed by restructuring the free-boundary equilibrium evolution equations into a form where it is computationally efficient to optimize the entire dynamic sequence. GSPulse can solve for hundreds of equilibria simultaneously within a few minutes. GSPulse has been validated against NSTX-U and MAST-U experiments and against SPARC feedback control simulations, and is being used to perform scenario design for SPARC. The computed trajectories can be used as feedforward inputs that are connected to the feedback controller to inform and improve feedback performance. The code for GSPulse is available open-source at github.com/jwai-cfs/GSPulse_public.

equilibrium↗

A review on machine learning-guided design of energy materials

Abstract The development and design of energy materials are essential for improving the efficiency, sustainability, and durability of energy systems to address climate change issues. However, optimizing and developing energy materials can be challenging due to large and complex search spaces. With the advancements in computational power and algorithms over the past decade, machine learning (ML) techniques are being widely applied in various industrial and research areas for different purposes. The energy material community has increasingly leveraged ML to accelerate property predictions and design processes. This article aims to provide a comprehensive review of research in different energy material fields that employ ML techniques. It begins with foundational concepts and a broad overview of ML applications in energy material research, followed by examples of successful ML applications in energy material design. We also discuss the current challenges of ML in energy material design and our perspectives. Our viewpoint is that ML will be an integral component of energy materials research, but data scarcity, lack of tailored ML algorithms, and challenges in experimentally realizing ML-predicted candidates are major barriers that still need to be overcome.

36 MATERIALS SCIENCE↗

Newton versus the machine: solving the chaotic three-body problem using deep neural networks

ABSTRACT Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system’s chaotic nature. We show that an ensemble of converged solutions for the planar chaotic three-body problem obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at a fixed computational cost and up to 100 million times faster than the numerical integrator. In addition, we demonstrate the importance of training an ANN using converged solutions from an arbitrary precise integrator, relative to solutions computed by a conventional fixed precision integrator, which can introduce errors in the training data, due to numerical round-off and time discretization, that are learned by the ANN. Our results provide evidence that, for computationally challenging regions of phase space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black hole binary systems or the origin of the core collapse in dense star clusters.

Breen, Philip G.↗

Primordial black hole dark matter: A quantitative parameter sensitivity comparison across formation mechanisms and particle candidates

Primordial black holes (PBHs) in the asteroid-mass window ( 10 17 – 10 22 g ) can account for all of the dark matter without violating any observational constraint, yet are routinely dismissed as fine-tuned. I put that dismissal to the test by applying three complementary sensitivity measures uniformly across a broad landscape: three noninflationary PBH production mechanisms, six classes of inflationary PBH models, and seven particle dark matter benchmarks, all evaluated against the same observable target. Three distinct naturalness universality classes emerge, determined entirely by the analytic structure of the abundance map rather than by the nature of the dark matter candidate. Biased-domain-wall PBHs, in their least model-dependent (free- V b ) form, have the same low sensitivity, Δ = 4.5 , as off-resonance weakly interacting massive particles and freeze-in particles ( Δ = 2 ), a sensitivity that, because it is constant over the entire parameter space of the construction, also coincides trivially with its own Wilson-normalized average within that parameter space (Section Definition and conventions), an equivalence that concerns only the space over which Δ is computed and is not a naturalness statement about the construction as a whole; a further reduction to Δ = 2 is possible only under the additional, independently motivated but not required, assumption that the domain-wall bias is generated by Planck-suppressed operators; early matter-domination PBHs occupy an intermediate tier alongside coannihilating weakly interacting massive particles (WIMPs), unified by a structural identity in which the sensitivity measure equals the logarithm of the ratio of the formation scale to the matter–radiation equality scale; first-order phase transition PBHs, once the more accurate super-exponential collapse probability is used in place of the single-exponential approximation, instead belong to the same highly sensitive tier as resonant WIMP annihilation and single-field inflationary collapse, for a structurally distinct reason; single-field ultraslow-roll inflationary collapse is severely tuned for a distinct reason: a double exponential in which the power spectrum amplitude is itself exponentially sensitive to the inflaton potential coefficients, on top of the exponential collapse sensitivity of the abundance map. My main conclusion is that the claim that PBH dark matter is generically fine-tuned conflates the worst case with a landscape spanning every naturalness tier. The Barbieri-Giudice sensitivity computed here and the Wilson-normalized measure of Iovino and Riotto answer distinct and mutually consistent questions about the same construction, a distinction I clarify within the two-layer decomposition.

Profumo, Stefano [University of California, Santa ↗

Off-shell effective energy theory: A unified treatment of the Hubbard model from $\textit{d} = 1$ to $\textit{d = ∞}$

Here we propose an exact formalism, off-shell effective energy theory (OET), which provides a thermodynamic description of a generic quantum Hamiltonian. The OET is based on a partitioning of the Hamiltonian and a corresponding density matrix ansatz constructed from an off-shell extension of the equilibrium density matrix; and there are dual realizations based on a given partitioning. To approximate OET, we introduce the central point expansion (CPE), which is an expansion of the density matrix ansatz, and we renormalize the CPE using a standard expansion of the ground-state energy. In this work we showcase the OET for the one-band Hubbard model in $\textit{d}$ = 1 , 2, and $\infty$, using a partitioning between kinetic and potential energy, yielding two realizations denoted as $\mathcal{K}$ and $\mathcal{X}$. OET shows favorable agreement with exact or state-of-the-art results over all parameter space, and has a negligible computational cost. Physically, $\mathcal{K}$ describes the Fermi liquid, while $\mathcal{X}$ gives an analogous description of both the Luttinger liquid and the Mott insulator. Our approach should find broad applicability in lattice model Hamiltonians, in addition to real materials systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Spatial correlation functions of paracrystals with radial symmetry

Here, we develop a phenomenological model to describe the structure of radially symmetric paracrystals whose long-range order are destroyed by propagation of particle fluctuations. General expressions are derived for the spatial correlation functions in one-, two-, and three-dimensional spaces. And the spatial correlation in paracrystals in reciprocal space is further discussed and clarified. The developed method can be used to quantitatively analyze the microstructure of paracrystalline materials in both real and reciprocal spaces via scattering experiments and computer simulations.

36 MATERIALS SCIENCE↗

An Online Approach to Solve the Dynamic Vehicle Routing Problem with Stochastic Trip Requests for Paratransit Services

Many transit agencies operating paratransit and microtransit services have to respond to trip requests that arrive in real-time, which entails solving hard combinatorial and sequential decision-making problems under uncertainty. To avoid decisions that lead to significant inefficiency in the long term, vehicles should be allocated to requests by optimizing a non-myopic utility function or by batching requests together and optimizing a myopic utility function. While the former approach is typically offline, the latter can be performed online. We point out two major issues with such approaches when applied to paratransit services in practice. First, it is difficult to batch paratransit requests together as they are temporally sparse. Second, the environment in which transit agencies operate changes dynamically (e.g., traffic conditions can change over time), causing the estimates that are learned offline to become stale. To address these challenges, we propose a fully online approach to solve the dynamic vehicle routing problem (DVRP) with time windows and stochastic trip requests that is robust to changing environmental dynamics by construction. We focus on scenarios where requests are relatively sparse—our problem is motivated by applications to paratransit services. We formulate DVRP as a Markov decision process and use Monte Carlo tree search to evaluate actions for any given state. Accounting for stochastic requests while optimizing a non-myopic utility function is computationally challenging; indeed, the action space for such a problem is intractably large in practice. To tackle the large action space, we leverage the structure of the problem to design heuristics that can sample promising actions for the tree search. Our experiments using real-world data from our partner agency show that the proposed approach outperforms existing state-of-the-art approaches both in terms of performance and robustness.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

GPU-Accelerated Drug Discovery with Docking on the Summit Supercomputer: Porting, Optimization, and Application to COVID-19 Research

Protein-ligand docking is an in silico tool used to screen potential drug compounds for their ability to bind to a given protein receptor within a drug-discovery campaign. Experimental drug screening is expensive and time consuming, and it is desirable to carry out large scale docking calculations in a high-throughput manner to narrow the experimental search space. Few of the existing computational docking tools were designed with high performance computing in mind. Therefore, optimizations to maximize use of high-performance computational resources available at leadership-class computing facilities enables these facilities to be leveraged for drug discovery. Here we present the porting, optimization, and validation of the AutoDock-GPU program for the Summit supercomputer, and its application to initial compound screening efforts to target proteins of the SARS-CoV-2 virus responsible for the current COVID-19 pandemic.

LeGrand, Scott↗

Memory Analysis Tool

This tool traces all memory accesses to stack (static allocation) and heap (dynamic allocation) on a trace run in a particular hardware and then estimates execution time on given arbitrary hardware configurations for hardware design space exploration. It can also compute memory access statistics, such as reuse distance.

Sato, Kento↗

FEDIS (The “Fresnel Equations” for Diffuse radiation on Inclined photovoltaic Surfaces) [SWR-22-64]

The well-known Fresnel equations solve for the reflection and transmission of light for precise incident angles. The transmission of diffuse radiation incident on a planar or domed surface is often needed for real-world applications. Due to the complexity of the Fresnel equations, the analytical solution of the integration has hitherto been unobtainable over the last centuries. Therefore, this problem was numerically solved by integrating the angular transmittances in space often leading to substantial computing burden and bias in the results. To efficiently estimate the solar energy resource for a glass-covered photovoltaic (PV) module, we derive an analytical solution of diffuse transmission based on the rigorous integration of an alternate form of the Fresnel equations. The approach leads to a simple yet accurate relative transmittance model that reconciles the solar energy sensed by pyranometers and PV panels. With limited and clearly stated approximations, the complex mathematical derivation resulted in an elegant solution.

Xie, Yu↗

Permanent Magnets Featuring Heavy Main Group Elements for Magnetic Anisotropy

Permanent magnets are the functional component of electric motors and generators found in numerous renewable energy applications. To improve energy conversion in such applications, we require fundamentally new magnets that generate higher magnetic flux per volume while retaining the properties conferred by rare-earth elements incorporated into current technologies. We hypothesize that by engendering a covalent interaction between two elements, we can access a new regime of magnetic materials where the two components of a magnetic moment—spin and orbital angular momentum—come from two separate atoms to form a complete magnetic moment. Our previous research utilized high-pressure conditions to discover two new candidate materials ideal for assessing this hypothesis. The first material, FeBi 2 , enables the study of an unprecedented solid state metal-metal bonding interaction. The second, MnBi 2 represents the second member of the promising Mn–Bi family known for its magnetic properties. Importantly, MnBi 2 is isostructural to FeBi 2 . Together these chemically simple but magnetically rich materials provide an elegant platform for elucidating fundamental design principals of magnetic anisotropy while inspiring the synthesis of new magnetic materials. More generally, solid-state chemistry remains a synthetic black box. To fully harness the potential of such new materials, it is vital to create a window into both the structure and the properties of our high-pressure materials. One missing area of importance to energy science is high-pressure magnetometry. To this end, we assessed the magnetic structure of these materials through two objectives: (1) we performed magnetometry measurements at high-pressures, including performing high-pressure X-ray Magnetic Circular Dichroism (XMCD) experiments on these materials that are only synthesized at high-pressure and cannot be recovered to ambient conditions and (2) we targeted the recovery of these materials using dynamic compression approaches, starting with laser shock techniques. Future directions include expanding our focus to ternary phase space, employing high-throughput computational approaches to direct our search for magnetic materials across the Periodic Table. This award was originated at Northwestern, and reopened as a new grant at MIT. This is the close out for Northwestern.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗