Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic 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 181 records · Page 10

Simulation of 100-300 GHz solid-state harmonic sources

Accurate and efficient simulations of the large-signal time-dependent characteristics of second-harmonic Transferred Electron Oscillators (TEO's) and Heterostructure Barrier Varactor (HBV) frequency triplers have been obtained. This is accomplished by using a novel and efficient harmonic-balance circuit analysis technique which facilitates the integration of physics-based hydrodynamic device simulators. The integrated hydrodynamic device/harmonic-balance circuit simulators allow TEO and HBV circuits to be co-designed from both a device and a circuit point of view. Comparisons have been made with published experimental data for both TEO's and HBV's. For TEO's, excellent correlation has been obtained at 140 GHz and 188 GHz in second-harmonic operation. Excellent correlation has also been obtained for HBV frequency triplers operating near 200 GHz. For HBV's, both a lumped quasi-static equivalent circuit model and the hydrodynamic device simulator have been linked to the harmonic-balance circuit simulator. This comparison illustrates the importance of representing active devices with physics-based numerical device models rather than analytical device models.

NONLINEAR CIRCUITS↗

IQ Collaboratory. II. The Quiescent Fraction of Isolated, Low-mass Galaxies across Simulations and Observations

We compare three major large-scale hydrodynamical galaxy simulations (EAGLE, Illustris-TNG, and SIMBA) by forward modeling simulated galaxies into observational space and computing the fraction of isolated and quiescent low-mass galaxies as a function of stellar mass. Using SDSS as our observational template, we create mock surveys and synthetic spectroscopic and photometric observations of each simulation, adding realistic noise and observational limits. All three simulations show a decrease in the number of quiescent, isolated galaxies in the mass range M * = 10 9−10 M ꙩ , in broad agreement with observations. However, even after accounting for observational and selection biases, none of the simulations reproduce the observed absence of quiescent field galaxies below M * = 10 9 M ꙩ . We find that the low-mass quiescent populations selected via synthetic observations have consistent quenching timescales, despite an apparent variation in the late-time star formation histories. The effect of increased numerical resolution is not uniform across simulations and cannot fully mitigate the differences between the simulations and the observations. The framework presented here demonstrates a path toward more robust and accurate comparisons between theoretical simulations and galaxy survey observations, while the quenching threshold serves as a sensitive probe of feedback implementations.

Claire M. Dickey↗

Solar Flare Heating with Turbulent Suppression of Thermal Conduction

During solar flares plasma is typically heated to very high temperatures, and the resulting redistribution of energy via thermal conduction is a primary mechanism transporting energy throughout the flaring solar atmosphere. The thermal flux is usually modeled using Spitzer’s theory, which is based on local Coulomb collisions between the electrons carrying the thermal flux and those in the background. However, often during flares, temperature gradients become sufficiently steep that the collisional mean free path exceeds the temperature gradient scale size, so that thermal conduction becomes inherently non-local. Further, turbulent angular scattering, which is detectable in nonthermal widths of atomic emission lines, can also act to increase the collision frequency and so suppress the heat flux. Recent work by Emslie & Bian (2018) extended Spitzer’s theory of thermal conduction to account for both non-locality and turbulent suppression. We have implemented their theoretical expression for the heat flux (which is a convolution of the Spitzer flux with a kernel function) into the RADYN flare-modeling code and performed a parameter study to understand how the resulting changes in thermal conduction affect flare dynamics and hence the radiation produced. We find that models with reduced heat fluxes predict slower bulk flows, less intense line emission, and longer cooling times. By comparing features of atomic emission lines predicted by the models with Doppler velocities and nonthermal line widths deduced from a particular flare observation, we find that models with suppression factors between 0.3 to 0.5 relative to the Spitzer value best reproduce observed Doppler velocities across emission lines forming over a wide range of temperatures. Interestingly, the model that best matches observed nonthermal line widths has a kappa-type velocity distribution function

Solar flares↗

X-Ray Spectra from MHD Simulations of Accreting Black Holes

We present new global calculations of X-ray spectra from fully relativistic magneto-hydrodynamic (MHO) simulations of black hole (BH) accretion disks. With a self consistent radiative transfer code including Compton scattering and returning radiation, we can reproduce the predominant spectral features seen in decades of X-ray observations of stellar-mass BHs: a broad thermal peak around 1 keV, power-law continuum up to >100 keV, and a relativistically broadened iron fluorescent line. By varying the mass accretion rate, different spectral states naturally emerge: thermal-dominant, steep power-law, and low/hard. In addition to the spectral features, we briefly discuss applications to X-ray timing and polarization.

Schnittman, Jeremy D.↗

Global MHD Simulation of Mesoscale Structures at the Magnetospheric Boundary

We have continued our modeling of the solar wind/magnetosphere interaction for steady solar Wind conditions by carrying out a series of 3 D global magneto-hydrodynamics (MHD) simulations using a set of predetermined solar wind parameters as input to the code. With the results from the simulations that we ran last year for two other dipole tilts (0 and 35 deg), we have 36 data sets to study the dynamics of mesoscale structures at the magnetospheric boundary. We have started the analysis of these runs to investigate geometrical properties and the topology of the magnetic and electric fields for the different solar wind regimes and orientations considered. Preliminary results indicate that merging sites are consistent with patterns proposed for antiparallel merging at the dayside magnetopause. Another goal of this investigation is to establish the displacement of the cusp region as a function of the solar wind dynamic pressure, IMF direction and magnetic field dipole tilt. One of the difficulties of that study is to locate precisely the cusp in the simulation results. We have used several case studies to establish a series of criteria in to permit a routine identification of that location. We need now to process the rest of the simulation results using these criteria to obtain statistical results.

Berchem, Jean↗

Simulated Time Lags of Hinode/XRT and SDO/AIA Lightcurves as an Indication of Loop Heating Scenario

The precise nature of the heating mechanism in coronal loops is still a matter of enormous research. We present the results from a 1D hydrodynamic loop simulation of a coronal loop which was run using different parameters such as loops length (50, 200 and 500 Mm), maximum temperature reached (3MK and 10 MK0, and abundances. For each scenario the model outputs were used to calculate the corresponding lightcurves as seen by XRT/Be-thin and various EUV AIA channels. The lag time between the peak of these lightcurves was computed and tested using cross-correlation and plotted as a function of loop length.

Loop heating↗

Simulated Time Lags of Hinode/XRT and SDO/AIA Lightcurves as an Indication of Loop Heating Scenario

The precise nature of the heating mechanism (location, duration) in coronal loops is still a matter of enormous research. We present results from a 1D hydrodynamic loop simulation of a coronal loop which was run using different parameters such as loops length (50, 200, and 500 Mm), maximum temperature reached (3MK and 10MK), and abundances. For each scenario the model outputs were used to calculate the corresponding lightcurves as seen by XRT/Be-thin and various EUV AIA channels. The lag time between the peak of these lightcurves was computed using cross-correlation and plotted as a function of loop length. Additional results were computed using the 0D EBTEL code in order to test the compatibility of the two codes and to investigate additional loop lengths. Initial results indicate that the long (greater than 5000s) lags observed in the approx. 100Mm loops of active regions can only be reproduced using photospheric abundances and much longer loop lengths.

Alexander, Caroline↗

Simulated Time Lags of Hinode/XRT and SDO/AIA Lightcurves as an Indication of Loop Heating Scenario

The precise nature of the heating mechanism (location, duration) in coronal loops is still a matter of enormous research. We present results from a one-dimensional (1D) hydrodynamic loop simulation of a coronal loop which was run using different parameters such as loops length (50, 200 and 500 light-seconds), maximum temperature reached (3 million degrees Kelvin and 10 million degrees Kelvin), and abundances. For each scenario the model outputs were used to calculate the corresponding lightcurves as seen by X-ray telescope/Be-thin filter and various Extreme Ultra Violet Atmospheric Imaging Assembly channels. The lag time between the peak of these lightcurves was computed using cross-correlation and plotted as a function of loop length. Additional results were computed using the zero-dimensional Enthalpy-Based Thermal Evolution of Loops (EBTEL) code in order to test the compatibility of the two codes and to investigate additional loop lengths. Initial results indicate that the long (greater than 5000 seconds) lags observed in the approximately 100 light-seconds loops of active regions can only be reproduced using photospheric abundances and much longer loop lengths. This result suggests that the observed time lags cannot be completely explained by impulsive heating.

nanoflares↗

Axial focusing of impact energy in the Earth's interior: Proof-of-principle tests of a new hypothesis

A causal link between major impact events and global processes would probably require a significant change in the thermal state of the Earth's interior, presumably brought about by coupling of impact energy. One possible mechanism for such energy coupling from the surface to the deep interior would be through focusing due to axial symmetry. Antipodal focusing of surface and body waves from earthquakes is a well-known phenomenon which has previously been exploited by seismologists in studies of the Earth's deep interior. Antipodal focusing from impacts on the Moon, Mercury, and icy satellites has also been invoked by planetary scientists to explain unusual surface features opposite some of the large impact structures on these bodies. For example, 'disrupted' terrains have been observed antipodal to the Caloris impact basis on Mercury and Imbrium Basin on the Moon. Very recently there have been speculations that antipodal focusing of impact energy within the mantle may lead to flood basalt and hotspot activity, but there has not yet been an attempt at a rigorous model. A new hypothesis was proposed and preliminary proof-of-principle tests for the coupling of energy from major impacts to the mantle by axial focusing of seismic waves was performed. Because of the axial symmetry of the explosive source, the phases and amplitudes are dependent only on ray parameter (or takeoff angle) and are independent of azimuthal angle. For a symmetric and homogeneous Earth, all the seismic energy radiated by the impact at a given takeoff angle will be refocused (minus attenuation) on the axis of symmetry, regardless of the number of reflections and refractions it has experienced. Mantle material near the axis of symmetry will experience more strain cycles with much greater amplitude than elsewhere and will therefore experience more irreversible heating. The situation is very different than for a giant earthquake, which in addition to having less energy, has an asymmetric focal mechanism and a larger area. Two independent proof-of-principle approaches were used. The first makes use of seismic simulations, which are being performed with a realistic Earth model to determine the degree of focusing along the axis and to estimate the volume of material, if any, that experiences significant irreversible heating. The second involves two-dimensional hydrodynamic code simulations to determine the stress history, internal energy, and temperature rise as a function of radius along the axis.

Boslough, M. B.↗

Simulation of Benchmark Cases with the Terminal Area Simulation System (TASS)

The hydrodynamic core of the Terminal Area Simulation System (TASS) is evaluated against different benchmark cases. In the absence of closed form solutions for the equations governing atmospheric flows, the models are usually evaluated against idealized test cases. Over the years, various authors have suggested a suite of these idealized cases which have become standards for testing and evaluating the dynamics and thermodynamics of atmospheric flow models. In this paper, simulations of three such cases are described. In addition, the TASS model is evaluated against a test case that uses an exact solution of the Navier-Stokes equations. The TASS results are compared against previously reported simulations of these benchmark cases in the literature. It is demonstrated that the TASS model is highly accurate, stable and robust.

Ahmad, Nashat N.↗

Simulation of Benchmark Cases with the Terminal Area Simulation System (TASS)

The hydrodynamic core of the Terminal Area Simulation System (TASS) is evaluated against different benchmark cases. In the absence of closed form solutions for the equations governing atmospheric flows, the models are usually evaluated against idealized test cases. Over the years, various authors have suggested a suite of these idealized cases which have become standards for testing and evaluating the dynamics and thermodynamics of atmospheric flow models. In this paper, simulations of three such cases are described. In addition, the TASS model is evaluated against a test case that uses an exact solution of the Navier-Stokes equations. The TASS results are compared against previously reported simulations of these banchmark cases in the literature. It is demonstrated that the TASS model is highly accurate, stable and robust.

Ahmad, Nash'at↗

Hierarchical Statistical 3D ' Atomistic' Simulation of Decanano MOSFETs: Drift-Diffusion, Hydrodynamic and Quantum Mechanical Approaches

When MOSFETs are scaled to deep submicron dimensions the discreteness and randomness of the dopant charges in the channel region introduces significant fluctuations in the device characteristics. This effect, predicted 20 year ago, has been confirmed experimentally and in simulation studies. The impact of the fluctuations on the functionality, yield, and reliability of the corresponding systems shifts the paradigm of the numerical device simulation. It becomes insufficient to simulate only one device representing one macroscopical design in a continuous charge approximation. An ensemble of macroscopically identical but microscopically different devices has to be characterized by simulation of statistically significant samples. The aims of the numerical simulations shift from predicting the characteristics of a single device with continuous doping towards estimating the mean values and the standard deviations of basic design parameters such as threshold voltage, subthreshold slope, transconductance, drive current, etc. for the whole ensemble of 'atomistically' different devices in the system. It has to be pointed out that even the mean values obtained from 'atomistic' simulations are not identical to the values obtained from continuous doping simulations. In this paper we present a hierarchical approach to the 'atomistic' simulation of aggressively scaled decanano MOSFETs. A full scale 3D drift-diffusion'atomostic' simulation approach is first described and used for verification of the more economical, but also more restricted, options. To reduce the processor time and memory requirements at high drain voltage we have developed a self-consistent option based on a thin slab solution of the current continuity equation only in the channel region. This is coupled to the Poisson's equation solution in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison with the full self-consistent solution. At low drain voltage only single solution of the nonlinear Poisson equation is sufficient to extract the current with satisfactory accuracy. A pilot version of a hydrodynamic 'atomistic' simulator has been developed in order to study the effect of the nonequilibrium, non local transport in decanano MOSFETs on the random dopant induced current fluctuations. For the first time we have also applied the density gradient approach in 3D to investigate the effect of the quantum confinement on the threshold voltage fluctuations. The developed 'atomistic' simulation techniques have been applied to study various fluctuation resistant MOSFET architectures including epitaxial and delta doped devices.

Asenov, Asen↗

Active Region Simulation with EBTEL

Active regions (ARs) are areas in the sun’s upper atmosphere that are bright in the extreme ultraviolet (EUV) and X-ray spectrum. Coronal loops are one key feature of an AR. They can be described as arcs of plasma resulting from loop-shaped magnetic fields on the sun. The average temperature of these loops is over one million kelvin, which is significantly hotter than the sun’s surface. The mechanism for the heat transfer of these high temperatures remains unknown and is referred to as the coronal heating problem. One leading hypothesis regarding this mystery is that small, impulsive heating events called nanoflares are a major contributor. Here, we explore the contribution of nanoflares in the heating of AR plasma. We utilize the Enthalpy-Based Thermal Evolution of Loops (EBTEL) program to simulate NOAA 12846, as observed on July 25, 2021. We then use that simulation to analyze how the frequency of nanoflare heating events affects the EUV and X-ray observations in existing instruments, including SDO’s AIA and Hinode’s XRT. We also predict how this AR would appear in recently developed X-ray instruments, such as NASA’s MaGIXS.

Active region↗

Models of the Cartwheel ring galaxy: Spokes and starbursts

Recent observations of this famous ring galaxy, including optical and near-infrared CCD surface photometry, and VLA radio continuum and 21 cm line mapping (Higdon 1992b, in prep.), have inspired a renewed modeling effort. Toomre's (1978, in The Large-scale Structure of the Universe, eds. Longair and Einasto) series of restricted three-body simulations demonstrated how the multiple rings could be produced in a nearly head-on galaxy collision. New models with a halo-dominated potential based on the 21 cm rotation curve are able to reproduce such details as the spacing between rings, ring widths, offset of the nucleus, and several kinematical features, thus providing strong support for the collisional theory. The new observations have shown there are little or no old stars in Cartwheel; it may consist almost entirely of gas and stars produced as a result of compression in the ring wave. To model this process Smooth Particle Hydrodynamics (SPH) simulations of the Cartwheel disk have been performed. Fixed gravitational potentials were used to represent the Cartwheel and a roughly 30 percent mass collision partner. The interaction dynamics was treated as in the usual restricted three-body approximation, and the effects of local self-gravity between disk particles were calculated. We are particularly interested in testing the theory that enhanced star formation in waves is the result of gravitational instability in the compressed region (see e.g. Kennicutt 1989, ApJ 344, 685). The gas surface density in a number of simulations was initialized to a value slightly below the threshold for local gravitational instability throughout most of the disk. The first ring wave produces relatively modest compressions (a factor of order a few), triggering instability in a narrow range of wavelengths. Self-gravity in the disk is calculated over a comparable range of scales. Simulations were run with isothermal, adiabatic, and adiabatic with radiative cooling characterized by a relatively short timescale. The isothermal approximation is good except in the vicinity of the strong second (inner) ring, and several snapshots from one case are shown in the figure below. Flocculent spiral segments are present before the collision, and these are compressed into dense knots in the ring wave. These knots are likely to be sites of vigorous star formation. In the strong rarefaction behind the outer ring most of the knots are radially stretched and sheared, giving rise to spoke-like features. A few dense knots are evidently very tightly bound, because they retain their coherence and are stretched relatively little through the rarefaction. This is in accord with evidence for continuing star formation in some spokes (Marcum, Appleton and Higdon 1992). The number and spacing of spokes is a direct function of the scale of the gravitational instability in the disk. Thus, the gravitational instability theory, together with the hypothesis that massive stars are only formed in dense knots of gas, can account for most of the distinct morphology of the Cartwheel.

Struck-Marcell, Curtis↗

High-Latitude Thermosphere Neutral Density Response to Solar Wind Dynamic Pressure Enhancement

We examine the response of the thermosphere to the impact of solar wind dynamic pressure enhancements using observations and global magneto-hydrodynamics (MHD) simulations by the OpenGGCM model. Combining neutral density observations from the Challenging Minisatellite Payload (CHAMP) and the Gravity Recovery and Climate Experiment (GRACE) satellites with simultaneous Poynting flux measurements from the Defense Meteorological Satellite Program (DMSP) F16, we find that thermospheric density as well as downward Poynting flux intensified shortly after a sudden enhancement of the solar wind dynamic pressure. The intensification manifested mostly on the dayside high-latitude region with peak intensity in the vicinity of the noon and prenoon cusp. OpenGGCM modeling results show that the ionospheric Joule heating increased abruptly in response to the sudden enhancement of the dynamic pressure in the same region as the observed Poynting flux and neutral density enhancements. The modeling results show that the enhanced Joule heating coincides, both in time and location, with the appearance of a pair of high-latitude localized field-aligned currents (FACs) in the cusp region. The FACs intensified and extended azimuthally. Coincidental with the solar wind dynamic pressure enhancement, the y component of the interplanetary magnetic field (IMF) By became strongly positive and, in addition, had some large fluctuations. We explore the separate and combined effects of the dynamic pressure and IMF By perturbations, with specifically designed simulation experiments that isolate the effect of each solar wind parameter. We find that the dynamic pressure enhancement is the primary source for the Joule heating and neutral density enhancements, but the IMF By modulates the level of enhancement.

magneto-hydrodynamics↗

Diagnostics of electron-heated solar flare models. III - Effects of tapered loop geometry and preheating

A series of hydrodynamic numerical simulations of nonthermal electron-heated solar flare atmospheres and their corresponding soft X-ray Ca XIX emission-line profiles, under the conditions of tapered flare loop geometry and/or a preheated atmosphere, is presented. The degree of tapering is parameterized by the magnetic mirror ratio, while the preheated atmosphere is parameterized by the initial upper chromospheric pressure. In a tapered flare loop, it is found that the upward motion of evaporated material is faster compared with the case where the flare loop is uniform. This is due to the diverging nozzle seen by the upflowing material. In the case where the flare atmosphere is preheated and the flare geometry is uniform, the response of the atmosphere to the electron collisional heating is slow. The upward velocity of the hydrodynamic gas is reduced due not only to the large coronal column depth, but also to the increased inertia of the overlying material. It is concluded that the only possible electron-heated scenario in which the predicted Ca XIX line profiles agree with the BCS observations is when the impulsive flare starts in a preheated dense corona.

Emslie, A. G.↗

Flows in Enthalpy-based Thermal Evolution of Loops

Plasma-filled loop structures are common in the solar corona. Because detailed modeling of the dynamical evolution of these structures is computationally costly, an efficient method for computing approximate but quick physics-based solutions is to rely on space-integrated 0D simulations. The enthalpy-based thermal evolution of loops (EBTEL) framework is a commonly used method to study the exchange of mass and energy between the corona and transition region. EBTEL solves for density, temperature, and pressure, averaged over the coronal part of the loop, velocity at coronal base, and the instantaneous differential emission measure distribution in the transition region. The current single-fluid version of the code, EBTEL2, assumes that at all stages the flows are subsonic. However, sometimes the solutions show the presence of supersonic flows during the impulsive phase of heat input. It is thus necessary to account for this effect. Here, we upgrade EBTEL2 to EBTEL3 by including the kinetic energy term in the Navier–Stokes equation. We compare the solutions from EBTEL3 with those obtained using EBTEL2, as well as the state-of-the-art field-aligned hydrodynamics code HYDRAD. We find that the match in pressure between EBTEL3 and HYDRAD is better than that between EBTEL2 and HYDRAD. Additionally, the velocities predicted by EBTEL3 are in close agreement with those obtained with HYDRAD when the flows are subsonic. However, EBTEL3 solutions deviate substantially from HYDRAD’s when the latter predicts supersonic flows. Using the mismatches in the solution, we propose a criterion to determine the conditions under which EBTEL can be used to study flows in the system.

Solar coronal heating↗

Origin of Mars's Moons By Disruptive Partial Capture of an Asteroid

The origin of Mars’s small moons, Phobos and Deimos, remains unknown. They are typically thought either to be captured asteroids or to have accreted from a debris disk produced by a giant impact. Here, we present an alternative scenario wherein fragments of a tidally disrupted asteroid are captured and evolve into a collisional proto-satellite disk. We simulate the initial disruption and the fragments’ subsequent orbital evolution. We find that tens of percent of an unbound asteroid’s mass can be captured and survive beyond collisional timescales, across a broad range of periapsis distances, speeds, masses, spins, and orientations in the Sun–Mars frame. Furthermore, more than one percent of the asteroid’s mass could evolve to circularise in the moons’ accretion region. This implies a lower mass requirement for the parent body than that for a giant impact, which could increase the likelihood of this route to forming a proto-satellite disk that, unlike direct capture, could also naturally explain the moons’ orbits. These three formation scenarios each imply different properties of Mars’s moons to be tested by upcoming spacecraft missions.

Impact phenomena↗