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 811 records · Page 45

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Instability of baroclinic flows with horizontal shear along topography

The stability of baroclinic flows with horizontal shear over sloping topography is analyzed with special emphasis on the structure and energetics of the unstable perturbations. The study is conducted by using a linearized two-layer quasi-geostrophic channel model for different topography profiles and distributions of the basic velocity field. Interactions between the two fluid layers and the energy conversions by the unstable perturbations are described. It is found that topography sloping as (opposed to) the fluid interface contributes to enhance the perturbation amplitude in the upper (lower) layer relative to the lower (upper) layer. The results for bottom topography with differing characteristics across the flow indicate pronounced localized effects on the energy conversions over the slopes and the meridional scale of the perturbations in the lower layer.

Mechoso, C. R.↗

NASA NTRS1985

Fluid surface behavior in low gravity. Center discretionary fund no. 83-21

Measurements of rotating equilibrium bubble shapes in the low-gravity environment of a free-falling aircraft are presented. Emphasis is placed on bubbles which intersect the container boundaries. These data are compared with theoretical profiles derived from Laplace's formula and are in good agreement with the measurements. Two types of instability are explored. The first occurs when the baffle spacing is too large for the bubble to intersect both the top and bottom boundaries. The second occurs when the hydrostatic pressure beneath a displaced free surface does not compensate for pressure change due to capillary forces. The interface shape depends on the contact angle, the radius of intersection with container, and the parameter F which is a measure of the relative importance of centrifugal force to surface tension. For isolated bubbles, F has a maximum value of 1/2. A further increase in F causes the bubble to break contact with the axis of rotation. For large values of F, the bubble becomes more cylindrical and the capillary rise occurs over a thinner layer so that the small radius of curvature can generate enough pressure drop to balance the increased hydrostatic contribution.

Leslie, F.↗

NASA NTRS1991

Laboratory and theoretical models of planetary-scale instabilities and waves

Meteorologists and planetary astronomers interested in large-scale planetary and solar circulations recognize the importance of rotation and stratification in determining the character of these flows. The two outstanding problems of interest are: (1) the origins and nature of chaos in baroclinically unstable flows; and (2) the physical mechanisms responsible for high speed zonal winds and banding on the giant planets. The methods used to study these problems, and the insights gained, are useful in more general atmospheric and climate dynamic settings. Because the planetary curvature or beta-effect is crucial in the large scale nonlinear dynamics, the motions of rotating convecting liquids in spherical shells were studied using electrohydrodynamic polarization forces to generate radial gravity and centrally directed buoyancy forces in the laboratory. The Geophysical Fluid Flow Cell (GFFC) experiments performed on Spacelab 3 in 1985 were analyzed. The interpretation and extension of these results have led to the construction of efficient numerical models of rotating convection with an aim to understand the possible generation of zonal banding on Jupiter and the fate of banana cells in rapidly rotating convection as the heating is made strongly supercritical. Efforts to pose baroclinic wave experiments for future space missions using a modified version of the 1985 instrument have led us to develop theoretical and numerical models of baroclinic instability. Some surprising properties of both these models were discovered.

Hart, John E.↗

DOE OSTI2021

Modeling Snowflake Divertors in MAST-U Tokamak

We report that in a snowflake (SF) divertor, two magnetic field nulls are placed close to each other, creating four strike points (SPs) compared to two in a standard X-point divertor. In preparation for MAST-U experiments, magnetic configurations with the standard and SF divertors with various locations and separation distances of the nulls were modeled using the two-dimensional multi-fluid code UEDGE with a full plasma transport model featuring charge-state-resolved sputtered carbon impurities. The complex interplay of the plasma transport and magnetic configurations was comprehensively studied using a simple model for the theoretically predicted fast plasma mixing driven by the 'churning' mode instability in the two-null SF region. The modeling results show that (1) all SF-plus configurations and SF-minus configuration with closely located nulls produce the same plasma parameters and heat fluxes at the same SPs; (2) SF divertors approach the outer and inner SP detachment conditions at lower upstream density w.r.t. the standard divertor; (3) heat flux profiles at primary SPs are substantially broadened and peak values are reduced in SF configurations w.r.t. SN divertors; this broadening becomes more pronounced with the fast plasma mixing increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1982

Fluid flow in crystallization processes

Investigations of the interactions of heat, mass, and momentum transport in crystal growth from the melt are described. Studies of the mall-scale floating zone process and on a prototype of the vertical Bridgman growth system were emphasized. In both systems detailed numerical calculations are used to dissect the interplay between fluid convection and dopant segregation. These calculations are based on finite element techniques that make feasible the complete solution of solidification problems which include convection. Other numerical methods were developed for solving thermal models of crystal growth processes with melt/gas menisci and for simulating the microscale instabilities in solidification interfaces.

Brown, R. A.↗

NASA NTRS1991

The wavelength of supercritical surface tension driven Benard convection

The size or the wavelength of moderately supercritical surface tension driven Benard convection has been investigated experimentally in a thin fluid layer of large aspect ratio. It has been found that the number of the hexagonal convection cells increases with increased temperature differences, up to 1.3 times the critical temperature difference. That means that the wavelength of surface tension driven convection decreases after onset of the instability for moderately nonlinear conditions. This result is in striking contrast to the well-known increase of the wavelength of buoyancy driven Rayleigh-Benard convection.

Koschmieder, E. L.↗

NASA NTRS1993

Assessment of Total Variation Diminishing schemes in combustion instability problems

Numerical simulations of acoustic waves in a shear layer and in an idealized combustion chamber using high resolution Total Variation Diminishing (TVD) schemes has been carried out to study the effects of inherent scheme dissipation and dispersion errors of this class of problems. The numerical results are compared against available exact solutions to quantify these errors. Several popular TVD limiters widely used in the Computational Fluid Dynamics (CFD) community have been assessed. Osher-Chakravarthy limiters are modified so that they can be used in explicit schemes. Among all the limiters investigated, Osher-Chakravarthy third-order limiter is identified as having performed the best. It is also found that all TVD schemes have exceptionally small dispersive errors.

Wang, Z. J.↗

DOE OSTI2025

Generalizing the compressible pairwise interaction extended point-particle model

Ejecta physics plays an important role in material interfaces that are impacted by a strong shock wave. When a shock impacts a rough surface of solid material and melts it, the Richtmyer–Meshkov instability grows perturbations on the surface, which can eject particles. After release, the ejecta travel through the post-shock compressible flow. To accurately simulate a large number of ejecta particles, an Euler–Lagrange approach is preferred, which requires modeling the subgrid-scale physics involved with fluid–particle interactions. We generalize the previous work from Hsiao et al. (2023) to consider systems of moving particles subject to any loading shock. The following improvements were made: (1) Particles are allowed to move relative to each other (2) Non-planar shocks are accounted for along with allowing for variable shock speeds. As a result, the generalized algorithm was tested with particle-resolved simulations for canonical test cases. The results of these tests are discussed and analyzed.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1989

Dendrites, viscous fingers, and the theory of pattern formation

Recent developments in the theory of pattern formation in dendritic crystal growth and viscous fingering in fluids are reviewed. Consideration is given to the discovery that weak capillary forces act as singular perturbations which lead to selection mechanisms in dendritic crystal growth and fingering patterns. Other topics include the conventional thermodynamic model of the solidification of a pure substance from its melt, fingering instability, pattern selection, the solvability theory, dendritic growth rates, the bubble effect discovered by Couder et al. (1986), the dynamics of pattern-forming systems, and snowflake formation.

Langer, J. S.↗

NASA NTRS1990

Transition to hard turbulence in thermal convection at infinite Prandtl number

Direct numerical simulations of two-dimensional high Rayleigh (Ra) number, base-heated thermal convection in large aspect-ratio boxes are presented for infinite Prandtl number fluids, as applied to the earth's mantle. A transition is characaterized in the flow structures in the neighborhood of Ra between 10 to the 7th and 10 to the 8th. These high Ra flows consist of large-scale cells with strong intermittent, boundary-layer instabilities. For Ra exceeding 10 to the 7th it is found that the heat-transfer mechanism changes from one characterized by mushroom-like plumes to one consisting of disconnected ascending instabilities, which do not carry with them all the thermal anomaly from the bottom boundary layer. Plume-plume collisions become much more prominent in high Ra situations and have a tendency of generating a pulse-like behavior in the fixed plume. This type of instability represents a distinct mode of heat transfer in the hard turbulent regime. Predictions of this model can be used to address certain issues concerning the mode of time-dependent convection in the earth's mantle.

Hansen, Ulrich↗

NASA NTRS2019

Preliminary Results from an Experimental Assessment of a Natural Laminar Flow Design Method

A 5.2% scale semispan model of the new Common Research Model with Natural Laminar Flow (CRM-NLF) was tested in the National Transonic Facility (NTF) at the NASA Langley Research Center. The model was tested at transonic cruise flight conditions with Reynolds numbers based on mean aerodynamic chord ranging from 10 to 30 million. The goal of the test was to experimentally validate a new design method, referred to as Crossflow Attenuated NLF (CATNLF), which shapes airfoils to have pressure distributions that delay transition on wings with high sweep and Reynolds numbers. Additionally, the test aimed to characterize the NTF laminar flow testing capabilities, as well as establish best practices for laminar flow wind tunnel testing. Preliminary results regarding the first goal of validating the new design method are presented in this paper. Experimental data analyzed in this assessment include surface pressure data and transition images. The surface pressure data acquired during the test agree well with computational fluid dynamics (CFD) results. Transition images at a variety of Reynolds numbers and angles of attack are presented and compared to computational transition predictions. The experimental data are used to assess transition due to a turbulent attachment line, as well as crossflow and Tollmien-Schlichting modal instabilities. Preliminary results suggest the CATNLF design method is successful at delaying transition on wings with high sweep. Initial analysis of the transition front images showed transition Reynolds numbers that exceed historic experimental values at similar sweep angles. , section lift

Lynde, Michelle N.↗

NASA NTRS1997

Gamma-Ray Bursts and Relativistic MHD around Coalescing Neutron Stars

In this proposal we performed the first relativistic hydrodynamic simulations of the binary neutron stars near coalescence. Calculations were performed in three spatial dimensions and were supplemented with magnetohydrodynamic simulations of the magnetic reconnection as the neutron star fluid responds to the relativistic forces of the binary. The ultimate goal of the proposed work was to develop a theoretical model for the temporal and spectral evolution of cosmological gamma-ray bursts produced during the magnetic reconnection. This is the first time such hydrodynamic calculations were performed and some surprising results were obtained. The neutrons stars were observed to compress and heat in the simulations. They could even collapse to black holes many orbits before coalescence. The question then arose as to whether this strange phenomenon could power a gamma-ray burst. The bulk of computing effort in this proposal was therefore devoted to an analysis of the implications of this possible new gamma-ray burst paradigm. The initial focus of the work performed under this proposal was therefore to utilize (3+1) dimensional and spherical numerical general relativistic hydrodynamics to study the origin, evolution, and parametric sensitivity of the collapse instability. We derived modified conditions of hydrostatic equilibrium for stars in the curved space of quasi-static orbits. We followed the magnetic field evolution for approximately one orbit (10 msec). The magnetic field energy was observed to exponentiate with an e-folding time of about 1 msec. Thus, the field should build up very quickly to a magnitude such that reconnection and back reaction of the fluid inhibits further growth.

Mathews, Grant J.↗

NASA NTRS2002

Microphysics of Waves and Instabilities in the Solar Wind and Their Macro Manifestations in the Corona and Interplanetary Space

A breakthrough resulted from the investigation of plasma heating by ion-cyclotron waves in rapidly expanding flow tubes, often referred to as coronal funnels and expected to originate in the low transition region where the temperature gradient in the solar atmosphere rises very sharply. The investigation was made using a three-fluid plasma consisting of protons, electrons and alpha particles. It was found that these waves heat the solar wind plasma by directly heating the heavier species, namely the alphas. Although only alpha particles dissipate the waves, the strong Coulomb coupling between alpha particles and protons, and between protons and electrons, makes it possible for protons and electrons to be heated also to more than one million degrees, i.e. to coronal temperatures. Interestingly, the extreme heating of the alpha particles, however, is such that they end up being hotter and faster, and are no longer in thermal equilibrium with the protons and electrons. once the rapid expansion of the flow tube is complete, the particles return to thermal equilibrium. The observational signatures of these results are such that if spectral lines formed below a million degrees are observed at different heights in the corona, the inferred outflow velocities may vary by a factor of 5 to 6. In addition, if minor ions are indeed much faster than protons and electrons at temperatures below a million degrees, then one cannot reliably determine the bulk outflow velocity of the solar wind in that region from inferences of minor ion outflow velocities. A detailed parameter study on the resonant interaction between ion cyclotron waves and alpha particles and other minor ions is also under way. It is found that the highest frequency of the waves, or the location of the heating determines whether a thermal equilibrium can be established between minor ions and protons. As the maximum frequency is increased, alpha particles and protons become farther away from thermal equilibrium. As the maximum frequency becomes smaller, alpha particles and protons are more likely to be in thermal equilibrium. However, if the maximum frequency is below a critical value, ion cyclotron resonance occurs at a greater height and a hot corona cannot be produced.

Habbal, Shadia Rifai↗

NASA NTRS1979

Attenuation of sound in a low Mach number nozzle flow

This paper examines the energy conversion mechanisms which govern the emission of low frequency sound from an axisymmetric jet pipe of arbitrary nozzle contraction ratio in the case of low Mach number nozzle flow. The incident acoustic energy which escapes from the nozzle is partitioned between two distinct disturbances in the exterior fluid. The first of these is the free-space radiation, whose directivity is equivalent to that produced by monopole and dipole sources. Second, essentially incompressible vortex waves are excited by the shedding of vorticity from the nozzle lip, and may be associated with the large-scale instabilities of the jet. Two linearized theoretical models are discussed. One of these is an exact linear theory in which the boundary of the jet is treated as an unstable vortex sheet. The second assumes that the finite width of the mean shear layer of the real jet cannot be neglected. The analytical results are shown to compare favourably with recent attenuation measurements.

Howe, M. S.↗

NASA NTRS2016

Preferential Heating of Oxygen 5+ Ions by Finite-Amplitude Oblique Alfven Waves

Minor ions in the fast solar wind are known to have higher temperatures and to flow faster than protons in the interplanetary space. In this study we combine previous research on parametric instability theory and 2.5D hybrid simulations to study the onset of preferential heating of Oxygen 5+ ions by large-scale finite-amplitude Alfven waves in the collisionless fast solar wind. We consider initially non-drifting isotropic multi-species plasma, consisting of isothermal massless fluid electrons, kinetic protons and kinetic Oxygen 5+ ions. The external energy source for the plasma heating and energization are oblique monochromatic Alfven-cyclotron waves. The waves have been created by rotating the direction of initial parallel pump, which is a solution of the multi-fluid plasma dispersion relation. We consider propagation angles theta less than or equal to 30 deg. The obliquely propagating Alfven pump waves lead to strong diffusion in the ion phase space, resulting in highly anisotropic heavy ion velocity distribution functions and proton beams. We discuss the application of the model to the problems of preferential heating of minor ions in the solar corona and the fast solar wind.

Maneva, Yana G.↗

DOE OSTI2023

Energy transport analysis of NSTX plasmas with the TGLF turbulent and NEO neoclassical transport models

This work presents a study of plasma transport at low aspect ratio on the National Spherical Torus Experiment tokamak, where the turbulent and neoclassical energy fluxes calculated by the quasilinear Trapped Gyro Landau Fluid (TGLF) model and the multi species drift-kinetic Neoclassical solver (NEO) are validated against experimental data. The turbulent energy transport of two plasma discharges, one in the L-mode confinement regime and another in the H-mode regime, is dominated by electrostatic drift-wave instabilities, while the ion heat transport has a significant neoclassical contribution. The data analysis workflow is described in detail to understand how the variations of mapping and fitting of experimental data affect the power balance solution and subsequent flux-matching plasma profile predictions with the TGYRO solver. On average, the predicted plasma profiles are consistent with experimental data. However, the solutions are sensitive to various input parameters, including boundary conditions, and the electron-ion coupling. Linear gyrokinetic stability analysis demonstrates close agreement of the real frequencies of unstable modes between TGLF and CGYRO gyrokinetic simulations, but higher growth rates are predicted by TGLF, especially for the H-mode case. Estimates of the low-k modes' contributions to the total flux are consistent with linear stability analysis and the E × B suppression of turbulence in TGLF simulations with the SAT1 saturation model, while the SAT2 saturation model over-predicts the low-k modes' contribution in the H-mode case. Moreover, the results with SAT1 model are consistent with power balance analysis, which indicates only neoclassical ion energy fluxes inside ρ < 0.4 in the L-mode case and $\rho \unicode{x2A7D} 0.7$ in the H-mode case. The presence of multi-scale turbulence and ion-scale driven zonal flow mixing effects are also observed in TGLF scans of the electron turbulent heat flux over a range of temperature gradients and the electron-ion temperature ratio, which could explain the strong model sensitivity to variations of input parameters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1978

Evidence for two-dimensional inertial turbulence in a cosmic-scale low-beta plasma

Magnetospheric electric-field power spectral densities measured from high-altitude balloons and the Hawkeye 1 satellite are combined in wavenumber space and exhibit a power-law spectrum with index -1.6 + or -0.3 below a wavenumber of 0.05 reciprocal km and index -2.8 + or -0.3 above 2 reciprocal km. This spectrum agrees with that predicted for the inertial subrange of an isotropic homogeneous two-dimensional fluid and for a two-dimensional plasma when energy enters the system at the spectral knee and there exists a viscous dissipation mechanism at large wavenumbers. The wavelength of the spectral knee corresponds to that of folds and curls in the aurora and may be due to an instability of the auroral particle beams.

Kelley, M. C.↗

NASA NTRS1980

Soliton solutions and their stability for the flow of relativistic fluids through channels

The flow of a perfect relativistic fluid through channels of various cross-sections is considered with reference to models of radio galaxies. Soliton-like solutions are found and their topologies are discussed. The calculations show that these solutions are unstable. It is suggested that under realistic astrophysical conditions the growth rate of the instabilities is so slow that soliton-type blobs may persist for a significant time.

Lerche, I.↗