Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D 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 1,675 records · Page 93
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Benchmark microgravity experiments and computations for 3D dendritic-array stability in directional solidification

In this study, we present a comprehensive quantitative analysis of stability bands for dendritic arrays during directional solidification of a transparent succinonitrile-0.46 wt % camphor alloy, spanning a broad range of pulling velocities. Taking advantage of the microgravity environment aboard the International Space Station where most convection effects are suppressed, we obtain unique measurements that quantify the stable primary spacing range of spatially extended three-dimensional dendritic array structures under purely diffusive growth conditions. Through carefully designed velocity jump experiments and detailed examination of sub-grain boundary dynamics, we characterize key instabilities, including elimination and tertiary branching, shedding new light on the mechanisms governing dynamic dendritic spacing selection in extended 3D arrays. Phase field simulations are performed to characterize the stability limits of dendritic array structures for quantitative comparison with the flight experiments. Although the simulations capture general trends, significant deviations are noted at the upper stability boundary, indicating the influence of additional, unexplored factors. These findings contribute to a deeper understanding of dendritic growth dynamics and offer valuable benchmark data that could aid in refining predictive models and improving control of dendritic microstructures in metallurgical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Enhanced MPM framework with multipatch isogeometric analysis for geotechnical applications

Achieving stable stress solutions at large strains using the Material Point Method (MPM) is challenging due to the accumulation of errors associated with geometry discretization, cell-crossing noise, and volumetric locking. Several simplified attempts exist in the literature to mitigate these errors, including higher-order frameworks. However, the stability of the MPM solution in such frameworks has been limited to simple geometries and the single-phase formulation (i.e., neglecting pore fluid). Although never explored, multipatch isogeometric analysis offers desirable qualities to simulate complex geometries while mitigating errors in the MPM. The degree of required high-order spatial integration has also never been investigated to infer a minimum limit for the stability of the stress solution in MPM. This paper presents a general-purpose numerical framework for simulating stable stresses in porous media, capturing both near incompressibility and multiphase interactions. First, the numerical framework is presented considering Non-Uniform Rational B-splines (NURBS) to perform isogeometric analysis (IGA) in MPM. Additionally, a volumetric strain smoothing algorithm is used to alleviate errors associated with volumetric locking. Second, the manifestation of cell-crossing errors is assessed via a series of problems with orders ranging from linear to cubic interpolation functions. Third, the use of NURBS is investigated and verified for problems with circular geometries. Finally, multipatch analysis is deployed to simulate plane strain and 3D penetration in soils, considering nearly incompressible elastoplastic (total stress) analysis and fully-coupled hydro-mechanical (effective stress) analysis. The stability of the solution is also analyzed for different constitutive models. From the results, it can be concluded that the framework using cubic interpolation functions with strain smoothing is the most convenient, presenting stable stress solutions for a broad range of multiphase geotechnical applications.

58 GEOSCIENCES↗
DOE OSTI2025

Influence of Initial fabric and water-wetting on particle fracture and force chain evolution in natural silica sand

Particle fracture has a significant influence on the engineering behavior of granular materials. However, the combined influence of the initial fabric and wetting on particle fracture and force chain evolution in granular materials, particularly in sands, remains insufficiently explored in the current literature. Here, in this paper, one-dimensional (1D) confined compression experiments were conducted on dry and wet specimens, and particle-level fracture was captured using three-dimensional (3D) in-situ Synchrotron Micro-Computed Tomography (SMT). A quantitative assessment of particle fracture revealed that wet specimens, which have disturbed fabric, exhibited a higher fracture percentage in comparison to dry specimens. Complementary 3D finite element (FE) simulations were performed using dry and wet (only solid fabric was considered and the influence of pore water is not considered) sand assemblies to assess interparticle contact forces and particle-scale stress distributions within the specimens. Force chain analysis was subsequently conducted, encompassing quantification of interparticle forces, characterization of force network structures, and monitoring the dynamic evolution of force chains under different strain levels. The results show that the specimens with disturbed fabric led to a more dynamic and less persistent force network, more fabric instability, and thus more reorganization of force chain structures. In addition, the frequent rearrangement of the force network in the presence of water (with reduced inter-particle contact friction) likely exacerbates localized stress concentrations, promoting failure in previously unengaged or weakly connected particles. The results reported in this paper offer a new insight into how the initial fabric and wetting cause different fracture behavior. The findings can also pave the way for more in-depth future investigations into the mechanics governing particle fracture in granular assemblies.

Finite element analysis↗
DOE OSTI2017

An unstructured-mesh finite-volume MPDATA for compressible atmospheric dynamics

An advancement of the unstructured-mesh finite-volume MPDATA (Multidimensional Positive Definite Advection Transport Algorithm) is presented that formulates the error-compensative pseudo-velocity of the scheme to rely only on face-normal advective fluxes to the dual cells, in contrast to the full vector employed in previous implementations. This is essentially achieved by expressing the temporal truncation error underlying the pseudo-velocity in a form consistent with the flux-divergence of the governing conservation law. The development is especially important for integrating fluid dynamics equations on non-rectilinear meshes whenever face-normal advective mass fluxes are employed for transport compatible with mass continuity—the latter being essential for flux-form schemes. In particular, the proposed formulation enables large-time-step semi-implicit finite-volume integration of the compressible Euler equations using MPDATA on arbitrary hybrid computational meshes. Furthermore, it facilitates multiple error-compensative iterations of the finite-volume MPDATA and improved overall accuracy. The advancement combines straightforwardly with earlier developments, such as the nonoscillatory option, the infinite-gauge variant, and moving curvilinear meshes. A comprehensive description of the scheme is provided for a hybrid horizontally-unstructured vertically-structured computational mesh for efficient global atmospheric flow modelling. The proposed finite-volume MPDATA is verified using selected 3D global atmospheric benchmark simulations, representative of hydrostatic and non-hydrostatic flow regimes. Besides the added capabilities, the scheme retains fully the efficacy of established finite-volume MPDATA formulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Maximum bound principle preserving integrating factor Runge–Kutta methods for semilinear parabolic equations

A large class of semilinear parabolic equations satisfy the maximum bound principle (MBP) in the sense that the time-dependent solution preserves for any time a uniform pointwise bound imposed by its initial and boundary conditions. Here, the MBP plays a crucial role in understanding the physical meaning and the well-posedness of the mathematical model. Investigation on numerical algorithms with preservation of the MBP has attracted increasingly attentions in recent years, especially for the temporal discretizations, since the violation of MBP may lead to nonphysical solutions or even blow-ups of the algorithms. In this paper, we study high-order MBP-preserving time integration schemes by means of the integrating factor Runge-Kutta (IFRK) method. Beginning with the space-discrete system of semilinear parabolic equations, we present the IFRK method in general form and derive the sufficient conditions for the method to preserve the MBP. In particular, we show that the classic four-stage, fourth-order IFRK scheme is MBP preserving for some typical semilinear systems although not strong stability preserving, which can be instantly applied to the Allen-Cahn type of equations. To our best knowledge, this is the first time to present a fourth-order linear numerical method preserving the MBP. In addition, convergence of these numerical schemes is proved theoretically and verified numerically, as well as their efficiency by simulations of 2D and 3D long-time evolutional behaviors. Numerical experiments are also carried out for a model which is not a typical gradient flow as the Allen-Cahn type of equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Phase contrast imaging measurements and numerical simulations of turbulent density fluctuations in gas-fuelled ECRH discharges in Wendelstein 7-X

The fundamental nature of turbulent density fluctuations in standard Wendelstein 7-X (W7-X) stellarator discharges is investigated experimentally via phase contrast imaging (PCI) in combination with gyrokinetic simulations with the code GENE. We find that density fluctuations are ion-temperature-gradient-driven and radially localised in the outer half of the plasma. It is shown that the line-integrated PCI measurements cover the right range of wavenumbers and a favourable toroidal and poloidal location to capture some of the strongest density fluctuations in W7-X. Due to the radial localisation of fluctuations, measured wavenumber–frequency spectra exhibit a dominant phase velocity, which can be related to the E × B rotation velocity at the radial position of a well in the neoclassical radial electric field. The match is robust against variations of heating power and line-integrated density, which is partly due to the localisation of fluctuations and partly due to effects of the radial gradient in the E × B velocity profile on the wavenumber–frequency spectrum. The latter effect is studied with a newly built synthetic PCI diagnostic and global gyrokinetic simulations with GENE-3D.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Linking Silica Support Morphology to the Dynamics of Aminopolymers in Composites

A combined computational and experimental approach is used to elucidate the effect of silica support morphology on polymer dynamics and CO 2 adsorption capacities in aminopolymer/silica composites. Furthermore, simulations are based on coarse-grained molecular dynamics simulations of aminopolymer composites where a branched aminopolymer, representing poly(ethylenimine) (PEI), is impregnated into different silica mesoporous supports. The morphology of the mesoporous supports varies from hexagonally packed cylindrical pores representing SBA-15, double gyroids representing KIT-6 and MCM-48, and cagelike structures representing SBA-16. In parallel, composites of PEI and the silica supports SBA-15, KIT-6, MCM-48, and SBA-16 are synthesized and characterized, including measuring their CO 2 uptake. Simulations predict that a 3D pore morphology, such as those of KIT-6, MCM-48, and SBA-16, will have faster segmental mobility and have lower probability of primary amine and surface silanol associations, which should translate to higher CO 2 uptake in comparison to a 2D pore morphology such as that of SBA-15. We found that KIT-6 has higher CO 2 uptake than SBA-15 at equivalent PEI loading, even though both supports have similar surface area and pore volume. But, this is not the case for the MCM-48 support, which has smaller pores, and SBA-16, whose pore structure rapidly degrades after PEI impregnation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

X-ray Micro-Computed Tomography for Structural Analysis of All-Solid-State Battery at Pouch Cell Level

Characterizing the microstructure of all-solid-state batteries (ASSBs) during fabrication and operation is vital for their advancement, particularly as scaling to pouch cell levels introduces challenges in probing large-scale microstructural evolution. This work highlights the potential of synchrotron X-ray micro-computed tomography (sXCT) as a nondestructive, rapid (<30 min), and high-resolution technique for visualizing and quantifying key microstructural features, including overhang, porosity, contact loss, active surface area, and tortuosity, in all-solid-state pouch cells. The large field of view (up to millimeters) of sXCT enables detailed analysis at an industry-relevant scale, bridging the gap between laboratory research and commercial applications. Furthermore, integrating realistic sXCT-derived 3D models into multiphysics simulations could provide insights into chemo-mechanical degradation, particularly at the edges of the pouch cells, offering a pathway for designing robust, high-performance ASSBs. This perspective establishes sXCT as an indispensable tool for advancing both the understanding and the engineering of next-generation energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2021

Stress Controls Rupture Extent and Maximum Magnitude of Induced Earthquakes

Abstract Seismic hazard forecasts of induced seismicity often require estimates of the maximum possible magnitude ( M max ). Empirical models suggest that maximum magnitudes, or expected number of earthquakes, are related to the volume of injected fluid. We perform a suite of 3D physics‐based earthquake simulations with rate‐ and state‐dependent friction, systematically varying the area of the pressurized region and the amplitude of the initial homogeneous or heterogeneous shear stress. Using the resulting catalog we explore the conditions that result in pressure‐controlled versus runaway ruptures that extend outside the pressurized zone. We find that proposed empirical scaling laws correctly predict M max when shear stresses are further from failure (≤90% of maximum shear stress) and for high amplitude stress fields. Runaway ruptures are observed for higher initial shear stresses and smoother stress fields. In these cases, runaway ruptures occur early after the onset of injection and rarely preceded by foreshock activity.

58 GEOSCIENCES↗
DOE OSTI2023

Three‐Dimensional Magnetic Reconnection Spreading in Current Sheets of Non‐Uniform Thickness

Abstract Magnetic reconnection in naturally occurring and laboratory settings often begins locally and elongates, or spreads, in the direction perpendicular to the plane of reconnection. Previous work has largely focused on current sheets with a uniform thickness, for which the predicted spreading speed for anti‐parallel reconnection is the local speed of the current carriers. We derive a scaling theory of three‐dimensional (3D) spreading of collisionless anti‐parallel reconnection in a current sheet with its thickness varying in the out‐of‐plane direction, both for spreading from a thinner to thicker region and a thicker to thinner region. We derive an expression for calculating the time it takes for spreading to occur for a current sheet with a given profile of its thickness. A key result is that when reconnection spreads from a thinner to a thicker region, the spreading speed in the thicker region is slower than both the Alfvén speed and the speed of the local current carriers by a factor of the ratio of thin to thick current sheet thicknesses. This is important because magnetospheric and solar observations have previously measured the spreading speed to be slower than previously predicted, so the present mechanism might explain this feature. We confirm the theory via a parametric study using 3D two‐fluid numerical simulations. We use the prediction to calculate the time scale for reconnection spreading in Earth's magnetotail during geomagnetic activity. The results are also potentially important for understanding reconnection spreading in solar flares and the dayside magnetopause of Earth and other planets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Atomic resolution coherent x-ray imaging with physics-based phase retrieval

Coherent x-ray imaging and scattering from accelerator based sources such as synchrotrons continue to impact biology, medicine, technology, and materials science. Many synchrotrons around the world are currently undergoing major upgrades to increase their available coherent x-ray flux by approximately two orders of magnitude. The improvement of synchrotrons may enable imaging of materials in operando at the atomic scale which may revolutionize battery and catalysis technologies. Current algorithms used for phase retrieval in coherent x-ray imaging are based on the projection onto sets method. These traditional iterative phase retrieval methods will become more computationally expensive as they push towards atomic resolution and may struggle to converge. Additionally, these methods do not incorporate physical information that may additionally constrain the solution. In this work, we present an algorithm which incorporates molecular dynamics into Bragg coherent diffraction imaging (BCDI). This algorithm, which we call PRAMMol (Phase Retrieval with Atomic Modeling and Molecular Dynamics) combines statistical techniques with molecular dynamics to solve the phase retrieval problem. We present several examples where our algorithm is applied to simulated coherent diffraction from 3D crystals and show convergence to the correct solution at the atomic scale.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Scintillation event imaging with a single photon avalanche diode camera

Position and time measurements of scintillation events encode information about the radiation source. Single photon avalanche diode (SPAD) arrays offer multiple-megapixel spatial resolution and tens of picoseconds temporal resolution for detecting single photons. Current lensless designs for measuring scintillation events use sensors that are lower in spatial resolution. Camera-based designs use sensors that are lower in temporal resolution or readout rate and cannot image individual interactions. Here we propose to image scintillation events in a thick, monolithic scintillator using a high-resolution SPAD camera. We demonstrate that a commercial SPAD camera is able to gather sufficient signal to image individual scintillation events and observe 3D shifts in their spatial distribution. Simulations show that a SPAD camera can localize individual scintillation events in 3D. We report direct imaging of gamma-ray interactions in a scintillator with a SPAD camera. The proposed design may allow to measure complex signatures of individual particles interacting in the scintillator.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Violent mergers can explain the inflated state of some of the fastest stars in the Galaxy

A significant number of hypervelocity stars with velocities between 1500 − 2500 km s −1 have recently been observed. The only plausible explanation so far is that they were produced through thermonuclear supernovae in white dwarf binaries. Since these stars are thought to be surviving donors of Type Ia supernovae, a surprising finding was that these stars are inflated, with radii an order of magnitude higher than expected for Roche-lobe-filling donors. Recent attempts at explaining them have combined 3D hydrodynamical supernova explosion simulations with 1D stellar modelling to explain the impact of supernova shocks on runaway white dwarfs. However, only the hottest and most compact of those runaway stars can so far marginally be reproduced by detailed models of runaways from supernova explosions. In this and a companion paper, we introduce a new AREPO simulation of two massive CO white dwarfs that explode via a violent merger. During the merger, the primary white dwarf ignites when the secondary is on its last orbit and plunging towards the primary. In the corresponding aftermath, the core of the secondary white dwarf of 0.16 M ⊙ remains bound, moving at a velocity of ∼2800 km s −1 . We mapped this object into MESA and show that this runaway star can explain the observations of two hypervelocity stars that were dubbed D6-1 and D6-3 based on their original discovery motivated by the D6 scenario, though the violent merger scenario presented here is somewhat distinct from the D6 scenario.

Astronomy and AstroPhysics↗
DOE OSTI2013

Sponge-like Si-SiO{sub 2} nanocomposite—Morphology studies of spinodally decomposed silicon-rich oxide

Sponge-like Si nanostructures embedded in SiO{sub 2} were fabricated by spinodal decomposition of sputter-deposited silicon-rich oxide with a stoichiometry close to that of silicon monoxide. After thermal treatment a mean feature size of about 3 nm was found in the phase-separated structure. The structure of the Si-SiO{sub 2} nanocomposite was investigated by energy-filtered transmission electron microscopy (EFTEM), EFTEM tomography, and atom probe tomography, which revealed a percolated Si morphology. It was shown that the percolation of the Si network in 3D can also be proven on the basis of 2D EFTEM images by comparison with 3D kinetic Monte Carlo simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Improvement of extraction system geometry with suppression of possible Penning discharge ignition

During the past two years, a new ECR 2.45 GHz type ion source has been developed especially dedicated to intense light ion injector project like IPHI (Injecteur Proton Haute Intensité), IFMIF (International Fusion Materials Irradiation Facility), to reduce beam emittance at RFQ entrance by shortening the length of the LEBT. This new ALISES concept (Advanced Light Ion Source Extraction System) is based on the use of an additional LEBT short length solenoid very close to the extraction aperture. The fringe field of this new solenoid produces the needed magnetic field to create the ECR resonance in the plasma chamber. Such geometry allows first putting the solenoid at ground potential, while saving space in front of the extraction to move the first LEBT solenoid closer and focus earlier the intense extracted beam. During the commissioning of the source in 2011–2012, ALISES has produced about 20 mA extracted from a 6 mm diameter plasma extraction hole at 23 kV. But the magnetic configuration combined to the new extraction system geometry led to important Penning discharge conditions in the accelerator column. Lots of them have been eliminated by inserting glass pieces between electrodes to modify equipotential lines with unfavorable ExB vacuum zones where particles were produced and trapped. To study Penning discharge location, several 3D calculations have been performed with OPERA-3D/TOSCA code to simulate the possible production and trapping of electrons in the extraction system. The results obtained on different sources already built have shown very good agreement with sparks location observed experimentally on electrodes. The simulations results as well as experimental measurements are presented and solutions to prevent possible Penning discharge in future source geometries are established.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Investigation of ionization-induced electron injection in a wakefield driven by laser inside a gas cell

Ionization-induced electron injection was investigated experimentally by focusing a driving laser pulse with a maximum normalized potential of 1.2 at different positions along the plasma density profile inside a gas cell, filled with a gas mixture composed of 99%H 2 +1%N 2 . Changing the laser focus position relative to the gas cell entrance controls the accelerated electron bunch properties, such as the spectrum width, maximum energy, and accelerated charge. Simulations performed using the 3D particle-in-cell code WARP with a realistic density profile give results that are in good agreement with the experimental ones. Lastly, we discuss the interest of this regime for optimizing the bunch charge in a selected energy window.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Particle rings and astrophysical accretion discs

Norman Rostoker had a wide range of interests and significant impact on the plasma physics research at Cornell during the time he was a Cornell professor. His interests ranged from the theory of energetic electron and ion beams and strong particle rings to the related topics of astrophysical accretion discs. We outline some of the topics related to rings and discs including the Rossby wave instability which leads to formation of anticyclonic vortices in astrophysical discs. These vorticies are regions of high pressure and act to trap dust particles which in turn may facilitate planetesimals growth in proto-planetary disks and could be important for planet formation. Analytical methods and global 3D magneto-hydrodynamic simulations have led to rapid advances in our understanding of discs in recent years.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Hydrodynamic dispersion in thin channels with micro-structured porous walls

Flowand transport within porous- and microtextured-walled channels is relevant to a number of natural and industrial processes. Designing and optimizing the topology of the pores and/or microstructure to achieve target performance at the system scale (or macroscale) is still an open question. In this work, we study whether hydrodynamic dispersion in microfluidic channels with walls structured by obstacles can be modeled by dispersion in channels with porous walls described as continuous porous media of zero or finite permeability. We perform single phase microfluidic non-reactive flow experiments in channels embedded in micropatterns with different topologies. Specifically, we focus on transverse riblets and arrays of pillars as examples of impermeable and permeable obstructions, respectively. We compare the experimental results with three models: 3D pore-scale simulations which resolve the micropattern geometry explicitly and two upscaled models which treat the micropattern as a continuum of zero or finite permeability. This study demonstrates that polydimethylsiloxane micromodels with appropriately patterned surfaces can be successfully employed to validate various continuum-scale modeling approximations in different physical regimes, identified by the order of magnitude of the Peclet number and the obstruction permeability.

Ling, Bowen↗