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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Production and Distribution of {sup 44}Ti and {sup 56}Ni in a Three-dimensional Supernova Model Resembling Cassiopeia A

The spatial and velocity distributions of nuclear species synthesized in the innermost regions of core-collapse supernovae can yield important clues about explosion asymmetries and the operation of the still disputed explosion mechanism. Recent observations of radioactive {sup 44}Ti with high-energy satellite telescopes ( Nuclear Spectroscopic Telescope Array [ NuSTAR ], INTEGRAL ) have measured gamma-ray line details, which provide direct evidence of large-scale explosion asymmetries in SN 1987A and in Cassiopeia A (Cas A) even by mapping of the spatial brightness distribution ( NuSTAR ). Here we discuss a 3D simulation of a neutrino-driven explosion, using a parameterized neutrino engine, whose {sup 44}Ti distribution is mostly concentrated in one hemisphere pointing opposite to the neutron star (NS) kick velocity. Both exhibit intriguing resemblance to the observed morphology of the Cas A remnant, although neither the progenitor nor the explosion was fine-tuned for a perfect match. Our results demonstrate that the asymmetries observed in this remnant can, in principle, be accounted for by a neutrino-driven explosion, and that the high {sup 44}Ti abundance in Cas A may be explained without invoking rapid rotation or a jet-driven explosion, because neutrino-driven explosions generically eject large amounts of high-entropy matter. The recoil acceleration of the NS is connected to mass ejection asymmetries and is opposite to the direction of the stronger explosion, fully compatible with the gravitational tugboat mechanism. Our results also imply that Cas A and SN 1987A could possess similarly “one-sided” Ti and Fe asymmetries, with the difference that Cas A is viewed from a direction with large inclination angle to the NS motion, whereas the NS in SN 1987A should have a dominant velocity component pointing toward us.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Gravitational Waves from Neutrino Emission Asymmetries in Core-collapse Supernovae

We present a broadband spectrum of gravitational waves (GWs) from core-collapse supernovae (CCSNe) sourced by neutrino emission asymmetries for a series of full 3D simulations. The associated GW strain probes the long-term secular evolution of CCSNe and small-scale turbulent activity and provides insight into the geometry of the explosion. For nonexploding models, both the neutrino luminosity and the neutrino gravitational waveform will encode information about the spiral SASI. The neutrino memory will be detectable for a wide range of progenitor masses for a galactic event. Our results can be used to guide near-future decihertz and long-baseline GW detection programs, including aLIGO, the Einstein Telescope, and DECIGO.

Astronomy & Astrophysics↗
DOE OSTI2021

Measuring an Off-center Detonation through Infrared Line Profiles: The Peculiar Type Ia Supernova SN 2020qxp/ASASSN-20jq

In this study, we present and analyze a near-infrared (NIR) spectrum of the underluminous Type Ia supernova SN 2020qxp/ASASSN-20jq obtained with NIRES at the Keck Observatory, 191 days after B-band maximum. The spectrum is dominated by a number of broad emission features, including the [Fe ιι] at 1.644 μm, which is highly asymmetric with a tilted top and a peak redshifted by ≈2000 km s -1 . In comparison with 2D non-LTE synthetic spectra computed from 3D simulations of off-center delayed-detonation Chandrasekhar-mass (M ch ) white dwarf (WD) models, we find good agreement between the observed lines and the synthetic profiles, and are able to unravel the structure of the progenitor’s envelope. We find that the size and tilt of the [Fe ιι] 1.644 μm profile (in velocity space) is an effective way to determine the location of an off-center delayed-detonation transition (DDT) and the viewing angle, and it requires a WD with a high central density of ~4 × 10 9 g cm -3 . We also tentatively identify a stable Ni feature around 1.9 μm characterized by a “pot-belly” profile that is slightly offset with respect to the kinematic center. In the case of SN 2020qxp/ASASSN-20jq, we estimate that the location of the DDT is ~0.3M WD off center, which gives rise to an asymmetric distribution of the underlying ejecta. We also demonstrate that low-luminosity and high-density WD SN Ia progenitors exhibit a very strong overlap of Ca and 56 Ni in physical space. This results in the formation of a prevalent [Ca ιι] 0.73 μm emission feature that is sensitive to asymmetry effects. Our findings are discussed within the context of alternative scenarios, including off-center C/O detonations in He-triggered sub-M Ch WDs and the direct collision of two WDs. Snapshot programs with Gemini/Keck/Very Large Telescope (VLT)/ELT-class instruments and our spectropolarimetry program are complementary to mid-IR spectra by the James Webb Space Telescope (JWST).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

What Sets the Star Formation Rate of Molecular Clouds? The Density Distribution as a Fingerprint of Compression and Expansion Rates

We use a suite of 3D simulations of star-forming molecular clouds, with and without stellar feedback, magnetic fields, and driven turbulence, to study the compression and expansion rates of the gas as functions of density. We show that, around the mean density, supersonic turbulence promotes rough equilibrium between the amounts of compressing and expanding gas, consistent with continuous gas cycling between high- and low-density states. We find that the inclusion of protostellar jets produces rapidly expanding and compressing low-density gas. We find that the gas mass flux peaks at the transition between the lognormal and power-law forms of the density probability distribution function (PDF). This is consistent with the transition density tracking the post-shock density, which promotes an enhancement of mass at this density (i.e., shock compression and filament formation). At high densities, the gas dynamics are dominated by self-gravity: the compression rate in all of our runs matches the rate of the run with only gravity, suggesting that processes other than self-gravity have little effect at these densities. The net gas mass flux becomes constant at a density below the sink formation threshold, where it equals the star formation rate. The density at which the net gas mass flux equals the star formation rate is one order of magnitude lower than our sink threshold density, corresponds to the formation of the second power-law tail in the density PDF, and sets the overall star formation rates of these simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

On a Spectral Method for β -particle Bound Excitation Collisions in Kilonovae

The interaction of β-particles with the weakly ionized plasma background is an important mechanism for powering the kilonova (KN) transient signal from neutron star mergers. For this purpose, we present an implementation of the approximate fast-particle collision kernel, described by Inokuti following the seminal formulation of Bethe, in a spectral solver of the Vlasov–Maxwell–Boltzmann equation. In particular, we expand the fast-particle plane-wave atomic excitation kernel into coefficients of the Hermite basis, and derive the relevant discrete spectral system. In this fast-particle limit, the approach permits the direct use of atomic data, including optical oscillator strengths, normally applied to photon–matter interaction. The resulting spectral matrix is implemented in the MASS-APP spectral solver framework, in a way that avoids full matrix storage per spatial zone. We numerically verify aspects of the matrix construction, and present a proof-of-principle 3D simulation of a 2D axisymmetric KN ejecta snapshot. Our preliminary numerical results indicate that a reasonable choice of Hermite basis parameters for β-particles in the KN is a bulk velocity parameter u = 0, a thermal velocity parameter α = 0.5c, and a 9 × 9 × 9 mode velocity basis set (Hermite orders of 0–8 in each dimension). For interior-ejecta sample zones, we estimate that the ratio of thermalization from large-angle (≳2fdg5) bound excitation scattering to total thermalization is ~0.002–0.003.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Modeling Solids in Nuclear Astrophysics with Smoothed Particle Hydrodynamics

Abstract Smoothed particle hydrodynamics (SPH) is a frequently applied tool in computational astrophysics to solve the fluid dynamics equations governing the systems under study. For some problems, for example when involving asteroids and asteroid impacts, the additional inclusion of material strength is necessary in order to accurately describe the dynamics. In compact stars, that is white dwarfs and neutron stars, solid components are also present. Neutron stars have a solid crust, which is the strongest material known in nature. However, their dynamical evolution, when modeled via SPH or other computational fluid dynamics codes, is usually described as a purely fluid dynamics problem. Here, we present the first 3D simulations of neutron star crustal toroidal oscillations including material strength with the Los Alamos National Laboratory SPH code FleCSPH. In the first half of the paper, we present the numerical implementation of solid material modeling together with standard tests. The second half is on the simulation of crustal oscillations in the fundamental toroidal mode. Here, we dedicate a large fraction of the paper to approaches that can suppress numerical noise in the solid. If not minimized, the latter can dominate the crustal motion in the simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1989

Field structure of collapsing wave packets in 3D strong Langmuir turbulence

A simple model is constructed for the electric fields in the collapsing wave packets found in 3D simulations of driven and damped isotropic strong Langmuir turbulence. This model, based on a spherical-harmonic decomposition of the electrostatic potential, accounts for the distribution of wave-packet shapes observed in the simulations, particularly the predominance of oblate wave packets. In contrast with predictions for undamped and undriven subsonic collapse of scalar fields, oblate vector-field wave packets do not flatten during collapse but, instead, remain approximately self-similar and rigid.

Newman, D. L.↗
NASA NTRS1990

Fluid dynamics and low gravity effects of chemical vapor deposition

Based on the comparison between experimental data and numerical results for the growth of GaAs from TMGa, it was shown that 3D simulations are necessary to simulate rectangular CVD reactors even when operated under subcritical (Ra) conditions. The important points found are summarized in the three attached reprints. The experimental studies of mixed convection in horizontal channels have shown three regimes of high Ra (22,220) number flows. At Re = 18.5, the rolls develop very quickly, significantly modulating the axial velocity even before it reaches the beginning of the hot plate. A few centimeters downstream, the velocities become asymmetric about the vertical centerplane and at x = 12 cm, become unsteady. These asymmetries were predicted theoretically, but experimental evidence has not been published prior to this work. At Re = 36, the axial velocity is only slightly modified at x = 0. Although the flow remains steady and symmetric about the vertical centerplane, there is a small spatial oscillation in the velocities over the length of the channel. The period of this oscillation was around 5 cm. At Re = 54, the longitudinal rolls developed smoothly over a length of 30 cm, with no asymmetries, unsteadiness, or spatial oscillations. Comparison of numerical simulations of these flows to experiments has revealed the importance and difficulty of setting proper thermal boundary conditions on the sidewalls. Calculated flows and experimentally measured flows showed very similar profiles, but at different axial locations, with the rolls developing more rapidly in the experiments. This is directly attributable to partially conducting sidewalls of the apparatus being hotter in the entrance section than the adiabatic walls of the simulations. A thorough comparison of the experimental data and numerical results for a variety of sidewall boundary conditions is in preparation.

Nyce, Thomas A.↗
NASA NTRS1991

Short-time Lyapunov exponent analysis and the transition to chaos in Taylor-Couette flow

The physical mechanism driving the weakly chaotic Taylor-Couette flow is investigated using the short-time Liapunov exponent analysis. In this procedure, the transition from quasi-periodicity to chaos is studied using direct numerical 3D simulations of axially periodic Taylor-Couette flow, and a partial Liapunov exponent spectrum for the flow is computed by simultaneously advancing the full solution and a set of perturbations. It is shown that the short-time Liapunov exponent analysis yields more information on the exponents and dimension than that obtained from the common Liapunov exponent calculations. Results show that the chaotic state studied here is caused by a Kelvin-Helmholtz-type instability of the outflow boundary jet of Taylor vortices.

Vastano, John A.↗
NASA NTRS1992

Timing considerations of Helmet Mounted Display performance

The Helmet Mounted Display (HMD) system developed in our lab should be a useful teleoperator systems display if it increases operator performance of the desired task; it can, however, introduce degradation in performance due to display update rate constraints and communication delays. Display update rates are slowed by communication bandwidth and/or computational power limitations. We used simulated 3D tracking and pick-and-place tasks to characterize performance levels for a range of update rates. Initial experiments with 3D tracking indicate that performance levels plateau at an update rate between 10 and 20 Hz. We have found that using the HMD with delay decreases performance as delay increases.

Tharp, Gregory↗
NASA NTRS1993

On the generation of sound by turbulent convection. I - A numerical experiment

Motivated by the problem of the origin of the solar p-modes, we study the generation of acoustic waves by turbulent convection. Our approach uses the results of high-resolution 3D simulations as the experimental basis for our investigation. The numerical experiment describes the evolution of a horizontally periodic layer of vigorously convecting fluid. The sound is measured by a procedure, based on a suitable linearization of the equations of compressible convection that allows the amplitude of the acoustic field to be determined. Through this procedure we identify unambiguously some 400 acoustic modes. The total energy of the acoustic field is found to be a fraction of a percent of the kinetic energy of the convection. The amplitudes of the observed modes depend weakly on (horizontal) wavenumber but strongly on frequency. The line widths of the observed modes typically exceed the natural linewidths of the modes as inferred from linear theory. This broadening appears to be related to the (stochastic) interaction between the modes and the underlying turbulence which causes abrupt, episodic events during which the phase coherence of the modes is lost.

Bogdan, Thomas J.↗
NASA NTRS1993

Line relaxation methods for the solution of 2D and 3D compressible flows

An implicit finite element based algorithm for the compressible Navier-Stokes equations is outlined, and the solution of the resulting equation by a line relaxation on general meshes of triangles or tetrahedra is described. The problem of generating and adapting unstructured meshes for viscous flows is reexamined, and an approach for both 2D and 3D simulations is proposed. An efficient approach appears to be the use of an implicit/explicit procedure, with the implicit treatment being restricted to those regions of the mesh where viscous effects are known to be dominant. Numerical examples demonstrating the computational performance of the proposed techniques are given.

Hassan, O.↗
NASA NTRS1992

Finite-element solution scheme for compressible viscous flow

Adaptivity offers the possibility of using FEM for problems that are otherwise intractable. The solution algorithm is combined with an adaptive approach in which the unstructured portion of the grid is adaptively regenerated as the solution proceeds; the structured portion of the mesh is modified accordingly, to ensure both compatibility along the interface between the implicit and explicit element groups and increased resolution of the viscous layer. The approach offers significant computational savings in 3D simulations.

Hassan, O.↗
NASA NTRS1993

Percolation technique for galaxy clustering

We study percolation in mass and galaxy distributions obtained in 3D simulations of the CDM, C + HDM, and the power law (n = -1) models in the Omega = 1 universe. Percolation statistics is used here as a quantitative measure of the degree to which a mass or galaxy distribution is of a filamentary or cellular type. The very fast code used calculates the statistics of clusters along with the direct detection of percolation. We found that the two parameters mu(infinity), characterizing the size of the largest cluster, and mu-squared, characterizing the weighted mean size of all clusters excluding the largest one, are extremely useful for evaluating the percolation threshold. An advantage of using these parameters is their low sensitivity to boundary effects. We show that both the CDM and the C + HDM models are extremely filamentary both in mass and galaxy distribution. The percolation thresholds for the mass distributions are determined.

Klypin, Anatoly↗
NASA NTRS1998

A Parameterization for the Triggering of Landscape Generated Moist Convection

A set of relatively high resolution three-dimensional (3D) simulations were produced to investigate the triggering of moist convection by landscape generated mesoscale circulations. The local accumulated rainfall varied monotonically (linearly) with the size of individual landscape patches, demonstrating the need to develop a trigger function that is sensitive to the size of individual patches. A new triggering function that includes the effect of landscapes generated mesoscale circulations over patches of different sizes consists of a parcel's perturbation in vertical velocity (nu(sub 0)), temperature (theta(sub 0)), and moisture (q(sub 0)). Each variable in the triggering function was also sensitive to soil moisture gradients, atmospheric initial conditions, and moist processes. The parcel's vertical velocity, temperature, and moisture perturbation were partitioned into mesoscale and turbulent components. Budget equations were derived for theta(sub 0) and q(sub 0). Of the many terms in this set of budget equations, the turbulent, vertical flux of the mesoscale temperature and moisture contributed most to the triggering of moist convection through the impact of these fluxes on the parcel's temperature and moisture profile. These fluxes needed to be parameterized to obtain theta(sub 0) and q(sub 0). The mesoscale vertical velocity also affected the profile of nu(sub 0). We used similarity theory to parameterize these fluxes as well as the parcel's mesoscale vertical velocity.

Lynn, Barry H.↗
NASA NTRS2004

Starting Conditions for Hydrothermal Systems Underneath Martian Craters: Hydrocode Modeling

Mars is the most Earth-like of the Solar System s planets, and the first place to look for any sign of present or past extraterrestrial life. Its surface shows many features indicative of the presence of surface and sub-surface water, while impact cratering and volcanism have provided temporary and local surface heat sources throughout Mars geologic history. Impact craters are widely used ubiquitous indicators for the presence of sub-surface water or ice on Mars. In particular, the presence of significant amounts of ground ice or water would cause impact-induced hydrothermal alteration at Martian impact sites. The realization that hydrothermal systems are possible sites for the origin and early evolution of life on Earth has given rise to the hypothesis that hydrothermal systems may have had the same role on Mars. Rough estimates of the heat generated in impact events have been based on scaling relations, or thermal data based on terrestrial impacts on crystalline basements. Preliminary studies also suggest that melt sheets and target uplift are equally important heat sources for the development of a hydrothermal system, while its lifetime depends on the volume and cooling rate of the heat source, as well as the permeability of the host rocks. We present initial results of two-dimensional (2D) and three-dimensional (3D) simulations of impacts on Mars aimed at constraining the initial conditions for modeling the onset and evolution of a hydrothermal system on the red planet. Simulations of the early stages of impact cratering provide an estimate of the amount of shock melting and the pressure-temperature distribution in the target caused by various impacts on the Martian surface. Modeling of the late stage of crater collapse is necessary to characterize the final thermal state of the target, including crater uplift, and distribution of the heated target material (including the melt pool) and hot ejecta around the crater.

Pierazzo, E.↗
NASA NTRS2009

Noise Radiation From a Leading-Edge Slat

This paper extends our previous computations of unsteady flow within the slat cove region of a multi-element high-lift airfoil configuration, which showed that both statistical and structural aspects of the experimentally observed unsteady flow behavior can be captured via 3D simulations over a computational domain of narrow spanwise extent. Although such narrow domain simulation can account for the spanwise decorrelation of the slat cove fluctuations, the resulting database cannot be applied towards acoustic predictions of the slat without invoking additional approximations to synthesize the fluctuation field over the rest of the span. This deficiency is partially alleviated in the present work by increasing the spanwise extent of the computational domain from 37.3% of the slat chord to nearly 226% (i.e., 15% of the model span). The simulation database is used to verify consistency with previous computational results and, then, to develop predictions of the far-field noise radiation in conjunction with a frequency-domain Ffowcs-Williams Hawkings solver.

Lockard, David P.↗
NASA NTRS2007

Detailed Studies of Pixelated CZT Detectors Grown with the Modified Horizontal Bridgman Method

The detector material Cadmium Zinc Telluride (CZT) achieves excellent spatial resolution and good energy resolution over a broad energy range, several keV up to some MeV. Presently, there are two main methods to grow CZT crystals, the Modified High-Pressure Bridgman (MHB) and the High-Pressure Bridgman (HPB) process. The study presented in this paper is based on MHB CZT substrates from the company Orbotech Medical Solutions Ltd. [Orbotech Medical Solutions Ltd., 10 Plaut St., Park Rabin, P.O. Box 2489, Rehovot, Israel, 76124]. Former studies have shown that high-work-function materials on the cathode side reduce the leakage current and, therefore, improve the energy resolution at lower energies. None of the studies have emphasized on the anode contact material. Therefore, we present in this paper the result of a detailed study in which for the first time the cathode material was kept constant and the anode material was varied. We used four different anode materials: Indium, Titanium, Chromium and Gold, metals with work-functions between 4.1 eV and 5.1 eV. The detector size was 2.0 x 2.0 x 0.5 cu cm with 8 x 8 pixels and a pitch of 2.46 mm. The best performance was achieved with the low-work-function materials Indium and Titanium with energy resolutions of 2.0 keV (at 59 keV) and 1.9 keV (at 122 keV) for Titanium and 2.1 keV (at 59 keV) and 2.9 keV (at 122 keV) for Indium. Taking into account the large pixel pitch of 2.46 mm, these resolutions are very competitive in comparison to those achieved with detectors made of material produced with the more expensive conventional HPB method. We present a detailed comparison of our detector response with 3D simulations. The latter comparisons allow us to determine the mobility-lifetime-products (mu tau-products) for electrons and holes. Finally, we evaluated the temperature dependency of the detector performance and ls-products. For many applications temperature dependence is important, therefore, we extended the scope of our study to temperatures as low as -30 C. There are two important results. The breakdown voltage increases with decreasing temperature, and electron mobility-lifetime-product decreases by about 30% over a range from 20 C to -30 C. The latter effect causes the energy resolution to deteriorate, but the concomitantly increasing breakdown voltage makes it possible to increase the applied bias voltage and restore the full performance

Jung, I.↗