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

3D modelling of MHD mixed convection flow in a vertical duct with transverse magnetic field and volumetric or surface heating

Understanding phenomena associated with the multiple effects/interactions of the fusion nuclear environment on liquid metal flow is required to correctly design liquid metal (LM) blankets for fusion facilities. These effects are investigated in the present work by numerically simulating 3D LM MHD flow. The simulated geometry consists of a straight, vertical duct which runs perpendicular to a strong, fringing applied magnetic field. There is also a region of applied heating as the primary goal is to explore buoyancy effects in MHD duct flows. Results are presented for both buoyancy assisted (upwards) and buoyancy opposed (downwards) flows in conducting and insulating ducts for a range of Hartmann numbers (Ha) up to 100, Reynolds numbers (Re) from 10 3 to 10 4 and Grashof (Gr) numbers from 10 7 to 10 8 . While increasing Gr or decreasing Re increases buoyancy effects, increasing Ha was shown to increase maximum temperature through turbulence reduction. The extent to which the MHD mixed convection flows are quasi-2D is analyzed and buoyant effects, in competition with electromagnetic forces, are shown to bring about 3D flow features not seen in purely MHD flows. Volumetric nuclear heating with steep gradients is applied to the vertical MHD flows for comparison to flows with surface heating only. Finally, surface heating generates stronger buoyancy effects than volumetric heating of the same total power; however, many of the same phenomena occur. Therefore, surface heating, the only option for lab experiments, can provide indication of the effects of volumetric heating in MHD flows.

42 ENGINEERING↗
DOE OSTI2025

Liquid crystal torons in Poiseuille-like flows

Three-dimensional (3D) simulations of the structure of liquid crystal (LC) torons, topologically protected distortions of the LC director field, under material flows are rare but essential in microfluidic applications. Here, we show that torons adopt a steady-state configuration at low flow velocity before disintegrating at higher velocities, in line with experimental results. Furthermore, we show that under partial slip conditions at the boundaries, the flow induces a reversible elongation of the torons, also consistent with the experimental observations. These results are in contrast with previous simulation results for 2D skyrmions under similar flow conditions, highlighting the need for a 3D description of this LC soliton in relation to its coupling to the material flow. These findings pave the way for future studies of other topological solitons, like hopfions and heliknotons, in flowing soft matter systems.

Science & Technology - Other Topics↗
DOE OSTI2017

Pedestal-to-wall 3D fluid transport simulations on DIII-D

The 3D fluid-plasma edge transport code EMC3-EIRENE is used to test several magnetic field models with and without plasma response against DIII-D experimental data for even and odd-parity n = 3 magnetic field perturbations. The field models include ideal and extended MHD equilibria, and the vacuum approximation. Plasma response is required to reduce the stochasticity in the pedestal region for even-parity fields, however too much screening suppresses the measured splitting of the downstream T-e profile. Odd-parity perturbations result in weak tearing and only small additional peaks in the downstream measurements. In this case plasma response is required to increase the size of the lobe structure. No single model is able to simultaneously reproduce the upstream and downstream characteristics for both odd and even-parity perturbations.

3D transport↗
DOE OSTI2019

Spatial decomposition of structured grids for nuclear reactor simulations

Spatial decomposition methods based on graph partitioning are developed and implemented in the high fidelity neutron transport code MPACT. These graph-based spatial decomposition methods are more general than previous decomposition methods and typically provide better load balance and reduced runtimes due to their improved parallel efficiency. Correlations are drawn between simulation runtime and the balance of the partition for 2D simulations. Comparisons are made using partition balance metrics for different decomposition schemes in 2D and 3D simulations. Finally, for typical ranges of subdomains, graph-based partitioning methods offer significant reductions to runtimes. However, for highly decomposed problems, these graph-based methods may decrease convergence rates, thus reducing parallel efficiency compared to older methods.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2017

Effects of dimensionality on kinetic simulations of laser-ion acceleration in the transparency regime

A particle-in-cell study of laser-ion acceleration mechanisms in the transparency regime illustrates how two-dimensional (2D) S and P simulations (laser polarization in and out of the simulation plane, respectively) capture different physics characterizing these systems, visible in their entirety in often cost-prohibitive three-dimensional (3D) simulations. The electron momentum anisotropy induced in the target by the laser pulse is dramatically different in the two 2D cases, manifested in differences in target expansion timescales, electric field strengths, and density thresholds for the onset of relativistically induced transparency. In particular, 2D-P simulations exhibit dramatically greater electron heating in the simulation plane, whereas 2D-S ones show a much more isotropic energy distribution, similar to 3D. An ion trajectory analysis allows one to isolate the fields responsible for ion acceleration and to characterize the acceleration regimes in time and space. The artificial longitudinal electron heating in 2D-P exaggerates the effectiveness of target-normal sheath acceleration into its dominant acceleration mechanism throughout the laser-plasma interaction, whereas 2D-S and 3D both have sizable populations accelerated preferentially during transparency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Non-equilibrium approach to prediction of microstructure evolution for metals undergoing severe plastic deformation

Highlights: • New consistent physically-based model of dynamic recrystallization. • The model is based on evolution of fraction of high-angle grain boundaries. • Possibility to predict the evolution of dislocation cells and grain boundaries. • The model is implemented for FEM simulations. • Results demonstrate inhomogeneous distribution of strain within the sample. - Abstract: Available models of dynamic recrystallization have a number of disadvantages that in most cases make them inapplicable for practical predictions of material microstructure evolution. Both the microstructural and the empirically based approaches do not reflect physical processes leading to evolution of material defect structure in the process of plastic deformation. This work presents an attempt to develop a consistent physically-based model of dynamic recrystallization. This model, accounting for physical nature of processes of material defect structure evolution, should provide a possibility to predict evolution of several different experimentally measurable parameters of material microstructure without introduction of big number of fitting parameters. It is suggested that such a model should be based on equation for evolution of fraction of high-angle grain boundaries (HAGBs) in the process of deformation. It is shown, that the new model gives a possibility to predict the evolution of dislocation cells and grain boundaries in copper-based alloys providing good coincidence with experimental observations. Full 3-dimensional numerical simulation of multidirectional forging of copper is performed utilizing the developed dynamic recrystallization model. The same 3D simulations demonstrate new noteworthy effects connected to inhomogeneous distribution of plastic strain within the bulk of the material and material strain hardening.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Modeling pulse-cleaning plasma mirrors from dielectric response to saturation: A particle-in-cell approach

Pulse-cleaning plasma mirrors are widely employed to improve ultraintense laser contrast, but much of the literature concerning their effect on the reflected pulse is empirical. A simulation study of pulse-cleaning plasma mirrors using the particle-in-cell code large scale plasma is presented here. The importance of capturing initial ionization from neutral atoms, collisional effects, and simulation dimensionality is considered. Excellent agreement with experimental data is obtained when a multiphoton ionization model is employed. Furthermore, a series of 2D simulations is shown to accurately replicate both the reflected light intensity and mode obtained from full 3D simulations at significantly reduced computation cost.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

The Role of Dimensionality in Understanding Granuloma Formation

Within the first 2–3 months of a Mycobacterium tuberculosis (Mtb) infection, 2–4 mm spherical structures called granulomas develop in the lungs of the infected hosts. These are the hallmark of tuberculosis (TB) infection in humans and non-human primates. A cascade of immunological events occurs in the first 3 months of granuloma formation that likely shapes the outcome of the infection. Understanding the main mechanisms driving granuloma development and function is key to generating treatments and vaccines. In vitro, in vivo, and in silico studies have been performed in the past decades to address the complexity of granuloma dynamics. This study builds on our previous 2D spatio-temporal hybrid computational model of granuloma formation in TB (GranSim) and presents for the first time a more realistic 3D implementation. We use uncertainty and sensitivity analysis techniques to calibrate the new 3D resolution to non-human primate (NHP) experimental data on bacterial levels per granuloma during the first 100 days post infection. Due to the large computational cost associated with running a 3D agent-based model, our major goal is to assess to what extent 2D and 3D simulations differ in predictions for TB granulomas and what can be learned in the context of 3D that is missed in 2D. Our findings suggest that in terms of major mechanisms driving bacterial burden, 2D and 3D models return very similar results. For example, Mtb growth rates and molecular regulation mechanisms are very important both in 2D and 3D, as are cellular movement and modulation of cell recruitment. The main difference we found was that the 3D model is less affected by crowding when cellular recruitment and movement of cells are increased. Overall, we conclude that the use of a 2D resolution in GranSim is warranted when large scale pilot runs are to be performed and if the goal is to determine major mechanisms driving infection outcome (e.g., bacterial load). To comprehensively compare the roles of model dimensionality, further tests and experimental data will be needed to expand our conclusions to molecular scale dynamics and multi-scale resolutions.

2D vs. 3D↗
DOE OSTI2016

SMR Re-Scaling and Modeling for Load Following Studies

This study investigates the creation of a new set of scaling parameters for the Oregon State University Multi-Application Small Light Water Reactor (MASLWR) scaled thermal hydraulic test facility. As part of a study being undertaken by Idaho National Lab involving nuclear reactor load following characteristics, full power operations need to be simulated, and therefore properly scaled. Presented here is the scaling analysis and plans for RELAP5-3D simulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

SMR Re-Scaling and Modeling for Load Following Studies

This study investigates the creation of a new set of scaling parameters for the Oregon State University Multi-Application Small Light Water Reactor (MASLWR) scaled thermal hydraulic test facility which modeled the NuScale power SMR design. As part of a study being undertaken by Idaho National Lab involving nuclear reactor load following characteristics, full power operations need to be simulated, and therefore properly scaled. Presented here is the scaling analysis and plans for RELAP5-3D simulation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1998

Virtual Reality Simulation of the Effects of Microgravity in Gastrointestinal Physiology

The ultimate goal of this research is to create an anatomically accurate three-dimensional (3D) simulation model of the effects of microgravity in gastrointestinal physiology and to explore the role that such changes may have in the pharmacokinetics of drugs given to the space crews for prevention or therapy. To accomplish this goal the specific aims of this research are: 1) To generate a complete 3-D reconstructions of the human GastroIntestinal (GI) tract of the male and female Visible Humans. 2) To develop and implement time-dependent computer algorithms to simulate the GI motility using the above 3-D reconstruction.

Compadre, Cesar M.↗
DOE OSTI2015

DYNAMO EFFECTS NEAR THE TRANSITION FROM SOLAR TO ANTI-SOLAR DIFFERENTIAL ROTATION

Numerical MHD simulations play an increasingly important role for understanding the mechanisms of stellar magnetism. We present simulations of convection and dynamos in density-stratified rotating spherical fluid shells. We employ a new 3D simulation code for obtaining the solution of a physically consistent anelastic model of the process with a minimum number of parameters. The reported dynamo simulations extend into a “buoyancy-dominated” regime where the buoyancy forcing is dominant while the Coriolis force is no longer balanced by pressure gradients, and strong anti-solar differential rotation develops as a result. We find that the self-generated magnetic fields, despite being relatively weak, are able to reverse the direction of differential rotation from anti-solar to solar-like. We also find that convection flows in this regime are significantly stronger in the polar regions than in the equatorial region, leading to non-oscillatory dipole-dominated dynamo solutions, and to a concentration of magnetic field in the polar regions. We observe that convection has a different morphology in the inner and the outer part of the convection zone simultaneously such that organized geostrophic convection columns are hidden below a near-surface layer of well-mixed highly chaotic convection. While we focus our attention on the buoyancy-dominated regime, we also demonstrate that conical differential rotation profiles and persistent regular dynamo oscillations can be obtained in the parameter space of the rotation-dominated regime even within this minimal model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Disorder effects in nitride semiconductors: impact on fundamental and device properties

Abstract Semiconductor structures used for fundamental or device applications most often incorporate alloy materials. In “usual” or “common” III–V alloys, based on the InGaAsP or InGaAlAs material systems, the effects of compositional disorder on the electronic properties can be treated in a perturbative approach. This is not the case in the more recent nitride-based GaInAlN alloys, where the potential changes associated with the various atoms induce strong localization effects, which cannot be described perturbatively. Since the early studies of these materials and devices, disorder effects have indeed been identified to play a major role in their properties. Although many studies have been performed on the structural characterization of materials, on intrinsic electronic localization properties, and on the impact of disorder on device operation, there are still many open questions on all these topics. Taking disorder into account also leads to unmanageable problems in simulations. As a prerequisite to address material and device simulations, a critical examination of experiments must be considered to ensure that one measures intrinsic parameters as these materials are difficult to grow with low defect densities. A specific property of nitride semiconductors that can obscure intrinsic properties is the strong spontaneous and piezoelectric fields. We outline in this review the remaining challenges faced when attempting to fully describe nitride-based material systems, taking the examples of LEDs. The objectives of a better understanding of disorder phenomena are to explain the hidden phenomena often forcing one to use ad hoc parameters, or additional poorly defined concepts, to make simulations agree with experiments. Finally, we describe a novel simulation tool based on a mathematical breakthrough to solve the Schrödinger equation in disordered potentials that facilitates 3D simulations that include alloy disorder.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Broadband (0–5 Hz) Fully Deterministic 3D Ground‐Motion Simulations of a Magnitude 7.0 Hayward Fault Earthquake: Comparison with Empirical Ground‐Motion Models and 3D Path and Site Effects from Source Normalized Intensities

We report on high‐performance computing (HPC) fully deterministic simulation of ground motions for a moment magnitude (⁠Mw⁠) 7.0 scenario earthquake on the Hayward fault resolved to 5 Hz using the SW4 finite‐difference code. We computed motions obeying physics‐based 3D wave propagation at a regional scale with an Mw 7.0 kinematic rupture model generated following Graves and Pitarka (2016). Both plane‐layered (1D) and 3D Earth models were considered, with 3D subsurface material properties and topography interpolated from a model of the U.S. Geological Survey (USGS). The resulting ground‐motion intensities cover a broader frequency range than typically considered in regional‐scale simulations, including higher frequencies relevant for engineering analysis of structures. Median intensities for sites across the domain are within the reported between‐event uncertainties (⁠τ⁠) of ground‐motion models (GMMs) across spectral periods 0.2–10 s (frequencies 0.1–5 Hz). The within‐event standard deviation ϕ of ground‐motion intensity measurement residuals range 0.2–0.5 natural log units with values consistently larger for the 3D model. Source‐normalized ratios of intensities (3D/1D) reveal patterns of path and site effects that are correlated with known geologic structure. These results demonstrate that earthquake simulations with fully deterministic wave propagation in 3D Earth models on HPC platforms produce broadband ground motions with median and within‐event aleatory variability consistent with empirical models. Systematic intensity variations for the 3D model caused by path and site effects suggest that these epistemic effects can be estimated and removed to reduce variation in site‐specific hazard estimates. This study motivates future work to evaluate the validity of the USGS 3D model and investigate the development of path and site corrections by running more scenarios.

Rodgers, Arthur J.↗
DOE OSTI2025

High-Resolution Simulations of Geological CO 2 Injection: Application to the SPE11 Benchmark

Geological carbon sequestration (GCS) will play a critical role in decarbonization and in facilitating the transition to clean energy systems. Because CO 2 is highly mobile, ensuring its safe and permanent injection into subsurface geological formations involves monitoring over larger spatial domains and longer time periods than is typical for hydrocarbon reservoirs. This can benefit from simulation tools capable of modeling key CO 2 trapping mechanisms, particularly those optimized for speed and scalability on high-performance computing systems. Using isothermal versions of the SPE11B and SPE11C benchmark cases, we conduct a mesh refinement study simulating CO 2 injection into kilometer-scale rock formations at centimeter resolution with the GEOS open-source simulation framework. We focus on how mesh refinement improves the accuracy of convective mixing in both 2D and 3D simulations. The computational costs associated with achieving a converged solution highlight the need for predictive upscaling techniques. A systematic performance scaling analysis—including both central processing unit (CPU) and graphics processing unit (GPU) architectures—complements the “Results” section.

Geosciences↗
DOE OSTI2020

Simulations of Laser Experiments to Study the Origin of Cosmic Magnetic Fields

Magnetic fields are ubiquitous in the universe; however, their origin is not fully understood. Cosmologists and astrophysicists have proposed a variety of ways in which small seed magnetic fields could be created. It is widely thought that the much larger values of the cosmic magnetic fields that we observe are a result of the amplification of these seed fields by the nonlinear turbulent dynamo mechanism. Such a mechanism had not yet been demonstrated in a controlled laboratory environment. We conceived experiments designed to demonstrate and study the nonlinear turbulent dynamo mechanism in the laboratory. These experiments characterize the distribution of turbulent energy among the velocity, magnetic field, and density fluctuations, providing a comprehensive picture of the energy cascade in a magnetized, turbulent plasma. The experiments build on the pathfinder experiments we have conducted on the Vulcan laser at the Rutherford-Appleton Laboratory in the UK and the OMEGA laser at Laboratory for Laser Energetics at Rochester. They utilize the high-intensity lasers at the Omega Laser Facility, the National Ignition Facility (NIF) at LLNL and the Laser Megajoule (LMJ) Facility in France – the largest laser facilities in the world. The goal of this project was to design and model the highly demanding experiments through simulation campaigns using FLASH, a highly capable radiation-MHD code we have developed, and large-scale 3D simulations on the Mira supercomputer at ANL. The simulations were vital to ensuring the experiments achieve the strong turbulence and large magnetic Reynolds numbers required for the nonlinear turbulent dynamo mechanism to operate. The simulations were also critical to determining when to fire the diagnostics, since the experiments last tens of nanoseconds but the signals of the strongly amplified magnetic fields last only a few nanoseconds. Finally, high-fidelity, validated FLASH simulations were crucial to interpreting the results of the experiments. By combining theory, experiments, and simulations, we were able to succeed in demonstrating and characterizing the nonlinear turbulent dynamo mechanism. The effort leveraged (1) the emergence of the Flash Center for Computational Science as a leader in numerical modeling of academic High Energy Density Physics (HEDP) experiments; (2) the FLASH code, which has been developed by the Flash Center and is a leadership-class, multi- physics, publicly available community code that has been applied successfully in many laser-driven experiments; (3) the expertise acquired through pathfinder experiments we performed on the Vulcan laser at the Rutherford-Appleton Laboratory in the UK and on the Omega laser facility; (4) the experience we have gained in designing and interpreting these experiments using validated simulations done with FLASH; and (5) the close collaborations we have established with outstanding experimental groups in the US and abroad. The broader impacts of the project are two-fold: (1) the work furthered the transformation of the academic community’s ability to design and analyze HEDP experiments at large laser facilities that is happening through the availability of FLASH and its widespread adoption by the academic community – a community that previously had had limited access to open and validated hydro/MHD simulation tools for HEDP experiments; and (2) the work trained young scientists to design and interpret HEDP experiments using validated simulations – a critical national need.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Multi-physics Preconditioning for Thermally Activated Batteries

Thermal batteries, also known as molten-salt batteries, are single-use reserve power systems activated by pyrotechnic heat generation, which transitions the solid electrolyte into a molten state. The simulation of these batteries relies on multiphysics modeling to evaluate performance and behavior under various conditions. This paper presents advancements in scalable preconditioning strategies for the Thermally Activated Battery Simulator (TABS) tool, enabling efficient solutions to the coupled electrochemical systems that dominate computational costs in thermal battery simulations. We propose a hierarchical block Gauss-Seidel preconditioner implemented through the Teko package in Trilinos, which effectively addresses the challenges posed by tightly coupled physics, including charge transport, porous flow, and species diffusion. The preconditioner leverages scalable subblock solvers, including smoothed aggregation algebraic multigrid (SA-AMG) methods and domain-decomposition techniques, to achieve robust convergence and parallel scalability. Strong and weak scaling studies demonstrate the solver’s ability to handle problem sizes up to 51.3 million degrees of freedom on 2048 processors, achieving near sub-second setup and solve times for the end-to-end electrochemical solve. These advancements significantly improve the computational efficiency and turnaround time of thermal battery simulations, paving the way for higher-resolution models and enabling the transition from 2D axisymmetric to full 3D simulations.

25 ENERGY STORAGE↗
DOE OSTI2016

Advanced Simulation and Computing FY17 Implementation Plan, Version 0

The Stockpile Stewardship Program (SSP) is an integrated technical program for maintaining the safety, surety, and reliability of the U.S. nuclear stockpile. The SSP uses nuclear test data, computational modeling and simulation, and experimental facilities to advance understanding of nuclear weapons. It includes stockpile surveillance, experimental research, development and engineering programs, and an appropriately scaled production capability to support stockpile requirements. This integrated national program requires the continued use of experimental facilities and programs, and the computational capabilities to support these programs. The Advanced Simulation and Computing Program (ASC) is a cornerstone of the SSP, providing simulation capabilities and computational resources that support annual stockpile assessment and certification, study advanced nuclear weapons design and manufacturing processes, analyze accident scenarios and weapons aging, and provide the tools to enable stockpile Life Extension Programs (LEPs) and the resolution of Significant Finding Investigations (SFIs). This requires a balance of resource, including technical staff, hardware, simulation software, and computer science solutions. ASC is now focused on increasing predictive capabilities in a three-dimensional (3D) simulation environment while maintaining support to the SSP. The program continues to improve its unique tools for solving progressively more difficult stockpile problems (sufficient resolution, dimensionality, and scientific details), and quantifying critical margins and uncertainties. Resolving each issue requires increasingly difficult analyses because the aging process has progressively moved the stockpile further away from the original test base. Where possible, the program also enables the use of high performance computing (HPC) and simulation tools to address broader national security needs, such as foreign nuclear weapon assessments and counter nuclear terrorism.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗