Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid wave”

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 595 records · Page 33

Numerical methods for turbulent flow

It has generally become accepted that the Navier-Strokes equations predict the dynamic behavior of turbulent as well as laminar flows of a fluid at a point in space away form a discontinuity such as a shock wave. Turbulence is also closely related to the phenomena of non-uniqueness of solutions of the Navier-Strokes equations. These second order, nonlinear partial differential equations can be solved analytically for only a few simple flows. Turbulent flow fields are much to complex to lend themselves to these few analytical methods. Numerical methods, therefore, offer the only possibility of achieving a solution of turbulent flow equations. In spite of recent advances in computer technology, the direct solution, by discrete methods, of the Navier-Strokes equations for turbulent flow fields is today, and in the foreseeable future, impossible. Thus the only economically feasible way to solve practical turbulent flow problems numerically is to use statistically averaged equations governing mean-flow quantities. The objective is to study some recent developments relating to the use of numerical methods to study turbulent flow.

Turner, James C., Jr.↗

Global convergence of inexact Newton methods for transonic flow

In computational fluid dynamics, nonlinear differential equations are essential to represent important effects such as shock waves in transonic flow. Discretized versions of these nonlinear equations are solved using iterative methods. In this paper an inexact Newton method using the GMRES algorithm of Saad and Schultz is examined in the context of the full potential equation of aerodynamics. In this setting, reliable and efficient convergence of Newton methods is difficult to achieve. A poor initial solution guess often leads to divergence or very slow convergence. This paper examines several possible solutions to these problems, including a standard local damping strategy for Newton's method and two continuation methods, one of which utilizes interpolation from a coarse grid solution to obtain the initial guess on a finer grid. It is shown that the continuation methods can be used to augment the local damping strategy to achieve convergence for difficult transonic flow problems. These include simple wings with shock waves as well as problems involving engine power effects. These latter cases are modeled using the assumption that each exhaust plume is isentropic but has a different total pressure and/or temperature than the freestream.

Young, David P.↗

Formation of intermediate shocks in both two-fluid and hybrid models

Intermediate shocks are studied in a two-fluid model that includes finite ion inertia dispersion and in a hybrid model in which the full ion dynamics is retained while the electrons are treated as a massless fluid. It is shown that in both models intermediate shocks can be formed through wave steepening, meaning that they are stable and possess shock structures.

Wu, C. C.↗

The flow of a thin liquid film on a stationary and rotating disk. I - Experimental analysis and flow visualization

The mean thickness of a thin liquid film of deionized water with a free surface on a stationary and rotating horizontal disk has been measured with a nonobtrusive capacitance technique. The measurements were taken when the rotational speed was 0-300 RPM and the flow rate was 7.0-15.0 LPM. A flow visualization study of the thin film was also performed to determine the characteristics of the waves on the free surface. When the disk was stationary, a circular hydraulic jump was present on the disk. Surface waves were found in the supercritical and subcritical regions at all flow rates studied. When the rotational speed of the disk is low, a standing wave at the edge of the disk was present. As the rotational speed increased, the surface waves changed from the wavy-laminar region to a region in which the waves ran nearly radially across the disk on top of a thin substrate of fluid.

Thomas, S.↗

A procedure for automating CFD simulations of an inlet-bleed problem

A procedure was developed to improve the turn-around time for computational fluid dynamics (CFD) simulations of an inlet-bleed problem involving oblique shock-wave/boundary-layer interactions on a flat plate with bleed into a plenum through one or more circular holes. This procedure is embodied in a preprocessor called AUTOMAT. With AUTOMAT, once data for the geometry and flow conditions have been specified (either interactively or via a namelist), it will automatically generate all input files needed to perform a three-dimensional Navier-Stokes simulation of the prescribed inlet-bleed problem by using the PEGASUS and OVERFLOW codes. The input files automatically generated by AUTOMAT include those for the grid system and those for the initial and boundary conditions. The grid systems automatically generated by AUTOMAT are multi-block structured grids of the overlapping type. Results obtained by using AUTOMAT are presented to illustrate its capability.

Chyu, Wei J.↗

Dynamic Magnetic Field Applications for Materials Processing

Magnetic fields, variable in time and space, can be used to control convection in electrically conducting melts. Flow induced by these fields has been found to be beneficial for crystal growth applications. It allows increased crystal growth rates, and improves homogeneity and quality. Particularly beneficial is the natural convection damping capability of alternating magnetic fields. One well-known example is the rotating magnetic field (RMF) configuration. RMF induces liquid motion consisting of a swirling basic flow and a meridional secondary flow. In addition to crystal growth applications, RMF can also be used for mixing non-homogeneous melts in continuous metal castings. These applied aspects have stimulated increasing research on RMF-induced fluid dynamics. A novel type of magnetic field configuration consisting of an axisymmetric magnetostatic wave, designated the traveling magnetic field (TMF), has been recently proposed. It induces a basic flow in the form of a single vortex. TMF may find use in crystal growth techniques such as the vertical Bridgman (VB), float zone (FZ), and the traveling heater method. In this review, both methods, RMF and TMF are presented. Our recent theoretical and experimental results include such topics as localized TMF, natural convection dumping using TMF in a vertical Bridgman configuration, the traveling heater method, and the Lorentz force induced by TMF as a function of frequency. Experimentally, alloy mixing results, with and without applied TMF, will be presented. Finally, advantages of the traveling magnetic field, in comparison to the more mature rotating magnetic field method, will be discussed.

Mazuruk, K.↗

Aeroelastic System Development Using Proper Orthogonal Decomposition and Volterra Theory

This research combines Volterra theory and proper orthogonal decomposition (POD) into a hybrid methodology for reduced-order modeling of aeroelastic systems. The out-come of the method is a set of linear ordinary differential equations (ODEs) describing the modal amplitudes associated with both the structural modes and the POD basis functions for the uid. For this research, the structural modes are sine waves of varying frequency, and the Volterra-POD approach is applied to the fluid dynamics equations. The structural modes are treated as forcing terms which are impulsed as part of the uid model realization. Using this approach, structural and uid operators are coupled into a single aeroelastic operator. This coupling converts a free boundary uid problem into an initial value problem, while preserving the parameter (or parameters) of interest for sensitivity analysis. The approach is applied to an elastic panel in supersonic cross ow. The hybrid Volterra-POD approach provides a low-order uid model in state-space form. The linear uid model is tightly coupled with a nonlinear panel model using an implicit integration scheme. The resulting aeroelastic model provides correct limit-cycle oscillation prediction over a wide range of panel dynamic pressure values. Time integration of the reduced-order aeroelastic model is four orders of magnitude faster than the high-order solution procedure developed for this research using traditional uid and structural solvers.

Lucia, David J.↗

Application of Micro-ramp Flow Control Devices to an Oblique Shock Interaction

Tests are planned in the 15cm x 15cm supersonic wind tunnel at NASA Glenn to demonstrate the applicability of micro-ramp flow control to the management of shock wave boundary layer interactions. These tests will be used as a database for computational fluid dynamics (CFD) validation and Design of Experiments (DoE) design information. Micro-ramps show potential for mechanically simple and fail-safe boundary layer control.

Hirt, Stefanie↗

Plume and Shock Interaction Effects on Sonic Boom in the 1-foot by 1-foot Supersonic Wind Tunnel

The desire to reduce or eliminate the operational restrictions of supersonic aircraft over populated areas has led to extensive research at NASA. Restrictions are due to the disturbance of the sonic boom, caused by the coalescence of shock waves formed by the aircraft. A study has been performed focused on reducing the magnitude of the sonic boom N-wave generated by airplane components with a focus on shock waves caused by the exhaust nozzle plume. Testing was completed in the 1-foot by 1-foot supersonic wind tunnel to study the effects of an exhaust nozzle plume and shock wave interaction. The plume and shock interaction study was developed to collect data for computational fluid dynamics (CFD) validation of a nozzle plume passing through the shock generated from the wing or tail of a supersonic vehicle. The wing or tail was simulated with a wedgeshaped shock generator. This test entry was the first of two phases to collect schlieren images and off-body static pressure profiles. Three wedge configurations were tested consisting of strut-mounted wedges of 2.5- degrees and 5-degrees. Three propulsion configurations were tested simulating the propulsion pod and aft deck from a low boom vehicle concept, which also provided a trailing edge shock and plume interaction. Findings include how the interaction of the jet plume caused a thickening of the shock generated by the wedge (or aft deck) and demonstrate how the shock location moved with increasing nozzle pressure ratio.

Shock wave interaction↗

Baffle Effect on Sloshing-Modulated Torques Responsed to Orbital Accelerations in Microgravity

The behavior of sloshing-dynamics-modulated fluid systems driven by orbital accelerations including gravity gradient and jitter accelerations has been studied. Partly liquid-filled rotating Dewars applicable to a full-scale Gravity Probe-B Spacecraft container with and without baffle are considered. Results show that a slosh wave excited along the liquid-vapor interface induced by jitter-acceleration-dominated orbital accelerations provides a torsional moment with an up-and-down movement of bubble oscillations in the rotating Dewar. The fluctuations of fluid forces and fluid moment (torque) exerted on the Dewar wall of the container caused by slosh-wave excitation driven by orbital accelerations are also investigated. Since the viscous force across a liquid-solid interface and the surface-tension force across a liquid-vapor or liquid-solid interface can greatly contribute to the damping effect of slosh-wave excitation, a rotating Dewar with a baffle provides more areas of liquid-solid, and liquid-vapor, and solid-vapor interfaces than does a rotating Dewar without a baffle. Results show that the damping effect provided by the baffle reduces the amplitudes of force and moment feedback from the fluids to the container-in particular, the components of fluctuations transverse to the direction of the baffle boards.

Pan, H. L.↗

Hydrodynamics of coalescing binary neutron stars: Ellipsoidal treatment

We employ an approximate treatment of dissipative hydrodynamics in three dimensions to study the coalescence of binary neutron stars driven by the emission of gravitational waves. The stars are modeled as compressible ellipsoids obeying a polytropic equation of state; all internal fluid velocities are assumed to be linear functions of the coordinates. The hydrodynamics equations then reduce to a set of coupled ordinary differential equations for the evolution of the principal axes of the ellipsoids, the internal velocity parameters, and the binary orbital parameters. Gravitational radiation reaction and viscous dissipation are both incorporated. We set up exact initial binary equilibrium configurations and follow the transition from the quasi-static, secular decay of the orbit at large separation to the rapid dynamical evolution of the configurations just prior to contact. A hydrodynamical instability resulting from tidal interactions significantly accelerates the coalescence at small separation, leading to appreciable radial infall velocity and tidal lag angles near contact. This behavior is reflected in the gravitational waveforms and may be observable by gravitational wave detectors under construction. In cases where the neutron stars have spins which are not aligned with the orbital angular momentum, the spin-induced quadrupole moment can lead to precession of the orbital plane and therefore modulation of the gravitational wave amplitude even at large orbital radius. However, the amplitude of the modulation is small for typical neutron star binaries with spins much smaller than the orbital angular momentum.

Lai, Dong↗

Propagation of acoustic modes in the transitional ionosphere

Continuum theory is unsatisfactory for describing wave propagation in the transitional ionosphere, since the mean free path is of the same order as the characteristic length (i.e., Knudsen number in this region is of the order of unity). Therefore a transition model in which isotropic electrons behave as a fluid while the ions are governed by the kinetic equation modified by ion-electron collision effects is proposed to study acoustic wave propagation in the upper atmosphere in the altitude region from 500 to 2000 km. The results show that the dissipation of acoustic waves by electrons is through the viscous and thermal conduction effects, and the dissipation of acoustic waves by ions is through Landau damping. Comparisons of the characteristics of acoustic mode propagation and wave dissipation mechanisms in the collision-dominated, transitional from collisional to collisionless, and collisionless media are also discussed.

Hung, R. J.↗

Sloshing dynamics modulated fluid angular momentum and moment fluctuations driven by orbital gravity gradient and jitter accelerations in microgravity

The dynamical behavior of spacecraft propellant affected by the asymmetric combined gravity gradient and jitter accelerations, in particular the effect of surface tension on partially-filled rotating fluids applicable to a full-scale Gravity Probe-B Spacecraft dewar tank has been investigated. Three different cases of orbital accelerations: (1) gravity gradient-dominated, (2) equally weighted between gravity gradient and jitter, and (3) gravity jitter-dominated accelerations are studied. The results of slosh wave excitation along the liquid-vapor interface induced by gravity gradient-dominated accelerations provide a torsional moment with tidal motion of bubble oscillations in the rotating dewar. The results are clearly seen from the twisting shape of the bubble oscillations driven by gravity gradient-dominated acceleration. The results of slosh wave excitation along the liquid-vapor interface induced by gravity jitter-dominated acceleration indicate the results of bubble motion in a manner of down-and-up and leftward-and-rightward movement of oscillation when the bubble is rotating with respect to rotating dewar axis. Fluctuations of angular momentum, fluid moment and bubble mass center caused by slosh wave excitations driven by gravity gradient acceleration or gravity jitter acceleration are also investigated.

Hung, R. J.↗

Computational Assessment of the Impact of Wave Count on Rotating Detonation Engine Performance

The impact of the number of azimuthally propagating waves on the performance of a Rotating Detonation Engine (RDE) is investigated using a simplified two-dimensional computational fluid dynamic simulation. The basic RDE configuration examined has no exit throat. The inlet is assumed lossless and does not allow backflow. The adiabatic, inviscid, and premixed simulation utilizes a particularly simple finite rate reaction mechanism that allows user control over the relative amounts of deflagration and detonation that occur, the presence and extent of a reaction delay associated with fuel and oxidizer mixing, and the number of waves present in the domain. The simplifications and idealizations decouple the effect of wave count from other potential loss mechanisms. Performance is measured using exhaust flow total pressure gain relative to the inlet total pressure. One, two, and three wave solutions are computed under identical boundary conditions and grid resolution. For simulations with no mixing delay and minimal deflagration, the number of waves present has negligible impact on performance. With a mixing delay, performance decreases with increasing wave number. With increased deflagration, performance increases with increasing wave number. With both effects simulated it is found that the two wave solution performs better than either the one or three wave solutions. The causes of these trends are explored. They imply that for practical, RDE’s, the number of waves that are present can impact performance.

Detonation↗

Computational Assessment of the Impact of Wave Count on Rotating Detonation Engine Performance

The impact of the number of azimuthally propagating waves on the performance of a Rotating Detonation Engine (RDE) is investigated using a simplified two-dimensional computational fluid dynamic simulation. The basic RDE configuration examined has no exit throat. The inlet is assumed lossless and does not allow backflow. The adiabatic, inviscid, and premixed simulation utilizes a particularly simple finite rate reaction mechanism that allows user control over the relative amounts of deflagration and detonation that occur, the presence and extent of a reaction delay associated with fuel and oxidizer mixing, and the number of waves present in the domain. The simplifications and idealizations decouple the effect of wave count from other potential loss mechanisms. Performance is measured using exhaust flow total pressure gain relative to the inlet total pressure. One, two, and three wave solutions are computed under identical boundary conditions and grid resolution. For simulations with no mixing delay and minimal deflagration, the number of waves present has negligible impact on performance. With a mixing delay, performance decreases with increasing wave number. With increased deflagration, performance increases with increasing wave number. With both effects simulated it is found that the two wave solution performs better than either the one or three wave solutions. The causes of these trends are explored. They imply that for practical, RDE’s, the number of waves that are present can impact performance.

Detonation↗

The effect of a coronal shock wave on the solar wind ionization state

In connection with studies of solar wind ionization state freezing, it is shown that, by using a Lagrangian approach of following individual fluid parcels, the techniques used previously for calculating ionization state variations in a steady state case can be extended straightforwardly to time-varying flows. The specific ionization state calculations presented are for a relatively simple picture of time-dependent coronal flow, based on a well-known model of a self-similar shock wave propagating through the corona. Time-dependent ionization effects for the sudden transition between two otherwise steady flows are likely to be limited to a narrow range of gas parcels which, having been shocked within the coronal freezing-in radius, pass a fixed interplanetary observer in an interval of a few tens of minutes. The amplitude of any rise in interplanetary ionization temperature associated with the coronal shock is likely to be considerably smaller than the jump in electron temperature that actually occurs in the corona.

Owocki, S. P.↗

Flow tracing fidelity of scattering aerosol in laser Doppler velocimetry

An experimental method for determinating the flow tracing fidelity of a scattering aerosol used in laser Doppler velocimeters was developed with particular reference to the subsonic turbulence measurements. The method employs the measurement of the dynamic response of a flow seeding aerosol excited by acoustic waves. The amplitude and frequency of excitation were controlled to simulate the corresponding values of fluid turbulence components. Experimental results are presented on the dynamic response of aerosols over the size range from 0.1 to 2.0 microns in diameter and over the frequency range 100 Hz to 100 kHz. It was observed that unit density spherical scatterers with diameters of 0.2 microns followed subsonic air turbulence frequency components up to 100 kHz with 98 percent fidelity.

Mazumder, M. K.↗

Skylab fluid mechanics demostration - A study of water droplet oscillations in space

In the Skylab experiment a drop of water was caused to oscillate at its natural oscillation frequency. The drop was observed until its oscillation began to decay due to its internal damping. A determination of the change in amplitude with time made it possible to verify the applicability of a theoretical model. The fundamental concept of the model is that fluid surfaces tend to an equilibrium shape produced by the balance of the forces of fluid pressure and surface tension. The theoretical model is based on the assumption that the wave amplitude is small compared with the wavelength.

Vaughan, O. H.↗