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 991 records · Page 55

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Visualization of R410A Flow Boiling Inside Horizontal Smooth Tubes under Diabatic Conditions

This article presents a novel approach to visualize the flow boiling in round tubes under “diabatic conditions”. In the visualization test section, a 3D printed round tube is placed inside a glass tube with secondary fluid at a higher temperature flowing in the annulus between the two tubes thus providing diabatic conditions. Experimental work was conducted for evaporating R410A flow in a horizontal smooth round tube with the inner diameter of 6.3 mm. Compared to the traditional adiabatic visualization technique, the more realistic evaporative condition is simulated and the developing evaporation process is captured. New insights for the flow behavior at low vapor quality are brought by this diabatic method. The vapor plug occurs at very low vapor quality, and this is attributed to the coalescence of bubbles generated on the heated tube wall. For the intermittent flow, liquid slug triggered by Kelvin–Helmholtz instability affects the nucleation on the tube wall. It brings the fresh liquid to rewet the dry surface and activate the nucleation sites. Finally, heat flux condition affects the heat transfer mechanisms during the flow boiling, and the heat transfer through the nucleate boiling mechanism increases as the heat flux increases.

3D printed↗

NASA NTRS2003

Mixing and Transport in the Solar Nebula

Boss & Vanhala (2000, 2001) prepared reviews of triggered collapse and injection models, using Prudence Foster's finite differences code at very high spatial resolution (440 x 1440 cells) to demonstrate the convergence of the R-T fingers in triggered injection models. A two dimensional hydrodynamical calculation with unprecedentedly high spatial resolution (960 x 2880 zones, or almost 3 million grid points) demonstrated that it suitable shock front can both trigger the collapse of an otherwise stable presolar cloud, and inject shock front particles into the collapsing cloud through the formation of what become Rayleigh-Taylor fingers of compressed fluid layers falling into the gravitational potential well of the growing protostar. These calculations suggest that heterogeneity derived from these R-T fingers will persist down to the scale of their injection onto the surface of the solar nebula. Haghighipour developed a numerical code capable of calculating the orbital evolution of dust grains of varied sizes in a gaseous nebula, subject to Epstein and Stokes drag as well as the self-gravity of the disk. In collaboration with the PI and George W. Wetherill, Haghighipour has been involved in development of a new idea on the possibility of rapid formation of ice giant planets via the disk instability mechanism. Haghighipour studied the stability of a five-body system consisting of the Sun and four protoplanets by numerically integrating their equations of motions. Using Levison and Duncan s SWIFT integrator, Haghighipour showed that, depending on the orbital parameters of the bodies, such a system can be stable for 0.1-10 Myr. Time periods of 1 Myr or more are long enough to be consistent with the time scale proposed for the formation of giant planets by the disk instability mechanism and the photoevaporation of the gaseous envelopes of the outermost protoplanets by a nearby OB star, resulting in the formation of ice giant planets. The PI has used his three dimensional models of marginally gravitationally unstable disks to study the preservation of isotopic heterogeneity in evolving protoplanetary disks. Such heterogeneity might arise from the infall onto the disk s surface of solids processed in the X-wind region of the disk, or derived from stellar nucleosynthesis and injected by R-T fingers. The technique used consists of solving a color equation, identical to the gas continuity equation, which follows the time evolution in three space dimensions of an arbitrarily placed initial color field, i.e., a dye inserted the disk. The models show that significant concentrations of color could persist for time periods of about a thousand years or more, even in the most dynamically active region of such a disk. Such a time period might be long enough for solids to coagulate and grow to significant sizes while retaining the isotopic signature of their birth region in the nebula.

Boss, Alan P.↗

NASA NTRS2002

Damping of Quasi-stationary Waves Between Two Miscible Liquids

Two viscous miscible liquids with an initially sharp interface oriented vertically inside a cavity become unstable against oscillatory external forcing due to Kelvin-Helmholtz instability. The instability causes growth of quasi-stationary (q-s) waves at the interface between the two liquids. We examine computationally the dynamics of a four-mode q-s wave, for a fixed energy input, when one of the components of the external forcing is suddenly ceased. The external forcing consists of a steady and oscillatory component as realizable in a microgravity environment. Results show that when there is a jump discontinuity in the oscillatory excitation that produced the four-mode q-s wave, the interface does not return to its equilibrium position, the structure of the q-s wave remains imbedded between the two fluids over a long time scale. The damping characteristics of the q-s wave from the time history of the velocity field show overdamped and critically damped response; there is no underdamped oscillation as the flow field approaches steady state. Viscous effects serve as a dissipative mechanism to effectively damp the system. The stability of the four-mode q-s wave is dependent on both a geometric length scale as well as the level of background steady acceleration.

Duval, Walter M. B.↗

DOE OSTI2021

Stability-based optimization of ACRT for the growth of CZT by the traveling heater method

Numerical modeling is employed to study the effects of several accelerated crucible rotation technique (ACRT) rotation schedules on the growth of cadmium zinc telluride (CZT) by the traveling heater method (THM). In conjunction with these analyses, a measure is developed to estimate the amplitude of disturbances associated with cellular interface growth morphologies that develop over time and across the melt-solid interface. Furthermore, this disturbance amplitude measure is used to assess the ability of ACRT to stabilize the growth interface and thereby minimize inclusion formation processes. Notably, our analyses reveal that interfacial instability is not well correlated with melt undercooling, explaining why classical ACRT approaches to reduce undercooling by mixing may not stabilize growth. Computations show that interface stability is best achieved via rotation schedules that accentuate outward Ekman flows during spin-up and, during spin-down, eliminate Taylor-Görtler flows and minimize inward Ekman flows.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Moiré nematic phase in twisted double bilayer graphene

Graphene moiré superlattices display electronic flat bands. At integer fillings of these flat bands, energy gaps due to strong electron-electron interactions are generally observed. However, the presence of other correlation-driven phases in twisted graphitic systems at non-integer fillings is unclear. Here, we report the existence of threefold rotational symmetry breaking in twisted double bilayer graphene. Using spectroscopic imaging over large and uniform areas to characterize the direction and degree of C 3 symmetry breaking, we find it to be prominent only at energies corresponding to the flat bands and nearly absent in the remote bands. We demonstrate that the magnitude of the rotational symmetry breaking does not depend on the degree of the heterostrain or the displacement field, being instead a manifestation of an interaction-driven electronic nematic phase. We show that the nematic phase is a primary order that arises from the normal metal state over a wide range of doping away from charge. Furthermore, our modelling suggests that the nematic instability is not associated with the local scale of the graphene lattice, but is an emergent phenomenon at the scale of the moiré lattice.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Exploring implosion designs for increased compression on the National Ignition Facility using high density carbon ablators

It has long been recognized that high compression, and hence good confinement, is essential to achieving high yields in inertial confinement fusion implosions. In pursuit of multi-megajoule yields on the National Ignition Facility (NIF), a new campaign has begun aimed at testing the hypothesis that controlling hydrodynamic stability is key to achieving effective higher compression with the high density carbon ablators currently fielded on NIF. This campaign is built around a new implosion design, called SQ-n, that is derived from the uniquely stable Bigfoot design tested on NIF in 2016–2019. While very stable and with performance that was quite close to one-dimensional expectations, Bigfoot was a relatively high adiabat, and consequently lower compression design. The goal of SQ-n is then to evolve Bigfoot toward a higher compression design but without compromising its unique stability characteristics. Specifically, SQ-n adopts a ramped foot pulse shape to minimize early time Richtmyer–Meshkov instability growth and uses an ablator dopant distribution extending all of the way to the fuel–ablator interface that simulations suggest further reduces perturbation growth. In conclusion, this paper describes the design philosophy pursued with SQ-n, the results of instability modeling of the candidate design, and the experimental campaign planned to test these ideas in the near future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1982

Vortex simulation of an inviscid shear layer

The accuracy of the vortex-blob method was tested by simulating a free-shear-layer instability, Kirchhoff's elliptical vortex, and a circular vortex. The main numerical parameters in the vortex-blob method are the density of the vortices, and the distribution of vorticity within each vortex core. The growth rate of a periodic unstable mode of the shear layer was calculated numerically and compared with the exact result. The error is only a few percent for about 10 rows of vortex blobs. The error is reduced by decreasing the spacing between vortices and, correspondingly, the core size. In the simulation of the motion of the elliptical vortex, the rotation of the boundary, without change of shape, and the circular particle paths of the vortical fluid were well simulated. For the circular vortex, optimum sets of parameters were obtained by comparing them with the exact velocity. The results are consistent with convergence theories of the vortex-blob method. In particular, second-order convergence is observed with a Gaussian core from velocity calculation.

Nakamura, Y.↗

NASA NTRS1992

Fluid flow in a jet and the Ca XIX line profiles observed during solar flares

The present study considers two types of fluid models as methods for reproducing the blueshifts and line broadening observed in soft X-ray lines during the rise of solar flares. The fluid models representing the laminar flow of material in a jet and through a pipe were used to derive the velocity at each location in the flow. The optimal values for the free parameters in the jet model were within the ranges expected for solar flare conditions. The jet was much smaller than typical flare volumes estimated from spatial images, but the resulting densities were consistent with the upper limit of densities derived using density-sensitive line ratios. This supports previous reports of small filling factors and is interpreted as evidence for a chromospheric origin for the flow. Stability analysis indicates that, although the laminar jet can reproduce the observations, instabilities may cause the flow to become turbulent, and that this laminar treatment is not strictly valid.

Bornmann, P. L.↗

NASA NTRSDate not supplied

Full Stall Simulations of a Redesigned Ventilation Fan for the ISS

The concept of a stall is studied rigorously in the aerospace industry. From a design standpoint, instabilities such as stall are undesirable– operation in the stall regime has a tremendous impact on aerodynamic performance as well as structural integrity. In extreme cases, operating in stall conditions can cause failure. Further, stall can cause a loss of lift on aircraft wings or a loss of thrust in aircraft engines. In any case, stall continues to be a topic of interest in the aerospace industry. The present work aims to analyze the stall characteristics of a ventilation fan that was recently designed for the International Space Station (ISS). Although the ventilation fan has a rotor-stator design, this paper considers a rotor-only configuration. The FUN3D Computational Fluid Dynamic (CFD) solver developed by NASA Langley Research Center was used to simulate the operational characteristics of the ventilation fan. FUN3D solves the Unsteady Reynolds-Averaged Naiver-Stokes (URANS) equations using implicit time marching and a dynamic overset grid. The FUN3D solver was originally written for exterior flow fields; however, this work represents an extension of the FUN3D solver to turbomachinery or interior flow fields. The FUN3D results for the rotor-only ventilation fan accurately captured the operational characteristics inherent to compressors– the results shared similar performance trends when compared to the experimental results for the rotor-stator case. A peak adiabatic efficiency of 96% occurred at a MFR of 105.2 CFM, the minimum aerodynamically stable point. A computationally stable stall occurred at a mass flow rate of 43.8 CFM where the adiabatic efficiency dropped to 69%.

CFD↗

DOE OSTI2021

Simulation and flow physics of a shocked and reshocked high-energy-density mixing layer

This work describes a computational investigation of multimode instability growth and multimaterial mixing induced by multiple shock waves in a high-energy-density (HED) environment, where pressures exceed 1 Mbar. The simulations are based on a series of experiments performed at the National Ignition Facility (NIF) and designed as an HED analogue of non-HED shock-tube studies of the Richtmyer–Meshkov instability and turbulent mixing. A three-dimensional computational modelling framework is presented. It treats many complications absent from canonical non-HED shock-tube flows, including distinct ion and free-electron internal energies, non-ideal equations of state, radiation transport and plasma-state mass diffusivities, viscosities and thermal conductivities. The simulations are tuned to the available NIF data, and traditional statistical quantities of turbulence are analysed. Integrated measures of turbulent kinetic energy and enstrophy both increase by over an order of magnitude due to reshock. Large contributions to enstrophy production during reshock are seen from both the baroclinic source and enstrophy–dilatation terms, highlighting the significance of fluid compressibility in the HED regime. Dimensional analysis reveals that Reynolds numbers and diffusive Péclet numbers in the HED flow are similar to those in a canonical non-HED analogue, but conductive Péclet numbers are much smaller in the HED flow due to efficient thermal conduction by free electrons. It is shown that the mechanism of electron thermal conduction significantly softens local spanwise gradients of both temperature and density, which causes a minor but non-negligible decrease in enstrophy production and small-scale mixing relative to a flow without this mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2020

Flow Simulations for Martian and Lunar Lander Plume-Surface Interaction Prediction

Landing vehicle propulsive systems introduce several risks associated with plume/surface interactions (PSI) during the final stages of descent. For example, view obscuration during descent due to the liberation of dust driven by the rocket plume’s interaction with the lunar surface was reported on several Apollo missions. Plume-driven debris impacts on vehicle and nearby assets as well as landing instabilities due to crater formation are also potential risks associated with plume/surface interactions. The importance of studying these plume/surface interactions has been elevated due to NASA’s mission to land humans on the Moon by 2024. Existing and under-development predictive simulation tools are currently being brought to application readiness under a project awarded by the NASA Space Technology Mission Directorate (STMD) Game Changing Development (GCD) program. This tool development is divided into four tasks: (1) simulating plume structure in low-pressure environments, (2) simulating crater development and ejecta, (3) regolith particle phase modeling, and (4) gas-particle interaction modeling. Recent efforts under task 1 include demonstrating existing computational fluid dynamics (CFD) tools for simulating plume structure in low pressure environments similar to those found on Mars. The results of these demonstrative simulations are compared to available experimental data. Emerging capability being developed under task 2 is also presented, demonstrating simulations of crater formation in Lunar conditions.

Thomas P Shurtz↗

DOE OSTI2023

A Conservative Finite-Volume Based Interface-Tracking Algorithm Using the Signed Distance Function

Methods for tracking an interface between two fluid phases are developed to ensure desired fluid properties, conservation, and stability are preserved in a finitevolume (FV) discretization. Here, the interface is tracked using a level set method where the signed distance function implicitly defines the interface. Marching methods are used to evaluate the value of the signed distance function, including a novel initialization method to initialize any implicit function to the signed distance function around sharp corners in the level set. Global conservation and consistency with a set of governing equations is enforced by a compression coefficient that measures the volumetric compression or expansion due to inaccuracies in the level set evaluation. A redistribution method is integrated into the volume correction to eliminate the small-cell instability while maintaining global conservation. This suite of methods is implemented and tested using static uniform velocity, and potential flow cases with multiple interface geometries. Results show these methods achieve up to second order accuracy, and are conservative. The application for these methods is intended to track the interface of a 3D printing filament in a finite-volume discretization of the all-speed Navier-Stokes equations.

42 ENGINEERING↗

DOE OSTI2022

Ephemeral antibubbles: Spatiotemporal evolution from direct numerical simulations

Antibubbles, which consist of a shell of a low-density fluid inside a high-density fluid, have several promising applications. We show, via extensive direct numerical simulations (DNSs), in both two and three dimensions, that the spatiotemporal evolution of antibubbles can be described naturally by the coupled Cahn-Hilliard-Navier-Stokes (CHNS) equations for a binary fluid. Our DNSs capture elegantly the gravity-induced thinning and breakup of an antibubble via the time evolution of the Cahn-Hilliard scalar-order-parameter field $\phi$, which varies continuously across interfaces, so we do not have to enforce complicated boundary conditions at the moving antibubble interfaces. To ensure that our results are robust, we supplement our CHNS simulations with sharp-interface volume-of-fluid DNSs. We track the thickness of the antibubble and calculate the dependence of the lifetime of an antibubble on several parameters; we show that our DNS results agree with various experimental results; in particular, the velocity with which the arms of the antibubble retract after breakup scales as σ 1/2 , where σ is the surface tension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2015

Advancing the State-of-the-Practice for Liquid Rocket Engine Injector Design

Current shortcomings in both the overall injector design process and its underlying combustion stability assessment methodology are rooted in the use of empirically based or low fidelity representations of complex physical phenomena and geometry details that have first order effects on performance, thermal environments and combustion stability. The result is a design and analysis capability that is often inadequate to reliably arrive at a suitable injector design in an efficient manner. Specifically, combustion instability has been particularly difficult to predict and mitigate. Large hydrocarbon-fueled booster engines have been especially problematic in this regard. Where combustion instability has been a problem, costly and time-consuming redesign efforts have often been an unfortunate consequence. This paper presents an overview of a recently completed effort at NASA Marshall Space Flight Center to advance the state-of-the-practice for liquid rocket engine injector design. Multiple perturbations of a gas-centered swirl coaxial (GCSC) element that burned gaseous oxygen and RP-1 were designed, assessed for combustion stability, and tested. Three designs, one stable, one marginally unstable and one unstable, were used to demonstrate both an enhanced overall injector design process and an improved combustion stability assessment process. High-fidelity results from state-of-the-art computational fluid dynamics CFD simulations were used to substantially augment and improve the injector design methodology. The CFD results were used to inform and guide the overall injector design process. They were also used to upgrade selected empirical or low-dimensional quantities in the ROCket Combustor Interactive Design (ROCCID) stability assessment tool. Hot fire single element injector testing was used to verify both the overall injector designs and the stability assessments. Testing was conducted at the Air Force Research Laboratory and at Purdue University. Companion papers provide details of the overall injector design process, full- and sub-scale testing, ROCCID-based stability assessments and the CFD simulations.

Tucker, P. K.↗

NASA NTRS1982

Thresholds for the onset of fluid and magnetofluid turbulence

Linear stability calculations conducted in plasma physics are based on the theory of hydrodynamic stability of neutral fluids. The present investigation is concerned with the validity of procedures based on linear stability analysis in six much-studied situations. All have the property that the fluid goes from laminar to turbulent at critical values of some dimensionless number, such as the Reynolds number or the Rayleigh number. The question is, at what threshold values of the dimensionless number do the unstable motions set in, and can these thresholds be predicted by a linear analysis of the stability of the laminar state as the threshold is approached from the stable side. In three cases linear stability analysis clearly seems to fail. These cases include Plane Poiseuille Flow, Plane Couette Flow, and Cylindrical Pipe Flow (Hagen-Poiseuille Flow). Situations in which predictions provided by linear stability analysis are correct are related to Rotating Couette Flow, Thermally-Driven Convection, and Instability of Laminar Boundary Layers.

Montgomery, D.↗

NASA NTRS1985

A self-consistent linear-mode model of stellar convection

A normal-mode expansion of the linearized fluid equations in terms of small subset of spherical harmonics can provide a foundation for a physically motivated, self-consistent description of a solar-type convection zone. In the absence of dissipation, a second-order differential equation governs the radial dependence of the modes, so that interpretation of the effects on convection quantities of the normal-form 'potential well' is straightforward. The philosophy is quite different from the more recent work of Narasimha and Antia (1982): all envelopes presented here differ substantially from MLT envelopes, and therefore, from theirs, which are constructed to be consistent with MLT. The amplitude of all modes is set by a Kelvin-Helmholtz-('shear'-) instability argument unrelated to solar observations, with the result that the convection description may be considered to arise from 'first-hueristic-principles'. The thermodynamics modelled vaguely resemble the sun's, and more vigorously convective envelopes show some phenomena qualitatively like solar observations (e.g., atmospheric velocity spectra).

Macauslan, J.↗

NASA NTRS1989

A normal-mode approach to Jovian atmospheric dynamic

A nonlinear, quasi-geostrophic, baroclinic model of Jovian atmospheric dynamics is proposed, in which vertical variations of velocity are represented by a truncated sum over a complete set of orthogonal functions obtained by a separation of variables of the linearized quasi-geostrophic potential vorticity equation. A set of equations for the time variation of the mode amplitudes in the nonlinear case is then derived. It is shown that, for a planet with a neutrally stable, fluid interior instead of a solid lower boundary, the barotropic mode represents motions in the interior, and is not affected by the baroclinic modes. One consequence of this is that a normal-mode model with one baroclinic mode is dynamically equivalent to a one-layer model with solid lower topography. It is also shown that, for motions in Jupiter's cloudy lower troposphere, the stratosphere behaves nearly as a rigid lid, so that the normal-mode is applicable to Jupiter. The accuracy of the normal-mode model for Jupiter is tested using the following simple problems: (1) forced, vertically propagating Rossby waves, using two and three baroclinic modes, and (2) baroclinic instability, using two baroclinic modes. It is found that the normal-mode model provides qualitatively correct results, even with only a very limited number of vertical degrees of freedom.

Achterberg, Richard K.↗

DOE OSTI2021

Implosion dynamics of triple-nozzle gas-puff z pinches on COBRA

Experiments on the 1-MA, 220-ns COBRA generator at Cornell University were conducted to provide detailed measurements of structured cylindrical gas-puff z pinches. In the experiments, a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet initialize the z-pinch load with various working gases, radial density profiles, and externally applied axial magnetic fields. Planar laser-induced fluorescence provides a measure of the initial neutral gas density of the load, while three-frame laser shearing interferometry and multi-frame extreme ultraviolet (XUV) cameras reveal the formation and propagation of a magneto-Rayleigh–Taylor (MRT) unstable shock layer. Implosion trajectories are compared to simple, experimentally informed models and found to be in good agreement. Differences in the structure of the accelerating plasma sheath and evolution of the MRT instability are observed for different gas species and axial magnetic field strengths, correlating with differences in pinch uniformity and x-ray emission. Here, the average instability growth is compared to linear MRT theory predictions using the instantaneous acceleration of the best-fit implosion models and characteristic instability wavelength, with the effective Atwood number and seed perturbation size as fit parameters. For high density argon center jets, ionization prior to the arrival of the imploding plasma sheath suggests a heating mechanism consistent with photoionization by XUV self-emission.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗