Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational physics applications”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Storage-Induced Collapse of Lignin Macromolecular Structure and Its Impacts on the Biorefinery

Lignin plays a vital role in the economics of biorefineries, serving as a source of process energy and a feedstock for sustainable fuels and chemical production. While understanding lignin’s chemical composition is crucial, emerging evidence suggests that a more comprehensive understanding of its macromolecular structure is critical to explaining its complex behavior in the biorefinery. This study investigated the partial collapse of the lignin network in corn stover feedstock after harvest and storage as a result of the microbial digestion of hemicellulose. Fluorescence microscopy was used to detect the collapse of lignin in terms of lignin’s inter-molecular interaction and the re-orientation of lignin’s chromophores, by the changes in lignin’s fluorescence lifetime, anisotropy, and the number of effective emitters. With minimal sample perturbation, our in-situ microscopic results revealed lignin's coil-globule transition phenomena, which was only previously predicted by molecular dynamics modeling extracted lignin in solvent. This collapse of lignin macromolecular structure was confirmed by results from NMR, IR, Raman, and powder X-ray diffraction. We also investigated the impact of this storage-induced collapse on the downstream biorefinery processes. Our study revealed that the two major approaches for lignin valorization in the lignin-first biorefinery model, namely monomer extraction and milled wood lignin extraction, were negatively impacted by the lignin collapse. As changes during storage are a source of feedstock variability, our study highlights the importance of understanding the effect of feedstock handling on biorefinery operations and economics.

09 BIOMASS FUELS↗
DOE OSTI2021

Machine Learning Assisted Safety Modeling and Analysis of Advanced Reactors

With the advances in computational power and numerical methods, analysts can now rely on first-principle simulations to predict ultra-fine details in a variety of applications. Advances in machine learning (ML) have produced algorithms that can now learn high-level abstractions via hierarchical models. This project aims to leverage advances in ML techniques and the available high-resolution simulation data to develop a novel modeling and simulation (M\&S) methodology for reactor safety analysis. While application-agnostic ML techniques are available, complex physics constraints need to be incorporated into ML techniques to build ML-based closures for computationally efficient predictive simulations. This project intends to develop a physics-guided data-driven multi-scale methodology for M\&S of advanced reactors. The project focuses on thermal fluid (T/F) phenomena, which play major roles in advanced reactor safety. Specifically, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor (SFR). The model has a coarse-mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data with Reynolds-averaged Navier-Stokes (RANS) turbulence model to ensure accuracy. Three different neural networks are developed and tested for loss-of-flow transients in the hot pool of SFR, i.e. the densely connected convolutional neural network (DCNN), long-short-term-memory network based on proper orthogonal decomposition (POD-LSTM), and the DCNN informed by LSTM (DCNN-LSTM). The performances of these three neural networks are evaluated based on baseline models. The DCNN-LSTM model has been chosen for further hyperparameter optimization. Furthermore, based on a simplified two-dimensional case, uncertainty quantification (UQ) of the developed ML-based closure are investigated with three methods, i.e. Monte Carlo dropout, deep ensemble, and Bayesian neural network. The developed ML-based turbulent viscosity closure relation based on deep ensemble is then integrated into the system analysis module SAM and serves as a term in the conservation equations. Such a SAM-ML based procedure guarantees that the obtained results are consistent with the physical constraints of the thermal-fluid system. The SAM-ML simulation on the same loss-of-flow transient showed comparable accuracy with the CFD simulation but with a much coarser mesh setup. Last but not least, the ML-based closure improvement with the support of higher-fidelity data from large eddy simulation (LES) is discussed. As a first step towards this direction, a baseline LES simulation is performed to obtain comparable data with RANS results. Based on the early results, future investigation on further improving the ML-based closure is discussed. We believe the developed approach that combines scientific machine learning with nuclear system analysis code can benefit the advanced reactor community as more accurate safety analyses will better characterize reactor safety margins and reduce licensing efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

LAN: A Materials Notation for Two-Dimensional Layered Assemblies

Two-dimensional (2D) layered materials offer intriguing possibilities for novel physics and applications. Before any attempt at exploring the materials space in a systematic fashion, or combining insights from theory, computation, and experiment, a formal description of information about an assembly of arbitrary composition is required. Here, we introduce a domain-generic notation that is used to describe the space of 2D layered materials from monolayers to twisted assemblies of arbitrary composition, existent or not yet fabricated. The notation corresponds to a theoretical materials concept of stepwise assembly of layered structures using a sequence of rotation, vertical stacking, and other operations on individual 2D layers. Its scope is demonstrated with several example structures using common single-layer materials as building blocks. This work overall aims to contribute to the systematic codification, capture, and transfer of materials knowledge in the area of 2D layered materials.

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

Anti-$\mathscr{PT}$-symmetric qubit: Decoherence and entanglement entropy

We investigate the dynamics of a general two-level based anti-parity-time (anti-$\mathscr{PT}$)-symmetric qubit and study its decoherence as well as entanglement entropy properties. We compare our findings with that of the corresponding parity-time ($\mathscr{PT}$)-symmetric and Hermitian qubits. To begin, we consider the time-dependent Dyson map to find the exact analytical dynamics for a general non-Hermitian qubit system weakly coupled with a thermal bath for pure dephasing, before specializing it to the case of a general anti-$\mathscr{PT}$-symmetric qubit. Basing the comparison under the same coupling strength or increasing the non-Hermiticity, we observe that the decoherence function and entanglement entropy of the anti-$\mathscr{PT}$-symmetric qubit decays and grows more slowly, respectively, compared to the $\mathscr{PT}$-symmetric and Hermitian qubits. Similarly, the corresponding variance and area of Fisher information are much higher compared to the $\mathscr{PT}$-symmetric and Hermitian qubits. These results demonstrate that anti-$\mathscr{PT}$-symmetric qubits may be better suited for quantum computing and quantum information processing applications than conventional Hermitian or even $\mathscr{PT}$-symmetric qubits.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

WarpX: Delivering Unprecedented Particle-in-Cell Simulation Capability

Since 1939, more than 30 percent of all Nobel prizes in physics and four of the past 14 prizes in chemistry have been awarded for work with particle accelerators. These tools are now used to treat cancer and produce medical supplies, support research in applied fields such as pharmaceuticals, create microcircuits, and sterilize food—and the applications list continues to grow. However, the size of the facilities needed to house these machines—which are often dozens of kilometers long—has limited their future potential. Thus, researchers are motivated to create novel particle accelerator designs that provide new imaging capability but have a much smaller footprint. The Exascale Computing Project’s WarpX application captures the full complexity of the acceleration processes within plasma-based design, supporting the rapid and cost-efficient development of tens of thousands of new particle accelerators for various applications from improving human health to nuclear fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Solving reaction dynamics with quantum computing algorithms

The description of quantum many-body dynamics is extremely challenging on classical computers, as it can involve many degrees of freedom. However, the time evolution of quantum states is a natural application for quantum computers that are designed to efficiently perform unitary transformations. Here, in this paper, we study quantum algorithms for response functions, relevant for describing different reactions governed by linear response. We focus on nuclear-physics applications and consider a qubit-efficient mapping on the lattice, which can efficiently represent the large volumes required for realistic scattering simulations. For the case of a contact interaction, we develop an algorithm for time evolution based on the Trotter approximation that scales logarithmically with the lattice size and is combined with quantum phase estimation. We eventually focus on the nuclear two-body system and a typical response function relevant for electron scattering as an example. We also investigate ground-state preparation and examine the total circuit depth required for a realistic calculation and the hardware noise level required to interpret the signal.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Classroom aerosol dispersion modeling: experimental assessment of a low-cost flow simulation tool

The purpose of this study was to assess the utility of a low-cost flow simulation tool for an indoor air modeling application by comparing its outputs with the results of a physical experiment, as well as those from a more advanced computational fluid dynamics (CFD) software package. In this study, five aerosol dispersion tests were performed in two different classrooms by releasing a CO 2 tracer gas from six student locations. Resultant steady-state concentrations were monitored at 13 locations around the periphery of the room. Subsequently, the experiments were modeled using both a low-cost tool (SolidWorks Flow Simulation) and a more sophisticated tool (STAR-CCM+). Models were evaluated based on their ability to predict the experimentally measured concentrations at the 13 monitoring locations by calculating four performance parameters commonly used in the evaluation of dispersion models: fractional mean bias (FB), normalized mean-square error (NMSE), fraction of predicted value within a factor of two (FAC2), and normalized absolute difference (NAD). The more sophisticated model performed better in 15 of the 20 possible cases (five tests at four parameters each), with parameters meeting acceptance criteria in 19 of 20 cases. However, the lower-cost tool was only slightly worse, with parameters meeting acceptance criteria in 18 of 20 cases, and it performed better than the other tool in 3 of 20 cases. Because it provides useful results at a fraction of the monetary and training cost and is already widely accessible to many institutions, such a tool may be worthwhile for many indoor aerosol dispersion applications, especially for students or researchers just beginning CFD modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Digital Modeling on Large Kernel Metamaterial Neural Network

Deep neural networks (DNNs) utilized recently are physically deployed with computational units (e.g., CPUs and GPUs). Such a design might lead to a heavy computational burden, significant latency, and intensive power consumption, which are critical limitations in applications such as Internet of Things (IoT), edge computing, and usage of drones. Recent advances in optical computational units (e.g., metamaterial) have shed light on energy-free and light-speed neural networks. However, the digital design of the metamaterial neural network (MNN) is fundamentally limited by its physical limitations, such as precision, noise, and bandwidth during fabrication. Moreover, the unique advantages of MNN’s (e.g., light-speed computation) are not fully explored via standard 3×3 convolution kernels. In this paper, we propose a novel large kernel metamaterial neural network (LMNN) that maximizes the digital capacity of the state-of-the-art (SOTA) MNN with model re-parametrization and network compression, while also considering the optical limitation explicitly. The new digital learning scheme can maximize the learning capacity of MNN while modeling the physical restrictions of meta-optics. With the proposed LMNN, the computation cost of the convolutional front-end can be offloaded to fabricated optical hardware. The experimental results on two publicly available datasets demonstrate that the optimized hybrid design improved classification accuracy while reducing computational latency. In conclusion, the development of the proposed LMNN is a promising step towards the ultimate goal of energy-free and light-speed AI.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Quantum Computation of Dynamical Quantum Phase Transitions and Entanglement Tomography in a Lattice Gauge Theory

Strongly coupled gauge theories far from equilibrium may exhibit unique features that could illuminate the physics of the early universe and of hadron and ion colliders. Studying real-time phenomena has proven challenging with classical-simulation methods but is a natural application of quantum simulation. To demonstrate this prospect, we quantum compute nonequal-time correlation functions and perform entanglement tomography of nonequilibrium states of a simple lattice gauge theory, the Schwinger model, using a trapped-ion quantum computer by IonQ Inc. As an ideal target for near-term devices, a recently predicted dynamical quantum phase transition in this model is studied by preparing, quenching, and tracking the subsequent nonequilibrium dynamics in three ways: (i) overlap echos signaling dynamical transitions, (ii) nonequal-time correlation functions with an underlying topological nature, and (iii) the entanglement structure of nonequilibrium states, including entanglement Hamiltonians. These results constitute the first observation of a dynamical quantum phase transition in a lattice gauge theory on a quantum computer and are a first step toward investigating topological phenomena in nuclear and high-energy physics using quantum technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Emergent Properties from Three-Dimensional Assemblies of (Nano)particles in Confined Spaces

The assembly of (nano)particles into compact hierarchical structures yields emergent properties not found in the individual constituents. The formation of these structures relies on a profound knowledge of the nanoscale interactions between (nano)particles, which are often designed by researchers aided by computational studies. These interactions have an effect when the (nano)particles are brought into close proximity, yet relying only on diffusion to reach these closer distances may be inefficient. Recently, physical confinement has emerged as an efficient methodology to increase the volume fraction of (nano)particles, rapidly accelerating the time scale of assembly. Specifically, the high surface area of droplets of one immiscible fluid into another facilitates the controlled removal of the dispersed phase, resulting in spherical, often ordered, (nano)particle assemblies. In this review, we discuss the design strategies, computational approaches, and assembly methods for (nano)particles in confined spaces and the emergent properties therein, such as trigger-directed assembly, lasing behavior, and structural photonic color. Finally, we provide a brief outlook on the current challenges, both experimental and computational, and farther afield application possibilities.

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

Evaluating Application Characteristics for GPU Portability Layer Selection

GPUs have become the dominant source of computing power for high performance computing and are increasingly being used across the High Energy Physics computing landscape for a wide variety of tasks. Though NVIDIA is currently the main provider of GPUs, AMD and Intel are rapidly increasing their market share. As a result, programming using a vendor-specific language such as CUDA can significantly reduce deployment choices. There are a number of portability layers such as Kokkos, Alpaka, SYCL, OpenMP and std::par that permit execution on a broad range of GPU and CPU architectures, significantly increasing the flexibility of application programmers. However, each of these portability layers has its own characteristics, performing better at some tasks and worse at others, or placing limitations on aspects of the application. In this presentation, we report on a study of application and kernel characteristics that can influence the choice of a portability layer and show how each layer handles these characteristics. We have analyzed representative heterogeneous applications from CMS (patatrack and p2r), DUNE (Wire-Cell Toolkit), and ATLAS (FastCaloSim) to identify key application characteristics that have different behaviors for the various portability technologies. Using these results, developers can make more informed decisions on which GPU portability technology is best suited to their application.

Atif, Mohammad [Brookhaven]↗
DOE OSTI2022

Computational Design of a Simple Flyer Plate Launcher

The response of materials to shock loading is important to understand for a variety of applications. When shock physics emerged during and after WWII, direct explosive loading or explosively driven plate impact was the primary tool for these studies. Subsequent decades have seen the widespread use of large caliber guns for plate impact studies, laser-shock facilities and pulsed power facilities. INL currently lacks a gun suitable for plate impact or explosives casting and machining facilities; however, it does possess explosives use and handling capabilities. An option for performing plate impact experiments was needed, therefore continuum scale models were utilized to explore a few simple donor-acceptor explosive plane wave lens designs, one of which could be hand packed with plastic explosives to launch flyer plates. 2D simulations were performed to study different geometries in an effort to minimize the difference in shock arrival across the central portion of a small copper flyer plate. A shock wave arrival time difference under 50 ns across 50 percent of the center of the flyer was achieved with a few designs. This work summarizes the computational models and results.

36 MATERIALS SCIENCE↗
DOE OSTI2026

New approaches to Bayesian uncertainty quantification for Nuclear Science (Final Technical Report)

Inverse problems play a central role in experimentation and theory/data comparisons for many areas of modern Nuclear Physics (NP) and High-Energy Physics (HEP). Bayes’s Theorem is a powerful tool for solving Inverse Problems, providing conceptually transparent and unbiased constraints on theoretical parameters and their uncertainties (“Bayesian Inference”) and enabling the quantification of agreement or tension between models and data. However, analyses based on Bayesian Inference are often challenging for NP and HEP applications, either because of the large number of parameters in the problem, the high computational cost, or both. We propose a multi-institutional collaboration to develop and deploy novel Bayesian analysis tools that advance the scientific scope of a broad range of current and future NP experiments. This project brings together NP domain scientists working on several high-profile NP projects for which new, high-performance Bayesian Uncertainty Quantification (“Bayesian UQ”) methods are essential to carry out the science, and data scientists who are developing state-of-the-art methods applicable to these problems. The NP projects in this proposal comprise measurements of the mass and fundamental nature of the neutrino; study of the Quark-Gluon Plasma that filled the early universe; and mapping of natural and anthropogenic radiation environments. While these NP projects have very different scientific goals, with datasets and analysis approaches that differ significantly, they share common requirements for improving computationally intensive Bayesian analyses using advanced Machine Learning algorithms and will benefit strongly from a coherent effort to develop general solutions. This proposal brings together these projects and forefront ML-based data science algorithms to develop such general solutions. The methods developed in this project will also be more widely applicable, thereby advancing science in the larger Nuclear Physics portfolio.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Beyond Single-Reference Fixed-Node Approximation in Ab Initio Diffusion Monte Carlo Using Antisymmetrized Geminal Power Applied to Systems with Hundreds of Electrons

Diffusion Monte Carlo (DMC) is an exact technique to project out the ground state (GS) of a Hamiltonian. Since the GS is always bosonic, in Fermionic systems, the projection needs to be carried out while imposing antisymmetric constraints, which is a nondeterministic polynomial hard problem. In practice, therefore, the application of DMC on electronic structure problems is made by employing the fixed-node (FN) approximation, consisting of performing DMC with the constraint of having a fixed, predefined nodal surface. How do we get the nodal surface? The typical approach, applied in systems having up to hundreds or even thousands of electrons, is to obtain the nodal surface from a preliminary mean-field approach (typically, a density functional theory calculation) used to obtain a single Slater determinant. This is known as single reference. In this paper, we propose a new approach, applicable to systems as large as the C 60 fullerene, which improves the nodes by going beyond the single reference. In practice, we employ an implicitly multireference ansatz (antisymmetrized geminal power wave function constraint with molecular orbitals), initialized on the preliminary mean-field approach, which is relaxed by optimizing a few parameters of the wave function determining the nodal surface by minimizing the FN-DMC energy. We highlight the improvements of the proposed approach over the standard single-reference method on several examples and, where feasible, the computational gain over the standard multireference ansatz, which makes the methods applicable to large systems. We also show that physical properties relying on relative energies, such as binding energies, are affordable and reliable within the proposed scheme.

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

Tuning of oscillation modes by controlling dimensionality of spin structures

Abstract Harmonic oscillation of spin structures is a physical phenomenon that offers great potential for applications in nanotechnologies such as nano-oscillators and bio-inspired computing. The effective tuning of oscillations over wide frequency ranges within a single ferromagnetic nanoelement is a prerequisite to realize oscillation-based nanodevices, but it has not been addressed experimentally or theoretically. Here, utilizing a vortex core structure, one of spin structures, we report a drastic change of oscillation modes over the frequency range from MHz to sub-GHz in a 100 nm-thick permalloy circular disk. Oscillation mode was found to considerably depend on the shape and dimension of the vortex core structure and various oscillation modes over a wide range of frequencies appeared with dimensional change in the vortex core structure. This work demonstrates that oscillation modes of the vortex core structure can be effectively tuned and opens a way to apply spin structures to oscillation-based technology.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dense hydrogen layers for high performance MagLIF

We report that Magnetized Liner Inertial Fusion (MagLIF) experiments driven by the Z machine produce >10 13 deuterium-deuterium fusion reactions. Simulations indicate high yields and gains (1000) with increased current and deuterium-tritium layers for burn propagation. Such a coating also isolates the metal liner from the gaseous fuel, which should reduce mixing of liner material into the fuel. However, the vapor density at the triple point is only 0.3 kg/m3, which is not high enough for MagLIF operation. We present two solutions to this problem. First, a fuel wetted low-density plastic foam can be used to form a layer on the inside of the liner. The desired vapor density can be obtained by controlling the temperature. This does however introduce carbon into the layer which will enhance radiation losses. Simulations indicate that this wetted foam layer can significantly contribute to the fusion yield when the foam density is less than 35 kg/m 3 . Second, we show that a pure frozen fuel layer can first be formed on the inside of the liner and then low temperature gaseous fuel can be introduced just before the implosion without melting a significant amount of the ice layer. This approach is the most promising for MagLIF to produce high yield and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Fast inversion, preconditioned quantum linear system solvers, fast Green's-function computation, and fast evaluation of matrix functions

Preconditioning is the most widely used and effective way for treating ill-conditioned linear systems in the context of classical iterative linear system solvers. We introduce a quantum primitive called fast inversion, which can be used as a preconditioner for solving quantum linear systems. The key idea of fast inversion is to directly block encode a matrix inverse through a quantum circuit implementing the inversion of eigenvalues via classical arithmetics. We demonstrate the application of preconditioned linear system solvers for computing single-particle Green's functions of quantum many-body systems, which are widely used in quantum physics, chemistry, and materials science. We analyze the complexities in three scenarios: the Hubbard model, the quantum many-body Hamiltonian in the plane-wave-dual basis, and the Schwinger model. We also provide a method for performing Green's function calculation in second quantization within a fixed-particle manifold and note that this approach may be valuable for simulation more broadly. Aside from solving linear systems, fast inversion also allows us to develop fast algorithms for computing matrix functions, such as the efficient preparation of Gibbs states. Furthermore, we introduce two efficient approaches for such a task, based on the contour-integral formulation and the inverse transform, respectively.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

MrHyDE v.1.0

SAND2024-01324O MrHyDE, which stands for Multi-resolution Hybridized Differential Equations, is a general-purpose C++ package for the solution of coupled multiphysics and multiscale systems on massively parallel computing systems. MrHyDE is designed to enable moving beyond forward simulation for multiscale applications which includes optimization, control, uncertainty quantification, and stochastic inversion. The framework provides interfaces to several packages within the Trilinos framework and leverages automatic differentiation to enable adjoint capabilities for large-scale, gradient-based optimization. MrHyDE provides automated multiscale capabilities through a subgrid model interface and multiscale Dirichlet-to-Neumann maps. For extreme-scale applications, MrHyDE provides in situ data-compression algorithms to reduce memory requirements while maintaining performance. MrHyDE is a general-purpose, computational framework for the solution of multiscale and multiphysics applications. It uses a combination of structure-preserving, physics-compatible discretizations, fully implicit methods, multi-resolution schemes, or fully explicit methods. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

SciDAC↗