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 469 records · Page 26

A granular flow model for dense planetary rings

In the present study of the viscosity of a differentially rotating particle disk, in the limiting case where the particles are densely packed and their collective behavior resembles that of a liquid, the pressure tensor is derived from both the equations of hydrodynamics and a simple kinetic model of collisions due to Haff (1983). Density waves and narrow circular rings are unstable if the liquid approximation applies, and the consequent nonlinear perturbations may generate 'splashing' of the ring material in the vertical direction. These results are pertinent to the origin of the ellipticities of ringlets, the nonaxisymmetric features near the outer edge of the Saturn B ring, and unexplained residuals in kinematic models of the Saturn and Uranus rings.

Borderies, N.↗

A three-dimensional model of co-rotating streams in the solar wind. 2: Hydrodynamic streams

Theoretical aspects of corotating solar wind dynamics on a global scale are explored by means of numerical simulations executed with a nonlinear, inviscid, adiabatic, single-fluid, three-dimensional (3-D) hydrodynamic formulation. A simple, hypothetical 3-D stream structure is defined on a source surface located at 35 solar radius and carefully documents its evolution to 1 AU under the influence of solar rotation. By manipulating the structure of this prototype configuration at the source surface, it is possible to elucidate the factors most strongly affecting stream evolution: (1) the intrinsic correlations among density, temperature, and velocity existing near the source; (2) the amplitude of the stream; (3) the longitudinal breadth of the stream; (4) the latitudinal breadth of the stream; and (5) the heliographic latitude of the centroid of the stream.

Pizzo, V. J.↗

CNO abundances and hydrodynamic studies of the Nova outburst. V - 1.00-solar-mass models with small mass envelopes

The paper reports on an investigation into the consequences of thermonuclear runaways in accreted hydrogen envelopes of 100 millionths of a solar mass on 1-solar-mass white dwarfs. These evolutionary sequences predict that from 10 to 50 millionths of a solar mass will be ejected with speeds from 300 to 3800 km/s (kinetic energies of 10 to the 44th-45th power ergs). Absolute visual magnitudes as high as -8.1 are attained, well within the observed range for fast novae. In addition, the shapes of the theoretical light curves are more reminiscent of an observed fast-nova light curve than those in earlier studies. The ejected material is strongly enhanced in the products of incomplete CNO burning; the most abundant of the ejected nuclei is C-13, followed by N-14 and C-12. The differences from previous studies are attributable to the lower peak temperatures reached in these sequences. These models also produce a large overabundance of Li-7, suggesting that novae may represent significant contributors to the galactic enrichment of this nucleus.

Starrfield, S.↗

Numerical models for Halley dust environments

The algebra, parameter values, and results of a reference model for the dust environment of Comet P/Halley are presented. The model employs simplified treatments of the icy conglomerate nucleus and of the hydrodynamic outflow of gas and dust. The model includes a gas production rate derived from observations of Halley's visual magnitude in 1910, gas and dust flux nonuniformly distributed over the nuclear surface, a dust distribution function consistent with large particle inferences from comet antitails, dust density variable with mass, and an approximation to the distribution of particles near the dust envelope. Many parameter values are derived by analogy with the results for other comets. The features are incorporated in a computer program which calculates the concentration, flux, and fluence of dust particles as functions of time, position and particle mass.

Divine, N.↗

Tidal disruption of Periodic Comet Shoemaker-Levy 9 and a constraint on its mean density

The apparent tidal disruption of Periodic Comet Shoemaker-Levy 9 (1993e) during a close encounter within approximately 1.62 planetary radii of Jupiter can be used along with theoretical models of tidal disruption to place an upper bound on the density of the predisruption body. Depending on the theoretical model used, these upper bounds range from rho(sub c) less than 0.702 +/- 0.080 g/cu cm for a simple analytical model calibrated by numerical smoothed particle hydrodynamics (SPH) simulations to rho(sub c) less than 1.50 +/- 0.17 g/cu cm for a detailed semianalytical model. The quoted uncertainties stem from an assumed uncertainty in the perijove radius. However, the uncertainty introduced by the different theoretical models is the major source of error; this uncertainty could be eliminated by future SPH simulations specialized to cometary disruptions, including the effects of initially prolate, spinning comets. If the SPH-based upper bound turns out to be most appropriate, it would be consistent with the predisruption body being a comet with a relatively low density and porous structure, as has been asserted previously based on observations of cometary outgassing. Regardless of which upper bound is preferable, the models all agree that the predisruption body could not have been a relatively high-density body, such as an asteroid with rho approximately = 2 g/cu cm.

Boss, Alan P.↗

Modeling the Compression of Merged Compact Toroids by Multiple Plasma Jets

A fusion propulsion scheme has been proposed that makes use of the merging of a spherical distribution of plasma jets to dynamically form a gaseous liner. The gaseous liner is used to implode a magnetized target to produce the fusion reaction in a standoff manner. In this paper, the merging of the plasma jets to form the gaseous liner is investigated numerically. The Los Alamos SPHINX code, based on the smoothed particle hydrodynamics method is used to model the interaction of the jets. 2-D and 3-D simulations have been performed to study the characteristics of the resulting flow when these jets collide. The results show that the jets merge to form a plasma liner that converge radially which may be used to compress the central plasma to fusion conditions. Details of the computational model and the SPH numerical methods will be presented together with the numerical results.

Thio, Y. C. Francis↗

Similitude in hydrodynamic tests involving planing

The problems of using models in planing tests are addressed. If one passes from the model to a hull of linear dimensions n times greater, the speeds are connected by the law of mechanical similitude. The normal forces given by the hydrodynamic equations (perfect fluid) also follow the law of dynamic similitude (Reech's method) and are multiplied by n(exp 3). A series of tests were performed and the actual results were compared to theoretical results.

Gruson, M F↗

Differential rotation and turbulent convection: A new Reynolds stress model and comparison with solar data

In most hydrodynamic cases, the existence of a turbulent flow superimposed on a mean flow is caused by a shear instability in the latter. Boussinesq suggested the first model for the turbulent Reynolds stresses bar-(u(sub i)u(sub j)) in which the mean shear S(sub ij) is the cause (or source) of turbulence represented by the stress bar-(u(sub i)u(sub j)). In the case of solar differential rotation, exactly the reverse physical process occurs: turbulence (which must pre-exist) generates a mean flow which manifests itself in the form of differential rotation. Thus, the Boussinesq model is wholly inadequate because in the solar case, cause and effect are reversed. Since the Boussinesq model is inadequate, one needs an alternative model for the Reynolds stresses. We present a new dynamical model for the Reynolds stresses, convective fluxes, turbulent kinetic energy, and temperature fluctuations. The complete model requires the solution of 11 differential equations. We then introduce a set of simplifying assumptions which reduce the full dynamical model to a set of algebraic Reynolds stress models. We explicitly solve one of these models that entails only one differential equation. The overall agreement with the data is obtained with a model that is neither phenomenological nor one that requires a full numerical simulation, since it is algebraic in nature. The new model can play an important role in understanding the complex physics underlying the interplay between solar differential rotation and convection, as many physical processes can naturally be incorporated into the model.

Canuto, V. M.↗

The Global Modeling and Assimilation Office (GMAO) 4d-Var and its Adjoint-based Tools

The fifth generation of the Goddard Earth Observing System (GEOS-5) Data Assimilation System (DAS) is a 3d-var system that uses the Grid-point Statistical Interpolation (GSI) system developed in collaboration with NCEP, and a general circulation model developed at Goddard, that includes the finite-volume hydrodynamics of GEOS-4 wrapped in the Earth System Modeling Framework and physical packages tuned to provide a reliable hydrological cycle for the integration of the Modern Era Retrospective-analysis for Research and Applications (MERRA). This MERRA system is essentially complete and the next generation GEOS is under intense development. A prototype next generation system is now complete and has been producing preliminary results. This prototype system replaces the GSI-based Incremental Analysis Update procedure with a GSI-based 4d-var which uses the adjoint of the finite-volume hydrodynamics of GEOS-4 together with a vertical diffusing scheme for simplified physics. As part of this development we have kept the GEOS-5 IAU procedure as an option and have added the capability to experiment with a First Guess at the Appropriate Time (FGAT) procedure, thus allowing for at least three modes of running the data assimilation experiments. The prototype system is a large extension of GEOS-5 as it also includes various adjoint-based tools, namely, a forecast sensitivity tool, a singular vector tool, and an observation impact tool, that combines the model sensitivity tool with a GSI-based adjoint tool. These features bring the global data assimilation effort at Goddard up to date with technologies used in data assimilation systems at major meteorological centers elsewhere. Various aspects of the next generation GEOS will be discussed during the presentation at the Workshop, and preliminary results will illustrate the discussion.

Todling, Ricardo↗

Two-fluid simulations of galaxy formation

We investigate the formation of galaxies and larger structure with a simulation modeling two gravitationally coupled fluids representing dark matter and baryons. The baryon gas dynamics are calculated with a smoothed particle hydrodynamics (SPH) method, and the physics modeled includes thermal pressure, shock heating, and radiative cooling. We simulate a 16 Mpc periodic cube with 64(exp 3) particles in each fluid and 10% baryon mass fraction. We confirm, for the first time experimentally, disk formation as a natural consequence of hierarchical clustering in a large-scale cosmological environment. The majority of isolated galaxies exhibit centrifugally supported disks. A power-law relation between cold baryonic mass and maximum rotation velocity is found, M varies as nu(sub rot)(exp alpha) with alpha = 2.5 after correcting for differential numerical resolution. Both the spatial and velocity distributions of the simulated galaxies are biased with respect to the dark matter. A counts-in-cells analysis indicates that an unphysical degree of merging in the central cluster is likely responsible for the antibias signal in the correlation function. A robust, scale-dependent velocity bias is measured. The ratio of galaxy to dark matter pairwise velocity dispersions on a scale of 1 Mpc is 0.7. The amplitude is only mildly dependent on redshift or mass cutoff and scales with separation as r(exp 0.2). The degree to which these results depend on numerical parameters is discussed. Mass resolution plays a key role in controlling the resulting fraction of cold, dense baryons. The mass fraction associated with galaxies decreases by a factor of approximately greater than 3 when the mass per particle is increased by a factor 8. Photoionization and energy input from supernova will have to be included to determine more carefully the fraction of highly dissipated material and the characteristics of the stellar component of galaxies.

Evrard, August E.↗

A Systems Approach to the Solid Lubrication of Foil Air Bearings for Oil-Free Turbomachinery

Foil air bearings are self-acting hydrodynamic bearings which rely upon solid lubricants to reduce friction and minimize wear during sliding which occurs at start-up and shut-down when surface speeds are too low to allow the formation of a hydrodynamic air film. This solid lubrication is typically accomplished by coating the non-moving foil surface with a thin, soft polymeric film. The following paper introduces a systems approach in which the solid lubrication is provided by a combination of self lubricating shaft coatings coupled with various wear resistant and lubricating foil coatings. The use of multiple materials, each providing different functions is modeled after oil-lubricated hydrodynamic sleeve bearing technology which utilizes various coatings and surface treatments in conjunction with oil lubricants to achieve optimum performance. In this study, room temperature load capacity tests are performed on journal foil air bearings operating at 14,000 rpm. Different shaft and foil coating technologies such as plasma sprayed composites, ceramic, polymer and inorganic lubricant coatings are evaluated as foil bearing lubricants. The results indicate that bearing performance is improved through the individual use of the lubricants and treatments tested. Further, combining several solid lubricants together yielded synergistically better results than any material alone.

DellaCorte, Christopher↗

Formation of the protosolar nebula

Theoretical models are discussed of the collapse of a dense molecular cloud core to form the protosolar nebula that produce the sun and the planet. The theoretical models use the equations of hydrodynamics, gravitation, and radiative transfer to follow the time evolution of a cloud collapsing under its own self-gravity. Both semi-analytical and fully numerical solutions (in two and three spatial dimensions) were calculated by several workers, One challenge is to find a set of initial conditions that will lead to the formation of a suitable protosolar nebula. Detailed results are shown for 2-D models, both with and without turbulent viscosity for redistributing angular momentum, and for 3-D models investigating the strength of gravitational torques associated with nonaxisymmetry produced during the collapse phase.

Tscharnuter, Werner M.↗

Interrogating Solar Flare Loop Models With IRIS Observations 1: Overview of the Models, and Mass Flows

Solar flares are transient yet dramatic events in the atmosphere of the Sun, during which a vast amount of magnetic energy is liberated. This energy is subsequently transported through the solar atmosphere or into the heliosphere, and together with coronal mass ejections flares comprise a fundamental component of space weather. Thus, understanding the physical processes at play in flares is vital. That understanding often requires the use of forward modelling in order to predict the hydrodynamic and radiative response of the solar atmosphere. Those predictions must then be critiqued by observations to show us where our models are missing ingredients. While flares are of course 3D phenomenon, simulating the flaring atmosphere including an accurate chromosphere with the required spatial scales in 3D is largely beyond current computational capabilities, and certainly performing parameter studies of energy transport mechanisms is not yet tractable in 3D. Therefore, field-aligned 1D loop models that can resolve the relevant scales have a crucial role to play in advancing our knowledge of flares. In recent years, driven in part by the spectacular observations from the Interface Region Imaging Spectrograph (IRIS), flare loop models have revealed many interesting features of flares. For this review I highlight some important results that illustrate the utility of attacking the problem of solar flares with a combination of high quality observations, and state-of-the-art flare loop models, demonstrating: 1) how models help to interpret flare observations from IRIS, 2) how those observations show us where we are missing physics from our models, and 3) how the ever increasing quality of solar observations drives model improvements. Here in Paper one of this two part review I provide an overview of modern flare loop models, and of electron-beam driven mass flows during solar flares.

solar flares↗

NICER Observes a Secondary Peak in the Decay of a Thermonuclear Burst from 4U 1608-52

We report for the first time below 1.5 keV, the detection of a secondary peak in an Eddington-limited thermonuclear X-ray burst observed by the Neutron Star Interior Composition Explorer (NICER) from the low-mass X-ray binary 4U 1608–52. Our time-resolved spectroscopy of the burst is consistent with a model consisting of a varying-temperature blackbody, and an evolving persistent flux contribution, likely attributed to the accretion process. The dip in the burst intensity before the secondary peak is also visible in the bolometric flux. Prior to the dip, the blackbody temperature reached a maximum of ≈3 keV. Our analysis suggests that the dip and secondary peak are not related to photospheric expansion, varying circumstellar absorption, or scattering. Instead, we discuss the observation in the context of hydrodynamical instabilities, thermonuclear flame spreading models, and reburning in the cooling tail of the burst.

Gaurava K. Jaisawal↗

Two-dimensional numerical simulation of a Stirling engine heat exchanger

The first phase of an effort to develop multidimensional models of Stirling engine components is described; the ultimate goal is to model an entire engine working space. More specifically, parallel plate and tubular heat exchanger models with emphasis on the central part of the channel (i.e., ignoring hydrodynamic and thermal end effects) are described. The model assumes: laminar, incompressible flow with constant thermophysical properties. In addition, a constant axial temperature gradient is imposed. The governing equations, describing the model, were solved using Crank-Nicloson finite-difference scheme. Model predictions were compared with analytical solutions for oscillating/reversing flow and heat transfer in order to check numerical accuracy. Excellent agreement was obtained for the model predictions with analytical solutions available for both flow in circular tubes and between parallel plates. Also the heat transfer computational results are in good agreement with the heat transfer analytical results for parallel plates.

Ibrahim, Mounir B.↗

Electrochromatographic separation of proteins

We have developed a modified electrochromatography system which minimizes Joule heating at electric field strengths up to 125 V/cm. A non-linear equilibrium model is described which incorporates electrophoretic mobility, hydrodynamic flow velocity, and an electrically induced concentration polarization at the surface of the stationary phase. This model is able to provide useful estimates of protein retention time and velocity in a column packed with Sephadex gel and subjected to an electric field. A correlation of electrophoretic mobility of peptide and proteins with respect to their charge, molecular mass, and asymmetry enables the selection of solute target molecules for electrochromatographic separations. Good separation of protein mixtures have been obtained.

NASA Discipline Life Support Systems↗

Design of Fiber Optic Sensors for Measuring Hydrodynamic Parameters

The science of optical hydrodynamics involves relating the optical properties to the fluid dynamic properties of a hydrodynamic system. Fiber-optic sensors are being designed for measuring the hydrodynamic parameters of various systems. As a flowing fluid makes an encounter with a flat surface, it forms a boundary layer near this surface. The region between the boundary layer and the flat plate contains information about parameters such as viscosity, compressibility, pressure, density, and velocity. An analytical model has been developed for examining the hydrodynamic parameters near the surface of a fiber-optic sensor. An analysis of the conservation of momentum, the continuity equation and the Navier-Stokes equation for compressible flow were used to develop expressions for the velocity and the density as a function of the distance along the flow and above the surface. When examining the flow near the surface, these expressions are used to estimate the sensitivity required to perform direct optical measurements and to derive the shear force for indirect optical measurements. The derivation of this result permits the incorporation of better design parameters for other fiber-based sensors. Future work includes analyzing the optical parametric designs of fiber-optic sensors, modeling sensors to utilize the parameters for hydrodynamics and applying different mixtures of hydrodynamic flow. Finally, the fabrication of fiber-optic sensors for hydrodynamic flow applications of the type described in this presentation could enhance aerospace, submarine, and medical technology.

Lyons, Donald R.↗

An analysis of the high-latitude thermospheric wind pattern calculated by a thermospheric general circulation model. I - Momentum forcing

A diagnostic processor (DP) was developed for analysis of hydrodynamic and thermodynamic processes predicted by the NCAR thermospheric general circulation model (TGCM). The TGCM contains a history file on the projected wind, temperature and composition fields at each grid point for each hour of universal time. The DP assimilates the history file plus ion drag tensors and drift velocities, specific heats, coefficients of viscosity, and thermal conductivity and calculates the individual forcing terms for the momentum and energy equations for a given altitude. Sample momentum forcings were calculated for high latitudes in the presence of forcing by solar radiation and magnetospheric convection with a 60 kV cross-tail potential, i.e., conditions on Oct. 21, 1981. It was found that ion drag and pressure forces balance out at F region heights where ion drift velocities are small. The magnetic polar cap/auroral zone boundary featured the largest residual force or net acceleration. Diurnal oscillations were detected in the thermospheric convection, and geostrophic balance was dominant in the E layer.

Killeen, T. L.↗