Search NASA⌕ Search

SEARCH · Search NASA

Results for “3D Simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,369 records · Page 76

A New 2D-Transport, 1D-Diffusion Approximation of the Boltzmann Transport equation

The work performed in this project consisted of the derivation, implementation, and testing of a new, computationally advantageous approximation to the 3D Boltzmann transport equation. The solution of the Boltzmann equation is the neutron flux in nuclear reactor cores and shields, but solving this equation is difficult and costly. The new “2D/1D” approximation takes advantage of a special geometric feature of typical 3D reactors to approximate the neutron transport physics in a specific (ax- ial) direction, but not in the other two (radial) directions. The resulting equation is much less expensive to solve computationally, and its solutions are expected to be sufficiently accurate for many practical problems. In this project we formulated the new equation, discretized it using standard methods, developed a stable iteration scheme for solving the equation, implemented the new numerical scheme in the MPACT code, and tested the method on several realistic problems. All the hoped- for features of this new approximation were seen. For large, difficult problems, the resulting 2D/1D solution is highly accurate, and is calculated about 100 times faster than a 3D discrete ordinates simulation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Massively Parallel Solver for the Mechanical Harmonic Analysis of Accelerator Cavities

ACE3P is a 3D massively parallel simulation suite that developed at SLAC National Accelerator Laboratory that can perform coupled electromagnetic, thermal and mechanical study. Effectively utilizing supercomputer resources, ACE3P has become a key simulation tool for particle accelerator R and D. A new frequency domain solver to perform mechanical harmonic response analysis of accelerator components is developed within the existing parallel framework. This solver is designed to determine the frequency response of the mechanical system to external harmonic excitations for time-efficient accurate analysis of the large-scale problems. Coupled with the ACE3P electromagnetic modules, this capability complements a set of multi-physics tools for a comprehensive study of microphonics in superconducting accelerating cavities in order to understand the RF response and feedback requirements for the operational reliability of a particle accelerator. (auth)

43 PARTICLE ACCELERATORS↗

SPE-5 Ground-Motion Prediction at Far-Field Geophone and Accelerometer Array Sites and SPE-5 Moment and Corner-Frequency Prediction

This report offers predictions for the SPE-5 ground-motion and accelerometer array sites. These predictions pertain to the waveform and spectral amplitude at certain geophone sites using Denny&Johnson source model and a source model derived from SPE data; waveform, peak velocity and peak acceleration at accelerometer sites using the SPE source model and the finite-difference simulation with LLNL 3D velocity model; and the SPE-5 moment and corner frequency.

47 OTHER INSTRUMENTATION↗

Final Report for Large Scale Laboratory Experiments on Kinetic Magnetic Reconnection

Magnetic reconnection plays a fundamental role in nearly all magnetized plasmas as it enables magnetic energy to be converted into high speed flows, and thermal energy. It allows the magnetic field lines to change topology in collisionless plasmas, and it thereby controls the spatial and temporal evolution of explosive phenomena such as solar flares, coronal mass ejections, and internal disruptions in fusion devices. Because of its strong impact on the plasma dynamics, understanding the physics of magnetic reconnection is one of the most important topics in plasma physics. New frontiers have now emerged in reconnection research with an emphasis on particle energization, reconnection with multiple X-lines in large systems, and the role of kinetic effects in collisionless reconnection. These research topics are closely linked and large scale computer simulations show that 3D effects are likely important to explain the detailed observations of reconnection in solar, interplanetary and magnetospheric plasmas. For experiments to remain relevant in this maturing field of research and help push forward the new frontiers, new devices are needed which access the regimes of large system size and low normalized collisionallity. To address this need, the project the new Terrestrial Reconnection EXperiment (TREX), which is spatially the largest dedicated reconnection experiment constructed to date. With the use of the 3m diameter spherical vacuum vessel of the Madison Plasma Dynamo Experiment (MPDX), the TREX is highly leveraged by earlier NSF investments. MPDX was constructed with funds from NSF's MRI program and is now being operated as a multi-investigator facility, which we name the Wisconsin Plasma Physics Laboratory (WiPPL). The project was carried out mainly by two UW-Madison graduate students, for reconnection in parameter regimes and configurations in TREX not attainable by previous experiments. The implementation of TREX within WiPPL renders the upgraded facility unique in its capabilities to extend the frontiers of experimental reconnection research, including access to reconnection in the collisionless regime at large system sizes and compatible with three-dimensional configurations. The theoretical teams in reconnection research have been given full access to the experimental results. Our experimental results havel allowed us to test existing theoretical models, provide resolution to controversy, and stimulate new understanding of magnetic reconnection. The WiPPAL, with its large basic plasma physics device, is largely operated by graduate and undergraduate students. It has proven an ideal facility for students to become immersed in state of the art research while learning the fundamental principles of plasma physics. The program has also proven effective in attracting new outstanding students into plasma physics and providing them with exciting learning experiences. While our research program and experimental results successfully compete with those of professional scientists at national labs, the most important legacy of our research is, perhaps, the broad training in experimental plasma physics of the large number of students affiliated with our program.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Formation and Diagnostics of Collisionless Shocks in Laser-Produced Plasmas (Final Report)

Collisions of supersonic flows occur frequently in astrophysics, and the resulting shock waves are responsible for the properties of many astrophysical phenomena, such as supernova remnants, Gamma Ray Bursts and jets from Active Galactic Nuclei. Because of the low density of astrophysical plasmas, the mean free path due to Coulomb collisions is typically very large; therefore, most shock waves in astrophysics are “collisionless” and form due to plasma instabilities and self-generated magnetic fields. Laboratory experiments at DOE high energy density (HED) facilities can achieve the conditions necessary for the formation of collisionless shocks and will provide a unique avenue for studying the nonlinear physics of shock waves. This project investigates the parameters and experimental configurations needed for the creation of collisionless shocks in the laboratory experiments. This study is being performed with state of the art 3D kinetic plasma simulations. Recent advances in ab-initio particle-in-cell (PIC) simulations of astrophysical shocks enable us to simulate shock formation on timescales and in geometries relevant to the laboratory experiments. The main goal of the proposal is to ascertain which signatures of shock physics can be extracted from the realistic experiments, accounting for limitations imposed by finite spatial and time scales in the lab. To aid the interpretation of current and future experimental data, PIC simulation data have been analyzed using the same diagnostic techniques as used in the laboratory, such as proton radiography and neutronics, to unambiguously predict the experimental signatures of relevant collisionless shock phenomena.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Representative-cell-based Modeling of HTTF in SAM

Although the current nuclear power market is primarily occupied by the light water reactors (LWRs), the TRISO-fueled helium-cooled graphite-moderated high temperature gas reactors (HTGRs) are drawing growing attention as the nuclear power industry marches towards more advanced systems. Compared to that of conventional LWRs, the thermal hydraulics of HTGR cores show extra complexity from their multi-scale heat transfer mechanisms with varying importance during different operation modes or transient stages. The modeling of HTGR cores for the system analysis purpose therefore faces the challenge of reaching reasonable fidelity and accuracy while maintaining sufficient simplicity. Previous development and validation efforts have demonstrated that a “2-D ring model” with reasonable performance can be implemented in the System Analysis Module (SAM) for prismatic HTGR cores. In the current project, an alternative modeling methodology based on “representative cells” is proposed for the typical prismatic core of HTGRs. Different from the previous ring model approach, the proposed methodology first separately models and then combines the small- and large-scale thermal conductions, by connecting representative local heat transfer units (cells) with effective core- wise thermal resistance. A modeling practice for an integral high-temperature test facility (HTTF) elaborates the modeling details. Steady-state validation of the resultant candidate model is performed against a higher-resolution benchmark from a 3D-1D coupled simulation, which shows satisfactory prediction with reasonably captured global parameters and well-represented temperature fields. A postulated pressurized conduction cooldown (PCC) is also simulated and analyzed, demonstrating the model’s capability of transient prediction with physically captured phenomena resolved in both small and large temporal and spatial scales. In general, the current work achieves a preliminary success in proposing a methodology using representative cells to model prismatic HTGR cores in SAM. Future efforts are envisioned with more validation activities and with potential modeling improvements to eventually achieve the high confidence on a high-fidelity robust modeling methodology with reasonable accuracy.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

THE INSTABILITY AND NON-EXISTENCE OF MULTI-STRANDED LOOPS WHEN DRIVEN BY TRANSVERSE WAVES

In recent years, omni-present transverse waves have been observed in all layers of the solar atmosphere. Coronal loops are often modeled as a collection of individual strands in order to explain their thermal behavior and appearance. We perform three-dimensional (3D) ideal magnetohydrodynamics simulations to study the effect of a continuous small amplitude transverse footpoint driving on the internal structure of a coronal loop composed of strands. The output is also converted into synthetic images, corresponding to the AIA 171 and 193 Å passbands, using FoMo. We show that the multi-stranded loop ceases to exist in the traditional sense of the word, because the plasma is efficiently mixed perpendicularly to the magnetic field, with the Kelvin–Helmholtz instability acting as the main mechanism. The final product of our simulation is a mixed loop with density structures on a large range of scales, resembling a power-law. Thus, multi-stranded loops are unstable to driving by transverse waves, and this raises strong doubts on the usability and applicability of coronal loop models consisting of independent strands.

79 ASTRONOMY AND ASTROPHYSICS↗

AGE dependence of wind properties for solar-type stars: a 3D study

Young and rapidly rotating stars are known for intense, dynamo-generated magnetic fields. Spectropolarimetric observations of those stars in precisely aged clusters are key input for gyrochronology and magnetochronology. We use Zeeman Doppler imaging maps of several young K-type stars of similar mass and radius but with various ages and rotational periods to perform three-dimensional (3D) numerical MHD simulations of their coronae and follow the evolution of their magnetic properties with age. Those simulations yield the coronal structure as well as the instant torque exerted by the magnetized, rotating wind on the star. As stars get older, we find that the angular momentum loss decreases with Ω{sub ⋆}{sup 3}, which is the reason for the convergence on the Skumanich law. For the youngest stars of our sample, the angular momentum loss shows signs of saturation around 8Ω{sub ⊙}, which is a common value used in spin evolution models for K-type stars. We compare these results to semianalytical models and existing braking laws. We observe a complex wind-speed distribution for the youngest stars with slow, intermediate, and fast wind components, which are the result of interaction with intense and nonaxisymmetric magnetic fields. Consequently, in our simulations, the stellar wind structure in the equatorial plane of young stars varies significantly from a solar configuration, delivering insight about the past of the solar system interplanetary medium.

79 ASTRONOMY AND ASTROPHYSICS↗

Time-resolved Rotational Velocities in the Upper Atmosphere of WASP-33 b

While steady empirical progress has been made in understanding the structure and composition of hot-planet atmospheres, direct measurements of velocity signatures, including winds, rotation, and jets, have lagged behind. Quantifying atmospheric dynamics of hot planets is critical for a complete understanding of their atmospheres, and such measurements may even illuminate other planetary properties, such as magnetic field strengths. In this manuscript we present the first detection of the Balmer lines Hα and Hβ in the atmosphere of the ultra-hot Jupiter WASP-33 b. Using atmospheric models that include the effects of atmospheric dynamics, we show that the shape of the average Balmer line transmission spectrum is consistent with rotational velocities in the planet’s thermosphere of v{sub rot}=10.1{sub −1.0}{sup +0.8} km s{sup −1}. We also measure a low-significance day- to nightside velocity shift of −4.6{sub −3.4}{sup +3.4} km s{sup −1} in the transmission spectrum, which is naturally explained by a global wind across the planet’s terminator. In a separate analysis the time-resolved velocity centroids of individual transmission spectra show unambiguous evidence of rotation, with a best-fit velocity of 10.0{sub −2.0}{sup +2.4} km s{sup −1}, consistent with the value of v {sub rot} derived from the shape of the average Balmer line transmission spectrum. Our observations and analysis confirm the power of time-resolved transmission spectra with a high signal-to-noise ratio to measure the velocity structures in exoplanet atmospheres. The high rotational and wind velocities we measure highlight the need for more detailed 3D global climate simulations of the rarefied upper atmospheres of ultra-hot gas giants.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Bombs and Flares at the Surface and Lower Atmosphere of the Sun

A spectacular manifestation of solar activity is the appearance of transient brightenings in the far wings of the H α line, known as Ellerman bombs (EBs). Recent observations obtained by the Interface Region Imaging Spectrograph have revealed another type of plasma “bombs” (UV bursts) with high temperatures of perhaps up to 8 × 10{sup 4} K within the cooler lower solar atmosphere. Realistic numerical modeling showing such events is needed to explain their nature. Here, we report on 3D radiative magnetohydrodynamic simulations of magnetic flux emergence in the solar atmosphere. We find that ubiquitous reconnection between emerging bipolar magnetic fields can trigger EBs in the photosphere, UV bursts in the mid/low chromosphere and small (nano-/micro-) flares (10{sup 6} K) in the upper chromosphere. These results provide new insights into the emergence and build up of the coronal magnetic field and the dynamics and heating of the solar surface and lower atmosphere.

79 ASTRONOMY AND ASTROPHYSICS↗

Feedback by AGN Jets and Wide-angle Winds on a Galactic Scale

To investigate the differences in mechanical feedback from radio-loud and radio-quiet active galactic nuclei on the host galaxy, we perform 3D AMR hydrodynamic simulations of wide-angle, radio-quiet winds with different inclinations on a single, massive, gas-rich disk galaxy at a redshift of 2–3. We compare our results to hydrodynamic simulations of the same galaxy but with a jet. The jet has an inclination of 0° (perpendicular to the galactic plane), and the winds have inclinations of 0°, 45°, and 90°. We analyze the impact on the host’s gas, star formation, and circumgalactic medium. We find that jet feedback is energy-driven and wind feedback is momentum-driven. In all the simulations, the jet or wind creates a cavity mostly devoid of dense gas in the nuclear region where star formation is then quenched, but we find strong positive feedback in all the simulations at radii greater than 3 kpc. All four simulations have similar SFRs and stellar velocities with large radial and vertical components. However, the wind at an inclination of 90° creates the highest density regions through ram pressure and generates the highest rates of star formation due to its ongoing strong interaction with the dense gas of the galactic plane. With increased wind inclination, we find greater asymmetry in gas distribution and resulting star formation. Our model generates an expanding ring of triggered star formation with typical velocities of the order of 1/3 of the circular velocity, superimposed on the older stellar population. This should result in a potentially detectable blue asymmetry in stellar absorption features at kiloparsec scales.

79 ASTRONOMY AND ASTROPHYSICS↗

Precipitation and Release of Solar Energetic Particles from the Solar Coronal Magnetic Field

Most solar energetic particles (SEPs) are produced in the corona. They propagate through complex coronal magnetic fields subject to scattering and diffusion across the averaged field lines by turbulence. We examine the behaviors of particle transport using a stochastic 3D focused transport simulation in a potential field source surface model of coronal magnetic field. The model is applied to an SEP event on 2010 February 7. We study three scenarios of particle injection at (i) the compact solar flare site, (ii) the coronal mass ejection (CME) shock, and (iii) the EUV wave near the surface. The majority of particles injected on open field lines are able to escape the corona. We found that none of our models can explain the observations of wide longitudinal SEP spread without perpendicular diffusion. If the perpendicular diffusion is about 10% of what is derived from the random walk of field lines at the rate of supergranular diffusion, particles injected at the compact solar flare site can spread to a wide range of longitude and latitude, very similar to the behavior of particles injected at a large CME shock. Stronger pitch-angle scattering results in a little more lateral spread by holding the particles in the corona for longer periods of time. Some injected particles eventually end up precipitating onto the solar surface. Even with a very small perpendicular diffusion, the pattern of the particle precipitation can be quite complicated depending on the detailed small-scale coronal magnetic field structures, which could be seen with future sensitive gamma-ray telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗

Gamma Rays from Kilonova: A Potential Probe of r-process Nucleosynthesis

The mergers of compact binaries with at least one neutron star component are the potential leading sites of the production and ejection of r-process elements. Discoveries of galactic binary pulsars, short gamma-ray bursts, and gravitational-wave detections have all been constraining the rate of these events, while the gravitational wave plus broadband electromagnetic coverage of binary neutron star merger (GW170817) has also placed constraints on the properties (mass and composition) of the merger ejecta. But uncertainties and ambiguities in modeling the optical and infrared emission make it difficult to definitively measure the distribution of heavy isotopes in these mergers. In contrast, gamma rays emitted in the decay of these neutron-rich ejecta may provide a more direct measurement of the yields. We calculate the gamma production in remnants of neutron star mergers, considering two epochs: a kilonova epoch, lasting about two weeks, and a much later epoch of tens and hundreds of thousands of years after the merger. For the kilonova epoch, when the expanding ejecta is still only partially transparent to gamma radiation, we use 3D radiative transport simulations to produce the spectra. Futhermore, we show that the gamma-ray spectra associated with beta- and alpha-decay provide a fingerprint of the ejecta properties and, for a sufficiently nearby remnant, may be detectable, even for old remnants. We compare our gamma spectra with the potential detection limits of next generation detectors, including the Lunar Occultation Explorer (LOX), the All-sky Medium Energy Gamma-ray Observatory (AMEGO), and the Compton Spectrometer and Imager (COSI). We show that fission models can be discriminated via the presence of short-lived fission fragments in the remnant spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

Diagnostics of Magnetohydrodynamic Modes in the Interstellar Medium through Synchrotron Polarization Statistics

One of the biggest challenges in understanding magnetohydrodynamic (MHD) turbulence is identifying the plasma mode components from observational data. Previous studies on synchrotron polarization from the interstellar medium (ISM) suggest that the dominant MHD modes can be identified via statistics of Stokes parameters, which would be crucial for studying various ISM processes such as the scattering and acceleration of cosmic rays, star formation, and dynamo. In this paper, we present a numerical study of the synchrotron polarization analysis (SPA) method through systematic investigation of the statistical properties of the Stokes parameters. We derive the theoretical basis for our method from the fundamental statistics of MHD turbulence, recognizing that the projection of the MHD modes allows us to identify the modes dominating the energy fraction from synchrotron observations. Based on the discovery, we revise the SPA method using synthetic synchrotron polarization observations obtained from 3D ideal MHD simulations with a wide range of plasma parameters and driving mechanisms, and present a modified recipe for mode identification. We propose a classification criterion based on a new SPA+ fitting procedure, which allows us to distinguish between Alfvén mode and compressible/slow mode dominated turbulence. We further propose a new method to identify fast modes by analyzing the asymmetry of the SPA+ signature and establish a new asymmetry parameter to detect the presence of fast mode turbulence. Additionally, we confirm through numerical tests that the identification of the compressible and fast modes is not affected by Faraday rotation in both the emitting plasma and the foreground.

97 MATHEMATICS AND COMPUTING↗

MAGNETIC HELICITY REVERSALS IN A CYCLIC CONVECTIVE DYNAMO

We investigate the role of magnetic helicity in promoting cyclic magnetic activity in a global, 3D, magnetohydrodynamic (MHD) simulation of a convective dynamo. This simulation is characterized by coherent bands of toroidal field that exist within the convection zone, with opposite polarities in the northern hemisphere (NH) and southern hemisphere (SH). Throughout most of the cycle, the magnetic helicity in these bands is negative in the NH and positive in the SH. However, during the declining phase of each cycle, this hemispheric rule reverses. We attribute this to a global restructuring of the magnetic topology that is induced by the interaction of the bands across the equator. This band interaction appears to be ultimately responsible for, or at least associated with, the decay and subsequent reversal of both the toroidal bands and the polar fields. We briefly discuss the implications of these results within the context of solar observations, which also show some potential evidence for toroidal band interactions and helicity reversals.

79 ASTRONOMY AND ASTROPHYSICS↗

Multipoint Turbulence Analysis with HelioSwarm

Exploration of plasma dynamics in space, including turbulence, is entering a new era of multisatellite constellation measurements that will determine fundamental properties with unprecedented precision. Familiar but imprecise approximations will need to be abandoned and replaced with more-advanced approaches. We present a preparatory study of the evaluation of second- and third-order statistics, using simultaneous measurements at many points. Here, for specificity, the orbital configuration of the NASA Swarm mission is employed in conjunction with 3D magnetohydrodynamics numerical simulations of turbulence. The HelioSwarm nine-spacecraft constellation flies virtually through the turbulence to compare results with the exact numerical statistics. We demonstrate novel increment-based techniques for the computation of (1) the multidimensional spectra and (2) the turbulent energy flux. This latter increment-space estimate of the cascade rate, based on the third-order Yaglom–Politano–Pouquet theory, uses numerous increment-space tetrahedra. Our investigation reveals that HelioSwarm will provide crucial information on the nature of astrophysical turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Geometry of the Near-core Magnetic Field in Massive Stars

It is well-known that the cores of massive stars sustain a stellar dynamo with a complex magnetic field configuration. However, the same cannot be said for the field's strength and geometry at the convective–radiative boundary, which are crucial when performing asteroseismic inference. In this Letter, we present 3D magnetohydrodynamic (MHD) simulations of a 7 M ⊙ mid-main-sequence star, with particular attention given to the convective–radiative boundary in the near-core region. Our simulations reveal that the toroidal magnetic field is significantly stronger than the poloidal field in this region, contrary to recent assumptions. Moreover, the rotational shear layer, also important for asteroseismic inference, is specifically confined within the extent of the Brunt–Väisälä frequency peak. These results, which are based on the inferred properties of HD 43317, have widespread implications for asteroseismic studies of rotation, mixing, and magnetism in stars. While we expect our results to be broadly applicable across stars with similar Brunt–Väisälä frequency profiles and stellar masses, we also expect the MHD parameters (e.g., Re m ) and the initial stellar rotation rate to impact the geometry of the field and differential rotation at the convective–radiative interface.

79 ASTRONOMY AND ASTROPHYSICS↗

On the Statistical Uncertainty of Monte Carlo-Calculated Scattering Sensitivities

Sensitivity coefficients calculated with Monte Carlo codes are widely used for nuclear data uncertainty quantification in the modeling and simulation of complex 3D reactor systems. This study systematically compares sensitivity coefficients and associated statistical uncertainties for the multiplication factor and fuel temperature reactivity across multiple Monte Carlo codes (SCALE/KENO, SCALE/Shift, MCNP, and Serpent) using simple models representing light-water reactors and advanced reactor concepts. For multiplication factor sensitivities, statistical uncertainties are generally acceptable, although scattering sensitivities show significantly larger statistical uncertainties than, for example, fission and capture reactions. Fuel temperature reactivity sensitivities show significantly larger statistical uncertainties across all reactions. Elastic scattering sensitivities are the most problematic: all Monte Carlo codes fail to resolve energy-dependent coefficients, and they produce dramatically different energy-collapsed values. Critically, the use of these sensitivity coefficients in nuclear data uncertainty propagation leads to reduced statistical uncertainties in individual uncertainty contributions. This can lead to the masking of unusable sensitivity coefficients and producing misleading uncertainty results. The findings of this study show that new or enhanced methods are needed to improve Monte Carlo elastic scattering sensitivity calculations. Additionally, this study shows the relevance of verifying sensitivity coefficients through direct perturbation calculations for individual nuclide reactions, instead of only for total cross sections as commonly done.

Bostelmann, Rike [ORNL] (ORCID:0000000165968088)↗