Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “rotating fluid”

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 1,063 records · Page 59

Simulation of Atmospheric-Entry Capsules in the Subsonic Regime

The accuracy of Computational Fluid Dynamics predictions of subsonic capsule aerodynamics is examined by comparison against recent NASA wind-tunnel data at high-Reynolds-number flight conditions. Several aspects of numerical and physical modeling are considered, including inviscid numerical scheme, mesh adaptation, rough-wall modeling, rotation and curvature corrections for eddy-viscosity models, and Detached-Eddy Simulations of the unsteady wake. All of these are considered in isolation against relevant data where possible. The results indicate that an improved predictive capability is developed by considering physics-based approaches and validating the results against flight-relevant experimental data.

Simulation↗

Decadal Polar Motion of the Earth Excited by the Convective Outer Core From Geodynamo Simulations

Long time geodetic observation records show that the orientation of the Earth's rotation axis with respect to the terrestrial reference frame, or polar motion, changes on a broad range of timescales. Apart from external torques from the luni-solar tides, these changes are excited by interactions among different components of the Earth system. The convective fluid outer core has long been conjectured a likely contributor to the observed polar motion on timescales upward of decades, such as the 30 year Markowitz wobble. We investigated the electromagnetic coupling scenario across the core-mantle boundary via numerical geodynamo simulation for different geodynamo parameters (Rayleigh numbers and magnetic Rossby numbers). Our simulated polar motion varies strongly with the dynamo parameters, while its excitation on decadal timescales appear to converge asymptotically within the adopted range of numerical Rossby numbers. Three strongest asymptotic modes emerge from numerical results, with periods around 30, 40, and 60 years for the prograde excitation and around 24, 30, and 60 years for the retrograde excitation. Their amplitudes are all larger than 5 times 10 (sup minus 8), or approximately 10 milliseconds of arc.The results suggest that the electromagnetic core-mantle coupling could explain a substantial portion, if not all, of the observed decadal polar motion. In particular, the predicted 60 year polar motion deserves special attention for future observations and studies.

Kuang, W.↗

Design and testing of a unique randomized gravity, continuous flow bioreactor

A rotating, null gravity simulator, or Couette bioreactor was successfully used for the culture of mammalian cells in a simulated microgravity environment. Two limited studies using Lipomyces starkeyi and Streptomyces clavuligerus were also conducted under conditions of simulated weightlessness. Although these studies with microorganisms showed promising preliminary results, oxygen limitations presented significant limitations in studying the biochemical and cultural characteristics of these cell types. Microbial cell systems such as bacteria and yeast promise significant potential as investigative models to study the effects of microgravity on membrane transport, as well as substrate induction of inactive enzyme systems. Additionally, the smaller size of the microorganisms should further reduce the gravity induced oscillatory particle motion and thereby improve the microgravity simulation on earth. Focus is on the unique conceptual design, and subsequent development of a rotating bioreactor that is compatible with the culture and investigation of microgravity effects on microbial systems. The new reactor design will allow testing of highly aerobic cell types under simulated microgravity conditions. The described reactor affords a mechanism for investigating the long term effects of reduced gravity on cellular respiration, membrane transfer, ion exchange, and substrate conversions. It offers the capability of dynamically altering nutrients, oxygenation, pH, carbon dioxide, and substrate concentration without disturbing the microgravity simulation, or Couette flow, of the reactor. All progeny of the original cell inoculum may be acclimated to the simulated microgravity in the absence of a substrate or nutrient. The reactor has the promise of allowing scientists to probe the long term effects of weightlessness on cell interactions in plants, bacteria, yeast, and fungi. The reactor is designed to have a flow field growth chamber with uniform shear stress, yet transfer high concentrations of oxygen into the culture medium. The system described allows for continuous, on line sampling for production of product without disturbing fluid and particle dynamics in the reaction chamber. It provides for the introduction of substrate, or control substances after cell adaptation to simulated microgravity has been accomplished. The reactor system provides for the nondisruptive, continuous flow replacement of nutrient and removal of product. On line monitoring and control of growth conditions such as pH and nutrient status are provided. A rotating distribution valve allows cessation of growth chamber rotation, thereby preserving the simulated microgravity conditions over longer periods of time.

Lassiter, Carroll B.↗

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗

CAS2D: FORTRAN program for nonrotating blade-to-blade, steady, potential transonic cascade flows

An exact, full-potential-equation (FPE) model for the steady, irrotational, homentropic and homoenergetic flow of a compressible, homocompositional, inviscid fluid through two dimensional planar cascades of airfoils was derived, together with its appropriate boundary conditions. A computer program, CAS2D, was developed that numerically solves an artificially time-dependent form of the actual FPE. The governing equation was discretized by using type-dependent, rotated finite differencing and the finite area technique. The flow field was discretized by providing a boundary-fitted, nonuniform computational mesh. The mesh was generated by using a sequence of conforming mapping, nonorthogonal coordinate stretching, and local, isoparametric, bilinear mapping functions. The discretized form of the FPE was solved iteratively by using successive line overrelaxation. The possible isentropic shocks were correctly captured by adding explicitly an artificial viscosity in a conservative form. In addition, a three-level consecutive, mesh refinement feature makes CAS2D a reliable and fast algorithm for the analysis of transonic, two dimensional cascade flows.

Dulikravich, D. S.↗

Convective effects in float-zone and Czochralski melts

The hydrodynamics of crystal-growth melts is a relatively new research area. Numerical modeling of these processes is necessary. The work discussed herein is in two parts: numerical simulations of the flow in a Czochralski melt, and also of that in a float zone. In addition, for the float-zone case, energy stability theory will be used to determine stability bounds for the onset of oscillatory thermo-capillary flow. Convective effects in crystal-growth melts arise from a variety of mechanisms. Temperature gradients both in the direction of gravity and normal to it give rise to convection due to buoyancy effects. Rotation of the crucible and/or crystal causes a forced convection which may augment or oppose the buoyancy-driven flow. Finally, thermo-capillary forces (due to the variation of surface tension with temperature) drive surface motions which in turn generate convection in the bulk fluid. All of these mechanisms are present in either Czochralski or float-zone growth. The objective of the Czochralski modeling is to develop an accurate numerical simulation of the flow in a Czochralski silicon melt and to investigate the effects of various parameters on the flow properties. Like some earlier investigations, the intent is to simulate the effects of buoyancy, forced and thermo-capillary convection, including unsteady effects. Unlike earlier work, the aim is to include the effects of a variable free surface and freezing interface and, possibly incorporate nonaxisymmetric effects.

Neitzel, G. P.↗

An investigation of a contoured wall injector for hypervelocity mixing augmentation

An experimental and computational investigation of a contoured wall fuel injector is presented. The injector was aimed at enabling shock-enhanced mixing for the supersonic combustion ramjet engines currently envisioned for applications on hypersonic vehicles. Three-dimensional flow field surveys, and temporally-resolved planar Rayleigh scattering measurements are presented for Mach 1.7 helium injection into Mach 6 air. These experimental data are compared directly with a three-dimensional Navier-Stokes simulation of the flow about the injector array. Two dominant axial vorticity sources are identified and characterized. The axial vorticity produced strong convective mixing of the injectant with the freestream. Shock-impingement was particularly effective as it assured seeding of baroclinic vorticity directly on the helium/air interface. The vorticity coalesced into a counter-rotating vortex pair of a sense which produced migration of the helium away from the wall. The influences of spatial averaging on the representation of the flow field as well as the importance of the fluctuating component of the flow in producing molecularly-mixed fluid are addressed.

Waitz, Ian A.↗

Meteorological Data for Geodynamics Use: IERS Special Bureau for the Atmosphere

During the life of this contract, the PI set up and operated a data center known as the "Special Bureau for the Atmosphere" of the International Earth Rotation Service (IERS) at Atmospheric and Environmental Research, Inc., and in conjunction with the U.S. National Oceanic and Atmospheric Administration. The role of the center is to calculate, archive, analyze, and distribute atmospheric data related to Earth rotation, polar motion and other motions of the Earth. We have produced data from 4 different operational meteorological centers. We have also produced data from atmospheric reanalyzes, spanning back about a half a century. The center is one of several centers under the Global Geophysical Fluids Center (GGFC) of the IERS; the GGFC is centered at NASA's Goddard Space Flight Center. Our data center has made data available to a large variety of scientists worldwide.

Salstein, David A.↗

Centrifugal Compressor Aeroelastic Analysis Code

Centrifugal compressors are very widely used in the turbomachine industry where low mass flow rates are required. Gas turbine engines for tanks, rotorcraft and small jets rely extensively on centrifugal compressors for rugged and compact design. These compressors experience problems related with unsteadiness of flowfields, such as stall flutter, separation at the trailing edge over diffuser guide vanes, tip vortex unsteadiness, etc., leading to rotating stall and surge. Considerable interest exists in small gas turbine engine manufacturers to understand and eventually eliminate the problems related to centrifugal compressors. The geometric complexity of centrifugal compressor blades and the twisting of the blade passages makes the linear methods inapplicable. Advanced computational fluid dynamics (CFD) methods are needed for accurate unsteady aerodynamic and aeroelastic analysis of centrifugal compressors. Most of the current day industrial turbomachines and small aircraft engines are designed with a centrifugal compressor. With such a large customer base and NASA Glenn Research Center being, the lead center for turbomachines, it is important that adequate emphasis be placed on this area as well. Currently, this activity is not supported under any project at NASA Glenn.

Keith, Theo G., Jr.↗

Navier-Stokes analysis and experimental data comparison of compressible flow within ducts

Many aircraft employ ducts with centerline curvature or changing cross-sectional shape to join the engine with inlet and exhaust components. S-ducts convey air to the engine compressor from the intake and often decelerate the flow to achieve an acceptable Mach number at the engine compressor by increasing the cross-sectional area downstream. Circular-to-rectangular transition ducts are used on aircraft with rectangular exhaust nozzles to connect the engine and nozzle. To achieve maximum engine performance, the ducts should minimize flow total pressure loss and total pressure distortion at the duct exit. Changes in the curvature of the duct centerline or the duct cross-sectional shape give rise to streamline curvature which causes cross stream pressure gradients. Secondary flows can be caused by deflection of the transverse vorticity component of the boundary layer. This vortex tilting results in counter-rotating vortices. Additionally, the adverse streamwise pressure gradient caused by increasing cross-sectional area can lead to flow separation. Vortex pairs have been observed in the exit planes of both duct types. These vortices are due to secondary flows induced by pressure gradients resulting from streamline curvature. Regions of low total pressure are produced when the vortices convect boundary layer fluid into the main flow. The purpose of the present study is to predict the measured flow field in a diffusing S-duct and a circular-to-rectangular transition duct with a full Navier-Stokes computer program, PARC3D, and to compare the numerical predictions with new detailed experimental measurements. The work was undertaken to extend previous studies and to provide additional CFD validation data needed to help model flows with strong secondary flow and boundary layer separation. The S-duct computation extends the study of Smith et al, and Harloff et al, which concluded that the computation might be improved by using a finer grid and more advanced turbulence models. The present study compares results for both the Baldwin-Lomas and k-epsilon turbulence models and is conducted with a refined grid. For the transition duct, two inlet conditions were considered, the first with straight flow and the second with swirling flow. The first case permits examination of the effects of the geometric transition on the flow field, while the second case includes the rotational flow effect characteristic of a gas turbine engine.

Harloff, G. J.↗

An investigation into the potential for wind turbines to cause barotrauma in bats

The high rates of bat mortality caused by operating wind turbines is a concern for wind energy and wildlife stakeholders. One theory that explains the mortality is that bats are not only killed by impact trauma, but also by barotrauma that results from exposure to the pressure variations caused by rotating turbine blades. To date, no published research has calculated the pressure changes that bats may be exposed to when flying near wind turbines and then used these data to estimate the likelihood that turbines cause barotrauma in bats. To address this shortcoming, we performed computational fluid dynamics simulations of a wind turbine and analytical calculations of blade-tip vortices to estimate the characteristics of the sudden pressure changes bats may experience when flying near a utility-scale wind turbine. Because there are no data available that characterize the pressure changes that cause barotrauma in bats, we compared our results to changes in pressure levels that cause barotrauma and mortality in other mammals of similar size. This comparison shows that the magnitude of the low-pressures bats experience when flying near wind turbines is approximately 8 times smaller than the pressure that causes mortality in rats, the smallest mammal for which data are available. The magnitude of the high-pressures that bats may experience are approximately 80 times smaller than the exposure level that causes 50% mortality in mice, which have a body mass similar to several bat species that are killed by wind turbines. Further, our results show that for a bat to experience the largest possible magnitude of low- and high-pressures, they must take very specific and improbable flight paths that skim the surface of the blades. Even a small change in the flight path results in the bat being hit by the blade or experiencing a much smaller pressure change. Accordingly, if bats have a physiological response to rapid low- and high-pressure exposure that is similar to other mammals, we conclude that it is unlikely that barotrauma is responsible for a significant number of turbine-related bat fatalities, and that impact trauma is the likely cause of the majority of wind-turbine-related bat fatalities.

17 WIND ENERGY↗

Alfven waves in a two-fluid model of the solar wind.

A two-fluid model for the solar wind is presented that includes the presence of Alfven waves which originate at the sun. The effective pressure of these Alfven waves is considered, as well as a model representation for proton heating through nonlinear Alfven wave damping. The effects of rotation in the solar equatorial plane are allowed for. The implications of the computational results obtained are discussed.

Hollweg, J. V.↗

Determination of physical and chemical states of lubricants in concentrated contacts, part 3

Solid and liquid thin films were analyzed by infrared emission Fourier microspectrophotometry. The apparatus used is a commercial absorption instrument modified to an emission instrument, comprising a rotating polarizing device, a miniature blackbody temperature reference adjustable in temperature and radiant flux and a microscope lens with a high numerical aperture in the entrance system for increased sensitivity and resolution. Studies of lubricant behavior in a simulated ball bearing showed the alignment of the fluid molecules in the Hertzian area. Polyphenyl ether plus 1% 1,1,2-trichloroethane (TCE) required lower shear rates for the same degree of alignment than without TCE. The experiment was run with 440 C stainless steel balls coated with TiN, a chemically inert material. In both cases, the alignment was strongly influenced by the presence of TCE. The results showed (1) the dependence of alignment of fluid molecules on flow and not on adsorption at metallic surfaces, (2) phase separation between lubricant and additive under high pressure which results in two phase flow and (3) reduction in traction of torque transmitting (traction) fluids.

Lauer, J. L.↗

CARS system for turbulent flame measurements

Simultaneous nitrogen number density and rotational-vibrational temperatures were measured in a turbulent diffusion flame with a Coherent Anti-Stokes Raman Scattering (CARS) instrument. The fuel jet was diluted with nitrogen (20 percent by volume) to allow temperature measurements across the entire jet mixing region. These measurements were compared with fluid dynamics computations. The CARS system incorporated a neodymium YAG laser, an intensified silicon photodiode array detector, and unique dynamic range enhancement methods. Theoretical calculations were based on a parabolic Navier-Stokes computer code. The comparison of these techniques will aid their development in the study of complex flowfields.

Antcliff, R. R.↗

Convection experiments in a centrifuge and the generation of plumes in a very viscous fluid

In this paper, experiments are described for which inertial effects are negligible. A small aspect-ratio tank filled with a very viscous fluid (Pr = 10 to the 6th) is used to observe the behavior of convection for Rayleigh numbers up to 6.3 x 10 to the 5th. These high values are reached by conducting the experiment in a centrifuge which provides a 130-fold increase in apparent gravity. Rotational effects are small, but cannot be totally dismissed. In this geometry, thermal boundary layer instabilities are indeed observed, and are found to be very similar to their lower Prandtl number counterparts. It is tentatively concluded that once given a certain degree of 'vulnerability' convection can develop 'plume' like instabilities, even when the Prandtl number is infinite. The concept is applied to the earth's mantle and it is speculated that 'plumes' could well be the dominant mode of small-scale convection under the lithospheric plates.

Nataf, H.-C.↗

A magnetic damper for first mode vibration reduction in multimass flexible rotors

Many rotating machines such as compressors, turbines and pumps have long thin shafts with resulting vibration problems, and would benefit from additional damping near the center of the shaft. Magnetic dampers have the potential to be employed in these machines because they can operate in the working fluid environment unlike conventional bearings. An experimental test rig is described which was set up with a long thin shaft and several masses to represent a flexible shaft machine. An active magnetic damper was placed in three locations: near the midspan, near one end disk, and close to the bearing. With typical control parameter settings, the midspan location reduced the first mode vibration 82 percent, the disk location reduced it 75 percent and the bearing location attained a 74 percent reduction. Magnetic damper stiffness and damping values used to obtain these reductions were only a few percent of the bearing stiffness and damping values. A theoretical model of both the rotor and the damper was developed and compared to the measured results. The agreement was good.

Kasarda, M. E. F.↗

High-energy particle acceleration by explosive electromagnetic interaction in an accretion disk

By examining electromagnetic field evolution occurring in an accretion disk around a compact object, we arrive at an explosive mechanism of particle acceleration. Flux-freezing in the differentially rotating disk causes the seed and/or generated magnetic field to wrap up tightly, becoming highly sheared and locally predominantly azimuthal in orientation. We show how asymptotically nonlinear solutions for the electromagnetic fields may arise in isolated plasma blobs as a result of the driving of the fluid equations by the accretion flow. These fields are capable of rapidly accelerating charged particles from the disk. Acceleration through the present mechanism from AGN can give rise to energies beyond 10 exp 20 eV. Such a mechanism may present an explanation for the extragalactic origin of the most energetic observed cosmic rays.

Haswell, C. A.↗

Flow properties along field lines in a 3-D tilted-dipole geometry

A 3-D MHD simulation of a global, tilted-dipole solar wind flow pattern is analyzed to determine flow properties along individual magnetic field lines. In the model, flow conditions near the Sun are chosen to provide a reasonable match to the interplanetary configuration prevailing during the recent south polar passage by Ulysses, i.e., a streamer belt inclined approximately 30 deg to the solar equator and speeds ranging from 325-800 km/s. Field lines all across the stream pattern are traced from 1 to 10 AU by following the motion of marker particles embedded in the flow. It is found that those field lines threading the core of the interaction region are subject to significant latitudinal and relative longitudinal displacement over this range of heliocentric distance. Thus, observations taken at a fixed latitude in the inner solar system sample, over the course of a solar rotation, field lines which connect to a range of latitudes in the outer heliosphere. Maps of the field line displacements are presented to help visualize these connections. In addition, it is found that depending upon the location relative to the CIR structure, the radial evolution of fluid elements frozen to different field lines can deviate considerably from that of the canonical solar wind. That is, for selected subsets of field lines, large speed changes (not just at shocks) can be experienced; the density variation can be far from 1/r(exp 2), and the magnetic field intensity need not decay monotonically with distance.

Pizzo, V. J.↗