Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “harmonic 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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

MARS-Q modeling of kink-peeling instabilities in DIII-D QH-mode plasma

In quiescent H-mode (QH-mode) regime, edge harmonic oscillations (EHOs) are believed to provide necessary radial transport to prevent occurrence of large edge localized modes. A systematic modeling study is performed here on the low- n EHOs in a DIII-D QH-mode plasma (Chen et al 2016 Nucl. Fusion 56 076011), by utilizing the MARS-Q code (Liu et al 2013 Phys. Plasmas 20 042503). Both the n = 1 and n = 2 instabilities are found to be strongly localized near the plasma edge, exhibiting the edge-peeling characteristics. The DIII-D resistive wall is found to have minor effects on these instabilities. The plasma resistivity is found to strongly modify the mode growth rate. Assuming the Spitzer model for the plasma resistivity, the computed mode growth rate scales as S sup>-1/3 with S being the Lundquist number. Toroidal flow of the plasma slightly stabilizes these edge localized kink-peeling modes. Drift kinetic effects all have a destabilization effect on these modes. Non-perturbative magneto-hydrodynamic-kinetic hybrid computations find that the drift kinetic effects associated with thermal particle species push the peak location of the eigenmode radially inward but still in the pedestal region. The modeled plasma temperature and density fluctuations in the plasma edge region, as well as the poloidal magnetic field perturbations along both the low and high field sides of the plasma surface, are in good agreement with experimental measurements. Finally, the quasi-linear initial value simulations find a strong non-linear interplay between the kink-peeling instability and the toroidal flow near the plasma edge. The combined effect of the damping of the flow amplitude and change of the edge flow shear is found to be the stabilizing factor for the kink-peeling mode, leading to the mode saturation and thus EHOs.

Physics↗
NASA NTRS1978

Incompressible fluid ellipsoids in halos. I - The second-harmonic oscillations of the Maclaurin spheroids

The structure and stability of Maclaurin spheroids embedded in rigid uniform-density oblate spheroidal halos are determined by the tensor virial-equation method. These spheroid-halo systems can be thought of as crude fluid analogs of disk galaxies with halos. The halos are assumed to have the same center, the same axis of symmetry, and the same equatorial radius as the Maclaurin spheroids. Only halos with lower eccentricity than the Maclaurin spheroids are considered. The dynamic instability of the toroidal (barlike) modes is suppressed when m, the ratio of the halo mass to Maclaurin spheroid mass, is greater than 3 pi/8 for spherical halos and when m is greater than 1/2 for halos congruent to the Maclaurin spheroids. Intermediate halo flattenings yield intermediate critical m-values. On the other hand, a neutral point of the toroidal modes in the rotating and inertial frames occurs for all m and for all allowed halo flattenings. Growth rates for secular instability beyond the neutral point are calculated, and the eigenfrequencies of all second-harmonic modes are given for select cases. The Ostriker-Peebles (1973) conjecture concerning the stability of disk galaxies against barlike perturbations appears to be incorrect.

Durisen, R. H.↗
NASA NTRS1985

Generation of capillary instabilities by external disturbances in a liquid jet

The receptivity problem in a circular liquid jet is considered. A time harmonic axial pressure gradient is imposed on the steady, parallel flow of a jet of liquid emerging from a circular duct. Using a technique developed in plasma physics a casual solution to the forced problem is obtained over certain ranges of Weber number for a number of mean velocity profiles. This solution contains a term which grows exponentially in the downstream direction and can be identified with a capillary instability wave. Hence, it is found that the externally imposed disturbances can indeed trigger instability waves in a liquid jet. The amplitude of the instability wave generated relative to the amplitude of the forcing is computed numerically for a number of cases.

Leib, S. J.↗
DOE OSTI2022

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E × B) flow relative to the unmagnetized ions can cause the so-called electron cyclotron drift instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in, for example, collisionless shock ramps in space, and in E × B discharge devices such as Hall thrusters. A prominent feature of ECDI is its capability to induce an electron flow parallel to the background >E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport is important due to its role in particle thermalization at space shocks, and in causing plasma flows towards the walls of E × B devices, leading to unfavorable erosion and performance degradation, etc. The development of ECDI and anomalous transport is often considered requiring a fully kinetic treatment. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a collisionless two-fluid framework without any adjustable collision parameter. By treating both electron and ion species on an equal footing, the free energy due to the inter-species velocity shear allows the growth of an anomalous electron flow parallel to the background E field. We will first present linear analyses of the instability in the two-fluid five- and ten-moment models, and compare them against the fully-kinetic theory. At low temperatures, the two-fluid models predict the fastest-growing mode in good agreement with the kinetic result. Also, by including more (> = 10) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on ion-to-electron mass ratio, plasma temperature, and background B field strength is also thoroughly explored. We then carry out direct numerical simulations of the cross-field setup using the five-moment model. The development of the instability, as well as the anomalous transport, is confirmed and in excellent agreement with theoretical predictions. The force balance properties are also studied using the five-moment simulation data. Here this work casts new insights into the nature of ECDI and the associated anomalous transport and demonstrates the potential of the two-fluid moment model in efficient modeling of E × B plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

3D Simulations of the Richtmyer-Meshkov Instability with Re-Shock

We present results of inviscid simulations, in three dimensions, of Richtmyer-Meshkov instability for high incident shock Mach number. The growth rate of a single harmonic perturbation is quantified and compared with the results of a 2D calculation. Upon re-shock, the perturbation amplitude undergoes a phase reversal while the mean velocity of the interface is zero. Before re-shock the normalized growth rate of a 2D and 3D interface are nearly the same, but the growth rate after re-shock is significantly larger for the 3D than the 2D case. We also examine the evolution of multiple harmonic perturbations. Computational and parallelization issues of the simulation code will also be briefly discussed. The computations were done on the T3E at Pittsburgh Supercomputing Center.

Meiron, Daniel I.↗
NASA NTRSDate not supplied

Separation and Transition on a Cone-Cylinder-Flare: Computational Investigations

Base flow computation and stability analysis were conducted for hypersonic flow over a cone-cylinder-flare (CCF) geometry for conditions that correspond to the experimental runs carried out in three wind tunnels. Owing to the presence of an attached boundary layer, a separation bubble induced by a shock-boundary layer interaction, and a reattachment region, the chosen flow configuration is physically rich. The complexity of this flowfield encompasses a combination of convective instabilities developing on the cone, global instabilities in the separation bubble, and shear-layer modes and streaks in the reattachment region. Thus, the selected CCF configuration provides the opportunity for a comprehensive comparison of the currently available methodologies for analyzing boundary layer instabilities. Various tools are used for the analysis, including global stability codes as well as convective instability analyses based on a local theory, a weakly non-parallel analysis, and tools that are applicable to strongly non-parallel flows. The paper presents a comparison of the convective instability characteristics based on different methodologies, such as linear stability theory (LST), the harmonic form of linearized Navier-Stokes equations (HLNSE), and resolvent analysis. The CCF configuration provided an effective framework for conducting a detailed cross-validation of this type, which had not yet been addressed in existing literature. This document is accompanied by a companion paper that is focused on the experimental aspects of the CCF configuration. Both papers are being presented in a dedicated session that highlights the research activities of the NATO STO Research Task Group AVT-346.

Boundary Layer Transition↗
NASA NTRS1985

Electron cyclotron maser instability in the solar corona - The role of superthermal tails

The effect of a superthermal component of electrons on the loss-cone-driven electron cyclotron maser instability is analyzed. It is found that for a superthermal tail with temperature about 10 KeV, the first harmonic (X- and O-mode) is suppressed for n(t)/n(r) of about 1 (n/t/ and n/r/ are the densities of superthermal tail and loss-cone electrons) and the second harmonic (X- and O-modes) is suppressed for n(t)/n(r) less than about 0.1. A qualitative discussion on the formation of superthermal tails is presented and it is suggested that superthermal tails play an important role on the observed or available power, at microwave frequencies, from the electron cyclotron maser instability in the solar corona.

Vlahos, L.↗
NASA NTRS2022

Combined Bluntness and Roughness Effects on Cones at Hypersonic Speeds

This computational study investigates the effects of discrete roughness elements on a blunt cone at zero degrees angle of attack in a Mach 6 flow. Motivation was provided by experiments conducted in the Air Force Research Laboratory Mach 6 High Reynolds Number facility on a 7-degree half-angle cone with a roughness array located at 45 degrees from the apex on two nose tips of different blutness but equivalent roughness Reynolds number. Transition was only affected on the blunter cone, indicating that the transition onset is associated with the combined effects of bluntness and roughness. The present study investigates the 15.24 mm nose radius, 420 azimuthal wavenumber case via Navier-Stokes computations of the laminar base flow and instability analysis. Plane-marching parabolized stability equations (PSE) and inflow-resolvent analysis based on the three-dimensional, harmonic linearized Navier-Stokes equations (HLNSE) are used to calculate the amplification of disturbances along the roughness wake as well as over the roughness nearfield. Results show that the roughness shape can have a great impact on the characteristics of the most amplified wake instabilities. For the experimental configuration with cubic roughness elements of 15 μ m height, the flow is marginally unstable. For prismatic elements of 20 μ m height, the PSE predicts a logarithmic disturbance amplification ratio of N = 5.6 along its wake, but this ratio increases to N = 9.6 when the amplification over the roughness and separation regions is included in the inflow-resolvent analysis.

Boundary-layer transition↗
NASA NTRS2022

Combined Bluntness and Roughness Effects on Cones at Hypersonic Speeds

This computational study investigates the effects of discrete roughness elements on a blunt cone at zero degrees angle of attack in a Mach 6 flow. Motivation was provided by experiments conducted in the Air Force Research Laboratory Mach 6 High Reynolds Number facility on a 7-degree half-angle cone with a roughness array located at 45 degrees from the apex on two nose tips of different blutness but equivalent roughness Reynolds number. Transition was only affected on the blunter cone, indicating that the transition onset is associated with the combined effects of bluntness and roughness. The present study investigates the 15.24 mm nose radius, 420 azimuthal wavenumber case via Navier-Stokes computations of the laminar base flow and instability analysis. Plane-marching parabolized stability equations (PSE) and inflow-resolvent analysis based on the three-dimensional, harmonic linearized Navier-Stokes equations (HLNSE) are used to calculate the amplification of disturbances along the roughness wake as well as over the roughness nearfield. Results show that the roughness shape can have a great impact on the characteristics of the most amplified wake instabilities. For the experimental configuration with cubic roughness elements of 15 μ m height, the flow is marginally unstable. For prismatic elements of 20 μ m height, the PSE predicts a logarithmic disturbance amplification ratio of N = 5.6 along its wake, but this ratio increases to N = 9.6 when the amplification over the roughness and separation regions is included in the inflow-resolvent analysis.

Boundary-layer transition↗
NASA NTRS2017

Critical Fluctuations in Beam-Plasma Systems and Solar Type III Radio Bursts

The type III radio bursts are the most intense radio emissions from the sun. In Figure 1, we present a typical type III burst observed by the STEREO spacecraft. Ginzburg and Zheleznyakov (1958) were the first to suggest that Langmuir waves excited by the solar flare accelerated electrons are the source of these bursts. The in situ detection of electron beams and Langmuir waves in association with type III bursts together with tracking of type III burst sources in the interplanetary medium confirmed this plasma hypothesis. However, the dynamics of type III electron beams and their interaction with the ambient plasma through the excited Langmuir waves is an unsolved problem. In order that the electron beams don't lose their energy by resonantly interacting with Langmuir waves, various nonlinear processes, such as the induced scattering of Langmuir waves off the background ions in weak turbulence regime and oscillating two-stream instability in strong turbulence regime were invoked. However, the observed electron density fluctuations do not allow Langmuir waves to grow to very high intensities as pointed out by several authors. As far as the conversion of Langmuir wave energy into electromagnetic energy is concerned, there is no consensus either for the fundamental or for the second harmonic emission. In this paper, we address these issues in terms of critical fluctuations near the boundary of the beam-plasma instability which is analogous to the anomalously growing fluctuations near the phase transition points called critical points. We show that this phenomenon can account for the survival of electron beams over large distances as well as the intensities of Langmuir waves excited during this phase transition regime are sufficient to explain the peak intensities of the fundamental as well as harmonic emissions

Fluctuations↗
NASA NTRS2022

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the non-local, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, as roughness elements, steps, gaps, or corners. To solve the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 107 to 108. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear system. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE.

Boundary Layer Stability↗
NASA NTRS2022

The Harmonic Linearized Navier-Stokes Equations for Transition Prediction in Three-Dimensional Flows

The conventional method to predict the onset of laminar-turbulent transition in convectively unstable boundary-layer flows is based on the logarithmic amplification ratio, the so-called N-factor, of the linear instability waves. To calculate the N-factor, the flow variables are decomposed into a laminar basic state solution and the linear disturbances, which are assumed to be harmonic in time. The most commonly used linear stability analysis approaches include the locally parallel linear stability theory (LST) and the nonlocal, weakly nonparallel parabolized stability equations (PSE). However, these methods do not account for strong streamwise gradients that are encountered in several configurations of interest, such as those in the vicinity of roughness elements, steps, gaps, or corners. To compute the linear evolution of disturbances along such strongly nonparallel regions, the harmonic linearized Navier-Stokes equations (HLNSE) need to be solved. The discretization of the HLNSE for spanwise/azimuthally inhomogeneous laminar basic states yields a linear system of complex arithmetic with a leading dimension of the order of 10^(7) to 10^(8) even in relatively simple flows. A combined multithread and multiprocessor algorithm is implemented for the direct solution of such linear systems. Results for a supersonic boundary layer over a three-dimensional roughness patch show good agreement with experimental measurements when the evolution of the instability waves over the roughness patch is included via the HLNSE. Additionally, inflow-resolvent analysis based on the HLNSE for discrete-roughness-induced disturbances in the nose tip of a blunt cone at Mach 6 demonstrates the importance of including the disturbance amplification along the near vicinity of the roughness element and separation region.

Boundary Layer Stability↗
DOE OSTI2024

Center for Tokamak Transients Simulations

The major goals for the Center for Tokamak Transients Simulations (CTTS) were to develop detailed numerical simulations of disruptive transients in tokamak experiments and to use them in theoretical studies of the thermal-quench (TQ) and current-quench (CQ) phases of disruption. The University of Wisconsin-Madison component of CTTS contributed relevant development for the NIMROD code (https://nimrodteam.org) and applied it to disruptions that involve global vertical displacement. Code-development contributions include boundary conditions that represent sheath effects which occur when tokamak plasma is in contact with the surrounding vessel wall during vertical displacement events (VDEs). We also implemented a reduced model of energetic runaway electrons (REs) to simulate the effect of these energetic electrons on the CQ phase of tokamak disruptions. Our computations of VDEs include 2D and 3D verification through benchmarks with two other widely used macroscale plasma simulation codes, M3D-C1 and JOREK. Over nonlinear evolution to the final termination of plasma current, the 2D simulation results on magnetic axis location, toroidal current, and halo current track each other well. With the 3D benchmark computation, all three codes found the same qualitative behavior of the asymmetric instabilities having increasing growth rate after contact with the wall, followed by nonlinear excitation of other toroidal harmonics. Our NIMROD computations for an idealized configuration show agreement with an analytical assessment by V. D. Pustovitov [Nucl. Fusion 55, 113032 (2015)] that the plasma and non-ideal wall remain in force-balance such that net forces only result from magnetic stresses over the outer surface of the wall. Our study of the influence of boundary conditions on VDEs shows that extended-MHD simulations of VDEs are sensitive to electron energy transport to the wall and that boundary conditions on flow are only significant to the extent that they influence the electron energy transport. Resistive-MHD NIMROD simulations of the tokamak discharges in the Madison Symmetric Torus show that magnetic perturbations of poloidal wavenumber m=3 lead to chaotic magnetic topology over the edge region, which helps deconfine energetic electrons. A second area of RE study shows that the coupled resistive MHD/reduced RE model reproduces the resistive-hose particle-beam instability. With respect to educational opportunities, this effort involved and supported four graduate students at the University of Wisconsin-Madison and one postdoctoral associate. The results of this study have been disseminated through journal publications, conference presentations, technical reports, and PhD dissertations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Harmonic oscillators and resonance series generated by a periodic unstable classical orbit

The presence of an unstable periodic classical orbit allows one to introduce the decay time as a purely classical magnitude: inverse of the Lyapunov index which characterizes the orbit instability. The Uncertainty Relation gives the corresponding resonance width which is proportional to the Planck constant. The more elaborate analysis is based on the parabolic equation method where the problem is effectively reduced to the multidimensional harmonic oscillator with the time-dependent frequency. The resonances form series in the complex energy plane which is equidistant in the direction perpendicular to the real axis. The applications of the general approach to various problems in atomic physics are briefly exposed.

Kazansky, A. K.↗
NASA NTRS1990

Impingement tones of large aspect ratio supersonic rectangular jets

The frequency structure of large aspect ratio jet impingement tones is analyzed. Feedback loops consisting of downstream propagating instability waves of the jet flow and upstream propagating acoustic waves from the wall to the nozzle lip produce as many as 20 or more tones when a supersonic rectangular jet of large aspect ratio is directed normal on a wall. It is seen that there are only two basic tones. The basic tone of the lower frequency is associated with a symmetric instability feedback mode of the jet, while that of the higher frequency is associated with an antisymmetric instability feedback mode. The remaining tones are combination tones of the two basic frequencies and their harmonics. It is determined that the measured basic impingement tone frequencies fall in the theoretical permissible frequency bands over a wide range of jet Mach number.

Tam, Christopher K. W.↗
NASA NTRS1996

Hypersonic Boundary-Layer Stability Experiments on a Flared-Cone Model at Angle of Attack in a Quiet Wind Tunnel

An experimental investigation of the effects of angle of attack on hypersonic boundary-layer stability on a flared-cone model was conducted in the low-disturbance Mach-6 Nozzle-Test Chamber Facility at NASA Langley Research Center. This unique facility provided a 'quiet' flow test environment which is well suited for stability experiments because the low levels of freestream 'noise' minimize artificial stimulation of flow-disturbance growth. Surface pressure and temperature measurements documented the adverse-pressure gradient and transition-onset location. Hot-wire anemometry diagnostics were applied to identify the instability mechanisms which lead to transition. In addition, the mean flow over the flared-cone geometry was modeled by laminar Navier-Stokes computations. Results show that the boundary layer becomes more stable on the windward ray and less stable on the leeward ray relative to the zero-degree angle-of-attack case. The second-mode instability dominates the transition process at a zero-degree angle of attack, however, on the windward ray at an angle of attack this mode was completely stabilized. The less-dominant first-mode instability was slightly destabilized on the windward ray. Non-linear mechanisms such as saturation and harmonic generation are identified from the flow-disturbance bispectra.

Doggett, Glen P.↗