Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TAYLOR INSTABILITY”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Laboratory study of the failed torus mechanism in arched, line-tied, magnetic flux ropes

Coronal mass ejections (CMEs) are some of the most energetic and violent events in our solar system. The prediction and understanding of CMEs are of particular importance due to the impact that they can have on Earth-based satellite systems and, in extreme cases, ground-based electronics. CMEs often occur when long-lived magnetic flux ropes (MFRs) anchored to the solar surface destabilize and erupt away from the Sun. One potential cause for these eruptions is an ideal magnetohydrodynamic (MHD) instability, such as the kink or torus instability. Previous experiments on the magnetic reconnection experiment revealed a class of MFRs that were torus-unstable but kink-stable, which failed to erupt. These “failed-tori” went through a process similar to Taylor relaxation, where the toroidal current was redistributed before the eruption ultimately failed. Herein we have investigated this behavior through additional diagnostics that measure the current distribution at the foot points and the energy distribution before and after an event. These measurements indicate that ideal MHD effects are sufficient to explain the energy distribution changes during failed torus events. This excludes Taylor relaxation as a possible mechanism of current redistribution during an event. A new model that only requires non-ideal effects in a thin layer above the electrodes is presented to explain the observed phenomena. This work broadens our understanding of the stability of MFRs and the mechanism behind the failed torus through the improved prediction of the torus instability and through new diagnostics to measure the energy inventory and current profile at the foot points.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Parameter space mapping of the Princeton magnetorotational instability experiment

Extensive simulations of the Princeton Magnetorotational Instability (MRI) Experiment with the Spectral/Finite Element code for Maxwell and Navier-Stokes Equations (SFEMaNS) have been performed to map the MRI-unstable region as a function of inner cylinder angular velocity and applied vertical magnetic field. The angular velocities of the outer cylinder and the end-cap rings follow the inner cylinder in fixed ratios optimized for MRI. Here, we first confirm the exponential growth of the MRI linear phase using idealized conducting vertical boundaries (end caps) rotating differentially with a Taylor-Couette profile. Subsequently, we run a multitude of simulations to scan the experimental parameter space and find that the normalized volume-averaged mean-square radial magnetic field, our main instability indicator, rises significantly where MRI is expected. At various locations, the local radial components of fluid velocity and generated magnetic field are well correlated with the volume-averaged indicator. Based on this correlation, a diagnostic system that will measure the radial magnetic field at several locations on the inner cylinder is proposed as the main comparison between simulation and experiment. A detailed analysis of poloidal mode structures in the SFEMaNS code indicates that MRI, rather than Ekman circulation or Rayleigh instability, dominates the fluid behavior in the region where MRI is expected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Instabilities and transition in the wall boundary layers of low-disturbance supersonic nozzles

Linear stability analysis and experimental nozzles were employed in studying means to control Tollmein-Schlichting (TS) waves and Taylor-Goertler (TG) vortices at transition points along walls in a supersonic wind tunnel. It was hoped that control of the TS waves and TG vortices would cause the wall flows to remain laminar, thereby reducing the turbulence noise in the tunnel. The test nozzles injecting the waves and vertices had Mach numbers from 3 to 5. TG vortices triggered transition, while axisymmetric nozzle flows suppressed it. TS waves were controllable with favorable pressure graients at Mach 3.5. Finally, advances were made in the capability of predicting the locations of the transition points.

Chen, F.-J.↗
NASA NTRS1992

The conductive propagation of nuclear flames. I - Degenerate C + O and O + Ne + Mg white dwarfs

The paper determines the physical properties - speed, width, and density structure - of conductive burning fronts in degenerate carbon-oxygen (C + O) and oxygen-neon-magnesium (O + Ne + Mg) compositions for a grid of initial densities and compositions. The dependence of the physical properties of the flame on the assumed values of nuclear reaction rates, the nuclear reaction network employed, the thermal conductivity, and the choice of coordinate system are investigated. The occurrence of accretion-induced collapse of a white dwarf is found to be critically dependent on the velocity of the nuclear conductive burning front and the growth rate of hydrodynamic instabilities. Treating the expanding area of the turbulent burning region as a fractal whose tile size is identical to the minimum unstable Rayleigh-Taylor wavelength, it is found, for all reasonable values of the fractal dimension, that for initial C + O or O + Ne + Mg densities above about 9 x 10 exp 9 g/cu cm the white dwarf should collapse to a neutron star.

Timmes, F. X.↗
DOE OSTI2023

On the relationship between the multi-region relaxed variational principle and resistive inner-layer theory

We show that the variational energy principle of the multi-region relaxed magnetohydrodynamic (MRxMHD) model can be used to predict finite-pressure linear tearing instabilities. In this model, the plasma volume is sliced into sub-volumes separated by 'ideal interfaces', and in each volume the magnetic field relaxes to a Taylor state, where the pressure gradient $\nabla p = 0$. The MRxMHD model is implemented in the Stepped-Pressure Equilibrium Code (SPEC) so that the equilibrium solution in each region is computed while preserving the force balance across the interfaces. As SPEC computes the Hessian matrix (a discretized stability matrix), the stability of an MRxMHD equilibrium can also be computed with SPEC. In this article, using SPEC, we investigate the effect of local pressure gradients and the $\nabla p = 0$ in the vicinity of the resonant surface of a tearing mode. For low-beta plasma, we have been able to illustrate a relationship between the resistive singular-layer theory (Coppi et al 1966 Nucl. Fusion6 101; Glasser et al 1975 Phys. Fluids18 875–88) and the MRxMHD model. Within the singular layer, the volume-averaged magnetic helicity and the flux-averaged toroidal flux are shown to be the invariants for the linear tearing modes in SPEC simulations. Our technique to compute MRxMHD stability is first tested numerically in a cylindrical tokamak and its application in toroidal geometry is demonstrated. We demonstrate an agreement between the stability boundary obtained with SPEC simulation and the resistive inner-layer theories.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2011

Effects of 30-, 60-, and 90-Day Bed Rest on Postural Control in Men and Women

INTRODUCTION Head-down-tilt bed rest (HDT) has been used as a safe gr ound-based analog to mimic and develop countermeasures for the physiological effects of spaceflight, including decrements in postural stability. The purpose of this investigation was to characterize the effects of 30-, 60-, and 90-day bed rest on postural control in men and women. METHODS Twenty-nine subjects (18M,11F) underwent 13 days of ambula tory acclimatization and were placed in 6? HDT for 30 (n=12), 60 (n=8), or 90 (n=9) days, followed by 14 days of ambulatory recovery. Computerized dynamic posturography (CDP) was used to assess changes in sensory and motor components of postural control, and recovery after HDT. Sensory Organization Tests (SOTs) objectively evaluate one?s ability to effectively use or suppress visual, vestibular, and proprioceptive information for postural control. Stability during the SOTs was assessed using peak-to-peak sway and convergence toward stability limits to derive an equilibrium score. Motor Control Tests (MCTs) evaluate one?s ability to recover from unexpected support surface perturbations, with performance determined by center-of-pressure path length. Whole-body kinematic data were collected to determine body-sway strategy used to maintain stability during each condition. Baselines were determined pre-HDT. Recovery was tracked post-HDT on days 0, 1, 2, and 4. RESULTS Immediately after HDT, subjects showed decreased performance on most SOTs, primarily on sway-referenced support conditions, typically returning to baseline levels within 4 days. MCT performance was not significantly affected. There were no significant gender or duration differences in performance. Kinematic data revealed a tendency to use ankle strategy to maintain an upright stance during most SOT conditions. Interestingly, six subjects (2M,4F) experienced orthostatic intolerance and were unable to complete day 0 testing. CONCLUSION HDT mimics some un loading mechanisms of spaceflight and elicits orthostatic issues present post-spaceflight (contributing to instability); however, it does not sufficiently address the vestibular dysfunction which occurs post-spaceflight.

Esteves, Julie↗
DOE OSTI2021

Modeling and simulation of transitional Rayleigh–Taylor flow with partially averaged Navier–Stokes equations

In this work, the partially averaged Navier–Stokes (PANS) equations are used to predict the variable-density Rayleigh–Taylor (RT) flow at Atwood number 0.5 and maximum Reynolds number 500. This is a prototypical problem of material mixing, featuring laminar, transitional, and turbulent flow, instabilities and coherent structures, density fluctuations, and production of turbulence kinetic energy by both shear and buoyancy mechanisms. These features pose numerous challenges to modeling and simulation, making the RT flow ideal to develop the validation space of the recently proposed PANS Besnard–Harlow–Rauenzahn-linear eddy viscosity model closure. The numerical simulations are conducted at different levels of physical resolution and test three approaches to set the parameters $f_\phi$ defining the range of physically resolved scales. The computations demonstrate the efficiency (accuracy vs cost) of the PANS model predicting the spatiotemporal development of the RT flow. Results comparable to large-eddy simulations and direct numerical simulations are obtained at significantly lower physical resolution without the limitations of the Reynolds-averaged Navier–Stokes equations in these transitional flows. The data also illustrate the importance of appropriate selection of the physical resolution and the resolved fraction of each dependent quantity $\phi$ of the turbulent closure, $f_\phi$. These two aspects determine the ability of the model to resolve the flow phenomena not amenable to modeling by the closure and, as such, the computations’ fidelity.

Navier Stokes equations↗
NASA NTRS1995

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.↗
NASA NTRS1986

Analytic evaluation of two-center molecular integrals

By using the Fourier-transform technique, the explicit expressions for the one-electron - two-center overlap integrals of Slater-type atomic orbitals up to 3d are derived. The final expressions are analytic, simple, and independent of local coordinates. Furthermore, they do not contain the nonclosed-form of exponential integrals which were presented in expressions given in earlier work. It is shown that the two-electron - two-center Coulomb integrals, as well as the hybrid integrals, can simply be expressed in terms of these integrals. The numerical instability arising from the situation in which the exponents of the two orbitals are almost equal is discussed, and a solution for this problem based on a Taylor-series expansion of the integral is suggested.

Tai, H.↗
NASA NTRS1993

A mesoscale gravity-wave event observed during CCOPE. IV - Stability analysis and Doppler-derived wave vertical structure

A summary of the results of a detailed study of the vertical structure of mesoscale gravity waves conducted during the Cooperative Convective Precipitation Experiment (CCOPE) is presented. Pressure perturbation fields derived from the Doppler wind fields are compared with the vertical structure of eigenfunctions resulting from a solution to the Taylor-Goldstein linear wave equation for an atmosphere whose mean state is described by vertical profiles obtained from a representative CCOPE sounding. An analysis of the potential for shear instability is also performed on all of the soundings taken on this day to assess the representativeness of the one chosen for the linear theoretical analysis.

Koch, Steven E.↗
DOE OSTI2021

Simulation of dynamic crystal plasticity with a Lagrangian discontinuous Galerkin hydrodynamic method

Here we present a new Lagrangian modal discontinuous Galerkin (DG) hydrodynamic method that supports a dynamic dislocation based crystal plasticity model for simulating the mechanical behavior of crystallographic materials, both single crystal and polycrystalline, under dynamic conditions. A modal DG approach is used to evolve fields relevant to conservation laws. These fields are approximated by Taylor series polynomials of varying degree. These polynomials describe macro-scale hydrodynamic behavior while their evolution is determined by evaluating the dynamic crystal plasticity model at material points within the element. The dynamic crystal plasticity model is sensitive to the time increment size, with too large of time increments leading to instability in the model. To mitigate this, the temporal evolution of the dynamic crystal plasticity model is achieved with the combination of a sub-incrementing scheme with Heun’s third-order time integration scheme, which is also used to temporally evolve the governing equations. The implementation of the dynamic crystal plasticity model within the DG framework is tested using a 2D approximation of the Taylor impact test with a single crystal material, using quadratic elements that have faces that can bend. In addition to the standard continuous material modeling, we propose a new simulation method that would represent the heterogeneous behavior of polycrystalline microstructures within an element by varying the position and material properties of the material points within the element. This method is demonstrated using random orientation distributions on materials points that are arranged in both structured and random configurations.

42 ENGINEERING↗
NASA NTRS1987

Spacelab experiments on convection in a rotating spherical shell with radial gravity

Data from the geophysical fluid flow cell experiment on Spacelab 3 revealed the diverse forms of fluid motion that can occur in geometrically simple models of global convection with rotation and radial gravity. The following types of convection were observed: midlatitude waves interacting with the low-latitude columnar convection rolls or 'banana cells', spiral waves' near the poles when latitudinal heating gradients are present on the bounding surfaces, and 'triangular waves' coupling midlatitude and equatorial disturbances under similar differential heating. It is believed that concepts based on the Taylor-Proudman theorem have a central role in many flows realized in the rapidly rotating experiments when H is approximately 0; however, the resulting banana cells are subject to secondary instabilities that produce intricate time dependence and eventually turbulent flows in which only hints of the simpler patterns are recognizable.

Toomre, Juri↗
NASA NTRS1978

Formation of a rotational accretion column

We consider the problem of the flow produced by the accretion of gas from a rotating cloud onto a point gravitational source. It is shown that, for a cloud initially in uniform rotation, the effect of accretion is to produce a Taylor column in the cloud parallel to the rotation axis. Only gas within the column is accreted; the gas outside the column undergoes oscillatory flow, which is probably subject to shear instability. The column forms in approximately one-half rotation period from the time at which accretion begins. Its initial radius is determined by the accretion rate and the angular velocity of the cloud. Application of the results to the formation of Jupiter by accretion is discussed.

Cassen, P.↗
NASA NTRS2009

Integrity of the Plasma Magnetic Nozzle

This report examines the physics governing certain aspects of plasma propellant flow through a magnetic nozzle, specifically the integrity of the interface between the plasma and the nozzle s magnetic field. The injection of 100s of eV plasma into a magnetic flux nozzle that converts thermal energy into directed thrust is fundamental to enabling 10 000s of seconds specific impulse and 10s of kW/kg specific power piloted interplanetary propulsion. An expression for the initial thickness of the interface is derived and found to be approx.10(exp -2) m. An algorithm is reviewed and applied to compare classical resistivity to gradient-driven microturbulent (anomalous) resistivity, in terms of the spatial rate and time integral of resistive interface broadening, which can then be related to the geometry of the nozzle. An algorithm characterizing plasma temperature, density, and velocity dependencies is derived and found to be comparable to classical resistivity at local plasma temperatures of approx. 200 eV. Macroscopic flute-mode instabilities in regions of "adverse magnetic curvature" are discussed; a growth rate formula is derived and found to be one to two e-foldings of the most unstable Rayleigh-Taylor (RT) mode. After establishing the necessity of incorporating the Hall effect into Ohm s law (allowing full Hall current to flow and concomitant plasma rotation), a critical nozzle length expression is derived in which the interface thickness is limited to about 1 ion gyroradius.

Gerwin, Richard A.↗
NASA NTRS1990

Convective instabilities in SN 1987A

Following Bandiera (1984), it is shown that the relevant criterion to determine the stability of a blast wave, propagating through the layers of a massive star in a supernova explosion, is the Schwarzschild (or Ledoux) criterion rather than the Rayleigh-Taylor criterion. Both criteria coincide only in the incompressible limit. Results of a linear stability analysis are presented for a one-dimensional (spherical) explosion in a realistic model for the progenitor of SN 1987A. When applying the Schwarzschild criterion, unstable regions get extended considerably. Convection is found to develop behind the shock, with a characteristic growth rate corresponding to a time scale much smaller than the shock traversal time. This ensures that efficient mixing will take place. Since the entire ejected mass is found to be convectively unstable, Ni can be transported outward, even into the hydrogen envelope, while hydrogen can be mixed deep into the helium core.

Benz, Willy↗
NASA NTRS1974

A passive submicrosecond time dissemination system

A submicrosecond time dissemination system was developed. The system is composed of an experimental transit improvement satellite (TRIAD), TRANET Satellite Tracking System, and special satellite signal receivers each tied to separate cesium clocks. The 400 MHz signal from the satellite is phase modulated plus or minus 45 deg by a pseudo-random noise (PRN) digital code with a bit rate of 5/3 MHz. Evaluation of this system has shown that satellite oscillator frequency instabilities are the major source of error. Experiments indicate that Global time transfers can be made with errors less than 100 nanoseconds.

Taylor, R. J.↗
NASA NTRS1989

Low-dimensional chaos in turbulence

Direct numerical simulations are being performed on two different fluid flows in an attempt to discover the mechanism underlying the transition to turbulence in each. The first system is Taylor-Couette flow; the second, two-dimensional flow over an airfoil. Both flows exhibit a gradual transition to high-dimensional turbulence through low-dimensional chaos. The hope is that the instabilities leading to chaos will be easier to relate to physical processes in this case, and that the understanding of these mechanisms can then be applied to a wider array of turbulent systems.

Vastano, John A.↗