Search NASASearch

SEARCH · Search NASA

Results for “eigenmode analysis”

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 19 records

Dynamic mode decomposition for gyrokinetic eigenmode analysis

Dynamic mode decomposition (DMD) is a post-processing approach to decompose a complex time series into a set of modes via spectral analysis. DMD provides a new and powerful method to recover gyrokinetic drift-wave eigenfrequencies and eigenfunctions based only on the solution of the gyrokinetic-Maxwell initial value problem with almost no added cost to the initial value solver. In the present paper, DMD is applied to the CGYRO gyrokinetic code using a newly-developed CGYRO-DMD post-processor. CGYRO-DMD is numerically efficient, even on a single CPU. It does not set any restrictions on the plasma shape, beta (ratio of the plasma pressure to the magnetic field pressure), collisionality or number of species, and allows one to resolve numerous eigenmodes, even of comparable growth rates. In addition, DMD is not limited to unstable modes, but rather can capture stable and unstable branches simultaneously. In this work, we illustrate the accuracy of DMD through gyrokinetic analysis of mode transition for electromagnetic drift wave instabilities.

drift-wave eigenmodes

Eigenmode analysis of unsteady one-dimensional Euler equations

The initial boundary value problem describing the evolution of unsteady linearized perturbations of a steady, uniform subsonic flow is analyzed. The eigenmodes and eigenfrequencies of the system are derived and several examples are presented to illustrate the effect of different boundary conditions on the exponential decay rate of the eigenmodes. The resultant implications for the stability and convergence rates of finite difference computations are discussed.

Giles, M.

Eigenmode Analysis of Boundary Conditions for One-Dimensional Preconditioned Euler Equations

An analysis of the effect of local preconditioning on boundary conditions for the subsonic, one-dimensional Euler equations is presented. Decay rates for the eigenmodes of the initial boundary value problem are determined for different boundary conditions. Riemann invariant boundary conditions based on the unpreconditioned Euler equations are shown to be reflective with preconditioning, and, at low Mach numbers, disturbances do not decay. Other boundary conditions are investigated which are non-reflective with preconditioning and numerical results are presented confirming the analysis.

Darmofal, David L.

Radiative and continuum dampings of reversed shear Alfvén eigenmodes and perturbative analysis limitations for tokamaks

A careful theoretical analysis of the excitation of Alfvén eigenmodes (AEs), such as TAE (toroidicity-induced AE) and RSAE (reversed shear AE), by superalfvenic energetic particles is required for reliable predictions of energetic ion relaxation in present day fusion experiments. This includes the evaluation of different AE damping mechanisms including radiative and continuum dampings which are the focus of this study. A recent comprehensive benchmark of different eigenmode solvers including gyrokinetic, gyrofluid and hybrid magenetohydrodynamics (MHD) has shown that employed models may have deficiencies when addressing some of them (Taimourzadeh et al., Nucl. Fusion, vol. 59, 2019, 066006). Here, in this paper, we are studying the radiative and continuum dampings of RSAEs in details which were missing in hybrid NOVA/NOVA-C calculations to prepare a NOVA-C package with a substantial upgrade. Both dampings require the finite Larmor radius (FLR) corrections to AE mode structures to be accounted for. Accurately calculating different damping rates and understanding their parametric dependencies, we resolve the limitation coming out of the perturbative approach. In particular, here, the radiative damping is included perturbatively, whereas the continuum damping is computed non-perturbatively. Our comparison leads to the conclusion that the non-perturbative treatment of the unstable RSAE modes is needed to find the agreement with the gyrokinetic calculations. We expect that the RSAE mode structure modification plays a dominant role in determining the RSAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Impedance and wakefield studies of the EIC RCS 591 MHz five-cell cavity

The Electron-Ion Collider (EIC) is a next-generation accelerator complex designed to enable high-luminosity collisions between highly polarized electrons and light ions (e.g., He-3). A central component of its Electron Injection System (EIS) is the Rapid Cycling Synchrotron (RCS), which accelerates a single 28 nC electron bunch from 750 MeV to 5, 10, or 18 GeV using an array of 591 MHz five-cell superconducting RF (SRF) cavities—eight at the current design stage. To ensure stable acceleration of high-charge bunches, we conducted detailed impedance and wakefield studies of the SRF cavity structure using both frequency- and time-domain methods. Wakefield solvers (ECHO3D, ECHO1D, CST), eigenmode analysis, and multi-particle tracking with ELEGANT were employed to evaluate longitudinal and transverse impedance effects and to determine instability thresholds. These studies provide critical input for the cavity design and operating parameters required to preserve beam quality and stability in the RCS.

Accelerator Physics

Hall current effect on tearing mode instability

From a linear 2-D eigenmode analysis, it is found that the Hall current effect on collisional tearing mode instability becomes important for the thin magnetic reversal layer whose width is comparable to the ion inertia length; Hall currents produce a three-dimensional field structure and increase the reconnection (growth) rate. Since the magnetaic reversal layer widths both in the magnetopause and in the magnetotail are reported to become as thin as the ion inertial length (several hundred km) when the reconnection process is supposed to occur, the Hall current effect may explain the appearance of the dawn-dusk component of the magnetic field in the magnetotail reconnection region.

Terasawa, T.

Computation Of Flutter In Turbomachinery

"Direct" solution procedure supplants conventional eigenmode analysis. Iterative computational procedure yields critical mach number for onset of flutter and flutter angular frequency, of propfan, turbine with unshrouded blades, or other turbomachinery. Procedure applied to aeroelastic analysis of propfan, in which finite-element model of propfan structure combined with model of unsteady aerodynamics based on three-dimensional subsonic-lifting-surface theory. Particularly suitable for optimization in design because as optimal design evolves, flutter solution expected to change incrementally, so previous solution provides good estimates for start of current solution.

Murthy, Durbha V.

Linear theory of fast reconnection at an X-type neutral point

A complete analytic theory of dynamic linear reconnection at an X-type neutral point is developed. An eigenmode analysis, using cylindrical coordinates centered on the neutral point, extends the work of Craig & McClymont (1991) to include nonazimuthally symmetric perturbations on the flux function. It is demonstrated that all physically significant disturbances, both reconnective and nonreconnective, decay resistively on a 'fast' time scale. The significance of the linear theory is discussed within the general context of steady state and dynamic reconnection studies. Disturbances can be expected to focus explosively in the vicinity of the neutral point. This suggests the formation of a 'flux pile-up' current layer in which the bulk of the magnetic energy is released as heat rather than kinetic energy of mass motion.

Craig, I. J. D.

Gradient drift eigenmodes in the equatorial electrojet

The problem of kilometer-scale irregularities in the daytime equatorial electrojet is revisited by means of an eigenmode analysis of the gradient drift instability. Realistic physical parameters are used, including the modeled altitude variations of ion and electron collision frequencies and mobilities. The full fourth-order system of two coupled differential equations (each of second order) for the denisty and electrostatic potential perturbations is solved numerically by a relaxation technique. Under some approximations, the fourth-order system can be shown to reduce to a second-order differential equation for the perturbed potential or density. The latter is solved using a shooting technique and provides initial guesses for numerical solutions to the full problem. It is shown that the linear growth rate peaks for kilometer-scale waves, contrary to the findings of recent initial-value studies. This occurs because the equilibrium velocity shear is much more effective as a damping mechanism for short-wavelength modes than it is for the longer, kilometer-scale modes. These results provide a natural qualitative explanation for the observed dominance of kilometer-scale structures in the daytime electrojet spectrum.

Wang, X.-H.

Mode Profiles in Waveguide-Coupled Resonators

Surface acoustic wave (SAW) waveguide-coupled resonators are of considerable interest for narrow-band filter applications, though to date there has been very little published on the acoustic details of their operation. As in any resonator, one must fully understand its mode structure and herein we study the SAW mode profiles in these devices. Transverse mode profiles in the resonant cavity of the device were measured at various frequencies of interest using a knife-edge laser probe. In addition we predict the mode profiles for the device structure by two independent methods. One is a stack-matrix approach adapted from integrated optics and the other is a conventional analytical eigenmode analysis of the Helmholtz equation. Both modeling techniques are in good agreement with the measured results.

Hunt, William D.

Stability analysis of Alfvén eigenmodes excited by ion cyclotron resonance heating on EAST

Alfvén wave instabilities driven by energetic particles are common in fusion devices. Understanding their behavior is essential for the good confinement of fusion plasma. For this purpose, a series of experiments are performed on EAST to investigate the excitation of toroidal Alfvén eigenmodes (TAEs). Experimentally, it was found that AEs with frequencies around 80, 134, 157 kHz are excited by ion cyclotron resonance heating (ICRH) hydrogen minority heating scenario. Moreover, the excitation of AEs is independent of the wall-coating materials, but determined with the generation and confinement of the fast ions. Statistical analysis suggests that the TAEs can only be excited when the ICRH power is larger than 2 MW and the H 98 factor is larger than 1.15. In line with the experiments, a set of simulations using TORIC, ASCOT, and NOVA-C is performed. The simulation results show that with the experimental density profile and safety factor, the measured TAE with frequency of 134 kHz is well reproduced by taking into account the plasma toroidal rotation frequency. However, the appearance of the 80 kHz mode cannot be captured in the simulations. The ICRH-generated fast hydrogen ions mainly have banana orbits and that their perpendicular velocity is much larger than that of the parallel ions. The stability analysis with NOVA-C suggests that the radial gradient of the fast H ion distribution is the only driving source of TAEs. For TAEs locate in the center of Alfvén continuum, the D ion Landau damping is the dominant damping term; For most other TAEs which have intersections with the continuum boundary, the continuum damping becomes the dominant damping mechanism.

Alfvén eigenmodes

Dynamic analysis of evolutive conservative systems. Discussion of eigenmode crossings

After an analysis of the close connection between the symmetries of a dynamical system and the multiplicity of its vibrational natural frequencies, it is proved by variational arguments that for a system of invariable symmetry the eigenfrequencies associated with the eigenmodes of a given symmetry type do not cross, in general, during the evolution of this system. The theory is implemented by some numerical calculations applied to the analysis of the evolution of the axisymmetric hydroelastic modes of the Ariane launch vehicle during burning of the first stage.

Morand, H. J. P.

Accuracy of an approximate static structural analysis technique based on stiffness matrix eigenmodes

Use of the stiffness matrix eigenmodes, instead of the vibration eigenmodes, as generalized coordinates is proposed for condensation of static load deflection equations in finite element stiffness method. The modes are selected by strain energy criteria and the resulting fast, approximate analysis technique is evaluated by applications to idealized built-up wings and a fuselage segment. The best results obtained are a two-order of magnitude reduction of the number of degrees of freedom in a high aspect ratio wing associated with less than one percent error in prediction of the largest displacement.

Sobieszczanski-Sobieski, J.

Analysis of the linear and nonlinear stability of Alfven eigenmodes and fish-bones in JET DT discharges: mode identification and shear flows generation

The plasma in future nuclear fusion reactors will be heated by neutral beam injectors (NBIs) and high frequency electromagnetic waves as well as fusion born alpha particles. Energetic particles (EPs), with energies up to two orders of magnitude larger than the thermal plasma, can trigger EP driven modes and induce harmful EP losses, reducing the plasma heating efficiency and the economical viability of the reactor. The present study is dedicated to analyze the Alfven Eigenmode (AE) activity in JET D–T discharges, the closest experiment to reactor-like operation performed until now. There, EP driven modes are induced by the combined effect of tangential NBIs and ion cyclotron resonance heating (ICRH) driven EP. Linear and nonlinear simulations are performed with the gyro-fluid FAR3d code to analyze the AE activity observed in the discharge 99896. The linear simulations reproduce the unstable n = 3 to 5 toroidal AEs (TAE) at the inner plasma region observed in the experiment, triggered by highly energetic passing deuterium populations injected by the tangential NBIs, further accelerated by the effect of the ICRH up to 1 MeV. In addition, fish-bones triggered by energetic trapped hydrogen induced by the ICRH are also reproduced. On the other hand, the alpha particles density is too small to destabilize AEs in the experiment. Nonetheless, increasing artificially the alpha density by one order of magnitude, an n = 1 beta induced AE can be destabilized in the inner plasma region. Nonlinear simulations indicate the generation of zonal structures during the AE/fish-bone saturation phase. TAE and fish-bones causes a rather weak increase of the passing D and trapped H EP (around 2%), respectively. Shear flows and zonal currents are generated during the saturation of TAE and fish-bones. Nonlinear simulations performed for D–T and pure deuterium thermal plasma indicate AE/fish-bone activity is weaker and shear flows are less intense in the pure deuterium case, trends consistent with the experimental observations that also indicates a deterioration of the thermal plasma confinement. Therefore, both numerical studies and experimental evidence indicate the generation of shear flows by AE/fish-bones could be connected with an improvement of the thermal plasma confinement.

AE

An investigation of the relationship between sub-Saharan rainfall and global sea surface temperatures

The relationship between SST and rainfall index anomalies over sub-Saharan Africa for the 1970-1984 period is investigated. Results of an empirical orthogonal function analysis indicate that the most dominant eigenmode, EOF1, is characterized by warming over the central eastern Pacific, cooling over the eastern midlatitude Pacific, and warming over the entire Atlantic and Indian ocean basins. EOF1 is found to have statistically signifiant monthly correlations for the Sahel and Soudan regions, with the warm El Nino-like phases of SST EOF1 corresponding to drought conditions. These results suggest that the large-scale SST anomalies may be responsible for a large component of the observed vacillation of sub-Saharan rainfall.

Semazzi, F. H. M.

Near-symmetric instability for general Prandtl number

An analysis of the stability of a simple baroclinic flow to perturbations whose horizontal wave vector lies in and near the vertical plane containing the density gradient has been performed. The Ekman number (E) and 'azimuthal' wavenumber (alpha) were both assumed much less than 1, and expansions about these quantities were performed for the eigenmodes and Richardson number (R). Agreement with the previous analysis of Busse and Chen (1981) that the correction to the critical R was O(alpha) for Prandtl number (P) away from unity was obtained. An expression for the correction to the critical R for arbitrary P, and an approximate expression for P = 1 were obtained as new results. The correction for P = 1 has the same sign as that for P less than 1, and is O(alpha E exp 2/3). This result compares well with the numerical results of Miller and Antar (1986) for small E.

Reynolds, Nathaniel D.

Nonlocal analysis of finite-beam-driven instabilities

The fully kinetic integral eigenmode equation in wave-number space is used to describe the nonlocal behavior of electrostatic waves in an electron-beam plasma, which are studied in the low-temperature-beam regime and the warm-beam regime. The case of strongly magnetized electrons and unmagnetized ions, which corresponds to the waves in a frequency range from the lower-hybrid to the electron plasma frequency, is examined. Three wave modes are found. The first group consists of modes that have dispersive properties similar to the uniform, infinite beam-plasma system. Depending on the beam width, the growth rates are strongly reduced. The second group, surface modes, are localized at the periphery of the beam region and are less unstable than the unstable modes of the first group. The third group represents natural oscillations of the background plasma. These modes are virtually unaffected by the beam.

Serizawa, Y.