Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational studies”

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

Optimization of Energy Flow through Synthetic Metabolic Modules and Regulatory Networks in a Model Photosynthetic Eukaryotic Microbe

Photosynthetic organisms have recently gained considerable attention for a role in development of renewable energy sources. Genome-enabled systems biology methods, coupled with functional and synthetic genomics, present opportunities to develop sustainable and economical applications such as fuel production within the next 10 to 15 years. However, optimization of light-driven metabolism for biomass or biofuel production will require a detailed systems biology understanding of photosynthetic processes and cellular metabolism. Genome-scale metabolic models (GEMs) are at the core of systems analysis of cellular processes and form a common organizational framework for analyses of data resulting from functional genomics experimental work and computational studies. Therefore, there is a clear demand for high quality photosynthetic model organisms and the appropriate computational tools that enable systems analysis of light-driven metabolism. Through research conducted we expanded the currently available repertoire of photosynthetic GEMs to include the commercially valuable model diatom Phaeoctylum tricornutum. Diatoms have a peculiar and distinct evolutionary footprint and represent a major eukaryotic lineage that is taxonomically and functionally distinct from green and red algae and vascular plants. Therefore, the true potential for light-driven metabolism aimed at biofuel production remains poorly understood at a systems level for a large subset of the global diversity of photosynthetic organisms. The metabolic capabilities of P. tricornutum were comparatively modeled with those from other photosynthetic groups in order to elucidate the occurrence of metabolic traits within and between phototrophs. Additionally, this research resulted in significant extension of the COnstraints Based Reconstruction and Analysis (COBRA) Toolbox to accommodate the crucial need for infrastructure required for ‘omics data integration and analysis in the context of genome-scale models. Therefore, the proposed research achieved two important goals. First, within the broad scope of photosynthetic organisms, we functionally compared and, as a result, identified cellular processes that require optimization in order to enable deployment as biofuel feedstock. Second, the proposed research resulted in development of key computational infrastructure, which can be further extended to other biological systems, that is currently lacking but necessary for multiple ‘omics data integration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Mechanical and Thermal Properties of Low-Density Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y Alloys

Refractory high-entropy alloys (RHEAs) Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y ((x, y) = (0, 0), (0, 10), and (10, 15)) were computationally studied to obtain a low density and a better mechanical property. The density functional theory (DFT) method was employed to compute the structural and mechanical properties of the alloys, based on a large unit cell model of randomly distributed elements. Debye–Grüneisen theory was used to study the thermal properties of Al 20+x Cr 20-x Mo 20-y Ti 20 V 20+y . The phase diagram calculation shows that all three RHEAs have a single body-centered cubic (BCC) structure at high temperatures ranging from 1000 K to 2000 K. The RHEA Al 30 Cr 10 Mo 5 Ti 20 V 35 has shown a low density of 5.16 g/cm 3 and a hardness of 5.56 GPa. The studied RHEAs could be potential candidates for high-temperature application materials where high hardness, ductility, and low density are required.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Actinide Molten Salts: A Machine-Learning Potential Molecular Dynamics Study

We know that actinide molten salts represent a class of important materials in nuclear energy. Understanding them at a molecular level is critical to proper and optimal design of relevant technological applications. Yet, owing to the complexity of electronic structure due to the 5f orbitals, computational studies of heavy elements in condensed phases using ab initio potentials to study the structure and dynamics of these elements embedded in molten salts are difficult. This lack of efficient computational protocols makes it difficult to obtain information on properties that require extensive statistical sampling like transport. To tackle this problem, we adopted a machine-learning approach to study ThCl 4 -NaCl and UCl 3 -NaCl binary systems. The machine-learning potential, with the density functional theory accuracy, allows us to obtain long molecular dynamics trajectories (ns) for large systems (10 3 atoms) at a considerably low computing cost, thereby efficiently gaining information about their bonding structures, thermodynamics, and dynamics at a range of temperature. We observed a considerable change in the coordination environments of actinide elements and their characteristic coordination-sphere lifetime. Our study also suggests that actinides in molten salts may not follow well known entropy-scaling laws.

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

Parallel computing for power system climate resiliency: Solving a large-scale stochastic capacity expansion problem with mpi-sppy

Here we propose a nodal stochastic generation and transmission expansion planning model that incorporates the output from high-resolution global climate models through load and generation availability scenarios. We implement our model in Pyomo and perform computational studies on a realistically-sized test case of the California electric grid in a high performance computing environment. We propose model reformulations and algorithm tuning to efficiently solve this large problem using a variant of the Progressive Hedging Algorithm. We utilize the parallelization capabilities and overall versatility of mpi-sppy, exploiting its hub-and-spoke architecture to concurrently obtain inner and outer bounds on an optimal expansion plan. Initial results show that instances with 360 representative days on a system with over 8,000 buses can be solved to within 5% of optimality in under 4 h of wall clock time, a first step towards solving a large-scale power system expansion planning problem across a wide range of climate-informed operational scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Revealing the Influence of Salts on the Hydration Structure of Ionic SDS Micelles by Contrast-Variation Small-Angle Neutron Scattering

The influence of lithium chloride (LiCl) on the hydration structure of anionic micelles of sodium dodecyl sulfate (SDS) in water was studied using the contrast-variation small-angle neutron scattering (SANS) technique. In the past, extensive computational studies have shown that the distribution of invasive water plays a critical role in the self-organization of SDS molecules and the stability of the assemblies. However, in past scattering studies the degree of the hydration level was not examined explicitly. Here, a series of contrast-variation SANS data was analyzed to extract the intramicellar radial distributions of invasive water and SDS molecules from the evolving spectral lineshapes caused by the varying isotopic ratios of water. By addressing the intramicellar inhomogeneous distributions of water and SDS molecules, a detailed description of how the counterion association influences the micellization behavior of SDS molecules is provided. The extension of our method can be used to provide an in-depth insight into the micellization phenomenon, which is commonly found in many soft matter systems.

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

Understanding the performance of a bisphosphonate Ru water oxidation catalyst

Water oxidation catalysts (WOCs) create electrons that allow generating H 2 from water and sunlight, a promising process for the production of clean energy. Here, the mechanism of water oxidation mediated by Ru(2,2'-bipyridine-6,6'-diphosphonato)(4-picoline) 2 has been studied computationally to comprehend the results obtained in the experiments performed by the Concepcion and Grotjahn groups. Our study was performed at pH = 8 and 1. At pH = 8, the phosphonates are fully deprotonated and the catalysis, which is shown to be more energetically costly than that of the dicarboxylato Ru catalyst counterpart, takes place through a mechanism that involves a bimolecular interaction between two metal-oxo units (I2M). At pH = 1, only one of the phosphonates of the catalyst can deprotonate. After testing all possible pathways and comparing with experimental data, we conclude that the catalysis proceeds neither through a water nucleophilic attack nor via I2M, but rather we hypothesize that the rate-determining step involves both the catalyst and the ceric ammonium nitrate.

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

Solid face sheets enable lattice metamaterials to withstand high-amplitude impulsive loading without yielding

Owing to their ability to provide tunable mechanical responses, lattice materials are frequently studied to elucidate their response to static and dynamic loads. However, these roles are typically in opposition: static loads must be supported sufficiently far away from the onset of buckling or yielding, whereas dynamic loads are typically ameliorated by crushing of the lattice, which provides excellent energy-absorption due to the large plastic deformation accompanying densification. In contrast, this work considers the octet truss as an exemplar topology, in a structural role where it must simultaneously support static loads while enduring high-amplitude impulsive loads. This study focuses on the ability to withstand impulsive loads without yielding, an essential prerequisite to enduring dual loading. Computational studies using the ALE3D hydrocode were performed to examine the response of the octet truss under a short temporal width impulse shape associated with laser-driven shocks. A key finding was that covering the lattice with a solid face sheet and treating this face sheet thickness as a design variable allows the Taylor-like pulse to be attenuated prior to entering the weaker lattice, at the cost of added mass up front. Experimental validation was accomplished by laser-driven shock testing, using octet trusses printed out of Ti-5Al-5V-5Mo-3Cr. The results show that for a given quantity of mass, the attenuation is maximized when as much mass as possible is moved into the face sheet, leaving a more slender lattice structure. The effect of placing mass in the face sheet rather than lattice beams dominates the effect of relative density, to the point where a low-mass structure with most of the mass concentrated in the face sheet can outperform a high-mass structure with most of the mass in the lattice. Finally, by further understanding the propagation of short pulse width waves within under-dense structures, this study expand the domain of applicability of such structures, including lattice materials, to challenging dual-loading regimes spanning decades of strain rates.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Computational Analysis of the Energetic Stability of High-Entropy Structures of a Prototypical Lanthanide-Based Metal–Organic Framework

High-entropy materials are characterized by their complex compositions, typically comprising five or more elements in near-equiatomic proportions. Applying this concept to metal ions in metal−organic frameworks (MOFs) has paved the way for exploring a new class of high-entropy MOFs. While the compositional strategy of high-entropy materials leverages configurational entropy to aid thermodynamic stability, it also poses significant analytical challenges due to the vast compositional landscape and diverse phases that these materials can adopt. We present a computational study of several complexities associated with selecting potential high-entropy versions of a prototype lanthanidebased MOF. We compute the energetics of metal mixing of these heterometallic MOFs using density functional theory (DFT) and machine learning interatomic potential (MLIP) methods. The use of MLIP methods allows a systematic exploration of the convex hull of thermodynamically stable MOF structures containing up to 5 distinct metals.

Chemical structure↗
DOE OSTI2025

Homologous Ladder Cyclohexasilanes

Here, we report the synthesis of five new examples of ladder cyclohexasilanes, possessing up to three consecutive fused rings and differing in relative ring fusion configuration and side chain structure. By coupling a 1,4-dipotassiooligosilyl dianion to a cyclohexasilane, we obtained bi- and tricyclic ladder cyclohexasilanes. Combined experimental and theoretical studies suggested that annulation could favor the cis configuration under kinetic control, while the trans configuration predominates under thermodynamic control. Computational studies show that with each additional ring in the trans-diastereomeric series, the predicted onset of light absorption shifts to longer wavelengths.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Computers in neutron science

Abstract Neutron scattering provides a sensitive probe for materials structure and dynamics and has had a profound impact on research ranging from quantum topological materials to archaeology to aerospace alloys. These experiments are incredibly data-rich and often challenging to interpret. Computational studies have been integral in the growing success of neutron science. From data collection and reduction to analysis and interpretation, computers are used at every stage of modern neutron experiments. State-of-the-art software allows for a deeper understanding of these powerful unique experiments. This themed issue, guest-edited by Dr. Jeff Armstrong, Dr. Keith T Butler, and Dr. Matthew R Ryder, brings together a collection of research highlighting the state-of-the-art in the application of computers in neutron science.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Global Sensitivity Analysis Using the Ultra‐Low Resolution Energy Exascale Earth System Model

Abstract For decades, Arctic temperatures have increased twice as fast as average global temperatures. As a first step toward quantifying parametric uncertainty in Arctic climate, we performed a variance‐based global sensitivity analysis (GSA) using a fully coupled, ultra‐low resolution (ULR) configuration of version 1 of the U.S. Department of Energy's Energy Exascale Earth System Model (E3SMv1). Specifically, we quantified the sensitivity of six quantities of interests (QOIs), which characterize changes in Arctic climate over a 75 year period, to uncertainties in nine model parameters spanning the sea ice, atmosphere, and ocean components of E3SMv1. Sensitivity indices for each QOI were computed with a Gaussian process emulator using 139 random realizations of the random parameters and fixed preindustrial forcing. Uncertainties in the atmospheric parameters in the Cloud Layers Unified by Binormals (CLUBB) scheme were found to have the most impact on sea ice status and the larger Arctic climate. Our results demonstrate the importance of conducting sensitivity analyses with fully coupled climate models. The ULR configuration makes such studies computationally feasible today due to its low computational cost. When advances in computational power and modeling algorithms enable the tractable use of higher‐resolution models, our results will provide a baseline that can quantify the impact of model resolution on the accuracy of sensitivity indices. Moreover, the confidence intervals provided by our study, which we used to quantify the impact of the number of model evaluations on the accuracy of sensitivity estimates, have the potential to inform the computational resources needed for future sensitivity studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Th 2 O – , Th 2 Au – , and Th 2 AuO 1,2 – Anions: Photoelectron Spectroscopic and Computational Characterization of Energetics and Bonding

The observation and characterization of the anions: Th 2 O – , Th 2 Au – , and Th 2 AuO 1,2 – is reported. These species were studied through a synergetic combination of anion photoelectron spectroscopy and ab initio correlated molecular orbital theory calculations at the CCSD(T) level with large correlation-consistent basis sets. To better understand the energetics and bonding in these anions and their corresponding neutrals, a range of smaller diatomic to tetratomic species were studied computationally. Correlated molecular orbital theory calculations at the CCSD(T) level showed that in most of these cases, there are close-lying anions and neutral clusters with different geometries and spin states and are consistent with the experimentally observed spectra. Thus, comparison of experimentally determined and computationally predicted vertical detachment energies and electron affinities for different optimized geometries and spin states shows excellent agreement to within 0.1 eV. The structures for both the neutrals and anions have a significant ionic component to the bonding because of the large electron affinity of the Au atom and modest ionization potentials for Th 2 , Th 2 O, and Th 2 O 2 . The analysis of the bonding for the Th–Th bonds from the molecular orbitals is consistent with this ionic model. The results show that there is a wide variation in the bond distance from 2.7 to 3.5 Å for the Th–Th bonds all of which are less than twice the atomic radius of Th of 3.6 Å. The bond distances encompass bond orders from 4 to 0. Furthermore, there can be different bond orders for the same bond distance depending on the nature of the ionic bonding suggesting that one may not be able to correlate the bond order with the bond distance in these types of clusters. In addition, the presence of an Au atom may provide a unique probe of the bonding in such clusters because of its ability to accept an electron from clusters with modest ionization potentials.

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

A Combined DFT, Energy Decomposition, and Data Analysis Approach to Investigate the Relationship Between Noncovalent Interactions and Selectivity in a Flexible DABCOnium/Chiral Anion Catalyst System

Developing strategies to study reactivity and selectivity in flexible catalyst systems has become an important topic of research. Herein, we report a combined experimental and computational study aimed at understanding the mechanistic role of an achiral DABCOnium cofactor in a regio- and enantiodivergent bromocyclization reaction. In this study, it was found that electron-deficient aryl substituents enable rigidified transition states via an anion–π interaction with the catalyst, which drives the selectivity of the reaction. In contrast, electron-rich aryl groups on the DABCOnium result in significantly more flexible transition states, where interactions between the catalyst and substrate are more important. An analysis of not only the lowest-energy transition state structures but also an ensemble of low-energy transition state conformers via energy decomposition analysis and machine learning was crucial to revealing the dominant noncovalent interactions responsible for observed changes in selectivity in this flexible system.

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

Design of Ru-Ni diatomic sites for efficient alkaline hydrogen oxidation

Anion exchange membrane fuel cells are limited by the slow kinetics of alkaline hydrogen oxidation reaction (HOR). Here, we establish HOR catalytic activities of single-atom and diatomic sites as a function of *H and *OH binding energies to screen the optimal active sites for the HOR. As a result, the Ru-Ni diatomic one is identified as the best active center. Guided by the theoretical finding, we subsequently synthesize a catalyst with Ru-Ni diatomic sites supported on N-doped porous carbon, which exhibits excellent catalytic activity, CO tolerance, and stability for alkaline HOR and is also superior to single-site counterparts. In situ scanning electrochemical microscopy study validates the HOR activity resulting from the Ru-Ni diatomic sites. Furthermore, in situ x-ray absorption spectroscopy and computational studies unveil a synergistic interaction between Ru and Ni to promote the molecular H 2 dissociation and strengthen OH adsorption at the diatomic sites, and thus enhance the kinetics of HOR.

25 ENERGY STORAGE↗
DOE OSTI2024

Stable Pentagonal Layered Palladium Diselenide Enables Rapid Electrosynthesis of Hydrogen Peroxide

Electrosynthesis of hydrogen peroxide (H 2 O 2 ) via the two-electron oxygen reduction reaction (2e – ORR) is promising for various practical applications, such as wastewater treatment. However, few electrocatalysts are active and selective for 2e – ORR yet are also resistant to catalyst leaching under realistic operating conditions. Here, a joint experimental and computational study reveals active and stable 2e – ORR catalysis in neutral media over layered PdSe 2 with a unique pentagonal puckered ring structure type. Computations predict active and selective 2e – ORR on the basal plane and edge of PdSe 2 , but with distinct kinetic behaviors. Further, electrochemical measurements of hydrothermally synthesized PdSe 2 nanoplates show a higher 2e – ORR activity than other Pd–Se compounds (Pd 4 Se and Pd 17 Se 15 ). PdSe 2 on a gas diffusion electrode can rapidly accumulate H 2 O 2 in buffered neutral solution under a high current density. The electrochemical stability of PdSe 2 is further confirmed by long device operational stability, elemental analysis of the catalyst and electrolyte, and synchrotron X-ray absorption spectroscopy. This work establishes a new efficient and stable 2e – ORR catalyst at practical current densities and opens catalyst designs utilizing the unique layered pentagonal structure motif.

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

Effect of Ni deposition on CrI 3 monolayer: First-principles study

In our recent experimental work, we demonstrated that Ni thin film deposition can be used to tune the surface magnetic states in bulk CrI 3 . Here, we present the results of a computational study of Ni deposition on a CrI 3 monolayer. In particular, we compare structural, electronic, and magnetic properties of three different systems: a pristine CrI 3 monolayer, a CrI 3 monolayer with partial Ni coverage, and a CrI 3 monolayer with full Ni coverage. The corresponding cells used in our calculations are Cr 2 I 6 , Cr 2 I 6 Ni 1 , and Cr 2 I 6 Ni 3 . In Cr 2 I 6 Ni 1 , the Ni atom is positioned over the hollow site, which corresponds to the lowest energy configuration. In Cr 2 I 6 Ni 3 , the three Ni atoms are positioned over the single hollow site, and over the two Cr sites. Our calculations indicate that the partial coverage of CrI 3 monolayer with Ni retains the semiconducting nature of the parent compound, but results in a strong reduction of the energy band gap. At the same time, the full coverage of CrI 3 monolayer with Ni results in a metallic transition, with a highly spin-polarized (possibly half-metallic) electronic structure. The magnetic alignment of the CrI 3 monolayer with the deposited Ni is ferrimagnetic, due to the anti-aligned magnetic moments of Cr and Ni atoms. While the partial coverage of the CrI 3 monolayer with Ni has a moderate impact on the lattice parameter of the parent compound, full coverage by Ni results in a significant increase of the lattice constant. For both partial and full coverages, the lowest energy configuration corresponds to the Ni atoms being approximately aligned with the out-of-plane coordinate of the surface iodine layer. In addition, in the full coverage mode, the Cr layer exhibits a significant shift away from the Ni atoms, effectively “detaching” the Ni / I atomic layer from the rest of the cell.

2D magnetism↗
DOE OSTI2022

Predictive coarse-grained (CG) modeling of morphologies in polymer nanocomposites with specific and directional intermolecular interactions (Final Report)

The overarching goal of the proposed work was to develop predictive models for investigating structure and dynamics in soft materials with chemistries that have specific and directional molecular interactions. The motivation behind studying materials with specific and directional interactions lies in the many desirable features these interactions provide when designing novel soft materials. Soft materials with specific and direction interactions (such as hydrogen bonds or H-bonds) can have a) thermally reversible phase behavior with different functions with varying temperature, b) precisely tuned nanostructure with desirable geometries that afford unique physical properties (e.g., color response, mechanical properties) and c) well-mixed/blended morphologies that are useful for variety of applications in energy field (e.g., materials for batteries require use of blended polymers where one domain gives superior mechanical properties and one domain promotes electrical conduction). Notably, biology makes extensive use of specific and directional interactions, in many cases based on H-bonds, to construct materials with precisely defined architectures and properties. Engineered soft materials with precisely tuned nanostructures and improved processiblity through thermoresponsive phase behavior are useful in numerous applications that are relevant to the Department of Energy (DOE) including high efficiency electronic devices, light-weight high-strength composite materials for batteries and fuel cells, and polymer membranes for separations, etc. While past computational studies have been tremendously useful in understanding molecular phenomena and guiding synthesis of new macromolecular soft materials for a wide variety of applications, the inability to capture small scale specific and directional interactions alongside macromolecular length and time scales represented a key limitation of most studies to date. Our work in this project addressed this grand challenge in computational materials chemistry, i.e., ability to model the anisotropic, directional, and specific interactions that govern the behavior of many macromolecular soft matter systems of interest, thus, has greatly expanded the predictive potential of simulations. Specifically the key outcomes were: successful development of new coarse-grained (CG) polymer models to study generic and specific polymer chemistries in which hydrogen-bonding interactions are dominant. These CG models were then used in molecular simulations to study structure and thermodynamics in polymer nanocomposites and blends; some studies were conducted in collaboration with experimentalists. We also published a perspective and a viewpoint which included some of the work we completed in this DOE project; we believe these perspective and viewpoint articles guide other researchers in the soft materials community on how to extend the computational approaches and models we have developed for the purposes of their studies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Magnetization dynamics fingerprints of an excitonic condensate t 2 g 4 magnet

The competition between spin-orbit coupling (SOC) λ and electron-electron interaction U leads to a plethora of novel states of matter, extensively studied in the context of t 2 g 4 and t 2 g 5 materials, such as ruthenates and iridates. Excitonic magnets—the antiferromagnetic state of bounded electron-hole pairs-are prominent examples of phenomena driven by those competing energy scales. Interestingly, recent theoretical studies predicted that excitonic magnets can be found in the ground state of SOC t 2 g 4 Hubbard models. Here we present a detailed computational study of the magnetic excitations in that excitonic magnet, employing one-dimensional chains (via density matrix renormalization group) and small two-dimensional clusters (via Lanczos). Specifically, first we show that the low-energy spectrum is dominated by a dispersive (acoustic) magnonic mode, with extra features arising from the λ = 0 state in the phase diagram. Second, and more importantly, we found a novel magnetic excitation forming a high-energy optical mode with the highest intensity at wave-vector q → 0 . In the excitonic condensation regime at large U , we also have found a novel high-energy π mode composed solely of orbital excitations. Overall, these features do not appear all together in any of the neighboring states in the phase diagram and thus constitute unique fingerprints of the t 2 g 4 excitonic magnet, of importance in the analysis of neutron and resonant inelastic x-ray scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗