Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computational Methods”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Two-particle correlation functions in cluster perturbation theory: Hubbard spin susceptibilities

Cluster perturbation theory (CPT) is a computationally economic method commonly used to estimate the momentum- and energy-resolved single-particle Green's function. It has been used extensively in direct comparisons with experiments that effectively measure the single-particle Green's function, e.g., angle-resolved photoemission spectroscopy. However, many experimental observables are given by two-particle correlation functions. CPT can be extended to compute two-particle correlation functions by approximately solving the Bethe-Salpeter equation. We implement this method and focus on the transverse spin susceptibility, measurable via inelastic neutron scattering or with optical probes of atomic gases in optical lattices. Finally, we benchmark the method with the one-dimensional Fermi-Hubbard model by comparing with known results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Deep Learning for Full Waveform Inversion of Elastic Active-Source Seismic Data to Estimate P-Wave Velocity Models

Seismic imaging methods are critical for Global Security and Energy & Homeland Security missions and activities that rely on subsurface characterization, but traditional methods remain computationally expensive and require significant labor hours and expertise to execute. Within the past few years, machine learning (ML), namely deep learning (DL), has been used to develop data-driven end-to-end full waveform inversion (FWI) methods to estimate 2D P-wave velocity (Vp) models in a fraction of the time as conventional FWI. These methods, however, are trained on simplistic acoustic wave seismic data and Vp models that are not realistic nor representative of real-world observations, leaving a large gap between the state-of-the-art and deployable, feasible, and practical DL FWI methods. Here, we generate a synthetic active-source, 3D, elastic wave seismic data set and a variety of Vp models with realistic geologic structure for training DL FWI methods. We evaluate six different methods that have performed well for acoustic DL FWI or medical imaging tasks using our more realistic dataset. We find that these six trained models do not match the performance of published acoustic end-to-end DL FWI methods, indicating more training data may be needed, physics may need to be incorporated to achieve good accuracy at the sacrifice of the end-to-end advantage, and/or novel methods need to be developed to enable end-to-end DL FWI methods to perform well for real-world seismic data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Green's function methods for multiphysics simulations (Final Report)

Green's functions are important tools for analyzing mathematical properties of partial differential equations (PDEs), and for numerically solving PDEs, especially when equations for the same operator but with multiple right hand sides need to be solved simultaneously. This proposal aims at developing efficient and accurate numerical methods for computing Green's functions, which can be used to tackle a challenging question in multiphysics simulation of DOE-mission science problems: how to couple quantum physics with classical physics. The key mathematical difficulty is properly formulate a "boundary condition" for the region described by quantum physics, and conventional approaches often model such boundary conditions in an empirical way. The proposed Green's function methods use the Dirichlet-to-Neumann map to formulate a boundary condition that is non-empirical and can couple the quantum and classical regions in an in principle exact way. The key ingredient of the new methods is to construct the Dirichlet-to-Neumann map in an efficient, accurate and versatile manner. The new methods have provably low complexity and are ideally suited for massively parallel and emerging many-core computational systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Computational modeling of phononic pseudocrystal isolators

Methods for the efficient computational prediction of the performance of phononic pseudocrystals (structured materials capable of blocking extraordinary ranges of frequency) in COMSOL and other comparable finite element method codes are set forth. These methods include boundary conditions that make possible halving the size of the computational domain. Also included is an introduction of elastic energy density methods for assessing the extinction of elastic waves within the patterned region.

Swift, Stephen Hales↗
DOE OSTI2024

Bootstrapping gauge theories

We consider asymptotically free gauge theories with gauge group S U ( N c ) and N f quarks with mass m q ≪ Λ QCD that undergo chiral symmetry breaking and confinement. We propose a bootstrap method to compute the S matrix of the pseudo-Goldstone bosons (pions) that dominate the low-energy physics. For the important case of N c = 3 , N f = 2 , a numerical implementation of the method gives the phase shifts of the S 0 , P 1 and S 2 waves in good agreement with experimental results. The method incorporates gauge theory information ( N c , N f , m q , Λ QCD ) by using the form-factor bootstrap recently proposed by Karateev, Kuhn and Penedones together with a finite energy version of the Shifman-Vainshtein-Zakharov (SVZ) sum rules. At low energy we impose constraints from chiral symmetry breaking. The only low-energy numerical inputs are the pion mass m π and the quark and gluon condensates. Published by the American Physical Society 2024

He, Yifei (ORCID:0000000213666157)↗
DOE OSTI2025

Computational alchemy clarifies origins of alloy strengthening

Solid solution strengthening (SSS) is widely used to enhance mechanical properties of metals. Originally developed for dilute alloys, classical SSS theories are presently challenged by the rise of complex concentrated alloys (CCA) with nearly equiatomic compositions. Here, we propose and develop a method of “computational alchemy” in which interatomic interactions are modified to systematically vary two key physical parameters defining SSS - atomic size misfit and elastic stiffness misfit - over a maximally wide range of two misfits. The resulting alchemical alloys are subjected to massive (~10 8 atoms) molecular dynamics (MD) simulations reproducing full complexity of plastic strength response. At variance with prevailing views, stiffness misfit is observed to contribute to SSS on par if not more than size misfit. Furthermore, depending on exactly how two misfits are combined, they result in synergistic (amplification) or antagonistic (compensation) effect on alloy strengthening. Unlike real CCAs in which each component element comes with its own specific size and stiffness, our alchemical model alloys span the space of two misfits continuously revealing trends in alloy strengthening unrecognized so far. Our study demonstrates unique value of intentionally unrealistic models for gaining deep physical insights into material behaviors that are difficult to reveal otherwise.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Efficient near-field ptychography reconstruction using the Hessian operator

X-ray ptychography is a powerful and robust coherent imaging method providing access to the complex object and probe (illumination). Ptychography reconstruction is typically performed using first-order methods due to their computational efficiency. Higher-order methods, while potentially more accurate, are often prohibitively expensive in terms of computation. In this study, we present a mathematical framework for reconstruction using second-order information derived from an efficient computation of the bilinear Hessian and Hessian operator. The formulation is provided for Gaussian-based models, enabling the simultaneous reconstruction of the object, probe, and object positions. Synthetic data tests, along with experimental near-field ptychography data processing, demonstrate a ten-fold reduction in computation time compared to first-order methods. The derived formulas for computing the Hessians, along with the strategies for incorporating them into optimization schemes, are well-structured and easily adaptable to various ptychography problem formulations.

Carlsson, Marcus [Lund Univ. (Sweden)] (ORCID:0000↗
DOE OSTI2020

Conditional Generative Adversarial Networks for Solving Heat Transfer Problems

Generative Adversarial Networks (GANs) have been used as a deep learning approach to solving physics and engineering problems. Using deep learning for these problems is attractive in that reasonably accurate models can be inferred from only raw data, eliminating the need to define the exact physical equations governing a problem. We expand on previous work using GANs to generate steady-state solutions to the two-dimensional heat equation. Using a basic conditional GAN (cGAN), we generate accurate solutions for rectangular domains conditioned on four edge boundary conditions (MAE < 0.5%). For finding steady-state solutions over arbitrary two-dimensional domains (not constrained to rectangles), we use a cGAN designed for image-to-image translation. We train this GAN on various types of geometric domains (circles, squares, triangles, shapes with one circular or rectangular hole), achieving accurate results on test data made up of geometries similar to those in training (MAE < 1%). For both of these GANs, we experiment with different loss function terms, showing that a term using the gradients of solution images significantly improves the basic cGAN but not the image-to-image GAN. Lastly, we show that the image-to-image GAN performs poorly when applied to two-dimensional geometries that vary in structure from training data (MAE < 8% for shapes with multiple holes or different shaped holes). This demonstrates the cGAN's lack of generalizability. While the cGAN is an accurate and computationally efficient method when trained and tested on similarly structured data, it is a much less reliable method when applied to data that is slightly different in structure from the training data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Prediction for Pressure Differential Across HEPA Filter Media Based on Media Characteristics and Particle Size Distribution

A new method for predicting the pressure drop across High Efficiency Particulate Air (HEPA) filter media is proposed based upon mass deposited onto the filter and known physical characteristics of the filter media. Detailed are the methods used in conjunction with current filter loading models to predict the pressure drop, as well as tests conducted to validate the methods. The benefit of a prediction model for practical use lies in the manufacturing and service life of nuclear grade HEPA filters. HEPA filters for use in nuclear facilities have a prescribed expiration date and a maximum allowable operating pressure drop. Therefore, the ability to predict the pressure drop across a filter and relate it to an expected length of service time can enable a reduction in wasted filters. This will allow for more informed decisions to be made based upon the dictated life cycle of the filters. Also, understanding and predicting how the pressure drop of HEPA filter media behaves as a function of physical characteristics and loaded mass can assist in future design and manufacturing of filter media. Currently, most existing pressure drop models are either computationally based or analytical methods relying on data gathered during tests. Neither are practical for prediction; the computational methods are difficult to implement, and the current analytical methods are more useful as tools for analysis. The current analytical model developed by Bergman et al. is based upon the pressure drop from each media fiber and modeling the deposited particles as newly formed fibers. Using Bergman's model and the media properties, air properties, and known aerosol particle size distribution, a pressure drop curve as a function of loaded mass can be created. This curve implies initial loading in the depth of the filter media with a transition to pure surface loading on the filter face, and a sensitivity to an evolving particle size distribution as mass continues to load onto the filter. To implement the predictive model, knowledge of the mean filter fiber diameter and porosity of the filter media is required. Traditionally, the mean fiber diameter is calculated as an effective diameter from prior media testing, however, in this study a Scanning Electron Microscope (SEM) technique is used to acquire this variable. Validation of the predictive model is provided by flat sheet media tests under a controlled environment with a measured particle size distribution of the challenge aerosol and shows reasonable preliminary agreement.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Code interoperability extends the scope of quantum simulations

The functionality of many materials is critically dependent on the integration of dissimilar components and on the interfaces that arise between them. The description of such heterogeneous components requires the development and deployment of first principles methods, coupled to appropriate dynamical descriptions of matter and advanced sampling techniques, in order to capture all the relevant length and time scales of importance to the materials’ performance. It is thus essential to build simple, streamlined computational schemes for the prediction and design of multiple properties of broad classes of materials, by developing interoperable codes which can be efficiently coupled to each other to perform complex tasks. We discuss the use of interoperable codes to simulate the structural and spectroscopic characterization of materials, including chemical reactions for catalysis, the description of defects for quantum information science, and heat and charge transport.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Advances in geophysical forensic event monitoring

Forensic analysis of man-made, non-nuclear events (such as industrial accidents, explosion experiments and mine collapses) has become more frequent and detailed owing to advancements in geophysical monitoring. Here, in this Technical Review, we demonstrate how geophysical forensic monitoring using seismic, infrasound and hydroacoustic recordings provides insights on events in the solid earth, atmosphere and underwater. Advanced techniques, including machine-learning-based models, have been developed to detect, identify and investigate these events, providing information on location, subevents, sources and explosive yield. The increase in data availability, application of advanced methods and computation and the growth of multitechnology approaches have increased the accuracy of forensic event analysis and enabled more realistic characterization of uncertainties. For example, the 2020 Beirut explosion in Lebanon demonstrated that various seismic, acoustic and other methods could be used to estimate explosive yield (and yield uncertainties) of about 1 ktonne, providing confidence in the application of these methods to smaller events where data are available. However, forensic investigations remain largely limited to known events with identified sources. Increased access to data, sophisticated analysis methods and high-resolution earth models will improve forensic event analysis further, enabling civil and scientific applications, such as localization in the search for the lost ARA San Juan submarine.

geophysics↗
DOE OSTI2020

Time-sliced quantum circuit partitioning for modular architectures

Current quantum computer designs will not scale. To scale beyond small prototypes, quantum architectures will likely adopt a modular approach with clusters of tightly connected quantum bits and sparser connections between clusters. We exploit this clustering and the statically-known control flow of quantum programs to create tractable partitioning heuristics which map quantum circuits to modular physical machines one time slice at a time. Specifically, we create optimized mappings for each time slice, accounting for the cost to move data from the previous time slice and using a tunable lookahead scheme to reduce the cost to move to future time slices. We compare our approach to a traditional statically-mapped, owner-computes model. Our results show strict improvement over the static mapping baseline. Here, we reduce the non-local communication overhead by 89.8% in the best case and by 60.9% on average. Our techniques, unlike many exact solver methods, are computationally tractable.

Baker, Jonathan M.↗
DOE OSTI2022

An Improved Method for Coupling Hydrodynamics with Astrophysical Reaction Networks

Reacting astrophysical flows can be challenging to model, because of the difficulty in accurately coupling hydrodynamics and reactions. This can be particularly acute during explosive burning or at high temperatures where nuclear statistical equilibrium is established. We develop a new approach, based on the ideas of spectral deferred corrections (SDC) coupling of explicit hydrodynamics and stiff reaction sources as an alternative to operator splitting, that is simpler than the more comprehensive SDC approach we demonstrated previously. We apply the new method to a double-detonation problem with a moderately sized astrophysical nuclear reaction network and explore the time step size and reaction network tolerances, to show that the simplified-SDC approach provides improved coupling with decreased computational expense compared to traditional Strang operator splitting. This is all done in the framework of the Castro hydrodynamics code, and all algorithm implementations are freely available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Astromers: Nuclear Isomers in Astrophysics

In this work, we develop a method to compute thermally mediated transition rates between the ground state and long-lived isomers in nuclei. We also establish criteria delimiting a thermalization temperature above which a nucleus may be considered a single species and below which it must be treated as two separate species: a ground-state species and an astrophysical isomer ("astromer") species. Below the thermalization temperature, the destruction rates dominate the internal transition rates between the ground state and the isomer. If the destruction rates also differ greatly from one another, the nuclear levels fall out of or fail to reach thermal equilibrium. Without thermal equilibrium, there may not be a safe assumption about the distribution of occupation probability among the nuclear levels when computing nuclear reaction rates. In these conditions, the isomer has astrophysical consequences and should be treated as a separate astromer species which evolves separately from the ground state in a nucleosynthesis network. We apply our transition-rate methods and perform sensitivity studies on a few well-known astromers. We also study transitions in several other isomers of likely astrophysical interest.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Real-time spin systems from lattice field theory

We construct a lattice field theory method for computing the real-time dynamics of spin systems in a thermal bath. This is done by building on previous work of Takano with Schwinger-Keldysh and functional differentiation techniques. We derive a Schwinger-Keldysh path integral for generic spin Hamiltonians, then demonstrate the method on a simple system. Our path integral has a sign problem, which generally requires exponential run time in the system size, but requires only linear storage. The latter may place this method at an advantage over exact diagonalization, which is exponential in both. Our path integral is amenable to contour deformations, a technique for reducing sign problems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Parsing Weather Variability and Wildfire Effects on the Post-Fire Changes in Daily Stream Flows: A Quantile-Based Statistical Approach and Its Application

Determining wildland fire impacts on streamflow can be problematic as the hydrology in burned watersheds is influenced by post-fire weather conditions. Here, this study presents a quantile-based analytical framework for assessing fire impacts on low and peak daily flow magnitudes, while accounting for post-fire weather influences. This framework entails (a) the bootstrap method to compute the relative change in the post-fire annual flow and weather statistics, (b) double mass analysis to detect if post-fire baseflow and quick-flow yield ratios are significantly altered, and (c) a quantile regression method to parse fire effects on flow at a specific quantile. We illustrate the applicability of this analytical framework using 44 western US streams with at least 5% of their watershed area burned. Results indicate that large, high-severity burns in upland watersheds can raise the streamflow magnitude at the 0. 05 th and 0. 95 th quantiles for at least the five post-fire years. Quantile regression results show that the median fire-related increase in flow for the five post-fire years can be up to 5,000% (Standard Error; S.E. < 2%) at the 0. 05 th E quantile and 161% (S.E. < 10%) at the 0. 95 th quantile. The fire-related increase in flow was often pronounced at the 0. 05 th quantile for streams in the Pacific Northwest and California regions. The difference in fire effects on flow (at both quantiles) across streams was related to post-fire weather, pyrology, physiography, and land cover. The proposed analytical framework can be useful for detecting and quantifying fire effects on the low and peak stream flows in burned watersheds without overlapping disturbances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A deep neural network regressor for phase constitution estimation in the high entropy alloy system Al-Co-Cr-Fe-Mn-Nb-Ni

High Entropy Alloys (HEAs) are composed of more than one principal element and constitute a major paradigm in metals research. The HEA space is vast and an exhaustive exploration is improbable. Therefore, a thorough estimation of the phases present in the HEA is of paramount importance for alloy design. Machine Learning presents a feasible and non-expensive method for predicting possible new HEAs on-the-fly. A deep neural network (DNN) model for the elemental system of: Mn, Ni, Fe, Al, Cr, Nb, and Co is developed using a dataset generated by high-throughput computational thermodynamic calculations using Thermo-Calc. The features list used for the neural network is developed based on literature and freely available databases. A feature significance analysis matches the reported HEAs phase constitution trends on elemental properties and further expands it by providing so far-overlooked features. The final regressor has a coefficient of determination ( r 2 ) greater than 0.96 for identifying the most recurrent phases and the functionality is tested by running optimization tasks that simulate those required in alloy design. The DNN developed constitutes an example of an emulator that can be used in fast, real-time materials discovery/design tasks.

36 MATERIALS SCIENCE↗