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.
Rotating stall in an axial-flow compressor is suppressed by the positioning of a fixed inlet flow divider in the annular inlet flow passage upstream of the compressor. The inlet flow divider is aligned with the flow of fluid through the duct and acts to block or interfere with any rotating wave in the inlet and thereby suppresses rotating stall in the compressor.
Vortex ring collisions have attracted intense interest in both water and air studies (Baird in Proc R Soc Lond Ser Math Phys Sci 409:59–65, 1987, Poudel et al. in Phys Fluids 33:096105, 2021, Lim and Nickels in Nature 357:225, 1992, New et al. in Exp Fluids 57:109, 2016, Suzuki et al. in Geophys Res Lett 34, 2007, Yan et al. in J Fluids Eng 140:054502, 2018, New et al. in J Fluid Mech 899, 2020, Cheng et al. in Phys Fluids 31:067107, 2019, Hernández and Reyes in 29:103604, 2017, Mishra et al. in Phys Rev Fluids, 2021, Zednikova et al. in Chem Eng Technol 42:843–850, 2019, Kwon et al. in Nature 600:64–69, 2021). These toroidal structures spin around a central axis and travel in the original direction of impulse while spinning around the core until inertial forces become predominant causing the vortex flow to spontaneously decay to turbulence (Vortex Rings, https://projects.iq.harvard.edu/smrlab/vortex-rings). Previous studies have shown the collision of subsonic vortex rings resulting in reconnected vortex rings, but the production of a shock wave from the collision has not been demonstrated visibly (Lim and Nickels in Nature 357:225, 1992, Cheng et al. in Phys Fluids 31:067107, 2019). Here we present the formation of a shock wave due to the collision of explosively formed subsonic vortex rings. As the vortex rings travel at Mach 0.66 toward the collision point, they begin to trap high pressure air between them. Upon collision, high pressure air was imploded and released radially away from the axis of the collision, generating a visible shock wave traveling through and away from the colliding vortices at Mach 1.22. Our results demonstrate a pressure gradient with high pressure release creating a shock wave. We anticipate our study to be a starting point for more explosively formed vortex collisions. For example, explosives with different velocities of detonation could be tested to produce vortex rings of varying velocities.
Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.
Results from wave theory and numerical simulation of the nonlinear MHD equations are used to study the response of a conducting fluid containing an embedded magnetic field with beta less than 1 to the sudden injection of material along the field lines. It is shown that the injection produces slow shocks with configurations which are concave toward the ejecta driver. Fast-mode waves which have not steepened into the shock precede the slow shock and alter the ambient medium. When beta equals 0.1, the fast mode becomes a transverse wave for parallel propagation, while the slow wave approaches a longitudinal, or sound, wave.
While developing a measurement acquisition system to be used to retrofit aging aircraft with vehicle health monitoring capabilities, Langley Research Center developed an innovative wireless fluid-level measurement system. The NASA technology was of interest to Tidewater Sensors LLC, of Newport News, Virginia, because of its many advantages over conventional fuel management systems, including its ability to provide an accurate measurement of volume while a boat is experiencing any rocking motion due to waves or people moving about on the boat. These advantages led the company to license this novel fluid-level measurement system from NASA for marine applications.
Injecting fluids into underground geologic structures is crucial for the development of long-term strategies for managing captured carbon and facilitating sustainable energy extraction operations. Here, we have previously reported that the injection of metal–organic frameworks (MOFs) into the subsurface can enhance seismic monitoring tools to track fluids and map complex structures, reduce risk, and verify containment in carbon storage reservoirs because of their absorption capacity of low-frequency seismic waves. Here, we demonstrate that water-based Cr/Zn/Zr MOF colloidal suspensions (nanofluids) are multimodal geophysical contrast agents that enhance near-wellbore logging tools. Based on experimental fluid-only measurements, MIL-101(Cr), ZIF-8, and UiO-66 nanofluids have distinct complex conductivity and/or low-field nuclear magnetic resonance (NMR) signatures that are relevant to field-deployed technologies, implying the potential to enhance near-wellbore monitoring of CO 2 injection and associated processes with downhole logging tools. Small- and wide-angle X-ray scattering characterization of ~0.5 wt % MIL-101(Cr) suspensions confirmed phase stability and provided insight into the fractal nature of colloidal nanoparticles. Finally, low-field (2 MHz) NMR measurements of MIL-101(Cr) nanofluid injection into a prototypical Berea sandstone demonstrate how paramagnetic high-surface area MOFs may dominate the relaxation times of hydrogen-bearing fluids in porous geologic matrices, enhancing the mapping of near-surface and near-wellbore transport pathways and advancing sustainable subsurface energy technologies.
The prediction of force on an isolated particle, while a shock is passing over it, is an important problem in many natural and industrial applications. Although the flow monotonically changes from the pre-shock to the post-shock state, the particle's force has been observed to behave nonmonotonically with a sharp peak when the shock is located halfway across the particle. This nonmonotonic behavior is due to the unsteady nature of the compression and rarefaction waves that radiate as the shock diffracts around the particle and, therefore, cannot be predicted by a quasi-steady model. An accurate force model must account for the unsteady nature of the flow and the sharp discontinues in the flow properties across the shock. In this work, we test four different inviscid models and observe that the compressible Maxey–Riley–Gatignol (C-MRG) model is the most accurate based on comparison with results from particle-resolved inviscid simulations at two different Mach numbers for both water and air as the medium. The C-MRG model is first demonstrated to predict the force on a stationary particle accurately and then extended to capture the force on a moving particle. Numerical complexities regarding the implementation of the C-MRG model are also discussed.
A procedure for determining the scattered pressure field resulting from a monochromatic harmonic wave that is incident upon a layer energy absorbing structure is treated. The situation where the structure is modeled with finite elements and the surrounding acoustic medium (water or air) is represented with either acoustic finite elements, or some type of boundary integral formulation, is considered. Finite element modeling problems arise when the construction of the structure, at the fluid structure interface, are nonhomogeneous and in particular, when the inhomogeneities are small relative to the acoustic wave length. An approximate procedure is presented for replacing the detailed microscopic representation of the layered surface configuration with an equivalent simple surface impedance finite element, which is especially designed to work only at limited frequencies. An example problem is presented using NASTRAN. However, the procedure is general enough to adapt to practically any finite element code having a steady state option.
A plasma simulation model was used to study the resonant interactions between streaming cosmic-ray ions and a self-consistent spectrum of Alfven waves, such as might exist in the interstellar medium upstream of a supernova remnant shock wave. The computational model is a hybrid one, in which the background interstellar medium is an MHD fluid and the cosmic-rays are discrete kinetic particles. The particle sources for the electromagnetic fields are obtained by averaging over the fast cyclotron motions. When the perturbed magnetic field is larger than 10 percent of the background field, the macro- and microphysics are no longer correctly predicted by quasi-linear theory. The particles are trapped by the waves and show sharp jumps in their pitch-angles relative to the background magnetic field, and the effective ninety-degree scattering time for diffusion parallel to the background magnetic field is reduced to between 5 and 30 cyclotron periods. Simulation results suggest that Type 1 supernova remnants may be the principal sites of cosmic ray acceleration.
An ultrasonic method is presented for non-intrusively measuring hydraulic fluid level in aircraft struts in the field quickly and easily without modifying the strut or aircraft. The technique interrogates the strut with ultrasonic waves generated and received by a removable ultrasonic transducer hand-held on the outside of the strut in a fashion that is in the presence or absence of hydraulic fluid inside the strut. This technique was successfully demonstrated on an A-6 aircraft strut on the carriage at the Aircraft Landing Dynamics Research Facility at NASA Langley Research Center. Conventional practice upon detection of strut problem symptoms is to remove aircraft from service for extensive maintenance to determine fluid level. No practical technique like the method presented herein for locating strut hydraulic fluid level is currently known to be used.
Laser wakefield accelerators rely on the extremely high electric fields of nonlinear plasma waves to trap and accelerate electrons to relativistic energies over short distances. When driven strongly enough, plasma waves break, trapping a large population of the background electrons that support their motion. Aside from limiting the maximum electric field, this trapping can lead to accelerated electron bunches with large energy spreads. Here we introduce a novel regime of plasma wave excitation and wakefield acceleration that allows for arbitrarily high electric fields while avoiding the deleterious effects of unwanted trapping. Here, the regime, enabled by spatiotemporal shaping of laser pulses, exploits the property that nonlinear plasma waves with superluminal phase velocities cannot trap charged particles and are therefore immune to wave breaking.
One of the primary motivations behind this work is the attempt to understand the physics of a two-phase jet pump which constitutes part of a flow boiling test facility at NASA-Marshall. The flow boiling apparatus is intended to provide data necessary to design highly efficient two-phase thermal control systems for aerospace applications. The facility will also be capable of testing alternative refrigerants and evaluate their performance using various heat exchangers with enhanced surfaces. The test facility is also intended for use in evaluating single-phase performance of systems currently using CFC refrigerants. Literature dealing with jet pumps is abundant and covers a very wide array of application areas. Example application areas include vacuum pumps which are used in the food industry, power station work, and the chemical industry; ejector systems which have applications in the aircraft industry as cabin ventilators and for purposes of jet thrust augmentation; jet pumps which are used in the oil industry for oil well pumping; and steam-jet ejector refrigeration, to just name a few. Examples of work relevant to this investigation includes those of Fairuzov and Bredikhin (1995). While past researchers have been able to model the two-phase flow jet pump using the one-dimensional assumption with no shock waves and no phase change, there is no research known to the author apart from that of Anand (1992) who was able to account for condensation shocks. Thus, one of the objectives of this work is to model the dynamics of fluid interaction between a two-phase primary fluid and a subcooled liquid secondary fluid which is being injected employing atomizing spray injectors. The model developed accounts for phase transformations due to expansion, compression, and mixing. It also accounts for shock waves developing in the different parts of the jet pump as well as temperature and pressure dependencies of the fluid properties for both the primary two-phase mixture and the secondary subcooled liquid. The research effort on which this document partly reports described a relatively simple model capable of describing the performance of a two-phase flow jet pump. The model is based on the isentropic homogeneous expansion/compression hypothesis and is capable of fully incorporating the effects of shocks in both the mixing chamber and the throat/diffuser parts of the pump. The physical system chosen is identical to that experimentally tested by Fairuzov and Bredikhin (1995) and should therefore be relatively easy to validate.
One of the outstanding issues regarding the excitation of magnetosonic waves has been the observational evidence that obliquely propagating waves are dominant at comets and planetary foreshocks despite the predictions of linear theory that maximum growth occurs at parallel propagation. To address this issue, we have conducted a detailed linear theory using a beam-ring distribution function. The results have shown that such distribution functions are associated with four separate instabilities. Two of these instabilities are similar to the right hand resonant ion/ion and the non-resonant instabilities which are also present in the case of a field aligned beam. The other two instabilities are associated with the presence of ring part of the distribution function. One of these instabilities excites magnetosonic waves with maximum growth in the oblique directions. The other excites Alfven waves with maximum growth in the oblique directions. The importance of the former instability is that it may explain the oblique nature of the magnetosonic waves observed at planetary foreshocks and comets. In order to understand the nonlinear properties of these various instabilities, we have also conducted 2-D hybrid (particle ions, fluid electrons) simulations. We have found that when the beam density is sufficiently large so that the non-resonant instability has the largest growth rate, these waves dominate the initial wave power in the system. At later times, however, the obliquely propagating magnetosonic waves become dominant. Another important finding was that when the two instabilities which excite magnetosonic waves have the largest growth rates, the system is dominated by the obliquely propagating magnetosonic waves. This is despite the fact that the largest growth rate for one of the instabilities occurs in the parallel direction. The exact cause of this is not currently understood and is under investigation.
Pulsing reaction control system (RCS) thrusters, vent valves, or other propulsion devices can preserve propellant resources in space, but this operation also effectively introduces a vibration to the vehicle. When the vibration is perpendicular to the liquid propellant surface, Faraday waves may be generated at the liquid-vapor interface. These Faraday instabilities can perturb or break up the liquid surface of cryogenic tanks, leading to inefficiencies in thermal management or even ullage collapse. Drawing from theory and experiments, an engineering model defining the allowable design regions for pulsed settling in microgravity was assembled and verified with computational fluid dynamics (CFD) simulations. A traditional settling metric, the Bond number, was also overlaid in the model to indicate which duty cycles were insufficient to overcome surface tension and aggregate propellant. Mission planners and engineers can consult the tool to rapidly evaluate the stability of a liquid-vapor interface given the pulse frequency and the excitation acceleration. Expressions developed for Faraday waves induced by a sinusoidal forcing input at standard gravity were found to provide excellent predictive capabilities for pulsed, or rectangular, waveforms in the absence of a consistent gravitational acceleration. This study extends the usage of these equations to an alternative forcing function and microgravity environments for the purpose of estimating natural frequencies, surface mode shapes, surface wave amplitudes, and the onset of droplet ejection. CFD simulations with the Loci/STREAM-VoF (Volume of Fluid) solver were initially validated against experimental results in standard gravity. Discrete points on the design map were then investigated with CFD and confirmed that the engineering model reliably indicates surface stability and most Faraday wave characteristics without requiring higher-fidelity tools. The engineering model is highly extensible and can be adapted for various propellant fill fractions, fluid properties, and tank sizes.
An analytic model for predicting the effect of unsteady local surface injection on the flow separating from a streamlined body at angle of attack is proposed. The model uses the premise that separation control results from enhanced mixing along the shear layer that develops between the main stream and the fluid in the underlying recirculation zone. High-Reynolds-number asymptotic methods are used to connect the unsteady surface injection to an instability wave propagating on the separating shear layer and then to the large-scale coherent structures that produce the increased mixing. The results is a tool that can guide the choice of fluid-actuator parameters to maximize flow-control effectiveness and may also facilitate computer-based numerical experiments.
Mixing of two fluids generated by steady and particularly g-jitter acceleration is fundamental towards the understanding of transport phenomena in a microgravity environment. We propose to carry out flight and ground-based experiments to quantify flow fields due to g-jitter type of accelerations using Stereo Imaging Velocimetry (SIV), and measure the concentration field using laser fluorescence. The understanding of the effects of g-jitter on transport phenomena is of great practical interest to the microgravity community and impacts the design of experiments for the Space Shuttle as well as the International Space Station. The aim of our proposed research is to provide quantitative data to the community on the effects of g-jitter on flow fields due to mixing induced by buoyancy forces. The fundamental phenomenon of mixing occurs in a broad range of materials processing encompassing the growth of opto-electronic materials and semiconductors, (by directional freezing and physical vapor transport), to solution and protein crystal growth. In materials processing of these systems, crystal homogeneity, which is affected by the solutal field distribution, is one of the major issues. The understanding of fluid mixing driven by buoyancy forces, besides its importance as a topic in fundamental science, can contribute towards the understanding of how solutal fields behave under various body forces. The body forces of interest are steady acceleration and g-jitter acceleration as in a Space Shuttle environment or the International Space Station. Since control of the body force is important, the flight experiment will be carried out on a tunable microgravity vibration isolation mount, which will permit us to precisely input the desired forcing function to simulate a range of body forces. To that end, we propose to design a flight experiment that can only be carried out under microgravity conditions to fully exploit the effects of various body forces on fluid mixing. Recent flight experiments, by the P.I. through collaboration with the Canadian Space Agency (STS-85, August 1997), aimed at determining the stability of the interface between two miscible liquids inside an enclosure show that a long liquid column (5 cm) under microgravity isolation conditions can be stable, i.e. the interface remains sharp and vertical over a short time scale; thus transport occurs by molecular mass diffusion. On the other hand, when the two liquids were excited from a controlled vibration source (Microgravity Vibration Isolation Mount) two to four mode large amplitude quasi-stationary waves were observed. The data was limited to CCD recording of the dynamics of the interface between the two fluids. We propose to carry out flight experiments to quantify the dynamics of the flow field using Stereo Imaging Velocimetry and measure the concentration field using laser fluorescence. The results will serve as a basis to understand effects of g-jitter on transport phenomena, in this case mass diffusion. As the measurement of the kinematics of the flow field will shed light on the instability mechanism. The research will allow measurement of the flow field in microgravity environment to prove two hypotheses: (1) Maxwell's hypothesis: finite convection always exists in diffusing systems, and (2) Quasi-stationary waves inside a bounded enclosure in a microgravity environment is generated by Kelvin-Helmholtz instability; resonance of the interface which produces incipient mixing is due to Rayleigh-Taylor instability. The first hypothesis can be used as a benchmark experiment to illustrate diffusive mixing. The second hypothesis will lead to the understanding of g-jitter effects on buoyancy driven flow fields which occur in many situations involving materials processing, and other basic fluid physics phenomena. In addition, the second hypothesis will also provide insight in how Rayleigh-Taylor and Kelvin-Helmholtz instabilities propagate concentration fronts during mixing. Measurement of the flow field using SIV is important because it is the flow field which causes instability at the interface between the two fluids. Mixing driven by buoyancy induced flow fields will be addressed both experimentally and computationally. The experimental effort will address the kinematics of mixing: stretching, transport and chaos. Quantification of the mechanisms of mixing will consists of measuring the flow field using the SIV system at Glenn and capturing the dynamics of the interface, to measure mass transport, using a CCD camera. These experiments will be carried out within the framework of Earth's gravity and g-jitter microgravity acceleration as in a Space Shuttle environment or the International Space Station. The g-jitter will be induced and controlled using a tunable vibration isolation platform to isolate against vibration as well as input periodic and random vibration to the system. The parametric range of the microgravity experiment will be extended from the experiments on STS-85 to investigate higher mode quasi-stationary waves (8 to 12), as well as resonance regions which leads to chaos and turbulence. Ground-based experiments will focus on effects of vibration on stably stratified fluid layers in order to scale for possible scenarios in a microgravity environment. These vibrations will be subjected perpendicular to the concentration field on the ground since the parallel case can only be carried out in a microgravity environment. The concept of dynamical similarity will be applied to tune the experiments as closely as possible to a Space Shuttle environment or the International Space Station. The computational effort will take advantage of the Computational Laboratory at Glenn to corroborate the experimental findings with predictions of the dynamics of the flow field using the codes FLUENT (finite difference based) and FIDAP (finite element based). We will investigate two important cases, single-fluid model to address dilute systems with negligible jump in viscosity and the more general two-fluid model which accounts for finite jump in viscosity. Apart from its microgravity relevance, this experiment is well suited to study dynamics in nonlinear systems.
We have developed a suite of 3D electromagnetic field solvers, both FEM and FDTD based, that account for the RF antenna and vacuum vessel geometries with unprecedented accuracy. Workflows were developed that make it possible to translate CAD models for the antenna and vacuum vessel to physics meshes for RF wave simulation. Nonlinear RF sheath formation has been incorporated self-consistently as a boundary condition in these solvers. We have also carried out extensive studies of the impact of RF sheaths on the ion energy angle distribution at plasma-material interfaces, using high fidelity particle-in-cell codes. Comprehensive simulation models were developed to assess the impact of blob-like edge turbulence on RF wave propagation and the impact of the RF ponderomotive force on the plasma scrape-off layer (SOL). A fluid transport solver for the far-SOL was also developed which accounts for the high parallel to perpendicular heat anisotropy on an unstructured mesh, thus making it possible to precisely represent an antenna structure in the presence of edge transport. Finally we have developed a hierarchy of core wave propagation and absorption models that self-consistently combine continuum Fokker Planck and Monte Carlo treatments of fast ion evolution with ICRF full-wave field solvers and continuum Fokker Planck treatments of fast electron evolution with both full-wave and ray tracing models.
Electromagnetic wave generation and resulting cross-field diffusion of plasma are considered at a tangential discontinuity, which characterizes the magnetopause for northward interplanetary magnetic field. Two-dimensional hybrid (particle ions, massless fluid electrons) simulations, in which the tangential discontinuity is generated self-consistently via a stream-stream interaction, are used to show that wave growth occurs when the ambient magnetic field is predominantly perpendicular to the direction of the density gradient. Low-frequency (much less than ion gyrofrequency) waves, with amplitudes delta B/B less than or equal to 0.2 and anticorrelated density fluctuations delta n/n less than or equal to 0.6, are generated at the discontinuity, resulting in cross-field diffusion that is comparable to the Bohm rate. Both the fluctuation level and the lack of ion heating in the calculations are consistent with observations at the magnetopause. The magnitude of the diffusion is considered in the presence of numerical effects and in the context of the inferred diffusion rate at the magnetopause. The relation of the low-frequency waves and their consequences to faster growing, short-wavelength waves due to the lower hybrid drift instability is also addressed. The overall conclusion of this initial study is that diffusion due to low frequency waves is not likely to be a major effect at the magnetopause.