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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,351 records · Page 75

All sources NASA NTRS DOE OSTI

DOE OSTI2015

THE FORMATION OF IRIS DIAGNOSTICS. IV. THE Mg ii TRIPLET LINES AS A NEW DIAGNOSTIC FOR LOWER CHROMOSPHERIC HEATING

A triplet of subordinate lines of Mg ii exists in the region around the h and k lines. In solar spectra these lines are seen mostly in absorption, but in some cases can become emission lines. The aim of this work is to study the formation of this triplet, and investigate any diagnostic value they can bring. Using 3D radiative magnetohydrodynamic simulations of quiet Sun and flaring flux emergence, we synthesize spectra and investigate how spectral features respond to the underlying atmosphere. We find that emission in the lines is rare and is typically caused by a steep temperature increase in the lower chromosphere (above 1500 K, with electron densities above 10{sup 17} m{sup −3}). In both simulations the lines are sensitive to temperature increases taking place at column masses ≳5 · 10{sup −4} g cm{sup −2}. Additional information can also be inferred from the peak-to-wing ratio and shape of the line profiles. Using observations from NASA's Interface Region Imaging Spectrograph we find both absorption and emission line profiles with similar shapes to the synthetic spectra, which suggests that these lines represent a useful diagnostic that complements the Mg ii h and k lines.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2019

3D numerical investigation of a free-burning argon arc with metal electrodes using a novel sheath coupling procedure

In this study, a novel coupling procedure describing the complicated plasma-electrode interaction process has been developed and applied into the 3D finite volume simulation of a direct current tungsten inert gas (TIG) welding system. This is achieved by making the space charge layer (sheath) incorporated into the computation domain and interact with both bulk plasma and cathode through the effective electrical conductivity. Both chemical and thermal nonequilibrium phenomena as well as the self-induced magnetic fields have been taken into consideration by the model to ensure a realistic numerical description of a non-thermal arc. The applicability of this coupling procedure is further improved by calculating the real electric potential, which is capable of accounting for the effects of the complicated drift and diffusion processes. Numerical results of both 100 and 200A discharge currents are presented, field reversal is obtained at near-anode regions in both cases, which is followed by the negative anode sheath potential drop. The region of the strongest electron overpopulation appears at the intersection of plasma fringes and electrode surface. Finally, numerical results of plasma temperature and voltage show good agreement with experimental measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Novel signatures of radiation reaction in electron–laser sidescattering

In this article we investigate novel signatures of radiation reaction via the angular deflection of an electron beam colliding at 90 degrees with an intense laser pulse. Due to the radiation reaction effect, the electrons can be deflected towards the beam axis for plane wave backgrounds, which is not possible in the absence of radiation reaction effects. The magnitude and size of the deflection angle can be controlled by tailoring the laser pulse shapes. The effect is first derived analytically using the Landau–Lifshitz equation, which allows to determine the important scaling behavior with laser intensity and particle energy. We then move on to full scale 3D Monte Carlo simulations to verify the effect is observable with present day laser technology. We investigate the opportunities for an indirect observation of laser depletion in such side scattering scenarios.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Impact of fusion reactor neutronics modeling for transmutation and thermal feedback *

Abstract Fusion neutronics calculations provide important metrics pertinent to fusion device operations, such as tritium breeding ratios (TBRs) and data on heat deposition, material activation, and damage. Because of the high computational burden required to generate a high-fidelity Monte Carlo simulation of a 3D fusion device, various assumptions are made to reduce computational time by simplifying the reactor model or the calculation iteration. This paper explores the impact of fusion neutronics metrics such as the TBR and decay heat of structural materials based on assumptions of material composition in the fusion reactor and temperature modeling of materials. Results show that for compact tokamaks with high power and long operational cycles, the transmutation of structural materials is significant enough to cause a substantial change in the flux spectrum and decrease the TBR by 1.68% after 2 years of full power operation. Additionally, assuming a constant temperature and material density can impact the TBR calculations up to 3%.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

On the use of neural networks for energy reconstruction in high-granularity calorimeters

We contrasted the performance of deep neural networks — Convolutional Neural Network (CNN) and Graph Neural Network (GNN) — to current state of the art energy regression methods in a finely 3D-segmented calorimeter simulated by GEANT4. This comparative benchmark gives us some insight to assess the particular latent signals neural network methods exploit to achieve superior resolution. A CNN trained solely on a pure sample of pions achieved substantial improvement in the energy resolution for both single pions and jets over the conventional approaches. It maintained good performance for electron and photon reconstruction. We also used the Graph Neural Network (GNN) with edge convolution to assess the importance of timing information in the shower development for improved energy reconstruction. We implement a simple simulation based correction to the energy sum derived from the fraction of energy deposited in the electromagnetic shower component. This serves as an approximate dual-readout analogue for our benchmark comparison. Although this study does not include the simulation of detector effects, such as electronic noise, the margin of improvement seems robust enough to suggest these benefits will endure in real-world application. We also find reason to infer that the CNN/GNN methods leverage latent features that concur with our current understanding of the physics of calorimeter measurement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Lagrangian statistics of a shock-driven turbulent dynamo in decaying turbulence

Small-scale fluctuating magnetic fields of order nG are observed in supernova shocks and galaxy clusters, where its amplification is likely caused by the Biermann battery mechanism. However, these fields cannot be amplified further without the turbulent dynamo, which generates magnetic energy through the stretch-twist-fold (STF) mechanism. Thus, we present here novel 3D magnetohydrodynamic (MHD) simulations of a laser-driven shock propagating into a stratified, multiphase medium, to investigate the post-shock turbulent magnetic field amplification via the turbulent dynamo. The configuration used here is currently being tested in the shock tunnel at the National Ignition Facility (NIF). In order to probe the statistical properties of the post-shock turbulent region, we use 384 × 512 × 384 tracers to track its evolution through the Lagrangian framework, thus providing a high-fidelity analysis of the shocked medium. Our simulations indicate that the growth of the magnetic field, which accompanies the near-Saffman kinetic energy decay (E kin ∝ t -1.15 ) without turbulence driving, exhibits slightly different characteristics as compared to periodic box simulations. Seemingly no distinct phases exist in its evolution because the shock passage and time to observe the magnetic field amplification during the turbulence decay are very short (~0.3 of a turbulent turnover time). Yet, the growth rate is still consistent with those expected for compressive (curl-free) turbulence driving in subsonic, compressible turbulence. Phenomenological understanding of the dynamics of the magnetic and velocity fields are also elucidated via Lagrangian frequency spectra, which are consistent with the expected inertial range scalings in the Eulerian–Lagrangian bridge.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Stellar wind impact on early atmospheres around unmagnetized Earth-like planets

ABSTRACT Stellar rotation at early ages plays a crucial role in the survival of primordial atmospheres around Earth-mass exoplanets. Earth-like planets orbiting fast-rotating stars may undergo complete photoevaporation within the first few hundred Myr driven by the enhanced stellar XUV [X-rays and extreme ultraviolet (EUV)] radiation, while planets orbiting slow-rotating stars are expected to experience difficulty in losing their primordial envelopes. Besides the action of stellar radiation, stellar winds induce additional erosion on these primordial atmospheres, altering their morphology, extent, and causing supplementary atmospheric losses. In this paper, we study the impact of activity-dependent stellar winds on primordial atmospheres to evaluate the extent to which the action of these winds can be significant in the whole planetary evolution at early evolutionary stages. We performed 3D magnetohydrodynamical (MHD) simulations of the interaction of photoevaporating atmospheres around unmagnetized Earth-mass planets in the time span between 50 and 500 Myr, analysing the joint evolution of stellar winds and atmospheres for both fast- and slow-rotating stars. Our results reveal substantial changes in the evolution of primordial atmospheres when influenced by fast-rotating stars, with a significant reduction in extent at early ages. In contrast, atmospheres embedded in the stellar winds from slow-rotating stars remain largely unaltered. The interaction of the magnetized stellar winds with the ionized upper atmospheres of these planets allows us to evaluate the formation and evolution of different MHD structures, such as double bow shocks and induced magnetospheres. This work will shed light on the first evolutionary stages of Earth-like exoplanets, which are of crucial relevance in terms of planet habitability.

Astronomy & Astrophysics↗

DOE OSTI2025

Spatiotemporal Energy Cascade in Three-Dimensional Magnetohydrodynamic Turbulence

We present a new scale decomposition method to investigate turbulence in wavenumber-frequency space. Using 3D magnetohydrodynamic turbulence simulations, we show that magnetic fluctuations with time scales longer than the nonlinear time exhibit an inverse cascade toward even smaller frequencies. Low frequency magnetic fluctuations support turbulence, acting as an energy reservoir that is converted into plasma kinetic energy, the latter cascading toward large wavenumbers and frequencies, where it is dissipated. Our results shed new light on the spatiotemporal properties of turbulence, potentially explaining the origin and role of low frequency turbulent fluctuations in the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Efficient Nonthermal Ion and Electron Acceleration Enabled by the Flux-Rope Kink Instability in 3D Nonrelativistic Magnetic Reconnection

The relaxation of field-line tension during magnetic reconnection gives rise to a universal Fermi acceleration process involving the curvature drift of particles. However, the efficiency of this mechanism is limited by the trapping of energetic particles within flux ropes. Using 3D fully kinetic simulations, we demonstrate that the flux-rope kink instability leads to strong field-line chaos in weak-guide-field regimes where the Fermi mechanism is most efficient, thus allowing particles to transport out of flux ropes and undergo further acceleration. As a consequence, both ions and electrons develop clear power-law energy spectra that contain a significant fraction of the released energy. The low-energy bounds are determined by the injection physics, while the high-energy cutoffs are limited only by the system size. Furthermore, these results have strong relevance to observations of nonthermal particle acceleration in space and astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Ideal MHD Limited Electron Temperature in Spherical Tokamaks

It is well documented that the central electron temperature in the national spherical torus experiment (NSTX) remains largely unchanged as the external heating power, and hence the normalized volume averaged plasma pressure β increases [Stutman, Phys. Rev. Lett. 102, 115002 (2009)]. Herein we present a hypothesis that low n, pressure driven ideal magnetohydrodynamic (MHD) instabilities that are nondisruptive, can break magnetic surfaces in the central region and thereby flatten the electron temperature profiles. We demonstrate this mechanism in a 3D resistive MHD simulation of a NSTX discharge. By varying the toroidal magnetic field strength, and/or the heating power, we show that there is a critical value of β, above which the central temperature profile no longer peaks on axis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Plasmoid Formation and Strong Radiative Cooling in a Driven Magnetic Reconnection Experiment

Herein we present the first experimental study of plasmoid formation in a magnetic reconnection layer undergoing rapid radiative cooling, a regime relevant to extreme astrophysical plasmas. Two exploding aluminum wire arrays, driven by the Z machine, generate a reconnection layer (S L ≈ 120) in which the cooling rate far exceeds the hydrodynamic transit rate (τ hydro =τ cool > 100). The reconnection layer generates a transient burst of >1 keV x-ray emission, consistent with the formation and subsequent rapid cooling of the layer. Time-gated x-ray images show fast-moving (up to 50 km s –1 ) hotspots in the layer, consistent with the presence of plasmoids in 3D resistive magnetohydrodynamic simulations. X-ray spectroscopy shows that these hotspots generate the majority of Al K-shell emission (around 1.6 keV) prior to the onset of cooling, and exhibit temperatures (170 eV) much greater than that of the plasma inflows and the rest of the reconnection layer, thus providing insight into the generation of high-energy radiation in radiatively cooled reconnection events.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Quenching through the QCD chiral phase transition

We present a detailed numerical and analytical study of the out-of-equilibrium dynamics of Model G, the dynamical universality class relevant to the chiral phase transition. We perform numerical 3D stochastic (Langevin) simulations of the 𝑂⁡(4) critical point for large lattices in the chiral limit. We quench the system from the high-temperature unbroken phase to the broken phase and study the nonequilibrium dynamics of pion fields. Strikingly, the nonequilibrium evolution of the two-point functions exhibits a regime of growth, a parametrically large enhancement, and a subsequent slow relaxation to equilibrium. We analyze our numerical results using dynamic critical scaling and mean-field theory. The growth of the two-point functions is determined by the nonlinear dynamics of an ideal non-Abelian superfluid, which is a limit of Model G that reflects the broken chiral symmetry. We also relate the nonequilibrium two-point functions to a long-lived parametric enhancement of soft pion yields relative to thermal equilibrium following a quench.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2014

SU-E-J-190: Characterization of Radiation Induced CT Number Changes in Tumor and Normal Lung During Radiation Therapy for Lung Cancer

Purpose: To measure CT number (CTN) changes in tumor and normal lung as a function of radiation therapy (RT) dose during the course of RT delivery for lung cancer using daily IGRT CT images and single respiration phase CT images. Methods: 4D CT acquired during planning simulation and daily 3D CT acquired during daily IGRT for 10 lung cancer cases randomly selected in terms of age, caner type and stage, were analyzed using an in-house developed software tool. All patients were treated in 2 Gy fractions to primary tumors and involved nodal regions. Regions enclosed by a series of isodose surfaces in normal lung were delineated. The obtained contours along with target contours (GTVs) were populated to each singlephase planning CT and daily CT. CTN in term of Hounsfield Unit (HU) of each voxel in these delineated regions were collectively analyzed using histogram, mean, mode and linear correlation. Results: Respiration induced normal lung CTN change, as analyzed from single-phase planning CTs, ranged from 9 to 23 (±2) HU for the patients studied. Normal lung CTN change was as large as 50 (±12) HU over the entire treatment course, was dose and patient dependent and was measurable with dose changes as low as 1.5 Gy. For patients with obvious tumor volume regression, CTN within the GTV drops monotonically as much as 10 (±1) HU during the early fractions with a total dose of 20 Gy delivered. The GTV and CTN reductions are significantly correlated with correlation coefficient >0.95. Conclusion: Significant RT dose induced CTN changes in lung tissue and tumor region can be observed during even the early phase of RT delivery, and may potentially be used for early prediction of radiation response. Single respiration phase CT images have dramatically reduced statistical noise in ROIs, making daily dose response evaluation possible.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2014

WE-G-17A-09: Novel Magnetic Shielding Design for Inline and Perpendicular Integrated 6 MV Linac and 1.0 T MRI Systems

Purpose: The influence of fringe magnetic fields delivered by magnetic resonance imaging (MRI) on the beam generation and transportation in Linac is still a major challenge for the integration of linear accelerator and MRI (Linac-MRI). In this study, we investigated an optimal magnetic shielding design for Linac-MRI and further characterized the beam trajectory in electron gun. Methods: Both inline and perpendicular configurations were analyzed in this study. The configurations, comprising a Linac-MRI with a 100cm SAD and an open 1.0 T superconductive magnet, were simulated by the 3D finite element method (FEM). The steel shielding around the Linac was included in the 3D model, the thickness of which was varied from 1mm to 20mm, and magnetic field maps were acquired with and without additional shielding. The treatment beam trajectory in electron gun was evaluated using OPERA 3d SCALA with and without shielding cases. Results: When Linac was not shielded, the uniformity of diameter sphere volume (DSV) (30cm) was about 5 parts per million (ppm) and the fringe magnetic fields in electron gun were more than 0.3 T. With shielding, the magnetic fields in electron gun were reduced to less than 0.01 T. For the inline configuration, the radial magnetic fields in the Linac were about 0.02T. A cylinder steel shield used (5mm thick) altered the uniformity of DSV to 1000 ppm. For the perpendicular configuration, the Linac transverse magnetic fields were more than 0.3T, which altered the beam trajectory significantly. A 8mm-thick cylinder steel shield surrounding the Linac was used to compensate the output losses of Linac, which shifted the magnetic fields' uniformity of DSV to 400 ppm. Conclusion: For both configurations, the Linac shielding was used to ensure normal operation of the Linac. The effect of magnetic fields on the uniformity of DSV could be modulated by the shimming technique of the MRI magnet. NIH/NIGMS grant U54 GM104944, Lincy Endowed Assistant Professorship.

43 PARTICLE ACCELERATORS↗

DOE OSTI2015

SU-E-T-786: Utility of Gold Wires to Optimize Intensity Modulation Capacity of a Novel Directional Modulated Brachytherapy Tandem Applicator for Image Guided Cervical Cancer Brachytherapy

Purpose: To evaluate the impact of using gold wires to differentially fill various channels on plan quality compared with conventional T&R applicator, inside a novel directional modulated brachytherapy (DMBT) tandem applicator for cervical cancer brachytherapy. Materials and Methods: The novel DMBT tandem applicator has a 5.4-mm diameter MR-compatible tungsten alloy enclosed in a 0.3-mm thick plastic tubing that wraps around the tandem. To modulate the radiation intensity, 6 symmetric peripheral holes of 1.3-mm diameter are grooved along the tungsten alloy rod. These grooved holes are differentially filled with gold wires to generate various degrees of directional beams. For example, three different fill patterns of 1) all void, 2) all filled except the hole containing the 192-Ir source, and 3) two adjacent holes to the 192-Ir source filled were Monte Carlo simulated. The resulting 3D dose distributions were imported into an in-house-coded inverse optimization planning system to generate HDR brachytherapy clinical plans for 19 patient cases. All plans generated were normalized to the same D90 as the clinical plans and D2cc doses of OARs were evaluated. Prescription ranged between 15 and 17.5Gy. Results: In general, the plans in case 1) resulted in the highest D2cc doses for the OARs with 11.65±2.30Gy, 7.47±3.05Gy, and 9.84±2.48Gy for bladder, rectum, and sigmoid, respectively, although the differences were small. For the case 2), D2cc doses were 11.61±2.29Gy, 7.41±3.07Gy, and 9.75±2.45Gy, respectively. And, for the case 3), D2cc doses were 11.60±2.28Gy, 7.41±3.05Gy, and 9.74±2.45Gy, respectively. Difference between 1) and 2) cases were small with the average D2cc difference of <0.64%. Difference between 1) and 3) cases were even smaller with the average D2cc difference of <0.1%. Conclusions: There is a minimal clinical benefit by differentially filling grooved holes in the novel DMBT tandem applicator for image guided cervical cancer brachytherapy.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2022

Effects of stratification and charge cooling on combustion in a gasoline direct-injection compression ignition (GDCI) engine

With the development of low temperature engine combustion strategies, performance of gasoline-type fuels under compression ignition conditions has attracted extensive research interest. Meanwhile, for the sake of co-optimization of engines and fuels for future ground transportation, identification and evaluation of general fuel properties should be a core research priority instead of endless testing of specific fuels. In this study, the roles of fuel octane sensitivity in characterizing the ignition performance of gasoline surrogates have been systematically investigated under typical gasoline direct ignition compression ignition (GDCI) engine conditions using 3D combustion CFD simulation, especially considering the subsequent in-cylinder charge stratification and charge cooling. Two different operating conditions, high boost pressure low boost temperature (beyond-RON) case and low boost pressure high boost temperature (beyond-MON) case, were considered. By comparing with our previous zero-dimensional chemical kinetic study of gasoline surrogates in advanced compression ignition (ACI) engines, the effects of stratification and charge cooling on the combustion processes are investigated. It is found that different fuel octane sensitivities lead to slight difference in equivalence ratio stratification and charge cooling due to differences in volatility. However, fuel reactivity is still the more dominant factor than the stratification and charge cooling effects in determining combustion phasing. The present results help to justify the P-T domain framework for engine autoignition analysis of overlapping pressure-temperature trajectory with ignition delay iso-contour. Here, the results also provide useful guidance to the understanding of GCI combustion process, and to the evaluation of controlling fuel properties and the selection of alternative fuels in GCI engines.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2015

Comparison of the PHISICS/RELAP5-3D ring and block model results for phase I of the OECD/NEA MHTGR-350 benchmark

The PHISICS code system has been under development at INL since 2010. It consists of several modules providing improved coupled core simulation capability: INSTANT (3D nodal transport core calculations), MRTAU (depletion and decay heat generation) and modules performing criticality searches, fuel shuffling and generalized perturbation. Coupling of the PHISICS code suite to the thermal hydraulics system code RELAP5-3D was finalized in 2013, and as part of the verification and validation effort the first phase of the OECD/NEA MHTGR-350 Benchmark has now been completed. The theoretical basis and latest development status of the coupled PHISICS/RELAP5-3D tool are described in more detail in a concurrent paper. This paper provides an overview of the OECD/NEA MHTGR-350 Benchmark and presents the results of Exercises 2 and 3 defined for Phase I. Exercise 2 required the modelling of a stand-alone thermal fluids solution at End of Equilibrium Cycle for the Modular High Temperature Reactor (MHTGR). The RELAP5-3D results of four sub-cases are discussed, consisting of various combinations of coolant bypass flows and material thermophysical properties. Exercise 3 required a coupled neutronics and thermal fluids solution, and the PHISICS/RELAP5-3D code suite was used to calculate the results of two sub-cases. The main focus of the paper is a comparison of results obtained with the traditional RELAP5-3D “ring” model approach against a much more detailed model that include kinetics feedback on individual block level and thermal feedbacks on a triangular sub-mesh. The higher fidelity that can be obtained by this “block” model is illustrated with comparison results on the temperature, power density and flux distributions. Furthermore, it is shown that the ring model leads to significantly lower fuel temperatures (up to 10%) when compared with the higher fidelity block model, and that the additional model development and run-time efforts are worth the gains obtained in the improved spatial temperature and flux distributions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2012

BEAMPATH User Manual

BEAMPATH is a 2D and 3D code for simulation of space charge dominated beam dynamics in linear particle accelerators and beam transport. The program is used for particle-in-cell simulation of axial-symmetric, quadrupole-symmetric and z-uniform beams in a channel containing RF gaps, radio-frequency quadrupoles, multipole lenses, solenoids, bending magnets, accelerating waveguides.

Accelerator Design, Technology, and Operations↗