application of 1D stagnation solutions to 3D wire array simulation in absence of radiation.
Abstract not provided.
Yu, Edmund P.↗
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.
Effects of Vacuum Impedance on MITL Flow Using 3D Electromagnetic PIC Simulations.
Abstract not provided.
Powell, Troy C.↗
EFFECTS OF VACUUM IMPEDANCE CHANGES ON MITL FLOW USING 3D ELECTROMAGNETIC PIC SIMULATIONS.
Abstract not provided.
Powell, Troy C.↗
ElRock-Net: Assessing the Utility of Machine Learning to Initialize 3D Electric Potential Simulations
Explore the source record for details and available documents.
Earth Sciences↗
Application of one-dimensional stagnation solutions to three-dimensional simulation of compact wire array in absence of radiation
We investigate the stagnation phase of a three-dimensional (3D), magnetohydrodynamic simulation of a compact, tungsten wire-array Z pinch, under the simplifying assumption of negligible radiative loss. In particular, we address the ability of one-dimensional (1D) analytic theory to describe the time evolution of spatially averaged plasma properties from 3D simulation. The complex fluid flows exhibited in the stagnated plasma are beyond the scope of 1D theory and result in centrifugal force as well as enhanced thermal transport. Despite these complications, a 1D homogeneous (i.e., shockless) stagnation solution can capture the increase of on-axis density and pressure during the initial formation of stagnated plasma. Later, when the stagnated plasma expands outward into the imploding plasma, a 1D shock solution describes the decrease of on-axis density and pressure, as well as the growth of the shock accretion region.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications
In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.
Bravo, Ramiro H.↗
Dendritic Growth Under Natural and Forced Convection in Al-Cu Alloys: From Equiaxed to Columnar Dendrites and from 2D to 3D Phase-Field Simulations
The interaction between convection and solute transport during solidification has significant influence on the dendritic evolution. By employing the phase-field lattice-Boltzmann approach together with the parallel and adaptive-mesh-refinement algorithm, the dendritic evolution under convection is simulated in both 2D and 3D cases. The flow-induced redistribution of the solute alters both tip velocity and the development of dendritic arms. The effect of both convection and undercooling is quantified and compared using the length ratio of the dendritic arms. The effect of convection behavior (i.e., natural and forced) and domain dimension (i.e., 2D and 3D) on dendritic growth is discussed. Results show that the convection effect is mainly dominated by the convection mode, and the melt flow in 2D can produce biased results comparing with those in 3D.
36 MATERIALS SCIENCE↗
Toward 2D Dynamo Models Calibrated by Global 3D Relativistic Accretion Disk Simulations
Two-dimensional models assuming axisymmetry are an economical way to explore the long-term evolution of black hole accretion disks, but they are only realistic if the feedback of the nonaxisymmetric turbulence on the mean momentum and magnetic fields is incorporated. Dynamo terms added to the 2D induction equation should be calibrated to 3D MHD simulations. For generality, the dynamo tensors should be calibrated as functions of local variables rather than explicit functions of spatial coordinates in a particular basis. In this paper, we study the feedback of non-axisymmetric features on the 2D mean fields using a global 3D, relativistic, Cartesian simulation from the Illinois GRMHD code. We introduce new methods for estimating overall dynamo alpha and turbulent diffusivity effects as well as measures of the dominance of non-axisymmetric components of energies and fluxes within the disk interior. We attempt closure models of the dynamo EMF using least squares fitting, considering both models where coefficient tensors are functions of space and more global, covariant models. None of these models are judged satisfactory, but we are able to draw conclusions on what sorts of generalizations are and are not promising.
Matthew D Duez↗
3D Convective Urca Process in a Simmering White Dwarf
Abstract A proposed setting for thermonuclear (Type Ia) supernovae is a white dwarf that has gained mass from a companion to the point of carbon ignition in the core. In the early stages of carbon burning, called the simmering phase, energy released by the reactions in the core drive the formation and growth of a core convection zone. One aspect of this phase is the convective Urca process, a linking of weak nuclear reactions to convection, which may alter the composition and structure of the white dwarf. The convective Urca process is not well understood and requires 3D fluid simulations to properly model the turbulent convection, an inherently 3D process. Because the neutron excess of the fluid both sets and is set by the extent of the convection zone, the realistic steady state can only be determined in simulations with real 3D mixing processes. Additionally, the convection is relatively slow (Mach number less than 0.005) and thus a low Mach number method is needed to model the flow over many convective turnovers. Using the MAESTROeX low Mach number hydrodynamic software, we present the first full-star 3D simulations of the A = 23 convective Urca process, spanning hundreds of convective turnover times. Our findings on the extent of mixing across the Urca shell, the characteristic velocities of the flow, the energy-loss rates due to neutrino emission, and the structure of the convective boundary can be used to inform 1D stellar models that track the longer-timescale evolution.
Boyd, Brendan (ORCID:0000000254199751)↗
Three-dimensional Hydrodynamic Simulations of Convective Nuclear Burning in Massive Stars Near Iron Core Collapse
Nonspherical structure in massive stars at the point of iron core collapse can have a qualitative impact on the properties of the ensuing core-collapse supernova explosions and the multimessenger signals they produce. Strong perturbations can aid successful explosions by strengthening turbulence in the postshock region. Here we report on a set of 4π 3D hydrodynamic simulations of O- and Si-shell burning in massive star models of varied initial masses using MESA and the FLASH simulation framework. We evolve four separate 3D models for roughly the final 10 minutes prior to and including iron core collapse. We consider initial 1D MESA models with masses of 14, 20, and 25 M ⊙ to survey a range of O/Si-shell density and compositional configurations. Further, we characterize the convective shells in our 3D models and compare them to the corresponding 1D models. In general, we find that the angle-average convective speeds in our 3D simulations near collapse are three to four times larger than the convective speeds predicted by MESA at the same epoch for our chosen mixing length parameter of α MLT = 1.5. In three of our simulations, we observe significant power in the spherical harmonic decomposition of the radial velocity field at harmonic indices of ℓ = 1–3 near collapse. Our results suggest that large-scale modes are common in massive stars near collapse and should be considered a key aspect of presupernova progenitor models.
79 ASTRONOMY AND ASTROPHYSICS↗
Instabilities and turbulence in low-$β$ guide field reconnection exhausts with kinetic Riemann simulations
The role of turbulence in low-$β$, guide-field reconnection exhausts is explored in 2D reconnection and 2D and 3D Riemann simulations. The structure of the exhaust and associated turbulence is controlled by a pair of rotational discontinuities (RDs) at the exhaust boundary and a pair of slow shocks (SSs) that are generated by counterstreaming ions beams. In 2D, the exhaust develops large-amplitude striations at the ion Larmor radius scale that are produced by electron-beam-driven ion cyclotron waves. The electron beams driving the instability are injected into the exhaust from one of the RDs. However, in 3D Riemann simulations, the additional dimension (in the out-of-plane direction) results in strong Buneman and electron-electron streaming instabilities at the RD which suppress the electron beam formation and therefore the striations in the exhaust. The strength of the streaming instabilities at the RD is controlled by the ratio of the electron thermal speed to Alfvén speed, with the lower thermal speed being more unstable. In the 3D simulations, an ion-ion streaming instability acts to partially thermalize the counterstreaming ion beams at the SSs. This instability is controlled by the ratio of the sound speed to Alfvén speed and is expected to be stable in the low $β$ solar corona. The results suggest that in a guide field reconnection exhaust with 1 $\gg$ $β >$ $m_e/m_i$, the kinetic-scale turbulence that develops will be too weak to play a significant role in energy conversion and particle acceleration. Therefore, the energy conversion will be mostly controlled by laminar physics or multi-x-line reconnection.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A NEW MODEL FOR MIXING BY DOUBLE-DIFFUSIVE CONVECTION (SEMI-CONVECTION). III. THERMAL AND COMPOSITIONAL TRANSPORT THROUGH NON-LAYERED ODDC
Oscillatory double-diffusive convection (ODDC; also known as semi-convection) refers to a type of double-diffusive instability that occurs in regions of planetary and stellar interiors that have a destabilizing thermal stratification and a stabilizing mean molecular weight stratification. In this series of papers, we use an extensive suite of three-dimensional (3D) numerical simulations to quantify the transport of heat and chemical species by ODDC. Rosenblum et al. first showed that ODDC can either spontaneously form layers that significantly enhance the transport of heat and chemical species compared to microscopic transport or remain in a state dominated by large-scale gravity waves, in which there is a more modest enhancement of the turbulent transport rates. Subsequent studies in this series focused on identifying under what conditions layers form and quantifying transport through layered systems. Here we proceed to characterize transport through systems that are unstable to ODDC, but do not undergo spontaneous layer formation. We measure the thermal and compositional fluxes in non-layered ODDC from both two-dimensional (2D) and 3D numerical simulations, and show that 3D simulations are well approximated by similar simulations in a 2D domain. We find that the turbulent mixing rate in this regime is weak and can, to a first-level approximation, be neglected. We conclude by summarizing the findings of papers I through III into a single prescription for transport systems unstable to ODDC.
79 ASTRONOMY AND ASTROPHYSICS↗
Understanding asymmetries using integrated simulations of capsule implosions in low gas-fill hohlraums at the National Ignition Facility
Current capsule implosions at the National Ignition Facility (NIF) using high-density-carbon ablators and laser energies close to 2 MJ have shown neutron yields in excess of 50 kJ. Improving on this performance requires understanding of the different degradation mechanisms. For many NIF implosions, nuclear diagnostic signatures have inferred residual hot-spot velocities that correlate with fuel areal density variations consistent with a low-order mode-1 asymmetry. A current working hypothesis attributes these asymmetries to a combination of beam to beam variations in the laser delivery and possibly coupling to target features, such as the diagnostic holes needed for x-ray imaging. Recently, a new source has been identified, thickness variation in the ablator shell. To gain better understanding and eventually mitigate the causes of <ρR> asymmetries, 3D integrated simulations using the actual delivered laser powers are needed. To capture the effect of the diagnostic holes (DHs) using direct numerical simulation would require significantly large computational resources. Instead, our 3D simulations make use of a subgrid model developed using highly-resolved 2D simulations that include several details of the DH engineering complexity. Simulations of NIF shots using hohlraums without DHs, to isolate the effect of beam-to-beam variations, reproduce fairly well the observed nuclear diagnostic signatures. Similarly, reasonable agreement between data and simulations is also obtained in the presence of diagnostic holes. To account for the remaining discrepancies a sensitivity study of ablator thickness variation showed that 1% thickness asymmetries are comparable in effect to 1% peak drive mode-1 asymmetries. Furthermore, this study identified sensitivity to variations in the imbedded doped layer (needed to shield the DT ice from the high energy x-rays generated in the hohlraum) thickness even when the inner and outer surface of the ablator are perfectly spherical.
radiation hydrodynamics↗
Progress on the Multiphysics Capabilities of the Parallel Electromagnetic ACE3P Simulation Suite
ACE3P is a 3D parallel simulation suite that is being developed at SLAC National Accelerator Laboratory. Effectively utilizing supercomputer resources, ACE3P has become a key tool for the coupled electromagnetic, thermal and mechanical research and design of particle accelerators. Based on the existing finite-element infrastructure, a massively parallel eigensolver is developed for modal analysis of mechanical structures. It complements a set of the multiphysics tools in ACE3P and, in particular, can be used for the comprehensive study of microphonics in accelerating cavities ensuring the operational reliability of a particle accelerator.
Kononenko, Oleksiy↗
An Efficient Three-Dimensional CFD-Based Numerical Wave Tank for a Wave Energy Converter in Extreme Irregular Waves
A numerical wave tank approach for computational fluid dynamics (CFD) modelling of an extreme irregular seastate is presented. The technique couples a potential flow solution with a CFD solver for more efficient numerical predictions. This method has recently become attractive both for the research community and the industry working with offshore structures. The model is used to determine the response of a submerged pressure differential wave energy converter (WEC) in a fully nonlinear irregular wave condition using the high-fidelity CFD code, STAR-CCM+. Potential flow based numerical models are commonly used to predict motions and performance of wave energy converters. Wave kinematics can deviate from potential flow predictions for extreme wave conditions; the excitation loads on an absorber can also be increasingly influenced by viscous effects, not predicted by potential flow engineering level models. In these extreme conditions, a Reynolds-averaged Navier-Stokes CFD model can better predict motions and loads for a WEC. Long time series with varying random seed numbers can be used to identify singular extreme wave events from a stochastic irregular sea state. This approach simulates a more realistic wave series for a given sea state than a regular wave or a focused wave. However, it is computationally infeasible to run these long time series for three-dimensional (3D) CFD simulations. In this work, two-dimensional (2D) CFD simulations with a long domain allow the full development of an extreme nonlinear wave condition. The results are used to identify extreme events from a 50-year storm condition for the PacWave site off the coast of Oregon. A relatively short time window including this extreme event is then mapped to a 3D simulation using a user defined wave methodology. Convergence studies for domain length, wave forcing lengths, and time before the extreme event were conducted.
CFD↗
An Efficient Three-Dimensional CFD-Based Numerical Wave Tank for a Wave Energy Converter in Extreme Irregular Waves: Preprint
A numerical wave tank approach for computational fluid dynamics (CFD) modelling of an extreme irregular seastate is presented. The technique couples a potential flow solution with a CFD solver for more efficient numerical predictions. This method has recently become attractive both for the research community and the industry working with offshore structures. The model is used to determine the response of a submerged pressure differential wave energy converter (WEC) in a fully nonlinear irregular wave condition using the high fidelity CFD code, STAR-CCM+. Potential flow based numerical models are commonly used to predict motions and performance of wave energy converters. Wave kinematics can deviate from potential flow predictions for extreme wave conditions; the excitation loads on an absorber can also be increasingly influenced by viscous effects, not captured in engineering level models. In these extreme conditions, a Reynolds-averaged Navier-Stokes CFD model can better predict motions and loads for a WEC. Long time series with varying random seed numbers can be used to identify singular extreme wave events from a stochastic irregular sea state. This approach simulates a more realistic wave series for a given sea state than a regular wave or a focused wave. However, it is computationally infeasible to run these long time series for three-dimensional (3D) CFD simulations. In this work, two-dimensional (2D) CFD simulations with a long domain allow the full development of an extreme nonlinear wave condition. The results are used to identify extreme events from a 50-year storm condition for the PacWave site off the coast of Oregon. A relatively short time window including this extreme event is then mapped to a 3D simulation using a user defined wave methodology. Convergence studies for domain length, wave forcing lengths, and time before the extreme event were conducted.
cfd↗
Nucleosynthetic Analysis of Three-dimensional Core-collapse Supernova Simulations
Abstract We study in detail the ejecta conditions and theoretical nucleosynthetic results for 18 three-dimensional core-collapse supernova (CCSN) simulations done by F ornax . Most of the simulations are carried out to at least 3 s after bounce, which allows us to follow their longer-term behaviors. We find that multidimensional effects introduce many complexities into the ejecta conditions. We see a stochastic electron fraction evolution, complex peak temperature distributions and histories, and long-tail distributions of the time spent within nucleosynthetic temperature ranges. These all lead to substantial variation in CCSN nucleosynthetic yields and differences from 1D results. We discuss the production of lighter α -nuclei, radioactive isotopes, heavier elements, and a few isotopes of special interest. Comparing pre-CCSN and CCSN contributions, we find that a significant fraction of elements between roughly Si and Ge are generally produced in CCSNe. We find that 44 Ti exhibits an extended production timescale as compared to 56 Ni, which may explain its different distribution and higher than previously predicted abundances in supernova remnants such as Cas A and SN1987A. We also discuss the morphology of the ejected elements. This study highlights the high-level diversity of ejecta conditions and nucleosynthetic results in 3D CCSN simulations and emphasizes the need for additional long-term (∼10 s) 3D simulations to properly address such complexities.
Astronomy & Astrophysics↗