Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D radiation”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Non-linear MHD modelling of edge localized modes suppression by resonant magnetic perturbations in ITER

Edge localized modes (ELMs) suppression by resonant magnetic perturbations (RMPs) was studied with the non-linear magneto-hydro-dynamic (MHD) code JOREK for the ITER H-mode scenarios at 15 MA, 12.5 MA, 10 MA/5.3 T. In this study, the main aim was to demonstrate that ELMs can be suppressed by RMPs while the divertor 3D footprints of heat and particle fluxes remain within divertor material limits. The unstable peeling–ballooning modes responsible for ELMs without RMPs were modelled first for each scenario using numerically accessible parameters for ITER. Then the stabilization of ELMs by RMPs was modelled with the same parameters. RMP spectra, optimized by the linear MHD MARS-F code, with main toroidal harmonics N = 2, N = 3, N = 4 have been used as boundary conditions of the computational domain of JOREK, including realistic RMP coils, main plasma, scrape off layer (SOL) divertor and realistic first wall. The model includes all relevant plasma flows: toroidal rotation, two fluid diamagnetic effects and neoclassical poloidal friction. With RMPs, the main toroidal harmonic and the non-linearly coupled harmonics remain dominant at the plasma edge, producing saturated modes and a continuous MHD turbulent transport thereby avoiding ELM crashes in all scenarios considered here. The threshold for ELM suppression was found at a maximum RMP coils current of 45 kAt–60 kAt compared to the coils maximum capability of 90 kAt. In the high beta poloidal steady-state 10 MA/5.3 T scenario, a rotating QH-mode without ELMs was observed even without RMPs. In this scenario with RMPs N = 3, N = 4 at 20 kAt maximum current in RMP coils, similar QH-mode behaviour was observed however with dominant edge harmonic corresponding to the main toroidal number of RMPs. The present MHD modelling was limited in time by few tens of ms after RMPs were switched on until the magnetic energy of the modes saturates. As a consequence the thermal energy was still evolving on this time scale, far from the ITER confinement time scale and hence only the form of 3D footprints on the divertor targets can be indicated within this set-up. Also note, that the divertor physics was missing in this model, so realistic values of fluxes are out of reach in this modelling. However the stationary 3D divertor and particle fluxes could be simply extrapolated from these results to the stationary situation considering that a large power fraction should be radiated in the core and SOL and only about 50 MW power is going to the divertor, which is an arbitrary, but reasonable number used here. The 3D footprints with RMPs show the characteristic splitting with the main RMP toroidal symmetry. The maximum radial extension of the footprints typically was ~20 cm in inner divertor and ~40 cm in outer divertor with stationary heat fluxes decreasing further out from the initial strike point from ~5 MW m ₋2 to ~1 MW m ₋2 assuming a total power in the divertor and walls is 50 MW. The heat fluxes remain within the divertor target and baffle areas, however with rather small margin in the outer divertor which could be an issue for the first wall especially in transient regimes when part of the plasma thermal energy is released due to switching on the RMP coils. This fact should be considered when RMPs are applied with a more favorable application before or soon after the L–H transition, although optimization is required to avoid increasing the L–H power threshold with RMPs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Simulation of runaway electron production with CQL3D coupled to NIMROD

Abstract A coupling between two distinctly different codes—one magnetohydrodynamic (MHD) and another kinetic—is achieved and applied for simulation of runaway electron (RE) production. The 3D initial value MHD code NIMROD simulates a DIII-D pure neon shattered pellet injection plasma quench including the propagation and ablation of the fragments, ionization and recombination of the impurities, and the radiated and transported energies. The field data from NIMROD is then used by the bounce-averaged Fokker–Planck Collisional QuasiLinear 3D (CQL3D) kinetic code to simulate the production of REs and their radial transport. The coupling procedure involves mapping of data between different grids and adjustment of the NIMROD toroidal electric field when REs appear. It is shown that without the radial transport, a large RE current is generated, up to 30% of the pre-pellet ohmic current. However, when the radial transport is included in CQL3D, the RE current is reduced to undetectable level, consistent with experiment. Various forms of the radial diffusion are surveyed to determine conditions when the fast electrons would not have time to be accelerated to relativistic energies before they are lost to chamber wall.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Development of an integrated core–edge scenario using the super H-mode

An optimized pedestal regime called the super-H (SH) mode is leveraged to couple a fusion relevant core plasma with a high density scrape-off layer appropriate for realistic reactor power exhaust solutions. Recent DIII-D experiments have expanded the operating space of the SH regime using advanced control algorithms and investigated optimization of impurity seeding, deuterium gas puffing, and 3D magnetic perturbations. Simultaneous real-time control of the pedestal density and radiated power with in-vessel coils and nitrogen seeding enable optimal coupled divertor and pedestal conditions. Four case studies are analysed with varied levels of radiated power in the divertor volume ranging from 0 (no seeding) to 8.5 MW radiated from carbon and nitrogen emission. Plasmas with a 4.5 MW radiated power target establish a radiative mantle, leading to divertor temperatures of ~16 eV while maintaining SH-mode, and with only marginal impact on the pedestal and core performance. Increased levels of N 2 seeding with a 7.5 MW radiated power target facilitate detachment onset and divertor temperatures <5 eV, with no degradation in stored energy and the operational point remaining inside the SH-mode channel for $ > 2.5{\tau }_{\text{E}}$. Finally, a 8.5 MW radiated power target leads to partial detachment, which is so far associated with the loss of access to SH-mode pedestal conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Internal-Gelation Production of Uranium Oxide Sol-Gel Particles for Forensic Applications

The purpose of this project was to develop and demonstrate a novel method for the production of uranium oxide microsphere particles with tunable chemical compositions via a sol-gel process using a 3D-printer setup. These particles can serve several purposes in research and development as a forensic training tool or as standard reference materials. A key component of the project was to demonstrate the ability to control physical and chemical parameters of the particles created. First, we demonstrated the ability to employ an internal gelation sol-gel process to create individual uranium oxide particles. The particles were successfully dispensed using a unique 3D-printing setup onto a substrate to react and then were collected and thermally processed. A series of temperatures for the annealing process was tested on individual samples to investigate the effect on the sol-gel chemical composition and physical integrity. Next, we demonstrated the ability to control matrix composition of the particles by separately incorporating fission product isotopes as well as Np-237 into the sol-gel solution. It was shown by gamma-ray spectroscopy that these matrix elements were successfully retained during the gelation process. We studied the retention of the elements across a series of annealing temperatures. Additionally, we demonstrated the ability to quantitatively control the isotopic composition of the particles by altering the U-237/U-238 ratio to a controlled value. Finally, X-ray diffraction analysis (XRD) was used to investigate the oxidation state of the sol-gel after annealing at different temperatures.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2020

Bremsstrahlung in GRMHD Models of Accreting Black Holes

Here, the role of bremsstrahlung in the emission from hot accretion flows around slowly accreting supermassive black holes is not thoroughly understood. In order to appraise the importance of bremsstrahlung relative to other radiative processes, we compute spectral energy distributions (SEDs) of accretion disks around slowly accreting supermassive black holes including synchrotron radiation, inverse Compton scattering, and bremsstrahlung. We compute SEDs for (i) four axisymmetric radiative general relativistic magnetohydrodynamics (RadGRMHD) simulations of 10 8 M ⊙ black holes with accretion rates between 10 -8 $\dot{M}$ Edd and 10 -5 $\dot{M}$ Edd , (ii) four axisymmetric RadGRMHD simulations of M87* with varying dimensionless spin a$_*$ and black hole mass, and (iii) a 3D GRMHD simulation scaled for Sgr A*. At 10 -8 $\dot{M}$ Edd , most of the luminosity is synchrotron radiation, while at 10 -5 $\dot{M}$ Edd the three radiative processes have similar luminosities. In most models, bremsstrahlung dominates the SED near 512 keV. In the M87* models, bremsstrahlung dominates this part of the SED if a$_*$ = 0.5, but inverse Compton scattering dominates if a$_*$ = 0.9375. Since scattering is more variable than bremsstrahlung, this result suggests that 512 keV variability could be a diagnostic of black hole spin. In the appendix, we compare some bremsstrahlung formulæ found in the literature.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Global Simulations of Gravitational Instability in Protostellar Disks with Full Radiation Transport. I. Stochastic Fragmentation with Optical-depth-dependent Rate and Universal Fragment Mass

Fragmentation in a gravitationally unstable accretion disk can be an important pathway for forming stellar/planetary companions. To characterize quantitatively the condition for and outcome of fragmentation under realistic thermodynamics, we perform global 3D simulations of gravitationally unstable disks at various cooling rates and cooling types, including the first global simulations of gravitational instability that employ full radiation transport. We find that fragmentation is a stochastic process, with the fragment generation rate per disk area p frag showing an exponential dependence on the parameter β ≡ Ω K t cool , where Ω K is the Keplerian rotation frequency, and t cool is the average cooling timescale. Compared to a prescribed constant β, radiative cooling in the optically thin/thick regime makes p frag decrease slower/faster in β; the critical β corresponding to ∼1 fragment per orbit is ≈3, 5, and 2 for constant β, optically thin, and optically thick cooling, respectively. The distribution function of the initial fragment mass is remarkably insensitive to disk thermodynamics. Regardless of cooling rate and optical depth, the typical initial fragment mass is m frag ≈ 40M tot h 3 , with M tot being the total (star+disk) mass and h = H/R being the disk aspect ratio. Applying this result to typical Class 0/I protostellar disks, we find m frag ∼ 20M J , suggesting that fragmentation more likely forms brown dwarfs. Given the finite width of the m frag distribution, forming massive planets is also possible.

gravitational instability↗
DOE OSTI2021

Automatic 3D image registration for nano-resolution chemical mapping using synchrotron spectro-tomography

Nano-resolution synchrotron X-ray spectro-tomography has been demonstrated as a powerful tool for probing the three-dimensional (3D) structural and chemical heterogeneity of a sample. By reconstructing a number of tomographic data sets recorded at different X-ray energy levels, the energy-dependent intensity variation in every given voxel fingerprints the corresponding local chemistry. The resolution and accuracy of this method, however, could be jeopardized by non-ideal experimental conditions, e.g. instability in the hardware system and/or in the sample itself. Herein is presented one such case, in which unanticipated sample deformation severely degrades the data quality. Here, to address this issue, an automatic 3D image registration method is implemented to evaluate and correct this effect. The method allows the redox heterogeneity in partially delithiated Li x Ta 0.3 Mn 0.4 O 2 battery cathode particles to be revealed with significantly improved fidelity.

3D image registration↗
DOE OSTI2024

Monte Carlo Simulation with CAD Interface for Calculation of 3D Maps of Residual Dose (CRADA)

Objective: To develop an easy-to-use software application to predict and mitigate radiation effects in research environment, space instruments, nuclear plants and medical facilities and help nonproliferation and national security efforts. Tech-X will develop standalone software libraries and command-line tools for ( 1) translating CAD into tessellated surfaces and tetrahedral meshes in GDML (for Geant4 and MARS 15), ROOT (for MARS 15) and HDF5 (for compact representation and for the visualization) formats, (2) healing CAD geometries to make them suitable for Monte Carlo simulations; (3) creating uniform and variable Cartesian and cylindrical meshes for detailed scoring; and ( 4) efficient Monte Carlo navigation in CAD geometries. JLAB will finish automation of simulations of residual dose in CAD geometries and integrate Tech-X software into Geant4 and MARS15. Finally, Tech-X will develop a Graphical User Interface to set up and heal CAD geometries, create input files, run and visualize simulations for residual dose. This application will run on local desktops, local and remote clusters and supercomputers and will be made available through public clouds, such as Amazon Web Services.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Autonomous Radiation Mapping and Quantification using an Unmanned Ground Vehicle. Part I - Environment Mapping

When characterizing facilities, advanced autonomous systems are safe, efficient, and cost-effective tools, which can safely deploy state-of-the-art instrumentation without exposing workers to radiation risks. One main focus of the research is the visual mapping of the environment.This paper presents a framework for the mapping aspect that utilizes lidar and SLAM technologies coupled with a monocular camera to create a point cloud of the environment. Upon completion, the present research work aims to provide an autonomous mobile robot the capability to localize itself and aid the user with visualizing the radiation data. Challenges: C++ node needs to properly synchronize image and point cloud data. 3D lidar has resolution of ± 3 cm. Proper camera calibration needed. Transforms need to be precise for proper alignment of image and point-cloud. Results: Voxelated point clouds. Accurate colors from calibrated cameras including details smaller than lidar resolution. Developed Robotics Operating System (ROS) package for ease of transfer to other robot platforms.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Global Simulations of Gravitational Instability in Protostellar Disks with Full Radiation Transport. II. Locality of Gravitoturbulence, Clumpy Spirals, and Implications for Observable Substructure

Spiral perturbations in a gravitationally unstable accretion disk regulate disk evolution through angular momentum transport and heating and provide an observational signature of gravitational instability (GI). We use global 3D simulations to systematically characterize and understand these spiral perturbations. The spiral perturbations and the resulting transport are overall insensitive to the cooling type, with the exception that radiative cooling, especially in the optically thick regime, reduces the amplitude of temperature perturbations. Spiral perturbations are localized around corotation, allowing transport to be approximated by a local α viscosity to zeroth order in aspect ratio (H/R), but only after averaging over multiple orbits in time and/or multiple scale heights in space. Meanwhile, large-amplitude perturbations from strong gravitoturbulence can cause $\mathcal{O}$(α 1/2 ) deviation in the cooling rate of the disk. We develop empirical prescriptions for the angular momentum transport, heating, and cooling in a gravitoturbulent disk that capture the deviation from a viscous, unperturbed disk to first order in H/R and α 1/2 . The spiral perturbations in saturated gravitoturbulence are clumpy, with dense clumps forming through the nonlinear coupling between multiple modes at different m. Observationally, the clumpy gravitoturbulence produced by saturated GI can be mistaken with observational noise or embedded companions, especially under finite resolution. Meanwhile, grand-design spirals with m-fold symmetry may be uncommon among disks in saturated gravitoturbulence, and we speculate that they may instead be a signature of recently triggered or decaying GI.

Gravitational instability↗
DOE OSTI2024

How numerical simulations helped to achieve breakeven on the NIF

The inertial confinement fusion program relies upon detailed simulations with inertial confinement fusion (ICF) codes to design targets and to interpret the experimental results. These simulations treat as much physics from essential principles as is practical, including laser deposition, cross beam energy transfer, x-ray production and transport, nonlocal thermal equilibrium kinetics, thermal transport, hydrodynamic instabilities, thermonuclear burn, and transport of reaction products. Improvements in radiation hydrodynamic code capabilities and vast increases in computing power have enabled more realistic, accurate 3D simulations that treat all known asymmetry sources. We describe how numerical simulations helped to guide the program, assess the impediments to breakeven, and optimize every aspect of target design. A preshot simulation of the first National Ignition Facility experiment that surpassed breakeven predicted an increased yield that matches the experimental result, within the preshot predicted uncertainty, with a target gain of 1.5. We will cover the key developments in Lawrence Livermore National Laboratory ICF codes that enabled these simulations and give specific examples of how they helped to guide the program.

Marinak, M. M. (ORCID:0009000331127696)↗
DOE OSTI2023

PIXEL2020: Tenth International Workshop on Semiconductor Pixel Detectors for Particles and Imaging

This workshop covers various topics related to pixel detector technology. Development and applications will be discussed for charged particle tracking in High Energy Physics, Nuclear Physics, Astrophysics, and X-ray imaging in Astronomy, Biology, Medicine and Material Science. The conference program will also include reports on front and back end electronics, radiation effects, low mass mechanics and construction techniques, and new technologies such as monolithic and 3D integrated detectors.

Sally, Seidel↗
DOE OSTI2020

Effect of Confinement of SIA Cluster Diffusion by Impurities on Radiation Defect Accumulation due to 14MeV neutrons in Tungsten

Impurities acting as traps, reduce the self-interstitial atom (SIA) cluster mobility, thereby increasing recombination and in turn enhancing radiation-resistance of polycrystalline materials. However, at elevated temperatures, depending on SIA-trap (impurity) binding energy, SIA clusters may trap and detrap (detach) numerous times. In tungsten, SIA clusters of all sizes glide one-dimensionally (1D) along their Burgers vector direction. Hence, detrapped SIAs would retrace their path, confining their 1D-glide between traps. However, small SIA clusters can change the direction of their 1D-glide by overcoming a rotation barrier. In addition, impurities or solutes are known to reduce this barrier. A lower rotation barrier can result in detrapped SIA clusters that diffuse in a random 1D direction, effectively leading to 3D (net-3D) diffusion. Here we investigate the effect of two cases of SIA diffusion, namely confined-1D diffusion and net-3D diffusion, on the damage accumulation in polycrystalline tungsten under 14 MeV neutron irradiation using the object kinetic Monte Carlo method. Simulations were performed using cascade debris obtained from molecular dynamics simulations with primary knock-on atom (PKA) energies corresponding to 14 MeV neutrons. In the simulations with net-3D diffusion, SIA clusters of all sizes diffuse in a random 1D direction after the clusters detrap. While for the confined 1D diffusion case, SIA clusters larger than size 5 retain their original direction of 1D-glide. In both types of simulations, trapped vacancies are assumed to be permanently immobilized. Unexpectedly, the damage accumulation was lower with confined 1D diffusion than that with net-3D. We present a systematic comparison of the influence of the two cases of SIA diffusion on the radiation damage accumulation as a function of dose rate, detrapping barrier, and impurity concentration at 1025 K.

Accumulation and recovery of lattice damage, tungs↗
DOE OSTI2021

General relativistic MHD simulations of non-thermal flaring in Sagittarius A*

Sgr A* exhibits regular variability in its multiwavelength emission, including daily X-ray flares and roughly continuous near-infrared (NIR) flickering. The origin of this variability is still ambiguous since both inverse Compton and synchrotron emission are possible radiative mechanisms. The underlying particle distributions are also not well constrained, particularly the non-thermal contribution. In this work, we employ the GPU-accelerated general relativistic magnetohydrodynamics code H-AMR to perform a study of flare flux distributions, including the effect of particle acceleration for the first time in high-resolution 3D simulations of Sgr A*. For the particle acceleration, we use the general relativistic ray-tracing code bhoss to perform the radiative transfer, assuming a hybrid thermal+non-thermal electron energy distribution. We extract ~60 h light curves in the sub-millimetre, NIR and X-ray wavebands, and compare the power spectra and the cumulative flux distributions of the light curves to statistical descriptions for Sgr A* flares. Our results indicate that non-thermal populations of electrons arising from turbulence-driven reconnection in weakly magnetized accretion flows lead to moderate NIR and X-ray flares and reasonably describe the X-ray flux distribution while fulfilling multiwavelength flux constraints. These models exhibit high rms per cent amplitudes, $\gtrsim 150{{\ \rm per\ cent}}$ both in the NIR and the X-rays, with changes in the accretion rate driving the 230 GHz flux variability, in agreement with Sgr A* observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Islands and Page curves in 4d from Type IIB

Variants of the black hole information paradox are studied in Type IIB string theory setups that realize four-dimensional gravity coupled to a bath. The setups are string theory versions of doubly-holographic Karch/Randall brane worlds, with black holes coupled to non-gravitating and gravitating baths. The 10d versions are based on fully backreacted solutions for configurations of D3, D5 and NS5 branes, and admit dual descriptions as N = 4 SYM on a half space and 3d $T$$^{σ}_{ρ}$ [SU(N)] SCFTs. Island contributions to the entanglement entropy of black hole radiation systems are identified through Ryu/Takayanagi surfaces and lead to Page curves. Analogs of the critical angles found in the Karch/Randall models are identified in 10d, as critical parameters in the brane configurations and dual field theories.

79 ASTRONOMY AND ASTROPHYSICS↗