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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Azimuthal and axial structures in 3D particle-in-cell simulation of Penning discharge

We report the results of a 3D particle-in-cell simulation of cylindrical Penning discharge with the axial magnetic field in the so-called reflex configuration, where the cathode and anti-cathode are biased to the same negative potential. The discharge is supported by thermal electron emission from the cathode. Electron and ion collisions, including ionization, are fully accounted for. The emphasis is on a specific regime in which the plasma potential at the center of the discharge is positive with respect to the chamber walls, serving as an anode. Spatial and temporal scales of the observed azimuthal and axial fluctuations and structures are characterized. It is suggested that azimuthal structures are caused by the dissipative gradient-drift instability. We find that the axial fluctuations related to the plasma-beam instabilities are weakly correlated with the azimuthal perturbations of the density, so that the azimuthal modes rotate as a whole and do not show any axial shear. The behavior of the electric potential is more involved, demonstrating intermittent standing wave and propagating structures in the axial direction that modulate the electron transport, producing the standing wave pattern (along the z-direction) in the radial electron flux.

3D PIC simulation of discharge↗
DOE OSTI2016

QUIESCENT PROMINENCES IN THE ERA OF ALMA: SIMULATED OBSERVATIONS USING THE 3D WHOLE-PROMINENCE FINE STRUCTURE MODEL

We use the detailed 3D whole-prominence fine structure model to produce the first simulated high-resolution ALMA observations of a modeled quiescent solar prominence. The maps of synthetic brightness temperature and optical thickness shown in the present paper are produced using a visualization method for synthesis of the submillimeter/millimeter radio continua. We have obtained the simulated observations of both the prominence at the limb and the filament on the disk at wavelengths covering a broad range that encompasses the full potential of ALMA. We demonstrate here extent to which the small-scale and large-scale prominence and filament structures will be visible in the ALMA observations spanning both the optically thin and thick regimes. We analyze the relationship between the brightness and kinetic temperature of the prominence plasma. We also illustrate the opportunities ALMA will provide for studying the thermal structure of the prominence plasma from the cores of the cool prominence fine structure to the prominence–corona transition region. In addition, we show that detailed 3D modeling of entire prominences with their numerous fine structures will be important for the correct interpretation of future ALMA observations of prominences.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Quantum mechanical model of crossing and anti-crossing points in 3D full-band Monte Carlo simulations

This work presents a 3D quantum mechanics based model to address the physics at band structure crossing/anti-crossing points in full band Monte Carlo (FBMC) simulations. The model solves the Krieger and Iafrate (KI) equations in real time using pre-computed coefficients at k-points spatially sampled within the first Brillouin zone. Solving the KI equations in real time makes this model applicable for all electric fields, which enables its use in FBMC device simulations. In this work, a two-level refinement scheme is used to aggressively sample regions in proximity to band crossings for accurate solutions to the KI equations and coarsely sample everywhere else to limit the number of k-points used. The presented sampling method is demonstrated on the band structure of silicon but is effective for the band structure of any semiconductor material. Next, the adaptation of the fully quantum KI model into the semi-classical FBMC method is discussed. Finally, FBMC simulations of hole transport in 4H silicon carbide with and without the KI model are performed. Results along different crystallographic directions for a wide range of electric fields are compared to previously published simulation and experimental values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

3D Radiation Nonideal Magnetohydrodynamical Simulations of the Inner Rim in Protoplanetary Disks

Many planets orbit within 1 au of their stars, raising questions about their origins. Particularly puzzling are the planets found near the silicate sublimation front. We investigate conditions near the front in the protostellar disk around a young intermediate-mass star, using the first global 3D radiation nonideal MHD simulations in this context. We treat the starlight heating; the silicate grains’ sublimation and deposition at the local, time-varying temperature and density; temperature-dependent ohmic dissipation; and various initial magnetic fields. The results show magnetorotational turbulence around the sublimation front at 0.5 au. The disk interior to 0.8 au is turbulent, with velocities exceeding 10% of the sound speed. Beyond 0.8 au is the dead zone, cooler than 1000 K and with turbulence orders of magnitude weaker. A local pressure maximum just inside the dead zone concentrates solid particles, favoring their growth. Over many orbits, a vortex develops at the dead zone’s inner edge, increasing the disk’s thickness locally by around 10%. We synthetically observe the results using Monte Carlo transfer calculations, finding that the sublimation front is near-infrared bright. The models with net vertical magnetic fields develop extended, magnetically supported atmospheres that reprocess extra starlight, raising the near-infrared flux 20%. The vortex throws a nonaxisymmetric shadow on the outer disk. At wavelengths >2 μm, the flux varies several percent on monthly timescales. The variations are more regular when the vortex is present. The vortex is directly visible as an arc at ultraviolet through near-infrared wavelengths, given sub-au spatial resolution.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

How to Turn Jets into Cylinders near Supermassive Black Holes in 3D General Relativistic Magnetohydrodynamic Simulations

Abstract Accreting supermassive black holes (SMBHs) produce highly magnetized relativistic jets that tend to collimate gradually as they propagate outward. However, recent radio interferometric observations of the 3C 84 galaxy reveal a stunning, cylindrical jet already at several hundred SMBH gravitational radii,r≳ 350r g . We explore how such extreme collimation emerges via a suite of 3D general relativistic magnetohydrodynamic simulations. We consider an SMBH surrounded by a magnetized torus immersed in a constant-density ambient medium that starts at the edge of the SMBH sphere of influence, chosen to be much larger than the SMBH gravitational radius,r B = 10 3 r g . We find that radiatively inefficient accretion flows (e.g., M87) produce winds that collimate the jets into parabolas near the black hole. After the disk winds stop collimating the jets atr≲r B , they turn conical. Once outsider B , the jets run into the ambient medium and form backflows that collimate the jets into cylinders some distance beyondr B . Interestingly, for radiatively efficient accretion, as in 3C 84, the radiative cooling saps the energy out of the disk winds; at early times, they cannot efficiently collimate the jets, which skip the initial parabolic collimation stage, start out conical near the SMBH, and turn into cylinders already atr≃ 300r g , as observed in 3C 84. Over time, the jet power remains approximately constant, whereas the mass accretion rate increases; the winds grow in strength and start to collimate the jets, which become quasi-parabolic near the base, and the transition point to a nearly cylindrical jet profile moves outward while remaining insider B .

Astronomy & Astrophysics↗
NASA NTRSDate not supplied

3D Phase-Field Simulations of Pattern Formation During Freeze Casting

We present the results of a combined experimental and phase-field modeling study of pattern formation during freeze casting. Those studies use unidirectional freezing of simple binary liquid mixtures of water and sugars (sucrose and trehalose) in a temperature gradient, which suffice to produce hierarchical templated structures similar to those observed in more complex multi-component freeze-cast systems including lamellae, undulated ridges, and more exotic “jellyfish-like” substructures. Multiscale 3D phase-field simulations reproduce remarkably well those structures quantitatively and identify key properties of the ice-water interface that control their formation. They further reveal that lamellae form as a result of a novel symmetry-breaking secondary instability of partially faceted cellular structures and pinpoint additional secondary instability mechanisms giving rise to smaller-scale substructures.

Kaihua Ji↗
NASA NTRSDate not supplied

3D Phase-Field Simulations of Pattern Formation During Freeze Casting

We present the results of a combined experimental and phase-field modeling study of pattern formation during freeze casting. Those studies use unidirectional freezing of simple binary liquid mixtures of water and sugars (sucrose and trehalose) in a temperature gradient, which suffice to produce hierarchical templated structures similar to those observed in more complex multi-component freeze-cast systems including lamellae, undulated ridges, and more exotic “jellyfish-like” substructures. Multiscale 3D phase-field simulations reproduce remarkably well those structures quantitatively and identify key properties of the ice-water interface that control their formation. They further reveal that lamellae form as a result of a novel symmetry-breaking secondary instability of partially faceted cellular structures and pinpoint additional secondary instability mechanisms giving rise to smaller-scale substructures.

Kaihua Ji↗
NASA NTRS2000

Statistically Reliable 'Atomistic' Simulation of Sub 100 nm MOSFETs

A 3D 'atomistic' simulation technique to study random impurity induced threshold voltage lowering and fluctuations in sub 0. 1 micron MOSFETs is presented. It allows statistical analysis of random impurity effects down to the individual impurity level-Efficient algorithms based on a single solution of Poisson's equation, followed by the solution of a simplified current continuity equation are used in the simulations.

Asenov, Asen↗
DOE OSTI2019

Earthquake Ground Motion Simulations on the Sierra High-Performance Computing System

Physics-based three-dimensional (3D) numerical simulations of earthquake ground motions provide an alternative to reliance on limited and highly scattered empirical data and models for seismic hazard assessment. We performed fully 3D earthquake ground motion simulations on the Sierra and Lassen high-performance computing (HPC) platforms at Livermore Computing. These new HPC platforms feature graphics processing unit (GPU) acceleration for improved performance. Ground motions were computed with SW4, a summation-by-parts finite difference code for anelastic wave propagation. SW4 was recently enhanced to run efficiently on GPU-accelerated platforms. We found that SW4 runs almost fifty (50) times faster in terms of node-hours on new GPU-accelerated platforms compared to conventional central processing unit (CPU) platforms with agreement of computed ground motion time-series to the full single precision which they are written. We focused on moment magnitude, M W , 7.0 earthquakes on the Hayward Fault. We were able to run three simulations resolving ground motions from static displacements (0 Hz) to frequencies of 10 Hz (wavelengths as short as 50 meters) making these Sierra calculations the highest resolution regional-scale earthquake simulations ever made for California. Results show good agreement between simulated ground motion intensities (peak ground velocity and acceleration, spectral response) and empirical models, including frequencies above 10 Hz. On Lassen we were able to run a suite of 52 M W 7.0 Hayward Fault earthquake simulations resolved 0-5 Hz for both plane-layered, one-dimensional (1D) and 3D Earth models. These simulations were used to develop a methodology for characterizing and removing systematic path and site effects which bias ground motion intensities. These simulations using the newly enhanced SW4 code allow us to run higher resolution seismic simulations with shorter run times, providing a new capability for seismic hazard and risk studies. Results were presented at scientific conferences and workshops and two publications are being prepared based on work performed in this project.

58 GEOSCIENCES↗
DOE OSTI2018

An O(N logN) multi-layer boundary element method for direct computation of sound propagation in shallow water environments

Direct three-dimensional (3D) numerical simulations of acoustic fields in range-dependent shallow water environments remains a challenge due to environmental complexities and large computational cost.We develop an efficient 3D boundary element method (BEM) for shallow water acoustic propagation which utilizes a Pre-corrected Fast Fourier Transform (PFFT) approach to reduce the computational effort from O ( N 2~3 ) to O ( N log N ) where N is the total number of boundary unknowns. To account for inhomogeneous media, the method allows for arbitrary number of coupled multi-layer BEM sub-domains. With O ( N log N ) efficiency and the use of massively parallel high-performance computing platforms, we are able to conduct multi-layer 3D direct simulations of low-mid frequency acoustics over kilometer ranges. Furthermore we perform extensive validations of the method and provide two shallow water waveguide examples benchmarked against theoretical solutions. To illustrate the efficacy and usefulness of the PFFT-BEM method, we perform 3D large-scale direct simulations to assess the performance of two established canonical models: axisymmetric coupled mode model for 3D seamount; and Kirchhoff approximation and perturbation theory for 3D rough surface scattering.

42 ENGINEERING↗
DOE OSTI2021

The impact of r -process heating on the dynamics of neutron star merger accretion disc winds and their electromagnetic radiation

Neutron star merger accretion discs can launch neutron-rich winds of >10 -2 M⊙. This ejecta is a prime site for r-process nucleosynthesis, which will produce a range of radioactive heavy nuclei. The decay of these nuclei releases enough energy to accelerate portions of the wind by ~0.1c. Here, we investigate the effect of r-process heating on the dynamical evolution of disc winds. We extract the wind from a 3D general relativistic magnetohydrodynamic simulation of a disc from a post-merger system. This is used to create inner boundary conditions for 2D hydrodynamic simulations that continue the original 3D simulation. We perform two such simulations: one that includes the r-process heating, and another one that does not. We follow the hydrodynamic simulations until the winds reach homology (60 s). Using time-dependent multifrequency multidimensional Monte Carlo radiation transport simulations, we then calculate the kilonova light curves from the winds with and without dynamical r-process heating. We find that the r-process heating can substantially alter the velocity distribution of the wind, shifting the mass-weighted median velocity from 0.06c to 0.12c. The inclusion of the dynamical r-process heating makes the light curve brighter and bluer at $\sim 1\, \mathrm{d}$ post-merger. However, the high-velocity tail of the ejecta distribution and the early ($\lesssim 1\, \mathrm{d}$) light curves are largely unaffected.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Kinematic Ground-Motion Simulations on Rough Faults Including Effects of 3D Stochastic Velocity Perturbations

Here, we describe a methodology for generating kinematic earthquake ruptures for use in 3D ground–motion simulations over the 0–5 Hz frequency band. Our approach begins by specifying a spatially random slip distribution that has a roughly wavenumber–squared fall–off. Given a hypocenter, the rupture speed is specified to average about 75%–80% of the local shear wavespeed and the prescribed slip–rate function has a Kostrov–like shape with a fault–averaged rise time that scales self–similarly with the seismic moment. Both the rupture time and rise time include significant local perturbations across the fault surface specified by spatially random fields that are partially correlated with the underlying slip distribution. We represent velocity–strengthening fault zones in the shallow (<5 km) and deep (>15 km) crust by decreasing rupture speed and increasing rise time in these regions. Additional refinements to this approach include the incorporation of geometric perturbations to the fault surface, 3D stochastic correlated perturbations to the P– and S–wave velocity structure, and a damage zone surrounding the shallow fault surface characterized by a 30% reduction in seismic velocity. We demonstrate the approach using a suite of simulations for a hypothetical Mw 6.45 strike–slip earthquake embedded in a generalized hard–rock velocity structure. The simulation results are compared with the median predictions from the 2014 Next Generation Attenuation–West2 Project ground–motion prediction equations and show very good agreement over the frequency band 0.1–5 Hz for distances out to 25 km from the fault. Additionally, the newly added features act to reduce the coherency of the radiated higher frequency (f>1 Hz) ground motions, and homogenize radiation–pattern effects in this same bandwidth, which move the simulations closer to the statistical characteristics of observed motions as illustrated by comparison with recordings from the 1979 Imperial Valley earthquake.

58 GEOSCIENCES↗
DOE OSTI2013

THE ROLE OF A FLUX ROPE EJECTION IN A THREE-DIMENSIONAL MAGNETOHYDRODYNAMIC SIMULATION OF A SOLAR FLARE

We investigated the dynamic evolution of a three-dimensional (3D) flux rope eruption and magnetic reconnection process in a solar flare by simply extending the two-dimensional (2D) resistive magnetohydrodynamic simulation model of solar flares with low β plasma to a 3D model. We succeeded in reproducing a current sheet and bi-directional reconnection outflows just below the flux rope during the eruption in our 3D simulations. We calculated four cases of a strongly twisted flux rope and a weakly twisted flux rope in 2D and 3D simulations. The time evolution of a weakly twisted flux rope in the 3D simulation shows behaviors similar to those of the 2D simulation, while a strongly twisted flux rope in the 3D simulation clearly shows a different time evolution from the 2D simulation except for the initial phase evolution. The ejection speeds of both strongly and weakly twisted flux ropes in 3D simulations are larger than in the 2D simulations, and the reconnection rates in 3D cases are also larger than in the 2D cases. This indicates positive feedback between the ejection speed of a flux rope and the reconnection rate even in the 3D simulation, and we conclude that the plasmoid-induced reconnection model can be applied to 3D. We also found that small-scale plasmoids are formed inside a current sheet and make it turbulent. These small-scale plasmoid ejections have a role in locally increasing the reconnection rate intermittently as observed in solar flares, coupled with a global eruption of a flux rope.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

Simulations using a 3D-moment turbulence closure in a two-dimensional cloud model - Comparisons with observations

A method for deep cumulus cloud simulations is presented, in which the plumes or clouds are simulated using a two-dimensional axisymmetric numerical cloud model with third-moment turbulence closure. This method is based on the assumption that the only uncertainty in the turbulence closure is in the turbulent length scale specification. Therefore, the length scale is varied, and, using different length scales, the simulated plumes are compared to the observed ones, making it possible to optimize the length scale prescription.

Krueger, Steven K.↗
NASA NTRS2014

Extending 3D Near-Cloud Corrections from Shorter to Longer Wavelengths

Satellite observations have shown a positive correlation between cloud amount and aerosol optical thickness (AOT) that can be explained by the humidification of aerosols near clouds, and/or by cloud contamination by sub-pixel size clouds and the cloud adjacency effect. The last effect may substantially increase reflected radiation in cloud-free columns, leading to overestimates in the retrieved AOT. For clear-sky areas near boundary layer clouds the main contribution to the enhancement of clear sky reflectance at shorter wavelengths comes from the radiation scattered into clear areas by clouds and then scattered to the sensor by air molecules. Because of the wavelength dependence of air molecule scattering, this process leads to a larger reflectance increase at shorter wavelengths, and can be corrected using a simple two-layer model. However, correcting only for molecular scattering skews spectral properties of the retrieved AOT. Kassianov and Ovtchinnikov proposed a technique that uses spectral reflectance ratios to retrieve AOT in the vicinity of clouds; they assumed that the cloud adjacency effect influences the spectral ratio between reflectances at two wavelengths less than it influences the reflectances themselves. This paper combines the two approaches: It assumes that the 3D correction for the shortest wavelength is known with some uncertainties, and then it estimates the 3D correction for longer wavelengths using a modified ratio method. The new approach is tested with 3D radiances simulated for 26 cumulus fields from Large-Eddy Simulations, supplemented with 40 aerosol profiles. The results showed that (i) for a variety of cumulus cloud scenes and aerosol profiles over ocean the 3D correction due to cloud adjacency effect can be extended from shorter to longer wavelengths and (ii) the 3D corrections for longer wavelengths are not very sensitive to unbiased random uncertainties in the 3D corrections at shorter wavelengths.

aerosol optical depth↗
NASA NTRS2013

Using the Flow-3D General Moving Object Model to Simulate Coupled Liquid Slosh - Container Dynamics on the SPHERES Slosh Experiment: Aboard the International Space Station

The SPHERES Slosh Experiment (SSE) is a free floating experimental platform developed for the acquisition of long duration liquid slosh data aboard the International Space Station (ISS). The data sets collected will be used to benchmark numerical models to aid in the design of rocket and spacecraft propulsion systems. Utilizing two SPHERES Satellites, the experiment will be moved through different maneuvers designed to induce liquid slosh in the experiment's internal tank. The SSE has a total of twenty-four thrusters to move the experiment. In order to design slosh generating maneuvers, a parametric study with three maneuvers types was conducted using the General Moving Object (GMO) model in Flow-30. The three types of maneuvers are a translation maneuver, a rotation maneuver and a combined rotation translation maneuver. The effectiveness of each maneuver to generate slosh is determined by the deviation of the experiment's trajectory as compared to a dry mass trajectory. To fully capture the effect of liquid re-distribution on experiment trajectory, each thruster is modeled as an independent force point in the Flow-3D simulation. This is accomplished by modifying the total number of independent forces in the GMO model from the standard five to twenty-four. Results demonstrate that the most effective slosh generating maneuvers for all motions occurs when SSE thrusters are producing the highest changes in SSE acceleration. The results also demonstrate that several centimeters of trajectory deviation between the dry and slosh cases occur during the maneuvers; while these deviations seem small, they are measureable by SSE instrumentation.

Schulman, Richard↗
DOE OSTI2021

Large-Scale Welding Process Simulation by GPU Parallelized Computing

The computational design of industrially relevant welded structures is extremely time consuming due to coupled physics and high nonlinearity. Previously, most welding distortion and residual stress simulations have been limited to small coupons and reduced order (from three-dimensional [3D] to two-dimensional [2D]), or inherent strain approximations were used for large structures. In this current study, an explicit finite element code based on a graphics processing unit was utilized to perform 3D transient thermomechanical simulation of structural components during welding. Laser brazing of aluminum alloy panels as representative of automotive manufacturing scenarios was simulated to predict out-of-plane distortion under different clamping conditions. The predicted deformation pattern and magnitude were validated by laser scanning data of physical assemblies. In addition, the code was used to investigate residual stresses developed during multipass arc welding of a nuclear industry pressurizer surge nozzle and subsequent welding repair where a 3D simulation was necessary. Taking the experimental data as reference, the 3D model predicted better residual stress distribution than a typical 2D asymmetrical model. Stress evolution in welding repair was also presented and discussed in this study. Furthermore, the efficient numerical model made it feasible to use integrated computational welding engineering to simulate welding processes for large-scale structures.

97 MATHEMATICS AND COMPUTING↗