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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS2017

Realistic Modeling of Interaction of Quiet-Sun Magnetic Fields with the Chromosphere

High‐resolution observations and 3D MHD simulations reveal intense interaction between the convection zone dynamics and the solar atmosphere on subarcsecond scales. To investigate processes of the dynamical coupling and energy exchange between the subsurface layers and the chromosphere we perform 3D radiative MHD modeling for a computational domain that includes the upper convection zone and the chromosphere, and investigate the structure and dynamics for different intensity of the photospheric magnetic flux. For comparison with observations, the simulation models have been used to calculate synthetic Stokes profiles of various spectral lines. The results show intense energy exchange through small‐scale magnetized vortex tubes rooted below the photosphere, which provide extra heating of the chromosphere, initiate shock waves, and small‐scale eruptions.

quiet-sun↗

DOE OSTI2019

Adaptive inverse ray-tracing for accurate and efficient modeling of cross beam energy transfer in hydrodynamics simulations

Integrated hydrodynamics simulations of inertial confinement fusion rely on reduced physics models. To reproduce experimental trends, these models often feature tuning parameters, but this comes with a risk: the over-tuning of one model can hide physics inadequacies in another. The ray-based models of cross-beam-energy transfer (CBET) represent this risk. In this paper, we present an accurate and efficient model of CBET suitable for inline implementation in 3D hydrodynamics simulations. Inverse Ray Tracing (IRT) is used to compute the ray field in a 3D permittivity profile described on an unstructured tetrahedral mesh using the Inline Field Reconstruction and Interaction using Inverse Tracing framework. CBET is accounted for through perturbations to the permittivity associated with ion acoustic waves driven by the overlapped fields. Large gradients in the permittivity are resolved by coupling the IRT to a recursive Adaptive Mesh Refinement (AMR) algorithm. The use of AMR also allows for the resolution of caustics, with accurate field reconstruction performed using the Etalon integral method. Comparisons of the model with wave-based solutions from the Laser Plasma Simulation Environment demonstrate its ability to control energy conservation and gain convergence through the AMR depth only, without the use of ad hoc physical models or artificial tuning parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

WOBBLING AND PRECESSING JETS FROM WARPED DISKS IN BINARY SYSTEMS

We present results of the first ever three-dimensional (3D) magnetohydrodynamic (MHD) simulations of the accretion–ejection structure. We investigate the 3D evolution of jets launched symmetrically from single stars but also jets from warped disks in binary systems. We have applied various model setups and tested them by simulating a stable and bipolar symmetric 3D structure from a single star–disk–jet system. Our reference simulation maintains a good axial symmetry and also a bipolar symmetry for more than 500 rotations of the inner disk, confirming the quality of our model setup. We have then implemented a 3D gravitational potential (Roche potential) due by a companion star and run a variety of simulations with different binary separations and mass ratios. These simulations show typical 3D deviations from axial symmetry, such as jet bending outside the Roche lobe or spiral arms forming in the accretion disk. In order to find indications of precession effects, we have also run an exemplary parameter setup, essentially governed by a small binary separation of only ≃200 inner disk radii. This simulation shows a strong indication that we observe the onset of a jet precession caused by the wobbling of the jet-launching disk. We estimate the opening angle of the precession cone defined by the lateral motion of the jet axis to be about 4° after about 5000 dynamical time steps.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1992

Simulation of automatic rotorcraft nap-of-the-earth flight in graphics workstation environment

This paper describes a three-dimensional (3D) helicopter flight simulation system. The simulation is designed to be a readily available tool for concept verification and tuning of automatic obstacle-avoidance guidance algorithms. The system has been implemented on networked workstations capable of interactive 3D graphics simulation. The simulation uses realistic terrain and obstacle models. The dynamics of the rotorcraft and the functional capabilities of the range sensors are simulated to provide all the components required to evaluate the guidance function. Standard graphics hardware available on the workstation is utilized to accelerate the range-data calculations for sensor simulation at the guidance rate. An example is given to demonstrate the performance of the obstacle-avoidance capability.

Lam, T.↗

DOE OSTI2024

Development and Refinement of Bison Models for Higher Burnup Application

Electric Power Research Institute (EPRI) is an international research organization funded by organizations in the energy and government sectors. The company conducts research to tackle the energy challenges of its member organizations. In support of the U.S. nuclear industry, EPRI and Idaho National Laboratory were awarded a U.S. Department of Energy Technology Commercialization Fund to develop a loss of coolant accident (LOCA) analysis tool. This tool was developed to support burnup extension activities by supplementing the existing capabilities of operators of existing large light-water reactor plants. The tool couples INL developed codes Bison and RELAP5-3D for advanced simulation of LOCA behavior. The project consists of seven major tasks, one of which is the development and refinement of Bison models for higher burnup applications. This task is broken down into six targeted areas of development and refinement including: formation of high burnup structure, fuel fracture and pulverization, transient fission gas release, axial gas communication within a fuel rod, cladding ballooning and burst, and axial fuel relocation. This report summarizes the models for both the fuel and cladding in Bison for high burnup application, identifies the areas of development and refinement necessary to enhance the Bison capabilities to be used in the coupled simulation tool with RELAP5-3D. Integration and regression testing to demonstrate proper implementation of the new features in Bison is highlighted.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2019

Simulated Phase-dependent Spectra of Terrestrial Aquaplanets in M Dwarf Systems

Orbital phase-dependent variations in thermal emission and reflected stellar energy spectra can provide meaningful constraints on the climate states of terrestrial extrasolar planets orbiting M dwarf stars. Spatial distributions of water vapor, clouds, and surface ice are controlled by climate. In turn, water, in each of its thermodynamic phases, imposes significant modulations to thermal and reflected planetary spectra. Here we explore these characteristic spectral signals, based on 3D climate simulations of Earth-sized aquaplanets orbiting M dwarf stars near the habitable zone. By using 3D models, we can self-consistently predict surface temperatures and the location ofwater vapor, clouds, and surface ice in the climate system. Habitable zone planets in M dwarf systems are expected to be in synchronous rotation with their host star and thus present distinct differences in emitted and reflected energy fluxes depending on the observed hemisphere. Here we illustrate that icy, temperate, and incipient runaway greenhouse climate states exhibit phase-dependent spectral signals that enable their characterization.

Wolf, E. T.↗

NASA NTRS1991

Numerical simulation of swept-wing flows

Efforts of the last six months to computationally model the transition process characteristics of flow over swept wings are described. Specifically, the crossflow instability and crossflow/Tollmien-Schlichting wave interactions are analyzed through the numerical solution of the full 3D Navier-Stokes equations including unsteadiness, curvature, and sweep. This approach is chosen because of the complexity of the problem and because it appears that linear stability theory is insufficient to explain the discrepancies between different experiments and between theory and experiment. The leading edge region of a swept wing is considered in a 3D spatial simulation with random disturbances as the initial conditions.

Reed, Helen L.↗

DOE OSTI2017

Magnetohydrodynamic Simulations of the Formation of Molecular Clouds toward the Stellar Cluster Westerlund 2: Interaction of a Jet with a Clumpy Interstellar Medium

The formation mechanism of CO clouds observed with the NANTEN2 and Mopra telescopes toward the stellar cluster Westerlund 2 is studied by 3D magnetohydrodynamic simulations, taking into account the interstellar cooling. These molecular clouds show a peculiar shape composed of an arc-shaped cloud on one side of the TeV γ -ray source HESS J1023-575 and a linear distribution of clouds (jet clouds) on the other side. We propose that these clouds are formed by the interaction of a jet with clumps of interstellar neutral hydrogen (H i). By studying the dependence of the shape of dense cold clouds formed by shock compression and cooling on the filling factor of H i clumps, we found that the density distribution of H i clumps determines the shape of molecular clouds formed by the jet–cloud interaction: arc clouds are formed when the filling factor is large. On the other hand, when the filling factor is small, molecular clouds align with the jet. The jet propagates faster in models with small filling factors.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

MODELING THE ROSSITER–MCLAUGHLIN EFFECT: IMPACT OF THE CONVECTIVE CENTER-TO-LIMB VARIATIONS IN THE STELLAR PHOTOSPHERE

Observations of the Rossiter–McLaughlin (RM) effect provide information on star–planet alignments, which can inform planetary migration and evolution theories. Here, we go beyond the classical RM modeling and explore the impact of a convective blueshift that varies across the stellar disk and non-Gaussian stellar photospheric profiles. We simulated an aligned hot Jupiter with a four-day orbit about a Sun-like star and injected center-to-limb velocity (and profile shape) variations based on radiative 3D magnetohydrodynamic simulations of solar surface convection. The residuals between our modeling and classical RM modeling were dependent on the intrinsic profile width and v sin i; the amplitude of the residuals increased with increasing v sin i and with decreasing intrinsic profile width. For slowly rotating stars the center-to-limb convective variation dominated the residuals (with amplitudes of 10 s of cm s{sup −1} to ∼1 m s{sup −1}); however, for faster rotating stars the dominant residual signature was due a non-Gaussian intrinsic profile (with amplitudes from 0.5 to 9 m s{sup −1}). When the impact factor was 0, neglecting to account for the convective center-to-limb variation led to an uncertainty in the obliquity of ∼10°–20°, even though the true v sin i was known. Additionally, neglecting to properly model an asymmetric intrinsic profile had a greater impact for more rapidly rotating stars (e.g., v sin i = 6 km s{sup −1}) and caused systematic errors on the order of ∼20° in the measured obliquities. Hence, neglecting the impact of stellar surface convection may bias star–planet alignment measurements and consequently theories on planetary migration and evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Neutrino-driven Winds in Three-dimensional Core-collapse Supernova Simulations

Abstract In this paper, we analyze the neutrino-driven winds that emerge in 12 unprecedentedly long-duration 3D core-collapse supernova simulations done using the code F ornax . The 12 models cover progenitors with zero-age main-sequence mass between 9 and 60 solar masses. In all our models, we see transonic outflows that are at least 2 times as fast as the surrounding ejecta and that originate generically from a proto−neutron star surface atmosphere that is turbulent and rotating. We find that winds are common features of 3D simulations, even if there is anisotropic early infall. We find that the basic dynamical properties of 3D winds behave qualitatively similarly to those inferred in the past using simpler 1D models, but that the shape of the emergent wind can be deformed, very aspherical, and channeled by its environment. The thermal properties of winds for less massive progenitors very approximately recapitulate the 1D stationary solutions, while for more massive progenitors they deviate significantly owing to aspherical accretion. The Y e temporal evolution in winds is stochastic, and there can be some neutron-rich phases. Though no strong r -process is seen in any model, a weak r -process can be produced, and isotopes up to 90 Zr are synthesized in some models. Finally, we find that there is at most a few percent of a solar mass in the integrated wind component, while the energy carried by the wind itself can be as much as 10%–20% of the total explosion energy.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Nonlinear excitation of the ablative Rayleigh-Taylor instability for all wave numbers

Small-scale perturbations in the ablative Rayleigh-Taylor instability (ARTI) are often neglected because they are linearly stable when their wavelength is shorter than a linear cutoff. Using 2D and 3D numerical simulations, it is shown that linearly stable modes of any wavelength can be destabilized. This instability regime requires finite amplitude initial perturbations and linearly stable ARTI modes are more easily destabilized in 3D than in 2D. In conclusion, it is shown that for conditions found in laser fusion targets, short wavelength ARTI modes are more efficient at driving mixing of ablated material throughout the target since the nonlinear bubble density increases with the wave number and small scale bubbles carry a larger mass flux of mixed material.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2014

Radiation asymmetries during disruptions on DIII-D caused by massive gas injection

One of the major challenges that the ITER tokamak will have to face during its operations are disruptions. During the last few years, it has been proven that the global consequences of a disruption can be mitigated by the injection of large quantities of impurities. But one aspect that has been difficult to study was the possibility of local effects inside the torus during such injection that could damage a portion of the device despite the global heat losses and generated currents remaining below design parameter. 3D MHD simulations show that there is a potential for large toroidal asymmetries of the radiated power during impurity injection due to the interaction between the particle injection plume and a large n = 1 mode. Another aspect of 3D effects is the potential occurrence of Vertical Displacement Events (VDE), which could induce large poloidal heat load asymmetries. This potential deleterious effect of 3D phenomena has been studied on the DIII-D tokamak, thanks to the implementation of a multi-location massive gas injection (MGI) system as well as new diagnostic capabilities. This study showed the existence of a correlation between the location of the n = 1 mode and the local heat load on the plasma facing components but shows also that this effect is much smaller than anticipated (peaking factor of ∼1.1 vs 3-4 according to the simulations). There seems to be no observable heat load on the first wall of DIII-D at the location of the impurity injection port as well as no significant radiation asymmetries whether one or 2 valves are fired. This study enabled the first attempt of mitigation of a VDE using impurity injection at different poloidal locations. The results showed a more favorable heat deposition when the VDE is mitigated early (right at the onset) by impurity injection. No significant improvement of the heat load mitigation efficiency has been observed for late particle injection whether the injection is done “in the way” of the VDE (upward VDE mitigated by injection from the upper part of the vessel vs the lower part) or not.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Integrating static PRA information with risk informed safety margin characterization (RISMC) simulation methods

The overall objective of the project was to develop a computationally feasible and user-friendly mechanized process to integrate traditional probabilistic risk assessment (PRA) and dynamic PRA (DPRA) results. Starting with the systematic identification of items in an existing PRA that need dynamic augmentation, the project used a generic 4-loop pressurized reactor (PWR) and 3-loop PWR as example plants. Station blackout (SBO) and large break loss of coolant accident SBLOCA) were selected as the example initiating events. Using the traditional event-tree (ET)/fault-tree (FT) methodology augmented by dynamic evet tree approach, the potential consequences of the initiating events were simulated with RELAP-3D and MELCOR/RASCAL codes to cover Level 1 through Level 3 of PRA. RAVEN and ADAPT software were used to generate Level 1 simulations with RELAP-3D and Level 2/3 simulations with MELCOR (Level 2)/RASCAL (Level 3), respectively. Example branching conditions (BCs) for SBO included AC power recovery time, valve repair failure time, reactor coolant pump leak time/break size and emergency power supply duration to a total of 9. Example BCs for LOCA included off-site power recovery time, diesel generator power recovery time, auxiliary feed water system operation time, safety relief valve failure to open upon demand, reactor coolant pump seal break time and size to a total of 21. Each RELAP-3D simulation (9,587 scenarios) was labelled OK or Core Damage based on the maximum allowed peak clad temperature (2,100oF). Each MELCOR simulation (4610 scenarios) was labeled as Bin over 10rem or Bin 0-10rem based on the dose at the site boundary. The scenarios were clustered based on the criteria above using the mean shift methodology. Classical PRA (CPRA) and DPRA results were compared to identify the ET sequences that need DPRA augmentation. Several approaches were proposed for the incorporation of these sequences into CPRA using clustering with the mean shift methodology, restructuring the CPRA ETs by adding new BCs/sequences, and using the concept of a limit surface. Procedures for decision making regarding the possible consequences of an initiating event (e.g. core damage or not, site evacuation or not) were developed using a convolutional neural network (CNN), a recurrent neural network (RNN) and a transformer neural network (TNN). The project has led to two PhD degrees, three archival journal papers and five refereed conference proceedings.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2017

Nonthermal Particle Acceleration in 3D Relativistic Magnetic Reconnection in Pair Plasma

As a fundamental process converting magnetic to plasma energy in high-energy astrophysical plasmas, relativistic magnetic reconnection is a leading explanation for the acceleration of particles to the ultrarelativistic energies that are necessary to power nonthermal emission (especially X-rays and gamma-rays) in pulsar magnetospheres and pulsar wind nebulae, coronae and jets of accreting black holes, and gamma-ray bursts. An important objective of plasma astrophysics is therefore the characterization of nonthermal particle acceleration (NTPA) effected by reconnection. Reconnection-powered NTPA has been demonstrated over a wide range of physical conditions using large 2D kinetic simulations. However, its robustness in realistic 3D reconnection—in particular, whether the 3D relativistic drift-kink instability (RDKI) disrupts NTPA—has not been systematically investigated, although pioneering 3D simulations have observed NTPA in isolated cases. Here, we present the first comprehensive study of NTPA in 3D relativistic reconnection in collisionless electron–positron plasmas, characterizing NTPA as the strength of 3D effects is varied systematically via the length in the third dimension and the strength of the guide magnetic field. We find that, while the RDKI prominently perturbs 3D reconnecting current sheets, it does not suppress particle acceleration, even for zero guide field; fully 3D reconnection robustly and efficiently produces nonthermal power-law particle spectra closely resembling those obtained in 2D. This finding provides strong support for reconnection as the key mechanism powering high-energy flares in various astrophysical systems. We also show that strong guide fields significantly inhibit NTPA, slowing reconnection and limiting the energy available for plasma energization, yielding steeper and shorter power-law spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2018

Sandia's Current Energy Conversion module for Flexible-Mesh Delft3D flow solver v. 1.0

The DOE has funded Sandia National Labs (SNL) to develop an open-source modeling tool to guide the design and layout of marine hydrokinetic (MHK) arrays to maximize power production while minimizing environmental effects. This modeling framework simulates flows through and around a MHK arrays while quantifying environmental responses. As an augmented version of the Dutch company, Deltares’s, environmental hydrodynamics code, Delft3D, SNL-Delft3D-CEC-FM includes a new module that simulates energy conversion (momentum withdrawal) by MHK current energy conversion devices with commensurate changes in the turbulent kinetic energy and its dissipation rate. SNL-Delft3D-CEC-FM modified the Delft3D flexible mesh flow solver, DFlowFM. The Following is a description of Deltares’s open-source code Delft3D from which Delft3D-CEC is built upon. “Delft3D is a world leading 3D modeling suite to investigate hydrodynamics, sediment transport and morphology and water quality for fluvial, estuarine and coastal environments. As per 1 January 2011, the Delft3D flow (FLOW), morphology (MOR) and waves (WAVE) modules are available in open source. The software is used and has proven his capabilities on many places around the world, like the Netherlands, USA, Hong Kong, Singapore, Australia, Venice, etc. The software is continuously improved and developed with innovating advanced modelling techniques as consequence of the research work of our institute and to stay world leading. The FLOW module is the heart of Delft3D and is a multi-dimensional (2D or 3D) hydrodynamic (and transport) simulation program which calculates non-steady flow and transport phenomena resulting from tidal and meteorological forcing on a curvilinear, boundary fitted grid or spherical coordinates. In 3D simulations, the vertical grid is defined following the so-called sigma coordinate approach or Z-layer approach. The MOR module computes sediment transport (both suspended and bed total load) and morphological changes for an arbitrary number of cohesive and non-cohesive fractions. Both currents and waves act as driving forces and a wide variety of transport formulae have been incorporated. For the suspended load this module connects to the 2D or 3D advection-diffusion solver of the FLOW module; density effects may be taken into account. An essential feature of the MOR module is the dynamic feedback with the FLOW and WAVE modules, which allow the flows and waves to adjust themselves to the local bathymetry and allows for simulations on any time scale from days (storm impact) to centuries (system dynamics). It can keep track of the bed composition to build up a stratigraphic record. The MOR module may be extended to include extensive features to simulate dredging and dumping scenarios. For over 30 years Deltares has been in the forefront of these types of combined morphological simulation techniques.”

Chartand, Chris↗

DOE OSTI2022

Kinetic Simulations of Instabilities and Particle Acceleration in Cylindrical Magnetized Relativistic Jets

Relativistic magnetized jets, such as those from AGN, GRBs, and XRBs, are susceptible to current- and pressure-driven MHD instabilities that can lead to particle acceleration and nonthermal radiation. Here, we investigate the development of these instabilities through 3D kinetic simulations of cylindrically symmetric equilibria involving toroidal magnetic fields with electron–positron pair plasma. Generalizing recent treatments by Alves et al. and Davelaar et al., we consider a range of initial structures in which the force due to toroidal magnetic field is balanced by a combination of forces due to axial magnetic field and gas pressure. We argue that the particle energy limit identified by Alves et al. is due to the finite duration of the fast magnetic dissipation phase. We find a rather minor role of electric fields parallel to the local magnetic fields in particle acceleration. In all investigated cases, a kink mode arises in the central core region with a growth timescale consistent with the predictions of linearized MHD models. In the case of a gas-pressure-balanced (Z-pinch) profile, we identify a weak local pinch mode well outside the jet core. We argue that pressure-driven modes are important for relativistic jets, in regions where sufficient gas pressure is produced by other dissipation mechanisms.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Towards Exascale Astrophysics of Mergers and Supernovae (TEAMS)

The TEAMS project brought together cutting-edge simulations, theoretical insights, and collaborative efforts to deepen our understanding of some of the universe’s most extreme phenomena—supernovae, neutron star mergers, and the powerful signals they emit. Using one of the largest suites of 3D supernova simulations ever conducted, researchers uncovered new insights into how massive stars explode, how those explosions vary by stellar mass, and what conditions lead to the birth of neutron stars or black holes. They also studied the radiation and gravitational wave signals emitted during these events, revealing how future observations can be used to uncover what happens deep inside collapsing stars. The team developed improved tools for modeling how light and neutrinos behave in such explosive environments, enabling more accurate predictions of what astronomers might observe. Work also explored how the chemical composition and geometry of kilonovae—the visible explosions that follow neutron star mergers—influence their signals and can reveal the origins of heavy elements like gold. These efforts not only advanced scientific knowledge, but also trained a new generation of researchers at the intersection of astrophysics, computational science, and nuclear theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2019

Independent Effects of Reynolds Number and Mach Number on the Aerodynamics of an Iced Swept Wing

Aerodynamic assessment of icing effects on swept wings is an important component of a larger effort to improve three-dimensional (3D) icing simulation capabilities. An understanding of ice-shape geometric fidelity and Reynolds and Mach number effects on the iced-wing aerodynamics is needed to guide the development and validation of ice-accretion simulation tools. To this end, wind tunnel testing was carried out for a 13.3-percent-scale semispan wing based upon the Common Research Model airplane configuration. The wind tunnel testing was conducted at the Office National d’Etudes et de Recherches Aérospatiales (ONERA) F1 pressurized wind tunnel with Reynolds numbers of 1.6 x 10(exp 6) to 11.9 x 10(exp 6 ) and Mach numbers of 0.09 to 0.34. Five different configurations were investigated using fully 3D, high-fidelity artificial ice shapes that maintain nearly all of the 3D ice-accretion features documented in prior icing wind tunnel tests. These large, leading-edge ice shapes were nominally based upon airplane holding in icing conditions scenarios. For three of these configurations, lower fidelity simulations were also built and tested. The results presented in this paper show that while Reynolds and Mach number effects are important for quantifying the clean-wing performance, there is very little to no effect for an iced wing with 3D, high-fidelity artificial ice shapes or 3D smooth ice shapes with grit roughness. These conclusions are consistent with the large volume of past research on iced airfoils. However, some differences were also noted for the associated stalling angle of the iced swept wing and for various lower fidelity versions of the leading-edge ice accretion. More research is planned to further investigate the key features of ice-accretion geometry that must be simulated in lower fidelity versions in order to capture the essential aerodynamics.

Aircraft icing↗