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 541 records · Page 30

CFD Code Calibration and Inlet-Fairing Effects On a 3D Hypersonic Powered-Simulation Model

A three-dimensional (3D) computational study has been performed addressing issues related to the wind tunnel testing of a hypersonic powered-simulation model. The study consisted of three objectives. The first objective was to calibrate a state-of-the-art computational fluid dynamics (CFD) code in its ability to predict hypersonic powered-simulation flows by comparing CFD solutions with experimental surface pressure data. Aftbody lower surface pressures were well predicted, but lower surface wing pressures were less accurately predicted. The second objective was to determine the 3D effects on the aftbody created by fairing over the inlet; this was accomplished by comparing the CFD solutions of two closed-inlet powered configurations with a flowing- inlet powered configuration. Although results at four freestream Mach numbers indicate that the exhaust plume tends to isolate the aftbody surface from most forebody flow- field differences, a smooth inlet fairing provides the least aftbody force and moment variation compared to a flowing inlet. The final objective was to predict and understand the 3D characteristics of exhaust plume development at selected points on a representative flight path. Results showed a dramatic effect of plume expansion onto the wings as the freestream Mach number and corresponding nozzle pressure ratio are increased.

Huebner, Lawrence D.↗

3D Geolocation of Simulated Lightning Sources from Low-Earth Orbit

The recent removal of the Lightning Imaging Sensor from the International Space Station has left an observational gap in lightning detection from low-Earth orbit (LEO). However, new studies have demonstrated the potential for 3D geolocation of lightning sources using orbiting sensors. The Cubespark mission concept aims to take advantage of these developments by deploying a constellation of satellites with radio frequency (RF) sensors and optical imagers to not only map lightning locations, but also to collect bi-spectral flash images. These new capabilities include mapping storm charge structure, flash channel structure, and distinguishing microphysical processes throughout flash development, helping link microphysics and convective processes with overall flash and storm structure around the globe from LEO. In this study, we simulate lightning RF sources in the very high frequency (VHF) band, extrapolate their signals to space-based detection using an improved ionospheric model, and reconstruct their 3D locations using a time-of-arrival (TOA) minimization algorithm. Various constellation configurations, locations, and atmospheric conditions are considered in order to identify and quantify the three main sources of geolocation error: geometric, ionospheric, and instrumental effects. The promising results of this study emphasize the potential of space-based 3D lightning mapping under diverse conditions. 3D resolution is shown to be better than 1-2 km in many cases, enabling new global applications in meteorology and climate sciences. Here we present a selection of these geolocation results as seen from space alongside recent advancements, paving the way for a future generation of LEO lightning mappers.

CubeSpark↗

Particle acceleration and plasma dynamics during magnetic reconnection in the magnetically dominated regime

Magnetic reconnection is thought to be the driver for many explosive phenomena in the universe. The energy release and particle acceleration during reconnection have been proposed as a mechanism for producing high-energy emissions and cosmic rays. We carry out two- and three-dimensional (3D) kinetic simulations to investigate relativistic magnetic reconnection and the associated particle acceleration. The simulations focus on electron–positron plasmas starting with a magnetically dominated, force-free current sheet (σ≡B{sup 2}/(4πn{sub e}m{sub e}c{sup 2})≫1). For this limit, we demonstrate that relativistic reconnection is highly efficient at accelerating particles through a first-order Fermi process accomplished by the curvature drift of particles along the electric field induced by the relativistic flows. This mechanism gives rise to the formation of hard power-law spectra f∝(γ−1){sup −p} and approaches p = 1 for sufficiently large σ and system size. Eventually most of the available magnetic free energy is converted into nonthermal particle kinetic energy. An analytic model is presented to explain the key results and predict a general condition for the formation of power-law distributions. The development of reconnection in these regimes leads to relativistic inflow and outflow speeds and enhanced reconnection rates relative to nonrelativistic regimes. In the 3D simulation, the interplay between secondary kink and tearing instabilities leads to strong magnetic turbulence, but does not significantly change the energy conversion, reconnection rate, or particle acceleration. This study suggests that relativistic reconnection sites are strong sources of nonthermal particles, which may have important implications for a variety of high-energy astrophysical problems.

79 ASTRONOMY AND ASTROPHYSICS↗

Interactive 3D Models and Simulations for Nuclear Security Education, Training, and Analysis.

By providing examples of products that have been produced in the past, it is the hopes of the authors that the audience will have a more thorough understanding of 3D modeling tools, potential applications, and capabilities that they can provide. Truly the applications and capabilities of these types of tools are only limited by one’s imagination. The future of three-dimensional models lies in the expansion into the world of virtual reality where one will experience a fully immersive first-person environment. The use of headsets and hand tools will allow students and instructors to have a more thorough spatial understanding of facilities and scenarios that they will encounter in the real world.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Quantitative Comparison of Ternary Eutectic Phase-Field Simulations with Analytical 3D Jackson–Hunt Approaches

For the analytical description of the relationship between undercoolings, lamellar spacings and growth velocities during the directional solidification of ternary eutectics in 2D and 3D, different extensions based on the theory of Jackson and Hunt are reported in the literature. Besides analytical approaches, the phase-field method has been established to study the spatially complex microstructure evolution during the solidification of eutectic alloys. The understanding of the fundamental mechanisms controlling the morphology development in multiphase, multicomponent systems is of high interest. For this purpose, a comparison is made between the analytical extensions and three-dimensional phase-field simulations of directional solidification in an ideal ternary eutectic system. Based on the observed accordance in two-dimensional validation cases, the experimentally reported, inherently three-dimensional chain-like pattern is investigated in extensive simulation studies. The results are quantitatively compared with the analytical results reported in the literature, and with a newly derived approach which uses equal undercoolings. A good accordance of the undercooling–spacing characteristics between simulations and the analytical Jackson–Hunt apporaches are found. The results show that the applied phase-field model, which is based on the Grand potential approach, is able to describe the analytically predicted relationship between the undercooling and the lamellar arrangements during the directional solidification of a ternary eutectic system in 3D.

36 MATERIALS SCIENCE↗

Towards an understanding of the resolution dependence of Core-Collapse Supernova simulations

ABSTRACT Using our new state-of-the-art core-collapse supernova (CCSN) code Fornax, we explore the dependence upon spatial resolution of the outcome and character of three-dimensional (3D) supernova simulations. For the same 19 M⊙ progenitor star, energy and radial binning, neutrino microphysics, and nuclear equation of state, changing only the number of angular bins in the θ and ϕ directions, we witness that our lowest resolution 3D simulation does not explode. However, when jumping progressively up in resolution by factors of two in each angular direction on our spherical-polar grid, models then explode, and explode slightly more vigorously with increasing resolution. This suggests that there can be a qualitative dependence of the outcome of 3D CCSN simulations upon spatial resolution. The critical aspect of higher spatial resolution is the adequate capturing of the physics of neutrino-driven turbulence, in particular its Reynolds stress. The greater numerical viscosity of lower resolution simulations results in greater drag on the turbulent eddies that embody turbulent stress, and, hence, in a diminution of their vigor. Turbulent stress not only pushes the temporarily stalled shock further out, but bootstraps a concomitant increase in the deposited neutrino power. Both effects together lie at the core of the resolution dependence we observe.

79 ASTRONOMY AND ASTROPHYSICS↗

Supo Thermal Model Development IV

This report describes the continuation of the Computational Fluid Dynamics (CFD) model of the Supo cooling system described in the reports, Supo Thermal Model Development, II, and III. In the previous reports, the solution natural convection of Supo was simulated using a 2D, axisymmetric model that altered the actual geometry. A 3D geometry of Supo will more accurately estimate the natural convection heat transfer coefficient (HTC) of the system since Supo contains several non-symmetric internal components that can only be modeled in 3D. This, in turn, will improve the natural convection HTC results for comparison purposes with past operational data and previous 2D CFD results. Unfortunately, the complexities of performing multiphase 3D CFD simulations has proven too large of a task for our current capabilities. Instead, this report will focus on the 3D model development for the 2.39 kW/L (31.3 kW) case and limitations and challenges associated with the 3D simulations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Plane shear flows of frictionless spheres: Kinetic theory and 3D soft-sphere discrete element method simulations

We use existing 3D Discrete Element simulations of simple shear flows of spheres to evaluate the radial distribution function at contact that enables kinetic theory to correctly predict the pressure and the shear stress, for different values of the collisional coefficient of restitution. Then, we perform 3D Discrete Element simulations of plane flows of frictionless, inelastic spheres, sheared between walls made bumpy by gluing particles in a regular array, at fixed average volume fraction and distance between the walls. The results of the numerical simulations are used to derive boundary conditions appropriated in the cases of large and small bumpiness. Those boundary conditions are, then, employed to numerically integrate the differential equations of Extended Kinetic Theory, where the breaking of the molecular chaos assumption at volume fraction larger than 0.49 is taken into account in the expression of the dissipation rate. We show that the Extended Kinetic Theory is in very good agreement with the numerical simulations, even for coefficients of restitution as low as 0.50. When the bumpiness is increased, we observe that some of the flowing particles are stuck in the gaps between the wall spheres. As a consequence, the walls are more dissipative than expected, and the flows resemble simple shear flows, i.e., flows of rather constant volume fraction and granular temperature.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Global seamless tidal simulation using a 3D unstructured-grid model (SCHISM v5.10.0)

We present a new 3D unstructured-grid global ocean model to study both tidal and nontidal processes, with a focus on the total water elevation. Unlike existing global ocean models, the new model resolves estuaries and rivers down to ~8 m without the need for grid nesting. The model is validated with both satellite and in situ observations for elevation, temperature, and salinity. Tidal elevation solutions have a mean complex root-mean-square error (RMSE) of 4.2 cm for M2 and 5.4 cm for all five major constituents in the deep ocean. The RMSEs for the other four constituents, S2, N2, K1, and O1, are, respectively, 2.05, 0.93, 2.08, and 1.34 cm). The nontidal residual assessed by a tide gauge dataset (GESLA) has a mean RMSE of 7 cm. For the first time ever, we demonstrate the potential for seamless simulation on a single mesh from the global ocean into several estuaries along the US West Coast. The model is able to accurately capture the total elevation, even at some upstream stations. The model can therefore potentially serve as the backbone of a global tide surge and compound flooding forecasting framework.

54 ENVIRONMENTAL SCIENCES↗

SU-E-T-330: To Analyze the Calculation Error of Live Dose-Volume Indices Applying 4D-CT in Radiotherapy for PTVs Within the Liver Completely

Purpose: To study the variation rule of normal liver dose-volume indices calculation for the liver malignancy patients whose plan target volumes were in the liver completely in all breath phases. Methods: Ten patients who accepted radiotherapy for malignant tumor were selected in our study. All patients underwent 4D-CT simulation and 3D-CT simulation in free breathing(FB). 4D-CT was sorted into 10 series CT images according to breath phase, named CT0, CT10 to CT90, respectively. And GTVs were contoured on different CT series, and the individual target volume(ITV) was obtained by merging 10 GTVs from 4D-CT. The PTVs were obtained from ITV applying margins. The PTVs were not beyond the boundary of liver in all breath phase observed by dynamic 4D-CT. The radiotherapy plans were designed and irradiation dose was calculated on 3D-CT images, and the livers were contoured on different series CT images and mapped to 3D-CT images applying rigid registration. To compare the dose-volume difference of livers based on distinct CT images. Results: (1)The liver volumes were similar on 4D-CT and 3D-CT images(CTFB 1485±500cm{sup 3}, CT0 1413±377cm{sup 3}, CT10 1409±396cm{sup 3}, CT20 1419±418cm{sup 3},CT30 1485±500cm{sup 3}, CT40 1438±392cm{sup 3}, CT50 1437±404cm{sup 3}, CT60 1439±409cm{sup 3}, CT70 1408±393cm{sup 3}, CT80 1384±397cm{sup 3}, CT90 1398±397cm{sup 3}; F=0.064,p=1.00) (2) The PTVs volume were 30.17±14.62cm{sup 3};(3) The mean dose and V5 to V10 of liver were similar among 4D-CT different series CT images(p>0.05), and the indices varied less than ±4% refer to liver on CT50. Conclusion: The calculation affection of liver dose-volume indices induced by breath motion were not significant for the PTV within liver completely as estimation before. And more objective prediction indices for radiation induced l.

60 APPLIED LIFE SCIENCES↗

A Physics-Constrained Deep Learning Model for Simulating Multiphase Flow in 3D Heterogeneous Porous Media

Physics-based simulators for multiphase flow in porous media emulate nonlinear processes with coupled physics, and usually require extensive computational resources for software development, maintenance and simulation execution. As a result, a huge demand exists for fast modeling of coupled processes in a wide range of subsurface applications including geological sequestration, hydrocarbon recovery and geothermal energy extraction. In this work, an efficient physics-constrained deep learning model is developed for solving multiphase flow in 3-Dimensional (3D) heterogeneous porous media. The model fully leverages the spatial topology predictive capability of convolutional neural networks, specifically U-Net with successive contracting and expansive steps, and is coupled with an efficient continuity-based smoother to predict flow responses that need spatial continuity. Furthermore, the transient regions are penalized to steer the training process such that the model can accurately capture flow in these regions. The model takes inputs including properties of porous media, fluid properties and well controls, and predicts the temporal-spatial evolution of the state variables (pressure and saturation). While maintaining the continuity of fluid flow, the 3D spatial domain is decomposed into 2D images for reducing training cost, and the decomposition results in an increased number of training data samples and better training efficiency. Additionally, a surrogate model is separately constructed as a postprocessor to calculate well flow rate based on the predictions of state variables from the deep learning model. We use the example of CO 2 injection into saline aquifers, and apply the physics-constrained deep learning model that is trained from physics-based simulation data and emulates the physics process. The model performs prediction with a speedup of ~ 1400 times compared to physics-based simulations, and the average temporal errors of predicted pressure and saturation plumes are 0.27% and 0.099% respectively. Furthermore, water production rate is efficiently predicted by a surrogate model for well flow rate, with a mean error less than 5%. Therefore, with its unique scheme to cope with the fidelity in fluid flow in porous media, the physics-constrained deep learning model can become an efficient predictive model for computationally demanding inverse problems or other coupled processes.

58 GEOSCIENCES↗

Seismic H2: Version 1.0

Seismic-H2 is an integrated software package for geological hydrogen reservoir simulation, optimization, and leakage monitoring. The package includes multiple components: (1) code used for modeling seismic wave propagation in 3D heterogeneous elastic media based on finite-difference method to support detection of geological hydrogen storage reservoir leakage; (2) 3D reservoir simulations of leaks from an underground reservoir and 3D simulations of saline aquifers and depleted gas reservoirs; (3) seismic monitoring costs of passive and active seismic monitoring required for UHS; (4) rock physics calculations and interpolations for converting the reservoir simulations from part (2) into the elastic media models in part (1); (5) pre-processing seismic data; and lastly (6), a GUI interface that combines these different components.

Creasy, Neala↗

Grand potential sintering simulations of doped UO 2 accident-tolerant fuel concepts

Chromium-doped UO is a widely-studied near-term deployable accident-tolerant fuel concept because it results in a dense, large-grain structure that increases the fuel resistance to densification, swelling, and fission gas release. A new charged-interstitial mechanism was recently proposed to describe the behavior of dopants like chromium in sintered UO 2 . Based on that mechanism, manganese was suggested as an even stronger dopant than chromium. In the current work we use mesoscale sintering simulations in an effort to validate the new mechanism. We compare the relative behavior of Cr-doped and undoped UO 2 against experimental data in the literature. We also make predictions of the relative behavior of Mn-doped UO 2 . Simulations are done using the phase field-based grand potential sintering model. Dopants have two effects on sintered UO 2 . They marginally increase the densification rate and greatly increase the average grain size. Both of these effects are individually tested in small-scale simulations. Then large 3D sintering simulations are used to test the combined behavior of both effects. The results for the densification rate simulations are consistent with experiments for Cr-doped fuels. However, we note the grain growth rates are lower than what is found in the literature for Cr-doped fuels.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Lower-hybrid drift waves and their interaction with plasmas in a 3D symmetric reconnection simulation with zero guide field

We investigate lower-hybrid drift waves (LHDW) in symmetric magnetic reconnection with zero guide field using three-dimensional particle-in-cell simulations. The long-wavelength mode with kρiρe∼1 develops in the bifurcated electron current layer around the X-line within the width of the electron meandering motion from the mid-plane, where ρi(e) is the ion (electron) gyroradius. The short-wavelength mode with kρe∼1 develops in the separatrix region downstream of the electron outflow jet, producing electron vortices in the background flow frame. Electrons follow the E × B drift with corrections from the diamagnetic drift and are heated inside the vortices with diverging electric fields. In the vortices, ions have comparable E × B and inertia drifts, which together mostly cancel the diamagnetic drift. Toward the center of diverging field vortices, ions are decelerated, leading to a decrease in the perpendicular temperature, while the loss of low-energy ions results in an increase in the parallel temperature. Parallel electric fields exist as a combination of the LHDW wave field projected to the magnetic field direction and the penetration of whistler waves that are mainly outside of the LHDW layer. The magnetic flux tube is twisted in the vortices. The twist may potentially lead to slippage reconnection, as indicated by the non-uniform parallel potential variation across field lines, while the periodic variations of the twisting directions are a limiting factor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulating MADMAX in 3D: Requirements for Dielectric Axion Haloscopes

We present 3D calculations for dielectric haloscopes such as the currently envisioned MADMAX experiment. For ideal systems with perfectly flat, parallel and isotropic dielectric disks of finite diameter, we find that a geometrical form factor reduces the emitted power by up to 30% compared to earlier 1D calculations. We derive the emitted beam shape, which is important for antenna design. We show that realistic dark matter axion velocities of 10 -3 c and inhomogeneities of the external magnetic field at the scale of 10% have negligible impact on the sensitivity of MADMAX. We investigate design requirements for which the emitted power changes by less than 20% for a benchmark boost factor with a bandwidth of 50 MHz at 22 GHz, corresponding to an axion mass of 90 μeV. We find that the maximum allowed disk tilt is 100 μm divided by the disk diameter, the required disk planarity is 20 μm (min-to-max) or better, and the maximum allowed surface roughness is 100 μm (min-to-max). We show how using tiled dielectric disks glued together from multiple smaller patches can affect the beam shape and antenna coupling.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

MERIDIONAL CIRCULATION DYNAMICS FROM 3D MAGNETOHYDRODYNAMIC GLOBAL SIMULATIONS OF SOLAR CONVECTION

The form of solar meridional circulation is a very important ingredient for mean field flux transport dynamo models. However, a shroud of mystery still surrounds this large-scale flow, given that its measurement using current helioseismic techniques is challenging. In this work, we use results from three-dimensional global simulations of solar convection to infer the dynamical behavior of the established meridional circulation. We make a direct comparison between the meridional circulation that arises in these simulations and the latest observations. Based on our results, we argue that there should be an equatorward flow at the base of the convection zone at mid-latitudes, below the current maximum depth helioseismic measures can probe (0.75 R{sub ⊙}). We also provide physical arguments to justify this behavior. The simulations indicate that the meridional circulation undergoes substantial changes in morphology as the magnetic cycle unfolds. We close by discussing the importance of these dynamical changes for current methods of observation which involve long averaging periods of helioseismic data. Also noteworthy is the fact that these topological changes indicate a rich interaction between magnetic fields and plasma flows, which challenges the ubiquitous kinematic approach used in the vast majority of mean field dynamo simulations.

79 ASTRONOMY AND ASTROPHYSICS↗