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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Reduction of asymmetric wall force in ITER disruptions with fast current quench

One of the problems caused by disruptions in tokamaks is the asymmetric electromechanical force produced in conducting structures surrounding the plasma. The asymmetric wall force in ITER asymmetric vertical displacement event (AVDE) disruptions is calculated in nonlinear 3D MHD simulations. It is found that the wall force can vary by almost an order of magnitude, depending on the ratio of the current quench time to the resistive wall magnetic penetration time. In ITER, this ratio is relatively low, resulting in a low asymmetric wall force. In JET, this ratio is relatively high, resulting in a high asymmetric wall force. Previous extrapolations based on JET measurements have greatly overestimated the ITER wall force. Here, it is shown that there are two limiting regimes of AVDEs, and it is explained why the asymmetric wall force is different in the two limits.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A new explanation of the sawtooth phenomena in tokamaks

The ubiquitous sawtooth phenomena in tokamaks are so named because the central temperature rises slowly and falls rapidly, similar to the blades of a saw. First discovered in 1974, it has so far eluded a theoretical explanation that is widely accepted and consistent with experimental observations. We propose here a new theory for the sawtooth phenomena in auxiliary heated tokamaks, which is motivated by our recent understanding of “magnetic flux pumping.” In this theory, the role of the (m, n) = (1, 1) mode is to generate a dynamo voltage, which keeps the central safety factor, q 0 , just above 1.0 with low central magnetic shear. When central heating is present, the temperature on axis will increase until at some point, and the configuration abruptly becomes unstable to ideal MHD interchange modes with equal poloidal and toroidal mode numbers, m = n > 1. It is these higher order modes and the localized magnetic stochasticity they produce that cause the sudden crash of the temperature profile, not magnetic reconnection. In conclusion, long time 3D MHD simulations demonstrate these phenomena, which appear to be consistent with many experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Plasma expansion and relativistic filamentation in intense laser-irradiated cone targets

We report that compound parabolic concentrator (CPC) cone targets have been shown to produce increased MeV photons on the NIF-ARC by 10× over flat targets. Multiple x-ray frames can potentially be generated by firing the NIF-ARC's beamlets into distinct cone targets at few nanosecond relative delays. This requires that the cone targets with delayed beams are not degraded by their proximity to previous targets. One concern is that the spatial wings of a beam fired into one target can fall on neighboring targets, producing a preformed plasma that may interfere with laser light reaching the tip of the cone. In this work, 3D hydra simulations of realistic targets and beam parameters show that hundreds of micrometer scale length preplasmas are produced in cones within 1 mm of the laser spot. 2D particle-in-cell simulations of the intense main pulse in this preplasma indicate a density threshold for the onset of relativistic filamentation in our conditions. Applying our modeling approach to a NIF-ARC shot with an intentional 15 J prepulse yields good agreement with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Effect of void positioning on the detonation sensitivity of a heterogeneous energetic material

We show although it is well-established that voids profoundly influence the initiation and reaction behaviors of heterogeneous energetic materials such as polymer-bonded explosives (PBX) and propellants, there has been little study of how void location in different constituents in the microstructures of such materials affect the macroscale behavior. Here, we use three-dimensional (3D) mesoscale simulations to study how void placement within the reactive grains versus the polymer binder influences the shock-to-detonation transition (SDT) in a polymer-bonded explosive. The material studied here has a microstructure comprised of 75% PETN (pentaerythritol tetranitrate) grains and 25% HTPB (hydroxyl-terminated polybutadiene) polymer binder by volume. Porosities up to 10% in the form of spherical voids distributed in both the grains and polymer are considered. An Arrhenius reactive burn relation is used to model the chemical kinetics of the PETN grains under shock loading, thereby resolving the heterogeneous detonation behavior of the PBX. The influence of void location on the shock initiation sensitivity of the material is quantitatively ranked by comparing the predicted run distance to detonation (RDD) for each sample. The analysis includes inherent quantification of uncertainties arising from the stochastic variations in the microstructure morphologies and void distributions by using statistically equivalent microstructure sample sets (SEMSS), leading to probabilistic formulations for the RDD as a function of shock pressure. The calculations reveal that the location of voids in the composite microstructure significantly affects the RDD. Specifically, voids exclusively within the grains cause the PBX to be more sensitive (having shorter RDD) than voids in the polymer binder. Unique probabilistic relationships are derived to map the probability of observing RDD for each void location material case, allowing for prediction of initiation behavior anywhere in the shock pressure – RDD space. These findings agree with trends reported in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Multi-fidelity modeling to predict the rheological properties of a suspension of fibers using neural networks and Gaussian processes

Unveiling the rheological properties of fiber suspensions is of paramount interest to many industrial applications. There are multiple factors, such as fiber aspect ratio and volume fraction, that play a significant role in altering the rheological behavior of suspensions. Three-dimensional (3D) numerical simulations of coupled differential equations of the suspension of fibers are computationally expensive and time-consuming. Machine learning algorithms can be trained on the available data and make predictions for the cases where no numerical data are available. However, some widely used machine learning surrogates, such as neural networks, require a relatively large training dataset to produce accurate predictions. Multi-fidelity models, which combine high-fidelity data from numerical simulations and less expensive lower fidelity data from resources such as simplified constitutive equations, can pave the way for more accurate predictions. Here, we focus on neural networks and the Gaussian processes with two levels of fidelity, i.e., high and low fidelity networks, to predict the steady-state rheological properties, and compare them to the single-fidelity network. High-fidelity data are obtained from direct numerical simulations based on an immersed boundary method to couple the fluid and solid motion. The low-fidelity data are produced by using constitutive equations. Multiple neural networks and the Gaussian process structures are used for the hyperparameter tuning purpose. Results indicate that with the best choice of hyperparameters, both the multi-fidelity Gaussian processes and neural networks are capable of making predictions with a high level of accuracy with neural networks demonstrating marginally better performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Diagnosing magnetic fields in cylindrical implosions with oblique proton radiography

Two experiments at the OMEGA Laser System used oblique proton radiography to measure magnetic fields in cylindrical implosions with and without an applied axial magnetic field. Although the goal of both experiments was to measure the magnitude of the compressed axial magnetic field in the core of the implosion, this field was obfuscated by two features in the coronal plasma produced by the compression beams: an azimuthal self-generated magnetic field and small length scale, high-amplitude structures attributed to collisionless effects. In order to understand these features, synthetic radiographs are generated using fields produced by 3D HYDRA simulations. These synthetic radiographs reproduce the features of the experimental radiographs with the exception of the small-scale structures. A direct inversion algorithm is successfully applied to a synthetic radiograph but is only partially able to invert the experimental radiographs in part because some protons are blocked by the field coils. The origins of the radiograph features and their dependence on various experimental parameters are explored. Furthermore, the results of this analysis should inform future measurements of compressed axial magnetic fields in cylindrical implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Exploration of cross-beam energy transfer mitigation constraints for designing an ignition-scale direct-drive inertial confinement fusion driver

The compression of direct-drive inertial confinement fusion (ICF) targets is strongly impacted by cross-beam energy transfer (CBET), a laser-plasma instability that limits ablation pressure by redirecting laser energy outward and that is projected to be mitigated by laser bandwidth. Here, we explore various CBET mitigation constraints to guide the design of future ICF facilities. First, we find that the flat, Gaussian, and Lorentzian spectral shapes have similar CBET mitigation properties, and a flat shape with nine spectral lines is a good surrogate for what can be obtained with other spectral shapes. Then, we conduct a comprehensive study across energy scales and ignition designs. 3D hydrodynamic simulations are used to derive an analytical model for the expected CBET mitigation as a function of laser and plasma parameters. From this model, we study the bandwidth requirements of conventional and shock ignition designs across four different energy scales and find that they require between 0.5 and 3±0.2% relative bandwidth. Further, best mitigation is achieved when the beam radius over critical radius R b /R c is kept low during the drive while the plasma temperature is kept high. In a steady state, we find that the bandwidth required to mitigate 85% of CBET scales as (R b /R c ) 2.15 L n –0.58 I 0.7 , where L n is the density scale length, and I the laser intensity. Finally, we find that the chamber beam port layout does not influence CBET mitigation. In the case of a driver using many monochromatic beamlets, we find that ~10 beamlets per port is required, with diminishing returns above ~20.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

On discriminating tokamak sawtooth crash models via localized density and temperature measurements

The core electron temperature drops rapidly during the sawtooth crash in tokamak plasmas, which causes heat loss and may lead to fast particle losses or even a disruption. Several models have been proposed for the periodic crash, including the Kadomtsev model with magnetic reconnection and the quasi-interchange model with the growth of higher-mode-number pressure-driven instabilities. 3D MHD simulations were performed for these two models with a goal to develop intuition and to predict qualitatively how different types of sawtooth will appear in various diagnostics. The structures of electron density ne and electron temperature Te show a dominant (1, 1) mode for the Kadomtsev case and a dominant (4, 4) mode for the quasi-interchange case. The oscillations of ne and Te have a positive correlation near the inversion layer for both cases, while their frequencies and amplitudes are different depending on the dominant modes. Particularly, for the Kadomtsev case, we find a relation between the amount of flux reconnected during a sawtooth event and ne or Te oscillations. Therefore, we connect recently developed measurement capabilities for ne and Te to the internal sawtooth behavior. We propose that this method of analysis can help in identifying the type of sawtooth in future experiments augmented by simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Explosively driven Richtmyer–Meshkov instability jet suppression and enhancement via coupling machine learning and additive manufacturing

The ability to control the behavior of fluid instabilities at material interfaces, such as the shock-driven Richtmyer–Meshkov instability, is a grand technological challenge with a broad number of applications ranging from inertial confinement fusion experiments to explosively driven shaped charges. In this work, we use a linear-geometry shaped charge as a means of studying methods for controlling material jetting that results from the Richtmyer–Meshkov instability. A shaped charge produces a high-velocity jet by focusing the energy from the detonation of high explosives. The interaction of the resulting detonation wave with a hollowed cavity lined with a thin metal layer produces the unstable jetting effect. By modifying the characteristics of the detonation wave prior to striking the lined cavity, the kinetic energy of the jet can be enhanced or reduced. Modifying the geometry of the liner material can also be used to alter jetting properties. We apply optimization methods to investigate several design parameterizations for both enhancing or suppressing the shaped-charge jet. This is accomplished using 2D and 3D hydrodynamic simulations to investigate the design space that we consider. We also apply new additive manufacturing methods for producing the shaped-charge assemblies, which allow for the experimental testing of complicated design geometries obtained through computational optimization. We present a direct comparison of our optimized designs with experimental results carried out at the High Explosives Application Facility at Lawrence Livermore National Laboratory.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Toward an integrated multi-gigahertz ionizing particle diagnostic

Needs arising at both current and future accelerator facilities call for the development of radiation-hardened position-sensing diagnostics that can operate with multi-GHz repetition rates. Such instruments are likely to also have applications in the diagnosis of rapid plasma behavior. Here, building on the recent work of our Advanced Accelerator Diagnostics Collaboration, we are exploring the development of integrated multi-GHz ionizing particle detection systems based on chemical-vapor deposition diamond sensors, with the initial goal of producing a quadrant detector that can determine the intensity and centroid position of a particle beam at a repetition rate between 5 and 10 GHz. Results from our initial high-speed characterization work are presented, including those from a single-channel sensor with a GHz response. Approaches to achieving multi-GHz (5–10 GHz) rate capability, including the design of a dedicated Application Specific Integrated Circuit and the use of 3D RF-solver computer aided design software, are presented and discussed in more detail. 3D RF simulations suggest clean pulses of duration less than 250 ps (FWHM less than 125 ps) can be achieved with the approaches developed by this work.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Fast ion stabilization of tilt in large radius FRCs

The field reversed configuration (FRC) has been a curious case in plasma physics research in that early MHD analysis suggested FRCs should be grossly unstable, while experimental results contradicted that prediction. Later, this theory was able to resolve this contradiction by understanding that finite Larmor radius effects largely negated the MHD predictions. Similarly, previous theoretical studies of beam driven FRCs predicted that such system would be unstable to beam driven modes while, again, experimental results indicated the contradiction. In this paper, we reconcile the theoretical understanding of beam driven modes with experimental observations of stability in these systems. By self-consistently capturing fast ion generation from neutral beam injection and its impact on the plasma equilibrium, we show that low amplitude perturbations in the magnetic field, driven by betatron particles, modify the precession frequencies of the betatron particles such that the drive for compressional Alfvén waves in the thermal plasma is reduced. Finally, we are able to demonstrate, for the first time, stable beam driven FRC evolution at high S*/E in 3D kinetic simulations.

Betatrons↗
DOE OSTI2025

Rotation of electrothermal-instability-driven overheating structure due to helically oriented surface magnetic field on a high-current-density aluminum rod

Experiments on the 1-MA, 100-ns-rise-time Mykonos Facility demonstrate rotation of electrothermal instability (ETI)-driven overheating structure on 1.00-mm-diameter, 10-nm-surface roughness, 99.999%-pure aluminum rods, which are pulsed with helically polarized surface magnetic field. Rods are machined to include pairs of 10-micron-scale quasi-hemispherical voids or “engineered defects (ED)” which provide the dominant current density perturbation from which ETI grows most rapidly. Experiments include an axial magnetic field component through the addition of a helically wound return-current electrode or “helical return can (HRC).” For a given HRC design, azimuthal field (B ɵ ) and axial field (B z ) components rise at a prescribed and fixed ratio, driving an increasing magnetic field of constant polarization at the rod's surface; most experiments generated surface magnetic field at a 15-degree field polarization angle (from horizontal) defined as ɸ B = arctan(B z /B ɵ ). ETI-driven emission patterns from individual ED are observed to rotate along ɸ B , while emission patterns from dielectric-coated ED pairs are shown to elongate and preferentially merge along ɸ B , in qualitative agreement with 3D-magnetohydrodynamic simulations. These data strongly support that for a randomized distribution of current density perturbations on a high-current density conductor, nearby perturbations will favorably merge about ɸ B , with the degree of merging increasing with current. Such observations offer fundamental new understanding of the seeding mechanisms of the helical magneto-Rayleigh Taylor (MRT) instabilities observed from axially magnetized magnetically driven imploding liners.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

A NEW MODEL FOR MIXING BY DOUBLE-DIFFUSIVE CONVECTION (SEMI-CONVECTION). II. THE TRANSPORT OF HEAT AND COMPOSITION THROUGH LAYERS

Regions of stellar and planetary interiors that are unstable according to the Schwarzschild criterion, but stable according to the Ledoux criterion, are subject to a form of oscillatory double-diffusive (ODD) convection often called ''semi-convection''. In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and composition by ODD convection, and ultimately propose a new 1D prescription that can be used in stellar and planetary structure and evolution models. The first paper in this series demonstrated that under certain conditions ODD convection spontaneously transitions from an initial homogeneous state of weak wave-breaking turbulence into a staircase of fully convective layers, which results in a substantial increase in the transport of heat and composition. Here, we present simulations of ODD convection in this layered regime, we describe the dynamical behavior of the layers, and we derive empirical scaling laws for the transport through layered convection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Cosmic ray transport in heliospheric magnetic structures. I. Modeling background solar wind using the CRONOS magnetohydrodynamic code

The transport of energetic particles such as cosmic rays is governed by the properties of the plasma being traversed. While these properties are rather poorly known for galactic and interstellar plasmas due to the lack of in situ measurements, the heliospheric plasma environment has been probed by spacecraft for decades and provides a unique opportunity for testing transport theories. Of particular interest for the three-dimensional (3D) heliospheric transport of energetic particles are structures such as corotating interaction regions, which, due to strongly enhanced magnetic field strengths, turbulence, and associated shocks, can act as diffusion barriers on the one hand, but also as accelerators of low energy CRs on the other hand as well. In a two-fold series of papers, we investigate these effects by modeling inner-heliospheric solar wind conditions with a numerical magnetohydrodynamic (MHD) setup (this paper), which will serve as an input to a transport code employing a stochastic differential equation approach (second paper). In this first paper, we present results from 3D MHD simulations with our code CRONOS: for validation purposes we use analytic boundary conditions and compare with similar work by Pizzo. For a more realistic modeling of solar wind conditions, boundary conditions derived from synoptic magnetograms via the Wang-Sheeley-Arge (WSA) model are utilized, where the potential field modeling is performed with a finite-difference approach in contrast to the traditional spherical harmonics expansion often utilized in the WSA model. Our results are validated by comparing with multi-spacecraft data for ecliptical (STEREO-A/B) and out-of-ecliptic (Ulysses) regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

FREQUENCY SHIFTS OF RESONANT MODES OF THE SUN DUE TO NEAR-SURFACE CONVECTIVE SCATTERING

Measurements of oscillation frequencies of the Sun and stars can provide important independent constraints on their internal structure and dynamics. Seismic models of these oscillations are used to connect structure and rotation of the star to its resonant frequencies, which are then compared with observations, the goal being that of minimizing the difference between the two. Even in the case of the Sun, for which structure models are highly tuned, observed frequencies show systematic deviations from modeled frequencies, a phenomenon referred to as the “surface term.” The dominant source of this systematic effect is thought to be vigorous near-surface convection, which is not well accounted for in both stellar modeling and mode-oscillation physics. Here we bring to bear the method of homogenization, applicable in the asymptotic limit of large wavelengths (in comparison to the correlation scale of convection), to characterize the effect of small-scale surface convection on resonant-mode frequencies in the Sun. We show that the full oscillation equations, in the presence of temporally stationary three-dimensional (3D) flows, can be reduced to an effective “quiet-Sun” wave equation with altered sound speed, Brünt–Väisäla frequency, and Lamb frequency. We derive the modified equation and relations for the appropriate averaging of 3D flows and thermal quantities to obtain the properties of this effective medium. Using flows obtained from 3D numerical simulations of near-surface convection, we quantify their effect on solar oscillation frequencies and find that they are shifted systematically and substantially. We argue therefore that consistent interpretations of resonant frequencies must include modifications to the wave equation that effectively capture the impact of vigorous hydrodynamic convection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

THERMAL DIAGNOSTICS WITH THE ATMOSPHERIC IMAGING ASSEMBLY ON BOARD THE SOLAR DYNAMICS OBSERVATORY: A VALIDATED METHOD FOR DIFFERENTIAL EMISSION MEASURE INVERSIONS

We present a new method for performing differential emission measure (DEM) inversions on narrow-band EUV images from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory. The method yields positive definite DEM solutions by solving a linear program. This method has been validated against a diverse set of thermal models of varying complexity and realism. These include (1) idealized Gaussian DEM distributions, (2) 3D models of NOAA Active Region 11158 comprising quasi-steady loop atmospheres in a nonlinear force-free field, and (3) thermodynamic models from a fully compressible, 3D MHD simulation of active region (AR) corona formation following magnetic flux emergence. We then present results from the application of the method to AIA observations of Active Region 11158, comparing the region's thermal structure on two successive solar rotations. Additionally, we show how the DEM inversion method can be adapted to simultaneously invert AIA and Hinode X-ray Telescope data, and how supplementing AIA data with the latter improves the inversion result. The speed of the method allows for routine production of DEM maps, thus facilitating science studies that require tracking of the thermal structure of the solar corona in time and space.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

Growth and evolution of thermal instabilities in idealized galaxy cluster cores

Heat input roughly balances radiative cooling in the gaseous cores of galaxy clusters even when the central cooling time is short, implying that cooling triggers a feedback loop that maintains thermal balance. Furthermore, cores with short cooling times tend to have multiphase structure, suggesting that the intracluster medium (ICM) becomes locally thermally unstable for cooling times ≲1 Gyr. Both observations and theoretical models have linked the condensation of cold gas with heating by an active galactic nucleus (AGN) through a cycle in which cooling gas fuels the AGN and drives energetic outbursts that reheat the ICM and maintain a state of approximate thermal balance. In this work, we use 2D and 3D hydrodynamic simulations to study the onset of condensation in idealized galaxy-cluster cores. In particular, we look at how the condensation process depends on the ratio of cooling time to freefall time and on the geometry of the gravitational potential. We conclude that the ICM can always evolve to a state in which condensation occurs if given enough time, but that an initial timescale ratio t{sub cool}/t{sub ff}≲10 is needed for thermal instability to grow quickly enough to affect realistic cluster cores within a timescale that is relevant for cosmological structure formation. We find that instability leads to convection and that perturbations continue to grow while the gas convects. Condensation occurs when the timescale ratio in the low-entropy tail of the perturbation distribution drops below t{sub cool}/t{sub ff}≲3, even if the volume-averaged timescale ratio is substantially greater. In our simulations, the geometry of the gravitational potential does not have a strong effect on thermal stability. Finally, we find that if condensation is powering feedback, a conversion efficiency of around 10{sup −3} for converting the condensed mass into thermal energy is sufficient to maintain thermal balance in the ICM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

SMALL-SCALE HEATING EVENTS IN THE SOLAR ATMOSPHERE. I. IDENTIFICATION, SELECTION, AND IMPLICATIONS FOR CORONAL HEATING

We present a comprehensive method to analyze small-scale heating events in detail in a 3D magnetohydrodynamics simulation for quiet-Sun conditions. The method determines the number, volume, and some general geometric properties of the small-scale heating events at different instants in a simulation with a volume of 16 × 8 × 16 Mm{sup 3}, spanning from the top of the convection zone to the corona. We found that there are about 10{sup 4} small-scale heating events at any instant above the simulated area of 128 Mm{sup 2}. They occur mainly at heights between 1.5 and 3.0 Mm. We determine the average value of their projected vertical extent, which ranges from 375 to 519 km over time, and we show that height, volume, and energy distribution of the events at any instant resemble power laws. Finally, we demonstrate that larger heating structures are a combination of much smaller heating events and that small-scale heating events dissipate enough energy to maintain the coronal energetic balance at any instant.

79 ASTRONOMY AND ASTROPHYSICS↗