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 109 records · Page 6

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.↗

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.↗

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.↗

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↗

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↗

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.↗

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.↗

Core plasma in the magnetosphere

The paper presents a review of new findings related to core plasma (which includes ions and electrons with energies less than 50 eV) identified in studies of magnetosphere during 1983-1986. These findings include the cleft ion fountain; the plasmaspheric heavy ion torus; core molecular ions; high core-plasma densities in the plasma-sheet boundary layer; intense transverse heating of core plasmas in the equatorial regions; the supersonic polar wind (H/+/ and He/+/); toroidal or ring heavy-ion distribution functions in the auroral region; conic-to-field-aligned ion signatures of inverted-V auroral events; and the large-scale electron density structure of the magnetosphere. The progress made on modeling core plasma includes simulation of thermal helium heating at the equator; kinetic models of the plasma transport from localized ionospheric regions; time-dependent hydrodynamic models of the polar wind and plasmasphere-ionosphere coupling; and kinetic models of plasmasphere refilling.

Horwitz, J. L.↗

Understanding Heating in Active Region Cores through Machine Learning. I. Numerical Modeling and Predicted Observables

To adequately constrain the frequency of energy deposition in active region cores in the solar corona, systematic comparisons between detailed models and observational data are needed. In this paper, we describe a pipeline for forward modeling active region emission using magnetic field extrapolations and field-aligned hydrodynamic models. We use this pipeline to predict time-dependent emission from active region NOAA 1158 for low-, intermediate-, and high-frequency nanoflares. In each pixel of our predicted multi-wavelength, time-dependent images, we compute two commonly used diagnostics: the emission measure slope and the time lag. We find that signatures of the heating frequency persist in both of these diagnostics. In particular, our results show that the distribution of emission measure slopes narrows and the mean decreases with decreasing heating frequency and that the range of emission measure slopes is consistent with past observational and modeling work. Furthermore, we find that the time lag becomes increasingly spatially coherent with decreasing heating frequency while the distribution of time lags across the whole active region becomes more broad with increasing heating frequency. In a follow-up paper, we train a random forest classifier on these predicted diagnostics and use this model to classify real observations of NOAA 1158 in terms of the underlying heating frequency.

UV radiation↗

Complete Galilean-Invariant Lattice BGK Models for the Navier-Stokes Equation

Galilean invariance has been an important issue in lattice-based hydrodynamics models. Previous models concentrated on the nonlinear advection term. In this paper, we take into account the nonlinear response effect in a systematic way. Using the Chapman-Enskog expansion up to second order, complete Galilean invariant lattice BGK models in one dimension (theta = 3) and two dimensions (theta = 1) for the Navier-Stokes equation have been obtained.

Qian, Yue-Hong↗

The application of single particle hydrodynamics in continuum models of multiphase flow

A review of the application of single particle hydrodynamics in models for the exchange of interphase momentum in continuum models of multiphase flow is presented. Considered are the equations of motion for a laminar, mechanical two phase flow. Inherent to this theory is a model for the interphase exchange of momentum due to drag between the dispersed particulate and continuous fluid phases. In addition, applications of two phase flow theory to de-mixing flows require the modeling of interphase momentum exchange due to lift forces. The applications of single particle analysis in deriving models for drag and lift are examined.

Decker, Rand↗

Chromospheric heating by acoustic shocks - A confrontation of GHRS observations of Alpha Tauri (K5 III) with ab initio calculations

We compare ab initio calculations of semiforbidden C II line profiles near 2325 A with recently published observations of the inactive red giant Alpha Tau (K5 III) obtained using the GHRS on board the Hubble Space Telescope. Our one-dimensional, time-dependent calculations assume that the chromosphere is heated by stochastic acoustic shocks generated by photospheric convection. We calculate various models using results from traditional (mixing length) convection zone calculations as input to hydrodynamical models. The semiforbidden C II line profiles and ratios provide sensitive diagnostics of chromospheric velocity fields, electron densities, and temperatures. We identify major differences between observed and computed line profiles which are related to basic gas dynamics and which are probably not due to technical modeling restrictions. If the GHRS observations are representative of chromospheric conditions at all epochs, then one (or more) of our model assumptions must be incorrect. Several possibilities are examined. We predict time variability of semiforbidden C II lines for comparison with observations. Based upon data from the IUE archives, we argue that photospheric motions associated with supergranulation or global pulsation modes are unimportant in heating the chromosphere of Alpha Tau.

Judge, P. G.↗

Orbital, Rotational and Climatic Interactions: Energy Dissipation and Angular Momentum Exchange in the Earth-Moon System

A numerical ocean tide model has been developed and tested using highly accurate TOPEX/Poseidon (T/P) tidal solutions. The hydrodynamic model is based on time stepping a finite difference approximation to the non-linear shallow water equations. Two novel features of our implementation are a rigorous treatment of self attraction and loading (SAL), and a physically based parameterization for internal tide (IT) radiation drag. The model was run for a range of grid resolutions, and with variations in model parameters and bathymetry. For a rational treatment of SAL and IT drag, the model run at high resolution (1/12 degree) fits the T/P solutions to within 5 cm RMS in the open ocean. Both the rigorous SAL treatment and the IT drag parameterization are required to obtain solutions of this quality. The sensitivity of the solution to perturbations in bathymetry suggest that the fit to T/P is probably now limited by errors in this critical input. Since the model is not constrained by any data, we can test the effect of dropping sea-level to match estimated bathymetry from the last glacial maximum (LGM). Our results suggest that the 100 m drop in sea-level in the LGM would have significantly increased tidal amplitudes in the North Atlantic, and increased overall tidal dissipation by about 40%. However, details in tidal solutions for the past 20 ka are sensitive to the assumed stratification. IT drag accounts for a significant fraction of dissipation, especially in the LGM when large areas of present day shallow sea were exposed, and this parameter is poorly constrained at present.

Egbert, Gary D.↗

The shock process and light element production in supernovae envelopes

Detailed hydrodynamic modeling of the passage of supernova shocks through the hydrogen envelopes of blue and red progenitor stars was carried out to explore the sensitivity to model conditions of light element production (specifically Li-7 and B-11) which was noted by Dearborn, Schramm, Steigman and Truran (1989) (DSST). It is found that, for stellar models with M is less than or approximately 100 M solar mass, current state of the art supernova shocks do not produce significant light element yields by hydrodynamic processes alone. The dependence of this conclusion on stellar models and on shock strengths is explored. Preliminary implications for Galactic evolution of lithium are discussed, and it is suspected that intermediate mass red giant stars may be the most consistent production site for lithium.

Brown, Lawrence E.↗

The shock process and light-element production in supernova envelopes

Detailed hydrodynamic modeling of the passage of supernova shocks through the hydrogen enevlopes of blue and red progenitor stars was carried out to explore the sensitivity to model conditions of light element production (specifically Li7 and B-11) which was noted by Dearborn, Schramm, Steigman and Truran (1989) (DSST). It is found that, for stellar models with M is less than or approximately 100 M solar mass, current state of the art supernova shocks do not produce significant light element yields by hydrodynamic processes alone. The dependence of this conclusion on stellar models and on shock strengths is explored. Preliminary implications for Galactic evolution of lithium are discussed, and it is suspected that intermediate mass red giant stars may be the most consistent production site for lithium.

Brown, Lawrence E.↗

The Influence of Dynamic Contact Angle on Wetting Dynamics

When surface tension forces dominate, and regardless of whether the situation is static or dynamic, the contact angle (the angle the interface between two immiscible fluids makes when it contacts a solid) is the key parameter that determines the shape of a fluid-fluid interface. The static contact angle is easy to measure and implement in models predicting static capillary surface shapes and such associated quantities as pressure drops. By contrast, when the interface moves relative to the solid (as in dynamic wetting processes) the dynamic contact angle is not identified unambiguously because it depends on the geometry of the system Consequently, its determination becomes problematic and measurements in one geometry cannot be applied in another for prediction purposes. However, knowing how to measure and use the dynamic contact angle is crucial to determine such dynamics as a microsystem throughput reliably. In this talk we will present experimental and analytical efforts aimed at resolving modeling issues present in dynamic wetting. We will review experiments that show the inadequacy of the usual hydrodynamic model when a fluid-fluid meniscus moves over a solid surface such as the wall of a small tube or duct. We will then present analytical results that show how to parametrize these problems in a predictive manner. We will illustrate these ideas by showing how to implement the method in numerical fluid mechanical calculations.

Rame, Enrique↗

The origin of planets orbiting millisecond pulsars

A model for the formation of planets around millisecond pulsar which no longer have stellar companions is suggested. Detailed hydrodynamical models are presented which suggest that planet formation can occur either in a low-mass X-ray binary progenitor to a progenitor of a star-vaporizing millisecond pulsar when the neutron star is accreting material driven off its companion by X-ray irradiation or after a pulsar has formed and is vaporizing its companion. In both cases a circumbinary disk is created in which planets can form on a timescale of 10 exp 5 to 10 exp 6 yrs and the planets can survive a second phase in which the companion star moves toward the pulsar and is completely vaporized.

Tavani, Marco↗