SEARCH · Search NASA

Results for “Computational methods”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

All sources NASA NTRS DOE OSTI

DOE OSTI2018

Theoretical Studies on Magnetic Structures, Hysteresis Loops and Size Effects of a Pair of Frustrated Double-Walled Nanotubes

We propose a theoretical quantum model and derive a set of analytic formulas to study the physical properties of a pair of double-walled magnetic nanotubes. The Heisenberg exchange parameters between the two walls of the nanotubes are assumed to differ only in sign. Thus, in the absence of external magnetic field, our calculated macroscopic properties of this pair of nanotubes are almost precisely identical, exhibiting fascinating duality of the nanosystems and demonstrating the correctness of our theoretical model. The two spin systems are all frustrated, so that sudden changes in the macroscopic properties are observed around T{sub M2} that is well below the transition temperature T{sub M1}. However, only the inner shell consisting of smaller A-type spins has been obviously affected. In the temperature range T{sub M2} < T < T{sub M1}, this shell becomes semi-antiferromagnetic and its magnetization is considerably suppressed, whereas as temperature falls below T{sub M2} the shell gradually restores its ferromagnetic nature. The longitudinal hysteresis behavior of such a double-waled nanotube is ferromagnetic-like below T{sub M2}, but antiferromagnetic-like in the temperature interval T{sub M2} < T < T{sub M1}. Moreover, we find that the diameter of the nanotube has seemly no effects on its physical properties, whereas its length does affects the two temperatures slightly, and also its spin configuration at very low temperatures if the tube is sufficiently long. More importantly, the theoretical results presented in this paper can be precisely reproduced with the quantum computational method we develop in recent years, justifying the validity of the numerical approach once again.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Modeling in High Temperature Corrosion: A Review and Outlook

Realizing higher operating temperatures to increase efficiency of future applications for energy conversion and storage while minimizing cost is a challenge for development of high-temperature materials. Simultaneous optimization of mechanical properties and corrosion resistance continues to be a difficult task but is essential due to the need to significantly accelerate the transition between technology readiness levels in the future. Furthermore, oxidation-induced degradation will be a critical life-limiting mechanism at increased operating temperatures. Suitable high-temperature materials cannot be solely identified by time-consuming experiments and reliable computational methods incorporating the relevant physics of processes must be considered to complement the experimental efforts. In the present work, a review of the methods employed to model oxidation-induced material degradation described in literature will be discussed. Furthermore, their capability to predict lifetime and aid in material selection will be evaluated.

36 MATERIALS SCIENCE↗

DOE OSTI2019

Characterization and Modeling of Hot Deformation Behavior of a Copper-Bearing High-Strength Low-Carbon Steel Microalloyed with Nb

This study investigates hot deformation behavior of a newly developed Cu-bearing high-strength low-carbon steel microalloyed with Nb (Nb-HSLC). A computational method based on experimental data was employed to design the chemical composition of the alloy. Compression tests were carried out in the temperature range of 850-1100 °C as well as strain rates of 0.001-10 s{sup −1} using BAHR Dil 805 A/D thermo-analyzer equipment. The Arrhenius-type constitutive equations were used to model the hot working behavior of the designed steel. Effects of friction and temperature rise during deformation were corrected to obtain the actual stresses. The results showed that the peak flow stress was increased with increasing Zener–Hollomon parameter. The obtained flow curves at strain rates lower than 0.1 s{sup −1} and temperatures above 950 °C represented the typical dynamic recrystallization (DRX) behavior, while the flow curves at temperatures lower than 950 °C at all strain rates were associated with continuous strain hardening. This feature is in good agreement with the precipitation temperature range of Nb(C, N) particles, i.e., 800-1000 °C. Moreover, the flow curves showed the serrations during hot deformation at strain rates of 0.001 and 0.01 s{sup −1}, indicating that the dynamic strain aging (DSA) phenomenon occurred at low strain rates. The best fit between “peak stress” and “deformation conditions” was obtained by a hyperbolic sine-type equation (R{sup 2} = 0.993). Therefore, the average activation energy was determined as 348 kJ mol{sup −1}. The agreement between the achieved model and experimental flow data was verified using the results of additional tests at a strain rate of 5 s{sup −1}. The maximum difference between the measured and predicted “peak stresses” was calculated as 5 Mpa.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Combined Experimental and Computational Study of Molybdenum and Niobium for Nuclear Sensor Application

Due to their novel electromagnetic and thermal properties, molybdenum (Mo) and niobium (Nb) become optimal temperature sensor materials for nuclear energy applications. We leveraged voltage recorded during a heat ramp to tune a computational method to predict the Seebeck electromotive force (EMF) of Mo and Nb. Using a combined Density Functional Theory (DFT) and Boltzmann Transport Equations (BTE) method the voltage was predicted but did not include the effects of temperature on atomic structure. Combining Ab Initio Molecular Dynamics (AIMD) and BTE included temperature effects on structure optimization and yielded voltages in a good agreement with experiment. Lanthanum (La) and Phosphorus (P) additives in Mo and Nb, respectively, could increase the EMF compared to those of the pure metals. The presence of oxygen (O) in Mo increases the EMF while O in Nb slightly reduces the EMF. Furthermore, our studies suggested that heat treatment-induced structural changes that lead to a reduction in voltage occur not only at the mesoscale as previously understood but also at the atomic scale.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2020

Machine learning a molecular Hamiltonian for predicting electron dynamics

We develop a computational method to learn a molecular Hamiltonian matrix from matrix-valued time series of the electron density. As we demonstrate for three small molecules, the resulting Hamiltonians can be used for electron density evolution, producing highly accurate results even when propagating 1,000 time steps beyond the training data. As a more rigorous test, we use the learned Hamiltonians to simulate electron dynamics in the presence of an applied electric field, extrapolating to a problem that is beyond the field-free training data. We find that the resulting electron dynamics predicted by our learned Hamiltonian are in close quantitative agreement with the ground truth. Our method relies on combining a reduced-dimensional, linear statistical model of the Hamiltonian with a time-discretization of the quantum Liouville equation within time-dependent Hartree Fock theory. In conclusion, we train the model using a least-squares solver, avoiding numerous, CPU-intensive optimization steps. For both field-free and field-on problems, we quantify training and propagation errors, highlighting areas for future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Chlorophylls as primary electron acceptors in reaction centers

Photosynthetic Reaction Centers (RCs) can be considered blueprints for highly efficient energy transfer. Embedded with an array of cofactors, including (bacterio)chlorophyll ((B)Chl) and (B)pheophytin ((B)Pheo) molecules, RCs function with a high quantum yield that spans a wide spectral range. Understanding the principles that underlie their function can influence the design of the next generation of artificial photosynthetic devices. We are particularly interested in the factors that influence the early stages of light-driven charge separation in RCs. With the recent publication of several highly anticipated RC structures and advanced computational methods available, it is possible to probe both the geometric and electronic structures of an array of RCs. In this chapter, we review the electronic and geometric structures of the (B)Chl and (B)Pheo primary electron acceptors from fiveRCs, comprising both Type I and Type II RCs and representing both heterodimeric and homodimeric systems. We showcase the dimeric A0•- state of Type I RCs, whereby the unpaired electron is delocalized, to various extents, over two (B)Chl molecules (B)Chl2 and (B)Chl3. This delocalization is controlled by several factors, including the structure of the (B)Chls, interactions with the surrounding protein matrix, and the orientation and distances of the cofactors themselves. In contrast, the primary acceptors of Type II RCs are entirely monomeric, with electron density residing solely on the (B)Pheo. We compare the natural design of the primary acceptors of the Type I and Type II RCs from both an evolutionary and application based perspective.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Interpreting and designing microbial communities for bioprocess applications, from components to interactions to emergent properties

Interest in microbial communities for bioprocessing has surged in recent years based on the potential to optimize multiple tasks simultaneously and to enhance process productivity and stability. The presence and magnitude of these desirable system properties often result from interactions between functionally distinct community members. The importance of interactions, while appreciated by some disciplines for decades, has gained interest recently due to the development of ‘omics techniques, polymicrobial culturing approaches, and computational methods which has made the systems-level analysis of interacting components more tractable. This review defines and categorizes natural and engineered system components, interactions, and emergent properties, as well as presents three ecological theories relevant to microbial communities. Case studies are interpreted to illustrate components, interactions, emergent properties and agreement with theoretical concepts. A general foundation is laid to facilitate interpretation of current systems and to aid in future design of microbial systems for the next generation of bioprocesses.

community↗

DOE OSTI2016

Molecular Dynamics Studies of Proton Transport in Hydrogenase and Hydrogenase Mimics

Protons are used throughout the biological world for a number of functions, from charge balance to energy carriers. Metalloenzymes use protons as energy carriers and control proton movement both temporally and spatially. Despite the interest and need for controlled proton movement in other systems, the scientific community has not been able to develop extensive general rules for developing synthetic proton pathways. In part this is due to the challenging nature of studying these large and complex molecules experimentally, although experiments have gleaned extensive critical insight. While computational methods are also challenging because of the large size of the molecules, they have been critical in advancing our knowledge of proton movement through pathways, but even further, they have advanced our knowledge in how protonation and proton movement is correlated with large and small scale molecular motions and electron movement. These studies often complement experimental studies but provide insight and depth simply not possible in many cases in the absence of theory. In this chapter, we will discuss advances and methods used in understanding proton movement in hydrogenases.

Electron Transport↗

DOE OSTI2024

A semantics-driven framework to enable demand flexibility control applications in real buildings

Decarbonising and digitalising the energy sector requires scalable and interoperable Demand Flexibility (DF) applications. Semantic models are promising technologies for achieving these goals, but existing studies focused on DF applications exhibit limitations. These include dependence on bespoke ontologies, lack of computational methods to generate semantic models, ineffective temporal data management and absence of platforms that use these models to easily develop, configure and deploy controls in real buildings. This paper introduces a semantics-driven framework to enable DF control applications in real buildings. The framework supports the generation of semantic models that adhere to Brick and SAREF while using metadata from Building Information Models (BIM) and Building Automation Systems (BAS). The work also introduces a web platform that leverages these models and an actor and microservices architecture to streamline the development, configuration and deployment of DF controls. The paper demonstrates the framework through a case study, illustrating its ability to integrate diverse data sources, execute DF actuation in a real building, and promote modularity for easy reuse, extension, and customisation of applications. The paper also discusses the alignment between Brick and SAREF, the value of leveraging BIM data sources, and the framework's benefits over existing approaches, demonstrating a 75% reduction in effort for developing, configuring, and deploying building controls.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

A Low-Enriched Uranium (LEU) option for the conversion of FRM II

The Forschungs-Neutronenquelle Heinz Maier-Leibnitz (FRM II) is Germany’s most powerful research reactor and uses Highly Enriched Uranium (HEU) fuel enriched at 93 %. The Technical University of Munich (TUM) operates the reactor and has been mandated to convert the FRM II to a lower-enrichment fuel. In preparation for the lower-enrichment fuel downselection, planned for 2023, TUM is evaluating multiple conversion scenarios. In this paper, it is demonstrated that a conversion of FRM II to Low-Enriched Uranium (LEU) is scientifically possible when using the novel monolithic U-10Mo fuel system, which is currently in the qualification process. As a research reactor, safety criteria, neutron flux distribution and cycle length are key metrics. To begin, neutronic and thermal-hydraulic models were created based on current technical drawings and relevant material properties. Specifically tailored computational methods and coupling schemes have also been developed to properly evaluate all relevant characteristics for the conversion of the FRM II reactor. Using these models and methods, a systematic parameter study was performed to explore the defined design space. The results of that study indicate that many LEU designs using monolithic U-10Mo appear to be viable to convert FRM II to LEU. One example of an LEU solution for the FRM II reactor that fulfills the stated safety, compatibility and scientific performance requirements is discussed in more detail.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2014

Real-time Feynman path integral with Picard–Lefschetz theory and its applications to quantum tunneling

Picard–Lefschetz theory is applied to path integrals of quantum mechanics, in order to compute real-time dynamics directly. After discussing basic properties of real-time path integrals on Lefschetz thimbles, we demonstrate its computational method in a concrete way by solving three simple examples of quantum mechanics. It is applied to quantum mechanics of a double-well potential, and quantum tunneling is discussed. We identify all of the complex saddle points of the classical action, and their properties are discussed in detail. However a big theoretical difficulty turns out to appear in rewriting the original path integral into a sum of path integrals on Lefschetz thimbles. We discuss generality of that problem and mention its importance. Real-time tunneling processes are shown to be described by those complex saddle points, and thus semi-classical description of real-time quantum tunneling becomes possible on solid ground if we could solve that problem. - Highlights: • Real-time path integral is studied based on Picard–Lefschetz theory. • Lucid demonstration is given through simple examples of quantum mechanics. • This technique is applied to quantum mechanics of the double-well potential. • Difficulty for practical applications is revealed, and we discuss its generality. • Quantum tunneling is shown to be closely related to complex classical solutions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Mechanistic insights into mode of action of potent natural antagonists of BACE-1 for checking Alzheimer’s plaque pathology

Highlights: •Accumulation of Aβ plaques is one of the major pathology associated with Alzheimer’s disease. •Inhibition of β-Secretase or BACE-1 offers a viable prospect to check the growth of these plaques. •A large virtual dataset of natural compounds was screened against BACE-1. •Top two hits were analyzed for thermodynamic and structural stability using MD simulations. •Their detailed binding mode of actions were elucidated. -- Abstract: Alzheimer’s is a neurodegenerative disorder resulting in memory loss and decline in cognitive abilities. Accumulation of extracellular beta amyloidal plaques is one of the major pathology associated with this disease. β-Secretase or BACE-1 performs the initial and rate limiting step of amyloidic pathway in which 37–43 amino acid long peptides are generated which aggregate to form plaques. Inhibition of this enzyme offers a viable prospect to check the growth of these plaques. Numerous efforts have been made in recent years for the generation of BACE-1 inhibitors but many of them failed during the preclinical or clinical trials due to drug related or drug induced toxicity. In the present work, we have used computational methods to screen a large dataset of natural compounds to search for small molecules having BACE-1 inhibitory activity with low toxicity to normal cells. Molecular dynamics simulations were performed to analyze molecular interactions between the screened compounds and the active residues of the enzyme. Herein, we report two natural compounds of inhibitory nature active against β-secretase enzyme of amyloidic pathway and are potent lead molecules against Alzheimer’s disease.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

Systematic drug safety evaluation based on public genomic expression (Connectivity Map) data: Myocardial and infectious adverse reactions as application cases

Adverse drug reaction (ADR) is of great importance to both regulatory agencies and the pharmaceutical industry. Various techniques, such as quantitative structure–activity relationship (QSAR) and animal toxicology, are widely used to identify potential risks during the preclinical stage of drug development. Despite these efforts, drugs with safety liabilities can still pass through safety checkpoints and enter the market. This situation raises the concern that conventional chemical structure analysis and phenotypic screening are not sufficient to avoid all clinical adverse events. Genomic expression data following in vitro drug treatments characterize drug actions and thus have become widely used in drug repositioning. In the present study, we explored prediction of ADRs based on the drug-induced gene-expression profiles from cultured human cells in the Connectivity Map (CMap) database. The results showed that drugs inducing comparable ADRs generally lead to similar CMap expression profiles. Based on such ADR-gene expression association, we established prediction models for various ADRs, including severe myocardial and infectious events. Drugs with FDA boxed warnings of safety liability were effectively identified. We therefore suggest that drug-induced gene expression change, in combination with effective computational methods, may provide a new dimension of information to facilitate systematic drug safety evaluation. - Highlights: • Drugs causing common toxicity lead to similar in vitro gene expression changes. • We built a model to predict drug toxicity with drug-specific expression profiles. • Drugs with FDA black box warnings were effectively identified by our model. • In vitro assay can detect severe toxicity in the early stage of drug development.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2015

The Phylogenetic Signature Underlying ATP Synthase c-Ring Compliance

The proton-driven ATP synthase (F O F 1 ) is comprised of two rotary, stepping motors (F O and F 1 ) coupled by an elastic power transmission. The elastic compliance resides in the rotor module that includes the membrane-embedded F O c-ring. Proton transport by F O is firmly coupled to the rotation of the c-ring relative to other F O subunits (ab 2 ). It drives ATP synthesis. We used a computational method to investigate the contribution of the c-ring to the total elastic compliance. We performed principal component analysis of conformational ensembles built using distance constraints from the bovine mitochondrial c-ring x-ray structure. Angular rotary twist, the dominant ring motion, was estimated to show that the c-ring accounted in part for the measured compliance. Ring rotation was entrained to rotation of the external helix within each hairpin-shaped c-subunit in the ring. Ensembles of monomer and dimers extracted from complete c-rings showed that the coupling between collective ring and the individual subunit motions was independent of the size of the c-ring, which varies between organisms. Molecular determinants were identified by covariance analysis of residue coevolution and structural-alphabet-based local dynamics correlations. The residue coevolution gave a readout of subunit architecture. The dynamic couplings revealed that the hinge for both ring and subunit helix rotations was constructed from the proton-binding site and the adjacent glycine motif (IB-GGGG) in the midmembrane plane. IB-GGGG motifs were linked by long-range couplings across the ring, while intrasubunit couplings connected the motif to the conserved cytoplasmic loop and adjacent segments. The correlation with principal collective motions shows that the couplings underlie both ring rotary and bending motions. Noncontact couplings between IB-GGGG motifs matched the coevolution signal as well as contact couplings. The residue coevolution reflects the physiological importance of the dynamics that may link proton transfer to ring compliance.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Rigorous Representation of Transporters as Allosteric Machines Enables a Quantitative Understanding of their Functional Mechanisms

The sodium-coupled transporters are prototypical molecular machines of the cell that maintain communication with the environment by moving substances across the membrane. How they use energy to perform the transfer of substrates against their chemical gradient, i.e., by sodium symport along its gradient, is of special biophysical interest. New structural information, and powerful computational methods applied to the analysis of extensive molecular dynamics trajectories, make it possible to incorporate mechanistic elements of ion and substrate transport in structure-informed allosteric models that can be interrogated for specific and quantitative inferences.

Biophysics↗

DOE OSTI2025

Prediction of α $IIb$ $β$ 3 integrin structures along its minimum free energy activation pathway

The adhesion protein integrin is a transmembrane heterodimer that plays a pivotal role in cellular processes such as cell signaling and cell migration. To execute its function, integrin undergoes extensive conformational changes from a bent-closed to an extended-open state. Resolving the structures across these changes remains a challenge with both experimental and computational methods, but it is crucial for understanding the activation mechanism of integrin. We address this challenge for the platelet integrin α IIb β 3 by employing finite temperature string method with structures of the images along the initial guess path generated by a multiscale data-driven framework. The full-length all-atom structures along the resulting minimum free energy path between the inactive bent-closed and active extended-open states of α IIb β 3 integrin are consistent with a variety of experimentally resolved structures. Changes in these predicted structures along the path show that the extension and separation of the α and β subunits from the bent-closed to the extended-open state require correlated movements between the subdomain pairs in α IIb β 3 . Furthermore, these results provide new insights into integrin activation mechanism, and the predicted structures have potential applications in guiding the design of integrin-targeting therapeutics.

Dasetty, Siva [University of Chicago, IL (United S↗

DOE OSTI2023

Revisiting the SGTE lattice stability of bcc aluminum

The CALPHAD framework rests on the ability to assign well-defined free energies to all phases, including mechanically unstable ones. Given the experimental difficulties associated with the latter, researchers have turned to computational methods to obtain a theoretically justified definition of the free energy of mechanically unstable phases. Recently, a method called “inflection detection” has been shown to yield free energies that are broadly consistent with the current SGTE (Scientific Group Thermodata Europe) standard for pure elements, which were obtained through an extrapolation of stable solid solutions free energies towards unstable end members. Despite this general agreement, apparent failures for a few elements stand out, notably for the bcc–fcc enthalpy difference for the common element Al. Furthermore we revisit the definition of the bcc Al enthalpy and show that a careful CALPHAD-based analysis of enthalpy trends in five different mechanically stable binary solid solutions in fact does agree very well with the inflection-detection approach.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Structure and spectroscopy of graphite monofluoride

The structure of graphite monofluoride, (CF) n , has been debated since its discovery in 1934. In this work, we investigate a commercial graphite monofluoride by vibrational spectroscopy (infrared, Raman and the first inelastic neutron scattering spectra of this material). The spectroscopy shows that the material contains unreacted graphite and the partially fluorinated product dicarbon fluoride, (C 2 F) n , We evaluate the previously proposed $P\bar{6}m2$ and $P\bar{3}m1$ structures using computational methods and find F···F contacts render the $P\bar{6}m2$ structure dynamically unstable. We propose two alternative structures, $Cmc2_1$ $P6_3mc$, generated by displacement of one layer relative to another and find that $Cmc2_1$ is also dynamically unstable. The calculations are validated by comparison of calculated and observed INS spectra

Density functional theory↗