Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics 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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

On the X-ray spectrum of Kepler's supernova remnant

We have devised a method to do nonequilibrium ionization calculations on the results of two-dimensional hydrodynamical simulations, based on the algorithm of Hughes & Helfand (1985). We have calculated the ionization structure and X-ray emission for a two-dimensional numerical hydrodynamical simulation for the remnant of Kepler's supernova (SN); the hydrodynamical model was presented in a previous paper. In this model, the progenitor of Kepler's SN is assumed to have been a massive runaway star ejected from the Galactic plane. In its red supergiant stage, its dense stellar wind was distorted and compressed into a bow shock by the ram pressure of the tenuous interstellar medium. The subsequent interaction of the supernova ejecta with this asymmetric circumstellar matter produced a strongly asymmetric supernova remnant (SNR). In this paper, we present calculated X-ray spectra for this hydrodynamical model. A comparison with observations implies a moderate overabundance of Fe in Kepler's SNR (only 50% larger than its cosmic value), in contrast to a large (6 to 15) Fe overabundance derived previously. However, we confirm earlier conclusions that Si and S abundances are 2 to 3 times solar. These modest enhancements of Si, S, and Fe may be attributed either to heavy-element enriched SN ejecta or to the initial chemical abundances of the SN progenitor, which originated in the metal-rich inner Galaxy. The comparison of our models with the observed spectra confirm theoretical predictions that moderate electron heating occurs at strong collisionless shock fronts, with the implied electron/mean temperature ratio of approximately 0.5.

Borkowski, Kazimierz J.↗
NASA NTRS2020

Spectral Evidence for Heating at Large Column Mass in Umbral Solar Flare Kernels. I. IRIS Near-UV Spectra of the X1 Solar Flare of 2014 October 25

The GOES X1 flare SOL2014-10-25T17:08:00 was a three-ribbon solar flare observed with the Interface Region Imaging Spectrograph (IRIS)in the near-UV(NUV)and far-UV. One of the flare ribbons crossed a sunspot umbra, producing a dramatic, ∼1000% increase in the NUV continuum radiation. We comprehensively analyze the UV spectral data of the umbral flare brightenings, which provide new challenges for radiative−hydrodynamic modeling of the chromospheric velocity field and the white-light continuum radiation. The emission line profiles in the umbral flare brightenings exhibit redshifts and profile asymmetries, but these are significantly smaller than in another, well-studied X-class solar flare. We present a ratio of the NUV continuum intensity to the FeIIλ2814.45intensity. This continuum-to-line ratio is a new spectral diagnostic of significant heating at high column mass (log m/[gcm−2]>−2) during solar flares because the continuum and emission line radiation originate from relatively similar temperatures but moderately different optical depths. The full spectral readout of these IRIS data also allow for a comprehensive survey of the flaring NUV landscape: in addition to many lines of Fe II and Cr II, we identify a new solar flare emission line, HeIλ2829.91(as previously identified in laboratory and early-type stellar spectra). The Fermi/GBM hard X-ray data provide inputs to radiative−hydrodynamic models (which will be presented in Paper II) in order to better understand the large continuum-to-line ratios, the origin of the white-light continuum radiation, and the role of electron beam heating in the low atmosphere.

Adam F Kowalski↗
NASA NTRS1993

Ejecting basaltic achondrites from Vesta: Hydrodynamical impact models

Vesta is a large asteroid whose crust is mostly basaltic. Spectral heterogeneity suggests a sizable olivine feature which may be explained as impact excavation. The spectral data probably show a localized feature approximately 200 km in diameter or a diffuse feature approximately 400 km in diameter. Lightcurve irregularities suggest heterogeneity on a similar scale. This heterogeneity may represent the crater bowl, the extent of its ejecta deposit, or indeed something unrelated to cratering. In any case, drawing direct inferences about the state of Vesta's surface on the basis of these observations involves substantial speculation. We are still far from understanding the boundary between cratering and catastrophic disruption, particularly on targets for which strength and self-gravity both matter. But we are now able to model the underlying physical process -- dynamic fragmentation -- accurately with fragmentation hydrocodes such as SALE 2D and SPH3D. We present several impact scenarios for Vesta; our study is similar to a previous impact model for the formation of Stickney crater on Phobos. We illustrate the effects of impactors of different sizes and velocities, and the effect of gravity and structural inhomogeneity.

Asphaug, E.↗
NASA NTRS2021

Modeling the Hydrodynamics, Sediment Transport, and Valley Incision of Outlet-forming Floods from Martian Crater Lakes

Open-basin lakes that have outlets, formed where the stored water breached the basin-confining topography (e.g., crater rims), have been recognized on Mars for decades. However, the mechanics involved in formation of these outlets has not been studied in detail. Here, we study this outlet formation process through a series of numerical experiments accomplished by adding erosion and sediment transport routines to an existing hydrodynamic simulation, ANUGA. These experiments are consistent with outlets being carved in single flood events, at least if sediment was available throughout outlet formation as in a transport-limited regime. Peak discharges from these open-basin lakes, including from Jezero crater, were 10^4–10^6 m^3/s, and the outlet-forming floods lasted for days to weeks, depending on lake size. Moreover, we find that sediment transport by suspension was likely important during these floods – outlets could not have reached their observed sizes if flow was limited to bedload transport alone. This in turn may explain why small open-basin lakes are less commonly recognized on Mars, as small lakes had less inherent capacity to transport sediment in suspension and thus incise confining topography catastrophically. Geometry, topographic setting, and earlier fluvial valley development may have played a significant role in influencing the magnitude of outlet-forming floods and resulting outlet canyon incision.

Mars↗
NASA NTRS1986

A model of inner cometary ionospheres

A hydrodynamical model to describe the movement of the thermalized charged components in the inner ionosphere of comet Halley is presented. Photoelectrons are included by applying a two-stream-type approach. The numerical scheme describes shock transitions in a natural way. Solutions are obtained for a number of different assumptions concerning electron heating rates but all show that the electron temperature increases sharply where the collisional electron neutral coupling becomes unimportant. This temperature increase is accompanied by an increase in the plasma pressure and in its associated polarization electric field, and causes the plasma flow to go subsonic. In certain cases this transition occurs as an inner shock which may explain the observed ion pile-up in Halley's comet.

Korosmezey, A.↗
NASA NTRS1991

A new approach to plasmasphere refilling: Anomalous plasma effects

During the last 10 months of the grant, both laminar and anomalous plasma processes occurring during the refilling of the outer plasmasphere after magnetic storms are investigated. Theoretical investigations were based on two types of models: (1) two-stream hydrodynamic model in which plasma flows from the conjugate ionospheres are treated as separate fluids and the ion temperature anisotropies are treated self-consistently; and (2) large-scale particle-in-cell code.

Singh, N.↗
NASA NTRS2015

Hazard Assessment from Storm Tides and Rainfall on a Tidal River Estuary

Here, we report on methods and results for a model-based flood hazard assessment we have conducted for the Hudson River from New York City to Troy/Albany at the head of tide. Our recent work showed that neglecting freshwater flows leads to underestimation of peak water levels at up-river sites and neglecting stratification (typical with two-dimensional modeling) leads to underestimation all along the Hudson. As a result, we use a three-dimensional hydrodynamic model and merge streamflows and storm tides from tropical and extratropical cyclones (TCs, ETCs), as well as wet extratropical cyclone (WETC) floods (e.g. freshets, rain-on-snow events). We validate the modeled flood levels and quantify error with comparisons to 76 historical events. A Bayesian statistical method is developed for tropical cyclone streamflows using historical data and consisting in the evaluation of (1) the peak discharge and its pdf as a function of TC characteristics, and (2) the temporal trend of the hydrograph as a function of temporal evolution of the cyclone track, its intensity and the response characteristics of the specific basin. A k-nearest-neighbors method is employed to determine the hydrograph shape. Out of sample validation tests demonstrate the effectiveness of the method. Thus, the combined effects of storm surge and runoff produced by tropical cyclones hitting the New York area can be included in flood hazard assessment. Results for the upper Hudson (Albany) suggest a dominance of WETCs, for the lower Hudson (at New York Harbor) a case where ETCs are dominant for shorter return periods and TCs are more important for longer return periods (over 150 years), and for the middle-Hudson (Poughkeepsie) a mix of all three flood events types is important. However, a possible low-bias for TC flood levels is inferred from a lower importance in the assessment results, versus historical event top-20 lists, and this will be further evaluated as these preliminary methods and results are finalized. Future funded work will quantify the influences of sea level rise and flood adaptation plans (e.g. surge barriers). It would also be valuable to examine how streamflows from tropical cyclones and wet cool-season storms will change, as this factor will dominate at upriver locations.

Hazard assessment↗
NASA NTRS2021

Accuracy Assessment of Global Internal-Tide Models Using Satellite Altimetry Models Using Satellite Altimetry

Altimeter measurements are corrected for several geophysical parameters in order to access ocean signals of interest, like mesoscale or sub-mesoscale variability. The ocean tide is one of the most critical corrections due to the amplitude of the tidal elevations and to the aliasing phenomena of high-frequency signals into the lower-frequency band, but the internal-tide signatures at the ocean surface are not yet corrected globally. Internal tides can have a signature of several centimeters at the surface with wavelengths of about 50–250 km for the first mode and even smaller scales for higher-order modes. The goals of the upcoming Surface Water Ocean Topography (SWOT) mission and other high-resolution ocean measurements make the correction of these small-scale signals a challenge, as the correction of all tidal variability becomes mandatory to access accurate measurements of other oceanic signals. In this context, several scientific teams are working on the development of new internal-tide models, taking advantage of the very long altimeter time series now available, which represent an unprecedented and valuable global ocean database. The internal-tide models presented here focus on the coherent internal-tide signal and they are of three types: empirical models based upon analysis of existing altimeter missions, an assimilative model and a three-dimensional hydrodynamic model. A detailed comparison and validation of these internal-tide models is proposed using existing satellite altimeter databases. The analysis focuses on the four main tidal constituents: M2, K1, O1 and S2. The validation process is based on a statistical analysis of multi-mission altimetry including Jason-2 and Cryosphere Satellite-2 data. The results show a significant altimeter variance reduction when using internal-tide corrections in all ocean regions where internal tides are generating or propagating. A complementary spectral analysis also gives some estimation of the performance of each model as a function of wavelength and some insight into the residual non-stationary part of internal tides in the different regions of interest. This work led to the implementation of a new internal-tide correction (ZARON'one) in the next geophysical data records version-F (GDR-F) standards.

Loren Carrere↗
NASA NTRS1997

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗
NASA NTRS1990

Modeling the dust size distribution in comets with dust fragmentation

A hydrodynamic model was developed of a spherically symmetric dusty gas flow in a cometary atmosphere assuming a single fluid, inviscid, perfect gas. The hydrodynamics for gas and dust, which involves the gas drag force (momentum transfer), heat exchange between gas and dust, photodissociation energy for H2O gas, and radiative heating and cooling terms for dust particles, are solved using the Gear method for stiff, coupled differential equations. Calculations were done with a dust size distribution for radii alpha = 0.01 micron to 10 cm with densities variable with the size. A nucleus size of 4.0 km radius with a density of 0.5 g/cu cm and a total dust-to-gas mass ratio chi = 1 was adopted. There are indications from in situ observations that dust particle fragment into smaller ones. Fragmentation of dust particles was incorporated into the model. This is done by adding source and sink terms in the continuity equations for the dust. Lifetimes for the decay of dust particles were assumed as a function of particle size. It is also assumed that dust particles always fragment only into the next smaller size.

Konno, Ichishiro↗
NASA NTRS1980

The chromosphere and transition region

The physical processes occurring as a result of the transfer of energy and momentum from the primary solar flare energy release site in the corona to the underlying chromosphere and transition region during the course of the flare are investigated through a comparison of theoretical models and observational data. Static, dynamic and hydrodynamic models of the lower-temperature chromospheric flare are reviewed. The roles of thermal conduction, radiation, fast particles and mass motion in chromosphere-corona interactions are analyzed on the basis of Skylab UV, EUV and X-ray data, and empirical and synthetic models of the chromospheric and upper photospheric responses to flares are developed. The canonical model of chromospheric heating during flares as a result of primary energy release elsewhere is found to be justified in the chromosphere as a whole, although not entirely as the temperature minimum, and a simplified model of horizontal chromospheric flare structure based on results obtained is presented.

Canfield, R. C.↗
NASA NTRS1987

SIR-B ocean-wave enhancement with fast Fourier transform techniques

Shuttle Imaging Radar (SIR-B) imagery is Fourier filtered to remove the estimated system-transfer function, reduce speckle noise, and produce ocean scenes with a gray scale that is proportional to wave height. The SIR-B system response to speckled scenes of uniform surfaces yields an estimate of the stationary wavenumber response of the imaging radar, modeled by the 15 even terms of an eighth-order two-dimensional polynomial. Speckle can also be used to estimate the dynamic wavenumber response of the system due to surface motion during the aperture synthesis period, modeled with a single adaptive parameter describing an exponential correlation along track. A Fourier filter can then be devised to correct for the wavenumber response of the remote sensor and scene correlation, with subsequent subtraction of an estimate of the speckle noise component. A linearized velocity bunching model, combined with a surface tilt and hydrodynamic model, is incorporated in the Fourier filter to derive estimates of wave height from the radar intensities corresponding to individual picture elements.

Tilley, David G.↗
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↗
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.↗
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.↗
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.↗