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.

SourceResearch connectionFrom yearThrough year

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

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Contribution to Solar Brightness of Small-Size Magnetic Elements

Irradiance variability is mostly driven by surface magnetism, each magnetic feature contributing in a different manner. The contribution of small-size magnetic elements observed ubiquitously at high-spatial resolution on the solar photosphere, is still debated, as such features are mostly unresolved on full-disk images employed to model irradiance variability. Understanding the contribution to solar brightness of small-size magnetic elements, especially in quiet regions, will help understanding irradiance variability on the decadal and longer temporal scales, which, in turn, are fundamental to understand the role of the Sun on the Earth’s climate. We present a preliminary study of the brightness of small-size magnetic elements using high-spatial resolution spectro-polarimetric observations acquired with the National Science Foundations Daniel K. Inouye Solar Telescope (DKIST) during Commissioning Phase 1. We focus on the Intensity contrast vs magnetic field relation, which is a fundamental observable in irradiance studies. Previous studies mostly focused on the FeI 630.1/630.2 nm spectral range and so to contextualize our results we focus here on this spectral region as well. By comparing DKIST results with MURaM magneto-hydrodynamic simulations, and with previous results published in the literature, we conclude that our results are in agreement with studies conducted on HINODE observations. We plan to extend the analysis to UV and IR ViSP observations and VBI blue continuum in order to better constrain the quiet sun contribution to spectral irradiance variations.

SMD↗

NASA NTRS2023

Prospects of Detecting Nonthermal Protons in Solar Flares via Lyman Line Spectroscopy: Revisiting the Orrall-Zirker Effect

Solar flares are efficient particle accelerators, with a substantial fraction of the energy released manifesting as nonthermal particles. While the role that nonthermal electrons play in transporting flare energy is well studied, the properties and importance of nonthermal protons are rather less well understood. This is in large part due to the paucity of diagnostics, particularly at the lower-energy (deka-keV) range of nonthermal proton distributions in flares. One means to identify the presence of deka-keV protons is by an effect originally described by Orrall & Zirker. In the Orrall–Zirker effect, nonthermal protons interact with ambient neutral hydrogen, and via charge exchange produce a population of energetic neutral atoms (ENAs) in the chromosphere. These ENAs subsequently produce an extremely redshifted photon in the red wings of hydrogen spectral lines. We revisit predictions of the strength of this effect using modern interaction cross sections, and numerical models capable of self-consistently simulating the flaring nonequilibrium ionization stratification, and the nonthermal proton distribution (and, crucially, their feedback on each other). We synthesize both the thermal and nonthermal emission from Ly α and Ly β, the most promising lines that may exhibit a detectable signal. These new predictions are weaker and more transient than prior estimates, but the effects should be detectable in fortuitous circumstances. We degrade the Ly β emission to the resolution of the Spectral Imaging of the Coronal Environment (SPICE) instrument on board Solar Orbiter, demonstrating that though likely difficult, it should be possible to detect the presence of nonthermal protons in flares observed by SPICE.

Solar flares↗

NASA NTRSDate not supplied

Comparison of Ring Diagrams Based on the Doppler Shifts of Synthetic Data Obtained with Bisector Method and SDO/HMI Pipeline

Ring diagrams are cross-sections of three-dimensional spatiotemporal power spectra of solar oscillations. The rings reveal information about sub-surface flows and represent an important tool for helioseismology. Ring diagrams can be constructed using Doppler velocity or intensity maps of spectral lines. How the velocities are computed is an important factor for accuracy of information we can retrieve on subsurface flows. In our work, the ring diagrams are generated from Doppler shift data of synthesized Fe I 6173 Å line. We compare ring diagrams computed by two methods–HMI line-of-sight pipeline and the bisector of Fe line. Fe I line is synthesized for StellarBox 3D Radiative hydrodynamic simulations under LTE assumption. We aim to answer the following questions: 1.How do power spectra obtained from velocities computed with the HMI pipeline and bisector of the Fe I 6173Å compare? 2.How do the power spectral density retrieved with each method vary with heliocentric angle? 3.What is the effect of changing resolution on the power spectral density in ring diagrams obtained with the two methods?

SMD↗

DOE OSTI2026

Dynamical Dark Energy Imprints in the Lyman-Alpha Forest

The nature of dark energy (DE) remains elusive, even though it constitutes the dominant energy-density component of the Universe and drives the late-time acceleration of cosmic expansion. By combining measurements of the expansion history from baryon acoustic oscillations, supernova surveys, and cosmic microwave background data, the Dark Energy Spectroscopic Instrument (DESI) Collaboration has inferred that the DE equation of state may evolve over time. The profound implications of a time-variable, ``dynamical" DE (DDE) that departs from a cosmological constant motivate the need for independent observational tests. In this work, we use cosmological hydrodynamical simulations of structure formation to investigate how DDE affects the properties of the Lyman-Alpha ``forest'' of absorption features produced by neutral hydrogen in the cosmic web. We find that DDE models consistent with the DESI constraints induce a spectral tilt in the forest transmitted flux power spectrum, imprinting a scale- and redshift-dependent signature relative to standard Lambda-CDM cosmologies. These models also yield higher intergalactic medium temperatures and reduced Lyman-Alpha opacity compared to Lambda-CDM. We discuss the observational implications of these trends as potential avenues for independent confirmation of DDE.

Garza, Diego [UC, Santa Cruz, Inst. Part. Phys.] (↗

DOE OSTI2026

The Challenge of Observing Patchy Reionization with CMB Optical-Depth Fluctuations

Spatial fluctuations in the Thomson optical depth encode information about the inhomogeneous nature of cosmic reionization. We compute the optical-depth angular power spectrum, $C_\ell^{ττ}$, using past lightcones constructed from five Cosmic Reionization on Computers (CROC) radiation-hydrodynamical simulations. By decomposing the electron-density field into patchy and density components, we quantify the separate contributions of ionization-fraction and baryon-density fluctuations to the optical-depth anisotropy. Because the simulations end at $z\approx5$, we supplement the reionization-era signal with an analytic estimate of the fully ionized low-redshift contribution. We find that baryon-density fluctuations dominate the high-redshift signal over most angular scales, while the accumulated low-redshift contribution exceeds the high-redshift signal across the full multipole range considered. Our results demonstrate that a significant fraction of the optical-depth power is not uniquely associated with reionization morphology, implying that future interpretations of $C_\ell^{ττ}$ must account for the density contribution in addition to patchy ionization.

Takoudes, Nick [Chicago U., Astron. Astrophys. Ctr↗

NASA NTRS2012

Improved Finite-Volume Method for Radiative Hydrodynamics

Fully coupled simulations of hydrodynamics and radiative transfer are essential to a number of fields ranging from astrophysics to engineering applications. Of particular interest in this work are hypersonic atmospheric entries and associated experimental apparatus, e.g., shock tubes and high enthalpy testing facilities. The radiative transfer calculations must supply to the CFD a heating term in the energy equation in the form of the divergence of the radiative heat flux and the radiative heat fluxes to bounding surfaces. It is most efficient to solve the radiative transfer equation on the same grid as the CFD solution, and this work presents an algorithm with improved accuracy for such simulations on structured and unstructured grids compared to more conventional approaches. Results will be shown for shock radiation during hypersonic reentry. Issues of parallelization within a radiation sweep will also be discussed.

Wray, Alan↗

NASA NTRS2015

Launch Environment Water Flow Simulations Using Smoothed Particle Hydrodynamics

This paper describes the use of Smoothed Particle Hydrodynamics (SPH) to simulate the water flow from the rainbird nozzle system used in the sound suppression system during pad abort and nominal launch. The simulations help determine if water from rainbird nozzles will impinge on the rocket nozzles and other sensitive ground support elements.

Smoothed Particle Hyrodynamics↗

DOE OSTI2025

Simulations of Radiation Hydrodynamics in Strongly Coupled Plasmas (Final Report)

Radiation propagates throughout the Universe and is often the only information we have from distant astronomical objects. Additionally, radiation plays a key role in systems ranging from the interior of stars to the establishment of shocks at the edges of forming galaxies to neutron star mergers. Here radiative processes, absorption, emission, and scattering, are fundamental to transport of momentum and energy, as well as, the overall evolution of these and other astrophysical systems. In some cases, radiation-dominant systems are also sufficiently dense and ionized, such that, the ions are strongly coupled (ie electrostatic energy >> thermal energy or Γ >> 1), which will affect radiative processes and the transport of radiation. During this project, we developed a computational ray-tracing algorithm to study radiation transport through a simulated neutron star capsule implosion under conditions relevant to Omega-60 laser experiments.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1995

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↗

DOE OSTI2025

Spherical and sessile droplet dynamics by fluctuating hydrodynamics

We simulate the mesoscopic dynamics of droplets formed by phase-separated fluids at nanometer scales where thermal fluctuations are significant. Both spherical droplets fully immersed in a second fluid and sessile droplets which are also in contact with a solid surface are studied. Our model combines a Cahn–Hilliard formulation with incompressible fluctuating hydrodynamics; for sessile droplets, the fluid–solid contact angle is specified as a boundary condition. Deterministic simulations with an applied body force are used to measure the droplets' mobility from which a diffusion coefficient is obtained using the Einstein relation. Stochastic simulations are independently used to obtain a diffusion coefficient from a linear fit of the variance of a droplet's position with time. In some scenarios, these two measurements give the same value but not in the case of a spherical droplet initialized near a slip wall or in the case of sessile droplets with large contact angles (≥90°) on both slip and no-slip surfaces.

Bell, John B↗

DOE OSTI2026

Computational capacity in hydrodynamic real-time hybrid simulation applied to simulate the dynamic response of floating offshore wind turbines

Real-time hybrid simulation (RTHS) mitigates similitude distortions in model-scale tests of floating offshore wind turbines (FOWTs) by coupling physical experiments with numerical models in real time. The coupling requires faster-than-real-time numerical computations to satisfy temporal similitude with the physical experiment, presenting a bottleneck for using more complex numerical models in RTHS. This paper presents a hydrodynamic-RTHS (hydro-RTHS) framework for FOWTs that simulates the hydrodynamics physically and the aerodynamics numerically with sensor feedback from the physical testing. The framework adapts the three-loop hardware architecture to leverage greater computational resources and mitigate strict temporal requirements, enabling more computationally demanding numerical analyses in hydro-RTHS. The three-loop hardware architecture integrates multiple machines, each dedicated to either numerical analysis or RTHS controls, with a rate-transition algorithm to synchronize the tasks executed across the different machine processors. Virtual and physical tests verified and validated the hydro-RTHS framework, respectively. The ”virtual” tests, which approximates the physical domain numerically, verified the RTHS framework with respect to a numerical full-scale complete FOWT model simulated in the open-source software, OpenFAST. The virtual tests were able to maintain comparable control signals while enabling greater computational resources for the numerical calculations. Real-world physical tests demonstrated that the hydro-RTHS framework computes aerodynamic forces similar to the complete OpenFAST model, validating the hydro-RTHS framework using the three-loop hardware architecture. Findings show that the hydro-RTHS framework with the three-loop hardware architecture is computationally efficient, with reserve capacity to simulate more complex problems due to the customized software, hardware, and rate-transition algorithm.

17 WIND ENERGY↗

DOE OSTI2026

Optimizing time integration for accurate recovery of shockwave interface location in radiography

We present simulations and experiments of time integrated radiographic imaging of a moving 1D shock wave front and a quantitative method for determining the statistical error in locating the shock front as a function of integration time and noise in the radiograph. We discuss the trade-off between increasing motion blur, which leads to decreased shock front location certainty, and increasing signal-to-noise, which leads to improved image quality with increasing integration time. We find an optimum integration time between a short integration time, where noise limits the error, and a long integration time, where motion blurring limits the error. This methodology can be used to tune experimental configurations to obtain the highest quality radiograph for a given experimental configuration.

Bremsstrahlung↗

DOE OSTI2024

Using Hydrodynamic Similarity as a Verification Method for Impact Cratering Simulations in the FLAG Hydrocode

Hydrodynamic codes (hydrocodes) are common tools for modeling hypervelocity impacts to provide insight into the physical phenomenon. Hydrocodes can simulate impacts from micrometer to kilometer spatial scales and reach impact velocities difficult to achieve in experimental settings. However, numerical models are approximations, and demonstrating that a numerical method is capable of providing physical results for these models is essential. In this work, we employ a hydrocode verification technique that leverages hydrodynamic similarity, a mathematical property of the conservation equations of fluid mechanics that form the basis for hydrocode models. Using the FLAG hydrocode, we simulate aluminum (Al) and basalt projectiles and targets at spatial scales spanning 7 orders of magnitude (hundreds of micrometers to kilometers). These materials were chosen because Al-6061 is a common material in spacecraft and satellites and basalt is a useful approximation of rocky astronomical bodies. Our results show that hydrodynamic similarity holds for each material model used and across spatial scales. We show that under certain conditions hydrodynamic similarity can apply in the presence of gravity and that similarity does not hold in the presence of strength models. We conclude that the FLAG hydrocode preserves important mathematical properties of fluid dynamics in hypervelocity impacts of Al-6061 and basalt.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2021

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↗

DOE OSTI2024

Modeling Dust Production, Growth, and Destruction in Reionization-era Galaxies with the CROC Simulations. II. Predicting the Dust Content of High-redshift Galaxies

Abstract We model the interstellar dust content of the reionization era with a suite of cosmological, fluid-dynamical simulations of galaxies with stellar masses ranging from ∼10 5 to 10 9 M ⊙ in the first 1.2 Gyr of the Universe. We use a post-processing method that accounts for dust creation and destruction processes, allowing us to systematically vary the parameters of these processes to test whether dust-dependent observable quantities of galaxies at these epochs could be useful for placing constraints on dust physics. We then forward model observable properties of these galaxies to compare to existing data. We find that we are unable to simultaneously match existing observational constraints with any one set of model parameters. Specifically, the models that predict the largest dust masses D / Z ≳ 0.1 at z = 5—because of high assumed production yields and/or efficient growth via accretion in the interstellar medium—are preferred by constraints on total dust mass and infrared (IR) luminosities, but these models produce far too much attenuation in the ultraviolet (UV), preventing them from matching observations of β UV . To investigate this discrepancy, we analyze the relative spatial distribution of stars and dust as probed by IR and UV emission, which appear to exhibit overly symmetric morphologies compared to existing data, likely due to the limitations of the stellar feedback model used in the simulations. Our results indicate that the observable properties of the dust distribution in high redshift galaxies are a particularly strong test of stellar feedback.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2022

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↗

DOE OSTI2025

Observation of kinetic mix enhancement in thin-shell OMEGA implosions

Recent separated reactant experiments for thin-shell (6 µ⁢m) shock-driven implosions on OMEGA have demonstrated significant mix from a buried deuterated layer of the shell into the hot spot. Time resolved D 3 He-p reaction history data demonstrate a (50 ± 20)⁢ ps shift earlier in peak nuclear emission for separated reactant experiments relative to control, in contrast to past experimental data for thicker, 20 µ⁢m shells with no laser burn through that show a 75 ps delay due to the time required for hydrodynamic instabilities to develop. This contrast suggests that the mix mechanism was not hydrodynamic. Ion kinetic simulations utilizing fall line analyses show much closer agreement with mix yield and temperature than diffusion models, predicting a D 3 He-p mix yield of 1.7 × 10 9 as compared to the experimental value of 9.3⁢ (±2.1) × 10 8 . This is three orders of magnitude closer than the fall line analysis from a hydrodynamic simulation with an inline diffusive mix model, which suggests minimal mix and D 3 He-p yields of 5×10 5 . This makes kinetic mechanisms the only feasible explanation for the mix seen, demonstrating impact of a non-standard mix mechanism. An analytical model of this kinetic mix mechanism suggests that it can remain significant in situations when the shell expands significantly to low densities, and diffusive models predict negligible mix. Finally, kinetic mix will impact multiple types of high energy density, laser-driven fusion experiments including high-adiabat direct drive cryoexperiments, nuclear cross section experiments, and thin-shell polar direct drive experiments used to tune heat conduction models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Weak lensing combined with the kinetic Sunyaev–Zel’dovich effect: a study of baryonic feedback

ABSTRACT Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmology Telescope DR5 kinematic Sunyaev–Zel’dovich (kSZ) to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model, ‘baryonification’. First, using WL only, we compare the $S_8$ constraints using baryonification to a simulation-calibrated halo model, a simulation-based emulator model, and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of $S_8$ and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters, with the joint WL + kSZ analysis constraining $S_8=0.823^{+0.019}_{-0.020}$. We measure the suppression of the non-linear matter power spectrum using WL + kSZ and constrain a mean feedback scenario that is more extreme than the predictions from most hydrodynamical simulations. We constrain the baryon fractions and the gas mass fractions and find them to be generally lower than inferred from X-ray observations and simulation predictions. We conclude that the WL + kSZ measurements provide a new and complementary benchmark for building a coherent picture of the impact of gas around galaxies across observations.

79 ASTRONOMY AND ASTROPHYSICS↗