Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Hydrodynamic models”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Modeling of Optoelectronic Devices

Ultrafast modulation of semiconductor quantum well (QW) laser is of technological importance for information technology. Improvement by order(s) of magnitude in data transfer rate is possible as terahertz (THz) radiation is available for heating the laser at picosecond time scale. Optical gain modulation in the QW is achieved via temperature modulation of electron-hole plasma (EHP). Applications include free-space THz communication, optical switching, and pulse generation. The EHP in the semiconductor QW is described with a two-band model. Semiconductor Bloch equations with many-body effects are used to derive a hydrodynamical model for the active QW region. Because of ultrafast carrier-carrier scatterings in the order of 50 fs, EHP follows quasiequilibrium Fermi-Dirac distributions and THz field interacts incoherently with it. Carrier-longitudinal optical (LO) phonon scatterings and coherent laser-EHP interaction are treated microscopically in our physical model. A set of hydrodynamical equations for plasma density, temperature, and laser envelop amplitude are derived and Runge-Kutta method is adopted for numerical simulation. A typical 8 nm GaAs/Al(0.3)Ga(0.7) As single QW at 300 K is used. Additional information is contained in the original extended abstract.

Li, Jian-Zhong↗
DOE OSTI2025

Optimal Uses of Magnetic Fields for Indirect-Drive Inertial Fusion

This project explored how applied magnetic fields can improve inertial confinement fusion (ICF), specifically the indirect-drive approach that uses a hohlraum. This has been proposed for several decades as potentially beneficial, due to thermal insulation (reduced losses) from the imploded hotspot. We performed the most advanced radiation-magneto-hydrodynamic modeling to date of magnetized ICF designs in the ignition regime. We found that adding technologically feasible fields up to 60 – 70 Tesla could increase the fusion yield of current igniting designs for the National Ignition Facility (NIF) by up to 8x. Also, simulations show that in certain cases relatively small fields of 3 – 5 Tesla could double the yield, and be implemented at much lower cost. Early work on re-optimizing NIF designs with magnetic fields, namely by using a thicker ablator with more mass remaining, could increase the yield of a sub-ignition target by 18x and bring it into the ignition regime. A separate benefit of magnetization besides reduced thermal loss is reduced hydrodynamic instability. Modeling work under this project shows this could be significant, though early experiments at NIF and the Omega Laser proved inconclusive. We designed and proposed an improved NIF experiment on magnetized mix, based on a large-amplitude imposed perturbation. The project also supported basic physics research into magnetized laser-plasma interactions, namely cross-beam energy transfer, both with experiments at Omega and theory / modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Collapse models for dark interstellar clouds

Properties of self-consistent numerical hydrodynamic models are compared with observed properties of several dark clouds. The results are consistent with the view that these clouds are undergoing gravitational collapse. The clouds appear to have evolved from similar initial states and to have ages comparable to their free-fall times. Derived cloud masses range between 10 and 1000 solar masses, correlating with cloud size. The models reproduce observed cloud-to-cloud variations in the (C-13)O abundance, and they offer additional evidence that the (C-13)O abundance varies within individual clouds.

Villere, K. R.↗
DOE OSTI2026

A hybrid CNN-LSTM surrogate model for hyper-resolution spatiotemporal flood forecasting in Norfolk, Virginia

Study region: Norfolk, Virginia, United States Study focus: Accurate and timely flood forecasting is essential for enhancing resilience in coastal urban areas in the context of increasing frequency and intensity of rainfall, sea level rise and urbanization. This study presents a hybrid deep learning-based surrogate model that integrates Convolutional Neural Networks (CNN) and Long Short-Term Memory (LSTM) networks to enable real-time spatiotemporal flood forecasting. The model leverages CNN to capture spatial features from inputs such as elevation and Topographic Wetness Index (TWI), while LSTM processes time-series inputs of rainfall and tide data to capture temporal features. New hydrologic insights for the region: The hybrid CNN-LSTM model was trained using the physics-based hydrodynamic model simulations obtained from the Two-dimensional Unsteady FLOW (TUFLOW) model for Norfolk, Virginia, and achieved high predictive accuracy across diverse flood-prone areas. The reduced computational time from four to six hours using TUFLOW to 3.2 min per event using CNN-LSTM enables rapid flood inundation mapping and early warning applications. The model effectively captured both spatial flood extents and their temporal evolution across different flooding scenarios, providing forecasts at a 2.5-m spatial resolution and 15-min temporal resolution and a one-hour-ahead prediction horizon. While challenges remain in terms of transferability to new regions and real-time data assimilation, this approach demonstrates strong potential for supporting operational flood risk management in coastal urban environments.

Coastal urban flooding↗
NASA NTRS1995

What is the source of observed annual variations in plasmaspheric density?

Plasmaspheric densities have been observed previously to be higher in December than in June, with the ratio varying between 1.5 and 3.0 and with larger variations at lower L shells. In order to search for the cause of the observed annual variations, we have modeled plasmaspheric density, using a time-dependent hydrodynamic model. On an L = 2 field line with geomagnetic longitude equal to 300 deg, the modeled plasmaspheric densities were a factor of 1.5 times higher in December than in June. The modeled December to June density ratio was found to increase slightly with L shell, in contrast to observations; this discrepancy may be due to the fact that outer plasmaspheric flux tubes are never completely full. In addition, for an L = 2 field line with geomagnetic longitude equal to 120 deg, the modeled plasmaspheric density was higher in June than in December by a factor of about 1.2. Various numerical tests were also performed in order to examine the sensitivity of plasmaspheric density to various parameters. In particular, a large vertical neutral wind was applied in order to raise the O(+) profile which had the effect of raising plasmaspheric density by a factor of 6. This in conjunction with a theoretical analysis suggests that plasmaspheric density levels are very sensitive to O(+) levels in the upper ionosphere. We conclude that annual variations in plasmaspheric density are due to similar variations in ionospheric O(+).

Guiter, S. M.↗
NASA NTRS1995

Plasmasphere Modeling with Ring Current Heating

Coulomb collisions between ring current ions and the thermal plasma in the plasmasphere will heat the plasmaspheric electrons and ions. During a storm such heating would lead to significant changes in the temperature and density of the thermal plasma. This was modeled using a time- dependent, one-stream hydrodynamic model for plasmaspheric flows, in which the model flux tube is connected to the ionosphere. The model simultaneously solves the coupled continuity, momentum, and energy equations of a two-ion (H(+) and O(+) quasineutral, currentless plasma. Heating rates due to collisions with ring current ions were calculated along the field line using a kinetic ring current model. First, diurnally reproducible results were found assuming only photoelectron heating of the thermal electrons. Then results were found with heating of the H(+) ions by the ring current during the recovery phase of a magnetic storm.

Guiter, S. M.↗
DOE OSTI2026

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

A three-dimensional time-dependent model of the plasmasphere

A three-dimensional, time-dependent, nonlinear, hydrodynamic model of the plasmasphere has been developed which includes the self-consistent coupling of conjugate hemispheres within the plasmasphere as well as the effects of cross-L drifts due to convection electric fields. The inner region of the plasmasphere was modeled from an L of 1.5 to an L of 3.5 to 4.5, depending on local time. In this first study of the global plasmasphere, the model was run for solar minimum conditions until diurnally reproducible results were obtained, indicating that the plasmasphere was fully filled. The results of the model were then compared with applicable measurements of the plasmasphere. It was found that the model tends to overestimate densities somewhat, although it is within a factor of 2 of an average of June and December whistler observations. The model predicts densities within a factor of 2 of most of the satellite observations as well. As anticipated, the diurnal variation of the plasmasphere was found to depend on magnetic latitude. At low latitudes the diurnal variation in density was relatively small, with the largest densities occurring in the afternoon time sector. However, near the plasmapause, the effects of changes in volume of drifting tubes of plasma due to cross-L drifts led to a factor of 3 variation in equatorial density, with the highest densities occurring near local midnight where the volume was lowest. Various limitations of the current model and suggested improvements are also discussed.

Rasmussen, C. E.↗
NASA NTRS1993

Numerical modeling of the interstellar medium in galactic disks

We have been developing detailed hydrodynamic models of the global interstellar medium in the hope of understanding the mass and volume occupied by various phases, as well as their structure and kinematics. In our model, the gas is modeled by one fluid while representative Pop 1 stars are modeled by a second fluid. The two fluids are coupled in that the gas forms into stars at a rate given by a Schmidt law while stellar mass loss returns matter into the gas phase (on a time scale of 100 Myr). Also, the stars heat the gas through stellar winds and the gas cools through optically thin radiation. The time behavior of these two fluids is studied in two spatial dimensions with the Eulerian finite difference numerical hydrodynamic code Zen. The two spatial dimensions are along the plane of a disk (x, total length of 2 kpc) and perpendicular to the disk (z, total height of +/- 15 kpc) and a galactic gravitational field in the z direction, typical of that at the solar circle, is imposed upon the simulation; self-gravity and rotation are absent. For the boundary conditions, outflow is permitted at the top and bottom of the grid (z = +/- 15 kpc) while periodic boundary conditions are imposed upon left and right sides of the grid. As initial conditions, we assumed a gaseous distribution like that seen for the H1 by earlier researchers, although the results are insensitive to the initial conditions. We have run simulations in which the heating due to stars, parameterized as a stellar wind velocity, a, is varied from low (a = 150 km/s), to intermediate (a = 300 km/s), to high (a = 600 km/s). Since the intermediate case is roughly equivalent to the Galactic energy injection rate from supernovae, this summary will concentrate on results from this simulation.

Rosen, A.↗
DOE OSTI2025

Life Cycle Cost Analysis of Prestressed Concrete Poles Subjected to Wind, Surges, and Waves

Prestressed concrete (PC) poles are becoming popular choices to support coastal power transmission systems. However, the existing literature does not offer a detailed analysis of the effectiveness of PC poles in terms of long-term vulnerabilities and the direct and indirect costs. This is due to (1) lack of fragility models for PC poles and (2) lack of probabilistic wind, storm surge, and wave models in coastal settings. In this study, we address these gaps through a series of Monte Carlo simulations to estimate fragility of PC poles as a function of age and hazard (wind, surges, and waves) intensity, and the development of a probabilistic hazard model based on 10,000 years of synthetic tropical cyclone data. The probabilistic hazard model is used in conjunction with high-resolution hydrodynamic models to generate realizations of coastal wind, storm surges, and waves for the Louisiana and Mississippi coasts. A comprehensive life cycle cost analysis for a service life of 70 years considering direct and indirect losses is conducted to compare the performance of a transmission line located in Pascagoula, Mississippi, when wood poles are replaced by PC poles. Results showed that aging has a minor effect on the reliability of PC poles, highlighting the advantages of replacing wood poles with PC poles, especially in coastal areas. In addition, PC poles are significantly more cost-effective compared with wood poles over their life cycle, leading to an estimated saving of $11.55 million (68.17% reduction). The results of this study provide key insight to inform decision-making processes to keep the coastal power grids resilient and cost-effective against future storm hazards.

natural disasters↗
DOE OSTI2025

Ponderosa pine hydraulic stress predicts more extreme wildfire behavior under future conditions in Bandelier National Monument, New Mexico

Background Live fuel moisture contributes to wildfire spread and reflects plant stress and physiological traits. The anticipated change in live fuel moisture under future conditions is likely non-linear, owing to physiological plant thresholds in water hydraulics. We constructed a mechanistic model of live fuel moisture’s response to water stress to understand the impact of future climate on live fuel moisture. We first gathered data on plant physiology and live fuel moisture for Pinus ponderosa at Bandelier National Monument, NM, USA, and modeled their relationship. We then parameterized a mechanistic plant hydrodynamics model (FATES-HYDRO) to simulate changes in plant stress and a statistical model to simulate the resulting impact on live fuel moisture. We then simulated FATES-HYDRO under future climate anomalies (SSP2-4.5 and SSP5-8.5: 2080–2100) to understand the change in plant stress and estimate its impact on live fuel moisture. Results We found that the number of days below crucial thresholds of live fuel moisture (100% and 79%) increased from contemporary levels (< 100%: 72 days, < 79%: 1.4 days) under SSP2-4.5 (< 100%: 185 days, < 79%: 10.2 day) and increased exponentially under SSP5-8.5 (< 100%: 215 days, < 79%: 65 days). We found that gross primary productivity decreased under both future climate scenarios (contemporary: 336 g C m −2 , SSP2-4.5: 203 g C m −2 , SSP5-8.5: 243 g C m −2 ); however, spring productivity increased under SSP5-8.5, potentially altering fuel loading. We additionally see a potentially lethal loss of conductivity in hydraulic P. ponderosa under SSP5-8.5. Conclusions The overall increase in plant water stress (as represented by loss of hydraulic conductivity and leaf water potential) and lower live fuel moisture appear to be driven by reduced precipitation during late summer monsoons typical of the region, extending the fire season later in the year. We predict increasing variability in the P. ponderosa wildfire regime under both climate projections driven by changing productivity, rising mortality, and an overall decrease in live fuel moisture.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Kinetic modeling of O(+) upflows resulting from E x B convection heating in the high-latitude F region ionosphere

We report here the results of modeling work aimed at understanding the development of ionospheric O(+) field-aligned upflows that develop in response to high-latitude E x B drift induced frictional heating. The model used is a collisional semikinetic model which includes ion-neutral resonant charge exchange and polarization collisions as well as Coulomb self-collisions. It also includes the process of chemical removal of O(+) as well as all of the macroscopic forces: ambipolar electric, gravity, magnetic mirror, and centripetal. Model results show the development of several types of non-Maxwellian velocity distributions including toroids at low altitude, distributions with large heat flow in the perpendicular component at intermediate altitudes, and distributions with a separate upflowing population or upward superthermal tail at high altitudes. Whenever the convection electric field increases from a small value (less than 25 mV/m) to a large value (100-200 mV/m) in 6 min or less large upflows develop with parallel drift speeds which peak (below 1000 km) at values between 500 m/s and 2 km/s, parallel fluxes which peak between 6.0 x 10(exp 8) and 3.2 x 10(exp 9)/sq cm/s, and parallel per particle heat flows which peak between 8.0 x 10(exp -9) and 8.0 x 10(exp -8) ergs cm/s. The higher values in these ranges occur for a cooler neutral atmosphere, with a larger convection electric field that is turned on quickly. The model produces field-aligned O(+) flow speeds that are larger than those produced by a 20-moment generalized transport model but smaller then those produced by an isotropic hydrodynamic model for comparable values of the convection turn on times. The model results compare favorably with some topside satellite and radar data.

Wilson, G. R.↗
NASA NTRS1988

Model structure of a cosmic-ray mediated stellar or solar wind

An idealized hydrodynamic model is presented for the mediation of a free-streaming stellar wind by galactic cosmic rays or energetic particles accelerated at the stellar wind termination shock. The spherically-symmetric stellar wind is taken to be cold; the only body force is the cosmic ray pressure gradient. The cosmic rays are treated as a massless fluid with an effective mean diffusion coefficient k proportional to radial distance r. The structure of the governing equations is investigated both analytically and numerically. Solutions for a range of values of k are presented which describe the deceleration of the stellar wind and a transition to nearly incompressible flow and constant cosmic ray pressure at large r. In the limit of small k the transition steepens to a strong stellar wind termination shock. For large k the stellar wind is decelerated gradually with no shock transition. It is argued that the solutions provide a simple model for the mediation of the solar wind by interstellar ions as both pickup ions and the cosmic ray anomalous component which together dominate the pressure of the solar wind at large r.

Lee, M. A.↗
NASA NTRS2019

Estuarine Dissolved Organic Carbon Flux from Space: With Application to Chesapeake and Delaware Bays

This study uses a neural network model trained with in situ data, combined with satellite data and hydrodynamic model products, to compute the daily estuarine export of dissolved organic carbon (DOC) at the mouths of Chesapeake Bay (CB) and Delaware Bay (DB) from 2007 to 2011. Both bays show large flux variability with highest fluxes in spring and lowest in fall as well as interannual flux variability (0.18 and 0.27 Tg C/year in 2008 and 2010 for CB; 0.04 and 0.09 Tg C/year in 2008 and 2011 for DB). Based on previous estimates of total organic carbon (TOCexp) exported by all Mid-Atlantic Bight estuaries (1.2 Tg C/year), the DOC export (CB + DB) of 0.3 Tg C/year estimated here corresponds to 25% of the TOCexp. Spatial and temporal covariations of velocity and DOC concentration provide contributions to the flux, with larger spatial influence. Differences in the discharge of fresh water into the bays (74 billion m(exp3)/year for CB and 21 billion m9exp3)/year for DB) and their geomorphologies are major drivers of the differences in DOC fluxes for these two systems. Terrestrial DOC inputs are similar to the export of DOC at the bay mouths at annual and longer time scales but diverge significantly at shorter time scales (days to months). Future efforts will expand to the Mid-Atlantic Bight and Gulf of Maine, and its major rivers and estuaries, in combination with coupled terrestrial-estuarine-ocean biogeochemical models that include effects of climate change, such as warming and CO2 increase.

Signorini, Sergio↗
NASA NTRS2019

Estuarine Dissolved Organic Carbon Flux From Space: With Application to Chesapeake and Delaware Bays

This study uses a neural network model trained with in situ data, combined with satellite data and hydrodynamic model products, to compute the daily estuarine export of dissolved organic carbon (DOC) at the mouths of Chesapeake Bay (CB) and Delaware Bay (DB) from 2007 to 2011. Both bays show large flux variability with highest fluxes in spring and lowest in fall as well as interannual flux variability (0.18 and 0.27 Tg C/year in 2008 and 2010 for CB; 0.04 and 0.09 Tg C/year in 2008 and 2011 for DB). Based on previous estimates of total organic carbon (TOCexp) exported by all Mid-Atlantic Bight estuaries (1.2 Tg C/year), the DOC export (CB + DB) of 0.3 Tg C/year estimated here corresponds to 25% of the TOCexp. Spatial and temporal covariations of velocity and DOC concentration provide contributions to the flux, with larger spatial influence. Differences in the discharge of fresh water into the bays (74 billion cu. m/year for CB and 21 billion cu. m/year for DB) and their geomorphologies are major drivers of the differences in DOC fluxes for these two systems. Terrestrial DOC inputs are similar to the export of DOC at the bay mouths at annual and longer time scales but diverge significantly at shorter time scales (days to months). Future efforts will expand to the Mid-Atlantic Bight and Gulf of Maine, and its major rivers and estuaries, in combination with coupled terrestrial-estuarine-ocean biogeochemical models that include effects of climate change, such as warming and CO2 increase.

Estuarine DOC export↗
NASA NTRS1995

Late-time spectra and type Ia supernova models: New clues from the Hubble Space Telescope

Calculated late-time spectra of two classical hydrodynamical models for Type Ia supernovae (deflagration model W7 of Nomoto, Thielemann, & Yokoi, and delayed detonation model DD4 of Woosley & Weaver) are compared with observations of SN 1992A and other spectroscopically normal SNe Ia. An important new piece of information is provided by observations done with the Hubble Space Telescope (HST) which cover the ultraviolet range at the nebular phase of a SN Ia: SN 1992A in NGC 1380. For the first time a picture of SN Ia emission from the ultraviolet through the optical is obtained at these phases. Predictions of the classical model (W7 and DD4) are compared with the observed spectrum of SN 1992A and with the optical spectra of SN 1989M in NGC 4579 and SN 1990N in NGC 4639 at similar epochs. The absolute B and V magnitudes of the models are also estimated at these late phases. Taken at face value the nebular spectra of these 'classical' models are more consistent with the long extragalactic distance scale, pointing to distances to NGC 4579 around 21 +/- 3 Mpc and a slightly larger distance, 22 +/- 3 Mpc, to NGC 4639, on the back side of the Virgo Cluster. However, the calculated Fe(+3) luminosity as predicted from the models exceeds the observed limit from the HST data of SN 1992A. Other differences in the ratios of the line intensities between calculated and observed spectra, show some disagreement with the observed spectra at the nebular phases. They may not be the best choice for spectroscopically normal SNe Ia, and their use as an independent calibration of the extragalactic distance scale should be viewed with caution.

Ruiz-Lapuente, P.↗
NASA NTRS1978

X-ray analysis of a polar plume

Polar plumes have been detected in X-rays with the Skylab S-054 experiment. Observations of one well-observed plume are analyzed and the plasma pressure is determined. An isothermal model in hydrostatic equilibrium is found to require unacceptably low temperatures. A hydrodynamic model is proposed which is consistent with earlier white light and EUV observations. Calculations indicate that the total outward mass flux in polar plumes is comparable to that in high speed solar wind streams expected from a polar coronal hole.

Ahmad, I. A.↗