Engineering Papers⌕ Search

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.

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

At least 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

On the nonlinear interfacial instability of rotating core-annular flow

The interfacial stability of rotating core-annular flows is investigated. The linear and nonlinear effects are considered for the case when the annular region is very thin. Both asymptotic and numerical methods are used to solve the flow in the core and film regions which are coupled by a difference in viscosity and density. The long-term behavior of the fluid-fluid interface is determined by deriving its nonlinear evolution in the form of a modified Kuramoto-Sivashinsky equation. We obtain a generalization of this equation to three dimensions. The flows considered are applicable to a wide array of physical problems where liquid films are used to lubricate higher or lower viscosity core fluids, for which a concentric arrangement is desired. Linearized solutions show that the effects of density and viscosity stratification are crucial to the stability of the interface. Rotation generally destabilizes non-axisymmetric disturbances to the interface, whereas the centripetal forces tend to stabilize flows in which the film contains the heavier fluid. Nonlinear affects allow finite amplitude helically travelling waves to exist when the fluids have different viscosities.

Coward, Aidrian V.↗
NASA NTRS2018

Collisionless Relaxation of the Ion Ring Distribution in Space Plasma

Energetic processes often produce transversely-heated angular distributions of the magnetized core (lowest energy) plasma. This characteristic is found in solar wind ion pickup, resulting from cometary or interstellar gas ionization, in Earths' ionosphere, and with hot ions formed around the Space Transportation System during gas releases. We investigate the thermalization of O+ ion pickup using the 2.5D hybrid simulation method (with fluid electrons and kinetic ions) of the ion pickup (ring) distributions, formed in the auroral ionosphere, with a range of ring velocities and thermal to magnetic pressure ratios. We find that in the unstable collisonless regime the anisotropy of the non-thermal distribution produces the ion-cyclotron instability, and the nonlinear relaxation is accompanied by wave-particle scattering that results in an emitted power of EMIC waves. We conclude that ionospheric pickup thermalization is slow due to the small ring speed compared to the thermal and Alfven speeds, while in the solar wind and other space plasmas regions with larger ion-ring velocity the collisionless relaxation and thermalization is rapid in terms of O+ ion gyro-period.

Ofman, Leon↗
DOE OSTI2020

High pressure/high temperature multiphase simulations of dodecane injection to nitrogen: Application on ECN Spray-A

The present work investigates the complex phenomena associated with pressure/high temperature dodecane injection for the Engine Combustion Network (ECN) Spray-A case, employing more elaborate thermodynamic closures, to avoid well known deficiencies concerning density and speed of sound prediction using traditional cubic models. A tabulated thermodynamic approach is proposed here, based on log 10 (p)-T tables, providing very high accuracy across a large range of pressures, spanning from 0 to 2500 bar, with only a small number of interpolation points. The tabulation approach is directly extensible to any thermodynamic model, existing or to be developed in the future. Here NIST REFPROP properties are used, combined with PC-SAFT Vapor-Liquid-Equilibrium to identify the liquid in mixtures penetration, hence avoiding the use of an arbitrary threshold for mass fraction. Identified liquid and vapour penetration are compared against experimental data from the ECN database showing a good agreement, within approximately 3–8% for axial penetration of liquid, 2% for vapor axial penetration and within experimental uncertainty for radial distribution of mass fraction. Analysis of the vortex evolution indicates that driving mechanisms behind the jet break-up are vortex tilting/stretching, then baroclinic torque, leading to Rayleigh-Taylor instabilities, closely followed by vortex dilation and finally viscous effects.

42 ENGINEERING↗
DOE OSTI2023

Drift kinetic effects on plasma response to resonant magnetic perturbation for EU DEMO design

Abstract A systematic investigation of the plasma response to an externally applied resonant magnetic perturbation (RMP) field, for the purpose of controlling edge localized modes, is carried out for an EU DEMO reference plasma. Particular emphasis is placed on the role of kinetic effects associated with both thermal particles and fusion-born alphas. The single fluid, resistive model predicts a large peak amplification of the n = 1 ( n is the toroidal mode number) plasma response to the target equilibrium, which is found to be close to the Troyon no-wall limit. A more advanced response model, including kinetic resonances between the RMP perturbation and drift motions of thermal and energetic particles, on the other hand, finds a strong suppression of the n = 1 field amplification. A major role is played by the precessional drift resonance of fusion-born alphas. A strong parallel sound wave damping model is found to well reproduce the full kinetic response results for the DEMO plasma, in terms of both the resonant field response amplitude and the plasma displacement. Finally, both fluid and kinetic models produce similar responses for the n = 2 and 3 RMP fields for the considered DEMO plasma, whilst kinetic effects again become important for the n = 4 RMP due to proximity of the reference plasma to the no-wall limit for the n = 4 ideal kink instability.

Physics↗
NASA NTRS2002

Fundamentals of Plasma Sails Propulsion Concept

The Mini-Magnetospheric Plasma Propulsion (M2P2), originally proposed by Winglee et al. [2000], is based on the two-fluid plasma model and requires a 15-km frontal standoff distance (or 20-km cross-sectional diameter) in order for the magnetic bubble to absorb sufficient momentum from the SW to accelerate a spacecraft to the unprecedented speeds of 50-80 km/s after an acceleration period of about three months. Winglee et al. [2000] derived the above size requirement based on an extrapolation of their simulated results in which a system much smaller than a M2P2 was used (p. 21,074 of their study). We submit, however, that a fluid model has no validity for such a small scale size-even in the region near the plasma source! It is assumed in the MHD fluid model, normally applied to the magnetosphere, that the characteristic scale-size is much greater than the Larmor radius and ion skin depth of the SW. In the case of the M2P2, however, the size of the magnetic bubble is actually less than or, at best, comparable to, the scale of these characteristic parameters and, therefore, a kinetic approach, which addresses the smallscale physical mechanisms involved, must be used. A fully three-dimensional version of the hybrid code is used in our M2P2 (Plasma Sails) studies was originally developed by Delamere et al. [1999]. The M2P2 plasma sail is an excellent application for this hybrid code. The primary advantage of this code is the seamless interface between fluid and kinetic descriptions of the ion populations. A kinetic description is not necessary for the dense inner regions of the magnetic bubble and tremendous computational savings can be realized by treating this dense, magnetized ion population with the fluid description. It is essential, however, that the outer bubble regions be treated kinetically as well as the SW protons. Comparison of full size M2P2 simulation based on 3D MHD and kinetic models show that kinetic treatment introduces much more asymmetry to the considering problem and the possibility of kinetic instabilities development.

Khazanov, G. V.↗
NASA NTRS2016

Cold Flow Propulsion Test Complex Pulse Testing

When the propellants in a liquid rocket engine burn, the rocket not only launches and moves in space, it causes forces that interact with the vehicle itself. When these interactions occur under specific conditions, the vehicle's structures and components can become unstable. One instability of primary concern is termed pogo (named after the movement of a pogo stick), in which the oscillations (cycling movements) cause large loads, or pressure, against the vehicle, tanks, feedlines, and engine. Marshall Space Flight Center (MSFC) has developed a unique test technology to understand and quantify the complex fluid movements and forces in a liquid rocket engine that contribute strongly to both engine and integrated vehicle performance and stability. This new test technology was established in the MSFC Cold Flow Propulsion Test Complex to allow injection and measurement of scaled propellant flows and measurement of the resulting forces at multiple locations throughout the engine.

McDougal, Kris↗
NASA NTRS1981

Effects of spacecraft motions on fluids experiments

The equations of motion governing an incompressible fluid contained in an orbiting laboratory were examined to isolate various fictitious forces and their relative influence on the fluid. The forces are divided into those arising from the orbital motions and those arising from small local motions of the spacecraft about its center of mass. The latter dominate the nonrotating experiments. Both are important for rotating experiments. A brief discussion of the onset of time-dependence and violent instability in earth-based rotating and processing systems is given.

Gans, R. F.↗
NASA NTRS1991

AIAA/IKI Microgravity Science Symposium, 1st, Moscow, USSR, May 13-17, 1991, Proceedings

The state of the art of microgravity science is examined in reviews and reports. Consideration is given to fluid physics; modeling and numerical analysis; crystal growth of electronic materials, proteins, metals, alloys, glasses, and ceramics; biotechnology, combustion science, microgravity hardware, and fundamental phenomena; and education. Particular attention is given to containerless science and technology, a low-gravity two-phase flow with heat transfer, numerical investigation on Benard instability in a finite liquid layer, synthesis and crystallization of refractory compounds from solutions in metallic melts under microgravitation conditions, the long duration protein crystallization experiment aboard the Mir Space Station, directional solidification in immiscible systems, formation of periodic layer structure in oxide melts under different gravity values, opposed flow flame spread in normal, enhanced, and reduced gravity, automatic equipment for semiconductor production in space, and consort and joust sounding rocket missions.

Source record↗
NASA NTRS1993

On the unique structure of the magnetic fields of Uranus and Neptune

The magnetic fields of Uranus and Neptune, which have comparable dipole, quadrupole, and octupole harmonics, are unique in the present-day solar system, but they resemble the geomagnetic field at the epochs of excursions and reversals known from paleomagnetic data. The precession dynamo model, in which the dominant role in the generation of the planetary magnetic fields is played by external gravitational forces, allows us to propose two scenarios for the formation of the unique topology of the magnetic fields of Uranus and Neptune. In the first case, tidal flows in the 'oceans' of these two planets extend down to the depths where the matter has a noticeable electric conductivity and velocity. A hydromagnetic interaction of the moving conducting fluid with the planetary magnetic field outside the generation region results in the deformation of the field and the deceleration of the motion under the action of the radial magnetic field. In the second case, the deformation of the field facilitates drastic changes in cyclonic cells within the generation region causing instabilities that result in a multi-polar field structure, excursions, and inversions. This paper considers this problem in greater detail by using the Neptune-Triton system as an example.

Dolginov, Sh. SH.↗
NASA NTRS1983

Post instability in continuous systems. II - Post-instability models of continua

New models simulating postinstability behavior of continuous systems are discussed. Multivalued velocity fields of the motion of a medium are analyzed, deriving concepts which are applied to a new model of turbulence. Postinstability models of solids are examined, emphasizing the multivaluedness of the parameters. The consequences are addressed of the simplifying assumption that in solids the rate of dissipation of energy of the pulsations is much higher than in fluids. The principal stress formulation of postinstability models is considered, including the geometry, kinematics, and dynamics of the principal stress coordinates. The fundamental properties of enlarged models of a solid are illustrated by the theory of wrinkling films. Results are presented for the statics of such films, the effects of shocks and shock waves, and the behavior of singular effects and singular thick wrinkles.

Zak, M.↗
NASA NTRS1984

On the stability of the ionopause of Venus

The stability of the Venus ionopause is examined in light of the importance of gravitation and curvature. Using a one-fluid approximation for the equation of motion of the plasma, and ignoring the effects of neutrals, a dispersion relation is obtained that includes the effects of the magnetic field, sheared plasma flow, buoyancy, centrifugal force and magnetic tension due to boundary curvature. It is found that buoyancy acts to neutralize the flute instability. As expected, the Kelvin-Helmholtz mode is the dominant instability over most of the dayside ionopause. The expected growth times of this mode are short in comparison with the wave-convection time over the boundary; the waves can grow and saturate quickly, producing a turbulent boundary that may affect electrodynamic coupling between the solar wind and ionospheric plasmas.

Elphic, R. C.↗
NASA NTRS1993

Mixing and unmixedness in plasma jets 1: Near-field analysis

The flow characteristics in the near-field of a plasma jet are simulated with a two-fluid model. This model accounts for both gradient-diffusion mixing and uni-directional sifting motion resulting from pressure-gradient-body-force imbalance. This latter mechanism is believed to be responsible for the umixedness observed in plasma jets. The unmixedness is considered to be essentially a Rayleigh-Taylor kind instability. Transport equations are solved for the individual plasma and ambient gas velocities, temperatures and volume fractions. Empirical relations are employed for the interface transfers of mass, momentum and heat. The empirical coefficients are first established by comparison of predictions with available experimental data for shear flows. The model is then applied to an Argon plasma jet ejecting into stagnant air. The predicted results show the significant build-up of unmixed air within the plasma gas, even relatively far downstream of the torch. By adjusting the inlet condition, the model adequately reproduces the experimental data.

Ilegbusi, Olusegun J.↗
DOE OSTI2023

Predicting Wind Loading and Instability in Solar Tracking PV Arrays

Wind loading and the fluctuating pressure loads it creates on PV panel surfaces are associated with multiple degradation mechanisms and failures. Modest wind speeds create reversing loads that can initiate cell cracks and weather cracked cells. Stronger wind speeds and extreme weather events can lead to larger scale forces and the aerodynamic instability known as torsional galloping. All these effects are dependent on the complex coupling between wind speed, panel orientation, and a myriad of other hardware and site-specific factors. In this work, we present the latest developments from our work to build an open-source, high-performance computing (HPC) fluid dynamics solver to predict and mitigate these effects. This simulation package allows users to easily specify different array layouts, solar-tracking angles, panel geometries, and weather conditions before automatically generating a refined computational mesh and solving for the unsteady loading on each panel surface. Small domains (e.g., a single panel row in isolation) can be solved on a modern laptop, while larger domains or very high-fidelity studies can be solved on distributed or HPC resources with minimal modifications to the underlying problem specification. We present preliminary case studies obtained using this simulation package and highlight how increased wind speeds combined with sub-optimal tracking angles can exacerbate degradation drivers.

aerodynamics↗
DOE OSTI2024

The Onset of Nucleate Boiling Prediction in a Mini Channel with the Eulerian Multiphase Flow

The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing is developed using a general buoyancy–shear–drag model. Analytic solutions to the simplified model equations corresponding to each instability are derived separately, which are then used to calibrate the model coefficients to predict specific values of the mixing layer growth parameters and exponents. The buoyancy-shear–drag equations for the bubble and spike mixing layer widths are then solved numerically, and a turbulent diffusivity (or viscosity) is constructed dimensionally as the product of the mixing layer width, h, and its time-derivative, dh/dt. Surrogate turbulent fields are then constructed by multiplying a presumed, approximate self-similar spatial profile by appropriate functions of h and dh/dt. Using several simplifying approximations, the turbulent diffusion equations satisfied by the mean mass fraction and mean shear velocity are solved analytically in a reference frame moving with the mean advection velocity, and it is shown that the mean fields evolve in space and time as expected. The explicit modeling and solution of turbulent transport equations (e.g., for the turbulent kinetic energy and its dissipation rate or a lengthscale) are not required. By using separate equations for the bubble and spike evolution, the model naturally captures the asymmetry induced by an increasing density contrast between the heavy and light fluids. The model also includes molecular dissipation and diffusion, and therefore, can describe transition to fully-developed turbulence if the initial turbulent diffusivity/viscosity is much smaller than the molecular diffusivity/viscosity. This model is applied to constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and variable-density Kelvin–Helmholtz mixing layers to demonstrate its utility. It is shown that the numerical solutions of the model calibrated using specific values of the instability growth parameters and exponents: (1) produces mixing layer widths in agreement with the expected self-similar growth power laws; (2) gives spatiotemporal profiles of turbulent fields that are expected and consistent with previous results; (3) predicts the expected power-law growths and decays of the spatially-integrated turbulent fields, and; (4) gives spatiotemporal profiles of the mean fields that are expected and consistent with previous results. This relatively simple model can be implemented numerically, and is expected to be useful for fundamental as well as astrophysical and inertial confinement fusion applications.

42 ENGINEERING↗
NASA NTRS1990

Ground-based experiments on thermal and thermosolutal convection in inclined low-aspect-ratio enclosures

Thermal and thermosolutal convection driven by buoyancy forces are investigated experimentally in inclined low-aspect-ratio enclosures. Thermal and thermosolutal convection are important in crystal growth processes. The transport process in the fluid phase during the growth of a crystal has a profound influence on the structure and quality of the solid phase. The present work is concerned with natural convection in a shallow rectangular enclosure with imposed inclined temperature and/or concentration gradients. Various complex flow patterns are observed with different experimental conditions. A flow instability under certain conditions is reported. The flow instability is found to depend not only on the thermal and solutal Grashof numbers, the aspect ratios, but also the inclined angles.

Wang, L. W.↗
NASA NTRS1992

The effect of gravity modulation on thermosolutal convection in an infinite layer of fluid

The effect of time-periodic vertical gravity modulation on the onset of thermosolutal convection in an infinite horizontal layer with stress-free boundaries is investigated using Floquet theory for the linear stability analysis. Situations for which the fluid layer is stably stratified in either the fingering or diffusive regimes of double-diffusive convection are considered. Results are presented both with and without steady background acceleration. Modulation may stabilize an unstable base solution or destabilize a stable base solution. In addition to synchronous and subharmonic response to the modulation frequency, instability in the double diffusive system can occur via a complex conjugate mode. In the diffusive regime, where oscillatory onset occurs in the unmodulated system, regions of resonant instability occur and exhibit strong coupling with the unmodulated oscillatory frequency. The response to modulation of the fundamental instability of the unmodulated system is described both analytically and numerically; in the double-diffusive system this mode persists under subcritical conditions as a high-frequency lobe.

Saunders, B. V.↗
NASA NTRS1973

Turbulent heating of colliding streams in the solar wind.

Turbulent heating of colliding plasma streams has previously been observed in the solar wind. The original data were interpreted in terms of a fluid model. It is contended that a plasma-kinetic description is the more appropriate theoretical approach and is necessary in order to better understand the microscopic physical phenomena that underlie all fluid models. Microscopic solar-wind parameters characteristic of conditions during the observations were used, together with the quasi-linear plasma-kinetic theory, to compute the expected magnetic field and temperature enhancements in the interaction region between two counterstreaming plasma beams. The physical mechanism of excitation is the electromagnetic two-stream instability in which Alfven waves are unstable. A total field in the interaction region of about 8 gamma and a change in temperature of about 100,000 K are obtained.

Goldstein, M. L.↗