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 109 records · Page 6

HYSTERESIS BETWEEN DISTINCT MODES OF TURBULENT DYNAMOS

Nonlinear mean-field models of the solar dynamo show long-term variability, which may be relevant to different states of activity inferred from long-term radiocarbon data. This paper is aimed at probing the dynamo hysteresis predicted by the recent mean-field models of Kitchatinov and Olemskoy with direct numerical simulations. We perform three-dimensional (3D) simulations of large-scale dynamos in a shearing box with helically forced turbulence. As an initial condition, we either take a weak random magnetic field or we start from a snapshot of an earlier simulation. Two quasi-stable states are found to coexist in a certain range of parameters close to the onset of the large-scale dynamo. The simulations converge to one of these states depending on the initial conditions. When either the fractional helicity or the magnetic Prandtl number is increased between successive runs above the critical value for onset of the dynamo, the field strength jumps to a finite value. However, when the fractional helicity or the magnetic Prandtl number is then decreased again, the field strength stays at a similar value (strong field branch) even below the original onset. We also observe intermittent decaying phases away from the strong field branch close to the point where large-scale dynamo action is just possible. The dynamo hysteresis seen previously in mean-field models is thus reproduced by 3D simulations. Its possible relation to distinct modes of solar activity such as grand minima is discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Validation of 3D MHD simulations of Ne/D 2 mixed shattered pellet injection in JET

Nonlinear 3D magnetohydrodynamic (MHD) modeling of shattered pellet injection (SPI) in JET tokamak plasmas is performed with the JOREK code. The study focuses on the validation of simulation result with respect to experimental observations, addressing in particular figures of merit for the efficiency of the SPI, as a technique to mitigate thermal loads by injecting radiative impurities like neon during the thermal quench (TQ) phase of disruptions. A set of JET pulses with neon/deuterium atomic mixture ratio in the shattered pellet varying from 10% to 100% provides the experimental data to be compared with. Simulations using different models for the ablation of solid fragments and the radiation of neon impurities are considered. Synthetic diagnostics are employed for a direct quantitative comparison with key JET diagnostic systems providing in particular radiation, electron temperature and density, and magnetic measurements. 3D radiation structures from MHD simulations are analyzed to estimate the toroidal asymmetry of radiation, which is difficult to measure in JET due to the presence of only two toroidally displaced multi-channel bolometry systems. JOREK simulations show that, before the TQ, the bulk of the radiation is concentrated in a toroidal region enclosing the radiating fragments that is missed by both multi-channel bolometers, suggesting a possible underestimation of the total radiation by such diagnostic systems in the pre-TQ phase. On the other hand, JOREK predicts that during the TQ the radiation becomes more toroidally symmetric, with positive implications for the heat load to plasma facing components.

JOREK↗

Black hole to breakout: 3D GRMHD simulations of collapsar jets reveal a wide range of transients

We present a suite of the first 3D GRMHD collapsar simulations, which extend from the self-consistent jet launching by an accreting Kerr black hole (BH) to the breakout from the star. We identify three types of outflows, depending on the angular momentum, l, of the collapsing material and the magnetic field, B, on the BH horizon: (i) subrelativistic outflow (low l and high B), (ii) stationary accretion shock instability (SASI; high l and low B), (iii) relativistic jets (high l and high B). In the absence of jets, free-fall of the stellar envelope provides a good estimate for the BH accretion rate. Jets can substantially suppress the accretion rate, and their duration can be limited by the magnetization profile in the star. We find that progenitors with large (steep) inner density power-law indices (≳ 2), face extreme challenges as gamma-ray burst (GRB) progenitors due to excessive luminosity, global time evolution in the light curve throughout the burst and short breakout times, inconsistent with observations. Our results suggest that the wide variety of observed explosion appearances (supernova/supernova + GRB/low-luminosity GRBs) and the characteristics of the emitting relativistic outflows (luminosity and duration) can be naturally explained by the differences in the progenitor structure. Our simulations reveal several important jet features: (i) strong magnetic dissipation inside the star, resulting in weakly magnetized jets by breakout that may have significant photospheric emission and (ii) spontaneous emergence of tilted accretion disc-jet flows, even in the absence of any tilt in the progenitor.

79 ASTRONOMY AND ASTROPHYSICS↗

3D MHD SIMULATION OF FLARE SUPRA-ARCADE DOWNFLOWS IN A TURBULENT CURRENT SHEET MEDIUM

Supra-arcade downflows (SADs) are sunward, generally dark, plasma density depletions originated above posteruption flare arcades. In this paper, using 3D MHD simulations we investigate whether the SAD cavities can be produced by a direct combination of the tearing mode and Kelvin–Helmholtz instabilities leading to a turbulent current sheet (CS) medium or if the current sheet is merely the background where SADs are produced, triggered by an impulsive deposition of energy. We find that to give an account of the observational dark lane structures an addition of local energy, provided by a reconnection event, is required. We suggest that there may be a closed relation between characteristic SAD sizes and CS widths that must be satisfied to obtain an observable SAD.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermal Performance of Spandrel Assemblies in Glazed Wall Systems: Laboratory Test Design – Challenges and Test Results

Accurate thermal performance calculation procedures for opaque spandrel areas in curtain wall and window wall systems are essential for rating systems when comparing spandrel systems. However, there is a lack of consensus in thermal modeling needed for accurately characterizing heat transfer through spandrel assemblies due to the complex arrangement of materials and structural components. Several studies indicate that conventional 2D thermal simulations may overestimate R-values by 30% compared to physical testing and 3D simulations. Detailed simulations and well-curated laboratory test data are necessary to build confidence in simulation models, which will later be used to develop correlations to improve widely used conventional 2D thermal simulations. This study aims to experimentally test heat transfer through various spandrel assemblies to validate 3D simulation models. Also, the challenges of conducting a thorough testing design along with the solutions would be documented. The team developed a design for testing spandrel assemblies, making appropriate modifications to the existing heat, air, and moisture (HAM) chamber to accommodate the testing needs. Two moveable baffles were designed and fabricated to guide airflow direction parallel to the test article surface. The data acquisition capabilities in the chamber were upgraded to add more than two hundred sensors to the climate and indoor side of the chamber. The goal is to provide a quality dataset for validating complex 3D modeling simulations, which will be used to develop improved thermal simulation techniques that more accurately represent the thermal behavior of spandrel assemblies and their integration within the building envelope. This paper will summarize the results for the boundary conditions of the testing and the temperature variation across different locations of the spandrel assemblies.

Kunwar, Niraj [ORNL] (ORCID:0000000263457652)↗

The thermal regime of the Campi Flegrei magmatic system reconstructed through 3D numerical simulations

In this paper, we illustrate a quantitative conductive/convective thermal model incorporating a wide range of geophysical, petrological, geological, geochemical and isotopical observations that constrain the thermal evolution and present state of the Campi Flegrei caldera (CFc) magmatic system. The proposed model has been computed on the basis of the current knowledge of: (1) the volcanic and magmatic history of the volcano over the last 44 ka, (2) its underlying crustal structure, and (3) the physical properties of the erupted magmas. 3D numerical simulations of heat conduction and convection within heterogeneous rock/magma materials with evolving heat sources and boundary conditions that simulate magma rise from a deep (≥ 8 km depth) to shallow (2–6 km) reservoirs, magma chamber formation, magma extrusion, caldera collapse, and intra-caldera hydrothermal convection, have been carried out. The evolution of the CFc magmatic system through time has been simulated through different steps related to its changes in terms of depth, location and size of magma reservoirs and their replenishment. The thermal modeling results show that both heat conduction and convection have played an important role in the CFc thermal evolution, although with different timing. Finally, the simulated present heat distribution is in agreement with the measured geothermal profiles (Agip, 1987), reproduces the thermal gradient peaks at the CFc margins in correspondence to the anomalies in surface gradients (Corrado et al., 1998), and suggests temperatures of 700 °C at depth of 4 km in the central portion of the caldera, in agreement with the estimated temperature for the brittle-ductile transition (Hill, 1992).

3D conductive/convective model↗

3D xRAGE simulation of inertial confinement fusion implosion with imposed mode 2 laser drive asymmetry

Low-mode asymmetries represent an important obstacle to achieving high-gain inertial confinement fusion implosions. As a step in learning how to control such effects, an OMEGA experiment with imposed mode 2 laser drive asymmetries was done to study the expected signatures of this type of asymmetry [M. Gatu Johnson et al., PRE 2018]. In the present work, a 3D xRAGE simulation including the stalk mount has been brought to bear on the data from that experiment. Furthermore, comprehensive comparisons between simulated and measured observables are made. Good agreement between simulated and measured x-ray image-inferred shell trajectories, bang times and neutron emission widths are seen, showing that the hydrodynamics are well captured in the simulation. Asymmetries seen in simulated and measured time-resolved and time-integrated x-ray images and areal densities also compare well, showing impact of both stalk and mode 2. On the other hand, important differences in measured and simulated neutron emission histories, yield, and ion temperature (T ion ) asymmetries are seen, suggesting that the simulation is overestimating shock yield. Overall, the results clearly demonstrate the importance of considering all asymmetry sources when interpreting measured signatures of asymmetry.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

3D numerical simulation of the long range propagation of acoustical shock waves through a heterogeneous and moving medium

Many situations involve the propagation of acoustical shock waves through flows. Natural sources such as lightning, volcano explosions, or meteoroid atmospheric entries, emit loud, low frequency, and impulsive sound that is influenced by atmospheric wind and turbulence. The sonic boom produced by a supersonic aircraft and explosion noises are examples of intense anthropogenic sources in the atmosphere. The Buzz-Saw-Noise produced by turbo-engine fan blades rotating at supersonic speed also propagates in a fast flow within the engine nacelle. Simulating these situations is challenging, given the 3D nature of the problem, the long range propagation distances relative to the central wavelength, the strongly nonlinear behavior of shocks associated to a wide-band spectrum, and finally the key role of the flow motion. With this in view, the so-called FLHOWARD (acronym for FLow and Heterogeneous One-Way Approximation for Resolution of Diffraction) method is presented with three-dimensional applications. A scalar nonlinear wave equation is established in the framework of atmospheric applications, assuming weak heterogeneities and a slow wind. It takes into account diffraction, absorption and relaxation properties of the atmosphere, quadratic nonlinearities including weak shock waves, heterogeneities of the medium in sound speed and density, and presence of a flow (assuming a mean stratified wind and 3D turbulent ? flow fluctuations of smaller amplitude). This equation is solved in the framework of the one-way method. A split-step technique allows the splitting of the non-linear wave equation into simpler equations, each corresponding to a physical effect. Each sub-equation is solved using an analytical method if possible, and finite-differences otherwise. Nonlinear effects are solved in the time domain, and others in the frequency domain. Homogeneous diffraction is handled by means of the angular spectrum method. Ground is assumed perfectly flat and rigid. Due to the 3D aspect, the code was massively parallelized using the single program, multiple data paradigm with the Message Passing Interfaces (MPI) for distributed memory architectures. This allows us to handle problems in the order of a thousand billion mesh points in the four dimensions (3 dimensions of space plus time). The validity of the method has been thoroughly evaluated on many cases with known solutions: linear piston, scattering of plane wave by a heterogeneous sphere, propagation in a waveguide with a shear flow, scattering by a finite amplitude vortex and nonlinear propagation in a thermoviscous medium. This validation process allows for a detailed assessment of the advantages and limitations of the method. Finally, applications to atmospheric propagation of shock waves will be presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

3D Numerical Simulation of Kink-driven Rayleigh–Taylor Instability Leading to Fast Magnetic Reconnection

Fast magnetic reconnection involving non-MHD microscale physics is thought to underlie both solar eruptions and laboratory plasma current disruptions. While there is extensive research on both the MHD macroscale physics and the non-MHD microscale physics, the process by which large-scale MHD couples to the microscale physics is not well understood. An MHD instability cascade from a kink to a secondary Rayleigh–Taylor instability in the Caltech astrophysical jet laboratory experiment provides new insights into this coupling and motivates a 3D numerical simulation of this transition from large to small scale. A critical finding from the simulation is that the axial magnetic field inside the current-carrying dense plasma must exceed the field outside. In addition, the simulation verifies a theoretical prediction and experimental observation that, depending on the strength of the effective gravity produced by the primary kink instability, the secondary instability can be Rayleigh–Taylor or mini-kink. Finally, it is shown that the kink-driven Rayleigh–Taylor instability generates a localized electric field sufficiently strong to accelerate electrons to very high energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Regional-Scale 3D Ground-Motion Simulations of M w 7 Earthquakes on the Hayward Fault, Northern California Resolving Frequencies 0–10 Hz and Including Site-Response Corrections

Large earthquake ground-motion simulations in 3D Earth models provide constraints on site-specific shaking intensities but have suffered from limited frequency resolution and ignored site response in soft soils. In this work we report new regional-scale 3D simulations for moment magnitude 7.0 scenario earthquakes on the Hayward Fault, northern California with SW4. Simulations resolved significantly broader band frequencies (0–10 Hz) than previous studies and represent the highest resolution simulations for any such earthquake to date. Seismic waves were excited by a kinematic rupture following Graves and Pitarka (2016) and obeyed wave propagation in a 3D Earth model with topography from the U.S. Geological Survey (USGS) assuming a minimum shear wavespeed, VSmin, of 500 m/s. We corrected motions for linear and nonlinear site response for the shear wavespeed, VS, from the USGS 3D model, using a recently developed ground-motion model (GMM) for Fourier amplitude spectra . At soft soil locations subjected to strong shaking, the site-corrected intensities reflect the competing effects of linear amplification by low V S material, reduction of stiffness during nonlinear deformation, and damping of high frequencies. Sites with near-surface V S of 500 m/s or greater require no linear site correction but can experience amplitude reduction due to nonlinear response. Averaged over all sites, we obtained reasonable agreement with empirical ergodic median GMMs currently used for seismic hazard and design ground motions (epsilon less than 1), with marked improvement at soft sedimentary sites. At specific locations, the simulated shaking intensities show systematic differences from the GMMs that reveal path and site effects not captured in these ergodic models. Results suggest how next generation regional-scale earthquake simulations can provide higher spatial and frequency resolution while including effects of soft soils that are commonly ignored in scenario earthquake ground-motion simulations.

58 GEOSCIENCES↗

THREE-DIMENSIONAL CORE-COLLAPSE SUPERNOVA SIMULATED USING A 15 M{sub ⊙} PROGENITOR

We have performed ab initio neutrino radiation hydrodynamics simulations in three and two spatial dimensions (3D and 2D) of core-collapse supernovae from the same 15 M{sub ☉} progenitor through 440 ms after core bounce. Both 3D and 2D models achieve explosions; however, the onset of explosion (shock revival) is delayed by ∼100 ms in 3D relative to the 2D counterpart and the growth of the diagnostic explosion energy is slower. This is consistent with previously reported 3D simulations utilizing iron-core progenitors with dense mantles. In the ∼100 ms before the onset of explosion, diagnostics of neutrino heating and turbulent kinetic energy favor earlier explosion in 2D. During the delay, the angular scale of convective plumes reaching the shock surface grows and explosion in 3D is ultimately lead by a single, large-angle plume, giving the expanding shock a directional orientation not dissimilar from those imposed by axial symmetry in 2D simulations. We posit that shock revival and explosion in the 3D simulation may be delayed until sufficiently large plumes form, whereas such plumes form more rapidly in 2D, permitting earlier explosions.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling of lithium granule injection in NSTX using M3D-C1

In this paper, we present simulations of pedestal control by lithium granule injection (LGI) in NSTX. A model for small granule ablation has been implemented in the M3D-C1 code (Jardin et al 2012 Comput. Sci. Discovery 5 014002), allowing the simulation of realistic lithium granule injections. 2D and 3D simulations of Li injections in NSTX H-mode plasmas are performed and the effect of granule size, injection angle and velocity on the pedestal gradient increase is studied. The amplitude of the local pressure perturbation caused by the granules is found to be highly dependent on the solid granule size. Adjusting the granule injection velocity allows one to inject more particles at the pedestal top. 3D simulations show the destabilization of high order MHD modes whose amplitude is directly linked to the localized pressure perturbation, which is found to depend on the toroidal localization of the granule density source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A 2.5D Computational Method to Simulate Cylindrical Fluidized Beds

In this paper, the limitations of axisymmetric and Cartesian two-dimensional (2D) simulations of cylindrical gas-solid fluidized beds are discussed. A new method has been proposed to carry out pseudo-two-dimensional (2.5D) simulations of a cylindrical fluidized bed by appropriately combining computational domains of Cartesian 2D and axisymmetric simulations. The proposed method was implemented in the open-source code MFIX and applied to the simulation of a lab-scale bubbling fluidized bed with necessary sensitivity study. After a careful grid study to ensure the numerical results are grid independent, detailed comparisons of the flow hydrodynamics were presented against axisymmetric and Cartesian 2D simulations. Furthermore, the 2.5D simulation results have been compared to the three-dimensional (3D) simulation for evaluation. This new approach yields better agreement with the 3D simulation results than with axisymmetric and Cartesian 2D simulations.

Li, Tingwen↗

A computational model to generate simulated three-dimensional breast masses

Purpose: To develop algorithms for creating realistic three-dimensional (3D) simulated breast masses and embedding them within actual clinical mammograms. The proposed techniques yield high-resolution simulated breast masses having randomized shapes, with user-defined mass type, size, location, and shape characteristics. Methods: The authors describe a method of producing 3D digital simulations of breast masses and a technique for embedding these simulated masses within actual digitized mammograms. Simulated 3D breast masses were generated by using a modified stochastic Gaussian random sphere model to generate a central tumor mass, and an iterative fractal branching algorithm to add complex spicule structures. The simulated masses were embedded within actual digitized mammograms. The authors evaluated the realism of the resulting hybrid phantoms by generating corresponding left- and right-breast image pairs, consisting of one breast image containing a real mass, and the opposite breast image of the same patient containing a similar simulated mass. The authors then used computer-aided diagnosis (CAD) methods and expert radiologist readers to determine whether significant differences can be observed between the real and hybrid images. Results: The authors found no statistically significant difference between the CAD features obtained from the real and simulated images of masses with either spiculated or nonspiculated margins. Likewise, the authors found that expert human readers performed very poorly in discriminating their hybrid images from real mammograms. Conclusions: The authors’ proposed method permits the realistic simulation of 3D breast masses having user-defined characteristics, enabling the creation of a large set of hybrid breast images containing a well-characterized mass, embedded within real breast background. The computational nature of the model makes it suitable for detectability studies, evaluation of computer aided diagnosis algorithms, and teaching purposes.

60 APPLIED LIFE SCIENCES↗

Simulation-Guided 3D Nanomanufacturing via Focused Electron Beam Induced Deposition

Focused electron beam induced deposition (FEBID) is one of the few techniques that enables direct-write synthesis of free-standing 3D nanostructures. While the fabrication of simple architectures such as vertical or curving nanowires has been achieved by simple trial and error, processing complex 3D structures is not tractable with this approach. This is due, inpart, to the dynamic interplay between electron–solid interactions and the transient spatial distribution of absorbed precursor molecules on the solid surface. Here, we demonstrate the ability to controllably deposit 3D lattice structures at the micro/nanoscale, which have received recent interest owing to superior mechanical and optical properties. Moreover, a hybrid Monte Carlo–continuum simulation is briefly overviewed, and subsequently FEBID experiments and simulations are directly compared. Finally, a 3D computer-aided design (CAD) program is introduced, which generates the beam parameters necessary for FEBID by both simulation and experiment. In using this approach, we demonstrate the fabrication of various 3D lattice structures using Pt-, Au-, and W-based precursors.

3D nanoprinting↗

Edge ICRF simulations in 3D geometry: From MHD equilibrium to coupling determination

We present in this work a consistent numerical scheme that allows the computation of 3D magnetic fields a nd 3D density profiles and their usage in ion cyclotron range of frequencies (ICRF) coupling simulations. We first utilize the PARVMEC code to compute the 3D free-boundary plasma equilibrium in the ideal magnetohydrodynamic (MHD) approximation. Since the PARVMEC solution is only defined within the last closed flux surface (LCFS), the magnetic field domain is extended to the scrape-off layer (SOL) via the BMW code, which computes a divergence-free magnetic field solution arising from the external conductors’ vacuum field a nd t he PARVMEC flux surface currents. This magnetic reconstruction is then used in the EMC3-EIRENE transport code in order to compute 3D density profiles. In the last step, the RAPLICASOL code is utilized to compute the ICRF antenna S-matrices resulting from the 3D density profiles. We exemplify this scheme for the ASDEX Upgrade tokamak. A new implementation of a curved model for the ASDEX Upgrade ICRF 2-strap antenna in RAPLICASOL allows simulations in realistic geometry, without any coordinate transformations.

Lopez, Graham↗

High-energy synchrotron flares powered by strongly radiative relativistic magnetic reconnection: 2D and 3D PIC simulations

ABSTRACT The time evolution of high-energy synchrotron radiation generated in a relativistic pair plasma energized by reconnection of strong magnetic fields is investigated with 2D and 3D particle-in-cell (PIC) simulations. The simulations in this 2D/3D comparison study are conducted with the radiative PIC code OSIRIS, which self-consistently accounts for the synchrotron radiation reaction on the emitting particles, and enables us to explore the effects of synchrotron cooling. Magnetic reconnection causes compression of the plasma and magnetic field deep inside magnetic islands (plasmoids), leading to an enhancement of the flaring emission, which may help explain some astrophysical gamma-ray flare observations. Although radiative cooling weakens the emission from plasmoid cores, it facilitates additional compression there, further amplifying the magnetic field B and plasma density n, and thus partially mitigating this effect. Novel simulation diagnostics utilizing 2D histograms in the n-B space are developed and used to visualize and quantify the effects of compression. The n-B histograms are observed to be bounded by relatively sharp power-law boundaries marking clear limits on compression. Theoretical explanations for some of these compression limits are developed, rooted in radiative resistivity or 3D kinking instabilities. Systematic parameter-space studies with respect to guide magnetic field, system size, and upstream magnetization are conducted and suggest that stronger compression, brighter high-energy radiation, and perhaps significant quantum electrodynamic effects such as pair production, may occur in environments with larger reconnection-region sizes and higher magnetization, particularly when magnetic field strengths approach the critical (Schwinger) field, as found in magnetar magnetospheres.

Schoeffler, K. M. (ORCID:0000000276447381)↗

Black Hole to Photosphere: 3D GRMHD Simulations of Collapsars Reveal Wobbling and Hybrid Composition Jets

Long-duration γ-ray bursts (GRBs) accompany the collapse of massive stars and carry information about the central engine. However, no 3D models have been able to follow these jets from their birth via black hole (BH) to the photosphere. We present the first such 3D general-relativity magnetohydrodynamic simulations, which span over six orders of magnitude in space and time. The collapsing stellar envelope forms an accretion disk, which drags inwardly the magnetic flux that accumulates around the BH, becomes dynamically important, and launches bipolar jets. The jets reach the photosphere at ~10 12 cm with an opening angle θ j ~ 6° and a Lorentz factor Γ j ≲ 30, unbinding ≳90% of the star. We find that (i) the disk–jet system spontaneously develops misalignment relative to the BH rotational axis. As a result, the jet wobbles with an angle θ t ~ 12°, which can naturally explain quiescent times in GRB lightcurves. The effective opening angle for detection θ j + θ t suggests that the intrinsic GRB rate is lower by an order of magnitude than standard estimates. This suggests that successful GRBs are rarer than currently thought and emerge in only ~0.1% of supernovae Ib/c, implying that jets are either not launched or choked inside most supernova Ib/c progenitors. (ii) The magnetic energy in the jet decreases due to mixing with the star, resulting in jets with a hybrid composition of magnetic and thermal components at the photosphere, where ~10% of the gas maintains magnetization σ ≳ 0.1. This indicates that both a photospheric component and reconnection may play a role in the prompt emission.

79 ASTRONOMY AND ASTROPHYSICS↗