Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Flow dynamics”

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 685 records · Page 38

Dynamic localization and second-order subgrid-scale models in large eddy simulations of channel flow

The objective here is to test the Dynamic Localization (DL) model in a wall-bounded channel flow for numerical stability and accuracy of results. Algebraic stress models suggest that the model for the residual subgrid-scale (SGS) Reynolds stress and scalar flux should generally have terms comprising most of the unique products of the resolved strain (S) and rotation (R) tensors with S and the resolved scalar gradient. The standard dynamic SGS model uses a simple (Smagorinsky) base model for the residual Reynolds stress, which is made proportional to S, and down-gradient base models for residual scalar fluxes; these correspond to the lowest, 'first-order' terms in algebraic stress models. Temporal scaling terms in these base models are formed from the magnitude of the resolved strain rate. While this is appropriate for simple shear flows, it may not be appropriate for more complicated flows (relevant to geophysical and astrophysical problems) that include any combination of shear, rotation, buoyancy, etc. On the other hand, the coefficient in the dynamic SGS model readily adjusts itself to different flow conditions and may adequately take account of these effects without the need for more complicated base models. Cabot (1993) has begun to test the dynamic SGS model in buoyant flows (Rayleigh-Benard and internally heated convection) with and without buoyancy terms explicitly included in the scaling terms of the base model; no great differences were found in large eddy simulation (LES) results for the different base model scalings. The second objective in this work is to test base models with additional, 'second-order' terms (e.g., S(sup 2) and RS for the residual Reynolds stress). These terms have been found to improve large-scale flow predictions by kappa-epsilon models in the presence of rotation and shear. Second-order base models will be tested here in the LES of channel flow with and without solid-body rotation and compared with results from the standard first-order base models to determine if there are significant differences or improvements in results that would warrant the added complexity of the second-order base models.

Cabot, William H.↗

Transitioning of power flow in beam models with bends

The propagation of power flow through a dynamically loaded beam model with 90 degree bends is investigated using NASTRAN and McPOW. The transitioning of power flow types (axial, torsional, and flexural) is observed throughout the structure. To get accurate calculations of the torsional response of beams using NASTRAN, torsional inertia effects had to be added to the mass matrix calculation section of the program. Also, mass effects were included in the calculation of BAR forces to improve the continuity of power flow between elements. The importance of including all types of power flow in an analysis, rather than only flexural power, is indicated by the example. Trying to interpret power flow results that only consider flexural components in even a moderately complex problem will result in incorrect conclusions concerning the total power flow field.

Hambric, Stephen A.↗

Effect of Air Swirler Configuration on Lean Direct Injector Flow Structure and Combustion Performance with a 7-Point Lean Direct Injector Array

Studies of various injector configurations in a 7-point Lean Direct Injector (LDI) array are reported for both non-reacting (cold) flow and for Jet-A/air reacting flows. For cold flow, central recirculation zone (CRZ) formation is investigated and for reacting flows, combustor operability and dynamics are of interest. 2D Particle Image Velocimetry (PIV) measurements are described for the cold flow experiments and flame chemiluminescence imaging and dynamic pressure results are discussed for the reacting flow cases. PIV results indicate that for this configuration the close spacing between swirler elements leads to strong interaction that affects whether a CRZ forms, and pilot recess and counter-swirl helps to isolate swirlers from one another. Dynamics results focus on features identified near 500-Hz.

aircraft combustors↗

Determination of Heritage SSME Pogo Suppressor Resistance and Inertance from Waterflow Pulse Testing

Waterflow tests of a heritage Space Shuttle Main Engine pogo suppressor were performed to experimentally quantify the resistance and inertance provided by the suppressor. Measurements of dynamic pressure and flow rate in response to pulsing flow were made throughout the test loop. A unique system identification methodology combined all sensor measurements with a one-dimensional perturbational flow model of the complete water flow loop to spatially translate physical measurements to the device under test. Multiple techniques were then employed to extract the effective resistance and inertance for the pogo suppressor. Parameters such as steady flow rate, perturbational flow rate magnitude, and pulse frequency were investigated to assess their influence on the behavior of the pogo suppressor dynamic response. These results support validation of the RS-25 pogo suppressor performance for use on the Space Launch System Core Stage.

McDougal, Chris↗

Quasi One-Dimensional Unsteady Modeling of External Compression Supersonic Inlets

The AeroServoElasticity task under the NASA Supersonics Project is developing dynamic models of the propulsion system and the vehicle in order to conduct research for integrated vehicle dynamic performance. As part of this effort, a nonlinear quasi 1-dimensional model of an axisymmetric external compression supersonic inlet is being developed. The model utilizes compressible flow computational fluid dynamics to model the internal inlet segment as well as the external inlet portion between the cowl lip and normal shock, and compressible flow relations with flow propagation delay to model the oblique shocks upstream of the normal shock. The external compression portion between the cowl-lip and the normal shock is also modeled with leaking fluxes crossing the sonic boundary, with a moving CFD domain at the normal shock boundary. This model has been verified in steady state against tunnel inlet test data and it s a first attempt towards developing a more comprehensive model for inlet dynamics.

Kopasakis, George↗

Rayleigh Scattering Diagnostic for Simultaneous Measurements of Dynamic Density and Velocity

A flow diagnostic technique based on the molecular Rayleigh scattering of laser light is used to obtain dynamic density and velocity data in turbulent flows. The technique is based on analyzing the Rayleigh scattered light with a Fabry-Perot interferometer and recording information about the interference pattern with a multiple anode photomultiplier tube (PMT). An artificial neural network is used to process the signals from the PMT to recover the velocity time history, which is then used to calculate the velocity power spectrum. The technique is illustrated using simulated data. The results of an experiment to measure the velocity power spectrum in a low speed (100 rn/sec) flow are also presented.

Seasholtz, Richard G.↗

Magnectic Probing of Core Geodynamics

To better understand geomagnetic theory and observation, we can use spatial magnetic spectra for the main field and secular variation to test core dynamical hypotheses against seismology. The hypotheses lead to theoretical spectra which are fitted to observational spectra. Each fit yields an estimate of the radius of Earth s core and uncertainty. If this agrees with the seismologic value, then the hypotheses pass the test. A new way to obtain theoretical spectra extends the hydromagnetic scale analysis of Benton to scale-variant field and flow. For narrow scale flow and a dynamically weak field by the top of Earth s core, this yields a JGR-PI, and a secular variation spectrum modulated by a cubic polynomial in spherical harmonic degree n. The former passes the tests. The latter passes many tests, but does not describe rapid dipole decline and quadrupole rebound; some tests suggest it is a bit hard, or rich in narrow scale change.In a core geodynamo, motion of the fluid conductor does work against the Lorentz force. This converts kinetic into magnetic energy which, in turn, is lost to heat via Ohmic dissipation. In the analysis at length- scale l/k, if one presumes kinetic energy is converted in either eddy- overturning or magnetic free-decay time-scales, then Kolmogorov or other spectra in conflict with observational spectra can result. Instead, the rate work is done roughly balances the dissipation rate, which is consistent with small scale flow. The conversion time-scale depends on dynamical constraints. These are summarized by the magneto- geostrophic vertical vorticity balance by the top of the core, which includes anisotropic effects of rotation, the magnetic field, and the core-mantle boundary. The resulting theoretical spectra for the core- source field and its SV are far more compatible with observation. The conversion time-scale of order 120 years is pseudo-scale-invariant. Magnetic spectra of other planets may differ; however, if a transition to non-conducting fluid hydrogen in Jupiter acts as barrier to vertical flow, as well as current, then the shape of the jovi-magnetic spectrum could be remarkably Earth-like.

Voorhies, Coerte↗

Magnetic Probing of Core Geodynamics

To better understand geomagnetic theory and observation, we can use spatial magnetic spectra for the main field and secular variation to test core dynamical hypotheses against seismology. The hypotheses lead to theoretical spectra which are fitted to observational spectra. Each fit yields an estimate of the radius of Earth's core and uncertainty. If this agrees with the seismologic value, then the hypothesis passes the test. A new way to obtain theoretical spectra extends the hydromagnetic scale analysis of Benton to scale-variant field and flow. For narrow scale flow and a dynamically weak field by the top of Earth's core, this yields a generalized Stevenson-McLeod spectrum for the core-source field, and a secular variation spectrum modulated by a cubic polynomial in spherical harmonic degree n. The former passes the tests. The latter passes many tests, but does not describe rapid dipole decline and quadrupole rebound; some tests suggest it is a bit hard, or rich in narrow scale range. In a core geodynamo, motion of the fluid conductor does work against the Lorentz force. This converts kinetic into magnetic energy which, in turn, is lost to heat via Ohmic dissipation. In the analysis at length-scale 1/k, if one presumes kinetic energy is converted in either eddy-overturning or magnetic free-decay time-scales, then Kolmogorov or other spectra in conflict with observational spectra can result. Instead, the rate work is done roughly balances the dissipation rate, which is consistent with small-scale flow. The conversion time-scale depends on dynamical constraints. These are summarized by the magnetogeostrophic vertical vorticity balance by the top of the core, which includes anisotropic effects of rotation, the magnetic field, and the core-mantle boundary. The resulting theoretical spectra for the core-source field and its SV are far more compatible with observation. The conversion time-scale of order 120 years is pseudo-scale-invariant. Magnetic spectra of other planets may differ; however, if a transition to non-conducting fluid hydrogen in Jupiter acts as a barrier to vertical flow, as well as current, then the shape of the jovi-magnetic spectrum could be remarkably Earth-like.

Voorhies, Coerte V.↗

Reconstruction of the dynamics of the 1800-1801 Hualalai eruption: Implications for planetary lava flows

The 1800-1801 eruption of alkalic basalt from the Hualalai volcano, Hawaii provides a unique opportunity for investigating the dynamics of lava flow emplacement with eruption rates and compositions comparable to those that have been suggested for planetary eruptions. Field observations suggest new considerations must be used to reconstruct the emplacement of these lava flows. These observations are: (1) the flow traversed the 15 km from the vent to the sea so rapidly that no significant crust formed and an observation of the eruption reported that the flow reach the sea from the vent in approximately 1 hour; (2) the drainage of beds of xenolith nodules indicates a highly fluid, low viscosity lava; (3) overspills and other morphologic evidence for a very low viscosity host fluid; (4) no significant longitudinal increase in flow thickness that might be associated with an increase in the rheological properties of the lava; and (5) the relatively large size of channels associated with the flow, up to 80 meters across and several km long. Models for many geologic mass movements and fast moving fluids with various loadings and suspensions are discussed.

Baloga, Stephen↗

Geometry of thin liquid sheet flows

Incompresible, thin sheet flows have been of research interest for many years. Those studies were mainly concerned with the stability of the flow in a surrounding gas. Squire was the first to carry out a linear, invicid stability analysis of sheet flow in air and compare the results with experiment. Dombrowski and Fraser did an experimental study of the disintegration of sheet flows using several viscous liquids. They also detected the formulation of holes in their sheet flows. Hagerty and Shea carried out an inviscid stability analysis and calculated growth rates with experimental values. They compared their calculated growth rates with experimental values. Taylor studied extensively the stability of thin liquid sheets both theoretically and experimentally. He showed that thin sheets in a vacuum are stable. Brown experimentally investigated thin liquid sheet flows as a method of application of thin films. Clark and Dumbrowski carried out second-order stability analysis for invicid sheet flows. Lin introduced viscosity into the linear stability analysis of thin sheet flows in a vacuum. Mansour and Chigier conducted an experimental study of the breakup of a sheet flow surrounded by high-speed air. Lin et al. did a linear stability analysis that included viscosity and a surrounding gas. Rangel and Sirignano carried out both a linear and nonlinear invisid stability analysis that applies for any density ratio between the sheet liquid and the surrounding gas. Now there is renewed interest in sheet flows because of their possible application as low mass radiating surfaces. The objective of this study is to investigate the fluid dynamics of sheet flows that are of interest for a space radiator system. Analytical expressions that govern the sheet geometry are compared with experimental results. Since a space radiator will operate in a vacuum, the analysis does not include any drag force on the sheet flow.

Chubb, Donald L.↗

Computational Assessment of Inlet Backflow Effects on Rotating Detonation Engine Performance and Operability

The performance impact of flow reversal at the inlet of an airbreathing rotating detonation engine (RDE) is investigated using 2 and 3-dimensional computational fluid dynamic (CFD) simulations. Flow reversal, or backflow, occurs in RDE inlets in the high-pressure region directly behind the rotating detonation front. This is also where most of the engine thrust or pressure gain is produced. The amount of backflow relative to throughflow depends on the inlet design. For the present work, a simple annular ‘slit’ design is used. The simulations are idealized in several ways, including that fuel and air are premixed, but prevented from reacting when within the inlet region. The results indicate that even with idealizations, the impact of inlet backflow on pressure gain can be substantial. The simulations also reveal an intriguing instability that develops in certain configurations. The mass flow rate into the RDE begins to oscillate at a regular frequency that is substantially less than the detonation rotational frequency. This is accompanied by oscillations in the detonation height. The oscillation amplitude grows over time until the detonation ultimately fails. Both the performance and instability results emphasize the need for carefully designed RDE inlets that provide low loss when flow is in the forward direction, but high resistance when the flow is reversed. Development of such high-diodicity inlets is critical to achieving pressure gain in airbreathing RDE’s.

Detonation↗

Computational Assessment of Inlet Backflow Effects on Rotating Detonation Engine Performance and Operability

The performance impact of flow reversal at the inlet of an airbreathing rotating detonation engine (RDE) is investigated using 2 and 3-dimensional computational fluid dynamic (CFD) simulations. Flow reversal, or backflow, occurs in RDE inlets in the high-pressure region directly behind the rotating detonation front. This is also where most of the engine thrust or pressure gain is produced. The amount of backflow relative to throughflow depends on the inlet design. For the present work, a simple annular ‘slit’ design is used. The simulations are idealized in several ways, including that fuel and air are premixed, but prevented from reacting when within the inlet region. The results indicate that even with idealizations, the impact of inlet backflow on pressure gain can be substantial. The simulations also reveal an intriguing instability that develops in certain configurations. The mass flow rate into the RDE begins to oscillate at a regular frequency that is substantially less than the detonation rotational frequency. This is accompanied by oscillations in the detonation height. The oscillation amplitude grows over time until the detonation ultimately fails. Both the performance and instability results emphasize the need for carefully designed RDE inlets that provide low loss when flow is in the forward direction, but high resistance when the flow is reversed. Development of such high-diodicity inlets is critical to achieving pressure gain in airbreathing RDE’s.

detonation↗

Power flows and Mechanical Intensities in structural finite element analysis

The identification of power flow paths in dynamically loaded structures is an important, but currently unavailable, capability for the finite element analyst. For this reason, methods for calculating power flows and mechanical intensities in finite element models are developed here. Formulations for calculating input and output powers, power flows, mechanical intensities, and power dissipations for beam, plate, and solid element types are derived. NASTRAN is used to calculate the required velocity, force, and stress results of an analysis, which a post-processor then uses to calculate power flow quantities. The SDRC I-deas Supertab module is used to view the final results. Test models include a simple truss and a beam-stiffened cantilever plate. Both test cases showed reasonable power flow fields over low to medium frequencies, with accurate power balances. Future work will include testing with more complex models, developing an interactive graphics program to view easily and efficiently the analysis results, applying shape optimization methods to the problem with power flow variables as design constraints, and adding the power flow capability to NASTRAN.

Hambric, Stephen A.↗

A coupled fluid-structure dynamic analysis for water impact loads

In this paper a nonlinear transient, hydroelastic method is developed for response analysis of the Space Shuttle recoverable booster on water impact. At each time step, coupled equations of fluid flow and structure dynamics are solved with an iterative, self-correcting process. Rigid body motions are represented by large angle equations. Local deformations are represented by a NASTRAN-generated model of the three-dimensional structure. The fluid is represented by series solutions of the potential flow equation which include the effects of the local structure motions. The resulting water loads and structural stresses obtained will aid in qualification of the final structural design. Comparisons are made to experimental impact data to validate the method.

Herting, D. N.↗

The source of the Earth's long wavelength geoid anomalies: Implications for mantle and core dynamics

The long wavelength components of the Earth's gravity field result mainly from density contrasts associated with convection in the mantle. Direct interpretation of the geoid for mantle convection is complicated by the fact that convective flow results in dynamically maintained deformation of the surface of the Earth, the core mantle boundary (CMB), and any interior chemical boundaries which might exist. These boundary deformations effect the geoid opposite in sign and are comparable in magnitude to those of the interior density contrasts driving the flow. The total difference of two relatively large quantities.

Hager, B. H.↗

Experimental study of vortex breakdown in a cylindrical, swirling flow

The stability of a steady, vortical flow in a cylindrical container with one rotating endwall has been experimentally examined to gain insight into the process of vortex breakdowwn. The dynamics of the flow are governed by the Reynolds number (Re) and the aspect ratio of the cylinder. Re is given by Omega R(sup 2)/nu, where Omega is the speed of rotation of the endwall, R is the cylinder radius, and nu is the kinematic viscosity of the fluid filling the cylinder. The aspect ratio is H/R, where H is the height of the cylinder. Numerical simulation studies disagree whether or not the steady breakdown is stable beyond a critical Reynolds number, Re(sub c). Previous experimental researches have considered the steady and unsteady flows near Re(sub c), but have not explored the stability of the steady breakdown structures beyond this value. In this investigation, laser induced fluorescence was utilized to observe both steady and unsteady vortex breakdown at a fixed H/R of 2.5 with Re varying around Re(sub c). When the Re of a steady flow was slowly increased beyond Re(sub c), the breakdown structure remained steady even though unsteadiness was possible. In addition, a number of hysteresis events involving the oscillation periods of the unsteady flow were noted. The results show that both steady and unsteady vortex breakdown occur for a limited range of Re above Re(sub c). Also, with increasing Re, complex flow transformations take place that alter the period at which the unsteady flow oscillates.

Stevens, J. L.↗

Preliminary efforts toward development of data handling and analysis software for unsteady flow measurements: An application for aeroelastic transonic flow configurations

A few years ago the Structural Dynamics Division at LaRC started ambitious experimental research efforts known as the Benchmark Models Program. The primary objective of this program was to provide experimental data that may serve as a calibration source for computational fluid dynamics (CFD) efforts that deal with aeroelastic unsteady flow configurations. It also focuses on the understanding of complex flow phenomenon associated with unsteady flow developments. The overall plan for the program has been described by Bennett, including a presentation of initial test results of flutter of a rigid wing mounted on flexible supports. An example of a test model employed to measure the dynamic response along with corresponding pressure distributions is shown. This model incorporates eighty pressure transducers distributed along two spanwise stations. In addition, the models are equipped with four accelerometers and two strain gages. The data handling system for the Benchmark Model Program is under development. Several interactive computer routines designed for the user interface, dynamic memory allocation, unsteady flow measurements data extraction, and further data processing were developed. To present a few examples of measured data, the unsteady pressure distributions and the wing model dynamic response were plotted.

Finaish, Fathi↗