Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer systems organization”

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 649 records · Page 36

Foreword: Selected papers from the 2020 Nuclear and Emerging Technologies for Space Topical Meeting (NETS 2020)

Humanity’s curiosity of and intrigue for exploring space has been with our species since we first looked up into the night sky. This age-old curiosity has driven humankind to innovate technologies for millennia, including efforts ranging from the construction of Stonehenge to track the annual movements of the sun to the investigation of advanced propulsion systems capable of transporting people to distant planets. Today’s technical communities dedicated to researching and developing innovative technologies for space exploration have seen a recent reinvigoration. Over the past few years, impressive progress made by industry, with exceptional support from government agencies, to reestablish and develop new capabilities has renewed confidence in humanity’s ability to reach the stars—a confidence that has not been present since the days of the Apollo program. Moreover, history shows that innovation and collaboration within the scientific and engineering communities reaches unparalleled pinnacles during periods of astrophilia. Finally, given the advancements in computational capacities for advanced modeling and simulation, unprecedented breakthroughs in commercializing advanced manufacturing techniques, and powerful collaborative relationships between industry and national scientific organizations, today we are limited only by our own imaginations.

99 GENERAL AND MISCELLANEOUS↗

Interplay between London Dispersion, Hubbard U , and Metastable States for Uranium Compounds

High-throughput computational studies of lanthanide and actinide chemistry with density-functional theory are complicated by the need for Hubbard U corrections, which ensure localization of the f-electrons, but can lead to metastable states. This work presents a systematic investigation of the effects of both Hubbard U value and metastable states on the predicted structural and thermodynamic properties of four uranium compounds central to the field of nuclear fuels: UC, UN, UO 2 , and UCl 3 . We also assess the impact of the exchange-hole dipole moment (XDM) dispersion correction on the computed properties. Overall, the choice of Hubbard U value and inclusion of a dispersion correction cause larger variations in the computed geometric properties than result from metastable states. Further, the weak dependence of structure optimization on metastable states should simplify future high-throughput calculations on actinides. Conversely, addition of the dispersion correction is found to offset the repulsion introduced by the Hubbard U term and provides greatly improved agreement with experiment for both cell volumes and heats of formation. The XDM dispersion correction is largely invariant to the chosen U value, making it a robust dispersion correction for actinide systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exascale Computing and Data Handling: Challenges and Opportunities for Weather and Climate Prediction

The emergence of exascale computing and artificial intelligence offer tremendous potential to significantly advance Earth system prediction capabilities. However, enormous challenges must be overcome to adapt models and prediction systems to use these new technologies effectively. A 2022 WMO report on exascale computing recommends “urgency in dedicating efforts and attention to disruptions associated with evolving computing technologies that will be increasingly difficult to overcome, threatening continued advancements in weather and climate prediction capabilities.” Further, the explosive growth in data from observations, model and ensemble output, and postprocessing threatens to overwhelm the ability to deliver timely, accurate, and precise information needed for decision-making. Artificial intelligence (AI) offers untapped opportunities to alter how models are developed, observations are processed, and predictions are analyzed and extracted for decision-making. Given the extraordinarily high cost of computing, growing complexity of prediction systems, and increasingly unmanageable amount of data being produced and consumed, these challenges are rapidly becoming too large for any single institution or country to handle. This paper describes key technical and budgetary challenges, identifies gaps and ways to address them, and makes a number of recommendations.

Atmosphere↗

Multi-Scale Modeling Framework for Mercury Biogeochemistry

Multi-Scale modeling of mercury (Hg) geochemical speciation and reactions has been performed by integrating atomistic quantum chemical calculations with continuum scale speciation models. Major progress has been made in the improvement of quantum chemical models to calculate critical thermodynamic data for Hg complexes in aquatic environments. Rapid and reliable quantum chemical approaches have been developed for calculating acid dissociation constants (pK a ) and stability constants (log K), with calculated mean unsigned errors of 0.5 and 1.5 log units, respectively for ligand molecules and Hg complexes. At the continuum scale, systematic analysis of uncertainty propagation in mercury (Hg) speciation modeling has been conducted and was used to identify environmental conditions under which thermodynamic constant uncertainties are significant and recommended to be accounted for. The integrated framework for multi-scale modeling of mercury geochemistry is open to the research community through the web-based multiscale modeling aqueous speciation resource, AQUA-MER. The improved quantum chemical approaches for thermodynamic constant calculations are accessible through AQUA-MER and can be used to provide the missing constants in the continuum scale speciation calculations. In addition to low molecular mass Hg complex speciation, modeling natural aquatic environments also involve the transport of high molecular weight dissolved organic matter (DOM) in reactive flows simultaneously with equilibrium and kinetic reactions. To this end, atomistic MD simulations were performed to capture the details of aggregation, mechanisms and distribution of functional groups in DOM at the molecular level. The elemental composition and calculated bulk properties of the DOM models are in close agreement with experimental measurements. A travel-time based reactive transport model in the hyporheic zone of stream corridors was established for the multicomponent Hg-DOM-S system and implemented through PFLOTRAN.

54 ENVIRONMENTAL SCIENCES↗

A simulation-based parametric study of CLOU chemical looping reactor performance

Chemical looping with oxygen uncoupling (CLOU) is a variant on chemical looping combustion in which the oxygen carrier releases gaseous O 2 in the fuel reactor, making it well-suited for solid fuels, since the released gaseous O 2 readily reacts with solid char. This work presents several computational fluid dynamic (CFD) simulations of copper-based CLOU in a dual fluidized bed system, each with different operating conditions. The modeling predicted that coal particle sizes as large as 1000 μm did not significantly affect performance. Increased oxygen carrier copper loading resulted in an excess of gaseous oxygen in the product gas stream. Decreasing the oxygen carrier bed mass as well as reducing the air reactor fluidizing velocity did not supply enough oxygen to the fuel reactor to complete combustion of the coal. This generated a failure state in which the temperature continued to decrease in the fuel reactor from the lack of combustion, which in turn reduced the O 2 equilibrium partial pressure, further lowering the amount of combustion possible. Sufficient O 2 can be maintained in the fuel reactor by ensuring a high enough air reactor velocity and a large enough supply of oxygen carrier inventory to handle the chosen coal feed rate.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Combining Force Fields and Neural Networks for an Accurate Representation of Chemically Diverse Molecular Interactions

A key goal of molecular modeling is the accurate reproduction of the true quantum mechanical potential energy of arbitrary molecular ensembles with a tractable classical approximation. The challenges are that analytical expressions found in general purpose force fields struggle to faithfully represent the intermolecular quantum potential energy surface at close distances and in strong interaction regimes; that the more accurate neural network approximations do not capture crucial physics concepts, e.g., nonadditive inductive contributions and application of electric fields; and that the ultra-accurate narrowly targeted models have difficulty generalizing to the entire chemical space. We therefore designed a hybrid wide-coverage intermolecular interaction model consisting of an analytically polarizable force field combined with a short-range neural network correction for the total intermolecular interaction energy. Here, we describe the methodology and apply the model to accurately determine the properties of water, the free energy of solvation of neutral and charged molecules, and the binding free energy of ligands to proteins. The correction is subtyped for distinct chemical species to match the underlying force field, to segment and reduce the amount of quantum training data, and to increase accuracy and computational speed. For the systems considered, the hybrid ab initio parametrized Hamiltonian reproduces the two-body dimer quantum mechanics (QM) energies to within 0.03 kcal/mol and the nonadditive many-molecule contributions to within 2%. Simulations of molecular systems using this interaction model run at speeds of several nanoseconds per day.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating Combustion and Surface Chemistry Simulations

Design of modern combustion systems relies on computer models to predict how changes in design will affect performance. These models have largely displaced previous methods that rely on the designer’s intuition or costly and time-consuming physical testing. By using improved models, design cycles can be shortened, and cleaner and more efficient combustion devices can be created. This project aims to improve computer simulations of transportation fuels with the goal of making these simulations faster and more accurate for predicting combustion in vehicles.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Structurally Constrained Evolutionary Algorithm for the Discovery and Design of Metastable Phases

Metastable materials are abundant in nature and technology, showcasing remarkable properties that inspire innovative materials design. However, traditional crystal structure prediction methods, which rely solely on energetic factors to determine a structure’s fitness, are not suitable for predicting the vast number of potentially synthesizable phases that represent a local minimum corresponding to a state in thermodynamic equilibrium. Here, we present a new approach for the prediction of metastable phases with specific structural features, and interface this method with the XTALOPT evolutionary algorithm. Our method relies on structural features that include the local crystalline order (e.g., the coordination number or chemical environment), and symmetry (e.g., Bravais lattice and space group) to filter the breeding pool of an evolutionary crystal structure search. The effectiveness of this approach is benchmarked on three known metastable systems: XeN 8 , with a two-dimensional polymeric nitrogen sublattice, brookite TiO 2 , and a high pressure BaH 4 phase that was recently characterized. Additionally, a newly predicted metastable melaminate salt, P1¯WC 3 N 6 , was found to possess an energy that is lower than two phases proposed in a recent computational study. Here, the method presented here could help in identifying the structures of compounds that have already been synthesized, and developing new synthesis targets with desired properties.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Selectivity-Determining Steps in O 2 Reduction Catalyzed by Iron(tetramesitylporphyrin)

The oxygen reduction reaction (ORR) is the cathode reaction in fuel cells and its selectivity for water over hydrogen peroxide production is important for these technologies. Iron porphyrin catalysts have long been studied for the ORR, but the origins of their selectivity are not well understood because the selectivity-determining step(s) usually occur after the rate-determining step. We report here the effects of acid concentration, as well as other solution conditions such as acid pK a , on the H 2 O 2 /H 2 O selectivity in electrocatalytic ORR by iron(tetramesitylporphyrin) (Fe(TMP)) in DMF. The results show that selectivity reflects a kinetic competition in which the dependence on [HX] is one order greater for the production of H 2 O than H 2 O 2 . Based on such experimental results and computational studies, we propose that the selectivity is governed by competition between protonation of the hydroperoxo intermediate, Fe III (TMP)(OOH), to produce water versus dissociation of the HOO – ligand to yield H 2 O 2 . The data rule out a bifurcation based on the regioselectivity of protonation of the hydroperoxide, as suggested in the enzymatic systems. Furthermore, the analysis developed in this paper should be generally valuable to the study of selectivity in other multi-proton/multi-electron electrocatalytic reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Modified Embedded-Atom Method Potential for a Quaternary Fe-Cr-Si-Mo Solid Solution Alloy

Ferritic-martensitic steels, such as T91, are candidate materials for high-temperature applications, including superheaters, heat exchangers, and advanced nuclear reactors. Considering these alloys’ wide applications, an atomistic understanding of the underlying mechanisms responsible for their excellent mechano-chemical properties is crucial. Here, we developed a modified embedded-atom method (MEAM) potential for the Fe-Cr-Si-Mo quaternary alloy system—i.e., four major elements of T91—using a multi-objective optimization approach to fit thermomechanical properties reported using density functional theory (DFT) calculations and experimental measurements. Elastic constants calculated using the proposed potential for binary interactions agreed well with ab initio calculations. Furthermore, the computed thermal expansion and self-diffusion coefficients employing this potential are in good agreement with other studies. This potential will offer insightful atomistic knowledge to design alloys for use in harsh environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Maximizing long-term biohydrogen production with Clostridium thermocellum for high solids conversion of lignocellulosic biomass

Biological hydrogen production from lignocellulosic biomass sustainably couples organic waste reduction with renewable energy generation. Efficient conversion is challenged by the structural complexity of lignocellulose and resulting recalcitrance to enzymatic degradation. Clostridium thermocellum natively breaks down biomass with highly effective hemi-/cellulases systems (i.e., cellulosomes) and generates hydrogen in anaerobic cultivation, creating a compelling platform for lignocellulosic biohydrogen production. Achieving commercially viable production rates requires balancing high biomass loading and throughput against uniform mixing conditions required for enzyme dispersion, pH and temperature control, and efficient hydrogen and metabolite removal in continuous operation. To address these barriers to process intensification, we implemented novel reactor and process designs for high-solids lignocellulosic biomass fermentations using the C. thermocellum KJC19-9 strain, genetically engineered for co-utilization of cellulose and hemicellulose sugars (i.e., xylose). Via computational fluid dynamics (CFD) modeling and experimental validation, we achieved a >50% improvement in biohydrogen production with an improved anchor-type impeller morphology, coupled to a threefold reduction in agitation rate. To further reduce rheological constraints and accumulation of toxic metabolites, we then transitioned the process to sequencing fed-batch operation. The resulting process generated 24.87 L H 2 L −1 from 160 g L −1 of deacetylated and mechanically refined (DMR)-pretreated corn stover biomass over 16 days while solubilizing >95% of influent cellulose and hemicellulose, setting a new performance benchmark for continuous production of biohydrogen from lignocellulose.

08 HYDROGEN↗

Hybrid quantum-classical simulations of magic angle spinning dynamic nuclear polarization in very large spin systems

Solid-state nuclear magnetic resonance can be enhanced using unpaired electron spins with a method known as dynamic nuclear polarization (DNP). Fundamentally, DNP involves ensembles of thousands of spins, a scale that is difficult to match computationally. This scale prevents us from gaining a complete understanding of the spin dynamics and applying simulations to design sample formulations. We recently developed an ab initio model capable of calculating DNP enhancements in systems of up to ~1000 nuclei; however, this scale is insufficient to accurately simulate the dependence of DNP enhancements on radical concentration or magic angle spinning (MAS) frequency. We build on this work by using ab initio simulations to train a hybrid model that makes use of a rate matrix to treat nuclear spin diffusion. We show that this model can reproduce the MAS rate and concentration dependence of DNP enhancements and build-up time constants. Finally, we then apply it to predict the DNP enhancements in core–shell metal-organic-framework nanoparticles and reveal new insights into the composition of the particles’ shells.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DFT-based QM/MM with Particle-Mesh Ewald for Direct, Long-Range Electrostatic Embedding

In this work, we present a DFT-based, QM/MM implementation with long-range electrostatic embedding achieved by direct real-space integration of the particle mesh Ewald (PME) computed electrostatic potential. The key transformation is the interpolation of the electrostatic potential from the PME grid to the DFT quadrature grid, from which integrals are easily evaluated utilizing standard DFT machinery. We provide benchmarks of the numerical accuracy with choice of grid size and real-space corrections, and demonstrate that good convergence is achieved while introducing nominal computational overhead. Furthermore, the approach requires only small modification to existing software packages, as is demonstrated with our implementation in the OpenMM and Psi 4 software. After presenting convergence benchmarks, we evaluate the importance of long-range electrostatic embedding in three solute/solvent systems modeled with QM/MM. Water and BMIM/BF 4 ionic liquid were considered as "simple" and "complex" solvents respectively, with water and p-phenylenediamine (PPD) solute molecules treated at QM level of theory. While electrostatic embedding with standard real-space truncation may introduce negligible error for simple systems such as water solute in water solvent, errors become more significant when QM/MM is applied to complex solvents such as ionic liquids. An extreme example is the electrostatic embedding energy for oxidized PPD in BMIM/BF 4 for which real-space truncation produces severe error even at 2-3 nm cutoff distances. This latter example illustrates that utilization of QM/MM to compute redox potentials within concentrated electrolytes/ionic media requires carefully chosen long-range electrostatic embedding algorithms, with our presented algorithm providing a general and robust approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Piecewise interaction picture density matrix quantum Monte Carlo

The density matrix quantum Monte Carlo (DMQMC) set of methods stochastically samples the exact N-body density matrix for interacting electrons at finite temperature. We introduce a simple modification to the interaction picture DMQMC (IP-DMQMC) method that overcomes the limitation of only sampling one inverse temperature point at a time, instead allowing for the sampling of a temperature range within a single calculation, thereby reducing the computational cost. At the target inverse temperature, instead of ending the simulation, we incorporate a change of picture away from the interaction picture. The resulting equations of motion have piecewise functions and use the interaction picture in the first phase of a simulation, followed by the application of the Bloch equation once the target inverse temperature is reached. We find that the performance of this method is similar to or better than the DMQMC and IP-DMQMC algorithms in a variety of molecular test systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational and experimental study of different brines in temperature swing solvent extraction desalination with amine solvents

Rapid global urbanization and high-salinity wastewater disposal from industrial activities have exerted significant pressure on water resources. Over the past few years, temperature swing solvent extraction (TSSE) has been identified as a promising technique to desalinate hypersaline brines. Despite its potential, the TSSE desalination literature has been mainly based on empirical insights, and the limited molecular simulation studies have primarily focused on NaCl brines. Herein, we use molecular dynamics (MD) simulations to study the TSSE desalination of four different brines, namely, KCl, KBr, NaCl, and NaBr using diisopropylamine as the solvent. Based on both bulk and interfacial brine-diisopropylamine MD results, here we investigate the qualitative and quantitative performance of the simulations by benchmarking these results against our experimental evaluations of these same systems. MD results provide satisfactory qualitative agreement with the experimental data of water solubilization in the organic phase and amine solubilization in the aqueous phase for the KBr, KCl, and NaBr brines. Also, the molecular mechanism of solvation of ionic species by water molecules over diisopropylamine suggested by the MD simulations is in agreement with our experimental data. However, larger qualitative and quantitative deviations were observed for the NaCl brines, and this is likely due to polarization and charge transfer effects, as quantified by our quantum chemical calculations.

42 ENGINEERING↗

Combined First-Principles Calculations of Electron–Electron and Electron–Phonon Self-Energies in Condensed Systems

In this work, we present a method to efficiently combine the computation of electron–electron and electron–phonon self-energies, which enables the evaluation of electron–phonon coupling at the G 0 W 0 level of theory for systems with hundreds of atoms. In addition, our approach, which is a generalization of a method recently proposed for molecules [J. Chem. Theory Comput.2018, 14, 6269–6275], enables the inclusion of nonadiabatic and temperature effects at no additional computational cost. We present results for diamond and defects in diamond and discuss the importance of numerically accurate G 0 W 0 band structures to obtain robust predictions of zero point renormalization (ZPR) of band gaps, and of the inclusion of nonadiabatic effects to accurately compute the ZPR of defect states in the band gap.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonideal mixing effects in warm dense matter studied with first-principles computer simulations

Here, we study nonideal mixing effects in the regime of warm dense matter (WDM) by computing the shock Hugoniot curves of BN, MgO, and MgSiO 3 . First, we derive these curves from the equations of state (EOS) of the fully interacting systems, which were obtained using a combination of path integral Monte Carlo calculations at high temperature and density functional molecular dynamics simulations at lower temperatures. We then use the ideal mixing approximation at constant pressure and temperature to rederive these Hugoniot curves from the EOS tables of the individual elements. We find that the linear mixing approximation works remarkably well at temperatures above ~2 × 10 5 K, where the shock compression ratio exceeds ~3.2. The shape of the Hugoniot curve of each compound is well reproduced. Regions of increased shock compression, which emerge because of the ionization of L and K shell electrons, are well represented, and the maximum compression ratio of the Hugoniot curves is reproduced with high precision. Some deviations are seen near the onset of the L shell ionization regime, where ionization equilibrium in the fully interacting system cannot be well reproduced by the ideal mixing approximation. This approximation also breaks down at lower temperatures, where chemical bonds play an increasingly important role. However, the results imply that the equilibrium properties of binary and ternary mixtures in the regime of WDM can be derived from the EOS tables of the individual elements. This significantly simplifies the characterization of binary and ternary mixtures in the WDM and plasma phases, which otherwise requires large numbers of more computationally expensive first-principles computer simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

CheKiPEUQ Intro 2: Harnessing Uncertainties from Data Sets, Bayesian Design of Experiments in Chemical Kinetics**

When choosing experimental conditions, Bayesian statistical tools can predict the experimental choices which will yield the highest information gain. Experimental choices could be temperature, pressure, reaction time, number of measurements, reactor volume, etc.. Three example analyses are presented here, each using the software Chemical Kinetics Parameter Estimation and Uncertainty Quantification (CheKiPEUQ). Information gain is a measure of reduction of uncertainty in a model's parameters. The three chemical system examples presented each illustrate Bayesian Design of Experiments using information gain. In the first chemical example, temperature selection impacts the information gain for the free energy of reaction in a two-component equilibrium reaction. In the second example, temperature and pressure are explored for a competitive adsorption Langmuir replacement reaction system. Finally, the third example is a catalytic membrane reactor which is a culmination of the previous examples. The catalytic membrane reactor has a complex and nonlinear response in the observables which is solved by numerical evaluation. In the three examples, the experimental conditions are treated as design variables for maximizing information gain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗