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 1,063 records · Page 59

Magnetic flux emergence in binary neutron star remnants

Using high-resolution AthenaK simulations of a twisted toroidal flux tube, we study the flux emergence of magnetic structures in the shear layer of a hot massive neutron star typical of a binary neutron star remnant. High-resolution simulations demonstrate that magnetic buoyant instabilities allow for emergence only for extremely large magnetic fields significantly exceeding 10 17 G, and more typical fields around 10 16 G are instead dominated by hydrodynamic effects. Because merger remnants tend to be stable against hydrodynamic convection, our work places strong limitations on the mechanisms by which massive binary neutron star remnants can produce the magnetically-driven outflows needed to power jets.

79 ASTRONOMY AND ASTROPHYSICS↗

Predictive Evaluations of Oxygen-Rich Hydrocarbon Combustion Gas-Centered Swirl Coaxial Injectors using a Flamelet-Based 3-D CFD Simulation Approach

NASA Marshall Space Flight Center (MSFC) has embarked upon a joint project with the Air Force to improve the state-of-the-art of space application combustion device design and operational understanding. One goal of the project is to design, build and hot-fire test a 40,000 pound-thrust Oxygen/Rocket Propellant-2 (RP-2) Oxygen-Rich staged engine at MSFC. The overall project goals afford the opportunity to test multiple different injector designs and experimentally evaluate the any effect on the engine performance and combustion dynamics. To maximize the available test resources and benefits, pre-test, combusting flow, Computational Fluid Dynamics (CFD) analysis was performed on the individual injectors to guide the design. The results of the CFD analysis were used to design the injectors for specific, targeted fluid dynamic features and the analysis results also provided some predictive input for acoustic and thermal analysis of the main Thrust Chamber Assembly (TCA). MSFC has developed and demonstrated the ability to utilize a computationally efficient, flamelet-based combustion model to guide the pre-test design of single-element Gas Centered Swirl Coaxial (GCSC) injectors. Previous, Oxygen/RP-2 simulation models utilizing the Loci-STREAM flow solver, were validated using single injector test data from the EC-1 Air Force test facility. The simulation effort herein is an extension of the validated, CFD driven, single-injector design approach applied to single injectors which will be part of a larger engine array. Time-accurate, Three-Dimensional, CFD simulations were performed for five different classes of injector geometries. Simulations were performed to guide the design of the injector to achieve a variety of intended performance goals. For example, two GCSC injectors were designed to achieve stable hydrodynamic behavior of the propellant circuits while providing the largest thermal margin possible within the design envelope. While another injector was designed to purposefully create a hydrodynamic instability in the fuel supply circuit as predicted by the CFD analysis. Future multi-injector analysis and testing will indicate what if any changes occur in the predicted behavior for the single-element injector when the same injector geometry is placed in a multi-element array.

Richardson, Brian R.↗

Fractional Laplacian spectral approach to turbulence in a dusty plasma monolayer

This work presents an analytical investigation of anomalous diffusion and turbulence in a dusty plasma monolayer, where energy transport across scales leads to the spontaneous formation of spatially disordered patterns. Many-body simulations of 10,000-particle dusty plasma monolayers are used to demonstrate how the global dynamics depend on the statistical properties of the dust assembly for realistic laboratory conditions. We find that disorder due to variations in the dust size distribution and charge-driven nonlocal interactions resulting in anomalous dust diffusion are key factors for the onset of instabilities. The resulting dynamics exhibit features of inertial turbulence over slightly more than half a decade of scales proportional or smaller than the Debye shielding length. These processes are examined analytically using a recently developed Fractional Laplacian Spectral (FLS) technique, which identifies the active energy channels as a function of scale, disorder concentration, and features of the nonlocal-interactions. The predictions from the theoretical (spectral) analysis demonstrate agreement with the results from the many-body (kinetic) simulations, thus providing a powerful tool for the study of active turbulence.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the nonlinear stability of the unsteady, viscous flow of an incompressible fluid in a curved pipe

The stability of the flow of an incompressible, viscous fluid through a pipe of circular cross-section curved about a central axis is investigated in a weakly nonlinear regime. A sinusoidal pressure gradient with zero mean is imposed, acting along the pipe. A WKBJ perturbation solution is constructed, taking into account the need for an inner solution in the vicinity of the outer bend, which is obtained by identifying the saddle point of the Taylor number in the complex plane of the cross-sectional angle co-ordinate. The equation governing the nonlinear evolution of the leading order vortex amplitude is thus determined. The stability analysis of this flow to periodic disturbances leads to a partial differential system dependent on three variables, and since the differential operators in this system are periodic in time, Floquet theory may be applied to reduce this system to a coupled infinite system of ordinary differential equations, together with homogeneous uncoupled boundary conditions. The eigenvalues of this system are calculated numerically to predict a critical Taylor number consistent with the analysis of Papageorgiou. A discussion of how nonlinear effects alter the linear stability analysis is also given, and the nature of the instability determined.

Shortis, Trudi A.↗

Toroidal modeling of anisotropic thermal transport and energetic particle effects on stability of resistive plasma resistive wall mode

Effects of anisotropic thermal transport on the linear stability of the resistive plasma resistive wall mode (RPRWM) are investigated by the magnetohydrodynamic-kinetic hybrid code MARS-K [Liu et al., Phys. Plasmas 15, 112503 (2008)], including the kinetic contribution from energetic particles (EPs). It is found that thermal transport can further stabilize the RPRWM in the presence of drift kinetic contributions from EPs. This is different from the fluid model, which always predicts destabilization of the mode by thermal transport. Furthermore, the thermal transport effect is found to amplify the role played by an adiabatic term, associated with the radial distribution of EPs’ birth energy, in modifying the mode stability as well as the mode eigenfunction. The shape of the equilibrium profiles of EPs, in particular that of the temperature, is also found to strongly affect the mode stabilization. This profile effect is more pronounced in the absence of thermal transport. MARS-K computations show that the stabilizing effect by thermal transport is more likely to occur at slower plasma rotation and lower EP energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Investigation of Dynamic Oxygen Adsorption in Molten Solder Jetting Technology

Surface tension forces play a critical role in fluid dynamic phenomena that are important in materials processing. The surface tension of liquid metals has been shown to be very susceptible to small amounts of adsorbed oxygen. Consequently, the kinetics of oxygen adsorption can influence the capillary breakup of liquid-metal jets targeted for use in electronics assembly applications, where low-melting-point metals (such as tin-containing solders) are utilized as an attachment material for mounting of electronic components to substrates. By interpreting values of surface tension measured at various surface ages, adsorption and diffusion rates of oxygen on the surface of the melt can be estimated. This research program investigates the adsorption kinetics of oxygen on the surface of an atomizing molten-metal jet. A novel oscillating capillary jet method has been developed for the measurement of dynamic surface tension of liquids, and in particular, metal melts which are susceptible to rapid surface degradation caused by oxygen adsorption. The experimental technique captures the evolution of jet swells and necks continuously along the jet propagation axis and is used in conjunction with an existing linear, axisymmetric, constant-property model to determine the variation of the instability growth rate, and, in turn, surface tension of the liquid as a function of surface age measured from the exit orifice. The conditions investigated so far focus on a time window of 2-4ms from the jet orifice. The surface properties of the eutectic 63%Sn-37%Pb solder alloy have been investigated in terms of their variation due to O2 adsorption from a N2 atmosphere containing controlled amounts of oxygen (from 8 ppm to 1000 ppm). The method performed well for situations where the oxygen adsorption was low in that time window. The value of surface tension for the 63Sn-37Pb solder in pure nitrogen was found to be 0.49 N/m, in good agreement with previously published work. A characteristic time of O(1ms) or less was determined for the molten-metal surface to be saturated by oxygen at 1000 ppm concentration in N2.

Megaridis, Constantine M.↗

Linear stability analysis of two fluid columns of different densities and viscosities in a gravity field

The linear stability of a vertical interface separating two miscible fluid columns of different densities and viscosities under the influence of gravity is investigated. This flow possesses a time-dependent reference state (each column accelerates at different rates owing to their different densities) and the interface thickness grows as the square root of time (by diffusion). Numerical integration of the linear initial-value problem is carried out and discussed in detail as a function of vertical and spanwise wavenumbers and the flow parameters. Adjoint-based optimization is performed in order to determine initial conditions that lead to maximum growth of disturbances in finite time. Results indicate that the rate of growth of the perturbation energy at small wavenumbers (less affected by viscosity initially) is dominated by two-dimensional modes (no spanwise variation). Substantial transient growth is observed at higher wave modes initially, followed by asymptotic decay of the perturbations at large time. Sensitivity of perturbation growth with respect to initial time, density and viscosity ratios is investigated. This work is complementary to previous inviscid analysis of this configuration, which showed that the interface was unconditionally unstable at all wave modes, even in the presence of surface tension, and that instability grew as the exponential of time squared.

Mechanics↗

Absolute vs Convective Instabilities and Front Propagation in Lipid Membrane Tubes

Here we analyze the stability of biological membrane tubes, with and without a base flow of lipids. Membrane dynamics are completely specified by two dimensionless numbers: the well-known Föppl-von Kármán number Γ and the recently introduced Scriven-Love number SL, respectively quantifying the base tension and base flow speed. For unstable tubes, the growth rate of a local perturbation depends only on Γ, whereas SL governs the absolute versus convective nature of the instability. Furthermore, nonlinear simulations of unstable tubes reveal an initially localized disturbance result in propagating fronts, which leave a thin atrophied tube in their wake. Depending on the value of Γ, the thin tube is connected to the unperturbed regions via oscillatory or monotonic shape transitions-reminiscent of recent experimental observations on the retraction and atrophy of axons. We elucidate our findings through a weakly nonlinear analysis, which shows membrane dynamics may be approximated by a model of the class of extended Fisher-Kolmogorov equations. Our study sheds light on the pattern selection mechanism in axonal shapes by recognizing the existence of two Lifshitz points, at which the front dynamics undergo steady-to-oscillatory bifurcations.

59 BASIC BIOLOGICAL SCIENCES↗

Time-dependent Computational Studies of Premixed Flames in Microgravity

This report describes the research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of NASA Microgravity Science and Applications Program. The primary focus of this research is on investigating fundamental questions concerning the propagation and extinction of premixed flames in earth gravity and in microgravity environments. Our approach is to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. Recently more advanced wall boundary conditions such as isothermal and no-slip have been added to the model. This enables the simulation of flames propagating in more practical systems than before. We have used the numerical simulations to investigate the effects of heat losses and buoyancy forces on the structure and stability of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the interaction between the various processes leading to flame instabilities and extinguishment, and to study the dynamics of cell formation and splitting. Our studies have been able to bring out the differences between upward- and downward-propagating flames and predict the zero-gravity behavior of these flames. The simulations have also highlighted the dominant role of wall heat losses in the case of downward-propagating flames. The simulations have been able to qualitatively predict the formation of multiple cells and the cessation of cell-splitting. Our studies have also shown that some flames in a microgravity environment can be extinguished due to a chemical instability and without any external losses. However, further simulations are needed to more completely understand upward-propagating and zero-gravity flames as well as to understand the potential effect of radiative heat losses.

Kailasanath, K.↗

Techniques for studying materials under extreme states of high energy density compression

The properties of materials under extreme conditions of pressure and density are of key interest to a number of fields, including planetary geophysics, materials science, and inertial confinement fusion. In geophysics, the equations of state of planetary materials, such as hydrogen and iron, under ultrahigh pressure and density provide a better understanding of their formation and interior structure [Celliers et al., “Insulator-metal transition in dense fluid deuterium,” Science 361, 677–682 (2018) and Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)]. The processes of interest in these fields occur under conditions of high pressure (100 GPa–100 TPa), high temperature (>3000 K), and sometimes at high strain rates (>103 s−1) depending on the process. With the advent of high energy density (HED) facilities, such as the National Ignition Facility (NIF), Linear Coherent Light Source, Omega Laser Facility, and Z, these conditions are reachable and numerous experimental platforms have been developed. To measure compression under ultrahigh pressure, stepped targets are ramp-compressed and the sound velocity, measured by the velocity interferometer system for any reflector diagnostic technique, from which the stress-density of relevant materials is deduced at pulsed power [M. D. Knudson and M. P. Desjarlais, “High-precision shock wave measurements of deuterium: Evaluation of exchange-correlation functionals at the molecular-to-atomic transition,” Phys. Rev. Lett. 118, 035501 (2017)] and laser [Smith et al., “Equation of state of iron under core conditions of large rocky exoplanets,” Nat. Astron. 2, 591–682 (2018)] facilities. To measure strength under high pressure and strain rates, experimenters measure the growth of Rayleigh–Taylor instabilities using face-on radiography [Park et al., “Grain-size-independent plastic flow at ultrahigh pressures and strain rates,” Phys. Rev. Lett. 114, 065502 (2015)]. The crystal structure of materials under high compression is measured by dynamic x-ray diffraction [Rygg et al., “X-ray diffraction at the national ignition facility,” Rev. Sci. Instrum. 91, 043902 (2020) and McBride et al., “Phase transition lowering in dynamically compressed silicon,” Nat. Phys. 15, 89–94 (2019)]. Medium range material temperatures (a few thousand degrees) can be measured by extended x-ray absorption fine structure techniques, Yaakobi et al., “Extended x-ray absorption fine structure measurements of laser-shocked V and Ti and crystal phase transformation in Ti,” Phys. Rev. Lett. 92, 095504 (2004) and Ping et al., “Solid iron compressed up to 560 GPa,” Phys. Rev. Lett. 111, 065501 (2013), whereas more extreme temperatures are measured using x-ray Thomson scattering or pyrometry. This manuscript will review the scientific motivations, experimental techniques, and the regimes that can be probed for the study of materials under extreme HED conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron Temperature Gradient Driven Transport Model for Tokamak Plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the Multi-Mode anomalous transport module [T. Rafiq \textit{et al.,} Phys Plasmas \textbf{20}, 032506 (2013)] to predict a time dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, $\beta$ and gradient of $\beta$ $(\betap)$, collisionality, safety factor, and toroidicity, where $\beta$ is the ratio of plasma pressure to the magnetic pressure. The electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

ETG↗

Instability and treatments of the coupled discrete element and lattice Boltzmann method by the immersed moving boundary scheme

The immersed moving boundary (IMB) scheme has been extensively used to couple the discrete element method (DEM) with the lattice Boltzmann method (LBM). In the literature, only the formulation of IMB for lattice nodal cells covered by a single-solid particle was given. The treatment of situations where a nodal cell is covered by two or more solid particles is seldom discussed. It is found that some numerical instability can occur for such situations due to an inappropriate computation of the weighting function in the IMB formulation. This work presents an enhanced treatment that can resolve the issue and validates it using some benchmark tests. Furthermore, to avoid the extra costs associated with the treatment and simplify the complicated procedure introduced, a simplified IMB scheme is proposed. The accuracy of both enhanced and simplified IMB schemes are validated by test cases including single-particle sedimentation, two-particle drafting-kissing-tumbling phenomenon, and multiple-particle sedimentation. Then, the robustness of both schemes is examined and discussed using a specially designed flow past cylinders test. Overall, the simplified IMB scheme is proved to be robust and sufficiently accurate and simpler and more effective than the enhanced scheme.

42 ENGINEERING↗

Isotope effects on energy transport in the core of ASDEX-Upgrade tokamak plasmas: Turbulence measurements and model validation

Design and operation of future tokamak fusion reactors using a deuterium–tritium 50:50 mix requires a solid understanding of how energy confinement properties change with ion mass. This study looks at how turbulence and energy transport change in L-mode plasmas in the ASDEX Upgrade tokamak when changing ion species between hydrogen and deuterium. For this purpose, both experimental turbulence measurements and modeling are employed. Local measurements of ion-scale (with wavevector of fluctuations perpendicular to the B-field k⊥< 2 cm−1, k⊥ρs< 0.2, where ρs is the ion sound Larmor radius using the deuterium ion mass) electron temperature fluctuations have been performed in the outer core (normalized toroidal flux ρTor=0.65−0.8) using a multi-channel correlation electron cyclotron emission diagnostic. Lower root mean square perpendicular fluctuation amplitudes and radial correlation lengths have been measured in hydrogen vs deuterium. Measurements of the cross-phase angle between a normal-incidence reflectometer and an ECE signal were made to infer the cross-phase angle between density and temperature fluctuations. The magnitude of the cross-phase angle was found larger (more out-of-phase) in hydrogen than in deuterium. TRANSP power balance simulations show a larger ion heat flux in hydrogen where the electron-ion heat exchange term is found to play an important role. These experimental observations were used as the basis of a validation study of both quasilinear gyrofluid trapped gyro-Landau fluid-SAT2 and nonlinear gyrokinetic GENE codes. Linear solvers indicate that, at long wavelengths (k⊥ρs<1), energy transport in the deuterium discharge is dominated by a mixed ion-temperature-gradient (ITG) and trapped-electron mode turbulence while in hydrogen transport is exclusively and more strongly driven by ITG turbulence. The Ricci validation metric has been used to quantify the agreement between experiments and simulations taking into account both experimental and simulation uncertainties as well as four different observables across different levels of the primacy hierarchy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Collisionless ion-electron energy exchange in magnetized shocks

Energy partition between ions and electrons in collisionless shocks has been a long-standing unsolved fundamental physical question. Here, we show that kinetic simulations of moderate Alfv´enic Mach number, magnetized, collisionless shocks reveal rapid, faster-than-Coulomb, energy exchange between ions and electrons when the plasma is sufficiently magnetized. Using kinetic and multi-fluid models with counter-streaming ions, we identify resonances between electron whistler and ion magnetohydrodynamic waves that account for this rapid energy exchange.

High-energy-density plasmas↗

Fully Kinetic Shearing-box Simulations of Magnetorotational Turbulence in 2D and 3D. I. Pair Plasmas

The magnetorotational instability (MRI) is a fundamental mechanism determining the macroscopic dynamics of astrophysical accretion disks. In collisionless accretion flows around supermassive black holes, MRI-driven plasma turbulence cascading to microscopic (i.e., kinetic) scales can result in enhanced angular-momentum transport and redistribution, nonthermal particle acceleration, and a two-temperature state where electrons and ions are heated unequally. However, this microscopic physics cannot be captured with standard magnetohydrodynamic (MHD) approaches typically employed to study the MRI. In this work, we explore the nonlinear development of MRI turbulence in a pair plasma, employing fully kinetic particle-in-cell (PIC) simulations in two and three dimensions. First, we thoroughly study the axisymmetric MRI with 2D simulations, explaining how and why the 2D geometry produces results that differ substantially from 3D MHD expectations. We then perform the largest (to date) 3D simulations, for which we employ a novel shearing-box approach, demonstrating that 3D PIC models can reproduce the mesoscale (i.e., MHD) MRI dynamics in sufficiently large runs. With our fully kinetic simulations, we are able to describe the nonthermal particle acceleration and angular-momentum transport driven by the collisionless MRI. Since these microscopic processes ultimately lead to the emission of potentially measurable radiation in accreting plasmas, our work is of prime importance to understand current and future observations from first principles, beyond the limitations imposed by fluid (MHD) models. While in this first study we focus on pair plasmas for simplicity, our results represent an essential step toward designing more realistic electron–ion simulations, on which we will focus in future work.

79 ASTRONOMY AND ASTROPHYSICS↗

CFD Integrated Transition Modeling for High-Speed Flows via Coupled OVERFLOW-LASTRAC Analysis

This work details ongoing efforts at the NASA Langley Research Center to develop and validate a general-use CFD tool that includes built-in predictions of boundary-layer transition in high-speed flows. Existing tools for the coupling of the NASA OVERFLOW structured overset RANS solver and LASTRAC stability analysis code have been extended to capture boundary-layer transition in high-speed flows driven by either Mack’s first-mode or secondmode instabilities and/or crossflow instabilities. The efficacy of this coupled approach has been demonstrated by examining a variety of supersonic and hypersonic test cases for which experimental validation data is available. The intermittency prescription parameters inherited from prior low-speed applications are found to be suitable for some, but not all high-speed flow scenarios. The method is found to be robust for straight cone configurations and improvements necessary to accurately capture transition on geometries with spatially complex amplification factor envelopes with intermediate regions of slow N-factor variation are examined. The automated, coupled analysis is demonstrated for a 3d supersonic test case and found to perform well within the limits of the linear stability analysis on which it relies. Finally, a preliminary investigation of the method’s robustness to non-ideal CFD meshes is conducted.

High-speed flow↗

A Transient Hydrodynamic Model of Screen Channel Liquid Acquisition Devices for In-Space Cryogenic Propellant Management

Screen channel liquid acquisition devices (LADs) will play a crucial role in future deep space travel. It is essential that vapor-free delivery of propellants during tank-to-tank transfer is ensured to maximize yield from storage tanks and prevent potential combustion instabilities. The screen channel LAD utilizes a fine screen wire mesh that can separate phases in a low Bond number (i.e. microgravity) environment using surface tension forces. This study presents the development and verification of a new model for transient screen compliance, one of the influential factors for screen channel LAD design. Screen compliance is crucial during LAD channel outflow transients because the slight deflection of the screen can provide needed mass to satisfy rapid outflow demands and reduce the pressure difference across the screen. The model is successfully verified against CFD simulations. In addition, the characteristic speed for the governing screen compliance equations is derived which allows for numerical stability criteria to be established. As shown in this study, the transient maximum pressure difference across the screen can greatly exceed the steady state maximum pressure difference across the screen in many cases.

Hydrodynamics Simulations↗

A learned closure method applied to phase mixing in a turbulent gradient-driven gyrokinetic system in simple geometry

We present a new method for formulating closures that learn from kinetic simulation data. We apply this method to phase mixing in a simple gyrokinetic turbulent system – temperature-gradient-driven turbulence in an unsheared slab. The closure, called the learned multi-mode (LMM) closure, is constructed by, first, extracting an optimal basis from a nonlinear kinetic simulation using singular value decomposition. Subsequent nonlinear fluid simulations are projected onto this basis and the results are used to formulate the closure. We compare the closure with other closures schemes over a broad range of the relevant two-dimensional parameter space (collisionality and gradient drive). We find that the turbulent kinetic system produces phase-mixing rates much lower than the linear expectations, which the LMM closure is capable of capturing. We also compare radial heat fluxes. A Hammett–Perkins closure, generalized to include collisional effects, is quite successful throughout the parameter space, producing ${\sim }14\,\%$ root-mean-square (r.m.s.) error. The LMM closure is also very effective: when trained at three (two) points (in a 35 point parameter grid), the LMM closure produces $8\,\%$ ( $12\,\%$ ) r.m.s. errors. The LMM procedure can be readily generalized to other closure problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗