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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

The Nature of the Long-Lived Excited State in a Ni II Phthalocyanine Complex Investigated by X-Ray Transient Absorption Spectroscopy

The nature of the photoexcited state of octabutoxy nickel(II) phthalocyanine (NiPcOBu8) with a 500ps lifetime was investigated by X-ray transient absorption (XTA) spectroscopy. Previous optical, vibrational, and computational studies have suggested that this photoexcited state has a ligand-to-metal charge transfer (LMCT) nature. By using XTA, which provides unambiguous information on the local electronic and nuclear configuration around the Ni center, the nature of the excited state of NiPcOBu8 was reassessed. Using X-ray probe pulses from a synchrotron source, the ground- and excited-state X-ray absorption spectra of NiPcOBu 8 were measured. Based on the results, we identified that the excited state exhibits spectral features that are characteristic of a Ni-1,Ni-3(3d(z2),3d(x2-y2)) state rather than a LMCT state with a transiently reduced Ni center. This state resembles the (d,d) state of nickel(II) tetramesitylphorphyrin. The XTA features are rationalized based on the inherent cavity sizes of the macrocycles. These results may provide useful guidance for the design of photocatalysts in the future.

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

Calculations of Solid-State 43Ca NMR Parameters: A Comparison of Periodic and Cluster Approaches and an Evaluation of DFT Functionals

We present a computational study of magnetic-shielding and quadrupolar-coupling tensors of 43Ca sites in crystalline solids. A comparison between periodic and cluster-based approaches for modeling solid-state interactions demonstrates that clusterbased approaches are suitable for predicting 43Ca NMR parameters. Several model chemistries, including Hartree–Fock theory and 17 DFT approximations (SVWN, CA-PZ, PBE, PBE0, PW91, B3PW91, rPBE, PBEsol, WC, PKZB, BMK, M06-L, M06, M06-2X, M06-HF, TPSS, and TPSSh), are evaluated for the prediction of 43Ca NMR parameters. Convergence of NMR parameters with respect to basis sets of the form cc-pVXZ (X5D, T, Q) is also evaluated. All DFT methods lead to substantial, and frequently systematic, overestimations of experimental chemical shifts. Hartree–Fock calculations outperform all DFT methods for the prediction of 43Ca chemical-shift tensors. VC 2017 Wiley Periodicals, Inc. DOI: 10.1002/jcc.24763

Holmes, Sean T.↗
DOE OSTI2024

Shortest path network interdiction with asymmetric uncertainty

Abstract This paper considers an extension of the shortest path network interdiction problem that incorporates robustness to account for parameter uncertainty. The shortest path interdiction problem is a game of two players with conflicting agendas and capabilities: an evader, who traverses the arcs of a network from a source node to a sink node using a path of shortest length, and an interdictor, who maximizes the length of the evader's shortest path by interdicting arcs on the network. It is usually assumed that the parameters defining the network are known exactly by both players. We consider the situation where the evader assumes the nominal parameter values while the interdictor uses robust optimization techniques to account for parameter uncertainty or sensor degradation. We formulate this problem as a nonlinear mixed‐integer semi‐infinite bilevel program and show that it can be converted into a mixed‐integer linear program with a second order cone constraint. We use random geometric networks and transportation networks to perform computational studies and demonstrate the unique decision strategies that our variant produces. Solving the shortest path interdiction problem with asymmetric uncertainty protects the interdictor from investing in a strategy that hinges on key interdictions performing as promised. It also provides an alternate strategy that mitigates the risk of these worst‐case possibilities.

Punla‐Green, She'ifa Z.↗
DOE OSTI2020

Sorption-induced polymer rearrangement: approaches from molecular modeling

With the growing need for chemical separation and chemical storage solutions, polymeric adsorbents have emerged as a promising class of candidate materials because of their potentially tunable sorption properties, membrane structure and relatively cost consciousness. Moreover, the developing field of polymeric membrane materials has shown particular success at integrating both experimental and computational studies. However, these material systems are known to suffer from varying degrees of induced membrane structural rearrangement upon adsorbate uptake, and thus many polymeric membrane performance metrics are often considered to degrade with an increasing number of ‘guest’ species. In this mini–review, we highlight methodology tradeoffs and provide insights into atomistic molecular simulations used to study adsorption with flexible frameworks, which have the potential to predict separation, storage or catalytic capabilities a priori to experimental efforts. Specifically, molecular simulation methods that have been applied to provide predictions of polymeric membrane properties that have included consideration for sorbate–induced polymer chain rearrangement, swelling and/or plasticization are reviewed. Here, the examples and methodologies described provide demonstrations of the applicability of simulations as an approach to understand adsorption–based phenomena at an atomistic/molecular level, and as a tool to carry out screening studies aimed at efficiently providing analysis for a diversity of polymeric adsorbent–adsorbate systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Elucidating the Role of Hydrogen at c‐In 2 O 3 /a‐In 2 O 3− x Interface

Ab initio molecular dynamics simulations and hybrid functional electronic structure calculations are employed to determine the formation, the structural and electronic properties, and the dynamics of covalent (In—OH) and ionic (In—H—In) hydrogen defects at crystalline-In 2 O 3 /amorphous-In 2 O 3−x interface. This comprehensive computational study considers (i) various interstitial and substitutional hydrogen site locations within the crystalline, amorphous, and interfacial regions; (ii) several oxygen-to-hydrogen ratios; and (iii) possible defect charged states. The results reveal hydrogen's inability to fully passivate the undercoordinated under-shared in atoms in amorphous highly substoichiometric oxide, giving rise to the formation of deep electron traps even in net-charge neutral structures. These trap defects are found to be sensitive to photoexcitation, in contrast to In—OH with H's electronic states located below the valence band and to In—H—In, where H is found to maintain its charge state upon illumination. Nevertheless, H plays a critical role in photoinduced conductivity and its relaxation by promoting In—O coordination transformations at the interfacial region, including deterioration of the crystalline layer. Finally, the results help identify mobile H species and metastable H defect complexes (such as In—H—H—In, In—OH–H—In, and In—OH–O—In) that are responsible for long relaxation times of the conductivity decay.

H defects↗
DOE OSTI2023

Underlying Substrate Effect on Electrochemical Activity for Hydrogen Evolution Reaction with Low‐Platinum‐Loaded Catalysts

Platinum is known as the best catalyst for the hydrogen evolution reaction (HER) but the scarcity and high cost of Pt limit its widespread applicability. Herein, the role of the underlying substrate on the HER activity of dispersed Pt atoms is uncovered. A direct current magnetron sputtering technique is utilized to deposit transition metal (TM) thin films of W, Ti, and Ta as underlying substrates for extremely low loading of Pt (<1.5 at%). The electrocatalytic performance of as‐synthesized samples for the HER is examined in both alkali and acidic media. The results show that despite the low loading of Pt, the Pt/TM catalysts produce hydrogen at a rate comparable to that of pristine bulk Pt. Pt/TM catalysts also display good stability with less than 5% decay in performance after 10 h of continuous HER operation. Based on the computational study, the excellent performance is attributed to the modified electronic properties of the Pt atoms, offering ideal binding energy for HER due to interaction with the underlying substrates. This work provides a robust and industry‐friendly route toward designing efficient catalytic systems for important electrochemical reactions such as HER and others.

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

Quantifying modeling uncertainty in simplified beam models for building response prediction

The use of simple models for response prediction of building structures is preferred in earthquake engineering for risk evaluations at regional scales, as they make computational studies more feasible. The primary impediment in their gainful use presently is the lack of viable methods for quantifying (and reducing upon) the modeling errors/uncertainties they bear. This study presents a Bayesian calibration method wherein the modeling error is embedded into the parameters of the model. Here, the method is specifically described for coupled shear-flexural beam models here, but it can be applied to any parametric surrogate model. The major benefit the method offers is the ability to consider the modeling uncertainty in the forward prediction of any degree-of-freedom or composite response regardless of the data used in calibration. The method is extensively verified using two synthetic examples. In the first example, the beam model is calibrated to represent a similar beam model but with enforced modeling errors. In the second example, the beam model is used to represent the detailed finite element model of a 52-story building. Both examples show the capability of the proposed solution to provide realistic uncertainty estimation around the mean prediction.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

When finite becomes infinite: convergence properties of vibrational spectra of oligomer chains

We present a computational study of convergence properties of vibrational IR and Raman spectra for a series of increasingly long units of polyethylene, cis- and trans-polyacetylenes, and polyynes. Convergent behavior to the spectra of infinitely long polymers was observed in all cases when chains reached lengths of approximately 60 carbon atoms, both with respect to the positions of the bands and to their intensities. The vibrational spectra of longer chains are practically indistinguishable. The convergence rate depends on the degree of the π conjugation in a studied system: Vibrational spectra for oligoethylenes converge noticeably faster than the spectra for the conjugated systems. Here, the slowest convergence is observed for skeletal motions of the oligomer chains, which may require more than a hundred carbon atoms in the chain to show deviations smaller than 1 cm –1 to the corresponding solid-state calculations. The results suggest that the boundary between the properties of finite and infinite molecular systems fades away for a surprisingly small number of atoms.

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

Experimental and Computational Investigation of the Emission of Radioactive Products from Irradiated Fuel of Fast Sodium-Cooled Reactor

The key experimental and computational studies of the emission of long-lived radioactive fission products and actinides from uranium dioxide irradiated in BN-600 into sodium and the gas phase in the range 150–800°C, performed on the experimental bench at the Leipunskii Institute of Physics and Power Engineering, are presented. The results of measurements of the properties of aerosols of certain radioactive fission products formed upon heating of irradiated nuclear fuel in sodium are presented. The experimental results are of interest for evaluating the radiation impact of reactor facilities on workers and the public during beyond design basis accidents.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Accident Resistance of Molten-Salt Nuclear Reactor

The aim of the present work is to examine the possibilities of the MSR-I molten-salt nuclear actinide-incinerator reactor with a cavity-type core as possessing internal and passive safety characteristics and meeting the requirements imposed on next-generation reactors. Starting from the premises that the basic principles of nuclear safety must not differ from all other types of reactors the MSR-I reactor was analyzed in transient processes without activation of the accident protection system in situations where fuel could not be discharged into subcritical drainage tanks because of operator error or failure of the emergency discharge system. A computational study showed that the basic types of emergency situations can be overcome without destroying the protective barriers or actuating accident protection systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

Evolution of Microstructure and Creep Behavior in an Fe-Ni-Cr-Nb-C Alloy during Service in Hydrocarbon Cracker Tubes

Heat resistant, cast austenitic Fe-Ni-Cr-Nb-C alloy tubes are used by the petrochemical industry to transport hydrocarbons at temperatures between 850 and 1050 °C. Exposure to carbon from the process stream and nitrogen from the air causes severe carburization and nitridation in the tube walls, which influences creep resistance and degrades their ductility. A better understanding of these processes would enable their service life to be extended. To this end, we report the results of an experimental and computational study of the changes in microstructure and creep behavior in a representative alloy during service. The microstructure of the fresh material is compared to that of 8- and 10-year-old tubes. The results confirm the presence of M{sub 23}C{sub 6}, M{sub 7}C{sub 3}, and M{sub 2}(CN) phases reported in prior studies, and also reveal the nucleation and growth of the carbo-nitride phase M{sub 6}(CN). We show that the spatial variations of carbide and nitride particle volume fractions in aged tubes can be predicted by a simple reaction–diffusion equation. In addition, the uniaxial compressive creep response of the fresh and aged tubes was obtained as functions of temperature and strain rate. The materials deform by power-law creep, with strain rate versus stress relation as ε-dot =Aσ{sup n}. Aging reduces the stress exponent from n = 8 in fresh tubes to n = 5 in aged tubes, and also reduces the pre-exponential coefficient A. Tests on model alloys with composition identical to that of the matrix combined with finite element simulations of creep in model microstructures reveal that Cr depletion is responsible for the change in stress exponent, while the change in A is caused by the increase in particle volume fraction. By combining the microstructure-scale simulations of creep with the predictions of the reaction–diffusion model of microstructure evolution, the creep behavior of tubes left in service beyond 10 years is estimated.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Ligation Kinetics as a Probe for Non-Covalent Electrostatic Bonding and Electron Solvation of Alkali and Alkaline Earth Cations with Ammonia

Mono-ligation kinetics were measured for ammonia reacting with atomic cations in the first two groups of the periodic table (K{sup +}, Rb{sup +}, Cs{sup +} and Ca{sup +}, Sr{sup +}, Ba{sup +}). Also, mono-ligation energies were computed using density functional theory (DFT) in an attempt to assess the role of non-covalent electrostatic interactions in these chemical reactions. The measurements were performed at room temperature in helium bath gas at 0.35 Torr using an inductively coupled plasma/selected-ion flow tube (ICP/SIFT) tandem mass spectrometer. Rate coefficients are reported for ammonia addition, the only reaction channel that was observed with all these cations. A systematic decrease in the rate of addition of NH{sub 3} was observed for both group 1 and 2 cations going down the periodic table. The computational studies predict a decrease in the adduct binding energy and an increase in the bond separation going down groups 1 and 2 of the periodic table and provide some insight into the role of the extra selectron in the group 2 radical cations in ligand bonding. A correlation is seen between the efficiency of ligation and the binding energy of the adduct ion and attributed to the lifetime of the intermediate encounter complex against back dissociation which is dependent on its well depth. Higher-order additions of ammonia were also observed. Remarkable differences in the extent and kinetics were seen between the group 1 and 2 cations, and these were attributed to the occurrence of ammonia solvation of the extra s electron in the higher-order adducts of the alkaline earth cations. .

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2018

GeV ion irradiation of NiFe and NiCo: Insights from MD simulations and experiments

Concentrated solid solution alloys have attracted rapidly increasing attention due to their potential for designing materials with high tolerance to radiation damage. To tackle the effects of chemical complexity in defect dynamics and radiation response, we present in this paper a computational study on swift heavy ion induced effects in Ni and equiatomic Ni -based alloys (Ni 50 Fe 50 , Ni 50 Co 50 ) using two-temperature molecular dynamics simulations (2T-MD). The electronic heat conductivity in the two-temperature equations is parameterized from the results of first principles electronic structure calculations. A bismuth ion (1.542 GeV) is selected and single impact simulations performed in each target. We study the heat flow in the electronic subsystem and show that alloying Ni with Co or Fe reduces the heat dissipation from the impact by the electronic subsystem. Simulation results suggest no melting or residual damage in pure Ni while a cylindrical region melts along the ion propagation path in the alloys. In Ni 50 Co 50 the damage consists of a dislocation loop structure (d = 2 nm) and isolated point defects, while in Ni 50 Fe 50 , a defect cluster (d = 4 nm) along the ion path is, in addition, formed. The simulation results are supported by atomic-level structural and defect characterizations in bismuth-irradiated Ni and Ni 50 Fe 50 . Finally, the significance of the 2T-MD model is demonstrated by comparing the results to those obtained with an instantaneous energy deposition model without consideration of e-ph interactions in pure Ni and by showing that it leads to a different qualitative behavior.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nanotwinned alloys under high pressure

Nanotwinned alloys are of interest due to their high strength and ductility, but twin boundaries may not be stable under shear. Computational studies indicate that high hydrostatic pressure may suppress detwinning mechanisms. Here, in this study, we investigate the microstructural changes of nanotwinned-nanocrystalline copper-nickel and Inconel 725 alloys under quasi-hydrostatic pressures up to 50 gigapascals (GPa). The alloys are compressed in a diamond anvil cell. In-situ x-ray diffraction (XRD) and ex-situ transmission electron microscopy (TEM) were employed to monitor microstructural changes. Twin boundary deformation and grain growth occur at 11.4 GPa quasi-hydrostatic pressure in the copper-nickel alloy. Molecular dynamics (MD) simulations reveal that hydrostatic pressure causes elevated local shear stress at grain boundaries, which leads to atomic rearrangements. A superposition of hydrostatic and deviatoric pressures lead to partial dislocation mediated twin boundary migration. In contrast, the Inconel 725 alloy showed stable twin and grain boundaries up to a quasi-hydrostatic pressure of 12.7 GPa. Texture, high solid solution strengthening, and low stacking fault energy are hypothesized to the enhanced microstructural stability in Inconel 725.

36 MATERIALS SCIENCE↗
DOE OSTI2015

At least 10% shorter C–H bonds in cryogenic protein crystal structures than in current AMBER forcefields

High resolution protein crystal structures resolved with X-ray diffraction data at cryogenic temperature are commonly used as experimental data to refine forcefields and evaluate protein folding simulations. However, it has been unclear hitherto whether the C–H bond lengths in cryogenic protein structures are significantly different from those defined in forcefields to affect protein folding simulations. This article reports the finding that the C–H bonds in high resolution cryogenic protein structures are 10–14% shorter than those defined in current AMBER forcefields, according to 3709 C–H bonds in the cryogenic protein structures with resolutions of 0.62–0.79 Å. Also, 20 all-atom, isothermal–isobaric, 0.5-μs molecular dynamics simulations showed that chignolin folded from a fully-extended backbone formation to the native β-hairpin conformation in the simulations using AMBER forcefield FF12SB at 300 K with an aggregated native state population including standard error of 10 ± 4%. However, the aggregated native state population with standard error reduced to 3 ± 2% in the same simulations except that C–H bonds were shortened by 10–14%. Furthermore, the aggregated native state populations with standard errors increased to 35 ± 3% and 26 ± 3% when using FF12MC, which is based on AMBER forcefield FF99, with and without the shortened C–H bonds, respectively. These results show that the 10–14% bond length differences can significantly affect protein folding simulations and suggest that re-parameterization of C–H bonds according to the cryogenic structures could improve the ability of a forcefield to fold proteins in molecular dynamics simulations. - Highlights: • Cryogenic crystal structures are commonly used in computational studies of proteins. • C–H bonds in the cryogenic structures are shorter than those defined in forcefields. • A survey of 3709 C–H bonds shows that the cryogenic bonds are 10–14% shorter. • The 10–14% differences affect molecular dynamics simulations of protein folding. • Shortening C–H bonds by 10–14% can improve protein folding simulations.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Quantifying Shear-Induced Deformation and Detachment of Individual Adherent Sickle Red Blood Cells

According to the article, vaso-occlusive crisis, a common painful complication of sickle cell disease, is a complex process triggered by intercellular adhesive interactions among blood cells and the endothelium in all human organs (e.g., the oxygen-rich lung as well as hypoxic systems such as liver and kidneys). We present a combined experimental-computational study to quantify the adhesive characteristics of sickle mature erythrocytes (SMEs) and irreversibly sickled cells (ISCs) under flow conditions mimicking those in postcapillary venules. We employed an in vitro microfluidic cell adherence assay, which is coated uniformly with fibronectin. We investigated the adhesion dynamics of SMEs and ISCs in pulsatile flow under well-controlled hypoxic conditions, inferring the cell adhesion strength by increasing the flow rate (or wall shear stress (WSS)) until the onset of cell detachment. In parallel, we performed simulations of individual SMEs and ISCs under shear. We introduced two metrics to quantify the adhesion process, the cell aspect ratio (AR) as a function of WSS and its rate of change (the dynamic deformability index). We found that the AR of SMEs decreases significantly with the increase of WSS, consistent between the experiments and simulations. In contrast, the AR of ISCs remains constant in time and independent of the flow rate. The critical WSS value for detaching a single SME in oxygenated state is in the range of 3.9–5.5 Pa depending on the number of adhesion sites; the critical WSS value for ISCs is lower than that of SMEs. Our simulations show that the critical WSS value for SMEs in deoxygenated state is above 6.2 Pa (multiple adhesion sites), which is greater than their oxygenated counterparts. We investigated the effect of cell shear modulus on the detachment process; we found that for the same cell adhesion spring constant, the higher shear modulus leads to an earlier cell detachment from the functionalized surface. These findings may aid in the understanding of individual roles of sickle cell types in sickle cell disease vaso-occlusion.

Biophysics↗