Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer experiments”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Thermal Experiments for Fractured Rock Characterization: Theoretical Analysis and Inverse Modeling

Abstract Field‐scale properties of fractured rocks play a crucial role in many subsurface applications, yet methodologies for identification of the statistical parameters of a discrete fracture network (DFN) are scarce. We present an inversion technique to infer two such parameters, fracture density and fractal dimension, from cross‐borehole thermal experiments data. It is based on a particle‐based heat‐transfer model, whose evaluation is accelerated with a deep neural network (DNN) surrogate that is integrated into a grid search. The DNN is trained on a small number of the heat‐transfer model runs and predicts the cumulative density function of the thermal field. The latter is used to compute fine posterior distributions of the (to be estimated) parameters. Our synthetic experiments reveal that fracture density is well constrained by data, while fractal dimension is harder to determine. Adding nonuniform prior information related to the DFN connectivity improves the inference of this parameter.

Zhou, Zitong↗
DOE OSTI2022

Project Description for Student Symposium 2022: Modeling Approach to Critical in the Upcoming CERBERUS Experiment

The Critical Experiment Reflected By copper to BettEr Understand Scattering [CERBERUS] is an integral experiment that will be executed at the end of FY23. The experiment is designed to study the intermediate energy range, in particular elastic and inelastic scattering of neutrons in copper. Approaching criticality safely is of utmost importance, and having a computer mode is a useful tool when determining changes that can safely be made to the experiment. Multiple simulations are needed to analyze the approach to critical, and manually creating MCNP input files is both tedious and error prone. To solve this problem, a python script was created to develop and run these input files.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Boson sampling with Gaussian input states: Toward efficient scaling and certification

A universal quantum computer of large scale is not available yet, however, intermediate models of quantum computation would still permit demonstrations of a quantum computational advantage over classical computing and could challenge the Extended Church-Turing Thesis. One of these models based on single photons interacting via linear optics is called Boson Sampling. Although Boson Sampling was demonstrated and the threshold to claim quantum computational advantage was achieved, the question of how to scale up Boson Sampling experiments remains. To make progress with this problem, here we present a practically achievable pathway to scale Boson Sampling experiments by combining continuous-variable quantum information and temporal encoding. Here, we propose the combination of switchable dual-homodyne and single-photon detections, the temporal loop technique, and scattershot-based Boson Sampling. We detail the required assumptions for concluding computational hardness for this configuration. Furthermore, this particular combination of techniques moves towards an efficient scaling and certification of Boson Sampling, all in a single experimental setup.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Chapter 14: Dynamic Flux Analysis: An Experimental Approach of Fluxomics

Metabolic flux analysis represents an essential perspective to understand cellular physiology and offers quantitative information to guide pathway engineering. A valuable approach for experimental elucidation of metabolic flux is dynamic flux analysis, which estimates the relative or absolute flow rates through a series of metabolic intermediates in a given pathway. It is based on kinetic isotope labeling experiments, liquid chromatography-mass spectrometry (LC-MS), and computational analysis that relate kinetic isotope trajectories of metabolites to pathway activity. Herein, we illustrate the mathematic principles underlying the dynamic flux analysis and mainly focus on describing the experimental procedures for data generation. This protocol is exemplified using cyanobacterial metabolism as an example, for which reliable labeling data for central carbon metabolites can be acquired quantitatively. This protocol is applicable to other microbial systems as well and can be readily adapted to address different metabolic processes.

BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Inside the hologram: reconstructing the bulk observer’s experience

We develop a holographic framework for describing the experience of bulk observers in AdS/CFT, that allows us to compute the proper time and energy distribution measured along any bulk worldline. Our method is formulated directly in the CFT language and is universal: it does not require knowledge of the bulk geometry as an input. When used to propagate operators along the worldline of an observer falling into an eternal black hole, our proposal resolves a conceptual puzzle raised by Marolf and Wall. Notably, the prescription does not rely on an external dynamical Hamiltonian or the AdS boundary conditions and is, therefore, outlining a general framework for the emergence of time.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Jet-images — deep learning edition

Building on the notion of a particle physics detector as a camera and the collimated streams of high energy particles, or jets, it measures as an image, we investigate the potential of machine learning techniques based on deep learning architectures to identify highly boosted W bosons. Modern deep learning algorithms trained on jet images can out-perform standard physically-motivated feature driven approaches to jet tagging. We develop techniques for visualizing how these features are learned by the network and what additional information is used to improve performance. Finally, this interplay between physically-motivated feature driven tools and supervised learning algorithms is general and can be used to significantly increase the sensitivity to discover new particles and new forces, and gain a deeper understanding of the physics within jets.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Anticipating gelation and vitrification with medium amplitude parallel superposition (MAPS) rheology and artificial neural networks

Abstract Anticipating qualitative changes in the rheological response of complex fluids (e.g., a gelation or vitrification transition) is an important capability for processing operations that utilize such materials in real-world environments. One class of complex fluids that exhibits distinct rheological states are soft glassy materials such as colloidal gels and clay dispersions, which can be well characterized by the soft glassy rheology (SGR) model. We first solve the model equations for the time-dependent, weakly nonlinear response of the SGR model. With this analytical solution, we show that the weak nonlinearities measured via medium amplitude parallel superposition (MAPS) rheology can be used to anticipate the rheological aging transitions in the linear response of soft glassy materials. This is a rheological version of a technique called structural health monitoring used widely in civil and aerospace engineering. We design and train artificial neural networks (ANNs) that are capable of quickly inferring the parameters of the SGR model from the results of sequential MAPS experiments. The combination of these data-rich experiments and machine learning tools to provide a surrogate for computationally expensive viscoelastic constitutive equations allows for rapid experimental characterization of the rheological state of soft glassy materials. We apply this technique to an aging dispersion of Laponite ® clay particles approaching the gel point and demonstrate that a trained ANN can provide real-time detection of transitions in the nonlinear response well in advance of incipient changes in the linear viscoelastic response of the system.

Lennon, Kyle R. (ORCID:0000000212515461)↗
DOE OSTI2017

Analysis of thermal history effects on mechanical anisotropy of 3D-printed polymer matrix composites via in situ X-ray tomography

The tensile behavior of an additively manufactured (AM) polymer matrix composite (PMC) is studied with in situ X-ray computed microtomography (CT) and digital volume correlation (DVC). In this experiment, the effects of recycled material content and print direction on the selective laser-sintered (SLS) material's mechanical response are explored. The PMC samples are printed in a tensile specimen geometry with gage lengths parallel to all three orthogonal, primary sintering directions. In situ tensile-CT experiments are conducted at Argonne National Laboratory's Advanced Photon Source 2-BM beamline. Analysis of the AM PMC's tensile response, failure, and strain evolution is analyzed both from a conventional standpoint, using the load-displacement data recorded by the loading fixture, and from a microstructural standpoint by applying DVC analysis to the reconstructed volumes. Significant variations on both strength and ductility are observed from both vantages with respect to print direction and the recycled material content in the printed parts. As a result, it is found that the addition of recycled source material with a thermal history reduces the tensile strength of the SLS composite for all directions, but the effect is drastic on the strength in the layering direction.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermomechanical conversion in metals: dislocation plasticity model evaluation of the Taylor-Quinney coefficient

Using a partitioned-energy thermodynamic framework which assigns energy to that of atomic configurational stored energy of cold work and kinetic-vibrational, in this study we derive an important constraint on the Taylor-Quinney coefficient, which quantifies the fraction of plastic work that is converted into heat during plastic deformation. Associated with the two energy contributions are two separate temperatures – the ordinary temperature for the thermal energy and the effective temperature for the configurational energy. We show that the Taylor-Quinney coefficient is a function of the thermodynamically defined effective temperature that measures the atomic configurational disorder in the material. Finite-element analysis of recently published experiments on the aluminum alloy 6016-T4 [1], using the thermodynamic dislocation theory (TDT), shows good agreement between theory and experiment for both stress-strain behavior and temporal evolution of the temperature. The simulations include both conductive and convective thermal energy loss during the experiments, and significant thermal gradients exist within the simulation results. Computed values of the differential Taylor-Quinney coefficient are also presented and suggest a value which differs between materials and increases with increasing strain.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Assessment of critical flaw sizes and crack driving forces during additive manufacturing of metallic materials

Additive manufacturing (AM) of complex engineering components is often plagued by a high susceptibility to cracking, particularly in high-strength metallic materials. While alloy design efforts have made progress in mitigating solidification defects, there remains a need for mechanistic guidelines to predict susceptibility to solid-state cracking. To address this gap, driving forces for the growth of melt pool cracks are calculated across a wide range of alloys using an efficient computational framework. Calculations are coupled with rapid single track laser experiments to elucidate trends in cracking from laser melting. The analyses conducted here highlight the important role of material properties in susceptibility to cracking, notably fracture toughness and elastic modulus. An important finding is that residual stresses that are limited in magnitude to the yield stress of the material are likely insufficient to drive cracking during cooling. Furthermore, the implications of these results are discussed in the context of alloy design for AM and residual stress accumulation during AM.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Frequency-independent approach to calculate physical optics radiations with the quadratic concave phase variations

In this work, we develop the numerical steepest descent path (NSDP) method to calculate the physical optics (PO) radiations with the quadratic concave phase variations. With the surface integral equation method, the physical optics (PO) scattered fields are formulated and further reduced to the surface integrals. The high frequency physical critical points contributions, including the stationary phase points, the boundary resonance points and the vertex points are comprehensively studied via the proposed NSDP method. The key contributions of this work are twofold. One is that together with the PO integrals taking the quadratic parabolic and hyperbolic phase terms, this work makes the NSDP theory be complete for treating the PO integrals with quadratic phase variations. Another is that, in order to illustrate the transition effect of the high frequency physical critical points, in this work, we consider and further extend the NSDP method to calculate the PO integrals with the coalescence of the high frequency critical points. Numerical results for the highly oscillatory PO integral with the coalescence of the critical points are given to verify the efficiency of the proposed NSDP method. The NSDP method could achieve the frequency independent computational workload and error controllable accuracy in all the numerical experiments, especially for the case of the coalescence of the high frequency critical points.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

DRIPS: A framework for dimension reduction and interpolation in parameter space

Reduced-order models are often used to describe the behavior of complex systems, whose simulation with a full model is too expensive, or to extract salient features from the full model’s output. We introduce a new model-reduction framework DRIPS (dimension reduction and interpolation in parameter space) that combines the offline local model reduction with the online parameter interpolation of reduced-order bases (ROBs). The offline step of this framework relies on dynamic mode decomposition (DMD) to build a low-rank linear surrogate model, equipped with a local ROB, for quantities of interest derived from the training data generated by repeatedly solving the (nonlinear) high-fidelity model for multiple parameter points. The online step consists of the construction of a parametric reduced-order model for each target/test point in the parameter space, with the interpolation of ROBs done on a Grassman manifold and the interpolation of reduced-order operators done on a matrix manifold. The DMD component enables DRIPS to model (typically low-dimensional) quantities of interest directly, without having to access the (typically high-dimensional and possibly nonlinear) operators in a high-fidelity model that governs the dynamics of the underlying high-dimensional state variables, as required in projection-based reduced-order modeling. A series of numerical experiments suggests that DRIPS yields a model reduction, which is computationally more efficient than the commonly used projection-based proper orthogonal decomposition; it does so without requiring a prior knowledge of the governing equation for quantities of interest. Furthermore, for the nonlinear systems considered, DRIPS is more accurate than Gaussian-process interpolation (Kriging).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Experimental Insights into the Coupling of Methane Combustion and Steam Reforming in a Catalytic Plate Reactor in Transient Mode

The microstructured reactor concept is very promising technology to develop a compact reformer for distributed hydrogen generation. In this work, a catalytic plate reactor (CPR) is developed and investigated for the coupling of methane combustion (MC) and methane steam reforming (MSR) over Pt/Al 2 O 3 -coated microchannels in cocurrent and counter-current modes in transient experiments during start-up. A three-dimensional (3D) computational fluid dynamics (CFD) simulation shows uniform velocity and pressure distribution profiles in microchannels. For a channel velocity from 5.1 to 57.3 m/s in the combustor, the oxidation of methane is complete and self-sustainable without explosion, blow-off, or extinction; nevertheless, flashbacks are observed in counter-current mode. In the reformer, the maximum methane conversion is 84.9% in cocurrent mode, slightly higher than that of 80.2% in counter-current mode at a residence time of 33 ms, but at the cost of three times higher energy input in the combustor operating at ~1000 °C. Nitric oxide (NO) is not identified in combustion products, but nitrous oxide (N 2 O) is a function of coupling mode and forms significantly in cocurrent mode. Furthermore, this research would be helpful to establish the start-up strategy and environmental impact of compact reformers on a small scale.

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

Simulating X-ray Absorption Spectra with Linear-Response Density Cumulant Theory

We present a new approach for simulating X-ray absorption spectra based on linear-response density cumulant theory (LR-DCT) [Copan, A. V.; Sokolov, A. Yu. J. Chem. Theory Comput. 2018, 14, 4097–4108]. Our new method combines the LR-ODC-12 formulation of LR-DCT with core- valence separation approximation (CVS) that allows us to efficiently access high-energy core-excited states. We describe our computer implementation of the CVS-approximated LRODC-12 method (CVS-ODC-12) and benchmark its performance by comparing simulated X-ray absorption spectra to those obtained from experiment for several small molecules. Our results demonstrate that the CVS-ODC-12 method shows good agreement with experiment for relative spacings between transitions and their intensities, but the excitation energies are systematically overestimated. When compared to results from excited-state coupled cluster methods with single and double excitations, the CVS-ODC-12 method shows a similar performance for intensities and peak separations, while coupled cluster spectra are less shifted, relative to experiment. An important advantage of CVS-ODC-12 is that its excitation energies are computed by diagonalizing a Hermitian matrix, which enables efficient computation of transition intensities.

Peng, Ruojing↗
DOE OSTI2018

Structure, Polarization, and Sum Frequency Generation Spectrum of Interfacial Water on Anatase TiO 2

The photocatalytic activity of TiO 2 for water splitting has been known for decades, yet the adsorption structure and hydrogen bonding of water at the interface with TiO 2 have remained controversial. We investigate the prototypical aqueous interface with anatase TiO 2 (101) using ab initio molecular dynamics (AIMD) with the strongly constrained and appropriately normed (SCAN) density functional, recently shown to provide an excellent description of the properties of bulk liquid water. We find that water forms a stable bilayer of intact molecules with ice-like dynamics and enhanced dipole moment and polarizability on the anatase surface. Here, the orientational order and H-bond environment of interfacial water are reflected in the computed sum frequency generation (SFG) spectrum, which agrees well with recent measurements in the OH stretching frequency range (3000–3600 cm –1 ). Additional AIMD simulations for a model interface with 66% of dissociated water in the contact layer show that surface hydroxyls disrupt the order in the bilayer and lead to a much faster orientational dynamics of interfacial water. Nonetheless, the computed SFG spectrum for the hydroxylated surface also agrees with experiment, suggesting that SFG measurements in a wider frequency range would be necessary to unambiguously identify the character of interfacial water on anatase.

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

An analysis and evaluation of the WeFold collaborative for protein structure prediction and its pipelines in CASP11 and CASP12

Every two years groups worldwide participate in the Critical Assessment of Protein Structure Prediction (CASP) experiment to blindly test the strengths and weaknesses of their computational methods. CASP has significantly advanced the field but many hurdles still remain, which may require new ideas and collaborations. In 2012 a web-based effort called WeFold, was initiated to promote collaboration within the CASP community and attract researchers from other fields to contribute new ideas to CASP. Members of the WeFold coopetition (cooperation and competition) participated in CASP as individual teams, but also shared components of their methods to create hybrid pipelines and actively contributed to this effort. We assert that the scale and diversity of integrative prediction pipelines could not have been achieved by any individual lab or even by any collaboration among a few partners. The models contributed by the participating groups and generated by the pipelines are publicly available at the WeFold website providing a wealth of data that remains to be tapped. Here, we analyze the results of the 2014 and 2016 pipelines showing improvements according to the CASP assessment as well as areas that require further adjustments and research.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

Using Wannier functions to improve solid band gap predictions in density functional theory

Enforcing a straight-line condition of the total energy upon removal/addition of fractional electrons on eigen states has been successfully applied to atoms and molecules for calculating ionization potentials and electron affinities, but fails for solids due to the extended nature of the eigen orbitals. Here we have extended the straight-line condition to the removal/addition of fractional electrons on Wannier functions constructed within the occupied/unoccupied subspaces. It removes the self-interaction energies of those Wannier functions, and yields accurate band gaps for solids compared to experiments. It does not have any adjustable parameters and the computational cost is at the DFT level. This method can also work for molecules, providing eigen energies in good agreement with experimental ionization potentials and electron affinities. Our approach can be viewed as an alternative approach of the standard LDA+U procedure.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Two mechanisms of nanoparticle generation in picosecond laser ablation in liquids: the origin of the bimodal size distribution

The synthesis of chemically clean and environmentally friendly nanoparticles through pulsed laser ablation in liquids has shown a number of advantages over conventional chemical synthesis methods and has evolved into a thriving research field attracting laboratory and industrial applications. The fundamental understanding of processes leading to the nanoparticle generation, however, still remains elusive. In particular, the origin of bimodal nanoparticle size distributions in femto- and picosecond laser ablation in liquids, where small nanoparticles (several nanometers) with narrow size distribution are commonly observed to coexist with larger (tens to hundreds of nanometers) ones, has not been explained so far. In this paper, joint computational and experimental efforts are applied to understand the mechanisms of nanoparticle formation in picosecond laser ablation in liquids and to explain the bimodal nanoparticle size distributions. The results of a large-scale atomistic simulation reveal the critical role of the dynamic interaction between the ablation plume and the liquid environment, leading to the generation of large nanoparticles through a sequence of hydrodynamic instabilities at the plume-liquid interface and a concurrent nucleation and growth of small nanoparticles in an expanding metal-liquid mixing region. The computational predictions are supported by a series of stroboscopic videography experiments showing the emergence of small satellite bubbles surrounding the main cavitation bubble generated in single pulse experiments. Lastly, carefully timed double pulse irradiation triggers expansion of secondary cavitation bubbles indicating, in accord with the simulation results, the presence of localized sites of laser energy deposition (possibly large nanoparticles) injected into the liquid at the early stage of the bubble formation.

36 MATERIALS SCIENCE↗