Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Study”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Computational Study of Triphosphine-Ligated Cu(I) Catalysts for Hydrogenation of CO 2 to Formate

The catalyzed hydrogenation of CO 2 to formate via a triphosphine-ligated Cu(I) was studied computationally. Two bases, DBU and TBD, were studied in the context of two proposed mechanisms in MeCN solvent. Of the four functionals benchmarked, M06 was generally in the best agreement with the experimentally estimated values. Activation of H 2 was explored by using LCu(DBU) + to form LCuH. Dissociation of a ligand arm results in higher barriers to form the key hydride complex, LCuH. There is no significant difference between the choice of base, DBU or TBD with respect to the proposed mechanisms. We propose that the experimentally observed differences between DBU and TBD reactivity for this mechanism are due to off-pathway changes.

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

Photoinduced charge transfer in Zn(II) and Au(III)-ligated symmetric and asymmetric bacteriochlorin dyads: A computational study

The excited-state properties and photoinduced charge-transfer (CT) kinetics in a series of symmetrical and asymmetrical Zn- and Au-ligated meso–meso-connected bacteriochlorin (BChl) complexes are studied computationally. BChl derivatives, which are excellent near-IR absorbing chromophores, are found to play a central role in bacterial photosynthetic reaction centers but are rarely used in artificial solar energy harvesting systems. The optical properties of chemically linked BChl complexes can be tuned by varying the linking group and involving different ligated metal ions. We investigate charge transfer in BChl dyads that are either directly linked or through a phenylene ring (1,4- phenylene) and which are ligating Zn or Au ions. The directly linked dyads with a nearly perpendicular arrangement of the BChl units bear markedly different properties than phenylene linked dyads. In addition, we find that the dielectric dependence of the intramolecular CT rate is very strong in neutral Zn-ligated dyads, whereas cationic Au-ligated dyads show negligible dielectric dependence of the CT rate. Rate constants of the photo induced CT process are calculated at the semiclassical Marcus level and are compared to fully quantum mechanical Fermi’s golden rule based values. The rates are calculated using a screened range separated hybrid functional that offers a consistent framework for addressing environment polarization. We study solvated systems in two solvents of a low and a high scalar dielectric constant.

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

Computational study of the ECN spray C via one-way coupling of internal nozzle flow and ensuing spray

Here, in this study, computational fluid dynamics simulations of the Engine Combustion Network Spray C injector were performed to investigate the effect of the internal nozzle flow on the development of the ensuing spray in automotive fuel injectors. A set of simulation best practices was first developed to achieve an accurate representation of the internal nozzle flow. Detailed, spatiotemporally defined information was extracted at the injector’s orifice exit and stored in maps that were subsequently employed to initialize a Lagrangian spray by means of a static, one-way coupling (OWC) approach. To highlight the impact of the in-nozzle flow on the spray, the OWC approach was compared with the conventional rate-of-injection (ROI) approach, for which the spray parcels are commonly initialized using a blob injection model. The results indicated that, using the ROI approach, considerable evaporation was observed within the plume in the near-nozzle region. On the other hand, the OWC approach showed evidence of evaporation farther downstream and at the plume periphery. The different behavior observed with the two approaches was mostly due to substantial differences in fuel mass distribution. The liquid fuel concentration predicted by the OWC approach matched well with X-ray tomography measurements of projected fuel mass per unit volume, in terms of both magnitude and distribution. In contrast, an unphysical distribution of fuel mass characterized by an extremely low liquid concentration was predicted by the ROI approach. In the far-field, the effect of the in-nozzle flow was found to be weaker as the spray penetrated and interacted with the ambient gas. In this region, the OWC and ROI approaches provided similar predictions of vapor penetration, spray morphology, and mixing-related characteristics.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

A Computational Study of RNA Tetraloop Thermodynamics, Including Misfolded States

An important characteristic of RNA folding is the adoption of alternative configurations of similar stability, often referred to as misfolded configurations. These configurations are considered to compete with correctly folded configurations, although their rigorous thermodynamic and structural characterization remains elusive. Tetraloop motifs found in large ribozymes are ideal systems for an atomistically detailed computational quantification of folding free energy landscapes and the structural characterization of their constituent free energy basins, including nonnative states. In this work, we studied a group of closely related 10-mer tetraloops using a combined parallel tempering and metadynamics technique that allows a reliable sampling of the free energy landscapes, requiring only knowledge that the stem folds into a canonical A-RNA configuration. Here we isolated and analyzed unfolded, folded, and misfolded populations that correspond to different free energy basins. We identified a distinct misfolded state that has a stability very close to that of the correctly folded state. This misfolded state contains a predominant population that shares the same structural features across all tetraloops studied here and lacks the noncanonical A-G base pair in its loop portion. Further analysis performed with biased trajectories showed that although this competitive misfolded state is not an essential intermediate, it is visited in most of the transitions from unfolded to correctly folded states. Moreover, the tetraloops can transition from this misfolded state to the correctly folded state without requiring extensive unfolding.

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

Enabling Efficient and Accurate Computational Studies of MOF Reactivity via QM/MM and QM/QM Methods

Electronic structure calculations can provide unique insight into metal–organic framework (MOF) reactivity and defect formation. Such calculations can be broadly categorized as utilizing either periodic or cluster models, each with their respective advantages and disadvantages. In the present study, we demonstrate how multiscale methods can leverage the advantages of both approaches to enable high levels of accuracy and computational efficiency in studies of MOF reactivity. Using defect formation in a zeolitic imidazolate framework (ZIF) as a prototypical example, we benchmark a quantum mechanics/molecular mechanics (QM/MM) scheme that enables highly efficient cluster-based calculations on MOFs. We demonstrate the importance of correctly accounting for the influence of both dative bond cleavage in the QM cluster and long-range mechanical coupling to the bulk to achieving accurate QM/MM studies of MOFs. We subsequently leverage these cluster models in a QM/QM scheme that goes beyond standard DFT to yield “gold-standard” correlated wave function results on MOFs at modest computational cost. Crucially, we find several cases in which the incorporation of these correlated corrections yields qualitatively important corrections over conventional DFT values.

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

Computational study of the effect of core–skin structure on the mechanical properties of carbon nanofibers

The effect of the core–skin structure on the mechanical properties of carbon nanofibers is investigated in large-scale molecular dynamics simulations of tensile deformation of carbon nanofibers with the core–skin and homogeneous structures. Contrary to an established notion of the deleterious effect of the skin layer on the strength of carbon fibers, the presence of a high-quality skin layer is found to increase both the Young’s modulus and tensile strength of the nanofiber. A detailed analysis of the fracture process indicates that the nanofiber strengthening is related to the ability of skin layer to suppress crack nucleation at the core–skin interface. Furthermore, the computational predictions suggest that the design of new approaches to carbon fiber manufacturing and processing leading to the generation of a high-quality skin layer while avoiding the introduction of structural defects at the core–skin interface may yield a significant enhancement of the mechanical properties of carbon fibers.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Catalytic decomposition of methane into hydrogen and high-value carbons: combined experimental and DFT computational study

Thermocatalytic decomposition (TCD) of methane can produce hydrogen and valuable nanocarbon co-products with low to near-zero CO2 emission. In this study, a series of Pd promoted Ni catalysts, prepared with various Ni/Pd ratios on a CNT support, were evaluated for methane TCD performance. Characterization and calculations using density functional theory (DFT) were carried out to elucidate the activity–structure relationship and growth mechanism of carbon nanomaterials. It was found that the methane conversion and stability of the catalysts were highly dependent on the Ni/Pd ratio and reaction temperature. DFT calculations revealed that the diffusion of carbon in the metal sublayer required for CNT growth was more favorable in the Ni–Pd alloy lattice suggesting that the buildup of carbon in the metal alloy sublayer facilitated the formation of CNTs and CNFs. A cyclic reaction–regeneration process for self-sustained TCD was experimentally demonstrated. In each cycle, a portion of the separated CNT product was used to re-synthesize the Ni–Pd/CNT catalyst for use in the next reaction cycle. After five cycles of operation, the CH4 conversion, morphology and crystallinity of the carbon product remained unchanged.

Wang, I-Wen↗
DOE OSTI2021

X-ray Crystal Structure of Thorium Tetrahydroborate, Th(BH 4 ) 4 , and Computational Studies of An(BH 4 ) 4 (An = Th, U)

Here the crystal structure of Th(BH 4 ) 4 is described. Two of the four BH 4 – ions are terminal and tridentate (κ 3 ), whereas the other two bridge between neighboring Th IV centers in a κ 2 ,κ 2 (i.e., bis-bidentate) fashion. Thus, each thorium center is bound to six BH 4 – groups by 14 Th–H bonds. The six boron atoms describe a distorted octahedron in which the κ 3 -BH 4 – ions are mutually cis; the 14 ligating hydrogen atoms define a highly distorted bicapped hexagonal antiprism. The thorium centers are linked into a polymer consisting of interconnected helical chains wound about 4-fold screw axes. The structures of An(BH 4 ) 4 (An = Th, U) were also investigated by DFT. The geometries of [An(BH 4 ) 6 ] 2– , [An3(BH 4 ) 16 ] 4– , and [An 5 (BH 4 ) 26 ] 6– fragments of the polymeric structures were optimized at the B3LYP and/or PBE levels. Most calculated geometries are 14-coordinate and agree with the experimental structures, but isolated [Th(BH 4 ) 6 ] 2– units are predicted to feature 16-coordinate Th IV centers.

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

Experimental and Computational Study of Weld Drop-Through Hydraulic Effects in Multi-Layer Concentric Pipes Servicing Cryogenic Liquid Hydrogen

The Cryogenic Moderator System (CMS) at the Oak Ridge National Laboratory (ORNL)-Spallation Neutron Source (SNS) supplies liquid hydrogen, through vacuum insulated transfer lines to the moderators. The welding process used to make the transfer lines produces weld drop-through, which can constrict the flow of hydrogen and/or cause thermal bridges between adjacent layers. To ensure proper operation of the CMS it is therefore necessary to conduct a flow test to verify the pressure drop caused by the welds. Measurements in as-built piping sections intended to provide liquid hydrogen service to the SNS-CMS system were performed to test the effect of weld penetrations into the fluid flow area focusing on pressure drop effects. The reference test sections were designed using prototypical piping dimensions, materials and expected penetration welds. A numerical study using Computational Fluid Dynamics (CFD) tools in the transfer lines provided a conversion factor from the pressure drop measured in the lines with water to the expected pressure drop using liquid cryogenic hydrogen. The results show evidence of a maximum pressure drop modification in water of 78 KPa in the tested parts, which scales to liquid hydrogen as PHydrogen = 5.2 KPa. The maximum measured equivalent blockage area percentage was found in the as-built TDM transfer line spool 2 with a value of 31%. The measured values obtained in the as-built transfer lines provided accurate values to estimate the expected total pressure drop in the CMS system providing operational limits for the current recirculatory capacity. The effect of the weld drop-through in hydrogen evacuation capacity scenarios was better evaluated as a result of this study.

Dominguez-Ontiveros, Elvis [ORNL] (ORCID:000000018↗
DOE OSTI2021

Unsaturated Sulfur Crown Ethers Can Extract Mercury(II) and Show Promise for Future Copernicium(II) Studies: A Combined Experimental and Computational Study

The unsaturated hexathia-18-crown-6 (UHT18C6) molecule was investigated for the extraction of Hg(II) in HCl and HNO3 media. This extractant can be directly compared to the recently studied saturated hexathia-18-crown-6 (HT18C6). The default conformation of the S lone pairs in UHT18C6 is endodentate, where the pocket of the charge density, according to the crystal structures, is oriented toward the center of the ring, which should allow better extraction for Hg(II) compared to the exodentate HT18C6. Batch study experiments showed that Hg(II) had better extraction at low acid molarity (ca. 99% in HCl and ca. 95% in HNO 3 ), while almost no extraction was observed above 0.4 M HCl and 4 M HNO 3 (<5%). Speciation studies were conducted with the goal of delineating a plausible extraction mechanism. Density functional theory calculations including relativistic effects were carried out on both Hg(II)-encapsulated HT18C6 and UHT18C6 complexes to shed light on the binding strength and the nature of bonding. Our calculations offer insights into the extraction mechanism. In addition to Hg(II), calculations were performed on the hypothetical divalent Cn(II) ion, and showed that HT18C6 and UHT18C6 could extract Cn(II). Finally, the extraction kinetics were explored to assess whether this crown can extract the short-lived Cn(II) species in a future online experiment.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2020

Computational study of TATB

We have performed a computational study of TATB. The study is composed of two parts, one where we perform static T = 0 K calculations with several different DFT functionals, to investigate structural properties and obtain a cold curve. Even though the functionals used in this study give poor results at ambient pressure, they perform much better at compression, when the dominance of the van der Waals’ forces in the binding is replaced by more normal interactions. In particular AM05 gives good results, mainly because it doesn’t include any “faulty” van der Waals’. Indeed the van der Waals’ corrected functional (PBE+D2) results shows the worst performance on structural properties. For the cold curve the PBE+D2 functional gives the best results compared to experiments (room temperature) but since the structural properties were bad with this functional we cannot for sure say that this is in any way a validation. It is shown that scaling the AM05 results with a common factor reproduces the experimental data. We leave the question on why to a follow up study but speculate that maybe the experiment and our calculations do not represent the same system set-up. From the structural properties we identify a site in the TATB lattice that contains a high concentration of oxygen and hydrogen. The consequences of this finding, if any, are left to a future project. In the second part of the study we perform extensive high-quality DFT-MD (molecular dynamics with DFT forces) calculations to obtain points on the Hugoniot. This part of the study did not get finished within this project and a more extensive summary of this part will be given in a follow up report at a future date. However, we show that obtaining calculated data is essential for equations of state development since the existing experimental Hugoniot data barely deviates from room temperature isotherm data and thus give no information of the high temperature dependency of materials properties. Our conclusion so far in the project is that being able to calculate materials properties is essential for high explosives since experiments are unable to give information in the high temperature part of phase space. While blind use of DFT cannot (yet) give unambiguous answers, we can still explore qualitatively and use such information to guide us in EOS development. However, the main conclusion is that we need better functionals in DFT to easily extract needed properties from calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

EXPERIMENTAL AND COMPUTATIONAL STUDY OF DIRECT CONTACT MEMBRANE DISTILLATION

Experimental and computational studies of direct contact membrane distillation (DCMD) are conducted. The permeate flux is measured in a module containing a flat sheet membrane. Polytetrafluoroethylene (PTFE) hydrophobic membrane (Membrane Solutions, FPB045A16) with a nominal pore size of 0.45 μm, the thickness of 179~239 μm, the porosity of 80% is used in the experiments. The feed solution is a binary mixture of water and NaCl. The inlet feed temperature is fixed at 70°C, and the feed Reynolds number is varied between 500 and 2000. Two inlet feed concentrations, 0 and 35g/L, are considered. The permeate flow rate was kept constant at Re of 1300 and temperature of 20°C. Three-dimensional steady-state computational fluid dynamics (CFD) simulations are performed using the laminar model. The Navier’s Stokes, energy, and mass transport equations in each channel coupled with flux boundary conditions imposed at the membrane surface are solved. The predicted module-averaged flux is compared to measured data for all operating conditions. Measured and predicted flux agrees well; validating both experiments and model, Temperature polarization in the feed and permeate channel and concentration polarization in the feed channel are characterized by validated CFD simulations.

Oztekin, Alparslan↗
DOE OSTI2025

Electronic, magnetic, and structural properties of V2CoAl: Experimental and computational study

Here, we present results of combined experimental and computations study of V2CoAl, a Heusler alloy that exhibits nearly perfect spin-polarization. Our calculations indicate that this material maintains a high degree of spin-polarization (over 90%) in the wide range of lattice parameters, except at the largest considered unit cell volume. The magnetic alignment of V2CoAl is ferrimagnetic, due to the antialignment of the magnetic moments of vanadium atoms in their two sublattices. The calculated total magnetic moment per formula unit is nearly integer at the optimal lattice parameter and at the smaller volumes of the unit cell, but it deviated from the integer values as the unit cell expands. This is consistent with the calculated variation in the degree of spin polarization with lattice constant. The expected ferrimagnetic behavior has been observed in the arc-melted V2CoAl sample, with a Curie temperature of about 80 K. However, the saturation magnetization is significantly smaller than the theoretical prediction of ∼2 μB/f.u., most likely due to the observed B2-type atomic disorder. The samples exhibit metallic electron transport across the measurement range of 2 K to 300 K.

Kharel, Parashu (ORCID:0000000171330718)↗
DOE OSTI2025

Computational Study of a Cryopump Design for Fusion Exhaust Gas Purification Using Direct Simulation Monte Carlo

Cryopump-based direct internal recycling (DIR) of fusion fuel is an attractive prospect because it provides pumping of helium ash from a reactor along with separation of helium and other impurities from the fuel. Previous studies have demonstrated a continuously regenerating cryopump, referred to as the Snail pump, that separates helium ash from fusion fuel for reactor-relevant flow rates. Here, in this study, a conceptual cryopump was designed as a proof-of-principle study to target other impurities, besides helium ash, from the fusion exhaust that would complement the Snail pump. The impurity-removal cryopump consists of four sets of chevron fins, two sets operating at 80 K and the other two at 30 K. A computational study using direct simulation Monte Carlo (DSMC) was performed with a gas mixture of 96.5% D2, 2.0% He, and 0.5% each of CO2, N2, and CH4. In this configuration, 80 K chevron fin sets are capable of capturing CO2 and 30 K sets capable of capturing CO2, N2, and CH4. The computational study showed that the cryopump is capable of reducing the impurity content by more than two orders of magnitude from the flow.

Adhikari, Nirajan [Oak Ridge National Laboratory (↗
DOE OSTI2021

Reviewing computational studies of defect formation and behaviors in carbon fiber structural units

Solid-state carbon materials, such as graphite and graphene, are at the forefront of materials research because of their unique electronic, vibrational, and mechanical properties, leading to a broad range of potential and realized applications. One key application is their role as basic structural units of carbon fiber (CF), a lightweight alternative to steel. In CF, a delicate relationship exists between ultimate material strength and the atomic-scale density of defects contained at the inter- and intra-subunit levels. Computational studies provide insight into the stability of various types of defects that can form in these systems and connect with experimental observables such as bandgap and spectroscopic measurements. Therefore, the literature contains many computational studies that focus on changes induced by defects, including vacancies and Dienes transformations (Stone-Wales and Thrower defects). However, wide-ranging methods and cell sizes have been used, and property-specific information is often lacking. Furthermore, this review summarizes key literature findings, paying particular attention to changes in the electronic, vibrational, and mechanical properties induced by defects in graphitic materials relevant to CF.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dissociating the phononic, magnetic and electronic contributions to thermal conductivity: a computational study in alpha-iron

Computational tools to study thermodynamic properties of magnetic materials have, until recently, been limited to phenomenological modeling or to small domain sizes limiting our mechanistic understanding of thermal transport in ferromagnets. Herein, we study the interplay of phonon and magnetic spin contributions to the thermal conductivity in a-iron utilizing non-equilibrium molecular dynamics simulations. It was observed that the magnetic spin contribution to the total thermal conductivity exceeds lattice transport for temperatures up to two-thirds of the Curie temperature after which only strongly coupled magnon-phonon modes become active heat carriers. Characterizations of the phonon and magnon spectra give a detailed insight into the coupling between these heat carriers, and the temperature sensitivity of these coupled systems. Comparisons to both experiments and ab initio data support our inferred electronic thermal conductivity, supporting the coupled molecular dynamics/spin dynamics framework as a viable method to extend the predictive capability for magnetic material properties.

36 MATERIALS SCIENCE↗