Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic modeling”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Gamma rays, X-rays, and optical light from the cobalt and the neutron star in SN 1987A

Recent developments in modeling the X-ray and gamma-ray emission from SN 1987A are discussed by taking into account both the decaying cobalt and the buried neutron star. The light curve and the spectra evolution of X-rays and gamma-rays are well modeled up to day of about 300 if mixing of Co-56 into hydrogen-rich envelope is assumed. However, the 16-28 keV flux observed by Ginga declines very slowly, whereas the spherical mixing model predicts that the flux should have decreased by a large factor at t greater than 300d. It is shown that this problem can be solved if the photoelectric absorption of X-rays is effectively reduced as a result of the formation of chemically inhomogeneous clumps. Based on the adopted hydrodynamical model and the abundance distribution, predictions are offered for future optical, X-ray, and gamma-ray light curves by taking into account other radioactive sources and various types of the central source, e.g., a buried neutron star accreting the reinfalling material or an isolated pulsar.

Kumagai, Shiomi↗
NASA NTRS1978

A three-dimensional model of corotating streams in the solar wind. I - Theoretical foundations

The paper is concerned with the development of the theoretical and mathematical background pertinent to the study of steady, corotating solar wind structure in all three spatial dimensions. The dynamical evolution of the plasma in interplanetary space (defined as the region beyond roughly 35 Rs where the flow is supersonic) is approximately described by the nonlinear, single-fluid, polytropic magnetohydrodynamic or hydrodynamic equations. Efficient numerical techniques are outlined for solving this complex system of coupled, hyperbolic partial differential equations. The present formulation is inviscid and nonmagnetic, but the methods used allow for the potential inclusion of both features with only modest modifications. A simple, highly idealized hydrodynamic model stream is examined to illustrate the fundamental processes involved in the three-dimensional dynamics of stream evolution. It is found that spatial variations in the rotational stream interaction mechanism produce small nonradial flows on a global scale that lead to the transport of mass, energy, and momentum away from regions of relative compression and into regions of relative rarefaction. Comparison with simpler models demonstrates the essential nonlinear, multidimensional nature of the interplanetary dynamics.

Pizzo, V.↗
DOE OSTI2023

Measurement of φ-meson production in Cu + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV and U + U collisions at $\sqrt{s_{NN}}$ = 193 GeV

The PHENIX experiment reports systematic measurements at the Relativistic Heavy Ion Collider of φ-meson production in asymmetric Cu + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV and in U + U collisions at $\sqrt{s_{NN}}$ = 193 GeV. Measurements were performed via the φ → K + K – decay channel at midrapidity |η| < 0.35. Features of φ- meson production measured in Cu + Cu, Cu + Au, Au + Au, and U + U collisions were found to not depend on the collision geometry, which was expected because the yields are averaged over the azimuthal angle and follow the expected scaling with nuclear-overlap size. The elliptic flow of the φ meson in Cu + Au, Au + Au, and U + U collisions scales with second-order-participant eccentricity and the length scale of the nuclear-overlap region (estimated with the number of participating nucleons). At moderate $p_T$ , φ-meson production measured in Cu + Au and U + U collisions is consistent with coalescence-model predictions, whereas at high $p_T$ the production is in agreement with expectations for in-medium energy loss of parent partons prior to their fragmentation. Finally, the elliptic flow for φ mesons measured in Cu + Au and U + U collisions is well described by a (2+1)-dimensional viscous-hydrodynamic model with specific-shear viscosity η/s = 1/4π.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Shell-model studies of the astrophysical r p -process reactions S 34 ( p , γ ) Cl 35 and Cl 34 g , m ( p , γ ) Ar 35

Background: Dust grains condensed in the outflows of presolar classical novae should have been present in the protosolar nebula. Candidates for such presolar nova grains have been found in primitive meteorites and can in principle be identified by their isotopic ratios, but the ratios predicted by state-of-the-art one-dimensional hydrodynamic models are uncertain due to nuclear-physics uncertainties. Purpose: To theoretically calculate the thermonuclear rates and uncertainties of the S 34 ( p , γ ) Cl 35 and Cl 34 g , m ( p , γ ) Ar 35 reactions and investigate their impacts on the predicted 34 S/ 32 S isotopic ratio for presolar nova grains. Method: A shell-model approach in a (0 +1 ) ℏ ω model space was used to calculate the properties of resonances in the 34 S(p ,γ ) 35 Cl and 34 g,m Cl (p ,γ ) 35 Ar reactions and their thermonuclear rates. Uncertainties were estimated using a Monte Carlo method. The implications of these rates and their uncertainties on sulfur isotopic nova yields were investigated using a postprocessing nucleosynthesis code. The rates for transitions from the ground state of 34 Cl as well as from the isomeric first excited state of 34 Cl were explicitly calculated. Results: At energies in the resonance region near the proton-emission threshold, many negative-parity states appear. Energies, spectroscopic factors, and proton-decay widths are reported. The resulting thermonuclear rates are compared with previous determinations. Conclusions: The shell-model calculations alone are sufficient to constrain the variation of the 34 S/ 32 S ratios to within about 30%. Uncertainties associated with other reactions must also be considered, but in general we find that the 34 S/ 32 S ratios are not a robust diagnostic to clearly identify presolar grains made from nova ejecta.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Model Advancements to Enable Impact Analysis of Climate Change on Streamflow Temperature

With support from the Department of Energy’s Water Power Technologies Office, Pacific Northwest National Laboratory (PNNL) has developed new tools that incorporate cutting-edge climate and hydrological science capabilities to assess the potential long-term impacts of future climate conditions on unregulated streamflow and water temperature within watershed-river-reservoir systems. The objectives of this project were achieved by enhancing key hydrologic and hydrodynamic models and transferring them to a high-performance computing environment to provide a high-spatiotemporal resolution, multi-scale modeling framework. The new modeling framework has the potential to quantify risks related climate change impacts on runoff, unregulated streamflow, and water temperature. Initial development and demonstration of the modeling framework was conducted under historical and future climate conditions in the Columbia River Basin in the Pacific Northwest and the Connecticut River Basin in New England.

13 HYDRO ENERGY↗
NASA NTRS2023

Towards an Indian Land Data Assimilation System (ILDAS): A Coupled Hydrologic-Hydraulic System for Water Balance Assessments

Effective management of water resources requires reliable estimates of land surface states and fluxes, including water balance components. But most land surface models run in uncoupled mode and do not produce river discharge at catchment scales to be useful for water resources management applications. Such integrated systems are also rare over India where hydrometeorological extremes have wreaked havoc on the economy and people. So, an Indian Land Data Assimilation System (ILDAS) with a coupled land surface and a hydrodynamic model has been developed and driven by multiple meteorological forcings (0.1°, daily) to estimate land surface states, channel discharge, and floodplain inundation. ILDAS benefits from an integrated framework as well as the largest suite of observation records collected over India and has been used to produce a reanalysis product for 1981–2021 using four forcing datasets, namely, Modern-Era Retrospective Analysis for Research and Applications, Version 2 (MERRA-2), Climate Hazards Group InfraRed Precipitation with Station data (CHIRPS), ECMWF’s ERA-5, and Indian Meteorological Department (IMD) gridded precipitation. We assessed the uncertainty and bias in these precipitation datasets and validated all major components of the terrestrial water balance, i.e., surface runoff, soil moisture, terrestrial water storage anomalies, evapotranspiration, and streamflow, against a combination of satellite and in situ observation datasets. Our assessment shows that ILDAS can represent the hydrological processes reasonably well over the Indian landmass with IMD precipitation showing the best relative performance. Evaluation against ESA-CCI soil moisture shows that MERRA-2 based estimates outperform the others, whereas ERA-5 performs best in simulating evapotranspiration when evaluated against MODIS ET. Evaluations against observed records show that CHIRPS-based estimates have the highest performance in reconstructing surface runoff and streamflow. Once operational, this system will be useful for supporting transboundary water management decision making in the region.

Indian Land Data Assimilation System (ILDAS)↗
DOE OSTI2021

Numerical Model of IProTech PIP WEC Device

iProTech PIP wave energy converter (WEC) is a slack moored, single hull device with no moving parts in the water, joints or bearings. This submission includes data of the simulation, reports, and code for the iProTech PIP (WEC) project. The organization of the data included in the provided archive is detailed below and in the data description of the archive. The data teamer-iprotech-nrel folder includes and explains matlab and python code developed to hydrodynamically model the PIP WEC device in WEC-Sim. The subfolders cover the following steps: 1) report: explanatory information on device geometry 2) pip_mesher: python code to generate mesh panels from device profile data 3) wec-sim_models: matlab code to run WEC-Sim The data uploaded is a snapshot as of 11/02/2121 of code residing in a Github repository administered by David Ogden of NREL.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Physics-driven Explosions of Stripped High-mass Stars: Synthetic Light Curves and Spectra of Stripped-envelope Supernovae with Broad Light Curves

Stripped-envelope supernovae (SESNe) represent a significant fraction of core-collapse supernovae, arising from massive stars that have shed their hydrogen and, in some cases, helium envelopes. The origins and explosion mechanisms of SESNe remain a topic of active investigation. In this work, we employ radiative-transfer simulations to model the light curves and spectra of a set of explosions of single, solar-metallicity, massive Wolf–Rayet stars with ejecta masses ranging from 4 to 11 M ⊙ , which were computed from a turbulence-aided and neutrino-driven explosion mechanism. We analyze these synthetic observables to explore the impact of varying ejecta mass and helium content on observable features. We find that the light curve shape of these progenitors with high ejecta masses is consistent with observed SESNe with broad light curves but not the peak luminosities. The commonly used analytic formula based on rising bolometric light curves overestimates the ejecta mass of these high-initial-mass progenitor explosions by a factor of up to 2.6. In contrast, the calibrated method by Haynie et al., which relies on late-time decay tails, reduces uncertainties to an average of 20% within the calibrated ejecta mass range. Spectroscopically, the He I 1.083 μm line remains prominent even in models with as little as 0.02 M ⊙ of helium. However, the strength of the optical He I lines is not directly proportional to the helium mass but instead depends on a complex interplay of factors such as the 56 Ni distribution, composition, and radiation field. Thus, producing realistic helium features requires detailed radiative transfer simulations for each new hydrodynamic model.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1993

Automata network models of galaxy evolution

Two ideas appear frequently in theories of star formation and galaxy evolution: (1) star formation is nonlocally excitatory, stimulating star formation in neighboring regions by propagation of a dense fragmenting shell or the compression of preexisting clouds; and (2) star formation is nonlocally inhibitory, making H2 regions and explosions which can create low-density and/or high temperature regions and increase the macroscopic velocity dispersion of the cloudy gas. Since it is not possible, given the present state of hydrodynamic modeling, to estimate whether one of these effects greatly dominates the other, it is of interest to investigate the predicted spatial pattern of star formation and its temporal behavior in simple models which incorporate both effects in a controlled manner. The present work presents preliminary results of such a study which is based on lattice galaxy models with various types of nonlocal inhibitory and excitatory couplings of the local SFR to the gas density, temperature, and velocity field meant to model a number of theoretical suggestions.

Chappell, David↗
DOE OSTI2021

Experimental observations of shock-wave-induced bubble collapse and hot-spot formation in nitromethane liquid explosive

We studied the collapse of individual helium gas bubbles in the homogeneous explosive nitromethane (NM) to investigate effects of hot-spot formation on the detonation process. A bubble was injected into a NM sample, and a shock wave from an explosive detonator compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and the pre- and post-bubble collapse processes. An infrared camera image showed the residual radiance temperature distribution after the bubble collapse, and an optical fiber pyrometer measured the time-resolved thermal radiance. We measured the optical spectra of light emitted from detonating NM without a bubble and from a collapsing bubble in shocked, undetonated NM. We estimated temperatures of the detonation fronts and of the hot spots formed by bubble collapse. To study the incipient detonation process, we performed all bubble collapse experiments at pressures below the threshold for creating a sustained detonation. Where the bubble collapsed, we observed an opaque, thermally emissive region believed to be chemical reaction products. Chemical reactions in NM can be produced with lower shock pressures (~1 GPa) when a helium bubble is present than without a bubble (~10 GPa). We used hydrodynamic modeling to predict shock wave propagation, extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. Simulations using a temperature-dependent Arrhenius burn model gave much better results than reactive burn models that depend only on pressure and density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Flow resistance of phloem sieve plates revisited using an experimental model

Sieve plate resistance has been identified as the key to understanding the efficiency of phloem transport of carbohydrates in trees. These plates with small holes connecting sieve tubes are responsible for the largest resistance to flow in the phloem. Their structure determines how fast sugars can be transported through the phloem with certain pressure differences and what the limits for phloem transport in different plants are. Because experiments with the phloem are very challenging, our understanding of sieve plate resistance is mostly based on anatomical studies and hydrodynamic modeling of flow through sieve plates. These models calculate the resistance of the entire sieve tube–sieve plate system using the Hagen–Poiseuille flow resistance for the sieve cell lumen and a combination of the Hagen–Poiseuille resistance and Sampson flow resistance through the sieve pores. The resistance of the entire sieve plate is calculated by summing the Hagen–Poiseuille and Sampson flow resistances of each pore. To test the validity of this model formulation, an experimental model with aspect parameters similar to phloem sieve tubes was built using polyvinyl chloride (PVC) piping and plastic straws of different diameters and lengths. This system was used to measure flow rates and calculate flow resistance at Reynolds numbers 0.5–300. The results suggest that the current models may significantly overestimate the flow resistance caused by sieve plates and that the resistance might be better described by formulations used for perforated plates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

QCD critical point from the Nambu–Jona-Lasino model with a scalar-vector interaction

We study the critical point in the QCD phase diagram in the Nambu–Jona-Lasino (NJL) model by including a scalar-vector coupled interaction. We find that varying the strength of this interaction, which has no effect on the vacuum properties of QCD, can significantly affect the location of the critical point in the QCD phase diagram, particularly the value of the critical temperature. This provides a convenient way to use the NJL-based transport or hydrodynamic model to extract information about the QCD phase diagram from relativistic heavy-ion collisions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Hydrological controls of a riparian wetland based on stable isotope data and model simulations

Isotopic evidence of groundwater and stream water is frequently used to investigate water exchanges with groundwater. Monthly sampling of rain, stream water, and groundwater was conducted at Tims Branch watershed in South Carolina for the oxygen and hydrogen stable isotope (δ 2 H and δ 18 O) measurement, as well as pH and oxidation–reduction potential (ORP). Together with a mass balance perspective, it was determined that it takes a few weeks to one month for groundwater in the hyporheic zone to fully exchange with stream water. From hydrodynamic modelling, we show that substantial (up to 70 %) groundwater exchange occurs at gaining and losing sites. Groundwater exfiltration, i.e. inflow into stream water, contributes up to 4 % to stream water, with the remainder from upstream exfiltration. A 2–4 % per day renewal rate of adjacent groundwater would indirectly indicate a groundwater residence time in the order of half a month to a full month (assuming either a well-mixed case or large dispersion rate in pulse flow case), in agreement with a greatly reduced variability of δ 2 H and δ 18 O of groundwater compared to stream water and rain. This reduced variability of stable isotope signal from groundwater confirms our hypothesis that riparian groundwater mixing at Tims Branch is more of a mixed type rather than a pulse flow type. As a result, a monthly time scale is sufficient for groundwater to become anoxic at exit points into stream water resulting in the episodic production of natural organic matter- and iron-rich flocs upon oxidation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1980

X-ray spectroscopy of supernova remnants

The information which can be obtained from X-ray spectroscopy of supernova remnants is considerd. The fitting of X-ray detector counts to models of the incident spectrum is discussed, and the types of thermal emission models generally employed are presented, including the power law, black body, and thermal bremsstrahlung models of the continua and models of the emission of a hot, optically thin plasma in collisional equilibrium. Observations of 12 supernova remnants made with the Solid State Spectrometer on board the Einstein Observatory are reported, and metal abundances inferred from the lines of the eight remnants showing thermal spectra are summarized. Questions raised by the failure to observe the overabundance of Fe predicted by stellar evolution and hydrodynamic modeling are discussed, and the need to develop more detailed models of the conditions in a supernova remnant in order to interpret the X-ray spectra is noted.

Szymkowiak, A. E.↗
NASA NTRS1987

Hydrodynamical constraints on cosmic-ray acceleration in relativistic shocks

A two-fluid hydrodynamical model governing the transport of cosmic rays in a relativistically moving background plasma is developed. The equations are used to discuss the time asymptotic structure of a relativistic, plane-parallel shock wave significantly modified by first-order Fermi acceleration of cosmic rays. The model allows for an anisotropic cosmic-ray pressure tensor with pressure components parallel and perpendicular to the shock normal. It is surmised that the well-known energy threshold problem for diffusive shock acceleration that downstream particles have a large enough velocity to overtake the shock and subsequently scatter in the upstream medium before returning to the downstream medium may be expressed hydrodynamically by the condition V less than V(cr), where V is the fluid velocity relative to the shock and V(cr) is the relativistic cosmic-ray sound speed. Astrophysical implications of the results are briefly discussed.

Webb, G. M.↗
DOE OSTI2022

Investigating laser ablated plume dynamics of carbon and aluminum targets

Recently acquired high-resolution images of nanosecond laser ablation plumes suggest a strong correlation between the internal plume structure and the type of material being ablated. However, the details of this relation are currently not well understood. In this work, we attempt to explore this correlation using a 2D radiation hydrodynamics model to study the dependence of internal plume structure formation on the ablation material. Here spatio-temporal emission maps and plume expansion velocities from experimental measurements are compared with the model predictions, including synthetic emission maps. The shape and expansion rate of an outer air plume region are found to be in good agreement for both carbon and aluminum, as are the inner material plume dynamics for carbon ablation. The largest disagreement is observed in the case of a polished aluminum target, where the chaotic inner plume features seen in the experimental images are not observed in the model. The possible physical mechanisms responsible for this discrepancy are discussed. This effort constitutes a continued development toward a predictive model of ablation plume dynamics and chemistry for various materials in extreme environments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Second-harmonic generation in nonlinear plasmonic lattices enhanced by quantum emitter gain medium

In this work, we report on a theoretical study of second-harmonic generation (SHG) in plasmonic nanostructures interacting with two-level quantum emitters (QEs) under incoherent energy pump. We generalize the driven-dissipative Tavis–Cummings model by introducing the anharmonic surface plasmon-polariton (SPP) mode coupled to QEs and examine physical properties of corresponding SPP-QE polariton states. Our calculations of the SHG efficiency for strong QE-SPP coupling demonstrate orders of magnitude enhancement facilitated by the polariton gain. We further discuss time-domain numerical simulations of SHG in a square lattice comprising Ag nanopillars coupled to QEs utilizing a fully vectorial nonperturbative nonlinear hydrodynamic model for conduction electrons coupled to Maxwell–Bloch equations for QEs. The simulations support the idea of gain enhanced SHG and show orders of magnitude increase in the SHG efficiency as the QEs are tuned in resonance with the lattice plasmon mode and brought above the population inversion threshold by incoherent pumping. By varying pump frequency and tuning QEs to a localized plasmon mode, we demonstrate further enhancement of the SHG efficiency facilitated by strong local electric fields. The incident light polarization dependence of the SHG is examined and related to the symmetries of participating plasmon modes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A HPC Theory-Guided Machine Learning Cyberinfrastructure for Communicating Hydrometeorological Data Across Scales

High-resolution predictions of hydrometeorological variables are critical for supporting hydropower generation decisions and flood control at hydroelectric power plants. Traditional climate and hydrologic models rely on the numerical simulation of detailed physical processes. Therefore, running these simulations is time-, labor-, and computation-intensive. Improving the spatial and temporal resolution in these modeling outputs could lead to cubic increases in both the simulation time and computational demands, rendering high-resolution hydrometeorological predictions expensive and impractical. Many past studies apply the super resolution (SR) technique to downscale climate models using deep learners. However, deep learners are deemed “black-boxes,” as their derivation processes from low-resolution outputs to high-resolution outputs are often hidden. Their results are difficult for domain scientists to interpret and validate. Thus, there is a need for an exploratory machine learning approach that can partially integrate domain-specific theory and knowledge into the data-driven mapping process between simulation outputs of different spatial scales. The domain-specific theory and knowledge can be incorporated into the data model through an inductive approach in which process-related environmental variables are used and analyzed as key drivers (i.e., environmental surrogates) to reflect the complex physical processes. Many of these variables, such as land use land cover, soil types, topography, digital elevation, air temperature, and various watershed characteristics, can be directly measured through sensors or remote sensing techniques. Additionally, SR applications that can downscale hydrological and hydrodynamics models to efficiently produce high-resolution (1 m) flood depth grids are still rare. Since the flood depth grid can be used to support critical decisions for flood control operation at hydroelectric power plants, it is crucial to enable an SR-based capability for interpolating high-resolution flood inundation maps.

13 HYDRO ENERGY↗