MAGNETOHYDRODYNAMIC WAKES IN A VISCOUS CONDUCTING FLUID
Magnetohydrodynamic wakes in a viscous conducting fluid
SEARCH · Engineering Papers
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.
Magnetohydrodynamic wakes in a viscous conducting fluid
Explore the source record for details and available documents.
The High Frequency Radiometer (HFR) is the only instrument on the Aeroassist Flight Experiment (AFE) with sufficient temporal resolution to discern the frequency of unsteady wake oscillations. Determining both the frequency and amplitude of wake unsteadiness during AFE atmospheric entry is essential for reliably predicting the geometry and motion of the wake of future Aeroassisted Space Transfer Vehicles (ASTV). These parameters directly affect the location and size of the pay- load and the weight of the required afterbody heat protection. The purpose of the AFE is to validate the technologies required for the design of ASTVs, which will be used to exploit Earth-lunar space. This validation will be conducted at a combination of vehicle size, altitude, and velocity not obtainable in ground-based facilities. The AFE will provide the experimental flight data needed to improve our understanding of hypersonic-wake physics and to validate computational predictions of the aerodynamic and heating loads, including afterbody radiative heating loads, on an ASTV. Reliable prediction of ASTV wake flows will ensure that payloads are located within the shear-layer envelope and will determine the amount of thermal protection the payloads require. Specifically, understanding the temporal nature of the wake unsteadiness is important for two reasons. Most importantly, analysis of ground-based experiments suggests that wake unsteadiness results in a variation of as much as +/- 5 deg in the shear-flow turning angle. This angle must be reliably predicted to avoid shear-layer impingement on the vehicle afterbody, which would result in heating rates of about 10 W/sq cm, of the same order as on the forebody stagnation point. Secondly, the energy associated with wake unsteadiness will reduce the static enthalpy of the wake fluid and cause an error of as much as 30% in the amount of predicted wake radiative heating. Therefore, the HFR flight data, which will quantify the frequency and amplitude of the wake unsteadiness, are required for the verification of computational models of ASTV flowfields. Measurement of wake unsteadiness has been endorsed by the Peer Science Steering Group and the AFE Computational Fluid Dynamics (CFD) Working Group. The Peer Science Steering Group has stressed that a reliable measurement of the wake unsteadiness is fundamental because incorporating unsteadiness represents a substantial challenge to the CFD community, and a reliable measurement will raise confidence in the computer simulation. The AFE CFD Working Group has stated that 44 any data which would prove the existence and effects of unsteady flow would be extremely valuable.
This work addresses the significant differences in compressor rotor wake mixing loss which exist in a stage environment relative to a rotor in isolation. The wake decay for a rotor in isolation is due solely to viscous dissipation which is an irreversible process and thus leads to a loss in both total pressure and efficiency. Rotor wake decay in the stage environment is due to both viscous mixing and the inviscid strain imposed on the wake fluid particles by the stator velocity field. This straining process, referred to by Smith (1993) as recovery, is reversible and for a 2D rotor wake leads to an inviscid reduction of the velocity deficit of the wake. A model for the rotor wake decay process is developed and used to quantify the viscous dissipation effects relative to those of inviscid wake stretching. The model is verified using laser anemometer measurements acquired in the wake of a transonic rotor operated in isolation and in a stage configuration at near peak efficiency and near stall operating conditions. Additional insight is provided by a time-accurate 3D Navier-Stokes simulation of the compressor stator flow field at the corresponding stage loading levels. Results from the wake decay model exhibit good agreement with the experimental data. Data from the model, laser anemometer measurements, and numerical simulations indicate that for the rotor/stator spacing used in this work, which is typical of core compressors, rotor wake straining (stretching) is the primary decay process in the stator passage with viscous mixing playing only a minor role. The implications of these results on compressor stage design are discussed.
The behavior of wake vortices in a convective boundary layer is investigated using a validated large eddy simulation model. Our results show that the vortices are largely deformed due to strong turbulent eddy motion while a sinusoidal Crow instability develops. Vortex rising is found to be caused by the updrafts (thermals) during daytime convective conditions and increases with increasing nondimensional turbulence intensity eta. In the downdraft region of the convective boundary layer, vortex sinking is found to be accelerated proportional to increasing eta, with faster speed than that in an ideal line vortex pair in an inviscid fluid. Wake vortices are also shown to be laterally transported over a significant distance due to large turbulent eddy motion. On the other hand, the decay rate of the, vortices in the convective boundary layer that increases with increasing eta, is larger in the updraft region than in the downdraft region because of stronger turbulence in the updraft region.
This project assembled a broad ensemble of modeling and experimentation tools to study the morphological and optical properties of detonation soots in explosive fireballs. A gram-scale hemispherical high explosive was studied in a low-pressure controlled environment using in-situ experimentation with diffusely illuminated visible absorption spectroscopy, particle sizing through light scattering techniques, and post-test collections with subsequent morphological analysis. Hydrocode modeling was performed to replicate the detonation flow observations, and subsequent aerosol kinetics models provided particle size distributions and extinction coefficients from the hydrocode results. Experimentally observed soot morphologies agreed with expectation from the literature - a bimodal distribution was found, brought upon by the particles growing to a size where their inertia and fluid wakes are non-negligible. The aerosol kinetics model did not replicate the observed bimodal size distribution for lack of a coagulation kernel to represent the behavior. To recover particulate optical properties, a spectrally resolved absorption spectroscopy method termed Spectral diffuse back-illuminated extinction imaging (SBI-EI) was developed and implemented on two explosive types. Inverting the absorption spectra using a Kramers-Kronig consistent method yielded the complex index of refraction for the soots produced by the explosives. This method resulted in an unrealistic index of refraction for one of the two explosives, and this is suspected to be due to the model neglecting scattering brought upon by the large particle sizes observed. In addition to the core work, three additional studies were performed in parallel. These investigated the impact of scattering on diffuse absorption spectroscopy, studied how soots oxidate and sublimate in a well-controlled shock tube, and laid the theoretical groundwork for a new collision kernel to replicate the bimodal size distribution from the observations. Summaries of these efforts are included at the end of this report.
The primary focus of this paper is to investigate the loss sources in an advanced GE transonic compressor design with high reaction and high stage loading. This advanced compressor has been investigated both experimentally and analytically in the past. The measured compressor efficiency is significantly lower than the efficiency calculated with various existing tools based on RANS and URANS. The general understanding is that some important flow physics in this modern compressor design are not represented in the current tools. To pinpoint the source of the efficiency miss, an advanced test with detailed flow traverse was performed for the front one and a half stage at the NASA Glenn Research Center. In the present paper, a Large Eddy Simulation (LES) is employed to determine whether a higher-fidelity simulation can pick up any additional flow physics that can explain past efficiency miss with RANS and URANS. The results from the Large Eddy Simulation were compared with the NASA test results and the GE interpretation of the test data. LES calculates lower total pressure and higher total temperature on the pressure side of the stator, resulting in large loss generation on the pressure side of the stator. On the other hand, existing tools based on the RANS and URANS do not calculate this high total temperature and low total pressure on the pressure side of the stator. The calculated loss through the stator from LES seems to match the measured data and the GE data interpretation. Detailed examination of the unsteady flow field from LES indicates that the accumulation of high loss near the pressure side of the stator is due to the interaction of the rotor wake with the stator blade. The strong rotor wake interacts quite differently with the pressure side of the stator than with the suction side of the stator blade. The concave curvature on the pressure side of the stator blade increases the mixing of the rotor wake with the pressure side boundary layer significantly. On the other hand, the convex curvature on the suction side of the stator blade decreases the mixing and the suction side blade boundary layer remains thin. The jet velocity in the rotor wake in the stator frame seems to magnify the curvature effect in addition to inviscid redistribution of wake fluid toward the pressure side of the blade.
Measurements of drag and wake structure in magnetofluid dynamic flow about a sphere
Measurements of drag and wake structure in magneto-fluid dynamic flow about a sphere
The plasma potential is measured near the edge of an electrically floating obstacle placed in a steady-state, supersonic, unmagnetized, neutral plasma flow. Equipotential contours show the sheath of the upstream side of the obstacle wrapping around the edge of the obstacle and fanning out into the near wake. Both fluid theory and the data find the near-wake plasma potential to be self-similar when ionization, charge exchange, and magnetic field can be neglected. The theory also finds that fluid velocity is self-similar, the near wake is nonneutral, and plasma density is not self-similar. Strong electric fields are found near the obstacle and equipotential contours are found to conform to all boundaries.
A major cost in the study of wind farm layout optimization is the repeated evaluation of the annual energy production (AEP). The current approach to estimating AEP requires a large set of flow simulations to be performed that cover each discrete wind speed and direction combination contained within the wind rose, followed by a probability-weighted sum of the power production resulting from each simulation. Even with inexpensive engineering wake models, this numerical integration scheme can lead to high computational costs. In this paper, we derive an analytical formulation for estimating farm AEP across every wind direction, based on a Gaussian wake velocity model, which reduces the number of wind farm simulations to a single function evaluation. As a result, we find that the Gaussian-FLOWERS approach reduces the time for AEP calculations by more than two orders of magnitude with a small trade-off in accuracy when compared to a conventional approach. This massive reduction in computation cost is useful to reduce overall costs in wind farm layout optimization studies.
The OSIRIS-REx sample return capsule's hypersonic re-entry into the atmosphere is a rare opportunity to test a variety of sonic boom source models since the projectile dimensions are well characterized. While the as-flown flight path is unknown, the predicted flight path enables a rough approximation of the source Mach number and location. Six infrasound microphones deployed in the boom carpet along the predicted flight path recorded impulsive signals from the OSIRIS-REx re-entry. Using a suite of atmosphere profiles and the geometric acoustics approximation, we estimate locations with uncertainty estimates along the flight path from which the signals were emitted. Acoustic overpressure and signal duration predictions from Whitham's far field theory, Carlson's simplified sonic boom prediction method, and a drag-dominated hypersonic model are analyzed with uncertainty estimates from the location estimate. While the Carlson simplified sonic boom prediction method could be accurate, our preference is for the drag-dominated source model. Using this source model with an inviscid Burgers's equation solver for propagation, we obtained an excellent match to the recorded data. In conclusion, these results will help better inform future sample return capsule re-entry observation campaigns as well as contribute to a better understanding of high altitude infrasonic sources.
Laboratory experiments and fluid theory find a stationary rarefaction shock in the near wake of an electrically grounded obstacle placed in a steady state, supersonic plasma flow. The shock is only found when two electron temperatures, differing by at least an order of magnitude, are present. These shocks are analogous to rarefaction shocks in plasma free expansions.
In the invention, propagating broad band and tonal acoustic components of noise characteristic of interaction of a turbomachine blade wake, produced by a turbomachine blade as the blade rotates, with a turbomachine component downstream of the rotating blade, are reduced. This is accomplished by injection of fluid into the blade wake through a port in the rotor blade. The mass flow rate of the fluid injected into the blade wake is selected to reduce the momentum deficit of the wake to correspondingly increase the time-mean velocity of the wake and decrease the turbulent velocity fluctuations of the wake. With this fluid injection, reduction of both propagating broad band and tonal acoustic components of noise produced by interaction of the blade wake with a turbomachine component downstream of the rotating blade is achieved. In a further noise reduction technique, boundary layer fluid is suctioned into the turbomachine blade through a suction port on the side of the blade that is characterized as the relatively low-pressure blade side. As with the fluid injection technique, the mass flow rate of the fluid suctioned into the blade is here selected to reduce the momentum deficit of the wake to correspondingly increase the time-mean velocity of the wake and decrease the turbulent velocity fluctuations of the wake; reduction of both propagating broad band and tonal acoustic components of noise produced by interaction of the blade wake with a turbomachine component downstream of the rotating blade is achieved with this suction technique. Blowing and suction techniques are also provided in the invention for reducing noise associated with the wake produced by fluid flow around a stationary blade upstream of a rotating turbomachine.
A solution method for problems that contain both continuum and rarefied flow regions is presented. The methodology is applied to flow about the 3-D Mars Sample Return Orbiter (MSRO) that has a highly compressed forebody flow, a shear layer where the flow separates from a forebody lip, and a low density wake. Because blunt body flow fields contain such disparate regions, employing a single numerical technique to solve the entire 3-D flow field is often impractical, or the technique does not apply. Direct simulation Monte Carlo (DSMC) could be employed to solve the entire flow field; however, the technique requires inordinate computational resources for continuum and near-continuum regions, and is best suited for the wake region. Computational fluid dynamics (CFD) will solve the high-density forebody flow, but continuum assumptions do not apply in the rarefied wake region. The CFD-DSMC approach presented herein may be a suitable way to obtain a higher fidelity solution.
From the theoretical fluid dynamics point of view, the wake region of a large horizontal-axis wind turbine has been defined and described, and numerical models of wake behavior have been developed. Wind tunnel studies of single turbine wakes and turbine array wakes have been used to verify the theory and further refine the numerical models. However, the effects of scaling, rotor solidity, and topography on wake behavior are questions that remain unanswered. In the wind tunnel studies, turbines were represented by anything from scaled models to tea strainers or wire mesh disks whose solidity was equivalent to that of a typical wind turbine. The scale factor compensation for the difference in Reynolds number between the scale model and an actual turbine is complex, and not typically accounted for. Though it is wise to study the simpler case of wakes in flat topography, which can be easily duplicated in the wind tunnel, current indications are that wind turbine farm development is actually occurring in somewhat more complex terrain. Empirical wake studies using large horizontal-axis wind turbines have not been thoroughly composited, and, therefore, the results have not been applied to the well-developed theory of wake structure. The measurement programs have made use of both in situ sensor systems, such as instrumented towers, and remote sensors, such as kites and tethered, balloonborne anemometers. We present a concise overview of the work that has been performed, including our own, which is based on the philosophy that the MOD-2 turbines are probably their own best detector of both the momentum deficit and the induced turbulence effect downwind. Only the momentum deficit aspects of the wake/machine interactions have been addressed. Both turbine power output deficits and wind energy deficits as measured by the onsite meteorological towers have been analyzed from a composite data set. The analysis has also evidenced certain topographic influences on the operation of spatially diverse wind turbines.
A three-dimensional, full-potential, quasi-steady code TFAR1 is proposed for calculating the transonic flow past a lifting helicopter rotor blade and oblique wing. The TFAR1 uses a two-dimensional nonlinear wake-model that allows a jump in velocity potential to propagate with the local fluid flow in the wake. Rotor calculations were made for a single blade at an advance ratio of 0.3, a rotational tip Mach number of 0.7, and at 0-degree incidence. A 1/7-scale model of the Cobra Operational Load Survey (OLS) rotor blade is calculated, and the pressure distributions are compared to the measurements for azimuth angles 0, 30, 60, 90, 120, and 150 degrees at the 95 percent spanwise station of the OLS blade. Furthermore, an oblique wing with Korn airfoil was calculated at the high transonic free-stream Mach number of 0.9791, zero incidence, and yaw angle of 40 degrees. The TFAR1, coupled with a helicopter performance code CAMRAD (Johnson, 1981), provides a full-potential code for calculating the entire flow field for a multiple-bladed rotor in transonic lifting forward flight.
A model is developed to predict the behavior of a thin fluid film in the wake of a tool in relative motion with respect to the table. Computational procedures are developed and limitations of the model are discussed. In general the fluid-interface temperature is controlled by conduction into the table. From the numerical results four regimes are identified: convective cooling, partial evaporation of the fluid film, extended evaporation due to a limiting evaporative heat flux, and surface dryout due to total evaporation of the fluid layer. These regimes are qualitatively illustrated in terms of the parameters film thickness, viscosity, and relative velocity.