Search NASA⌕ Search

SEARCH · Search NASA

Results for “Structural simulation”

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

Investigation of the Effect of Temperature on the Structure of SARS-CoV-2 Spike Protein by Molecular Dynamics Simulations

Statistical and epidemiological data imply temperature sensitivity of the SARS-CoV2 coronavirus. However, the molecular level understanding of the virus structure at different temperature is still not clear. Spike protein is the outermost structural protein of the SARS-CoV-2 virus which interacts with the Angiotensin Converting Enzyme 2 (ACE2), a human receptor, and enters the respiratory system. In this study, we performed an all atom molecular dynamics simulation to study the effect of temperature on the structure of the Spike protein. After 200 ns of simulation at different temperatures, we came across some interesting phenomena exhibited by the protein. We found that the solvent exposed domain of Spike protein, namely S1, is more mobile than the transmembrane domain, S2. Structural studies implied the presence of several charged residues on the surface of N-terminal Domain of S1 which are optimally oriented at 10– 30°C. Bioinformatics analyses indicated that it is capable of binding to other human receptors and should not be disregarded. Additionally, we found that receptor binding motif (RBM), present on the receptor binding domain (RBD) of S1, begins to close around temperature of 40°C and attains a completely closed conformation at 50°C. We also found that the presence of glycan moieties did not influence the observed protein dynamics. Nevertheless, the closed conformation disables its ability to bind to ACE2, due to the burying of its receptor binding residues. Our results clearly show that there are active and inactive states of the protein at different temperatures. This would not only prove beneficial for understanding the fundamental nature of the virus, but would be also useful in the development of vaccines and therapeutics.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Smoothened stimulation by membrane sterols drives Hedgehog pathway activity

Hedgehog signalling is fundamental to embryonic development and postnatal tissue regeneration. Aberrant postnatal Hedgehog signalling leads to several malignancies, including basal cell carcinoma and paediatric medulloblastoma. Hedgehog proteins bind to and inhibit the transmembrane cholesterol transporter Patched-1 (PTCH1), which permits activation of the seven-transmembrane transducer Smoothened (SMO) via a mechanism that is poorly understood. Here we report the crystal structure of active mouse SMO bound to both the agonist SAG21k and to an intracellular binding nanobody that stabilizes a physiologically relevant active state. Analogous to other G protein-coupled receptors, the activation of SMO is associated with subtle motions in the extracellular domain, and larger intracellular changes. In contrast to recent models, a cholesterol molecule that is critical for SMO activation is bound deep within the seven-transmembrane pocket. We propose that the inactivation of PTCH1 by Hedgehog allows a transmembrane sterol to access this seven-transmembrane site (potentially through a hydrophobic tunnel), which drives the activation of SMO. Overall, these results—combined with signalling studies and molecular dynamics simulations—delineate the structural basis for PTCH1–SMO regulation, and suggest a strategy for overcoming clinical resistance to SMO inhibitors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Elucidating the formation mechanisms of zeolites using data-driven modelling and in-situ characterization

Zeolites are the main solid catalysts used by the chemical industry. The use of zeolites in separations and as shape selective catalysts requires control of the width and connectivity of their pores. 235 distinct zeolite frameworks have been synthesized to date, of over 2 million that have been proposed. Recent work indicates that the limitation is in large part kinetic: new synthetic pathways are required to access new zeolites. Organic cations are used to direct the synthesis towards specific zeolites. However, the molecular mechanisms by which cations direct the nucleation towards specific zeolites is not known. Elucidating these mechanisms is key to realize new zeolites for catalysis and separations, and is the focus of this project. This project developed and implemented a synergistic, data-driven computational and experimental approach to resolve the molecular pathways of nucleation, growth, and polymorph selection of zeolites and the role of organic cations in directing their formation. The project developed computationally efficient and accurate models for the study of the nucleation and growth of pure silica zeolites in molecular simulations, using machine learning with data from experiments. Simulations with these models were integrated with scanning tunneling electron microscopy, computer vision, and deep learning to unveil the molecular pathways of formation of a zeolite. Of particular interest in this project was to elucidate the role of amorphous precursors in the nucleation of the zeolite. Previous experiments indicate that zeolites are born within non-crystalline aggregates in which the silicates and organic cations have local and medium range order similar to that of the zeolite. The organic cations that direct the formation of zeolites and those that direct the formation of ordered mesoporous silicas are similar. We hypothesized that the frustrated attraction that for large organic cations leads to the formation of stable mesophases that direct the synthesis of mesoporous silicas, could promote the formation of metastable mesophases that can assist in the nucleation and polymorph selection of zeolites. The simulations resolved how structure directing agents build crystalline order and showed that mesoscopic pre-ordering occurs is not required to facilitate the nucleation of zeolites, because the synthesis occurs at high driving forces, where the barriers for nucleation are negligible. This project unveiled that polymorph selection in zeolite synthesis occurs after nucleation, opening a distinct area of control through the kinetics of growth and not through nucleation barriers.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An Open-Source Parallel EMT Simulation Framework

As the integration level of inverter-based resources (IBRs) increases, ensuring the reliable operation of the bulk power systems requires the use of electromagnetic transient (EMT) simulation tools to identify and mitigate system-wide stability risks. Conducting EMT studies for large-scale, IBR-rich grids, however, is challenging due to the inherent computational bottleneck caused by the underlying high-fidelity models and required small time steps. This paper introduces ParaEMT: an open-source, generic EMT simulation framework designed to accelerate simulations by leveraging advanced parallel computational technologies, such as high-performance computers. This paper presents a comprehensive exposition of ParaEMT, covering its modeling library, simulation strategy, framework structure, operational procedures, and auxiliary features, alongside its extensible parallel computational architecture. Notably, ParaEMT is a publicly accessible and modularized framework written in Python, thereby facilitating future development and the integration of new models and algorithms. The accuracy and efficiency of ParaEMT are demonstrated by rigorous validations via multiple case studies.

electromagnetic transient simulation↗
DOE OSTI2023

An Open-Source Parallel EMT Simulation Framework: Preprint

As the integration level of inverter-based resources (IBR) increases, ensuring the reliable operation of the bulk power systems requires the use of electromagnetic transient (EMT) simulation tools to identify and mitigate system-wide stability risks. Conducting EMT studies for large-scale, IBR-rich grids, however, is challenging due to the inherent computational bottleneck caused by the underlying high-fidelity models and required small time steps. This paper introduces ParaEMT: an open-source, generic EMT simulation framework designed to accelerate simulations by leveraging advanced parallel computational technologies, such as high-performance computers. This paper presents a comprehensive exposition of ParaEMT, covering its modeling library, simulation strategy, framework structure, operational procedures, and auxiliary features, alongside its extensible parallel computational architecture. Notably, ParaEMT is a publicly accessible and modularized framework written in Python, thereby facilitating future development and the integration of new models and algorithms. The accuracy and efficiency of ParaEMT are demonstrated by rigorous validations via multiple case studies.

electromagnetic transient simulation↗
DOE OSTI2025

Numerical Simulation Studies of Ultrasonic De-Icing for Heating, Ventilation, Air Conditioning, and Refrigeration Structures

Ice accumulation on heating, ventilation, air conditioning, and refrigeration (HVACR) structures presents significant operational challenges. These challenges include reduced efficiency, increased energy consumption, and potential damage to equipment. Traditional de-icing methods, such as chemical treatments, mechanical scraping, or heating-based techniques, are often labor-intensive, costly, and environmentally harmful. Here, this study uniquely investigates ultrasonic de-icing as an energy-efficient alternative for HVACR applications, focusing on the specific structural geometries found in these systems. A comprehensive numerical simulation framework was developed using finite element analysis to explore ultrasonic wave propagation across four distinct HVACR structures. Key parameters such as ultrasonic frequency, power levels, and the number and placement of actuators were examined for their impact on ice detachment efficiency. Results from simulations on a plate structure reveal that ultrasonic excitation can propagate effectively across large areas (at least 150 × 150 mm), enhancing the de-icing coverage. Lower frequency (e.g., 30 to 45 kHz) excitation results in greater displacement, improving de-icing performance, while increased actuator numbers with the same total power input also enhance effectiveness. Two actuators seem sufficient for the de-icing of a 300 × 300 mm plate. For tube-and-fin structures, specific high-power ultrasonic frequencies selectively excite the fin plates, demonstrating efficient ice removal when actuated on the tube. However, optimal performance requires careful design of actuator placement and vibration modes to accommodate the irregular shapes of these structures.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

Acetylcholinesterase: Structure, dynamics, and interactions with organophosphorus compounds

Acetylcholinesterase (AChE) is an enzyme that hydrolyzes the neurotransmitter acetylcholine (ACh), removing it from the synaptic cleft after the transmission of an electrical signal, making it an essential component of chemical neurotransmission. AChE is a serine hydrolase, containing a catalytic triad of Ser/His/Glu. AChE is a prime target for pharmaceuticals treating a variety of neurological disorders. It is also the target of synthetic organophosphorus (OP) compounds that have been used as pesticides and chemical warfare agents. OP compounds contain a potent leaving group, such as fluorine, and act by forming a covalent adduct with the catalytic serine of the AChE active site. A wealth of structural information is available for AChE, including over 300 structures, including a subset of structures in complex with drugs as well as OP compounds. This review will highlight the interactions between OP compounds and AChE from a structural and computational perspective, with a discussion of access to the active site, as well as side reactions that lead to dealkylation of the OP-catalytic serine adduct, a process known as aging. We conclude that while the majority of the conformational changes needed to accommodate the OP compounds are localized to the acyl loop in the crystal structures, molecular dynamics simulations highlight the potential for a far more dynamic enzyme.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Structure of the Bastnäsite (001) Surface by Crystal Truncation Rod X-ray Diffraction and Ab Initio Molecular Dynamics: Implications for Separations of a Rare Earth Ore Mineral

Bastnäsite ((Ce,La)FCO 3 ) is the primary mineral source of light rare earth elements, but its surface structure is not well understood. This presents a major challenge in improving beneficiation strategies. In this work, a synergistic combination of X-ray scattering and ab initio molecular dynamics (AIMD) was used to gain atomistic insight into the interfacial structure of bastnäsite. Surface X-ray scattering was used to measure crystal truncation rods (CTRs) of the bastnäsite (001) surface, a significant crystal face with a previously unknown termination. The best-fit atomic-scale model of the CTR data features a carbonate layer at the surface, which is stabilized by the relaxation of carbonate groups from their bulk structural positions. AIMD simulations predict similar surface relaxations, which are shown to be influenced by the protonation of oxygen atoms at the surface. Evidence of ordered water at the interface is also observed in the best-fit model and AIMD simulations. Further, the presence of a carbonate layer at this dominant crystal surface is significant for improving separation technologies because most commonly used ligands utilize anionic functional groups to chelate metal cations at particle surfaces. Without modification, anionic ligands are expected to have poor affinity for the carbonate-terminated (001) surface.

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

Theoretical and experimental investigations on H2 sensing properties of flower-like titanium dioxide

Highlights: • Flower-like rutile TiO{sub 2} has been successfully synthesized. • The unique nanostructure show an great enhanced gas-sensing performance. • A sensing mechanism based on our calculation was proposal. - Abstract: Microstructure and gas sensing properties of partly agglomerating flower-like rutile TiO{sub 2} obtained using a one-pot hydrothermal method was investigated. Test results exhibited an excellent sensitivity to H{sub 2} gas and stable reproducibility. Further, we performed a first principles calculation of an established surface model to delve into the sensing mechanism. Adsorbed oxygen preferred to locate upon the Ti{sub 5c} site, while the electronic structure of surface had rather weak changes after hydrogen adsorption. This simulation may pave a new path to further discover the sensing mechanism of rutile TiO{sub 2}-based gas sensor.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Mol* Viewer: modern web app for 3D visualization and analysis of large biomolecular structures

Abstract Large biomolecular structures are being determined experimentally on a daily basis using established techniques such as crystallography and electron microscopy. In addition, emerging integrative or hybrid methods (I/HM) are producing structural models of huge macromolecular machines and assemblies, sometimes containing 100s of millions of non-hydrogen atoms. The performance requirements for visualization and analysis tools delivering these data are increasing rapidly. Significant progress in developing online, web-native three-dimensional (3D) visualization tools was previously accomplished with the introduction of the LiteMol suite and NGL Viewers. Thereafter, Mol* development was jointly initiated by PDBe and RCSB PDB to combine and build on the strengths of LiteMol (developed by PDBe) and NGL (developed by RCSB PDB). The web-native Mol* Viewer enables 3D visualization and streaming of macromolecular coordinate and experimental data, together with capabilities for displaying structure quality, functional, or biological context annotations. High-performance graphics and data management allows users to simultaneously visualise up to hundreds of (superimposed) protein structures, stream molecular dynamics simulation trajectories, render cell-level models, or display huge I/HM structures. It is the primary 3D structure viewer used by PDBe and RCSB PDB. It can be easily integrated into third-party services. Mol* Viewer is open source and freely available at https://molstar.org/.

Sehnal, David↗
DOE OSTI2023

Exploring the Performance Boundaries of a Small Reconfigurable Multi-Mission UAV through Multidisciplinary Analysis

The performance of a small reconfigurable unmanned aerial vehicle (UAV) is evaluated, combining a multidisciplinary approach in the computational analysis of additive manufactured structures, fluid dynamics, and experiments. Reconfigurable UAVs promise cost savings and efficiency, without sacrificing performance, while demonstrating versatility to fulfill different mission profiles. The use of computational fluid dynamics (CFD) in UAV design produces higher accuracy aerodynamic data, which is particularly important for complex aircraft concepts such as blended wing bodies. To address challenges relating to anisotropic materials, the Tsai–Wu failure criterion is applied to the structural analysis, using CFD solutions as load inputs. Aerodynamic performance results show the low-speed variant attains an endurance of 1 h, 48 min, whereas its high-speed counterpart is 29 min at a 66.7% higher cruise speed. Each variant serves different aspects of small UAS deployment, with low speed envisioned for high-endurance surveying, and high speed for long-range or time-critical missions such as delivery. The experimental and simulation results suggest room for design iteration, in wing area and geometry adjustments. Structural simulations demonstrated the need for airframe improvements to the low-speed configuration. This paper highlights the potential of reconfigurable UAVs to be useful across multiple industries, advocating for further research and design improvements.

42 ENGINEERING↗
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 OSTI2022

Using Machine Learning to Greatly Accelerate Path Integral Ab Initio Molecular Dynamics

Ab initio molecular dynamics (AIMD) has become one of the most popular and robust approaches for modeling complicated chemical, liquid, and material systems. However, the formidable computational cost often limits its widespread application in simulations of the largest-scale systems. Furthermore, the situation becomes even more severe in cases where the hydrogen nuclei may be better described as quantized particles using a path integral representation. Here, we present a computational approach that combines machine learning with recent advances in path integral contraction schemes, and we achieve a 2 orders of magnitude acceleration over direct path integral AIMD simulation while at the same time maintaining its accuracy.

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

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Critical Factors in Computational Characterization of Hydrogen Storage in Metal–Organic Frameworks

Inconsistencies in high-pressure H2 adsorption data and a lack of comparative experiment–theory studies have made the evaluation of both new and existing metal–organic frameworks (MOFs) challenging in the context of hydrogen storage applications. In this work, we performed grand canonical Monte Carlo (GCMC) simulations in nearly 500 experimentally refined MOF structures to examine the variance in simulation results because of the equation of state, H 2 potential, and the effect of density functional theory structural optimization. We find that hydrogen capacity at 77 K and 100 bar, as well as hydrogen 100-to-5 bar deliverable capacity, is correlated more strongly with the MOF pore volume than with the MOF surface area (the latter correlation is known as the Chahine’s rule). The tested methodologies provide consistent rankings of materials. In addition, four prototypical MOFs (MOF-74, CuBTC, ZIF-8, and MOF-5) with a range of surface areas, pore structures, and surface chemistries, representative of promising adsorbents for hydrogen storage, are evaluated in detail with both GCMC simulations and experimental measurements. Simulations with a three-site classical potential for H 2 agree best with our experimental data except in the case of MOF-5, in which H 2 adsorption is best replicated with a five-site potential. However, for the purpose of ranking materials, these two choices for H 2 potential make little difference. Here more significantly, 100 bar loading estimates based on more accurate equations of state for the vapor–liquid equilibrium yield the best comparisons with the experiment.

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

Lattice Structured Lightweight Structural Materials

The development of lightweight structural materials is crucial for enhancing the performance and deployment feasibility of fission batteries. This study aims to produce lightweight structural materials whose strength-to-weight ratios exceed those of current widely used structural materials. To achieve this, advanced modeling and simulation tools were employed to design lattice structures with different lattice parameters and different lattice types. A process was successfully developed for transforming lattice-structured models into Multiphysics Object Oriented Simulation Environment (MOOSE) inputs. Finite element modeling (FEM) was used to simulate the uniaxial tensile testing of the lattice-structured parts to investigate the stress distribution at a given displacement. The modeling results showed that the lattice-structured sample displayed a lower Young’s modulus in comparison to the solid material; the increase in solid shell thickness and blend radius enhances the mechanical performance; and the effect of unit cell size on macro scale stress is minimal. Tensile testing was conducted on the solid and lattice-structured materials fabricated by laser powder bed fusion (LPBF) additive manufacturing. The experimental results agreed well with the model prediction. The approach of using modeling as a guiding tool for preliminary material design can significantly save time and cost for new material development.

lightweight material↗
DOE OSTI2023

Structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquids close to Earth's liquid core conditions

Using an artificial neural-network machine learning interatomic potential, we have performed molecular dynamics simulations to study the structure and dynamics of $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid close to the Earth's liquid core conditions. Here, the simulation results reveal that the short-range structural order (SRO) in the $\mathrm{Fe_{90}}$$\mathrm{Si_3}$$\mathrm{O_7}$ liquid is very strong. About 80% of the atoms are arranged in crystallinelike SRO motifs. In particular, ~ $70$% of Fe-centered clusters can be classified as either hexagonal-close-packed-like or icosahedrallike SRO motifs. The SRO clusters centered on Fe, Si, or O atoms are strongly intermixed and homogenously distributed throughout the liquid. The atomic structure of the liquid and the fractions of dominant SRO clusters are not sensitive to pressure/temperature used in the simulations except that the SRO of the O-centered clusters is enhanced close to inner core pressures. The O-diffusion coefficient is about two to three times larger than the Fe and Si ions and increases more rapidly in the deeper core regions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗