SEARCH · Search NASA

Results for “Coarse grained 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 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗

DOE OSTI2024

CoarsenConf: Equivariant Coarsening with Aggregated Attention for Molecular Conformer Generation

Molecular conformer generation (MCG) is an important task in cheminformatics and drug discovery. The ability to efficiently generate low-energy 3D structures can avoid expensive quantum mechanical simulations, leading to accelerated virtual screenings and enhanced structural exploration. Several generative models have been developed for MCG, but many struggle to consistently produce high-quality conformers for meaningful downstream applications. To address these issues, we introduce CoarsenConf, which coarse-grains molecular graphs based on torsional angles and integrates them into an SE(3)-equivariant hierarchical variational autoencoder. Through equivariant coarse-graining, we aggregate the fine-grained atomic coordinates of subgraphs connected via rotatable bonds, creating a variable-length coarse-grained latent representation. Our model uses a novel aggregated attention mechanism to restore fine-grained coordinates from the coarse-grained latent representation, enabling efficient generation of accurate conformers. Furthermore, we evaluate the chemical and biochemical quality of our generated conformers on multiple downstream applications, including property prediction and large-scale oracle-based protein docking. Overall, CoarsenConf generates more accurate conformer ensembles compared to prior generative models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRSDate not supplied

Dust Particle Aeroheating Calculations for Mars Entry Hypersonic Flows

The purpose of the current research effort is to predict particle trajectories in a hypersonic flowfield while accounting for all relevant particle-fluid and particle-particle interactions. The resulting particle solver is titled DUst Simulation & Tracking (DUST) and works in conjunction with the US3D CFD flow solver. The key elements underpinning the current work are an efficient mesh-localiztion algorithm, time-driven hard-sphere model, point-to-point MPI framework, coarse-graining using computational parcels, and high-order Adams-Bashforth time stepping. The coupled US3D-DUST framework will be applied to conduct a multi-physics examination of dust-laden flows around the the Schiarparelli capsule.

Dusty Flows↗

NASA NTRSDate not supplied

Dust Particle Aeroheating Calculations for Mars Entry Hypersonic Flows

Dusty Flows↗

DOE OSTI2024

A Scale‐Dependent Analysis of the Barotropic Vorticity Budget in a Global Ocean Simulation

Abstract The climatological mean barotropic vorticity budget is analyzed to investigate the relative importance of surface wind stress, topography, planetary vorticity advection, and nonlinear advection in dynamical balances in a global ocean simulation. In addition to a pronounced regional variability in vorticity balances, the relative magnitudes of vorticity budget terms strongly depend on the length‐scale of interest. To carry out a length‐scale dependent vorticity analysis in different ocean basins, vorticity budget terms are spatially coarse‐grained. At length‐scales greater than 1,000 km, the dynamics closely follow the Topographic‐Sverdrup balance in which bottom pressure torque, surface wind stress curl and planetary vorticity advection terms are in balance. In contrast, when including all length‐scales resolved by the model, bottom pressure torque and nonlinear advection terms dominate the vorticity budget (Topographic‐Nonlinear balance), which suggests a prominent role of oceanic eddies, which are of km in size, and the associated bottom pressure anomalies in local vorticity balances at length‐scales smaller than 1,000 km. Overall, there is a transition from the Topographic‐Nonlinear regime at scales smaller than 1,000 km to the Topographic‐Sverdrup regime at length‐scales greater than 1,000 km. These dynamical balances hold across all ocean basins; however, interpretations of the dominant vorticity balances depend on the level of spatial filtering or the effective model resolution. On the other hand, the contribution of bottom and lateral friction terms in the barotropic vorticity budget remains small and is significant only near sea‐land boundaries, where bottom stress and horizontal viscous friction generally peak.

Khatri, Hemant↗

DOE OSTI2024

Complex orders and chirality in the classical Kitaev-Γ model

It is well recognized that the low-energy physics of many Kitaev materials is governed by two dominant energy scales, the Ising-type Kitaev coupling 𝐾 and the symmetric off-diagonal Γ coupling. An understanding of the interplay between these two scales is therefore the natural starting point toward a quantitative description that includes subdominant perturbations that are inevitably present in real materials. This study focuses on the classical 𝐾−Γ model on the honeycomb lattice, with a specific emphasis on the region 𝐾< 0 and Γ > 0 , which is the most relevant for the available materials and which remains enigmatic in both quantum and classical limits, despite much effort. We employ large-scale Monte Carlo simulations on specially designed finite-size clusters and unravel the presence of a complex multisublattice magnetic order in a wide region of the phase diagram, whose structure is characterized in detail. We show that this order can be quantified in terms of a coarse-grained scalar-chirality order, featuring a counterrotating modulation on the two spin sublattices. Here, we also provide a comparison to previous studies and discuss the impact of quantum fluctuations on the phase diagram.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2011

Software Testing and Verification in Climate Model Development

Over the past 30 years most climate models have grown from relatively simple representations of a few atmospheric processes to a complex multi-disciplinary system. Computer infrastructure over that period has gone from punch card mainframes to modem parallel clusters. Model implementations have become complex, brittle, and increasingly difficult to extend and maintain. Existing verification processes for model implementations rely almost exclusively upon some combination of detailed analysis of output from full climate simulations and system-level regression tests. In additional to being quite costly in terms of developer time and computing resources, these testing methodologies are limited in terms of the types of defects that can be detected, isolated and diagnosed. Mitigating these weaknesses of coarse-grained testing with finer-grained "unit" tests has been perceived as cumbersome and counter-productive. In the commercial software sector, recent advances in tools and methodology have led to a renaissance for systematic fine-grained testing. We discuss the availability of analogous tools for scientific software and examine benefits that similar testing methodologies could bring to climate modeling software. We describe the unique challenges faced when testing complex numerical algorithms and suggest techniques to minimize and/or eliminate the difficulties.

Clune, Thomas L.↗

NASA NTRS2002

Quench Crack Behavior of Nickel-base Disk Superalloys

There is a need to increase the temperature capability of superalloy turbine disks to allow higher operating temperatures in advanced aircraft engines. When modifying processing and chemistry of disk alloys to achieve this capability, it is important to preserve the ability to use rapid cooling during supersolvus heat treatments to achieve coarse grain, fine gamma prime microstructures. An important step in this effort is an understanding of the key variables controlling the cracking tendencies of nickel-base disk alloys during quenching from supersolvus heat treatments. The objective of this study was to investigate the quench cracking tendencies of several advanced disk superalloys during simulated heat treatments. Miniature disk specimens were rapidly quenched after solution heat treatments. The responses and failure modes were compared and related to the quench cracking tendencies of actual disk forgings. Cracking along grain boundaries was generally observed to be operative. For the alloys examined in this study, the solution temperature not alloy chemistry was found to be the primary factor controlling quench cracking. Alloys with high solvus temperatures show greater tendency for quench cracking.

Gayda, John↗

DOE OSTI2025

An adaptive, data-driven multiscale approach for dense granular flows

The accuracy of coarse-grained continuum models of dense granular flows is limited by the lack of high-fidelity closure models for granular rheology. One approach to addressing this issue, referred to as the hierarchical multiscale method, is to use a high-fidelity fine-grained model to compute the closure terms needed by the coarse-grained model. The difficulty with this approach is that the overall model can become computationally intractable due to the high computational cost of the high-fidelity model. In this work, we describe a multiscale modeling approach for dense granular flows that utilizes neural networks trained using high-fidelity discrete element method (DEM) simulations to approximate the constitutive granular rheology for a continuum incompressible flow model. Our approach leverages an ensemble of neural networks to estimate predictive uncertainty that allows us to determine whether the rheology at a given point is accurately represented by the neural network model. Additional DEM simulations are only performed when needed, minimizing the number of additional DEM simulations required when updating the rheology. This adaptive coupling significantly reduces the overall computational cost of the approach while controlling the error. In addition, the neural networks are customized to learn regularized rheological behavior to ensure well-posedness of the continuum solution. We first validate the approach using two-dimensional steady-state and decelerating inclined flows. We then demonstrate the efficiency of our approach by modeling three-dimensional sub-aerial granular column collapse for varying initial column aspect ratios, where our multiscale method compares well with the computationally expensive computational fluid dynamics (CFD)-DEM simulation.

Dense granular flows↗

DOE OSTI2024

Deciphering the Cofilin Oligomers via Intermolecular Disulfide Bond Formation: A Coarse-Grained Molecular Dynamics Approach to Understanding Cofilin’s Regulation on Actin Filaments

Cofilin, a key actin-binding protein, orchestrates the dynamics of the actomyosin network through its actin-severing activity and by promoting the recycling of actin monomers. Recent experimental work suggests that cofilin also forms functionally distinct oligomers through thiol post-translational modification (PTM) that encourages actin nucleation and assembly. Despite these advances, the structural conformations of cofilin oligomers that modulate actin activity remain elusive because there are combinatorial ways to oxidize thiols in cysteines to form disulfide bonds rapidly. This study employs molecular dynamics simulations to investigate human cofilin 1 as a case study for exploring cofilin dimers via disulfide bond formation. Using the free energy profiling, our simulations unveil a range of probable cofilin dimer structures not represented in current Protein Data Bank entries. These candidate dimers are characterized by their distinct population distributions and relative free energies. Of particular note is a dimer featuring an interface between cysteines 139 and 147 residues, which demonstrates stable free energy characteristics and intriguingly symmetrical geometry. In contrast, the experimentally proposed dimer structure exhibits a less stable free energy profile. Here, we also evaluate frustration quantification based on the energy landscape theory in the protein-protein interactions at the dimer interfaces. Notably, the 39-39 dimer configuration emerges as a promising candidate for forming cofilin tetramers, as substantiated by frustration analysis. Additionally, docking simulations with actin filaments further evaluate the stability of these cofilin dimer-actin complexes. Our findings thus offer a computational framework for understanding the role of thiol post-translational modification cofilin proteins in regulating oligomerization, and the subsequent cofilin-mediated actin dynamics in the actomyosin network.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2001

The Tensile Properties of Advanced Nickel-Base Disk Superalloys During Quenching Heat Treatments

There is a need to increase the temperature capabilities of superalloy turbine disks. This would allow full utilization of higher temperature combustor and airfoil concepts under development. One approach to meet this goal is to modify the processing and chemistry of advanced alloys, while preserving the ability to use rapid cooling supersolvus heat treatments to achieve coarse grain, fine gamma prime microstructures. An important step in this effort is to understand the key high temperature tensile properties of advanced alloys as they exist during supersolvus heat treatments. This could help in projecting cracking tendencies of disks during quenches from supersolvus heat treatments. The objective of this study was to examine the tensile properties of two advanced disk superalloys during simulated quenching heat treatments. Specimens were cooled from the solution heat treatment temperatures at controlled rates, interrupted, and immediately tensile tested at various temperatures. The responses and failure modes were compared and related to the quench cracking tendencies of disk forgings.

Gabb, Timothy P.↗

NASA NTRS1998

An Assessment of the Axial and Radial Dilation of a DPIMS Tantalum Cartridge for Space Shuttle Flight Experiments

Ground-based heat treatment tests are planned on an argon gas-filled tantalum cartridge developed as pan of a Diffusion Processes in Molten Semiconductors (DPIMS) experiment conducted on NASA's Space Shuttle. The possibility that the cartridge may creep during testing and touch the furnace walls is of real concern in this program. The present paper discusses the results of calculations performed to evaluate this possibility. Deformation mechanism maps were constructed using literature data in order to identify the creep mechanism dominant under the appropriate stresses and temperatures corresponding to the test conditions. These results showed that power-law creep was dominant when the grain size of the material exceeded 55 gm but Coble creep was the important mechanism below this value of grain size. Finite element analysis was used to analyze the heat treatment tWs assuming a furnace run away condition (which is a worst case scenario) using the appropriate creep parameters corresponding to grain sizes of 1 and 100 gm. Calculations were also conducted to simulate the effect of an initial 3 tilt of the cartridge assembly, the maximum possible tilt angle. The von Mises stress and su-ain distributions were calculated assuming that the cartridge was fixed at one end as it was heated from ambient temperature to 1823 K in 1.42 h, maintained at 1823 K for 9.5 h and then further heated to an over temperature condition of 2028 K in 0.3 h. The inelastic axial and radial displacements of the cartridge walls were evaluated by resolving the von Mises strain along the corresponding directions. These calculations reveal that the maximum axial and radial displacements are expected to be about 2.9 and 0.25 mm, respectively, for both fine and coarse-grained materials at 2028 K. It was determined that these displacements occur during heat-up to temperature and creep of the cartridge is likely to be relatively insignificant irrespective of grain size. Furthermore, with a 3' tilt of the cartridge, the deflection is increased by only 0.39 gm which is negligible. Since the gap between the furnace heating elements and the cartridge is about 7.5 mm and less than the maximum radial dilation of 0.25 mm at 2028 K, it is concluded that the cartridge is unlikely to touch the furnace walls during the experiments.

Raj, S.V.↗

DOE OSTI2024

VOTCA: multiscale frameworks for quantum and classical simulations in soft matter

Many physical phenomena in liquids and soft matter are multiscale by nature and can involve processes with quantum and classical degrees of freedom occurring over a vast range of length- and timescales. Examples range from structure formation processes of complex polymers or even polymer blends (Svaneborg & Everaers, 2023) on the classical side to charge and energy transport and conversion processes (Lee et al., 2019) involving explicit electronic and, therefore, quantum information. The Versatile Object-oriented Toolkit for Coarse-graining Applications (VOTCA) provides multiscale frameworks built on a comprehensive set of methods for the development of classical coarse-grained potentials (VOTCA-CSG) as well as state-of-the art excited state electronic structure methods based on density-functional and many-body Green’s function theories, coupled in mixed quantum-classical models and used in kinetic network models (VOTCA-XTP).

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Quantum real-time evolution using tensor renormalization group methods

We introduce an approach for approximate real-time evolution of quantum systems using tensor renormalization group (TRG) methods originally developed for imaginary time. We use higher-order TRG to generate a coarse-grained time evolution operator for a 1+1⁢D transverse Ising model with a longitudinal field. We show that the standard tensor norm used for the singular value decomposition-based truncation is degenerate and propose an alternate method to discriminate. We show that it is effective and efficient in evolving Gaussian wave packets for one and two particles in the disordered phase, while ordered phase behavior is more challenging to capture. We compare our algorithm with local simulators for universal quantum computers and discuss possible benchmarking in the near future.

lattice gauge theory↗

DOE OSTI2026

Can classical DEM simultaneously capture compressibility and flowability of milled biomass?

Accurate prediction of the rheological behavior of biomass is essential for the design and operation of hoppers, feeders, and storage systems in biorefineries. This study examines whether the classical, coarse-grained discrete element method (DEM) formulation can simultaneously reproduce the compressibility and flowability of milled herbaceous biomass, using Miscanthus × giganteus as a representative material. The model represents particles as rigid spheres interacting through Hertz-Mindlin elastic-frictional contacts augmented with an area-dependent cohesion term. Laboratory cyclic compression and wedge-shaped hopper discharge experiments were used as calibration benchmarks. Although the model can independently reproduce each behavior by appropriately tuning particle Young's modulus E and cohesion energy density k, an extensive parametric investigation comprising more than 600 simulations reveals that the optimal parameter regions for compression and hopper flow are distinct and non-overlapping in (E, k) space. Surrogate surface analysis further shows that the corresponding objective-function valleys exhibit similar trends but are approximately parallel and spatially offset, precluding a unified calibration within the explored domain. Sensitivity analysis indicates that compressibility is governed predominantly by stiffness and cohesion, whereas the slope of the mass flow rate-opening relation in hopper discharge is primarily controlled by tangential friction. Extensions incorporating particle size distribution and clumped-sphere representations do not eliminate the incompatibility. These results systematically reveal, for the first time, the structural limitation of simplified DEM formulations in representing biomass rheological behavior, underscoring the necessity for models incorporating additional physical mechanisms, such as particle deformability or enhanced interlocking, to achieve unified predictive capability for biomass handling behavior.

09 BIOMASS FUELS↗

NASA NTRS1992

Lunar resources: Oxygen from rocks and soil

The first set of hydrogen reduction experiments to use actual lunar material was recently completed. The sample, 70035, is a coarse-grained vesicular basalt containing 18.46 wt. percent FeO and 12.97 wt. percent TiO2. The mineralogy includes pyroxene, ilmenite, plagioclase, and minor olivine. The sample was crushed to a grain size of less than 500 microns. The crushed basalt was reduced with hydrogen in seven tests at temperatures of 900-1050 C and pressures of 1-10 atm for 30-60 minutes. A capacitance probe, measuring the dew point of the gas stream, was used to follow reaction progress. Experiments were also conducted using a terrestrial basalt similar to some lunar mare samples. Minnesota Lunar Simulant (MLS-1) contains 13.29 wt. percent FeO, 2.96 wt. percent Fe2O3, and 6.56 wt. percent TiO2. The major minerals include plagioclase, pyroxene, olivine, ilmenite, and magnetite. The rock was ground and seived, and experiments were run on the less than 74- and 500-1168-micron fractions. Experiments were also conducted on less than 74-micron powders of olivine, pyroxene, synthetic ilmenite, and TiO2. The terrestrial rock and mineral samples were reduced with flowing hydrogen at 1100 C in a microbalance furnace, with reaction progress monitored by weight loss. Experiments were run at atmospheric pressure for durations of 3-4 hr. Solid samples from both sets of experiments were analyzed by Mossbauer spectroscopy, petrographic microscopy, scanning electron microscopy, tunneling electron microscopy, and x-ray diffraction. Apollo 17 soil 78221 was examined for evidence of natural reduction in the lunar environment. This sample was chosen based on its high maturity level (I sub s/FeO = 93.0). The FeO content is 11.68 wt. percent and the TiO2 content is 3.84 wt. percent. A polished thin section of the 90-150 micron size fraction was analyzed by petrographic microscopy and scanning electron microscopy.

Allen, C. C.↗

NASA NTRSDate not supplied

Coarsening Up: Expanding the Particle Size Distributions of Lunar Simulants to Encompass A Comprehensive Range of Regolith Granularity

Due to a paucity of lunar regolith samples, simulated granular materials (simulants) have been developed and deployed for a variety of engineering and research purposes in an attempt to meet the requirements for lunar exploration. Broadly speaking, these manufactured materials aim to replicate the compositional (i.e., mineralogy and chemistry) and physical (e.g., particle size distribution, particle shape) profiles of lunar soil samples collected during the Apollo missions. While creating a “near-perfect” lunar simulant is prohibitive due to the unique formative processes (i.e., meteorite impacts, solar wind implantation) and conditions (anhydrous, reducing environment and exposure to radiation on the surface over billions of years) encountered on the moon, certain granular properties of the regolith can be recreated with a high level of fidelity by leveraging contemporary technologies. A fundamental physical property and descriptor, grain size can be correlated between samples and replicated through processing; however, a holistic understanding of the particle size distribution of lunar regolith must be considered when creating analog materials. This study produces particle size distributions for Apollo 16 returned samples, incorporating the sieved out >1 cm coarse-grained fraction, in order to provide a more comprehensive characterization of surficial lunar regolith grain sizes for the development of highlands-type granular analog testing materials.

R. N. Kovtun↗

NASA NTRS2003

Performance Analysis of a Hybrid Overset Multi-Block Application on Multiple Architectures

This paper presents a detailed performance analysis of a multi-block overset grid compu- tational fluid dynamics app!ication on multiple state-of-the-art computer architectures. The application is implemented using a hybrid MPI+OpenMP programming paradigm that exploits both coarse and fine-grain parallelism; the former via MPI message passing and the latter via OpenMP directives. The hybrid model also extends the applicability of multi-block programs to large clusters of SNIP nodes by overcoming the restriction that the number of processors be less than the number of grid blocks. A key kernel of the application, namely the LU-SGS linear solver, had to be modified to enhance the performance of the hybrid approach on the target machines. Investigations were conducted on cacheless Cray SX6 vector processors, cache-based IBM Power3 and Power4 architectures, and single system image SGI Origin3000 platforms. Overall results for complex vortex dynamics simulations demonstrate that the SX6 achieves the highest performance and outperforms the RISC-based architectures; however, the best scaling performance was achieved on the Power3.

Djomehri, M. Jahed↗