SEARCH · Engineering Papers

Results for “Fluid instabilities”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Nonclassical aileron buzz in transonic flow

A computational study of inviscid, transonic aileron and trailing-edge buzz instabilities is presented. A mixed Eulerian-Lagrangian formulation is used to model the fluid-structure system and to obtain a system of space-discretized equations that is time-marched to simulate the aeroelastic behavior of the wing-aileron system. Results obtained suggest that shock-induced separation may not be an essential driving force behind all buzz phenomena. Several examples are shown where the shock motion interacts with the aileron motion to extract energy from the flow. If the trailing-edge region is sufficiently flexible and the shocks are at the trailing edge, a trailing-edge buzz instability appears possible.

Bendiksen, Oddvar O.↗

DOE OSTI2023

Fluid dynamic mathematical aspects of supernova remnants

Supernovae—explosions of stars—are a central problem in astrophysics since they contain information on the entire process of stellar evolution and nucleosynthesis. Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities, developing during the supernova blast, lead to intense interfacial RT/RM mixing of the star's materials and couple astrophysical to atomic scales. This work analyzes some fluid dynamic mathematical aspects of the titanic task of supernova's blast. We handle mathematical challenges of RT/RM dynamics in supernova relevant conditions by directly linking the conservation laws governing RT/RM dynamics to symmetry-based momentum model, by exactly deriving the model parameters in the scale-dependent and scale-invariant regimes, and by exploring the special self-similar class for RT/RM interfacial mixing with variable accelerations. Here we reveal that RT/RM dynamics is strongly influenced by deterministic (the initial and the flow) conditions in the scale-dependent linear and nonlinear regimes and in the self-similar mixing regime. The theory outcomes are consistent with the observations of supernova remnants, explain the results of the scaled laboratory experiments in high energy density plasmas, and yield the design of future experiments for the accurate quantification of RT/RM dynamics in supernova relevant conditions. We find that from fluid dynamic mathematical perspectives, supernovae can be regarded as an astrophysical initial value problem. Along with the guidance of what explodes at microscopic scales, supernova remnants encapsulate information on the explosion hydrodynamics and the associated deterministic conditions at macroscopic scales. We urge such effects be considered in interpretations of the observational data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Simultaneous measurements of kinetic and scalar energy spectrum time evolution in the Richtmyer–Meshkov instability upon reshock

The Richtmyer–Meshkov instability (Richtmyer, Commun. Pure Appl. Maths , vol. 13, issue 2, 1960, pp. 297–319; Meshkov, Fluid Dyn. , vol. 4, issue 5, 1972, pp. 101–104) of a twice-shocked gas interface is studied using both high spatial resolution single-shot (SS) and lower spatial resolution, time-resolved, high-speed (HS) simultaneous planar laser-induced fluorescence and particle image velocimetry in the Wisconsin Shock Tube Laboratory's vertical shock tube. The initial condition (IC) is a shear layer with broadband diffuse perturbations at the interface between a helium–acetone mixture and argon. This IC is accelerated by a shock of nominal strength Mach number $M = 1.75$ , and then accelerated again by the transmitted shock that reflects off the end wall of the tube. An ensemble of experiments is analysed after reshock while the interface mixing width grows linearly with time. The kinetic and scalar energy spectra and the terms of their evolution equation are calculated and compared between SS and HS experiments. The inertial range scaling of the scalar power spectrum is found to follow Gibson's relation (Gibson, Phys. Fluids , vol. 11, issue 11, 1968, pp. 2316–2327) as a function of Schmidt number when the effective turbulent Schmidt number is used in place of the material Schmidt number that controls equilibrium scaling. Further, the spatially integrated scalar flux follows similar behaviour observed for the kinetic energy in large eddy simulation studies by Zeng et al. ( Phys. Fluids , vol. 30, issue 6, 2018, 064106) while the spatially varying scalar flux exhibits back scatter along the centre of the mixing layer and forward energy transfer in the spike and bubble regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1972

Hydrodynamic stability of viscoelastic fluids - Importance of fluid model, overstability, and form of disturbance.

A linearized stability analysis has been applied to a fluid flowing in a gravity field between horizontal planes in Couette flow under conditions such that the temperature of the bottom plane exceeds that of the top. It is shown that, under conditions likely to be encountered with polymer solutions, oscillatory instabilities will not be controlling. Criteria are offered for ascertaining when an analysis based upon a second-order fluid model may be expected to yield physically meaningful results. It is also shown that, for the fluid model considered, critical conditions for stability are not changed when disturbances which vary in the flow direction are substituted for those which are a function of the coordinate transverse to the flow.

Mcintire, L. V.↗

DOE OSTI2023

Intermediate Temperature Proton‐Conducting Solid Oxide Electrolysis Cells with Improved Performance and Durability

Technical understanding: A high-performance and stable intermediate temperature proton-conducting solid oxide electrolysis cell (PCEC) will add significantly to the hydrogen roadmap of EERE. However, the low electrocatalytic activity of the air electrode and the instability of the electrolyte are technical barriers. In this research, we utilized computational fluid dynamic modeling and materials/microstructure engineering to investigate the complex reaction chain and to improve the performance and reliability of the protonic ceramic electrochemical cells. An active and stable Pr 1.8 Ba 0.2 NiO 4+δ -BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3-δ (BZCYYb) scaffold conformal coating electrode was rationally designed where the robust Pr 1.8 Ba 0.2 NiO 4+δ functions both as a catalytic active and a protective layer. The PCECs with porosity-optimized scaffold conformal coating electrodes demonstrate 1.46 Acm -2 under 1.3 V at 600°C, which is among the highest in recent years. The final optimization by computational fluid dynamic modeling and conformal coating technique was applied to evolve the PCEC air electrode from Pr 2 NiO 4+δ to Pr 1.8 Ba 0.2 NiO 4+δ -BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3-δ scaffold conformal coating. The electrolysis current was boosted from 0.445A/cm 2 to 1.46A/cm 2 at 600°C and 1.3V. A superior low degradation rate of 1mV/kh was achieved in a 5000-hour of operation. Technical effectiveness and economic feasibility: The effectiveness of this TCO conformal coating on electrolyte backbone as air electrode has been proven by long-term demonstration. Based on the highly conductive BZCYYb electrolyte and the TCO steam electrode PNO, this project rationally took advantage of the existent materials to structure a PCEC device with high TRL and low risk. The solution infiltration method used to structure the conformal coating is a standard method used by the solid oxide cell community. Very minimal adoption is needed in the current SOC fabrication process. It is a cost-effective, scalable method for manufacturing, and brings up remarkable impact on the current technical status of PCEC. Benefit to the public: The findings in the research showed the industry that the severe challenge in PCEC can be overcome in a practical way. A various energy-hydrogen ecosystem is not beyond reach. The success of this technology provides sustainable and affordable H 2 to society. It will stimulate investment from the governmental and private sectors in response to the device manufacture and deployment, and the related infrastructure construction. It will generate high-pay domestic jobs and environmental benefits to the general public in the United States.

08 HYDROGEN↗

NASA NTRS1988

The strong nonlinear interaction of Tollmien-Schlichting waves and Taylor-Goertler vortices in curved channel flow

Viscous fluid flows with curved streamlines can support both centrifugal and viscous traveling wave instabilities. Here the interaction of these instabilities in the context of the fully developed flow in a curved channel is discussed. The viscous (Tollmein-Schlichting) instability is described asymptotically at high Reynolds numbers and it is found that it can induce a Taylor-Goertler flow even at extremely small amplitudes. In this interaction, the Tollmein-Schlichting wave can drive a vortex state with wavelength either comparable with the channel width or the wavelength of lower branch viscous modes. The nonlinear equations which describe these interactions are solved for nonlinear equilibrium states.

Bennett, J.↗

NASA NTRS1987

Interaction of Tollmien-Schlichting waves and Goertler vortices

There are many fluid flows of practical interest where transition can be caused by competing hydrodynamic instabilities. Thus in three-dimensional boundary-layer flows over curved walls, instability might be caused by Tollmien-Schlichting waves, Goertler vortices or crossflow vortices. If a particular type of instability is suppressed by some means, there is the possibility that another one might be stimulated. Hence it is important to understand the mechanisms by which these different instabilities interact. Some properties of the interaction which can take place between Goertler vortices and Tollmien-Schlichting waves are discussed.

Hall, Phil↗

NASA NTRS1992

Propagating confined states in phase dynamics

Theoretical treatment is given to the possibility of the existence of propagating confined states in the nonlinear phase equation by generalizing stationary confined states. The nonlinear phase equation is set forth for the case of propagating patterns with long wavelengths and low-frequency modulation. A large range of parameter values is shown to exist for propagating confined states which have spatially localized regions which travel on a background with unique wavelengths. The theoretical phenomena are shown to correspond to such physical systems as spirals in Taylor instabilities, traveling waves in convective systems, and slot-convection phenomena for binary fluid mixtures.

Brand, Helmut R.↗

DOE OSTI2026

Capturing Secondary Kinetic Instabilities in Three‐Dimensional Dayside Reconnection Using an Improved Gradient‐Based Closure

Magnetic reconnection is a highly dynamic process that excites a wide variety of kinetic waves and instabilities. Transverse current sheet instabilities such as the lower-hybrid drift and secondary drift-kink instabilities in particular have been shown by kinetic simulations to modify the reconnection and introduce significant turbulence and mixing to the reconnection layer. Past studies using the ten-moment fluid model to capture important kinetic physics such as the electron inertia and full representation of the pressure tensor proved advantageous to a two-fluid representation of reconnection, but the model struggled when using a local relaxation closure for the heat flux to replicate the current sheet instabilities and subsequent mixing seen in kinetic simulations. This work uses the Gkeyll software framework to perform simulations of asymmetric reconnection based on the 16 October 2015 MMS crossing of a diffusion region, the Burch event. An improved gradient-based heat flux closure is implemented, showing significant improvement in secondary kinetic instabilities that grow in the current sheet. These instabilities generate turbulence which leads to growth of secondary magnetic islands and flux ropes.

Bradshaw, K. [Princeton University, NJ (United Sta↗

NASA NTRSDate not supplied

Surface Instability of Liquid Propellants in Microgravity During Pulsed Settling Operations

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.

Faraday Waves↗

DOE OSTI2022

Full fluid moment modeling of rotating spokes in Penning-type configuration

Rotating spokes are observed in a partially magnetized plasma using a two-dimensional full fluid moment (FFM) model. In the present setup, where the radial electric field and plasma density gradient exist in opposite directions, it is observed that the spokes propagate in the direction of the diamagnetic drift and not the E × B drift. This is contrary to the modified Simon–Hoh instability, and the results suggest that the spokes can be driven by a strong diamagnetic drift. Different parameters, including magnetic field amplitude and physical domain size, influence the growth of the rotational instability as well as the dominant wave modes that arise. Furthermore, the propagation speed of the rotating spokes obtained from the FFM simulation are in good agreement with the observations in experimental and other computational work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2020

Cross-field Streaming Instabilities in the Lower Hybrid Frequency Near the Front Pole Regions of a Magnetically Shielded Hall Thruster

The mechanism that led to the erosion of the inner front pole cover in wear tests of the Hall-effect Rocket with Magnetic Shielding (HERMeS) at 300 V and 12.5 kW has been elusive. Several classical mechanisms have been considered yet none has explained conclusively the source of the discrepancy between wear measurements and simulations. In this article we consider a new mechanism. The electrostatic dispersion relation for an unbounded, homogeneous plasma in a fixed magnetic field is solved numerically in the near-plume region of this thruster. The plasma conditions have been obtained from Hall2De (r-z) multifluid-PIC simulations. We find that the modified two-stream instability can be excited in most regions of the domain, with the fastest growth approaching the fluid-limit value of ωLH/2 near the magnetic poles where ωpe/ωce~1. In this region the relative drift between magnetized electrons and beam ions perpendicular to the magnetic field can exceed 3uTi, while that between counter-streaming beam and cathode ions can exceed 4uTi. The latter suggests that the ion-ion cross-field lower hybrid instability may also be active here. Unlike the classical ion-acoustic instability, lower hybrid instabilities can heat ions to temperatures that are comparable to the electron temperature. Such anomalous heating would enhance the average energy of these ions and, in turn, the erosion of the nearby surfaces they sputter. An idealized model that scales the anomalous ion heating with ωLH has been incorporated in Hall2De and the simulations now yield significantly improved agreement with the measured wear rates and IVDFs, as reported in our companion paper.

Ortega, Alejandro Lopez↗

DOE OSTI2023

The Kelvin–Helmholtz Instability at the Boundary of Relativistic Magnetized Jets

We study the linear stability of a planar interface separating two fluids in relative motion, focusing on conditions appropriate for the boundaries of relativistic jets. The jet is magnetically dominated, whereas the ambient wind is gas-pressure-dominated. We derive the most general form of the dispersion relation and provide an analytical approximation of its solution for an ambient sound speed much smaller than the jet Alfvén speed v A , as appropriate for realistic systems. The stability properties are chiefly determined by the angle ψ between the wavevector and the jet magnetic field. For ψ = π/2, magnetic tension plays no role, and our solution resembles the one of a gas-pressure-dominated jet. Here, only sub-Alfvénic jets are unstable ($0\lt {M}_{e}\equiv (v/{v}_{{\rm{A}}})\cos \theta \lt 1$ , where v is the shear velocity and θ the angle between the velocity and the wavevector). For ψ = 0, the free energy in the velocity shear needs to overcome the magnetic tension, and only super-Alfvénic jets are unstable ($1\lt {M}_{e}\lt \sqrt{{(1+{{\rm{\Gamma }}}_{w}^{2})/[1+({v}_{{\rm{A}}}/c)}^{2}{{\rm{\Gamma }}}_{w}^{2}]}$, with Γ w the wind adiabatic index). Our results have important implications for the propagation and emission of relativistic magnetized jets.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2004

Multi-Element Unstructured Analyses of Complex Valve Systems

The safe and reliable operation of high pressure test stands for rocket engine and component testing places an increased emphasis on the performance of control valves and flow metering devices. In this paper, we will present a series of high fidelity computational analyses of systems ranging from cryogenic control valves and pressure regulator systems to cavitating venturis that are used to support rocket engine and component testing at NASA Stennis Space Center. A generalized multi-element framework with sub-models for grid adaption, grid movement and multi-phase flow dynamics has been used to carry out the simulations. Such a framework provides the flexibility of resolving the structural and functional complexities that are typically associated with valve-based high pressure feed systems and have been difficult to deal with traditional CFD methods. Our simulations revealed a rich variety of flow phenomena such as secondary flow patterns, hydrodynamic instabilities, fluctuating vapor pockets etc. In the paper, we will discuss performance losses related to cryogenic control valves, and provide insight into the physics of the dominant multi-phase fluid transport phenomena that are responsible for the choking like behavior in cryogenic control elements. Additionally, we will provide detailed analyses of the modal instability that is observed in the operation of the dome pressure regulator valve. Such instabilities are usually not localized and manifest themselves as a system wide phenomena leading to an undesirable chatter at high flow conditions.

Sulyma, Peter↗

NASA NTRSDate not supplied

Fluid System Stability Analysis Techniques

Fluid systems, or networks, consist of multiple components that work together to achieve some desired thermofluid state. For any generic application, this fluid state can be a combination of the fluid pressure, flow rate, enthalpy, or species concentration. Fluid system components, such as pumps and valves, are often governed by nonlinear differential equations, resulting in complex component-to-component interactions. System-level fluid network stability occurs when the flow through the system can maintain a steady-state solution in the presence of small perturbations, which depends on these component interactions. System instability, however, can go undetected until issues arise during integrated system testing. This presentation explores a method for system designers to think of the fluid network as an assembly of components, each with their own thermofluid surfaces of partial stability, called nullclines. The intersections of all nullclines yields system-level solutions, called equilibrium points. When designers define operating points, they are tuning system parameters so that these equilibrium points move to the desired location in the thermofluid state plane. However, linearization theory shows us that the dynamic behavior around these points can be unstable. The local stability of these equilibrium points can be assessed analytically with eigen-analysis, or numerically by propagating state-plane samples to construct a phase portrait. Investigating a phase portrait can help designers gain a qualitative understanding of a system’s dynamic performance. This understanding can then help inform requirement definitions, component selection, and operational procedures. This presentation includes an example of the phase portrait technique on a system featuring a centrifugal pump and a back-pressure regulator (BPR). Numerical modeling of this system suggests that equilibrium points on the left-hand side of the pump curve are dynamically unstable.

Fluid System↗

NASA NTRSDate not supplied

Fluid System Stability Analysis Techniques

Fluid System↗

DOE OSTI2020

Flow Z-pinch plasma production on the FuZE experiment

The Fusion Z Pinch Experiment (FuZE) investigates sheared-flow stabilization of classic m = 0 and m = 1 instabilities in Z pinches with an embedded axial flow. FuZE consists of a 100 cm coaxial plasma accelerator, where neutral gas is ionized and accelerated in a pulsed electrical discharge, followed by a 50 cm assembly region, where pinches are formed. Maintaining the pinch requires continuous plasma injection provided by a deflagration mode in the coaxial accelerator. Two discharge modes, with and without deflagration, are investigated on FuZE. Pinch formation is observed with deflagration only. Plasma velocities in the assembly region are found to match the E × B velocity estimated in the accelerator based on a 1D circuit model, indicating that a 1D MHD approximation may offer a valid description of the plasma in the accelerator channel. The velocity of magnetic field propagation is found to agree with the snowplow model based on momentum conservation, and the lifetime of the pinch is shown to be in agreement with constraints imposed by mass conservation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2004

Injection Induced Mixing in Flows Separating From Smooth Surfaces

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.

Adamczyk, John J.↗