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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Solvent-Mediated Control of Nanocellulose Dispersion: An Integrated Computational and Experimental Investigation

Fibrillated cellulose derived from forestry feedstocks represents a renewable and high-strength materials platform for circular bioeconomies. However, its practical implementation is hindered by the irreversible aggregation of nanocellulose architectures, including cellulose nanofibers (CNFs). Solvent-based dispersion offers a simple and practical route to prevent CNF aggregation. Here, in this work, we integrate classical and enhanced sampling molecular dynamics (MD) simulations with experimental suspension rheology and atomic force microscopy (AFM) to elucidate how solvent environments tune CNF–CNF interactions and dispersion stability. CNF–CNF contact free energies computed from MD simulations reveal reduced aggregation in acetone/water, γ-valerolactone (GVL)/water, and tetrahydrofuran (THF)/water and pure acetone compared with pure water, reflecting stronger CNF-solvent relative to inter-CNF interactions. Correspondingly, CNF-solvent suspensions in these solvent systems exhibit stronger inter-fibril network structures and enhanced recovery compared to water, indicating improved CNF-solvent affinity. Liquid cell AFM imaging in acetone–water mixtures and in pure acetone further confirm the presence of well-dispersed CNFs. By combining multiscale computation with targeted experiments, this study establishes a rational framework for solvent design to achieve stable nanocellulose dispersions for high-strength biobased materials and efficient bioenergy conversion.

cellulose↗
DOE OSTI2017

Progress in Machine Learning Studies for the CMS Computing Infrastructure

Here, computing systems for LHC experiments developed together with Grids worldwide. While a complete description of the original Grid-based infrastructure and services for LHC experiments and its recent evolutions can be found elsewhere, it is worth to mention here the scale of the computing resources needed to fulfill the needs of LHC experiments in Run-1 and Run-2 so far.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

JARVIS-Leaderboard: a large scale benchmark of materials design methods

Abstract Lack of rigorous reproducibility and validation are significant hurdles for scientific development across many fields. Materials science, in particular, encompasses a variety of experimental and theoretical approaches that require careful benchmarking. Leaderboard efforts have been developed previously to mitigate these issues. However, a comprehensive comparison and benchmarking on an integrated platform with multiple data modalities with perfect and defect materials data is still lacking. This work introduces JARVIS-Leaderboard, an open-source and community-driven platform that facilitates benchmarking and enhances reproducibility. The platform allows users to set up benchmarks with custom tasks and enables contributions in the form of dataset, code, and meta-data submissions. We cover the following materials design categories: Artificial Intelligence (AI), Electronic Structure (ES), Force-fields (FF), Quantum Computation (QC), and Experiments (EXP). For AI, we cover several types of input data, including atomic structures, atomistic images, spectra, and text. For ES, we consider multiple ES approaches, software packages, pseudopotentials, materials, and properties, comparing results to experiment. For FF, we compare multiple approaches for material property predictions. For QC, we benchmark Hamiltonian simulations using various quantum algorithms and circuits. Finally, for experiments, we use the inter-laboratory approach to establish benchmarks. There are 1281 contributions to 274 benchmarks using 152 methods with more than 8 million data points, and the leaderboard is continuously expanding. The JARVIS-Leaderboard is available at the website: https://pages.nist.gov/jarvis_leaderboard/

36 MATERIALS SCIENCE↗
DOE OSTI2023

PETSc/TAO Users Manual (Rev. 3.19)

This manual describes the use of the Portable, Extensible Toolkit for Scientific Computation (PETSc) and the Toolkit for Advanced Optimization (TAO) for the numerical solution of partial differential equations and related problems on high-performance computers. PETSc/TAO is a suite of data structures and routines that provide the building blocks for the implementation of large-scale application codes on parallel (and serial) computers. PETSc uses the MPI standard for all distributed memory communication. PETSc/TAO includes a large suite of parallel linear solvers, nonlinear solvers, time integrators, and opti mization that may be used in application codes written in Fortran, C, C++, and Python (via petsc4py; see Getting Started). PETSc provides many of the mechanisms needed within parallel application codes, such as parallel matrix and vector assembly routines. The library is organized hierarchically, enabling users to employ the level of abstraction that is most appropriate for a particular problem. By using techniques of object-oriented programming, PETSc provides enormous flexibility for users. PETSc is a sophisticated set of software tools; as such, for some users it initially has a much steeper learning curve than packages such as MATLAB or a simple subroutine library. In particular, for individuals without some computer science background, experience programming in C, C++, python, or Fortran and experience using a debugger such as gdb or lldb, it may require a significant amount of time to take full advantage of the features that enable efficient software use. However, the power of the PETSc design and the algorithms it incorporates may make the efficient implementation of many application codes simpler than “rolling them” yourself. For many tasks a package such as MATLAB is often the best tool; PETSc is not intended for the classes of problems for which effective MATLAB code can be written. There are several packages, built on PETSc, that may satisfy your needs without requiring directly using PETSc. We recommend reviewing these packages functionality before starting to code directly with PETSc. PETSc can be used to provide a “MPI parallel linear solver” in an otherwise sequential, or OpenMP parallel code. This approach cannot provide extremely large improvements in the application time by utilizing large numbers of MPI processes but can still improve the performance. Certainly all parts of a previously sequential code need not be parallelized but the matrix generation portion must be parallelized to expect true scalability to large numbers of MPI processes. See PCMPI for details on how to utilize the PETSc MPI linear solver server. Since PETSc is under continued development, small changes in usage and calling sequences of routines will occur. PETSc has been supported for twenty-five years; see mailing list information on our website for information on contacting support.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Linear permeability evolution of expanding conduits due to feedback between flow and fast phase change: Linear Permeability Evolution of Conduit

Conduits are ubiquitous and critical pathways for many fluids relevant for geophysical processes such as magma, water, and gases. Predicting flow through conduits is challenging when the conduit geometry coevolves with the flow. We theoretically show that the permeability (k) of a conduit whose walls are eroding due to fast phase change increases linearly with time because of a self-reinforcing mechanism. This simple result is surprising given complex feedbacks between flow, transport, and phase change. The theory is congruent with previous experimental observations of fracture dissolution in calcite. Supporting computational fracture dissolution experiments showed that k only slightly increases until the dissolution front reaches the narrowest conduit constriction, after which the linear evolution of k manifests. Finally, the theory holds across multiple scales and a broad range of Peclet and Damkohler numbers and thus advances the prediction of dynamic mass fluxes through expanding conduits in various geologic and environmental settings

58 GEOSCIENCES↗
DOE OSTI2020

Toward Higher Voltage Solid-State Batteries by Metastability and Kinetic Stability Design

The energy density of battery systems is limited largely by the electrochemical window of the electrolyte. Here, the combined thermodynamic and kinetic effects of mechanically induced metastability are shown to greatly widen the operational voltage window of solid-state batteries based on ceramic-sulfide electrolytes. Solid electrolyte voltage stability up to 10 V is achieved with minimal degradation, far beyond the capability of organic liquid electrolytes. Furthermore, combined experiment, ab initio computation, and theoretical modeling identify the nature of mechanically constrained Li 10 GeP 2 S 12 decomposition both within the bulk and at interfaces with cathode materials at very high voltages. Previously unclear kinetic processes are identified that, when properly implemented, can potentially allow solid-state full cells with remarkably high operational voltages.

25 ENERGY STORAGE↗
DOE OSTI2021

Practical choices for space-filling designs

Space-filling designs are now commonly used as a flexible model-free strategy for providing good coverage throughout an input space of interest for a variety of computer and physical experiment design scenarios. Some of the preliminary choices about how to frame the problem and which type of design to use can have a substantial impact on the success or failure of the experiment, and yet how to make these critical choices is often under-emphasized in the literature. In this paper, we explore several of the practical choices required by the experimenter and describe a sequence of steps to help create an ideal design that matches the goals and constraints of the experiment. These choices include the specification of the input space, the scaling of the variables, the degree of uniformity of the design points across the input space, and the space-filling characteristics. In addition, some new tools for defining the weights to implement non-uniform space-filling designs are provided. Finally, the methods are demonstrated with several illustrative examples and a real-world chemical engineering experiment for carbon capture.

42 ENGINEERING↗
DOE OSTI2018

Porting DMRG++ Scientific Application to OpenPOWER

With the rapidly changing microprocessor designs and architectural diversity (multi-cores, many-cores, accelerators) for the next generation HPC systems, scientific applications must adapt to the hardware, to exploit the different types of parallelism and resources available in the architecture. To get the benefit of all the in-node hardware threads, it is important to use a single programming model to map and coordinate the available work to the different heterogeneous execution units in the node (e.g., multi-core hardware threads (latency optimized), accelerators (bandwidth optimized), etc.).Our goal is to show that we can manage the node complexity of these systems by using OpenMP for in-node parallelization by exploiting different “programming styles” supported by OpenMP 4.5 to program CPU cores and accelerators. Finding out the suitable programming-style (e.g., SPMD style, multi-level tasks, accelerator programming, nested parallelism, or a combination of these) using the latest features of OpenMP to maximize performance and achieve performance portability across heterogeneous and homogeneous systems is still an open research problem.We developed a mini-application, Kronecker Product (KP), from the original DMRG++ application (sparse matrix algebra) computational motif to experiment with different OpenMP programming styles on an OpenPOWER architecture and present their results in this paper.

Chatterjee, Ronnie↗
DOE OSTI2021

Modeling in High Temperature Corrosion: A Review and Outlook

Realizing higher operating temperatures to increase efficiency of future applications for energy conversion and storage while minimizing cost is a challenge for development of high-temperature materials. Simultaneous optimization of mechanical properties and corrosion resistance continues to be a difficult task but is essential due to the need to significantly accelerate the transition between technology readiness levels in the future. Furthermore, oxidation-induced degradation will be a critical life-limiting mechanism at increased operating temperatures. Suitable high-temperature materials cannot be solely identified by time-consuming experiments and reliable computational methods incorporating the relevant physics of processes must be considered to complement the experimental efforts. In the present work, a review of the methods employed to model oxidation-induced material degradation described in literature will be discussed. Furthermore, their capability to predict lifetime and aid in material selection will be evaluated.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mechanistic understanding of support effect on the activity and selectivity of indium oxide catalysts for CO 2 hydrogenation

Herein we present a mechanistic study on the support effect (ZrO 2 and CeO 2 ) of In 2 O 3 catalysts in CO 2 hydrogenation by a combined experimental and computational approach. Kinetic experiments and surface characterization suggested that the activity of In 2 O 3 catalysts cannot be simply correlated with the abundance of surface oxygen vacancies (O v ) formed by either H 2 -reduction or thermal treatment, which has been frequently invoked in previous studies. The support effect should originate from the electronic interactions between In 2 O 3 and the support oxide, rather than geometric factors or the difference in the particle size of In 2 O 3 . Theoretical modelling revealed that surface O v facilitate the formation and stabilization of the formate (HCOO*) intermediate. While a carbonate-like structure is favored for CO 2 adsorption on CeO 2 -supported or unsupported In 2 O 3 catalysts, CO 2 tends to bind strongly in a bent configuration on the O v site at the In 2 O 3 -ZrO 2 interface. The distinct CO 2 adsorption structures on different supported In 2 O 3 catalysts may account for the different reaction energy profiles in the subsequent hydrogenation reactions, especially the rate-limiting step, i.e., hydrogenation of HCOO* to CH 2 O* and methoxy (CH 3 O*). The relatively higher methanol selectivity of In 2 O 3 catalyst supported on ZrO 2 with respect to that on CeO 2 are suggested to stem from the greater energy difference (Δ$E_a$) between the parallel hydrogenation and C-O bond cleavage of HCOO*, which leads to the formation of methanol and CO, respectively. We report this study underlines the important role of metal-oxide-interface in determining the catalytic behavior of oxide-supported In 2 O 3 catalysts in CO 2 conversion.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

End-gas autoignition fraction and flame propagation rate in laser-ignited primary reference fuel mixtures at elevated temperature and pressure

Knock in spark-ignited (SI) engines is initiated by autoignition of the unburned gasses upstream of spark-ignited, propagating, turbulent premixed flames. Knock propensity of fuel/air mixtures is typically quantified using research octane number (RON), motor octane number (MON), or methane number (MN; for gaseous fuels), which are measured using single-cylinder, variable compression ratio engines. In this study, knock propensity of SI fuels was quantified via observations of end-gas autoignition (EGAI) in unburned gasses upstream of laser-ignited, premixed flames at elevated pressures and temperatures in a rapid compression machine. Stoichiometric primary reference fuel (PRF; n-heptane/isooctane) blends of varying reactivity (50 ≤ PRF ≤ 100) were ignited using an Nd:YAG laser over a range of temperatures and pressures, all in excess of 545 K and 16.1 bar. Laser ignition produced outwardly-propagating premixed flames. High-speed pressure measurements and schlieren images indicated the presence of EGAI. The fraction of the total heat release attributed to EGAI (i.e., EGAI fraction) varied with fuel reactivity (i.e., octane number) and the time-integrated temperature of the end-gas prior to ignition. Flame propagation rates, which were measured using schlieren images, were only weakly correlated with octane number but were affected by turbulence caused by variation in piston timing. Under conditions of low turbulence, measured flame propagation rates approached one-dimensional premixed laminar flame speed computations performed at the same conditions. Experiments were simulated with a three-dimensional CONVERGE™ model using reduced chemical kinetics (121 species, 538 reactions). The simulations accurately captured the measured flame propagation rates, as well as the variation in EGAI fraction with fuel reactivity and time-integrated end-gas temperature. The simulations also revealed low-temperature heat release as well as formaldehyde and hydrogen peroxide formation in the end-gas upstream of the propagating flame, which increased the temperature and degree of chain branching in the end-gas, ultimately leading to EGAI.

42 ENGINEERING↗
DOE OSTI2023

A simple, validated approach for design of two-dimensional periodic particle patterns via acoustophoresis

Two-dimensional patterning of microparticles enables a wide range of functional materials, including patterned energy storage electrodes, flexible electronics, and sensor arrays. Particle patterning via acoustics offers an attractive path to generate a wide variety of 2D periodic patterns that introduce tailorable hierarchical porosity, useful for controlling surface area, transport distances, and other properties. This method is most effective with micron scale particles and patterns of tens to hundreds of microns. To enable systematic exploration of the broad design space for such patterns, this work develops a model of 2D and 3D assembly of particles at high loadings and validates the obtained patterns against both experiments and more computationally intensive modeling techniques. Using this simple model, connections are mapped between input parameters (like actuation conditions, particle volume fraction, material properties) and output geometrical features (like void size and shape, pattern connectivity, and surface area) so that they can be tailored to given applications. The utility of this simple model is illustrated by predicting and then experimentally demonstrating new hierarchical patterns resulting from multiple waves of different frequencies interacting. These multiscale patterns offer the potential to lift the limits on surface area, diffusion distances, and other features.

42 ENGINEERING↗
DOE OSTI2021

Data on erosion and hydrogen fuel retention in Beryllium plasma-facing materials

ITER will use beryllium as a plasma-facing material in the main chamber, covering a total surface area of about 620 m 2 . Given the importance of beryllium erosion and co-deposition for tritium retention in ITER, significant efforts have been made to understand the behaviour of beryllium under fusion-relevant conditions with high particle and heat loads. This paper provides a comprehensive report on the state of knowledge of beryllium behaviour under fusion-relevant conditions: the erosion mechanisms and their consequences, beryllium migration in JET, fuel retention and dust generation. The paper reviews basic laboratory studies, advanced computer simulations and experience from laboratory plasma experiments in linear simulators of plasma–wall interactions and in controlled fusion devices using beryllium plasma-facing components. A critical assessment of analytical methods and simulation codes used in beryllium studies is given. The overall objective is to review the existing set of data with a broad literature survey and to identify gaps and research needs to broaden the database for ITER.

Erosion–deposition↗
DOE OSTI2021

CUORE opens the door to tonne-scale cryogenics experiments

The past few decades have seen major developments in the design and operation of cryogenic particle detectors. This technology offers an extremely good energy resolution – comparable to semiconductor detectors – and a wide choice of target materials, making low temperature calorimetric detectors ideal for a variety of particle physics applications. Rare event searches have continued to require ever greater exposures, which has driven them to ever larger cryogenic detectors, with the CUORE experiment being the first to reach a tonne-scale, mK-cooled, experimental mass. CUORE, designed to search for neutrinoless double beta decay, has been operational since 2017 at a temperature of about 10 mK. This result has been attained by the use of an unprecedentedly large cryogenic infrastructure called the CUORE cryostat: conceived, designed and commissioned for this purpose. In this article the main characteristics and features of the cryogenic facility developed for the CUORE experiment are highlighted. In this work, a brief introduction of the evolution of the field and of the past cryogenic facilities are given. The motivation behind the design and development of the CUORE cryogenic facility is detailed as are the steps taken toward realization, commissioning, and operation of the CUORE cryostat. The major challenges overcome by the collaboration and the solutions implemented throughout the building of the cryogenic facility will be discussed along with the potential improvements for future facilities. The success of CUORE has opened the door to a new generation of large-scale cryogenic facilities in numerous fields of science. Broader implications of the incredible feat achieved by the CUORE collaboration on the future cryogenic facilities in various fields ranging from neutrino and dark matter experiments to quantum computing will be examined.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A solvent-targeted recovery and precipitation scheme for the recycling of up to ten polymers from post-industrial mixed plastic waste

Solvent-Targeted Recovery and Precipitation (STRAP TM ) separates polymers within a plastic waste stream by selective dissolution. In this work, the STRAP framework, which combines computational modeling and experiments, was applied to develop a series of steps to separate up to 10 polymers from post-industrial mixed plastic waste (MPW) and the main components recovered were LDPE, HDPE, and PET. The STRAP steps were initially demonstrated with a physical polymer mixture containing LDPE, HDPE, PS, PVC, EVOH, PET, PP, PA6, PA66, and PA66/6, in which recoveries of 89% or higher were achieved for each polymer. Furthermore, this paper demonstrates a solvent selection approach that can be applied to separate unknown plastics materials into purer components.

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

Extension of Solvent-Assisted Linker Exchange to Supported Metal–Organic Framework Thin Films

Solvent-assisted linker exchange (SALE) is a common postsynthetic approach used to accommodate metal–organic frameworks (MOFs) for targeted applications. As this technique is relatively untested in immobilized MOF structures, we repeated the generalized SALE approach with ZIF-8 deposited directly onto solid substrates and assessed the results with a variety of techniques, including time-of-flight secondary-ion mass spectrometry (ToF-SIMS). Our methods confirm that the supported films maintain their key structural properties while allowing for linker substitution with foreign species. Alongside experiments, we used computational modeling to quantify the relative energetics of SALE reactions with a variety of linkers. These techniques together reveal that SALE is translatable to immobilized MOFs. Additionally, we recognize how various aspects of the postsynthetic approach can bear consequences for the material’s final micro- and macroscopic properties. Furthermore, the results uncover basic yet critical details needed to inform future efforts to engineer MOF films with functional properties.

Ligands↗
DOE OSTI2019

Influence of Longer Range Transfer of Vapor Interface Modified Caging Constraints on the Spatially Heterogeneous Dynamics of Glass-Forming Liquids

Based on the elastically collective nonlinear Langevin equation (ECNLE) theory of bulk relaxation in glass-forming liquids, and our recent ideas of how interface-nucleated modification of caging constraints are spatially transferred into the interior of a thick film, we present a force-based theory for dynamical gradients in thick films with one vapor interface that includes collective elasticity effects. Quantitative applications to the foundational hard sphere fluid and polymer melts of diverse fragilities are presented. In this work, we predict a roughly exponential spatial variation of the total activation barrier which is nonperturbatively modified to similar to 10 particle diameters into the film. This leads, to leading order, to the prediction of a reduced alpha relaxation time gradient of a double-exponential form characterized by a nearly constant spatial decay length, in qualitative accord with simulations. The relaxation acceleration at the surface grows exponentially with increasing packing fraction or decreasing temperature. The rate of increase with cooling of the alpha time strongly weakens upon approaching the vapor interface, with the top two layers remaining liquid-like down to similar to ~80-85% of the bulk glass transition temperature. These spatially heterogeneous changes of the temperature variation of the alpha time result in a large and long-range gradient of the local glass transition temperature. An average interfacial layer thickness relevant to ensemble-averaged dielectric and other experiments is also computed. It is predicted to be rather large at the bulk glass transition temperature, decreasing roughly linearly with heating. Remarkably, to leading order the ratio of this layer-averaged interfacial relaxation time to its bulk analogue is invariant to chemistry, volume fraction, and temperature and of modest absolute magnitude. Spatially inhomogeneous power law decoupling of the alpha relaxation time from its bulk value is predicted, with an effective exponent that decays to zero with distance from the free surface in a nearly exponential manner, trends which are in qualitative accord with recent simulations. Finally, this behavior and the double-exponential alpha time gradient are related and can be viewed as consequences of an effective quasi-universal factorization of the total barrier in films into the product of its bulk-temperature-dependent value times a function solely of location in the film.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Effect of Lithium Doping on MgO Hydroxylation and Carbonation

Recovery of magnesium from brines can potentially be used to source MgO (periclase) as a CO 2 sorbent or for Mg-based cements. However, it is not clear how common impurities in brines, such as lithium, affect the resulting MgO reactivity. Here, to test the effect of lithium incorporation on MgO reactivity for hydration and carbonation, we combined computational simulations with experiments. Experimentally altered (Mg,Li)O with a low dopant concentration (0.012 ± 0.002% w/w Li) was characterized using synchrotron-based X-ray scattering and high-resolution electron microscopy to measure reaction layer formation on (Mg,Li)O. Single-crystal X-ray diffraction analysis of (Mg,Li)O demonstrates that the incorporation of lithium leads to the formation of oxygen vacancies. The presence of vacancies is likely causing faster hydroxylation rates as predicted by ab initio molecular dynamics simulations. However, the faster hydroxylation rates likely lead to faster passivation of the surface because we observe thinner reaction layers on (Mg,Li)O samples both over short time periods (30 days) and over long time periods (28 years). After 28 years, the reaction layer on the (Mg,Li)O sample was less than one-third of the thickness of that of the pure MgO sample. In addition, over 30 days, reaction layers on (Mg,Li)O samples primarily formed at steps rather than on terraces, in contrast to our previous observations on MgO. Based on our results, naturally occurring impurities in MgO modify its reactivity even at very low concentrations and need to be considered for accurate reaction rate prediction for application of MgO as a CO 2 sorbent or in cements.

MgO↗