Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid instabilities”

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 523 records · Page 29

Fluid Physical and Transport Phenomena Studies aboard the International Space Station: Planned Experiments

This paper provides an overview of the microgravity fluid physics and transport phenomena experiments planned for the International Spare Station. NASA's Office of Life and Microgravity Science and Applications has established a world-class research program in fluid physics and transport phenomena. This program combines the vast expertise of the world research community with NASA's unique microgravity facilities with the objectives of gaining new insight into fluid phenomena by removing the confounding effect of gravity. Due to its criticality to many terrestrial and space-based processes and phenomena, fluid physics and transport phenomena play a central role in the NASA's Microgravity Program. Through widely publicized research announcement and well established peer-reviews, the program has been able to attract a number of world-class researchers and acquired a critical mass of investigations that is now adding rapidly to this field. Currently there arc a total of 106 ground-based and 20 candidate flight principal investigators conducting research in four major thrust areas in the program: complex flows, multiphase flow and phase change, interfacial phenomena, and dynamics and instabilities. The International Space Station (ISS) to be launched in 1998, provides the microgravity research community with a unprecedented opportunity to conduct long-duration microgravity experiments which can be controlled and operated from the Principal Investigators' own laboratory. Frequent planned shuttle flights to the Station will provide opportunities to conduct many more experiments than were previously possible. NASA Lewis Research Center is in the process of designing a Fluids and Combustion Facility (FCF) to be located in the Laboratory Module of the ISS that will not only accommodate multiple users but, allow a broad range of fluid physics and transport phenomena experiments to be conducted in a cost effective manner.

Singh, Bhim S.↗

A Technique for Rapidly Deploying a Concentration Gradient with Applications to Microgravity

The latter half of the last century has seen rapid advancements in semiconductor crystal growth powered by the demand for high performance electronics in myriad applications. The reduced gravity environment of space has also been used for crystal growth tests, especially in instances where terrestrial growth has largely been unsuccessful. While reduced gravity crystal growth affords some control of the gravity parameter, many crystals grown in space, to date, have structural flaws believed to result from convective motions during the growth phase. The character of these instabilities is not well understood but is associated with thermal and solutal density variations near the solidification interface in the presence of residual gravity and g-jitter. In order to study these instabilities in a separate, controlled space experiment, a concentration gradient would first have to be artificially established in a timely manner as an initial condition. This is generally difficult to accomplish in a microgravity environment because the momentum of the fluid injected into a test cell tends to swirl around and mix in the absence of a restoring force. The use of magnetic fields to control the motion and position of liquids has received growing interest in recent times. The possibility of using the force exerted by a non-uniform magnetic field on a ferrofluid to not only achieve fluid manipulation but also to actively control fluid motion makes it an attractive candidate for space applications. This paper describes a technique for quickly establishing a linear or exponential fluid concentration gradient using a magnetic field in place of gravity to stabilize the deployment. Also discussed is a photometric technique for measuring the concentration profile using light attenuation. Results of the ground-based experiments indicate that the concentration distribution is within 3% of the predicted value. Although any range of concentations can be realized, photometric constraints are discussed which impose some limitations on measurements.

Leslie, Fred↗

The Use of Ferrofluids to Model Materials Processing (MSFC Center Director's Discretionary Fund)

Many crystals grown in space have structural flaws believed to result from convective motions during the growth phase. The character of these instabilities is not well understood but is associated with thermal and solutal density variations near the solidification interface in the presence of residual gravity and g-jitter. To study these instabilities in a separate, controlled space experiment, a concentration gradient would first have to be artificially established in a timely manner as an initial condition. This is generally difficult to accomplish in a microgravity environment because the momentum of the fluid injected into a test cell tends to swirl around and mix in the absence of a restoring force. The use of magnetic fields to control the motion and position of liquids has received recent, growing interest. The possibility of using the force exerted by a non-uniform magnetic field on a ferrofluid to not only achieve fluid manipulation but also to actively control fluid motion makes it an attractive candidate for space applications. This paper describes a technique for quickly establishing a linear or exponential fluid concentration gradient using a magnetic field in place of gravity to stabilize the deployment. Also discussed is a photometric technique for measuring the concentration profile using light attenuation. Although any range of concentrations can be realized, photometric constraints impose some limitations on measurements. Results of the ground-based experiments indicate that the species distribution is within 3 percent of the predicted value.

Leslie, F.↗

Numerical Simulation of Protoplanetary Vortices

The fluid dynamics within a protoplanetary disk has been attracting the attention of many researchers for a few decades. Previous works include, to list only a few among many others, the well-known prescription of Shakura & Sunyaev, the convective and instability study of Stone & Balbus and Hawley et al., the Rossby wave approach of Lovelace et al., as well as a recent work by Klahr & Bodenheimer, which attempted to identify turbulent flow within the disk. The disk is commonly understood to be a thin gas disk rotating around a central star with differential rotation (the Keplerian velocity), and the central quest remains as how the flow behavior deviates (albeit by a small amount) from a strong balance established between gravitational and centrifugal forces, transfers mass and momentum inward, and eventually forms planetesimals and planets. In earlier works we have briefly described the possible physical processes involved in the disk; we have proposed the existence of long-lasting, coherent vortices as an efficient agent for mass and momentum transport. In particular, Barranco et al. provided a general mathematical framework that is suitable for the asymptotic regime of the disk; Barranco & Marcus (2000) addressed a proposed vortex-dust interaction mechanism which might lead to planetesimal formation; and Lin et al. (2002), as inspired by general geophysical vortex dynamics, proposed basic mechanisms by which vortices can transport mass and angular momentum. The current work follows up on our previous effort. We shall focus on the detailed numerical implementation of our problem. We have developed a parallel, pseudo-spectral code to simulate the full three-dimensional vortex dynamics in a stably-stratified, differentially rotating frame, which represents the environment of the disk. Our simulation is validated with full diagnostics and comparisons, and we present our results on a family of three-dimensional, coherent equilibrium vortices.

Lin, H.↗

Flow instability in particle-bed nuclear reactors

The particle-bed core offers mitigation of some of the problems of solid-core nuclear rocket reactors. Dividing the fuel elements into small spherical particles contained in a cylindrical bed through which the propellant flows radially, may reduce the thermal stress in the fuel elements, allowing higher propellant temperatures to be reached. The high temperature regions of the reactor are confined to the interior of cylindrical fuel assemblies, so most of the reactor can be relatively cool. This enables the use of structural and moderating materials which reduce the minimum critical size and mass of the reactor. One of the unresolved questions about this concept is whether the flow through the particle-bed will be well behaved, or will be subject to destructive flow instabilities. Most of the recent analyses of the stability of the particle-bed reactor have been extensions of the approach of Bussard and Delauer, where the bed is essentially treated as an array of parallel passages, so that the mass flow is continuous from inlet to outlet through any one passage. A more general three dimensional model of the bed is adopted, in which the fluid has mobility in three dimensions. Comparison of results of the earlier approach to the present one shows that the former does not accurately represent the stability at low Re. The more complete model presented should be capable of meeting this deficiency while accurately representing the effects of the cold and hot frits, and of heat conduction and radiation in the particle-bed. It can be extended to apply to the cylindrical geometry of particle-bed reactors without difficulty. From the exemplary calculations which were carried out, it can be concluded that a particle-bed without a cold frit would be subject to instability if operated at the high temperatures desired for nuclear rockets, and at power densities below about 4 megawatts per liter. Since the desired power density is about 40 megawatts per liter, it can be concluded that operation at design exit temperature but at reduced power could be hazardous for such a reactor. But the calculations also show that an appropriate cold frit could very likely cure the instability. More definite conclusions must await calculations for specific designs.

Kerrebrock, Jack L.↗

Time-dependent thermocapillary convection in a Cartesian cavity - Numerical results for a moderate Prandtl number fluid

The present numerical simulation explores a thermal-convective mechanism for oscillatory thermocapillary convection in a shallow Cartesian cavity for a Prandtl number 6.78 fluid. The computer program developed for this simulation integrates the two-dimensional, time-dependent Navier-Stokes equations and the energy equation by a time-accurate method on a stretched, staggered mesh. Flat free surfaces are assumed. The instability is shown to depend upon temporal coupling between large scale thermal structures within the flow field and the temperature sensitive free surface. A primary result of this study is the development of a stability diagram presenting the critical Marangoni number separating steady from the time-dependent flow states as a function of aspect ratio for the range of values between 2.3 and 3.8. Within this range, a minimum critical aspect ratio near 2.3 and a minimum critical Marangoni number near 20,000 are predicted below which steady convection is found.

Peltier, L. J.↗

Exact Solution to Stationary Onset of Convection Due to Surface Tension Variation in a Multicomponent Fluid Layer With Interfacial Deformation

Stationary onset of convection due to surface tension variation in an unbounded multicomponent fluid layer is considered. Surface deformation is included and general flux boundary conditions are imposed on the stratifying agencies (temperature/composition) disturbance equations. Exact solutions are obtained to the general N-component problem for both finite and infinitesimal wavenumbers. Long wavelength instability may coexist with a finite wavelength instability for certain sets of parameter values, often referred to as frontier points. For an impermeable/insulated upper boundary and a permeable/conductive lower boundary, frontier boundaries are computed in the space of Bond number, Bo, versus Crispation number, Cr, over the range 5 x 10(exp -7) less than or equal to Bo less than or equal to 1. The loci of frontier points in (Bo, Cr) space for different values of N, diffusivity ratios, and, Marangoni numbers, collapsed to a single curve in (Bo, D(dimensional variable)Cr) space, where D(dimensional variable) is a Marangoni number weighted diffusivity ratio.

Skarda, J. Raymond Lee↗

Flap Edge Aeroacoustic Measurements and Predictions

An aeroacoustic model test has been conducted to investigate the mechanisms of sound generation on high-lift wing configurations. This paper presents an analysis of flap side-edge noise, which is often the most dominant source. A model of a main element wing section with a half-span flap was tested at low speeds of up to a Mach number of 0.17, corresponding to a wing chord Reynolds number of approximately 1.7 million. Results are presented for flat (or blunt), flanged, and round flap-edge geometries, with and without boundary-layer tripping, deployed at both moderate and high flap angles. The acoustic database is obtained from a Small Aperture Directional Array (SADA) of microphones, which was constructed to electronically steer to different regions of the model and to obtain farfield noise spectra and directivity from these regions. The basic flap-edge aerodynamics is established by static surface pressure data, as well as by Computational Fluid Dynamics (CFD) calculations and simplified edge flow analyses. Distributions of unsteady pressure sensors over the flap allow the noise source regions to be defined and quantified via cross-spectral diagnostics using the SADA output. It is found that shear layer instability and related pressure scatter is the primary noise mechanism. For the flat edge flap, two noise prediction methods based on unsteady-surface-pressure measurements are evaluated and compared to measured noise. One is a new causality spectral approach developed here. The other is a new application of an edge-noise scatter prediction method. The good comparisons for both approaches suggest that much of the physics is captured by the prediction models. Areas of disagreement appear to reveal when the assumed edge noise mechanism does not fully define, the noise production. For the different edge conditions, extensive spectra and directivity are presented. Significantly, for each edge configuration, the spectra for different flow speeds, flap angles, and surface roughness were successfully scaled by utilizing aerodynamic performance and boundary layer scaling method developed herein.

Brooks, Thomas F.↗

Flap Edge Aeroacoustic Measurements and Predictions

An aeroacoustic model test has been conducted to investigate the mechanisms of sound generation on high-lift wing configurations. This paper presents an analysis of flap side-edge noise, which is often the most dominant source. A model of a main element wing section with a half-span flap was tested at low speeds of up to a Mach number of 0.17, corresponding to a wing chord Reynolds number of approximately 1.7 million. Results are presented for flat (or blunt), flanged, and round flap-edge geometries, with and without boundary-layer tripping, deployed at both moderate and high flap angles. The acoustic database is obtained from a Small Aperture Directional Array (SADA) of microphones, which was constructed to electronically steer to different regions of the model and to obtain farfield noise spectra and directivity from these regions. The basic flap-edge aerodynamics is established by static surface pressure data, as well as by Computational Fluid Dynamics (CFD) calculations and simplified edge flow analyses. Distributions of unsteady pressure sensors over the flap allow the noise source regions to be defined and quantified via cross-spectral diagnostics using the SADA output. It is found that shear layer instability and related pressure scatter is the primary noise mechanism. For the flat edge flap, two noise prediction methods based on unsteady surface pressure measurements are evaluated and compared to measured noise. One is a new causality spectral approach developed here. The other is a new application of an edge-noise scatter prediction method. The good comparisons for both approaches suggest that much of the physics is captured by the prediction models. Areas of disagreement appear to reveal when the assumed edge noise mechanism does not fully define the noise production. For the different edge conditions, extensive spectra and directivity are presented. Significantly, for each edge configuration, the spectra for different flow speeds, flap angles, and surface roughness were successfully scaled by utilizing aerodynamic performance and boundary layer scaling methods developed herein.

Brooks, Thomas F.↗

Physics Based Model for Cryogenic Chilldown and Loading. Part I: Algorithm

We report the progress in the development of the physics based model for cryogenic chilldown and loading. The chilldown and loading is model as fully separated non-equilibrium two-phase flow of cryogenic fluid thermally coupled to the pipe walls. The solution follow closely nearly-implicit and semi-implicit algorithms developed for autonomous control of thermal-hydraulic systems developed by Idaho National Laboratory. A special attention is paid to the treatment of instabilities. The model is applied to the analysis of chilldown in rapid loading system developed at NASA-Kennedy Space Center. The nontrivial characteristic feature of the analyzed chilldown regime is its active control by dump valves. The numerical predictions are in reasonable agreement with the experimental time traces. The obtained results pave the way to the development of autonomous loading operation on the ground and space.

heat transfer↗

On the initiation of melt fracture.

Discussion of a problem encountered in the processing of viscoelastic materials that is caused by the presence of a hydrodynamic instability in the extrusion of polymer melts. The importance of the so-called Weissenberg number in determining the onset of the melt fracture is examined using classical linearized hydrodynamic stability analysis. It is shown that the simple shearing flow of a viscoelastic fluid becomes unstable at a critical value of the Weissenberg number. Implications for the extrusion processing of polymers are reviewed.

Mcintire, L. V.↗

Incompressible fluid ellipsoids in halos. I - The second-harmonic oscillations of the Maclaurin spheroids

The structure and stability of Maclaurin spheroids embedded in rigid uniform-density oblate spheroidal halos are determined by the tensor virial-equation method. These spheroid-halo systems can be thought of as crude fluid analogs of disk galaxies with halos. The halos are assumed to have the same center, the same axis of symmetry, and the same equatorial radius as the Maclaurin spheroids. Only halos with lower eccentricity than the Maclaurin spheroids are considered. The dynamic instability of the toroidal (barlike) modes is suppressed when m, the ratio of the halo mass to Maclaurin spheroid mass, is greater than 3 pi/8 for spherical halos and when m is greater than 1/2 for halos congruent to the Maclaurin spheroids. Intermediate halo flattenings yield intermediate critical m-values. On the other hand, a neutral point of the toroidal modes in the rotating and inertial frames occurs for all m and for all allowed halo flattenings. Growth rates for secular instability beyond the neutral point are calculated, and the eigenfrequencies of all second-harmonic modes are given for select cases. The Ostriker-Peebles (1973) conjecture concerning the stability of disk galaxies against barlike perturbations appears to be incorrect.

Durisen, R. H.↗

Kinetic effects on Alfven wave nonlinearity. II - The modified nonlinear wave equation

A previously developed Vlasov theory is used here to study the role of resonant particle and other kinetic effects on Alfven wave nonlinearity. A hybrid fluid-Vlasov equation approach is used to obtain a modified version of the derivative nonlinear Schroedinger equation. The differences between a scalar model for the plasma pressure and a tensor model are discussed. The susceptibilty of the modified nonlinear wave equation to modulational instability is studied. The modulational instability normally associated with the derivative nonlinear Schroedinger equation will, under most circumstances, be restricted to left circularly polarized waves. The nonlocal term in the modified nonlinear wave equation engenders a new modulational instability that is independent of beta and the sense of circular polarization. This new instability may explain the occurrence of wave packet steepening for all values of the plasma beta in the vicinity of the earth's bow shock.

Spangler, Steven R.↗

Tuning Neural Network Models for Improved Prediction of Boundary Layer Transition

Boundary layer transition can strongly impact flight vehicle performance as it influences surface skin friction and aerodynamic heating, making accurate transition prediction a key to designing next generation aircraft. Artificial neural networks (ANNs) have shown promise toward predicting laminar-turbulent transition based on linear stability correlations. The computational efficiency of ANNs and the substantially reduced user involvement in relation to direct computations based on the linear stability theory (LST) makes them an attractive methodology for integrating the LST based correlations in computational fluid dynamics codes. Tollmien-Schlichting (TS) waves correspond to the dominant transition mechanism in 2D or weakly 3D subsonic boundary layers, such as those encountered in general aviation applications. Improvements to neural network model accuracy in predicting the amplification rates of TS instability waves have been investigated by leveraging recent machine learning developments in conjunction with surrogate optimization techniques and via suitable augmentation of the data used to train the networks. The optimized models trained on the modified dataset reduced the average transition location errors on different airfoils at several flow conditions by 51% of the original manually-tuned network’s errors on the same flow cases. The actual transition locations were derived from the Langley Stability and Transition Analysis Code (LASTRAC).

Machine Learning↗

Ellipsoidal figures of equilibrium - Compressible models

The results of Chandrasekhar (1969) are generalized to polytropes, using a formalism based on ellipsoidal energy variational principle to construct approximate stellar equilibrium solutions and study their stability. After reviewing the energy variational method and describing the approach, several equivalent stability conditions are established and secular vs. dynamical instabilities are discussed. Then, the equilibrium structure equations are derived for isolated, rotating polytropes, and axisymmetric configurations (compressible Maclaurin spheroids) are considered. Particular attention is given to triaxial configurations, either in a state of uniform rotation (generalizing the classical Jacobi ellipsoids) or with internal fluid motions of uniform vorticity (the compressible analogues of Riemann-S ellipsoids) and to the stability of these single star configurations. The compressible generalizations of the Roche and Roche-Riemann problems for a polytrope in orbit about a point-mass companion are solved, and the generalized Darwin problem for two identical polytropes in a binary is considered.

Lai, Dong↗

Coefficient matrices for implicit finite difference solution of the inviscid fluid conservation law equations

Although the Navier-Stokes equations describe most flows of interest in aerodynamics, the inviscid conservation law equations may be used for small regions with viscous forces. Thus, Euler equations and several time-accurate finite difference procedures, explicit and implicit, are discussed. Although implicit techniques require more computational work, they permit larger time steps to be taken without instability. It is noted that the Jacobian matrices for Euler equations in conservation-law form have certain eigenvalue-eigenvector properties which may be used to construct conservative-form coefficient matrices. This reduces the computation time of several implicit and semiimplicit schemes. Extensions of the basic approach to other areas are suggested.

Steger, J. L.↗

Microstructure of the solar wind

In this review microstructure is understood to summarize the dynamical properties of the solar wind that enter the fluid equations, usually employed to describe the overall flow, only in one or the other averaged form. Thus, first a description of the kinetic picture of the solar wind is given; it mainly involves the structure of the particle distributions in velocity space and associated instabilities. Waves and discontinuities are discussed next. Two selected examples are mentioned where the effects of waves and discontinuities play a role. Specific physical properties of the interplanetary magnetic irregularity spectrum and its relevance for cosmic ray propagation are discussed at the end.

Voelk, H. J.↗