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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Self-organization in collisionless, high- β turbulence

The magnetohydrodynamic (MHD) equations, as a collisional fluid model that remains in local thermodynamic equilibrium (LTE), have long been used to describe turbulence in myriad space and astrophysical plasmas. Yet, the vast majority of these plasmas, from the solar wind to the intracluster medium (ICM) of galaxy clusters, are only weakly collisional at best, meaning that significant deviations from LTE are not only possible but common. Recent studies have demonstrated that the kinetic physics inherent to this weakly collisional regime can fundamentally transform the evolution of such plasmas across a wide range of scales. Here, we explore the consequences of pressure anisotropy and Larmor-scale instabilities for collisionless, β $\gg$ 1, turbulence, focusing on the role of a self-organizational effect known as ‘magneto-immutability’. We describe this self-organization analytically through a high-β, reduced ordering of the Chew–Goldberger–Low-MHD (CGL-MHD) equations, finding that it is a robust inertial-range effect that dynamically suppresses magnetic-field-strength fluctuations, anisotropic-pressure stresses and dissipation due to heat fluxes. As a result, the turbulent cascade of Alfvénic fluctuations continues below the putative viscous scale to form a robust, nearly conservative, MHD-like inertial range. These findings are confirmed numerically via Landau-fluid CGL-MHD turbulence simulations that employ a collisional closure to mimic the effects of microinstabilities. We find that microinstabilities occupy a small (~5%) volume-filling fraction of the plasma, even when the pressure anisotropy is driven strongly towards its instability thresholds. We discuss these results in the context of recent predictions for ion-vs-electron heating in low-luminosity accretion flows and observations implying suppressed viscosity in ICM turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Fluid Solver for Studying Torsional Galloping in Solar-Tracking PV Panel Arrays

Solar-tracking photovoltaic arrays are susceptible to aeroelastic fluttering during high-wind events. This dynamic fluttering behavior can grow in amplitude until the panels enter an unstable mode known as torsional galloping which can lead to panel failure or total array destruction. To better understand the physics of the torsional galloping phenomenon, and to inform the discussion around panel design and recommended panel stow positions during high wind events, a fluid-structure interaction solver composed of a simulated atmospheric boundary layer with simplified panel structural responses was designed. The simulation choices and features of this solver were informed by the geometry and physical properties of an experimental panel array known to exhibit torsional galloping behavior during hind-wind events. These simulations revealed that the torsional galloping instability is driven by a combination of cyclic vortex shedding from the sun-facing side of the panel and the elastic properties of the torque tube linking the panel assemblies. Testing different stow angles across a range of wind speeds indicates that panels are generally more stable when stowed at negative angles where the leading edge is closer to the ground, hypothesized to be due to ground-blocking effects.

41 EE - Solar Energy Technologies Office (EE-4S)↗
NASA NTRS2000

A Technique for Rapidly Deploying a Concentration Gradient with Applications to Microgravity

The latter half of the last century has seen rapid advancements in semiconductor crystal growth powered by the demand for high performance electronics in myriad applications. The reduced gravity environment of space has also been used for crystal growth tests, especially in instances where terrestrial growth has largely been unsuccessful. While reduced gravity crystal growth affords some control of the gravity parameter, many crystals grown in space, to date, have structural flaws believed to result from convective motions during the growth phase. The character of these instabilities is not well understood but is associated with thermal and solutal density variations near the solidification interface in the presence of residual gravity and g-jitter. In order to study these instabilities in a separate, controlled space experiment, a concentration gradient would first have to be artificially established in a timely manner as an initial condition. This is generally difficult to accomplish in a microgravity environment because the momentum of the fluid injected into a test cell tends to swirl around and mix in the absence of a restoring force. The use of magnetic fields to control the motion and position of liquids has received growing interest in recent times. The possibility of using the force exerted by a non-uniform magnetic field on a ferrofluid to not only achieve fluid manipulation but also to actively control fluid motion makes it an attractive candidate for space applications. This paper describes a technique for quickly establishing a linear or exponential fluid concentration gradient using a magnetic field in place of gravity to stabilize the deployment. Also discussed is a photometric technique for measuring the concentration profile using light attenuation. Results of the ground-based experiments indicate that the concentration distribution is within 3% of the predicted value. Although any range of concentations can be realized, photometric constraints are discussed which impose some limitations on measurements.

Leslie, Fred↗
NASA NTRS2000

The Use of Ferrofluids to Model Materials Processing (MSFC Center Director's Discretionary Fund)

Many crystals grown in space have structural flaws believed to result from convective motions during the growth phase. The character of these instabilities is not well understood but is associated with thermal and solutal density variations near the solidification interface in the presence of residual gravity and g-jitter. To study these instabilities in a separate, controlled space experiment, a concentration gradient would first have to be artificially established in a timely manner as an initial condition. This is generally difficult to accomplish in a microgravity environment because the momentum of the fluid injected into a test cell tends to swirl around and mix in the absence of a restoring force. The use of magnetic fields to control the motion and position of liquids has received recent, growing interest. The possibility of using the force exerted by a non-uniform magnetic field on a ferrofluid to not only achieve fluid manipulation but also to actively control fluid motion makes it an attractive candidate for space applications. This paper describes a technique for quickly establishing a linear or exponential fluid concentration gradient using a magnetic field in place of gravity to stabilize the deployment. Also discussed is a photometric technique for measuring the concentration profile using light attenuation. Although any range of concentrations can be realized, photometric constraints impose some limitations on measurements. Results of the ground-based experiments indicate that the species distribution is within 3 percent of the predicted value.

Leslie, F.↗
NASA NTRS2003

Numerical Simulation of Protoplanetary Vortices

The fluid dynamics within a protoplanetary disk has been attracting the attention of many researchers for a few decades. Previous works include, to list only a few among many others, the well-known prescription of Shakura & Sunyaev, the convective and instability study of Stone & Balbus and Hawley et al., the Rossby wave approach of Lovelace et al., as well as a recent work by Klahr & Bodenheimer, which attempted to identify turbulent flow within the disk. The disk is commonly understood to be a thin gas disk rotating around a central star with differential rotation (the Keplerian velocity), and the central quest remains as how the flow behavior deviates (albeit by a small amount) from a strong balance established between gravitational and centrifugal forces, transfers mass and momentum inward, and eventually forms planetesimals and planets. In earlier works we have briefly described the possible physical processes involved in the disk; we have proposed the existence of long-lasting, coherent vortices as an efficient agent for mass and momentum transport. In particular, Barranco et al. provided a general mathematical framework that is suitable for the asymptotic regime of the disk; Barranco & Marcus (2000) addressed a proposed vortex-dust interaction mechanism which might lead to planetesimal formation; and Lin et al. (2002), as inspired by general geophysical vortex dynamics, proposed basic mechanisms by which vortices can transport mass and angular momentum. The current work follows up on our previous effort. We shall focus on the detailed numerical implementation of our problem. We have developed a parallel, pseudo-spectral code to simulate the full three-dimensional vortex dynamics in a stably-stratified, differentially rotating frame, which represents the environment of the disk. Our simulation is validated with full diagnostics and comparisons, and we present our results on a family of three-dimensional, coherent equilibrium vortices.

Lin, H.↗
DOE OSTI2024

A numerical method for simulating variable density flows in membrane desalination systems

Here, we present a novel method for simulating unsteady, variable density, fluid flows in membrane desalination systems. By assuming the density varies only with concentration and temperature, the scheme decouples the solution of the governing equations into two sequential blocks. The first solves the governing equations for the temperature and concentration fields, which are used to compute all thermophysical properties. The second block solves the conservation of mass and momentum equations for the velocity and pressure. We show that this is computationally more efficient than schemes that iterate over the full coupled equations in one block. We verify that the method achieves second-order spatial-temporal accuracy, and we use the method to investigate buoyancy-driven convection in a desalination process called vacuum membrane distillation. Specifically, we show that with gravity properly oriented, variations in temperature and concentration can trigger a double-diffusive instability that enhances mixing and improves water recovery. We also show that the instability can be strengthened by providing external heating.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗
NASA NTRS1993

Flow instability in particle-bed nuclear reactors

The particle-bed core offers mitigation of some of the problems of solid-core nuclear rocket reactors. Dividing the fuel elements into small spherical particles contained in a cylindrical bed through which the propellant flows radially, may reduce the thermal stress in the fuel elements, allowing higher propellant temperatures to be reached. The high temperature regions of the reactor are confined to the interior of cylindrical fuel assemblies, so most of the reactor can be relatively cool. This enables the use of structural and moderating materials which reduce the minimum critical size and mass of the reactor. One of the unresolved questions about this concept is whether the flow through the particle-bed will be well behaved, or will be subject to destructive flow instabilities. Most of the recent analyses of the stability of the particle-bed reactor have been extensions of the approach of Bussard and Delauer, where the bed is essentially treated as an array of parallel passages, so that the mass flow is continuous from inlet to outlet through any one passage. A more general three dimensional model of the bed is adopted, in which the fluid has mobility in three dimensions. Comparison of results of the earlier approach to the present one shows that the former does not accurately represent the stability at low Re. The more complete model presented should be capable of meeting this deficiency while accurately representing the effects of the cold and hot frits, and of heat conduction and radiation in the particle-bed. It can be extended to apply to the cylindrical geometry of particle-bed reactors without difficulty. From the exemplary calculations which were carried out, it can be concluded that a particle-bed without a cold frit would be subject to instability if operated at the high temperatures desired for nuclear rockets, and at power densities below about 4 megawatts per liter. Since the desired power density is about 40 megawatts per liter, it can be concluded that operation at design exit temperature but at reduced power could be hazardous for such a reactor. But the calculations also show that an appropriate cold frit could very likely cure the instability. More definite conclusions must await calculations for specific designs.

Kerrebrock, Jack L.↗
DOE OSTI2021

Alpha particle driven Alfvénic instabilities in ITER post-disruption plasmas

Fusion-born alpha particles in ITER disruption simulations are investigated as a possible drive of Alfvénic instabilities. The ability of these waves to expel runaway electron (RE) seed particles is explored in the pursuit of a passive, inherent RE mitigation scenario. The spatiotemporal evolution of the alpha particle distribution during the disruption is calculated using the linearized Fokker–Planck solver CODION coupled to a fluid disruption simulation. These simulations are done in the limit of no alpha particle transport during the thermal quench, which can be seen as a most pessimistic situation where there is also no RE seed transport. Under these assumptions, the radial anisotropy of the resulting alpha population provides free energy to drive Alfvénic modes during the quench phase of the disruption. We use the linear gyrokinetic magnetohydrodynamic code LIGKA to calculate the Alfvén spectrum and find that the equilibrium is capable of sustaining a wide range of modes. The self-consistent evolution of the mode amplitudes and the alpha distribution is calculated utilizing the wave-particle interaction tool HAGIS. Intermediate mode number (n = 7–15, 22–26) toroidal Alfvén eigenmodes are shown to saturate at an amplitude of up to δB/B ≈ 0.1% in the spatial regimes crucial for RE seed formation. We find that the mode amplitudes are predicted to be sufficiently large to permit the possibility of significant radial transport of REs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗
NASA NTRS1993

Time-dependent thermocapillary convection in a Cartesian cavity - Numerical results for a moderate Prandtl number fluid

The present numerical simulation explores a thermal-convective mechanism for oscillatory thermocapillary convection in a shallow Cartesian cavity for a Prandtl number 6.78 fluid. The computer program developed for this simulation integrates the two-dimensional, time-dependent Navier-Stokes equations and the energy equation by a time-accurate method on a stretched, staggered mesh. Flat free surfaces are assumed. The instability is shown to depend upon temporal coupling between large scale thermal structures within the flow field and the temperature sensitive free surface. A primary result of this study is the development of a stability diagram presenting the critical Marangoni number separating steady from the time-dependent flow states as a function of aspect ratio for the range of values between 2.3 and 3.8. Within this range, a minimum critical aspect ratio near 2.3 and a minimum critical Marangoni number near 20,000 are predicted below which steady convection is found.

Peltier, L. J.↗
NASA NTRS1998

Exact Solution to Stationary Onset of Convection Due to Surface Tension Variation in a Multicomponent Fluid Layer With Interfacial Deformation

Stationary onset of convection due to surface tension variation in an unbounded multicomponent fluid layer is considered. Surface deformation is included and general flux boundary conditions are imposed on the stratifying agencies (temperature/composition) disturbance equations. Exact solutions are obtained to the general N-component problem for both finite and infinitesimal wavenumbers. Long wavelength instability may coexist with a finite wavelength instability for certain sets of parameter values, often referred to as frontier points. For an impermeable/insulated upper boundary and a permeable/conductive lower boundary, frontier boundaries are computed in the space of Bond number, Bo, versus Crispation number, Cr, over the range 5 x 10(exp -7) less than or equal to Bo less than or equal to 1. The loci of frontier points in (Bo, Cr) space for different values of N, diffusivity ratios, and, Marangoni numbers, collapsed to a single curve in (Bo, D(dimensional variable)Cr) space, where D(dimensional variable) is a Marangoni number weighted diffusivity ratio.

Skarda, J. Raymond Lee↗
DOE OSTI2023

The saturation mechanism of thermal instability

The literature on thermal instability (TI) reveals that even for a simple homogeneous plasma, the nonlinear outcome ranges from a gentle reconfiguration of the initial state to an explosive one, depending on whether the condensations that form evolve in an isobaric or nonisobaric manner. After summarizing the recent developments on the linear and nonlinear theory of TI, here we derive several general identities from the evolution equation for entropy that reveal the mechanism by which TI saturates; whenever the boundary of the instability region (the Balbus contour) is crossed, a dynamical change is triggered that causes the comoving time derivative of the pressure to change the sign. This event implies that the gas pressure force reverses direction, slowing the continued growth of condensation. For isobaric evolution, this “pressure reversal” occurs nearly simultaneously for every fluid element in condensation and a steady state is quickly reached. For nonisobaric evolution, the condensation is no longer in mechanical equilibrium and the contracting gas rebounds with greater force during the expansion phase that accompanies the gas reaching the equilibrium curve. The cloud then pulsates because the return to mechanical equilibrium becomes wave mediated. We show that both the contraction rebound event and subsequent pulsation behavior follow analytically from an analysis of the new identities. Our analysis also leads to the identification of an isochoric TI zone and makes it clear that unless this zone intersects the equilibrium curve, isochoric modes can only become unstable if the plasma is in a state of thermal non-equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

High resolution numerical simulations of methane pool fires using adaptive mesh refinement

The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. Here, in this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.

42 ENGINEERING↗
DOE OSTI2026

Axisymmetric hydrodynamics in numerical relativity: treating coordinate singularity, artificial heating and modeling MHD instabilities

Two-dimensional axisymmetric simulations of binary neutron star (BNS) merger remnant are a cheap alternative to 3D simulations. To maintain realism for secular timescales, simulations must avoid accumulated errors from drifts in conserved quantities and artificial heating, and they must model turbulent transport in a way that remains plausible throughout the evolution. It is also crucial to avoid numerical artifacts due to the polar coordinate axis singularity. Methods that behave well near the axis often break flux-conservative form of the hydrodynamic equations, resulting in significant drifts in conserved quantities. We present a flux-conservative scheme that maintains smoothness near the axis without sacrificing conservative formulation of the equations or incurring drifts in conserved global quantities. We compare the numerical performance of different treatments of the hydrodynamic equations when evolving a hypermassive neutron star resembling the remnant of a BNS merger. These simulations demonstrate that the new scheme combines the axis smoothness of non-conservative methods with the mass and angular momentum conservation of other conservative methods on $\sim 10^2$ ms timescales of viscous and neutrino-driven evolution. Because fluid profiles remain smooth in the remnant interior, it is possible to remove artificial heating by evolving the entropy density. We show how physical heating and cooling terms can be easily calculated from source terms of the conservative evolution variables and demonstrate our implementation. Finally, we discuss and implement improvements to the effective viscosity scheme to better model the effect of magnetohydrodynamic instabilities as the remnant evolves.

axisymmetry↗
NASA NTRS2000

Flap Edge Aeroacoustic Measurements and Predictions

An aeroacoustic model test has been conducted to investigate the mechanisms of sound generation on high-lift wing configurations. This paper presents an analysis of flap side-edge noise, which is often the most dominant source. A model of a main element wing section with a half-span flap was tested at low speeds of up to a Mach number of 0.17, corresponding to a wing chord Reynolds number of approximately 1.7 million. Results are presented for flat (or blunt), flanged, and round flap-edge geometries, with and without boundary-layer tripping, deployed at both moderate and high flap angles. The acoustic database is obtained from a Small Aperture Directional Array (SADA) of microphones, which was constructed to electronically steer to different regions of the model and to obtain farfield noise spectra and directivity from these regions. The basic flap-edge aerodynamics is established by static surface pressure data, as well as by Computational Fluid Dynamics (CFD) calculations and simplified edge flow analyses. Distributions of unsteady pressure sensors over the flap allow the noise source regions to be defined and quantified via cross-spectral diagnostics using the SADA output. It is found that shear layer instability and related pressure scatter is the primary noise mechanism. For the flat edge flap, two noise prediction methods based on unsteady-surface-pressure measurements are evaluated and compared to measured noise. One is a new causality spectral approach developed here. The other is a new application of an edge-noise scatter prediction method. The good comparisons for both approaches suggest that much of the physics is captured by the prediction models. Areas of disagreement appear to reveal when the assumed edge noise mechanism does not fully define, the noise production. For the different edge conditions, extensive spectra and directivity are presented. Significantly, for each edge configuration, the spectra for different flow speeds, flap angles, and surface roughness were successfully scaled by utilizing aerodynamic performance and boundary layer scaling method developed herein.

Brooks, Thomas F.↗
NASA NTRS2000

Flap Edge Aeroacoustic Measurements and Predictions

An aeroacoustic model test has been conducted to investigate the mechanisms of sound generation on high-lift wing configurations. This paper presents an analysis of flap side-edge noise, which is often the most dominant source. A model of a main element wing section with a half-span flap was tested at low speeds of up to a Mach number of 0.17, corresponding to a wing chord Reynolds number of approximately 1.7 million. Results are presented for flat (or blunt), flanged, and round flap-edge geometries, with and without boundary-layer tripping, deployed at both moderate and high flap angles. The acoustic database is obtained from a Small Aperture Directional Array (SADA) of microphones, which was constructed to electronically steer to different regions of the model and to obtain farfield noise spectra and directivity from these regions. The basic flap-edge aerodynamics is established by static surface pressure data, as well as by Computational Fluid Dynamics (CFD) calculations and simplified edge flow analyses. Distributions of unsteady pressure sensors over the flap allow the noise source regions to be defined and quantified via cross-spectral diagnostics using the SADA output. It is found that shear layer instability and related pressure scatter is the primary noise mechanism. For the flat edge flap, two noise prediction methods based on unsteady surface pressure measurements are evaluated and compared to measured noise. One is a new causality spectral approach developed here. The other is a new application of an edge-noise scatter prediction method. The good comparisons for both approaches suggest that much of the physics is captured by the prediction models. Areas of disagreement appear to reveal when the assumed edge noise mechanism does not fully define the noise production. For the different edge conditions, extensive spectra and directivity are presented. Significantly, for each edge configuration, the spectra for different flow speeds, flap angles, and surface roughness were successfully scaled by utilizing aerodynamic performance and boundary layer scaling methods developed herein.

Brooks, Thomas F.↗