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 73 records · Page 4

Calibration of Two-dimensional Floodplain Modeling in the Atchafalaya River Basin Using SAR Interferometry

Two-dimensional (2D) satellite imagery has been increasingly employed to improve prediction of floodplain inundation models. However, most focus has been on validation of inundation extent, with little attention on the 2D spatial variations of water elevation and slope. The availability of high resolution Interferometric Synthetic Aperture Radar (InSAR) imagery offers unprecedented opportunity for quantitative validation of surface water heights and slopes derived from 2D hydrodynamic models. In this study, the LISFLOOD-ACC hydrodynamic model is applied to the central Atchafalaya River Basin, Louisiana, during high flows typical of spring floods in the Mississippi Delta region, for the purpose of demonstrating the utility of InSAR in coupled 1D/2D model calibration. Two calibration schemes focusing on Manning s roughness are compared. First, the model is calibrated in terms of water elevations at a single in situ gage during a 62 day simulation period from 1 April 2008 to 1 June 2008. Second, the model is calibrated in terms of water elevation changes calculated from ALOS PALSAR interferometry during 46 days of the image acquisition interval from 16 April 2008 to 1 June 2009. The best-fit models show that the mean absolute errors are 3.8 cm for a single in situ gage calibration and 5.7 cm/46 days for InSAR water level calibration. The optimum values of Manning's roughness coefficients are 0.024/0.10 for the channel/floodplain, respectively, using a single in situ gage, and 0.028/0.10 for channel/floodplain the using SAR. Based on the calibrated water elevation changes, daily storage changes within the size of approx 230 sq km of the model area are also calculated to be of the order of 107 cubic m/day during high water of the modeled period. This study demonstrates the feasibility of SAR interferometry to support 2D hydrodynamic model calibration and as a tool for improved understanding of complex floodplain hydrodynamics

Hahn Chul Jung↗

Non-Parabolic Hydrodynamic Formulations for the Simulation of Inhomogeneous Semiconductor Devices

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models cannot fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations or the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship ((hk)(exp 2)/2m = W(1 + alphaW). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(exp y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships; parabolic, Kane dispersion and power law dispersion.

Smith, A. W.↗

Non-parabolic hydrodynamic formulations for the simulation of inhomogeneous semiconductor devices

Hydrodynamic models are becoming prevalent design tools for small scale devices and other devices in which high energy effects can dominate transport. Most current hydrodynamic models use a parabolic band approximation to obtain fairly simple conservation equations. Interest in accounting for band structure effects in hydrodynamic device simulation has begun to grow since parabolic models can not fully describe the transport in state of the art devices due to the distribution populating non-parabolic states within the band. This paper presents two different non-parabolic formulations of the hydrodynamic model suitable for the simulation of inhomogeneous semiconductor devices. The first formulation uses the Kane dispersion relationship (hk)(exp 2)/2m = W(1 + alpha(W)). The second formulation makes use of a power law ((hk)(exp 2)/2m = xW(sup y)) for the dispersion relation. Hydrodynamic models which use the first formulation rely on the binomial expansion to obtain moment equations with closed form coefficients. This limits the energy range over which the model is valid. The power law formulation readily produces closed form coefficients similar to those obtained using the parabolic band approximation. However, the fitting parameters (x,y) are only valid over a limited energy range. The physical significance of the band non-parabolicity is discussed as well as the advantages/disadvantages and approximations of the two non-parabolic models. A companion paper describes device simulations based on the three dispersion relationships: parabolic, Kane dispersion, and power low dispersion.

Smith, Arlynn W.↗

Energy Models for One-Carrier Transport in Semiconductor Devices

Moment models of carrier transport, derived from the Boltzmann equation, made possible the simulation of certain key effects through such realistic assumptions as energy dependent mobility functions. This type of global dependence permits the observation of velocity overshoot in the vicinity of device junctions, not discerned via classical drift-diffusion models, which are primarily local in nature. It was found that a critical role is played in the hydrodynamic model by the heat conduction term. When ignored, the overshoot is inappropriately damped. When the standard choice of the Wiedemann-Franz law is made for the conductivity, spurious overshoot is observed. Agreement with Monte-Carlo simulation in this regime required empirical modification of this law, or nonstandard choices. Simulations of the hydrodynamic model in one and two dimensions, as well as simulations of a newly developed energy model, the RT model, are presented. The RT model, intermediate between the hydrodynamic and drift-diffusion model, was developed to eliminate the parabolic energy band and Maxwellian distribution assumptions, and to reduce the spurious overshoot with physically consistent assumptions. The algorithms employed for both models are the essentially non-oscillatory shock capturing algorithms. Some mathematical results are presented and contrasted with the highly developed state of the drift-diffusion model.

Jerome, Joseph W.↗

Enthalpy-Based Thermal Evolution of Loops: III. Comparison of Zero-Dimensional Models

Zero dimensional (0D) hydrodynamic models, provide a simple and quick way to study the thermal evolution of coronal loops subjected to time-dependent heating. This paper presents a comparison of a number of 0D models that have been published in the past and is intended to provide a guide for those interested in either using the old models or developing new ones. The principal difference between the models is the way the exchange of mass and energy between corona, transition region and chromosphere is treated, as plasma cycles into and out of a loop during a heating-cooling cycle. It is shown that models based on the principles of mass and energy conservation can give satisfactory results at some, or, in the case of the Enthalpy Based Thermal Evolution of Loops (EBTEL) model, all stages of the loop evolution. Empirical models can lead to low coronal densities, spurious delays between the peak density and temperature, and, for short heating pulses, overly short loop lifetimes.

Cargill, P. J.↗

Improved Calibration of Modeled Discharge and Storage Change in the Atchafalaya Floodplain Using SAR Interferometry

This study focuses on the feasibility of using SAR interferometry to support 2D hydrodynamic model calibration and provide water storage change in the floodplain. Two-dimensional (2D) flood inundation modeling has been widely studied using storage cell approaches with the availability of high resolution, remotely sensed floodplain topography. The development of coupled 1D/2D flood modeling has shown improved calculation of 2D floodplain inundation as well as channel water elevation. Most floodplain model results have been validated using remote sensing methods for inundation extent. However, few studies show the quantitative validation of spatial variations in floodplain water elevations in the 2D modeling since most of the gauges are located along main river channels and traditional single track satellite altimetry over the floodplain are limited. Synthetic Aperture Radar (SAR) interferometry recently has been proven to be useful for measuring centimeter-scale water elevation changes over the floodplain. In the current study, we apply the LISFLOOD hydrodynamic model to the central Atchafalaya River Basin, Louisiana, during a 62 day period from 1 April to 1 June 2008 using two different calibration schemes for Manning's n. First, the model is calibrated in terms of water elevations from a single in situ gauge that represents a more traditional approach. Due to the gauge location in the channel, the calibration shows more sensitivity to channel roughness relative to floodplain roughness. Second, the model is calibrated in terms of water elevation changes calculated from ALOS PALSAR interferometry during 46 days of the image acquisition interval from 16 April 2008 to 1 June 2009. Since SAR interferometry receives strongly scatters in floodplain due to double bounce effect as compared to specular scattering of open water, the calibration shows more dependency to floodplain roughness. An iterative approach is used to determine the best-fit Manning's n for the two different calibration approaches. Results suggest similar floodplain roughness but slightly different channel roughness. However, application of SAR interferometry provides a unique view of the floodplain flow gradients, not possible with a single gauge calibration. These gradients, allow improved computation of water storage change over the 46-day simulation period. Overall, the results suggest that the use of 2D SAR water elevation changes in the Atchafalaya basin offers improved understanding and modeling of floodplain hydrodynamics.

Jung, Hahn Chul↗

Integrated Analysis of Interferometric SAR, Satellite Altimetry and Hydraulic Modeling to Quantify Louisiana Wetland Dynamics

Wetland loss in Louisiana has been accelerating due primarily to anthropogenic and nature processes, and is being advocated as a problem with national importance. Accurate measurement or modeling of wetland-wide water level changes, its varying extent, its storage and discharge changes resulting in part from sediment loads, erosion and subsidence are fundamental to assessment of hurricane-induced flood hazards and wetland ecology. Here, we use innovative method to integrate interferometric SAR (InSAR) and satellite radar altimetry for measuring absolute or geocentric water level changes and applied the methodology to remote areas of swamp forest in coastal Louisiana. Coherence analysis of InSAR pairs suggested that the HH polarization is preferred for this type of observation, and polarimetric analysis can help to identi:fy double-bonnce backscattering areas in the wetland. Envisat radar altimeter-measured 18- Hz (along-track sampling of 417 m) water level data processed with regional stackfile method have been used to provide vertical references for water bodies separated by levees. The high-resolution (approx.40 m) relative water changes measured from ALOS PALSAR L-band and Radarsat-l C-band InSAR are then integrated with Envisat radar altimetry to obtain absolute water level. The resulting water level time series were validated with in situ gauge observations within the swamp forest. Furthermore, we compare our water elevation changes with 2D flood modeling from LISFLOOD hydrodynamic model. Our study demonstrates that this new technique allows retrospective reconstruction and concurrent monitoring of water conditions and flow dynamics in wetlands, especially those lacking gauge networks.

Lee, Hyongki↗

Properties of large scale plasma flow during the early stage of the plasmaspheric refilling

The objective is to better characterize the macroscopic properties of the interhemisphere plasma flow by solving a more complete set of hydrodynamic equations than that solved previously. Specifically, the ion continuity, momentum and energy equations were solved for the plasma flow along the closed magnetic field lines. During the initial stage of the supersonic outflow in the equatorial region, the ions cool substantially. Using the hydrodynamic model for the large-scale plasma flow, the dynamics of shocks was examined which form in the geomagnetic flux tubes during the early stages of refilling. These shocks are more like those forming in neutral gases than the electrostatic shocks driven by microinstabilities involving ion-ion interaction. Therefore, the shocks seen in the hydrodynamic model are termed as hydrodynamic shocks. Such shocks are generally unsteady and therefore the usual shock jump conditions given by Rankine-Hugoniot relations are not strictly applicable to them. The density, flow velocity and temperature structures associated with the shocks are examined for both asymmetrical and symmetrical flows. In the asymmetrical flow the outflow from one of two conjugate ionospheres is dominant. On the other hand, in the symmetrical case outflows from the two ionospheric sources are identical. Both cases are treated by a two-stream model. In the late type of flow, the early-time refilling shows a relaxation type of oscillation, which is driven by the large-scale interactions between the two identical streams. After this early stage, the resulting temperature structure shows some interesting features. In the equatorial region the streams are isothermal, but in the off-equatorial regions the streams have quite different temperatures, and also densities and flow velocities. The dense and slow stream is found to be warmer than the low-density fast stream. In the late stage of refilling, the temperature is found to steadily increase from the conjugate ionospheres towards the equator; the equatorial temperature is found to be as high as about 8000 K compared to the ionospheric temperature of 3600 K.

Singh, Nagendra↗

A modified thermal conductivity for low density plasma magnetic flux tubes

In response to inconsistencies which have arisen in results from a hydrodynamic model in simulation of high ion temperature (1-2 eV) observed in low density, outer plasmasphere flux tubes, we postulate a reduced thermal conductivity coefficient in which only particles in the loss cone of the quasi-collisionless plasma contribute to the thermal conduction. Other particles are assumed to magnetically mirror before they reach the topside ionosphere and therefore not to remove thermal energy from the plasmasphere. This concept is used to formulate a mathematically simple, but physically limiting model for a modified thermal conductivity coefficient. When this modified coefficient is employed in the hydrodynamic model in a case study, the inconsistencies between simulation results and observations are largely resolved. The high simulated ion temperatures are achieved with significantly lower ion temperatures in the topside ionosphere. We suggest that this mechanism may be operative under the limited low density, refilling conditions in which high ion temperatures are observed.

Comfort, R. H.↗

Non-LTE model calculations for SN 1987A and the extragalactic distance scale

This paper presents model atmospheres for the first week of SN 1987A, based on the luminosity and density/velocity structure from hydrodynamic models of Woosley (1988). The models account for line blanketing, expansion, sphericity, and departures from LTE in hydrogen and helium and differ from previously published efforts because they represent ab initio calculations, i.e., they contain essentially no free parameters. The formation of the UV spectrum is dominated by the effects of line blanketing. In the absorption troughs, the Balmer line profiles were fit well by these models, but the observed emissions are significantly stronger than predicted, perhaps due to clumping. The generally good agreement between the present synthetic spectra and observations provides independent support for the overall accuracy of the hydrodynamic models of Woosley. The question of the accuracy of the Baade-Wesselink method is addressed in a detailed discussion of its approximations. While the application of the standard method produces a distance within an uncertainty of 20 percent in the case of SN 1987A, systematic errors up to a factor of 2 are possible, particularly if the precursor was a red supergiant.

Schmutz, W.↗

Gas flow in the barred spiral galaxy NGC 5383

Line-of-sight gas velocities calculated from a numerical hydrodynamical model of a barred spiral galaxy are compared with the observations and models made by Peterson, Rubin, Ford, and Thonnard of the gas flow in NGC 5383. The hydrodynamical model provides a somewhat better fit to the observations of NGC 5383 than does their symmetric, warped-disk model. It is argued that the basic model appropriate to understanding the gas motions in barred spiral galaxies is not circular rotation but rather a characteristic elliptical pattern of flow.

Huntley, J. M.↗

Development of an in-vitro circulatory system with known resistance and capacitance

An in-vitro (hydrodynamic) model of the circulatory system was developed. The model consisted of a pump, compliant tubing, and valves for resistance. The model is used to simulate aortic pressure and flow. These parameters were measured using a Konigsburg Pressure transducer and a Triton ART2 flow probe. In addition, venous pressure and flow were measured on the downstream side of the resistance. The system has a known compliance and resistance. Steady and pulsatile flow tests were conducted to determine the resistance of the model. A static compliance test was used to determine the compliance of the system. The aortic pressure and flow obtained from the hydrodynamic model will be used to test the accuracy of parameter estimation models such as the 2-element and 4-element Windkessel models and the 3-element Westkessel model. Verifying analytical models used in determining total peripheral resistance (TPR) and systemic arterial compliance (SAC) is important because it provides insight into hemodynamic parameters that indicate baroreceptor responsiveness to situations such as changes in gravitational acceleration.

NASA Discipline Cardiopulmonary↗

Water resources planning for rivers draining into mobile bay

A hydrodynamic model describing water movement and tidal elevation is formulated, computed, and used to provide basic data about water quality in natural systems. The hydrodynamic model is based on two-dimensional, unsteady flow equations. The water mass is considered to be reasonably mixed such that integration (averaging) in the depth direction is a valid restriction. Convective acceleration, the Coriolis force, wind and bottom interactions are included as contributing terms in the momentum equations. The solution of the equations is applied to Mobile Bay, and used to investigate the influence that river discharge rate, wind direction and speed, and tidal condition have on water circulation and holdup within the bay. Storm surge conditions, oil spill transport, artificial island construction, dredging, and areas subject to flooding are other topics which could be investigated using the mathematical modeling approach.

Ng, S.↗

Modelling time variable and total eclipses of the millisecond pulsar PSR 1744-24A

The 11 ms pulsar PSR 1744-24A in the globular clusters Terzan 5 is the first and only known example of a binary eclipsing pulsar showing time-variable eclipses with irregular and nonsymmetric eclipse patterns. Occasionally, the pulsar eclipse becomes total, and no radio emission is detected for an extended period of time. We present here numerical results of a hydrodynamical model aimed at explaining the remarkable eclipse properties of PSR 1744-24A. We show that the hydrodynamic time scale for the partial or total pulsar enshrouding is of order of the orbital period for a relatively large mass loss rate, in agreement with observations. The time-variable eclipse egress positions imply that the mass loss rate from the companion varies within one orbital period. Hydrodynamic modelling of the mass outflow of PSR 1744-24A is of crucial importance for the understanding of a class of millisecond pulsars whose radio emission might be occasionally or permanently quenched by mass outflows from their companions. The properties of PSR 1744-24A suggest a new way of studying millisecond pulsars involving radio, X-ray, gamma-ray, and optical observations.

Tavani, Marco↗

Evaluating the Impact of Land Use Change on Submerged Aquatic Vegetation Stressors in Mobile Bay

Alabama coastal systems have been subjected to increasing pressure from a variety of activities including urban and rural development, shoreline modifications, industrial activities, and dredging of shipping and navigation channels. The impacts on coastal ecosystems are often observed through the use of indicator species. One such indicator species for aquatic ecosystem health is submerged aquatic vegetation (SAV). Watershed and hydrodynamic modeling has been performed to evaluate the impact of land use change in Mobile and Baldwin counties on SAV stressors and controlling factors (temperature, salinity, and sediment) in Mobile Bay. Watershed modeling using the Loading Simulation Package in C++ (LSPC) was performed for all watersheds contiguous to Mobile Bay for land use scenarios in 1948, 1992, 2001, and 2030. Landsat-derived National Land Cover Data (NLCD) were used in the 1992 and 2001 simulations after having been reclassified to a common classification scheme. The Prescott Spatial Growth Model was used to project the 2030 land use scenario based on current trends. The LSPC model simulations provided output on changes in flow, temperature, and sediment for 22 discharge points into the Bay. Theses results were inputted in the Environmental Fluid Dynamics Computer Code (EFDC) hydrodynamic model to generate data on changes in temperature, salinity, and sediment on a grid with four vertical profiles throughout Mobile Bay. The changes in the aquatic ecosystem were used to perform an ecological analysis to evaluate the impact on SAV habitat suitability. This is the key product benefiting the Mobile Bay coastal environmental managers that integrates the influences of temperature, salinity, and sediment due to land use driven flow changes with the restoration potential of SAVs.

Al-Hamdan, Mohammad↗

Multi-wavelength Observations and Modelling of a Solar Flare

Aims: The aim of this work is to investigate the dynamic behavior of a C-class solar flare through the evolution of temperature, emission measure, energy loss and velocity. In particular, the variation of these properties with time are studied using multi-wavelength observations in combination with a recently developed 0-D hydrodynamic model. Methods: The temperature and emission measure evolution were studied using several instruments covering a wide range of temperatures - the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI, >5 MK), GOES-12 (5- 30 MK), the Transition Region and Coronal Explorer (TRACE 171 A, 1 MK) and the Coronal Diagnostic Spectrometer (CDS, 0.03-8 MK). The temperature and emission measure were analysed through the systematic cooling of flare plasma through the response functions of these instruments. These parameters were then investigated using the Enthalpy Based Thermal Evolution of Loops model (EBTEL). The Doppler shifts at both flare footpoints were analysed using five emission lines seen by CDS. Results: The flare began with clear evidence for pre-flare heating. Upflows of approx.90 km/s and low level emission, both observed in Fe XIX before the main impulsive phase were explained by pre-flare gentle chromospheric evaporation. During the main impulsive phase, the flare plasma was heated to a temperature of >13 MK in approximately 10 minutes. Explosive chromospheric evaporation was observed, driving upflows of approx.80 km/s in Fe XIX and simultaneous downflows of approx.20 km/s in He I and O v. At the peak of the Rare, conduction modelled by EBTEL was found to be the dominant loss mechanism, working efficiently to both lower the temperatures and drive gentle chromospheric evaporation. As the temperature fell below approx.8 MK, radiation became the dominant loss mechanism. During the final stages of the decay phase, downflowing plasma was observed at the footpoints in He I, O v and Mg x at velocities of up to approx.40 km/s, suggesting loop draining occurred. Conclusions. This is the first extensive study of the evolution of flare plasma using both spectroscopic and broad-band instruments in conjunction with a comprehensive hydrodynamic model. The flare began with pre-flare heating and then evolved following the predictions of the standard flare model. Detailed analysis of the plasma heating mechanisms was carried out and the heating function most consistent with observations was found to be Gaussian in shape. The simulations suggested that both direct heating and heating by a non-thermal beam played significant roles in this event.

Raftery, Claire L.↗

Development of the Packed Bed Reactor ISS Flight Experiment

Packed bed reactors are compact, require minimum power and maintenance to operate, and are highly reliable. These features make this technology a leading candidate as a potential unit operation in support of long duration human space exploration. On earth, this type of reactor accounts for approximately 80% of all the reactors used in the chemical process industry today. Development of this technology for space exploration is truly crosscutting with many other potential applications (e.g., in-situ chemical processing of planetary materials and transport of nutrients through soil). NASA is developing an ISS experiment to address this technology with particular focus on water reclamation and air revitalization. Earlier research and development efforts funded by NASA have resulted in two hydrodynamic models which require validation with appropriate instrumentation in an extended microgravity environment. The first model developed by Motil et al., (2003) is based on a modified Ergun equation. This model was demonstrated at moderate gas and liquid flow rates, but extension to the lower flow rates expected in many advanced life support systems must be validated. The other model, developed by Guo et al., (2004) is based on Darcy s (1856) law for two-phase flow. This model has been validated for a narrow range of flow parameters indirectly (without full instrumentation) and included test points where the flow was not fully developed. The flight experiment presented will be designed with removable test sections to test the hydrodynamic models. The experiment will provide flexibility to test additional beds with different types of packing in the future. One initial test bed is based on the VRA (Volatile Removal Assembly), a packed bed reactor currently on ISS whose behavior in micro-gravity is not fully understood. Improving the performance of this system through an accurate model will increase our ability to purify water in the space environment.

Patton, Martin O.↗