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 1,009 records · Page 56

All sources NASA NTRS DOE OSTI

NASA NTRS1996

CFD Analysis of the 24-inch JIRAD Hybrid Rocket Motor

A series of multispecies, multiphase computational fluid dynamics (CFD) analyses of the 24-inch diameter joint government industry industrial research and development (JIRAD) hybrid rocket motor is described. The 24-inch JIRAD hybrid motor operates by injection of liquid oxygen (LOX) into a vaporization plenum chamber upstream of ports in the hydroxyl-terminated polybutadiene (HTPB) solid fuel. The injector spray pattern had a strong influence on combustion stability of the JIRAD motor so a CFD study was initiated to define the injector end flow field under different oxidizer spray patterns and operating conditions. By using CFD to gain a clear picture of the flow field and temperature distribution within the JIRAD motor, it is hoped that the fundamental mechanisms of hybrid combustion instability may be identified and then suppressed by simple alterations to the oxidizer injection parameters such as injection angle and velocity. The simulations in this study were carried out using the General Algorithm for Analysis of Combustion SYstems (GALACSY) multiphase combustion codes. GALACSY consists of a comprehensive set of droplet dynamic submodels (atomization, evaporation, etc.) and a computationally efficient hydrocarbon chemistry package built around a robust Navier-Stokes solver optimized for low Mach number flows. Lagrangian tracking of dispersed particles describes a closely coupled spray phase. The CFD cases described in this paper represent various levels of simplification of the problem. They include: (A) gaseous oxygen with combusting fuel vapor blowing off the walls at various oxidizer injection angles and velocities, (B) gaseous oxygen with combusting fuel vapor blowing off the walls, and (C) liquid oxygen with combusting fuel vapor blowing off the walls. The study used an axisymmetric model and the results indicate that the injector design significantly effects the flow field in the injector end of the motor. Markedly different recirculation patterns are observed in the vaporization chamber as the oxygen velocity and/or spray pattern is varied. The ability of these recirculation patterns to stabilize the diffusion flame above the surface of the solid fuel gives a plausible explanation for the experimentally determined combustion stability characteristics of the JIRAD motor, and suggests how combustion stability can be assured by modifications to the injector design.

Liang, Pak-Yan↗

NASA NTRS1982

Three-dimensional stability of vortex arrays

The stability to three-dimensional disturbances of three classical steady vortex configurations in an incompressible inviscid fluid is studied in the limit of small vortex cross-sectional area and long axial disturbance wavelength. The configurations examined are the single infinite vortex row, the Karman vortex street of staggered vortices and the symmetric vortex street. It is shown that the single row is most unstable to a two-dimensional disturbance, while the Karman vortex street is most unstable to a three-dimensional disturbance over a significant range of street spacing ratios. The symmetric vortex street is found to be most unstable to three-dimensional or two-dimensional symmetric disturbances depending on the spacing ratio of the street. Short remarks are made concerning the relevance of the calculations to the observed instabilities in free shear layer, wake and boundary-layer type flows.

Robinson, A. C.↗

DOE OSTI2022

Cross-beam energy transfer in direct-drive ICF. I. Nonlinear and kinetic effects

Results are presented from a series of simulations examining the susceptibility of the cross-beam energy transfer (CBET) instability to nonlinear processes in the context of direct-drive inertial confinement fusion experiments on the OMEGA laser facility. These form the basis for the second paper of this series [A. G. Seaton, L. Yin, R. Follett, B. J. Albright, and A. Le, “Cross-beam energy transfer in direct-drive ICF. II. Theory and simulation of mitigation through increased laser bandwidth,” Phys. Plasmas 29, 042707 (2022)], where we examine the efficacy of increases in laser bandwidth at suppressing CBET. We choose laser and plasma conditions for the simulations that are favorable to CBET and promote nonlinearity. Through a comparison of outputs from the particle-in-cell code vector particle in cell (VPIC) and the linearized fluid code laser-plasma simulation environment (LPSE), a series of nonlinear effects have been identified in the kinetic simulations that include particle trapping, the two-ion wave decay, and ion-acoustic wave self-focusing. These effects produce time-dependent energy transfer, in contrast to the linearized fluid simulations in which a steady state is reached after an initial transient. Ion trapping is shown to allow for increased energy transfer relative to fluid simulations, with the remaining nonlinear processes acting to reduce the energy transfer. Nonlinear dynamics is contrasted for low- and high-intensity beams as well as between speckled and planar beams. For the parameters under consideration, beam profile has a significant effect on nonlinear dynamics, though the greatest sensitivity is to beam intensity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1992

Perspectives on hypersonic viscous and nonequilibrium flow research

An attempt is made to reflect on current focuses in certain areas of hypersonic flow research by examining recent works and their issues. Aspects of viscous interaction, flow instability, and nonequilibrium aerothermodynamics pertaining to theoretical interest are focused upon. The field is a diverse one, and many exciting works may have either escaped the writer's notice or been abandoned for the sake of space. Students of hypersonic viscous flow must face the transition problems towards the two opposite ends of the Reynolds or Knudsen number range, which represents two regimes where unresolved fluid/gas dynamic problems abound. Central to the hypersonic flow studies is high-temperature physical gas dynamics; here, a number of issues on modelling the intermolecular potentials and inelastic collisions remain the obstacles to quantitative predictions. Research in combustion and scramjet propulsion will certainly be benefitted by advances in turbulent mixing and new computational fluid dynamics (CFD) strategies on multi-scaled complex reactions. Even for the sake of theoretical development, the lack of pertinent experimental data in the right energy and density ranges is believed to be among the major obstacles to progress in aerothermodynamic research for hypersonic flight. To enable laboratory simulation of nonequilibrium effects anticipated for transatmospheric flight, facilities capable of generating high enthalpy flow at density levels higher than in existing laboratories are needed (Hornung 1988). A new free-piston shock tunnel capable of realizing a test-section stagnation temperature of 10(exp 5) at Reynolds number 50 x 10(exp 6)/cm is being completed and preliminary tests has begun (H. Hornung et al. 1992). Another laboratory study worthy of note as well as theoretical support is the nonequilibrium flow experiment of iodine vapor which has low activation energies for vibrational excitation and dissociation, and can be studied in a laboratory with modest resources (Pham-Van-Diep et al. 1992).

Cheng, H. K.↗

DOE OSTI2021

Modification of favorable average curvature effect by changing parallel sound wave behavior in tokamak plasmas

Abstract The favorable average curvature effect, also known as the GGJ effect (Glasser et al 1975 Phys. Fluids 18 875), is intrinsically associated with parallel sound wave propagation in a tokamak plasma. This work investigates how the GGJ effect is modified by changing the parallel sound wave behavior. Two physics models beyond the standard single fluid theory, i.e. an anisotropic thermal transport model and a parallel sound wave damping model, are employed to change parallel sound waves in a toroidal plasma, and the consequence on the GGJ effect is demonstrated for two important classes of problems, i.e. the resistive plasma response to the applied resonant magnetic perturbation and the stability of the tearing mode. Toroidal modeling reveals that the GGJ effect is significantly altered by both of the aforementioned physics effects. Compared to the thermal transport physics, which completely removes the GGJ effect, the sound wave damping effect only offers partial mitigation. The differences between these two models are further illustrated in terms of the radial structure of the shielding current and the eigenfunction of the tearing instability. In particular, a fundamental reason for complete suppression of the GGJ effect by the thermal transport is identified as an extra toroidal coupling of the poloidal harmonics.

Physics↗

DOE OSTI2026

Validation of the GFS model for gyrokinetic stability of NSTX pedestal data

This study presents a large database validation of the gyro fluid system (GFS) model for linear gyrokinetic stability for high-mode (H-mode) edge transport barrier conditions in the national spherical torus experiment (NSTX) tokamak. The database of linear stability calculations with the CGYRO gyrokinetic code was produced using plasma profile measurements from NSTX discharges to identify kinetic ballooning modes (KBM), trapped electron modes (TEM), and micro-tearing modes (MTM) that limit the pressure profile gradient in the H-mode barrier. A novel Bayesian optimization approach determines optimal resolution parameters for GFS specifically for spherical tokamak pedestal conditions. Our results demonstrate that GFS, with optimized resolution, can achieve accurate linear stability analysis in NSTX pedestal conditions for reduced resolution compared to CGYRO. GFS can accurately find the KBM, TEM, and MTM instability branches. Parametric analysis reveals that GFS accuracy in this extreme pedestal parameter range is degraded for low magnetic shear and near the separatrix conditions. These findings establish GFS as a fast linear eigenmode solver for spherical tokamak pedestal gyrokinetic stability and demonstrate a systematic methodology for determining the optimum resolution settings.

Yang, Minglei [Oak Ridge National Laboratory (ORNL↗

NASA NTRS1984

An analysis of the stability of the compressible Ekman boundary layer

The linear stability problem for the compressible Ekman boundary layer common to rotating fluids is formulated and the stability properties determined numerically. Three classes of unstable waves are identified (called class A, B, and C), their properties are described. The class C waves have only recently been reported in the literature and are present only in compressible Ekman boundary layers. Most of the calculations presented here are for uranium hexafluoride gas; however, critical Reynolds numbers are also computed for air and ammonia gas. Compressibility is generally found to decrease the critical Reynolds number for each class of wave. A comparison of results for the three different gases shows the stability to be largely unaffected by changes in the gas properties. Maximum growth rate calculations for each wave show the class A and B waves to be the dominant instabilities.

Spall, J. R.↗

NASA NTRS1987

Straight cylindrical seal for high-performance turbomachines

A straight cylindrical seal configuration representing the seal for a high-performance turbopump (e.g., the space shuttle main engine fuel pump) was tested under static (nonrotating) conditions. The test data included critical mass flux and pressure profiles over a wide range of inlet temperatures and pressures for fluid nitrogen and fluid hydrogen with the seal in concentric and fully eccentric positions. The critical mass fluxes (or leakage rates) for the concentric and fully eccentric configurations were nearly the same when based on stagnation conditions upstream of the seal. The fully eccentric configuration pressure profiles of the gas and liquid were different. Further, the pressure differences between the maximum and the minimum clearance positions were highly dependent on the geometric conditions, the temperature, and the absolute pressure at both the inlet and the exit. The pressure differences were greatest in the inlet region. The results, although complex, tend to follow the corresponding-states principles for critical flows. Gaseous injection near the seal exit plane significantly altered the pressure profiles and could be used to control turbomachine instabilities.

Hendricks, Robert C.↗

NASA NTRS1991

Simulation of high-speed (orbital) releases of electron attachment materials in the ionosphere

The dynamics of ionospheric plasma irregularities produced by the release of electron attachment materials at orbital velocities across the geomagnetic field is studied. A two-dimensional electrostatic fluid model which includes electron attachment and mutual neutralization chemistry, self-consistent electric fields, and three-species transport is developed. Numerical simulations are performed to study the behavior at early and at late times after the release. At early times, of the order of or less than the attachment material neutral collision time, the negative ion cloud produced by the release may structure owing to the shear in the E x B velocity within the cloud. At high altitudes the cloud may bifurcate and form vortices on the back. At lower altitudes where ion-neutral collisional effects dominate, this structuring is suppressed. At late times, after a plasma depletion has formed due to neutralization chemistry, the cloud structures by the E x B interchange instability. Depending on the release altitude, the depletion structures by the collisional or inertial limit of this instability.

Scales, W. A.↗

DOE OSTI2022

Higher Harmonics in Multipactor Induced Plasma Ionization Breakdown near a Dielectric Surface

In this Letter, the novel physics of higher harmonic (HH) generation of the normal electric field near a dielectric surface is reported in multipactor induced plasma ionization breakdown, as determined by kinetic particle-in-cell simulations. The observed HH frequency is around ten times the fundamental rf driving frequency, but lower than the electron plasma frequency. A theory is constructed which indicates that stream plasma interaction-induced instability is the mechanism of HH generation in the collisional regime. The HH frequency and its corresponding growth rate of the HH oscillation amplitude from the theory are in good agreement with kinetic particle-in-cell simulations.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2021

Simulations of ion heating due to ion-acoustic instabilities in presheaths

Particle-in-cell, direct simulation Monte Carlo simulations reveal that ion-acoustic instabilities excited in presheaths can cause significant ion heating. Ion-acoustic instabilities are excited by the ion flow toward a sheath when the neutral gas pressure is small enough and the electron temperature is large enough. A series of 1D simulations were conducted in which neutral plasma (electrons and ions) was uniformly sourced with an ion temperature of 0.026 eV and different electron temperatures (0.1 eV–50 eV). Ion heating was observed when the electron-to-ion temperature ratio exceeded the minimum value predicted by linear response theory to excite ion-acoustic instabilities at the sheath edge (T e /T i ≈ 28). When this threshold was exceeded, the temperature equilibration rate between ions and electrons rapidly increased near the sheath so that the local temperature ratio did not significantly exceed the threshold for instability. This resulted in significant ion heating near the sheath edge, which also extended back into the bulk plasma; presumably due to wave reflection from the sheath. Furthermore, this ion-acoustic wave heating mechanism was found to decrease for higher neutral pressures, where ion-neutral collisions damp the ion-acoustic waves and ion heating is instead dominated by inelastic collisions in the presheath.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1993

Analytical description of the breakup of liquid jets in air

A viscous or inviscid cylindrical jet with surface tension in a vacuum tends to pinch due to the mechanism of capillary instability. Similarity solutions are constructed which describe this phenomenon as a critical time is encountered, for two physically distinct cases: inviscid jets governed by the Euler equations and highly viscous jets governed by the Stokes equations. In both cases the only assumption imposed is that at the time of pinching the jet shape has a radial length scale which is smaller than the axial length scale. For the inviscid case, we show that our solution corresponds exactly to one member of the one-parameter family of solutions obtained from slender jet theories and the shape of the jet is locally concave at breakup. For highly viscous jets our theory predicts local shapes which are monotonic increasing or decreasing indicating the formation of a mother drop connected to the jet by a thin fluid tube. This qualitative behavior is in complete agreement with both direct numerical simulations and experimental observations.

Papageorgiou, Demetrios T.↗

NASA NTRS2013

Characterization of Aerodynamic Interactions with the Mars Science Laboratory Reaction Control System Using Computation and Experiment

On August 5 , 2012, the Mars Science Laboratory (MSL) entry capsule successfully entered Mars' atmosphere and landed the Curiosity rover in Gale Crater. The capsule used a reaction control system (RCS) consisting of four pairs of hydrazine thrusters to fly a guided entry. The RCS provided bank control to fly along a flight path commanded by an onboard computer and also damped unwanted rates due to atmospheric disturbances and any dynamic instabilities of the capsule. A preliminary assessment of the MSL's flight data from entry showed that the capsule flew much as predicted. This paper will describe how the MSL aerodynamics team used engineering analyses, computational codes and wind tunnel testing in concert to develop the RCS system and certify it for flight. Over the course of MSL's development, the RCS configuration underwent a number of design iterations to accommodate mechanical constraints, aeroheating concerns and excessive aero/RCS interactions. A brief overview of the MSL RCS configuration design evolution is provided. Then, a brief description is presented of how the computational predictions of RCS jet interactions were validated. The primary work to certify that the RCS interactions were acceptable for flight was centered on validating computational predictions at hypersonic speeds. A comparison of computational fluid dynamics (CFD) predictions to wind tunnel force and moment data gathered in the NASA Langley 31-Inch Mach 10 Tunnel was the lynch pin to validating the CFD codes used to predict aero/RCS interactions. Using the CFD predictions and experimental data, an interaction model was developed for Monte Carlo analyses using 6-degree-of-freedom trajectory simulation. The interaction model used in the flight simulation is presented.

Schoenenberger, Mark↗

NASA NTRS2022

On the Stability and Phase Behavior of Titan's Subsurface Liquid Columns

On Titan, liquid hydrocarbon may stay in the subsurface porous permeable crust known as the alkanofer, analogous to water in Earth’s aquifer. In addition to pressure gradient, the subsurface liquid in alkanofers is subject to vertical compositional grading due to the gravity and temperature gradient. The common wisdom is that the liquid would normally stay underground in a stability established by the pressure that increases with depth as observed in aquifers on Earth. However, Titan’s liquids consist of nitrogen and hydrocarbons, mainly methane and ethane, the behavior of which is very sensitive to temperature and pressure. Consequently, the liquid density does not always increase with depth, thus may introduce a reverse density profile that leads to vertical convective instability of the liquid column. If reverse density profiles are present, capillary pressures arising from liquid trapped within small pores in the crust can help with the column stability. The liquid held in the capillaries can seal the space below it thus helping with the stability, unless the overpressure built from underneath becomes larger than the capillary pressure, which causes leakage to allow the liquids to seep upward from the deep. This situation is analogous to hydrocarbon seeps on Earth, where oil and natural gas escape the reservoir and flow slowly through network of cracks to the surface. An algorithm based on an extended Gibbs equation commonly used in petroleum reservoir engineering is employed in this work to produce pressure, density, and compositional profiles for evaluating the stability and phase behavior of Titan’s subsurface fluids.

Sugata Tan↗

DOE OSTI2022

Coupling a Lagrangian–Eulerian Spark-Ignition (LESI) model with LES combustion models for engine simulations

In the United States transportation sector, Light-Duty Vehicles (LDVs) are the largest energy consumers and CO 2 emitters. Electrification of LDVs is posed as a potential solution, but SI engines can still contribute to decarbonization. Car manufacturers have turned to unconventional engine operation to increase the efficiency of Spark-Ignition (SI) engines and reduce the carbon emissions of their fleets. Dilute, lean, and stratified-charge engine operation has the potential for engine efficiency improvements at the expense of increased cyclic variability and combustion instability. At such demanding engine conditions, the spark ignition event is key for flame initiation and propagation and for enhanced combustion stability. Reliable and accurate spark ignition models can help design ignition systems that reduce cyclic variability. Multiple computational spark-ignition models exist that perform well under conventional conditions, but the underlying physics needs to be expanded, for unconventional engine operation. In this paper, a hybrid Lagrangian–Eulerian Spark-Ignition (LESI) model is coupled with different turbulent flame propagation models for engine simulations. LESI relies on Lagrangian arc tracking and Eulerian energy deposition. The LESI model is coupled with the Well-Stirred Reactor (WSR), Thickened Flame Model (TFM), and g-equation model and used to simulate several cycles of a Direct-Injection Spark-Ignition (DISI) engine using a commercial Computational Fluid Dynamics (CFD) engine solver. The results showcase the successful coupling of LESI with the combustion models. Global engine metrics, such as pressure and Apparent Heat Release Rate (AHRR), for each simulation setup are compared to experimental engine results, for validation. In addition, results highlight the successful prediction of spark channel movement by comparing simulation images to experimental optical engine images. Finally, the successful coupling of LESI to combustion models, making it a usable model in the engine modeling community, is emphasized and future development details are discussed.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2024

Core electron temperature turbulence and transport during sawtooth oscillations in the DIII-D tokamak

Abstract Sawteeth are one of the concerning instabilities in ITER and future burning plasma experiments. Sawtooth dynamics and its interaction with broadband plasma turbulence has been a challenge for predictive simulations of core transport in future fusion devices. This study provides new observations of core turbulence behavior during sawtooth oscillations in DIII-D hydrogen L-mode neutral beam injection heated plasmas in an inner wall limited configuration. A strong correlation of electron temperature and density turbulence levels with the sawtooth oscillation phase has been observed at locations inside the T e inversion radius and/or safety factor q = 1 magnetic surface. The T e turbulence amplitude in the core during the sawtooth ramp exhibits a critical T e gradient behavior inside but not near the T e inversion radius/ q = 1 magnetic surface. The most unstable mode calculated from the trapped gyro-landau fluid turbulence simulations reveal a change from low- k ion-type to low-k electron-type modes from pre- to post- sawtooth crash time periods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Operando synchrotron X-ray analysis of melt pool dynamics in an Al-Sn immiscible alloy

The melt flow in an Al-50vol.% Sn immiscible alloy, produced by single-track laser melting of Al and Sn elemental powders, was studied in real time. High-speed synchrotron X-ray imaging was used to track Al and Sn liquids’ movements, and also to examine elemental distributions in the laser tracks, complimented by electron microscopy after solidification. Key aspects, including melt pool geometry, keyhole instability, and flow dynamics (flow pattern and velocity), were examined using digital image analysis. Relatively deeper melt pools formed at 400 W and 300 mm/s exhibited greater stability, with smooth surfaces, consistent outward flow, and minor vortices near the keyhole. In contrast, shallower pools produced at higher scanning speeds (>500 mm/s) demonstrated greater instability with increased surface waviness, and stronger velocity fluctuations, leading to numerous micro-vortices and increased Al-Sn heterogeneity. Velocity scale estimations, supported by experimental observations, examined the roles of vapour pressure, Marangoni effect, buoyancy, inertial, and surface tension forces in the flow. The results revealed that vapour pressure and mechanical waves dominated at high scanning speeds (shallow pools), while Marangoni forces were equally significant in deep pools at lower speeds (300 mm/s). Buoyancy was found to have minimal impact in both cases. Furthermore, the interaction between inertial and surface tension forces played a critical role in determining the degree of waviness of the pools’ surfaces. These findings offer valuable insights into melt pool dynamics during laser processing of immiscible alloys and other metallic systems using elemental powders, and provide guidance for developing high-fidelity computational fluid dynamics models.

immiscible alloys↗

DOE OSTI2021

Probing nonlinear velocity profiles of shear-thinning, nematic platelet dispersions in Couette flow using x-ray photon correlation spectroscopy

In this study, we report experiments employing x-ray photon correlation spectroscopy (XPCS) to characterize the velocity profiles of complex fluids in Couette flow. The approach involves modeling the XPCS correlation functions obtained with the incident x-ray beam passing tangentially through the Couette cell gap at various distances from the inner wall. We first demonstrate the technique with measurements on a dilute colloidal dispersion in the Newtonian liquid glycerol, where the expected linear velocity profiles are recovered. We then employ the technique to map the shear-rate-dependent velocity profiles of a shear-thinning dispersion of nematically ordered Gibbsite platelets. The nonlinear velocity profiles of the Gibbsite dispersion include a narrow slip region adjacent to the outer wall and a band with a small velocity gradient in the interior of the gap that evolves into a region increasingly resembling plug flow with increasing shear rate. Variations in the velocity profile along the vorticity direction indicate an instability in the interface between this region of small velocity gradient and a region of high velocity gradient near the inner wall. The analysis of the small-angle scattering patterns provides information about the spatial and temporal variations in the nematic order of the Gibbsite dispersion and their coupling to the velocity profile. Additional potential applications of this XPCS-based technique and comparisons with established methods for characterizing velocity profiles are discussed.

42 ENGINEERING↗