Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dissipative 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 433 records · Page 24

Unusual Bow Shock and Magnetopause Locations During May 11-12, 1999

On May 11 and 12, 1999, when the solar wind proton density dropped to only about 3% of its nominal value, WIND and IMP 8 (and possibly Lunar Prospector) encountered the Earth's bow shock at unusually large distances. These observations, at 57-63 Re, correspond to the most distant crossings of the bow shock ever recorded. This paper will present observations from Interball-1 and Geotail, showing the bow shock near its nominal position, and from IMP 8 and WIND (and perhaps Lunar Prospector) crossing the bow shock at increasing radial distances. The bow shock is extremely weak by the time of its observation by WIND. However, a clear jump in both the magnetic field strength and solar wind bulk velocity is quite apparent, although with no significant changes in solar wind proton density. The nature of this dissipating shock will be discussed. Also, interball-1 has observed the magnetopause at an unusually large standoff distance. The observed motion of the magnetopause and bow shock will be compared with model predictions. This set of observations provide an excellent test case for the behavior of gas dynamic and MHD models under extreme conditions.

Szabo, Adam↗

A Self-Contained Mapping Closure Approximation for Scalar Mixing

Scalar turbulence exhibits interplays of coherent structures and random fluctuations over a broad range of spatial and temporal scales. This feature necessitates a probabilistic description of the scalar dynamics, which can be achieved comprehensively by using probability density functions (PDFs). Therefore, the challenge is to obtain the scalar PDFs (Lundgren 1967; Dopazo 1979). Generally, the evolution of a scalar is governed by three dynamical processes: advection, diffusion and reaction. In a PDF approach (Pope 1985), the advection and reaction can be treated exactly but the effect of molecular diffusion has to be modeled. It has been shown (Pope 1985) that the effect of molecular diffusion can be expressed as conditional dissipation rates or conditional diffusions. The currently used models for the conditional dissipation rates and conditional diffusions (Pope 1991) have resisted deduction from the fundamental equations and are unable to yield satisfactory results for the basic test cases of decaying scalars in isotropic turbulence, although they have achieved some success in a variety of individual cases. The recently developed mapping closure approach (Pope 1991; Chen, Chen & Kraichnan 1989; Kraichnan 1990; Klimenko & Pope 2003) provides a deductive method for conditional dissipation rates and conditional di usions, and the models obtained can successfully describe the shape relaxation of the scalar PDF from an initial double delta distribution to a Gaussian one. However, the mapping closure approach is not able to provide the rate at which the scalar evolves. The evolution rate has to be modeled. Therefore, the mapping closure approach is not closed. In this letter, we will address this problem.

He, Guo-Wei↗

Method Development for Experimental Characterization of Dynamic Strength of Aluminum Structures

Space Flight Hardware (SFH) experiences intense vibratory loading during flight, which only lasts a few minutes. When determining the appropriate size of these components to withstand such loading without becoming damaged, standard design practice is to assume that the peak dynamic loads are applied statically. In doing so, the resulting stress is compared against a material strength parameter obtained from a quasi-static experiment. Since the near-peak stresses are only experienced over a small fraction of time in reality, this approach leads to design conservatism that unnecessarily increases structural mass as well as the associated inefficiency and financial cost. In an effort to modernize engineering design standards to appropriately consider the higher practical strength of dynamically loaded structures, this overarching research project seeks to develop an experimental test procedure for quantifying the dynamic strength of metallic alloys as a function of excitation frequency. In the ideal case, the characterization test would include an in-situ method for monitoring the onset and progression of plastic deformation of the test specimen undergoing vibratory loading. These new tests are designed to be high intensity (forcing amplitude), short term (60 seconds at full amplitude), and cyclic in nature (sinusoidal excitation via attached stinger to an otherwise cantilevered beam). In addition, the initial alloy under investigation is 6061 aluminum, due to its wide use and applicability for SFH. Thus far, the primary candidate under development with live-capacity is to track hysteresis behavior of the beam from power dissipation trends, calculated via force (from transducer on excitation stinger) and velocity (measurements from Laser Doppler Vibrometer) data and work to distinguish between elastic and plastic features. As a key component of this ongoing method development, any pseudo-live indication of plasticity could be corroborated against the outcome of a pre-post assessment of damage via macroscopic evaluation of beam geometry (i.e. assessing any permanent change in the beam’s tip deflection). Support for experimental design decisions as well as dynamic strength data from tests with excitation frequencies of 10, 40, and 55 Hz will be discussed. This work contributes to the foundation for a new type of vibration-based characterization experiments and generates initial data on the functional strength of 6061 aluminum under the conditions considered.

vibrations, plastic deformation, aluminum, testing↗

Aeroacoustic diffraction and dissipation by a short propeller cowl in subsonic flight

This report develops and applies an aeroacoustic diffraction theory for a duct, or cowl, placed around modelled sources of propeller noise. The regime of flight speed is high subsonic. The modelled cowl's inner wall contains a liner with axially variable properties. Its exterior is rigid. The analysis replaces both sides with an unsteady lifting surface coupled to a dynamic thickness problem. The resulting pair of aeroacoustic governing equations for a lined 'ring wing' is valid both for a passive and for an active liner. Their numerical solution yields the effective dipole and monopole distributions of the shrouding system and thereby determines the cowl-diffracted component of the total radiated field. The sample calculations here include a preliminary parametric search for that liner layout which maximizes the cowl's shielding effectiveness. The main conclusion of the study is that a short cowl, passively lined, should provide moderate reductions in propeller noise.

Martinez, Rudolph↗

Characteristics of a Lunar Landing Configuration Having Various Multiple-Leg Landing-Gear Arrangements

An experimental investigation has been made of some lunar-landing characteristics of a 1/6-scale dynamic model of a landing module having multiple-leg landing-gear systems. Symmetric four-point and five-point systems were investigated. The landing-gear legs were inverted tripod arrangements having a telescoping main strut which incorporated a yielding-metal strap for energy dissipation, hinged V-struts, and circular pads. The landing tests were made by launching a free model onto an impenetrable hard surface (concrete) and onto a powdered-pumice overlay of various depths. Landing motion and acceleration data were obtained for a range of touchdown speeds, touchdown attitudes, and landing-surface conditions. Maximum normal acceleration experienced at the module center of gravity during landings on hard surface or pumice was 2g (full-scale lunar value in terms of earth's gravity) over a wide range of touchdown conditions. acceleration experienced was 12 1/2 radians/sec(exp 2) and maximum longitudinal acceleration was 1 3/4g. The module was very stable with all gear configurations during landings on hard surface (coefficient of friction, micron = 0.4) at all conditions tested. Some overturn instability occurred during landings on powdered pumice (micron = 0.7 to 1.0) depending upon flight path, pitch and yaw attitude, depth of pumice, surface topography, and landing-gear configuration. The effect on stability of roll attitude for the limited amount of roll-attitude landing data obtained was insignificant. Compared with the four-point landing gear, the five-point system with equal maximum gear radius increased landing stability slightly and improved the static stability for subsequent lunar launch. A considerable increase in landing stability in the direction of motion was obtained with an asymmetric four-point gear having two pads offset to increase gear radius by 33 percent in the direction of horizontal flight.

Blanchard, Ulysse J.↗

An experimental investigation of gas jets in confined swirling air flow

The fluid dynamics of jets in confined swirling flows which is of importance to designers of turbine combustors and solid fuel ramjets used to power missiles fired from cannons were examined. The fluid dynamics of gas jets of different densities in confined swirling flows were investigated. Mean velocity and turbulence measurements are made with a one color, one component laser velocimeter operating in the forward scatter mode. It is shown that jets in confined flow with large area ratio are highly dissipative which results in both air and helium/air jet centerline velocity decays. For air jets, the jet like behavior in the tube center disappears at about 20 diameters downstream of the jet exit. This phenomenon is independent of the initial jet velocity. The turbulence field at this point also decays to that of the background swirling flow. A jet like behavior in the tube center is noticed even at 40 diameters for the helium/air jets. The subsequent flow and turbulence field depend highly on the initial jet velocity. The jets are fully turbulent, and the cause of this difference in behavior is attributed to the combined action swirl and density difference. This observation can have significant impact on the design of turbine combustors and solid fuel ramjets subject to spin.

Mongia, H.↗

A study of numerical methods for hyperbolic conservation laws with stiff source terms

In the present study of the behavior of typical numerical methods in the case of a model advection equation having a parameter-dependent source term, two approaches to the incorporation of the source terms are used: MacCormack-type predictor-corrector methods with flux limiters, and splitting methods in which the fluid dynamics and chemistry are handled in separate steps. The latter are found to perform slightly better. The model scalar equation is used to show that the incorrectness of the propagation speeds of discontinuities observed in the stiff case is due to the introduction of nonequilibrium values through numerical dissipation in the advection step.

Leveque, R. J.↗

Method Development for Experimental Characterization of Dynamic Strength of Aluminum Structures

Space Flight Hardware (SFH) experiences intense vibratory loading during flight, which only lasts a few minutes. When determining the appropriate size of these components to withstand such loading without becoming damaged (plastically deforming or breaking), standard design practice is to assume that the peak dynamic loads are applied statically. In doing so, the resulting stress is compared against a material strength parameter obtained from a quasi-static experiment. Since the near-peak stresses are only experienced over a small fraction of time in reality, this approach leads to design conservatism that unnecessarily increases structural mass as well as the associated inefficiency and financial cost. In an effort to modernize engineering design standards to appropriately consider the higher practical strength of dynamically loaded structures, this overarching research project seeks to develop an experimental test procedure for quantifying the dynamic strength of metallic alloys as a function of excitation frequency. In the ideal case, the characterization test would include an in-situ method for monitoring the onset and progression of plastic deformation of the test specimen undergoing vibratory loading. These new tests are designed to be high intensity (forcing amplitude), short term (60 seconds at full amplitude), and cyclic in nature (sinusoidal excitation via attached stinger to an otherwise cantilevered beam). In addition, the initial metallic alloy under investigation is 6061 aluminum, due to its wide use and applicability for SFH. Thus far, the most promising monitoring method under development with live-capacity is to track hysteresis behavior of the beam from power dissipation trends, calculated via force (from transducer on excitation stinger) and velocity (measurements from Laser Doppler Vibrometer) data and work to distinguish between elastic and plastic features. As a key component of this ongoing method development, any pseudo-live indication of plasticity can be corroborated against the outcome of a pre-post assessment of damage via macroscopic evaluation of beam geometry (i.e. assessing any permanent change in the beam’s tip deflection). Support for experimental design decisions as well as dynamic strength data from tests with excitation frequencies of 10, 40, and 55 Hz will be discussed. This work contributes to the foundation for a new type of vibration-based characterization experiments and generates initial data on the functional strength of 6061 aluminum under the conditions considered. Keywords: vibrations, plastic deformation, aluminum

vibrations↗

Measurements of Acoustic Properties of Porous and Granular Materials and Application to Vibration Control

For application of porous and granular materials to vibro-acoustic controls, a finite dynamic strength of the solid component (frame) is an important design factor. The primary goal of this study was to investigate structural vibration damping through this frame wave propagation for various poroelastic materials. A measurement method to investigate the vibration characteristics of the frame was proposed. The measured properties were found to follow closely the characteristics of the viscoelastic materials - the dynamic modulus increased with frequency and the degree of the frequency dependence was determined by its loss factor. The dynamic stiffness of hollow cylindrical beams containing porous and granular materials as damping treatment was measured also. The data were used to extract the damping materials characteristics using the Rayleigh-Ritz method. The results suggested that the acoustic structure interaction between the frame and the structure enhances the dissipation of the vibration energy significantly.

Park, Junhong↗

A CFD Approach to Modeling Spacecraft Fuel Slosh

Energy dissipation and resonant coupling from sloshing fuel in spacecraft fuel tanks is a problem that occurs in the design of many spacecraft. In the case of a spin stabilized spacecraft, this energy dissipation can cause a growth in the spacecrafts' nutation (wobble) that may lead to disastrous consequences for the mission. Even in non-spinning spacecraft, coupling between the spacecraft or upper stage flight control system and an unanticipated slosh resonance can result in catastrophe. By using a Computational Fluid Dynamics (CFD) solver such as Fluent, a model for this fuel slosh can be created. The accuracy of the model must be tested by comparing its results to an experimental test case. Such a model will allow for the variation of many different parameters such as fluid viscosity and gravitational field, yielding a deeper understanding of spacecraft slosh dynamics. In order to gain a better understanding of the dynamics behind sloshing fluids, the Launch Services Program (LSP) at the NASA Kennedy Space Center (KSC) is interested in finding ways to better model this behavior. Thanks to past research, a state-of-the-art fuel slosh research facility was designed and fabricated at Embry Riddle Aeronautical University (ERAU). This test facility has produced interesting results and a fairly reliable parameter estimation process to predict the necessary values that accurately characterize a mechanical pendulum analog model. The current study at ERAU uses a different approach to model the free surface sloshing of liquid in a spherical tank using Computational Fluid Dynamics (CFD) methods. Using a software package called Fluent, a model was created to simulate the sloshing motion of the propellant. This finite volume program uses a technique called the Volume of Fluid (VOF) method to model the interaction between two fluids [4]. For the case of free surface slosh, the two fluids are the propellant and air. As the fuel sloshes around in the tank, it naturally displaces the air. Using the conservation of mass, momentum, and energy equations, as well as the VOF equations, one can predict the behavior of the sloshing fluid and calculate the forces, pressure gradients, and velocity field for the entire liquid as a function of time.

Marsell, Brandon↗

Physics-Based Approach to Predict the Solar Activity Cycles

Observations of the complex highly non-linear dynamics of global turbulent flows and magnetic fields are currently available only from Earth-side observations. Recent progress in helioseismology has provided us some additional information about the subsurface dynamics, but its relation to the magnetic field evolution is not yet understood. These limitations cause uncertainties that are difficult take into account, and perform proper calibration of dynamo models. The current dynamo models have also uncertainties due to the complicated turbulent physics of magnetic field generation, transport and dissipation. Because of the uncertainties in both observations and theory, the data assimilation approach is natural way for the solar cycle prediction and estimating uncertainties of this prediction. I will discuss the prediction results for the upcoming Solar Cycle 25 and their uncertainties and affect of Ensemble Kalman Filter parameters to resulting predictions.

Physics-Based↗

Machine-Learning-Based Adaptive Thinning of CrIS Radiances to Improve Global Tropical Cyclone Analysis and Forecasts

This work is focused on optimizing the assimilation of hyperspectral infrared (IR) radiances from the Cross-track Infrared Sounder (CrIS) with the goal of improving the representation of tropical cyclones (TCs) in global analyses and forecasts. Current operational assimilation systems rely on subsampling IR radiances on a regular thinning grid. A new and improved adaptive methodology based on machine learning (ML) recognizes TCs from geostationary satellite imagery and is implemented in the Goddard Earth Observing System (GEOS) model and data assimilation framework. The ML methodology is extensively trained on existing TC data sets and creates for each TC a dynamic mask, based on the evolving shape and life cycle of that specific event. Once a TC mask is created, a switch is then activated in the data assimilation system to alter the thinning, ingesting more CrIS radiances within the moving mask, thus increasing the TC sampling. After the TC dissipates, the assimilation of CrIS radiances reverts to normal data density. Results of TC segmentation provided by a state-of-the-art generative machine learning model known as the Denoising Diffusion Probabilistic Model (DDPM) are compared to the previously used U-Net model. The new approach surpasses the performance of the previously developed one. The methodology is applied to both clear-sky and cloud-cleared radiances. Benefits from the latter methodology, particularly in improving the structure of TCs and the intensity forecasts, are presented.

Oreste Reale↗

Origin and dynamics of emission line clouds in cooling flow environments

The author suggests that since clouds are born co-moving in a turbulent intra-cluster medium (ICM), the allowed parameter space can now be opened up to a more acceptable range. Large-scale motions can be driven in the central parts of cooling flows by a number of mechanisms including the motion of the central and other galaxies, and the dissipation of advected, focussed rotational and magnetic energy. In addition to the velocity width paradox, two other paradoxes (Heckman et al. 1989) can be solved if the ICM is turbulent. Firstly, the heating source for the emission line regions has always been puzzling - line luminosities are extremely high for a given (optical or radio) galaxy luminosity compared to those in non-cooling flow galaxies, therefore a mechanism peculiar to cooling flows must be at work. However most, if not all, previously suggested heating mechanisms either fail to provide enough ionization or give the wrong line ratios, or both. The kinetic energy in the turbulence provides a natural energy source if it can be efficiently converted to cloud heat. Researchers suggest that this can be done via magneto-hydrodynamic waves through plasma slip. Secondly, while the x ray observations indicate extended mass deposition, the optical line emission is more centrally concentrated. Since many of the turbulence-inducing mechanisms are strongest in the central regions of the ICM, so is the method of heating. In other words material is dropping out everywhere but only being lit up in the center.

Loewenstein, Michael↗

Polar thermospheric Joule heating, and redistribution of recombination energy in the upper mesosphere

Kellogg (1961), suggested that transport of atomic oxygen from the summer into the winter hemisphere and subsequent release of energy by three body recombination, O+O+N2 yields O2+N2+E, may contribute significantly to the so-called mesopause temperature anomaly. Earlier model calculations have shown that Kellogg's mechanism produces about a 10-percent increase in the temperature from summer to winter at 90 km. This process, however, is partly compensated by differential heating from absorption of UV radiation associated with dissociation of O2. In the auroral region of the thermosphere, there is a steady energy dissipation by Joule heating causing a redistribution and depletion of atomic oxygen due to wind-induced diffusion. With the removal of O, latent chemical energy normally released by three body recombination is also removed, and the result is that the temperature decreases by almost 2 percent near 90 km. Through dynamic feedback, this process reduces the depletion of atomic oxygen by about 25 percent and the temperature perturbation in the exosphere from 10 to 7 percent at polar latitudes. Under the influence of the internal dynamo interaction, the prevailing zonal circulation in the upper thermosphere changes direction when the redistribution of recombination energy is considered.

Mayr, H. G.↗

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗

MHD Turbulence and Magnetic Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V↗

Bulk-Flow Analysis of Hybrid Thrust Bearings for Advanced Cryogenic Turbopumps

A bulk-flow analysis and computer program for prediction of the static load performance and dynamic force coefficients of angled injection, orifice-compensated hydrostatic/hydrodynamic thrust bearings have been completed. The product of the research is an efficient computational tool for the design of high-speed thrust bearings for cryogenic fluid turbopumps. The study addresses the needs of a growing technology that requires of reliable fluid film bearings to provide the maximum operating life with optimum controllable rotordynamic characteristics at the lowest cost. The motion of a cryogenic fluid on the thin film lands of a thrust bearing is governed by a set of bulk-flow mass and momentum conservation and energy transport equations. Mass flow conservation and a simple model for momentum transport within the hydrostatic bearing recesses are also accounted for. The bulk-flow model includes flow turbulence with fluid inertia advection, Coriolis and centrifugal acceleration effects on the bearing recesses and film lands. The cryogenic fluid properties are obtained from realistic thermophysical equations of state. Turbulent bulk-flow shear parameters are based on Hirs' model with Moody's friction factor equations allowing a simple simulation for machined bearing surface roughness. A perturbation analysis leads to zeroth-order nonlinear equations governing the fluid flow for the thrust bearing operating at a static equilibrium position, and first-order linear equations describing the perturbed fluid flow for small amplitude shaft motions in the axial direction. Numerical solution to the zeroth-order flow field equations renders the bearing flow rate, thrust load, drag torque and power dissipation. Solution to the first-order equations determines the axial stiffness, damping and inertia force coefficients. The computational method uses well established algorithms and generic subprograms available from prior developments. The Fortran9O computer program hydrothrust runs on a Windows 95/NT personal computer. The program, help files and examples are licensed by Texas A&M University Technology License Office. The study of the static and dynamic performance of two hydrostatic/hydrodynamic bearings demonstrates the importance of centrifugal and advection fluid inertia effects for operation at high rotational speeds. The first example considers a conceptual hydrostatic thrust bearing for an advanced liquid hydrogen turbopump operating at 170,000 rpm. The large axial stiffness and damping coefficients of the bearing should provide accurate control and axial positioning of the turbopump and also allow for unshrouded impellers, therefore increasing the overall pump efficiency. The second bearing uses a refrigerant R134a, and its application in oil-free air conditioning compressors is of great technological importance and commercial value. The computed predictions reveal that the LH2 bearing load capacity and flow rate increase with the recess pressure (i.e. increasing orifice diameters). The bearing axial stiffness has a maximum for a recess pressure rati of approx. 0.55. while the axial damping coefficient decreases as the recess pressure ratio increases. The computer results from three flow models are compared. These models are a) inertialess, b) fluid inertia at recess edges only, and c) full fluid inertia at both recess edges and film lands. The full inertia model shows the lowest flow rates, axial load capacity and stiffness coefficient but on the other hand renders the largest damping coefficients and inertia coefficients. The most important findings are related to the reduction of the outflow through the inner radius and the appearance of subambient pressures. The performance of the refrigerant hybrid thrust bearing is evaluated at two operating speeds and pressure drops. The computed results are presented in dimensionless form to evidence consistent trends in the bearing performance characteristics. As the applied axial load increases, the bearing film thickness and flow rate decrease while the recess pressure increases. The axial stiffness coefficient shows a maximum for a certain intermediate load while the damping coefficient steadily increases. The computed results evidence the paramount of centrifugal fluid inertia at low recess pressures (i.e. low loads), and where there is actually an inflow through the bearing inner diameter, accompanied by subambient pressures just downstream of the bearing recess edge. These results are solely due to centrifugal fluid inertia and advection transport effects. Recommendations include the extension of the computer program to handle flexure pivot tilting pad hybrid bearings and the ability to calculate moment coefficients for shaft angular misalignments.

SanAndres, Luis↗

A New View of the Dynamics of Reynolds Stress Generation in Turbulent Boundary Layers

The structure of a numerically simulated turbulent boundary layer over a flat plate at Re(theta) = 670 was studied using the invariants of the velocity gradient tensor (Q and R) and a related scalar quantity, the cubic discriminant (D = 27R(exp 2)/4 + Q(exp 3)). These invariants have previously been used to study the properties of the small-scale motions responsible for the dissipation of turbulent kinetic energy. In addition, these scalar quantities allow the local flow patterns to be unambiguously classified according to the terminology proposed by Chong et al. The use of the discriminant as a marker of coherent motions reveals complex, large-scale flow structures that are shown to be associated with the generation of Reynolds shear stress -u'v'(bar). These motions are characterized by high spatial gradients of the discriminant and are believed to be an important part of the mechanism that sustains turbulence in the near-wall region.

Cantwell, Brian J.↗