Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Inflow/outflow boundary conditions and global dynamics of spatial mixing layers

The numerical simulation of incompressible spatially-developing shear flows poses a special challenge to computational fluid dynamicists. The Navier-Stokes equations are elliptic and boundary equations need to be specified at the inflow and outflow boundaries in order to compute the fluid properties within the region of interest. It is, however, difficult to choose inflow and outflow conditions corresponding to a given experimental situation. Furthermore the effects that changes in the boundary conditions or in the size of the computational domain may induce on the global dynamics of the flow are presently unknown. These issues are examined in light of recent developments in hydrodynamic stability theory. The particular flow considered is the spatial mixing layer but it was expected that similar phenomena were bound to occur in other cases such as channel flow, the boundary layer, etc. A short summary of local/global and absolute/convective instability concepts is given. The results of numerical simulations are presented which strongly suggest that global resonances may be triggered in domains of finite streamwise extent although the evolution of the perturbation vorticity field is everywhere locally convective. A relationship between finite domains and pressure sources which might help in devising a scheme to eliminate these difficulties is discussed.

Buell, J. C.↗
NASA NTRS1990

Monte Carlo simulation of nonadiabatic expansion in cometary atmospheres - Halley

Monte Carlo methods developed for the characterization of velocity-dependent collision processes and ballistic transports in planetary exospheres form the basis of the present computer simulation of icy comet atmospheres, which iteratively undertakes the simultaneous determination of velocity distribution for five neutral species (water, together with suprathermal OH, H2, O, and H) in a flow regime varying from the hydrodynamic to the ballistic. Experimental data from the neutral mass spectrometer carried by Giotto for its March, 1986 encounter with Halley are compared with a model atmosphere.

Hodges, R. Richard, Jr.↗
NASA NTRS1994

Synthesized spectra of turbulent clouds

A comparison has been made of the velocity structure of non-star-forming dense molecular clouds and the output of a 512(exp 3) hydrodynamic code for compressible gas displaying features of turbulence. A remarkable similarity exists between the interstellar line profiles and the velocity profiles found by integrating, with an optically thin assumption, through subsections of the simulation. This similarity is in terms of both the structure of individual spectra, showing skewness, multiple peaks, and wings, and also the variations of the spectra from point to point. It is possible to see the origin of the various velocity substructures and how they correlate with regions of high density and velocity divergence, or regions of high vorticity, by deprojecting individual spectra of the simulation. It is found that the interstellar medium spectra are best represented by the epoch of the simulation in which the energy spectrum is essentially Kolmogorov-like and initially developed shocks have been largely converted into vorticity structures. Further, it is found that the vorticity peaks ascribed to intermittency contribute to the non-Gaussian features and wings of the line profile, but at an intensity level that is weak compared to the contribution of the bulk of the flow.

Falgarone, E.↗
NASA NTRS1995

Unsteady numerical simulations of the stability and dynamics of flames

In this report we describe the research performed at the Naval Research Laboratory in support of the NASA Microgravity Science and Applications Program over the past three years (from Feb. 1992) with emphasis on the work performed since the last microgravity combustion workshop. The primary objective of our research is to develop an understanding of the differences in the structure, stability, dynamics and extinction of flames in earth gravity and in microgravity environments. Numerical simulations, in which the various physical and chemical processes can be independently controlled, can significantly advance our understanding of these differences. Therefore, our approach is to use detailed time-dependent, multi-dimensional, multispecies numerical models to perform carefully designed computational experiments. The basic issues we have addressed, a general description of the numerical approach, and a summary of the results are described in this report. More detailed discussions are available in the papers published which are referenced herein. Some of the basic issues we have addressed recently are (1) the relative importance of wall losses and gravity on the extinguishment of downward-propagating flames; (2) the role of hydrodynamic instabilities in the formation of cellular flames; (3) effects of gravity on burner-stabilized flames, and (4) effects of radiative losses and chemical-kinetics on flames near flammability limits. We have also expanded our efforts to include hydrocarbon flames in addition to hydrogen flames and to perform simulations in support of other on-going efforts in the microgravity combustion sciences program. Modeling hydrocarbon flames typically involves a larger number of species and a much larger number of reactions when compared to hydrogen. In addition, more complex radiation models may also be needed. In order to efficiently compute such complex flames recent developments in parallel computing have been utilized to develop a state-of-the-art parallel flame code. This is discussed below in some detail after a brief discussion of the numerical models.

Kailasanath, K.↗
NASA NTRS1993

Hydrodynamic models of the Cartwheel ring galaxy

A series of increasingly sophisticated models of the Cartwheel ring galaxy is studied in order to test the collisional model for the galaxy formation and examine the star formation processes in this unique environment, using new data acquired in the last decade. The simulations provided some possible answers to a number of questions about the Cartwheel. First, an explanation for the wide spacing between inner and outer rings is suggested by the simple epicyclic kinematics within the dark matter-dominated potential implied by H I rotation curve. These models and the kinematic model of Struck-Marcell and Lotan (1990) also predict that the outer ring should be relatively weak, while the second inner ring should be stronger, with a dense orbit-crossing region of significant width bounded by sharp, caustic edges. The collisional model is given support by the agreement between the observations and the morphological and kinematic properties of the numerical simulations presented.

Struck-Marcell, Curtis↗
NASA NTRS1995

Solar wind/local interstellar medium interaction including charge exchange with neural hydrogen

We present results from a hydrodynamic model of the interaction of the solar wind with the local interstellar medium (LISM), self-consistently taking into account the effects of charge exchange between the plasma component and the interstellar neutrals. The simulation is fully time dependent, and is carried out in two or three dimensions, depending on whether the helio-latitudinal dependence of the solar wind speed and number density (both giving rise to three dimensional effects) are included. As a first approximation it is assumed that the neutral component of the flow can be described by a single, isotropic fluid. Clearly, this is not the actual situation, since charge exchange with the supersonic solar wind plasma in the region of the nose results in a 'second' neutral fluid propagating in the opposite direction as that of the LISM neutrals.

Pauls, H. Louis↗
NASA NTRS2005

On the Origin of the Saturnian Satellite System: Did Iapetus Form In-Situ?

Current models of planet and satellite formation are marred by our lack of understanding regarding the turbulent state of accretion disks. According to the Rayleigh criterion, Keplerian disks are hydrodynamically stable. Indeed, it has been argued that a carefully designed Taylor Couette experiment shows stability in the case of positive radial gradients in specific angular momentum even for high Reynolds numbers [1], in agreement with numerical simulations which consistently show turbulence decay [2]. Other possible sources of turbulence may fail due to low ionization, may decay as the optical depth decreases due to dust coagulation, may involve unrealistic boundary conditions, or result in limited transport. The difficulty stems not only from the degree of turbulence, but also from the kind of turbulence, and whether it may be characterized by an alpha parameter.

Mosqueira, I.↗
NASA NTRS1989

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

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. Parallel plate and tubular heat exchanger models are described, with emphasis on the central part of the channel (i.e., ignoring hydrodynamic and thermal end effects). 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 have been solved using the Crack-Nicloson finite-difference scheme. Model predictions are compared with analytical solutions for oscillating/reversing flow and heat transfer in order to check numerical accuracy. Excellent agreement is obtained for flow both in circular tubes and between parallel plates. The computational heat transfer results are in good agreement with the analytical heat transfer results for parallel plates.

Ibrahim, Mounir↗
NASA NTRS2007

Core Physics and Kinetics Calculations for the Fissioning Plasma Core Reactor

Highly efficient, compact nuclear reactors would provide high specific impulse spacecraft propulsion. This analysis and numerical simulation effort has focused on the technical feasibility issues related to the nuclear design characteristics of a novel reactor design. The Fissioning Plasma Core Reactor (FPCR) is a shockwave-driven gaseous-core nuclear reactor, which uses Magneto Hydrodynamic effects to generate electric power to be used for propulsion. The nuclear design of the system depends on two major calculations: core physics calculations and kinetics calculations. Presently, core physics calculations have concentrated on the use of the MCNP4C code. However, initial results from other codes such as COMBINE/VENTURE and SCALE4a. are also shown. Several significant modifications were made to the ISR-developed QCALC1 kinetics analysis code. These modifications include testing the state of the core materials, an improvement to the calculation of the material properties of the core, the addition of an adiabatic core temperature model and improvement of the first order reactivity correction model. The accuracy of these modifications has been verified, and the accuracy of the point-core kinetics model used by the QCALC1 code has also been validated. Previously calculated kinetics results for the FPCR were described in the ISR report, "QCALC1: A code for FPCR Kinetics Model Feasibility Analysis" dated June 1, 2002.

Butler, C.↗
NASA NTRS1991

Accretion in wind-driven X-ray sources

The hydrodynamics of axisymmetric accretion flow in stellar wind-fed X-ray sources is investigated, including momentum deposition by radiation, as well as radiative heating and cooling by Compton and bremsstrahlung processes. The results of two-dimensional numerical simulations for mass accretion rates ranging from 0.02 to 0.4 of the Eddington value reveal a variety of behaviors. At low rates radiative effects are unimportant, and the accretion flow is nearly adiabatic. In this case the flow is steady. For intermediate accretion rates, radiative heating and cooling effects become important. The flow remains steady, but the accretion rate is significantly reduced from the Hoyle-Lyttleton estimate. At the highest mass accretion rates, radiative momentum transfer and energy exchange are important and lead to nonsteady behavior. For the parameters relevant to massive X-ray binary systems accreting at high rates, the results reveal that the accretion flows can be complex and time dependent.

Taam, Ronald E.↗
NASA NTRS1999

The Micromechanics of the Moving Contact Line

A transient moving contact line is investigated experimentally. The dynamic interface shape between 20 and 800 microns from the contact line is compared with theory. A novel experiment is devised, in which the contact line is set into motion by electrically altering the solid-liquid surface tension gamma(sub SL). The contact line motion simulates that of spontaneous wetting along a vertical plate with a maximum capillary number Ca approx. = 4 x 10(exp -2). The images of the dynamic meniscus are analyzed as a funtion of Ca. For comparison, the steady-state hydrodynamic equation based on the creeping flow model in a wedge geometry and the three-region uniform perturbation expansion of Cox (1986) is adopted. The interface shape is well depicted by the uniform solutions for Ca <= 10(exp -3). However, for Ca > 10(exp -3), the uniform solution over-predicts the viscous bending. This over-prediction can be accounted for by modifying the slip coefficient within the intermediate solution. With this correction, the measured interface shape is seen to match the theoretical prediction for all capillary numbers. The amount of slip needed to fit the measurements does not scale with the capillary number.

Lichter, Seth↗
NASA NTRS2012

Pros and Cons of ID vs. 3D Modeling

Advances in computing capability have led to tremendous improvements in 3D modeling. Entire active regions are being simulated in what might be described as a first principles way, in which plasma heating is treated self consistently rather than through the specification of heating functions. There are limitations to this approach, however, as actual heating mechanisms on the Sun involve spatial scales orders of magnitude smaller than what these simulations can resolve. Other simulations begin to resolve these scales, but they only treat a tiny volume and do not include the all important coupling with larger scales or with other parts of the atmosphere, and so cannot be readily compared with observations. Finally, ID hydrodynamic models capture the field-aligned evolution of the plasma extremely well and are ideally suited for data comparison, but they treat the heating in a totally ad hoc manner. All of these approaches have important contributions to make, but we must be aware of their limitations. I will highlight some of the strengths. and weaknesses of each.

Klimchuk, James A.↗
NASA NTRS1988

The simulation of the alternate turbopump development high pressure oxygen and fuel turbopumps for the space shuttle main engine using the Shaberth computer program

The Space Shuttle Main Engine (SSME) is basically comprised of a combustion chamber and nozzle, high and low pressure oxygen turbopumps and high and low pressure fuel turbopumps. In the current configuration, the high pressure fuel (HPTFP) and high pressure oxygen turbopumps (HPOTP) have experienced a history of ball bearing wear. The wear problem can be attributed to numerous factors including the hydrodynamic axial and radial loads caused by the flow of liquid oxygen and liquid hydrogen through the turbopump impellers and turbine. Also, friction effects between the rolling elements, races, and cage can create thermally induced bearing geometry changes. To alleviate some of the current configuration problems, an alternate turbopump development (ATD) was proposed. However, the ATD HPOTP and HPTFP are constrained to operate interchangeably with the current turbopumps, thus, the operation conditions must be similar. The ATD configuration features a major change in bearings used to support the integrated shaft, impeller, and turbine system. A single ball and single roller will replace the pump-end and turbine and duplex ball bearings. The Shaft-Bearing-Thermal (SHABERTH) computer code was used to model the ATD HPOTP and ATD HPFTP configurations. A two bearing model was used to simulate the HPOTP and HPFTP bearings and shaft geometry. From SHABERTH, a comparison of bearing reaction loads, frictional heat generation rates, and Hertz contact stresses will be attempted with analysis at the 109 percent and 65 percent power levels.

Mcdonald, Gary H.↗
NASA NTRS1990

Hydrodynamic instability and coherent structures in turbulence

The evolution of a coherent structure in the near wall region is examined in the context of linear theory. In accordance with this theory, a three-dimensional inviscid disturbance is shown to grow into a streamwise elongated streaky structures in a manner similar to what has been observed in experiments and numerical simulations. It is found that spanwise asymmetry is necessary for the structures to evolve into long streaks. The implications of the results for turbulent drag reduction through active or passive means are discussed.

Landahl, Marten T.↗
NASA NTRS1992

The three-dimensional evolution of a plane mixing layer - The Kelvin-Helmholtz rollup

The hydrodynamic evolution of an incompressible plane mixing layer is addressed to elucidate scalar mixing in free shear flows. A detailed description of the onset of three-dimensionality in a mixing layer before or in the absence of pairing is presented. Various simulations were performed to investigate the sensitivity of these results to variations in initial conditions. These variations included changes in amplitude, wavelength, functional form, and relative phasing of the initial low-wavenumber disturbances. Pierrehumber and Widnall's (1982) translative instability eigenfunctions are found to include rib vortices in the braid region and oppositely signed streamwise vorticity in the roller core. The translative instability is an instability of the late-time oversaturated flow. Three-dimensional perturbation growth similar to that of the translative instability can occur whenever spanwise vorticity is present in the braid region. The nonlinear effects that occur when the initial rib circulation is sufficiently high are discussed.

Rogers, Michael M.↗
NASA NTRS1992

Magellan observations of extended impact crater related features on the surface of Venus

The hydrodynamic evolution of an incompressible plane mixing layer is addressed to elucidate scalar mixing in free shear flows. A detailed description of the onset of three-dimensionality in a mixing layer before or in the absence of pairing is presented. Various simulations were performed to investigate the sensitivity of these results to variations in initial conditions. These variations included changes in amplitude, wavelength, functional form, and relative phasing of the initial low-wavenumber disturbances. Pierrehumber & Widnall's (1982) translative instability eigenfunctions are found to include rib vortices in the braid region and oppositely signed streamwise vorticity in the roller core. The translative instability is an instability of the late-time oversaturated flow. Three-dimensional perturbation growth similar to that of the translative instability can occur whenever spanwise vorticity is present in the braid region. The nonlinear effects that occur when the initial rib circulation is sufficiently high are discussed.

Campbell, D. B.↗
NASA NTRS1994

Studies on equatorial shock formation during plasmaspheric refilling

Investigations based on small-scale simulations of microprocesses occurring when a magnetic flux tube refills with a cold plasma are summarized. Results of these investigations are reported in the following attached papers: (1) 'Numerical Simulation of Filling a Magnetic Flux Tube with a Cold Plasma: The Role of Ion Beam-Driven Instabilities'; and (2) 'Numerical Simulation of Filling a Magnetic Flux Tube with a Cold Plasma: Effects of Magnetically Trapped Hot Plasma'. Other papers included are: 'Interaction of Field-Aligned Cold Plasma Flows with an Equatorially-Trapped Hot Plasma: Electrostatic Shock Formation'; and 'Comparison of Hydrodynamic and Semikinetic Treatments for a Plasma Flow along Closed Field Lines'. A proposal for further research is included.

Singh, N.↗