SEARCH · Search NASA

Results for “Cluster Dynamics Simulations”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Simulation of polymerization induced phase separation in model thermosets

Polymerization induced phase separation (PIPS) in a three component thermoset is studied using molecular dynamics simulations of a new coarse-grained thermoset model. The system includes two crosslinker molecules, which differ in their glass transition temperatures (T g ) and chain length and thus have the potential for phase separation. Here, one crosslinker has a high T g corresponding to a rubbery behavior, and simulations were performed for a short length (4 beads) and a long length (33 beads). The resin and other crosslinker have low T g . A coarse-grained model is developed with these features and with interaction parameters determined so that for either rubbery crosslinker length, the system is in the liquid state at the cure temperature. For sufficiently slow reaction rates, the long rubbery molecule exhibits PIPS into a bicontinuous array of nanoscale domains, but the short one does not, reproducing recent experimental results. The simulations demonstrate that the reaction rates must be slow enough to allow diffusion to yield phase separation. Particularly, the reaction rate corresponding to the secondary amine must be very slow, else the structure of crosslinked clusters and the substantially increased diffusion time will prevent PIPS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Molecular Design Principles for Photosystem I-Based Biohybrid Solar Fuel Catalysts

Direct solar-to-chemical conversion offers a compelling route to clean, dispatchable energy. Photosystem I (PSI), an evolutionarily optimized light-driven oxidoreductase, can be repurposed for solar-fuel production by coupling its photochemistry to catalytic interfaces. However, the molecular determinants that govern productive electron transfer to abiotic catalysts remain poorly understood. Here, we present molecular structures of active PSI-Pt nanoparticle (PtNP) biohybrids that reveal how protein architecture controls catalyst access, binding geometry, and photocatalytic efficiency. Removal of stromal subunits exposes the electron transfer chain and enables PtNP binding proximal to the F X cluster, demonstrating that steric occlusion limits access to native acceptor regions in PSI. In contrast, in trimeric PSI, PtNPs bind at multiple sites per monomer, but only a subset are positioned within electron transfer distance of terminal cofactors, resulting in a heterogeneous population of productive and nonproductive configurations. Structural analyses and molecular dynamics simulations define the interface topology, electrostatics, and cofactor-to-nanoparticle distances that govern catalyst binding and electron transfer. These results establish that catalytic inefficiency arises not only from intrinsic electron transfer constraints but also from the distribution of binding geometries imposed by the protein scaffold. Together, these findings provide a molecular framework linking protein structure to biohybrid function and define design principles for engineering PSI-based solar fuel systems and protein-nanomaterial interfaces for light-driven catalysis.

biohybrid↗

NASA NTRS1995

Towards Distributed Fluids/Controls Simulations

Computational prototyping of aircraft in nonlinear regimes will require the coupling of high-order disciplines on high performance computers. This effort investigates the coupling of the Reynolds-averaged Navier-Stokes equations with six degree-of-freedom nonlinear body dynamics and controls on distributed architectures. Towards this goal, solutions obtained on a wide-node IBM SP-2, an Ethernet-linked cluster of workstations, and a Cray C90 are compared in terms of accuracy and throughput.

Atwood, Christopher A.↗

DOE OSTI2020

Probing the nature of dark matter with accreted globular cluster streams

ABSTRACT The steepness of the central density profiles of dark matter (DM) in low-mass galaxy haloes (e.g. dwarf galaxies) is a powerful probe of the nature of DM. We propose a novel scheme to probe the inner profiles of galaxy subhaloes using stellar streams. We show that the present-day morphological and dynamical properties of accreted globular cluster (GC) streams – those produced from tidal stripping of GCs that initially evolved within satellite galaxies and later merged with the Milky Way (MW) – are sensitive to the central DM density profile and mass of their parent satellites. GCs that accrete within cuspy cold dark matter (CDM) subhaloes produce streams that are physically wider and dynamically hotter than streams that accrete inside cored subhaloes. A first comparison of MW streams ‘GD-1’ and ‘Jhelum’ (likely of accreted GC origin) with our simulations indicates a preference for cored subhaloes. If these results hold up in future data, the implication is that either the DM cusps were erased by baryonic feedback, or their subhaloes naturally possessed cored density profiles implying particle physics models beyond CDM. Moreover, accreted GC streams are highly structured and exhibit complex morphological features (e.g. parallel structures and ‘spurs’). This implies that the accretion scenario can naturally explain the recently observed peculiarities in some of the MW streams. We also propose a novel mechanism for forming ‘gaps’ in stellar streams when the remnant of the parent subhalo (which hosted the GC) later passes through the GC stream. This encounter can last a longer time (and have more of an impact) than the random encounters with DM subhaloes previously considered, because the GC stream and its parent subhalo are on similar orbits with small relative velocities. Current and future surveys of the MW halo will uncover numerous faint stellar streams and provide the data needed to substantiate our preliminary tests with this new probe of DM.

(Galaxy:) globular clusters: individual↗

DOE OSTI2023

Ab Initio Prediction of Ultra‐Wide Band Gap B x Al 1− x N Materials

Abstract Ultra‐wide bandgap (UWBG) materials are poised to play an important role in the future of power electronics. Devices made from UWBG materials are expected to operate at higher voltages, frequencies, and temperatures than current silicon and silicon‐carbide‐based devices and can even lead to significant miniaturization of such devices. In the UWBG field, aluminum nitride and boron nitride have attracted great interest; however, the B x Al 1− x N alloys are much less studied. In this work, using first‐principles simulations combining density‐functional theory and the cluster expansion method, the crystal structure of B x Al 1− x N alloys is predicted. Seventeen ground state structures of B x Al 1− x N with formation energies between 0.11 and 0.25 eV atom ‐1 are found. All of these structures are found to be dynamically stable. The B x Al 1− x N structures are found to have predominantly a tetrahedral bonding environment; however, some structures exhibit sp 2 bonds similar to hexagonal BN. This work expands knowledge of the structures, energies, and bonding in B x Al 1− x N which aids their synthesis, the innovation of lateral or vertical devices, and discovery of compatible dielectric and Ohmic contact materials.

36 MATERIALS SCIENCE↗

DOE OSTI2024

A Generative Model for Realistic Galaxy Cluster X-Ray Morphologies

Abstract The X-ray morphologies of clusters of galaxies display significant variations, reflecting their dynamical histories and the nonlinear dependence of X-ray emissivity on the density of the intracluster gas. Qualitative and quantitative assessments of X-ray morphology have long been considered a proxy for determining whether clusters are dynamically active or “relaxed.” Conversely, the use of circularly or elliptically symmetric models for cluster emission can be complicated by the variety of complex features realized in nature, spanning scales from megaparsecs down to the resolution limit of current X-ray observatories. In this work, we use mock X-ray images from simulated clusters from The Three Hundred project to define a basis set of cluster image features. We take advantage of the clusters’ approximate self-similarity to minimize the differences between images before encoding the remaining diversity through a distribution of high-order polynomial coefficients. Principal component analysis then provides an orthogonal basis for this distribution, corresponding to natural perturbations from an average model. This representation allows novel, realistically complex X-ray cluster images to be easily generated, and we provide code to do so. The approach provides a simple way to generate training data for cluster image analysis algorithms and could be straightforwardly adapted to generate clusters displaying specific types of features or selected by physical characteristics available in the original simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2007

Computational Plume Modeling of COnceptual ARES Vehicle Stage Tests

The plume-induced environment of a conceptual ARES V vehicle stage test at the NASA Stennis Space Center (NASA-SSC) was modeled using computational fluid dynamics (CFD). A full-scale multi-element grid was generated for the NASA-SSC B-2 test stand with the ARES V stage being located in a proposed off-center forward position. The plume produced by the ARES V main power plant (cluster of five RS-68 LOX/LH2 engines) was simulated using a multi-element flow solver - CRUNCH. The primary objective of this work was to obtain a fundamental understanding of the ARES V plume and its impingement characteristics on the B-2 flame-deflector. The location, size and shape of the impingement region were quantified along with the un-cooled deflector wall pressures, temperatures and incident heating rates. Issues with the proposed tests were identified and several of these addressed using the CFD methodology. The final results of this modeling effort will provide useful data and boundary conditions in upcoming engineering studies that are directed towards determining the required facility modifications for ensuring safe and reliable stage testing in support of the Constellation Program.

Allgood, Daniel C.↗

DOE OSTI2024

Quantum fluctuations drive nonmonotonic correlations in a qubit lattice

Abstract Fluctuations may induce the degradation of order by overcoming ordering interactions, consequently leading to an increase of entropy. This is particularly evident in magnetic systems characterized by nontrivial, constrained disorder, where thermal or quantum fluctuations can yield counterintuitive forms of ordering. Using the proven efficiency of quantum annealers as programmable spin system simulators, we present a study based on entropy postulates and experiments on a platform of programmable superconducting qubits to show that a low level of uncertainty can promote ordering in a system impacted by both thermal and quantum fluctuations. A set of experiments is proposed on a lattice of interacting qubits arranged in a triangular geometry with precisely controlled disorder, effective temperature, and quantum fluctuations. Our results demonstrate the creation of ordered ferrimagnetic and layered anisotropic disordered phases, displaying characteristics akin to the elegant order-by-disorder phenomenon. Extensive experimental evidence is provided for the role of quantum fluctuations in lowering the total energy of the system by increasing entropy and defect clustering. Our thorough and comprehensive application of an intentionally introduced noise on a quantum platform provides insight into the dynamics of defects and fluctuations in quantum devices, which may help to reduce the cost associated with quantum processing.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1994

Computational strategies for three-dimensional flow simulations on distributed computer systems

An increasing amount of research activity in computational fluid dynamics has been devoted to the development of efficient algorithms for parallel computing systems. The increasing performance to price ratio of engineering workstations has led to research to development procedures for implementing a parallel computing system composed of distributed workstations. This thesis proposal outlines an ongoing research program to develop efficient strategies for performing three-dimensional flow analysis on distributed computing systems. The PVM parallel programming interface was used to modify an existing three-dimensional flow solver, the TEAM code developed by Lockheed for the Air Force, to function as a parallel flow solver on clusters of workstations. Steady flow solutions were generated for three different wing and body geometries to validate the code and evaluate code performance. The proposed research will extend the parallel code development to determine the most efficient strategies for unsteady flow simulations.

Weed, Richard Allen↗

NASA NTRS1995

Computational strategies for three-dimensional flow simulations on distributed computer systems

This research effort is directed towards an examination of issues involved in porting large computational fluid dynamics codes in use within the industry to a distributed computing environment. This effort addresses strategies for implementing the distributed computing in a device independent fashion and load balancing. A flow solver called TEAM presently in use at Lockheed Aeronautical Systems Company was acquired to start this effort. The following tasks were completed: (1) The TEAM code was ported to a number of distributed computing platforms including a cluster of HP workstations located in the School of Aerospace Engineering at Georgia Tech; a cluster of DEC Alpha Workstations in the Graphics visualization lab located at Georgia Tech; a cluster of SGI workstations located at NASA Ames Research Center; and an IBM SP-2 system located at NASA ARC. (2) A number of communication strategies were implemented. Specifically, the manager-worker strategy and the worker-worker strategy were tested. (3) A variety of load balancing strategies were investigated. Specifically, the static load balancing, task queue balancing and the Crutchfield algorithm were coded and evaluated. (4) The classical explicit Runge-Kutta scheme in the TEAM solver was replaced with an LU implicit scheme. And (5) the implicit TEAM-PVM solver was extensively validated through studies of unsteady transonic flow over an F-5 wing, undergoing combined bending and torsional motion. These investigations are documented in extensive detail in the dissertation, 'Computational Strategies for Three-Dimensional Flow Simulations on Distributed Computing Systems', enclosed as an appendix.

Sankar, Lakshmi N.↗

DOE OSTI2020

Control of spin dynamics in artificial honeycomb spin-ice-based nanodisks

In this paper, we report the experimental and theoretical characterization of the angular-dependent spin dynamics in arrays of ferromagnetic nanodisks arranged on a honeycomb lattice. The magnetic field and microwave frequency dependence, measured by broadband ferromagnetic resonance, reveal a rich spectrum of modes that is strongly affected by the microstate of the network. Based on symmetry arguments with respect to the external field, we show that certain parts of the ferromagnetic network contribute to the detected signal. A comparison of the experimental data with micromagnetic simulations reveals that different subsections of the lattice predominantly contribute to the high-frequency response of the array. This is confirmed by optical characterizations using microfocused Brillouin light scattering. Furthermore, we find indications that nucleation and annihilation of vortexlike magnetization configurations in the low-field range affect the dynamics, which is different from clusters of ferromagnetic nanoellipses. Our work opens up new perspectives for designing magnonic devices that combine geometric frustration in gyrotropic vortex crystals at low frequencies with magnonic crystals at high frequencies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1990

Traffic placement policies for a multi-band network

Recently protocols were introduced that enable the integration of synchronous traffic (voice or video) and asynchronous traffic (data) and extend the size of local area networks without loss in speed or capacity. One of these is DRAMA, a multiband protocol based on broadband technology. It provides dynamic allocation of bandwidth among clusters of nodes in the total network. A number of traffic placement policies for such networks are proposed and evaluated. Metrics used for performance evaluation include average network access delay, degree of fairness of access among the nodes, and network throughput. The feasibility of the DRAMA protocol is established through simulation studies. DRAMA provides effective integration of synchronous and asychronous traffic due to its ability to separate traffic types. Under the suggested traffic placement policies, the DRAMA protocol is shown to handle diverse loads, mixes of traffic types, and numbers of nodes, as well as modifications to the network structure and momentary traffic overloads.

Maly, Kurt J.↗

NASA NTRS1991

Multiple-Bit Errors Caused By Single Ions

Report describes experimental and computer-simulation study of multiple-bit errors caused by impingement of single energetic ions on 256-Kb dynamic random-access memory (DRAM) integrated circuit. Studies illustrate effects of different mechanisms for transport of charge from ion tracks to various elements of integrated circuits. Shows multiple-bit errors occur in two different types of clusters about ion tracks causing them.

Zoutendyk, John A.↗

DOE OSTI2026

Protonation Dynamics of Confined Ethanol–Water Mixtures in H-ZSM-5 from Machine Learning-Driven Metadynamics

Zeolites are indispensable heterogeneous catalysts in industrial chemical processes, valued for their strong Brønsted acidity, well-defined microporous frameworks, and tunable pore structures. Their catalytic activity arises primarily from Brønsted acid sites (BAS), typically present as bridging hydroxyl groups (Si–OH–Al). Under aqueous reaction conditions, these protons interact dynamically with water and alcohol molecules, leading to complex solvation and protonation behavior within confined pores. In this study, we investigate the protonation equilibrium occurring between ethanol and water at the BAS of acidic zeolites under varying hydration levels, i.e., C2H5OH–(H2O)n, n=1–4. Local structure was analyzed through an adaptive-learning global optimization algorithm, while enhanced sampling molecular dynamics simulations with Well-Tempered Metadynamics (WMetaD) and machine learning interatomic potentials (MLPs) provide free-energy surfaces (FES) at variable hydration levels. The results reveal a strong dependence of proton localization on the degree of hydration. At low hydration (1 water molecule), the proton resides predominantly on ethanol; with 2 water molecules, it shifts toward water, and at higher hydration (3 or more water molecules), it becomes extensively delocalized over the water cluster. These findings underscore the critical role of solvation in modulating acid site behavior and suggest that a minimum of three water molecules is necessary to fully stabilize the proton on water within the zeolite framework. This solvation threshold has significant implications for catalytic processes, particularly in biomass conversion reactions where alcohol protonation is a key step in dehydration mechanisms.

machine learning↗

DOE OSTI2021

Molecular dynamics simulations of radiation cascade evolution near cellular dislocation structures in additively manufactured stainless steels

The dislocation dense cellular structures observed in additively manufactured (AM) stainless steels can allow these materials to have enhanced strength and ductility compared to conventionally manufactured materials. Effective design of radiation resistant materials often necessitates creating microstructures containing of a high density of dislocations that act as sinks for radiation induced defects. This work uses molecular dynamics simulations to study the impact of pre-existing defects on radiation damage in stainless steel 316L fabricated by the laser powder bed fusion process. The evolution of synthetically generated dislocations in a crystal system in response to multiple radiation collision cascades suggest that dislocation dense regions reduce the probability of surviving point defects forming clusters. This finding seems to result from the dislocation cores’ ability to absorb point defects, thereby limiting the formation of clusters during the recombination phase. However, after successive cascades, the ability of the dislocation entanglement to prevent defect formation is diminished and defect cluster formation trends to an equilibrium. Accordingly, it is hypothesized that in an experimental setting the dislocation cells will initially act as neutral sinks for point defects, potentially delaying the onset of radiation damage effects; but, these structures will be degraded due to radiation-enhanced diffusion and lose their effectiveness at doses above 1.5–2 displacements per atom. Nevertheless, the ability to produce complex and feature-specific microstructures suggests promise for the design of radiation tolerant materials using AM methods.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Metastable–solid phase diagrams derived from polymorphic solidification kinetics

Nonequilibrium processes during solidification can lead to kinetic stabilization of metastable crystal phases. A general framework for predicting the solidification conditions that lead to metastable-phase growth is developed and applied to a model face-centered cubic (fcc) metal that undergoes phase transitions to the body-centered cubic (bcc) as well as the hexagonal close-packed phases at high temperatures and pressures. Large-scale molecular dynamics simulations of ultrarapid freezing show that bcc nucleates and grows well outside of the region of its thermodynamic stability. An extensive study of crystal–liquid equilibria confirms that at any given pressure, there is a multitude of metastable solid phases that can coexist with the liquid phase. We define for every crystal phase, a solid cluster in liquid (SCL) basin, which contains all solid clusters of that phase coexisting with the liquid. Additionally, a rigorous methodology is developed that allows for practical calculations of nucleation rates into arbitrary SCL basins from the undercooled melt. It is demonstrated that at large undercoolings, phase selections made during the nucleation stage can be undone by kinetic instabilities amid the growth stage. On these bases, a solidification–kinetic phase diagram is drawn for the model fcc system that delimits the conditions for macroscopic grains of metastable bcc phase to grow from the melt. We conclude with a study of unconventional interfacial kinetics at special interfaces, which can bring about heterogeneous multiphase crystal growth. A first-order interfacial phase transformation accompanied by a growth-mode transition is examined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Spatio-Temporal Clustering of Seismicity Enabled by Off-Fault Plasticity

While significant progress has been made in understanding earthquake source processes in linear elastic domains, the effect of more realistic rheologies including plasticity is poorly understood. Here, in this work, we simulate the sequence of earthquake and aseismic slip of a 2D antiplane rate-and-state fault embedded in a full-space elastic-plastic bulk. We show that off-fault plasticity may lead to partial ruptures as well as temporal clustering of seismic events. Furthermore, the interaction of fault slip and off-fault plasticity results in pockets of slip deficit. While the energy dissipated through plastic deformation remains a small fraction of the total energy budget, its impact on the source characteristics is disproportionally large through the redistribution of stresses and viscous relaxation. Our results suggest a new mechanism of dynamic heterogeneity in earthquake physics that may have important implications on earthquake size distribution and energy budget.

58 GEOSCIENCES↗

DOE OSTI2026

Dynamics of iodine geminate recombination in supercritical xenon solvent: Caging effect

Understanding the dynamics of chemical reactions in solutions is vital, as their rates and kinetics are significantly affected by the solvent environment. Supercritical solvents offer extensive applications in chemical reactions by enabling the manipulation of the solution environment. Here, in this study, we investigate the geminate recombination of iodine in a supercritical xenon solvent by using ReaxFF-based molecular dynamics simulations. Our findings reveal that the highest iodine recombination rate occurs near supercritical conditions, while lower-pressure conditions lead to reduced collision rates and unstable recombination, and higher-pressure conditions hinder iodine diffusion, resulting in a lower recombination rate. Our analysis shows that the xenon local density at the time of recombination is at least 2.5 times higher than the global density, confirming the presence of xenon clusters surrounding the Iodine atoms. This observation is further supported by coordination number analysis, which confirms an elevated xenon local density during recombination. In addition, the correlation between the total energy of xenon atoms within a cluster and recombined iodine atoms underscores the kinetic energy transfer process, validating the occurrence of geminate recombination. The excess kinetic energy from the recombining iodine atoms is transferred to the surrounding xenon atoms. Our examination of geminate recombination demonstrates that iodine atoms confined within xenon clusters—whether through manual insertion of atoms or the fast dissociation of an iodine molecule within xenon clusters—are more likely to recombine as primary geminate recombination. However, extending the iodine molecule dissociation time allows iodine atoms to diffuse out of the cluster, and the recombination to shift toward secondary geminate recombination.

Cage effect↗