Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computer model”

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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Identifying mechanistic differences between co-fed CO 2 hydrogenation and reactive CO 2 capture using Ru and Pd dual function materials

Dual function materials (DFMs) enable reactive carbon capture (RCkeC), an intensified approach to carbon dioxide capture and utilization for cost and energy input reductions. Yet, there is a fundamental lack of understanding of mechanisms around CO 2 adsorption and subsequent conversion on these materials, hindering further development. Herein, we investigated several supported alkaline metal oxides for their CO 2 adsorption characteristics to find that Na/Al 2 O 3 had the highest CO 2 adsorption capacity, accompanied by a variety of CO 2 adsorption geometries as identified by in situ DRIFTS and computational modeling. The addition of catalytic metals (Ru, Pd) increased the adsorption capacity of Na/Al 2 O 3 without altering binding modes. In the subsequent reactive desorption step, acetate and formate intermediates were observed. Notably, this mechanistic investigation identified that the formation of acetate species was unique to RCC on a DFM, as these species were not observed in co-fed hydrogenation over the DFM or RCC over a Na-free catalyst.

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

Illuminating the mechanistic impacts of an Fe-quaterpyridine functionalized crystalline poly(triazine imide) semiconductor for photocatalytic CO 2 reduction

The strategy of incorporating earth-abundant catalytic centers into light-absorbing architectures is desirable from the viewpoint of low cost, low toxicity, and versatility for activating small molecules to produce solar-based fuels. Herein, we show that an Fe-quaterpyridine molecular catalyst can be anchored to a light-absorbing, crystalline, carbon nitride (PTI), to yield a molecular-catalyst/material hybrid, Fe-qpy-PTI, capable of facilitating CO 2 reduction to CO selectively (up to ∼97–98%) in aqueous solution under low-intensity light irradiation. This hybrid material leverages the ability of the Fe-qpy catalyst to bind CO 2 upon a one-electron reduction, as achieved by transfer of excited electrons from the carbon-nitride semiconductor. At a low incident power density of only 50 mW cm −2 , the catalytic activity of the hybrid material was measured across a range of catalyst loadings from 0.1–3.8 wt%, yielding CO rates of up to 596 μmol g −1 h −1 for a 3.8 wt% loading during a 3 h experiment. Over the course of 8 h, the hybrid material attained a CO evolution rate of 608 μmol g −1 h −1 and 305 turnovers for a TOF of ∼38 h −1 and an apparent quantum yield of 2.6%. Higher light intensities provided an initial increase in activity but negatively impacted photocatalytic rates with time, with an AQY of 0.6% at 150 mW cm −2 and 0.4% at 250 mW cm −2 . Transient absorption spectroscopy results showed electron survival probabilities consistent with the trends in observed product rates. Computational modeling was also used to evaluate and understand the mechanistic pathway of the high product selectivity for CO versus H 2 . These results thus help unveil key factors for leveraging the mechanistic understanding of molecular catalysts for CO 2 reduction for pairing with light absorbing semiconductors and establishing optimal conditions to attain maximal rates in aqueous solution.

McGuigan, Scott [Baylor University, Waco, TX (Unit↗
DOE OSTI2025

Subsurface hydrogen, curvature, and strain: lessons from electro-reduction of benzaldehyde on nano-structured Pd catalysts

The unique ability of palladium (Pd) to absorb hydrogen and form a bulk hydride is vital for chemical transformations that involve hydrogenation reactions. Nano-structured Pd catalysts offer a promise of tuning these reaction rates by exploiting variations of reactant binding energies depending on the surface structure and morphological constraints that result in inhomogeneous strain. However, the interplay between the nano-structure of Pd and the ability of Pd to adsorb (and absorb) hydrogen as well as other reactive species needs to be better understood for a rational understanding of competitive chemical transformations at Pd surfaces. We consider the effects of the surface corrugation, strain, and subsurface Pd hydride on the reduction of benzaldehyde to benzyl alcohol in two qualitatively different samples – Pd nanoparticles and Pd gels formed by quasi-one-dimensional chains of these nanoparticles. Our electrochemical measurements and computational modelling suggest that surface concave sites, inherent to Pd gels, facilitate hydrogen transfer to the Pd subsurface region, thus weakening benzaldehyde binding to the surface. This effect is further modulated by the strain, depending on the local coordination environment on the corrugated surface. Furthermore, these findings demonstrate how structurally complex samples in the form of gels provide degrees of freedom for controlling the behavior of metal catalysts that are not available in isolated nanoparticles, which paves the way for new approaches in the design of catalytic materials and synthesis of metal hydrides.

Padavala, Sri Krishna Murthy [University of Minnes↗
DOE OSTI2023

Outcomes of WPEC SG47 on "Use of Shielding Integral Benchmark Archive and Database for Nuclear Data Validation"

The Working Party on International Nuclear Data Evaluation Co-operation Subgroup 47 (WPECSG47) entitled "Use of Shielding Integral Benchmark Archive and Database for Nuclear Data Validation" was organised from 2019 and 2022 with the objectives to promote more systematic and wider use of shielding benchmark experiments in nuclear data (ND) and transport code validation and development, to provide feedback on the Shielding Integral Benchmark Archive and Database (SINBAD), and to promote its further development in coordination with the Expert Group on Physics of Reactor Systems (EGPRS). Altogether 9 meetings, the large majority (8) held remotely, were organised during the past 3 years to discuss the experience on the use of SINBAD, evaluation of new benchmarks and improvements to be contributed to the database which was severely neglected and lacking maintenance over the past ← 10+ years. Several proposals for new or updated benchmark evaluation were presented and discussed, such as FNG copper, LLNL pulsed spheres, CIAE iron sphere, KFK 1977 gamma measurements, Rez Fe sphere, ASPIS, ORNL Oxygen broomstick, TIARA and others. Complementing the database with new features was also discussed, for example providing the nuclear data sensitivity profiles more systematically would facilitate and better guide the use of data. Information on the geometry, (radiation source) and materials available in CAD format is expected to allow an easier and less error prone reference for computational model preparation and a potential input to CAD based workflows. Inputs for various transport codes and other benchmark data from participants have been shared via the NEA GitLab which could hopefully in the future evolve and form a bases for critically checked and validated benchmark data. Future development of SINBAD will be monitored by EGPRS and the newly created SINBAD Task Force.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Hydropower Scheduling Toolchains: Comparing Experiences in Brazil, Norway, and USA and Implications for Synergistic Research

While hydropower scheduling is a well-defined problem, there are institutional differences that need to be identified to promote constructive and synergistic research. We study how established toolchains of computer models are organized to assist operational hydropower scheduling in Brazil, Norway, and the United States’ Colorado River System (CRS). These three systems have vast hydropower resources, with numerous, geographically widespread, and complex reservoir systems. Although the underlying objective of hydropower scheduling is essentially the same, the systems are operated in different market contexts and with different alternative uses of water, where the stakeholders’ objectives clearly differ. This in turn leads to different approaches when it comes to the scope, organization, and use of models for operational hydropower scheduling and the information flow between the models. Here, we describe these hydropower scheduling toolchains, identify the similarities and differences, and shed light on the original ideas that motivated their creation. We then discuss the need to improve and extend the current toolchains and the opportunities to synergistic research that embrace those contextual differences.

13 HYDRO ENERGY↗
DOE OSTI2020

On the expansion of non-ideal copper plasma into vacuum

The expansion of non-ideal copper plasma into vacuum is analyzed for the conditions typical for explosive electron emission in vacuum arcs. The computational model solves the Euler equations with an equation of state (EoS) for weakly non-ideal plasma, taking into account the ionization energy correction due to electron–ion Coulomb coupling. We have ascertained that the EoS has insignificant influence on the plasma expansion while the plasma properties and composition change drastically. Based on our simulation results, we discuss the validity of the “frozen” state theory often used in the vacuum arc plasma diagnostics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Roles of bandgrading, lifetime, band alignment, and carrier concentration in high-efficiency CdSeTe solar cells

CdSeTe alloying has significantly increased the efficiency of CdTe-based solar technology. Here, computational modeling compares how different CdSeTe bandgrading, carrier lifetimes, band alignment, and carrier concentrations contribute to transport, recombination, and performance. We find that the gain in photocurrent caused by bandgap narrowing alone is insufficient to describe experimental efficiency gains. Performance can be increased by adjusting CdSeTe compositions and bandgrading depths. However, these performance gains are small relative to the contributions of enhanced lifetime by Se alloying, which can explain record cell efficiency gains with minimal open-circuit voltage loss despite significant bandgap narrowing. Similarly, CdSeTe band alignment shifts can significantly increase performance if front interface recombination is prevalent. For a wide range of CdSeTe grading profiles, the hole density is a critical component to achieve efficiencies exceeding 25%.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Interfacial water asymmetry at ideal electrochemical interfaces

Controlling electrochemical reactivity requires a detailed understanding of the charging behavior and thermodynamics of the electrochemical interface. Experiments can independently probe the overall charge response of the electrochemical double layer by capacitance measurements and the thermodynamics of the inner layer with potential of maximum entropy measurements. Relating these properties by computational modeling of the electrochemical interface has so far been challenging due to the low accuracy of classical molecular dynamics (MD) for capacitance and the limited time and length scales of ab initio MD. Here, we combine large ensembles of long-time-scale classical MD simulations with charge response from electronic density functional theory to predict the potential-dependent capacitance of a family of ideal aqueous electrochemical interfaces with different peak capacitances. We show that while the potential of maximum capacitance varies, this entire family exhibits an electrode charge of maximum capacitance (CMC) between -2.9 and -2.2 μC/cm 2 , regardless of the details in the electronic response. Simulated heating of the same interfaces reveals that the entropy peaks at a charge of maximum entropy (CME) of -5.1 ± 0.6 μC/cm 2 , in agreement with experimental findings for metallic electrodes. The CME and CMC both indicate asymmetric response of interfacial water that is stronger for negatively charged electrodes, while the difference between CME and CMC illustrates the richness in behavior of even the ideal electrochemical interface.

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

A laser powder bed fusion system for operando synchrotron x-ray imaging and correlative diagnostic experiments at the Stanford Synchrotron Radiation Lightsource

Laser powder bed fusion (LPBF) is a highly dynamic multi-physics process used for the additive manufacturing (AM) of metal components. Improving process understanding and validating predictive computational models require high-fidelity diagnostics capable of capturing data in challenging environments. Synchrotron x-ray techniques play a vital role in the validation process as they are the only in situ diagnostic capable of imaging sub-surface melt pool dynamics and microstructure evolution during LPBF-AM. In this article, a laboratory scale system designed to mimic LPBF process conditions while operating at a synchrotron facility is described. The system is implemented with process accurate atmospheric conditions, including an air knife for active vapor plume removal. Significantly, the chamber also incorporates a diagnostic sensor suite that monitors emitted optical, acoustic, and electronic signals during laser processing with coincident x-ray imaging. The addition of the sensor suite enables validation of these industrially compatible single point sensors by detecting pore formation and spatter events and directly correlating the events with changes in the detected signal. Experiments in the Ti–6Al–4V alloy performed at the Stanford Synchrotron Radiation Lightsource using the system are detailed with sufficient sampling rates to probe melt pool dynamics. X-ray imaging captures melt pool dynamics at frame rates of 20 kHz with a 2 µm pixel resolution, and the coincident diagnostic sensor data are recorded at 470 kHz. Furthermore, this work shows that the current system enables the in situ detection of defects during the LPBF process and permits direct correlation of diagnostic signatures at the exact time of defect formation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Vorticity screening by dense layers

We investigate variable density layers in the presence of vortices, by analogy to line charges in the presence of a dielectric slab. We adapt solutions for dielectric layers from the literature to the variable density fluid case, obtaining exact hypergeometric expressions for flow velocities induced everywhere by a vortex in the vicinity of a finite-thickness slab of different densities. The dimensionless Atwood number, which appears naturally in other variable density phenomena, such as the Rayleigh–Taylor instability, reappears here as the natural ordering parameter for the influence of the slab, which acts primarily to screen the influence of vorticity across the slab to a lower effective circulation. The solution takes the form of scaling correction factors that might be applied to computational models or evaluated directly. Extensions, where vorticity is localized within the dense layer itself and amplified outside of it, are considered, as well as some speculative applications such as inertial confinement fusion capsule designs.

42 ENGINEERING↗
DOE OSTI2023

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Chemistry↗
DOE OSTI2023

Continuum kinetic investigation of the impact of bias potentials in the current saturation regime on sheath formation

In this work, we examine sheath formation in the presence of bias potentials in the current saturation regime for pulsed power fusion experiments. It is important to understand how the particle and heat fluxes at the wall may impact the wall material and affect electrode degradation. Here, simulations are performed using the 1X-1V Boltzmann–Poisson system for a proton–electron plasma in the presence of bias potentials ranging from 0 to 10 kV. The results indicate that the sheath near the anode remains generally the same as that of a classical sheath without the presence of a bias potential. However, the sheath near the cathode becomes more prominent with a larger potential drop, a significant decrease in the electron density, and larger sheath lengths. The spatially constant current density increases to a saturation value with the increasing bias potential. For high bias potentials, the material choice needs to consider that the anode faces significantly larger particle and heat fluxes compared to the cathode. In general, the results trend with theory with differences attributed to the simplified assumptions in the theory and the kinetic effects considered in the simulations. Due to the significant computational cost of a well resolved 1X-2V simulation, only one such simulation is performed for the 5 kV case showing higher current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

The dynamics, mixing, and thermonuclear burn of compressed foams with varied gas fills

Inertial confinement fusion (ICF) implosions involve highly coupled physics and complex hydrodynamics that are challenging to model computationally. Due to the sensitivity of such implosions to small features, detailed simulations require accurate accounting of the geometry and dimensionality of the initial conditions, including capsule defects and engineering features such as fill tubes used to insert gas into the capsule, yet this is computationally prohibitive. It is therefore difficult to evaluate whether discrepancies between the simulation and experiment arise from inadequate fidelity to the capsule geometry and drive conditions, uncertainties in physical data used by simulations, or inadequate physics. We present results from detailed high-resolution three-dimensional simulations of ICF implosions performed as part of the MARBLE campaign on the National Ignition Facility [Albright et al., Phys. Plasmas 29, 022702 (2022)]. These experiments are foam-filled separated-reactant experiments, where deuterons reside in the foam and tritons reside in the capsule gas fill and deuterium–tritium (DT) fusion reactions only occur in the presence of mixing between these materials. Material mixing in these experiments is primarily seeded by shock interaction with the complex geometry of the foam and gas fill, which induces the Richtmyer–Meshkov instability. We compare results for experiments with two different gas fills (ArT and HT), which lead to significant differences in the hydrodynamic and thermodynamic developments of the materials in the implosion. Our simulation results show generally good agreement with experiments and demonstrate a substantial impact of hydrodynamic flows on measured ion temperatures. The results suggest that viscosity, which was not included in our simulations, is the most important unmodeled physics and qualitatively explains the few discrepancies between the simulation and experiment. The results also suggest that the hydrodynamic treatment of shocks is inadequate to predict the heating and yield produced during shock flash, when the shock converges at the center of the implosion. Alternatively, underestimation of the level of radiative preheat from the shock front could explain many of the differences between the experiment and simulation. Nevertheless, simulations are able to reproduce many experimental observables within the level of experimental reproducibility, including most yields, time-resolved X-ray self-emission images, and an increase in burn-weighted ion temperature and neutron down-scattered ratio in the line of sight that includes a jet seeded by the glue spot that joins capsule hemispheres.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Formation of a spherical plasma liner for plasma-jet-driven magneto-inertial fusion

Plasma-jet-driven magneto-inertial fusion is an alternative approach to controlled nuclear fusion, which aims to utilize a line-replaceable dense plasma liner as a repetitive spherical compression driver. In this experiment, first measurements of the formation of a spherical argon plasma liner formed from 36 discrete pulsed plasma jets are obtained on the Plasma Liner Experiment. Properties including liner uniformity and morphology, plasma density, temperature, and ram pressure are assessed as a function of time throughout the implosion process and indicate an apparent transition from initial kinetic inter-jet interpenetration to collisional regime near stagnation times, in accordance with theoretical expectation. A lack of primary shock structures between adjacent jets during flight implies that arbitrarily smooth liners may be formed by way of corresponding improvements in jet parameters and control. The measurements facilitate the benchmarking of computational models and understanding the scaling of plasma liners toward fusion-relevant energy density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An updated technique to obtain explosive kinetics data on microsecond timescales

There are few techniques available for chemists to obtain time-to-explosion data with known temperature inputs at the early stages of the design and synthesis of new explosives. In the 1960s, a technique was developed to rapidly heat milligram-quantities of confined explosives to ~1000 K on microsecond timescales. Wenograd loaded explosives inside stainless steel hypodermic needles, connected them to a fireset and rapidly discharged a capacitor through the steel. He obtained the temperature by measuring the needle resistance in a Wheatstone bridge arrangement and the time to explosion from a needle rupture. However, owing to the narrow-gauge needles used in the original research, the experiment was only possible with melt-castable explosives; it was never replicated, and modern diagnostics are now available with advances beyond the 1960s. Here, we report the development of the High Explosives Initiation Time (HEIT) test, which utilizes a 250 J pulsed power system to heat the needles. This work extends the Wenograd approach by using optical diagnostics, computational modeling, and advanced techniques to measure needle resistance and needle rupture. Preliminary rate information for pentaerythritol tetranitrate (PETN) will be presented.

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

Effect of particle shape on stratification in drying films of binary colloidal mixtures

The role of particle shape in evaporation-induced auto-stratification in polydisperse colloidal suspensions is explored with molecular dynamics simulations of mixtures of spheres and aspherical particles. A unified framework based on the competition between diffusion and diffusiophoresis is proposed to understand the effects of shape and size dispersity. In general, particles diffusing more slowly (e.g., larger particles) tend to accumulate more strongly at the evaporation front. However, larger particles have larger surface areas and therefore greater diffusiophoretic mobility. Hence, they are more likely to be driven away from the evaporation front via diffusiophoresis. For a rapidly dried bidisperse suspension containing small and large spheres, the competition leads to “small-on-top” stratification. Here, we employ a computational model in which the diffusion coefficient is inversely proportional to particle mass. For a mixture of spheres and aspherical particles with similar mass, the diffusion contrast is reduced, and the spheres are always enriched at the evaporation front as they have the smallest surface area for a given mass and, therefore, the lowest diffusiophoretic mobility. Furthermore, for a mixture of solid and hollow spheres that have the same outer radius and thus the same surface area, the diffusiophoretic contrast is suppressed, and the system is dominated by diffusion. Consequently, the solid spheres, which have a larger mass and diffuse more slowly, accumulate on top of the hollow spheres. Finally, for a mixture of thin disks and long rods that differ significantly in shape but have similar mass and surface area, both diffusion and diffusiophoresis contrasts are suppressed, and the mixture does not stratify.

Classical molecular dynamic simulations↗
DOE OSTI2020

Phase-shift speckle-shearing interferometry

We have numerically simulated the process of measuring stress – strain states by the method of speckle-shearing interferometry using the phase-shift technique. A computer model with the possibility of setting its strain and roughness is developed, which includes a model of a diffusely reflecting test object corresponding to the characteristics of a real membrane, as well as a speckle interferometer model that allows speckle interferograms to be obtained for different speckle sizes and angles between interfering beams. The process of reconstructing the object surface topogram from model speckle interferograms by the phase-shift technique is implemented. Using the developed models, a two-dimensional shearogram are obtained, which is a derivative of the strain field of a circular membrane. Comparison of the results of numerical simulation with experimental data shows that the differences (rms deviations) do not exceed 0.02 μm. It is also shown that the error of interferogram reconstruction by the phase-shift technique increases significantly when the test object strains exceed 12 μm. (paper)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Emergence of self-organized multivortex states in flocks of active rollers

Active matter, both synthetic and biological, demonstrates complex spatiotemporal self-organization and the emergence of collective behavior. A coherent rotational motion, the vortex phase, is of great interest because of its ability to orchestrate well-organized motion of self-propelled particles over large distances. However, its generation without geometrical confinement has been a challenge. Here, we show by experiments and computational modeling that concentrated magnetic rollers self-organize into multivortex states in an unconfined environment. We find that the neighboring vortices more likely occur with the opposite sense of rotation. Our studies provide insights into the mechanism for the emergence of coherent collective motion on the macroscale from the coupling between microscale rotation and translation of individual active elements. These results may stimulate design strategies for self-assembled dynamic materials and microrobotics.

36 MATERIALS SCIENCE↗