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 1,045 records · Page 58
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Kinetic stability of Chapman–Enskog plasmas

In this paper, we investigate the kinetic stability of classical, collisional plasma – that is, plasma in which the mean-free-path $\lambda$ of constituent particles is short compared with the length scale $L$ over which fields and bulk motions in the plasma vary macroscopically, and the collision time is short compared with the evolution time. Fluid equations are typically used to describe such plasmas, since their distribution functions are close to being Maxwellian. The small deviations from the Maxwellian distribution are calculated via the Chapman–Enskog (CE) expansion in $\lambda /L \ll 1$ , and determine macroscopic momentum and heat fluxes in the plasma. Such a calculation is only valid if the underlying CE distribution function is stable at collisionless length scales and/or time scales. We find that at sufficiently high plasma $\beta$ , the CE distribution function can be subject to numerous microinstabilities across a wide range of scales. For a particular form of the CE distribution function arising in strongly magnetised plasma ( viz. plasma in which the Larmor periods of particles are much smaller than collision times), we provide a detailed analytic characterisation of all significant microinstabilities, including peak growth rates and their associated wavenumbers. Of specific note is the discovery of several new microinstabilities, including one at sub-electron-Larmor scales (the ‘whisper instability’) whose growth rate in certain parameter regimes is large compared with other instabilities. Our approach enables us to construct the kinetic stability maps of classical, two-species collisional plasma in terms of $\lambda$ , the electron inertial scale $d_e$ and the plasma $\beta$ . This work is of general consequence in emphasising the fact that high- $\beta$ collisional plasmas can be kinetically unstable; for strongly magnetised CE plasmas, the condition for instability is $\beta \gtrsim L/\lambda$ . In this situation, the determination of transport coefficients via the standard CE approach is not valid.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Diffusion at the magnetopause: Hybrid simulations

Electromagnetic wave generation and resulting cross-field diffusion of plasma are considered at a tangential discontinuity, which characterizes the magnetopause for northward interplanetary magnetic field. Two-dimensional hybrid (particle ions, massless fluid electrons) simulations, in which the tangential discontinuity is generated self-consistently via a stream-stream interaction, are used to show that wave growth occurs when the ambient magnetic field is predominantly perpendicular to the direction of the density gradient. Low-frequency (much less than ion gyrofrequency) waves, with amplitudes delta B/B less than or equal to 0.2 and anticorrelated density fluctuations delta n/n less than or equal to 0.6, are generated at the discontinuity, resulting in cross-field diffusion that is comparable to the Bohm rate. Both the fluctuation level and the lack of ion heating in the calculations are consistent with observations at the magnetopause. The magnitude of the diffusion is considered in the presence of numerical effects and in the context of the inferred diffusion rate at the magnetopause. The relation of the low-frequency waves and their consequences to faster growing, short-wavelength waves due to the lower hybrid drift instability is also addressed. The overall conclusion of this initial study is that diffusion due to low frequency waves is not likely to be a major effect at the magnetopause.

Winske, D.↗
NASA NTRS2002

Lattice Boltzmann Method for 3-D Flows with Curved Boundary

In this work, we investigate two issues that are important to computational efficiency and reliability in fluid dynamics applications of the lattice, Boltzmann equation (LBE): (1) Computational stability and accuracy of different lattice Boltzmann models and (2) the treatment of the boundary conditions on curved solid boundaries and their 3-D implementations. Three athermal 3-D LBE models (D3QI5, D3Ql9, and D3Q27) are studied and compared in terms of efficiency, accuracy, and robustness. The boundary treatment recently developed by Filippova and Hanel and Met et al. in 2-D is extended to and implemented for 3-D. The convergence, stability, and computational efficiency of the 3-D LBE models with the boundary treatment for curved boundaries were tested in simulations of four 3-D flows: (1) Fully developed flows in a square duct, (2) flow in a 3-D lid-driven cavity, (3) fully developed flows in a circular pipe, and (4) a uniform flow over a sphere. We found that while the fifteen-velocity 3-D (D3Ql5) model is more prone to numerical instability and the D3Q27 is more computationally intensive, the 63Q19 model provides a balance between computational reliability and efficiency. Through numerical simulations, we demonstrated that the boundary treatment for 3-D arbitrary curved geometry has second-order accuracy and possesses satisfactory stability characteristics.

Mei, Renwei↗
NASA NTRS2021

Pretest Computational Assessment of Boundary Layer Transition in the NASA Juncture Flow Model with an NACA 0015-Based Wing

The first two phases of the NASA Juncture Flow experiment were carried out on a DLR-F6 swept-wing model and were designed to provide “CFD validation-quality” data toward the assessment and improvement of existing CFD turbulence models in predicting onset and extent of three-dimensional separated flow near the wing-juncture trailing-edge region. The next phase of experiments will involve an NACA 0015-based swept wing, as prior risk reduction experiments had indicated that this wing shape resulted in reduced separation near the juncture region than the DLR-F6 wing, thus providing a better option to evaluate the ability of CFD models to predict incipient turbulent separation. The NACA 0015 measurements will also include IR thermography to infer the variation of transition front with respect to an increasing angle of attack. The primary objective of this work is to computationally make a preliminary assessment of the transition front on both surfaces of the NACA 0015 wing at a crank-chord-based Reynolds number of 2.4 x 106 for four different angles of attack, (0°, 2.5°, 5°, and 7.5°) and to determine the dominant mechanisms responsible for transition. This assessment includes both RANS-based transition models from NASA’s OVERFLOW 2.3b flow solver and linear parabolized stability equations (PSE) stability analysis based on the Langley Stability and Transition Analysis code, LASTRAC. Linear PSE results indicate that the upper surface of the wing is dominated by Tollmien- Schlichting (TS) instabilities, and that the laminar flow region shrinks from about 50% chord to a very small region just downstream of the attachment line as the angle of attack is increased from 0° to 7.5°. Consequently, the transition fronts predicted by the Spalart- Allmaras-based amplification factor transport (AFT-2017b) equation model (which accounts for the TS instabilities alone) and the Menter’s shear-stress transport equation (SST2003)- based Langtry-Menter transition model with ability to account for both TS and crossflow effects (LM2015) compare well with those predicted using linear PSE. On the lower surface of the wing, stationary crossflow (CF) instabilities begin to appear on the inboard portion of the wing in addition to the TS-instabilities for the larger angles of attack (5° and 7.5°), further reducing the laminar flow extent within the inboard region. The LM2015 model that accounts for CF effects is able to replicate this trend but appears to predict a slightly earlier transition. The outcome of this effort will inform the experiment and, when the actual experimental data become available, provide further opportunity to assess and improve the various transition models.

CFD modeling↗
DOE OSTI2022

Sudden collapse of a pressure profile generated by off-axis heating in a linear magnetized plasma

Here, the features of an unexpected, large event that arises spontaneously during a basic heat transport experiment are presented. It consists of the sudden collapse of the radial plasma pressure profile, akin to disruption events observed in toroidal magnetic confinement devices. The experiment is performed on the Large Plasma Device at the University of California, Los Angeles (UCLA). It uses a LaB6 thermionic emitter of annular shape to induce off-axis heating of a cold, afterglow plasma, in a linear magnetic geometry. The temporal evolution consists of three regimes. During an early, quiescent period, classical heat transport along and across the magnetic field arises from Coulomb collisions. After significant pressure gradients develop, drift-Alfvén waves become unstable. Upon reaching large amplitude, they trigger avalanche events that flatten the outer part of the heated region, which, in turn, quenches the instability. Due to the sustained heating, the pressure profile rebuilds and the process repeats, leading to a relatively long, second regime that displays multiple avalanches, but suddenly, the annular pressure profile is observed to collapse. After this collapse, the system enters a third regime with large fluctuations. Before the collapse, a rapid, runaway heating environment arises whose time evolution exhibits a self-similar dependence on the applied voltage. The time evolution, morphology, and scaling of the collapse event are presented, and an examination is made of the underlying mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2020

Generation and Calibration of Linear Models of Aircraft with Highly Coupled Aeroelastic and Flight Dynamics

The lightweight structures and unconventional configurations being considered for the next generation of aircraft mean that any effort to predict or control the flight dynamics is impacted by the structural dynamics. One of the most severe forms of coupling between aeroelasticity and flight dynamics is an instability called body freedom flutter. The existing tools often assume a relatively weak effect of structural dynamics on the flight dynamics, and are therefore incapable of modeling strong interactions like body freedom flutter. A method of combining different sources of data traditionally used for aeroelasticity and flight dynamics is described by reconciling many of the differences between these models. By building upon past modeling efforts, a level of familiarity in the approach is achieved. Generally the differences from the traditional approaches are subtle but significant. The traditional frequency domain flutter model in a modal coordinate system is converted to a form consistent with a time domain flight dynamics model. The time domain rational function approximation about a non-inertial coordinate system and the unique constraints for the conversion between the inertial and non-inertial coordinate systems are discussed. A consistent transformation of the states of aeroelastic models to flight dynamics models is derived, which enables the integration of data from higher fidelity computational fluid dynamics models or wind-tunnel testing. The present method of integrating multidisciplinary data was used to create models that compare well with X-56A flight-test data, including conditions past the flutter speed.

Ouellette, Jeffrey↗
NASA NTRS1991

Turbulization of shear flows in astrophysics

The problem of instability in shear flows is treated in a number of works which present a global analysis of irrotational (potential) modes in a 2-D compressible shear flow, implying that, in the presence of reflecting boundaries, those modes can increase and that the characteristic time of the increase is several orders as high as the dynamic time in the medium. Due to the slowness of the increase of irrotational modes and to the problematic character of the theory's applicability for thin accretion disks, an attempt was made to develop an alternative theory of turbulization shear flows. It is quite contrary to the above theory. Instead of the global analysis, a local one is made far from the flow boundaries. Instead of irrotational perturbations, vortical ones are examined. Lastly, the medium is considered to be incompressible. The flow in the accretion disk can be modeled by a plane Couette flow, without regard for boundary conditions: it can be modeled by a plane shear flow in infinite incompressible fluid. According to this scenario, there is a critical perturbation level in the free shear flow. If the level of initial perturbations exceeds the critical one, the flow is turbulized.

Chagelishvili, G. D.↗
NASA NTRS1993

A time-accurate finite volume method valid at all flow velocities

A finite volume method to solve the Navier-Stokes equations at all flow velocities (e.g., incompressible, subsonic, transonic, supersonic and hypersonic flows) is presented. The numerical method is based on a finite volume method that incorporates a pressure-staggered mesh and an incremental pressure equation for the conservation of mass. Comparison of three generally accepted time-advancing schemes, i.e., Simplified Marker-and-Cell (SMAC), Pressure-Implicit-Splitting of Operators (PISO), and Iterative-Time-Advancing (ITA) scheme, are made by solving a lid-driven polar cavity flow and self-sustained oscillatory flows over circular and square cylinders. Calculated results show that the ITA is the most stable numerically and yields the most accurate results. The SMAC is the most efficient computationally and is as stable as the ITA. It is shown that the PISO is the most weakly convergent and it exhibits an undesirable strong dependence on the time-step size. The degenerated numerical results obtained using the PISO are attributed to its second corrector step that cause the numerical results to deviate further from a divergence free velocity field. The accurate numerical results obtained using the ITA is attributed to its capability to resolve the nonlinearity of the Navier-Stokes equations. The present numerical method that incorporates the ITA is used to solve an unsteady transitional flow over an oscillating airfoil and a chemically reacting flow of hydrogen in a vitiated supersonic airstream. The turbulence fields in these flow cases are described using multiple-time-scale turbulence equations. For the unsteady transitional over an oscillating airfoil, the fluid flow is described using ensemble-averaged Navier-Stokes equations defined on the Lagrangian-Eulerian coordinates. It is shown that the numerical method successfully predicts the large dynamic stall vortex (DSV) and the trailing edge vortex (TEV) that are periodically generated by the oscillating airfoil. The calculated streaklines are in very good comparison with the experimentally obtained smoke picture. The calculated turbulent viscosity contours show that the transition from laminar to turbulent state and the relaminarization occur widely in space as well as in time. The ensemble-averaged velocity profiles are also in good agreement with the measured data and the good comparison indicates that the numerical method as well as the multipletime-scale turbulence equations successfully predict the unsteady transitional turbulence field. The chemical reactions for the hydrogen in the vitiated supersonic airstream are described using 9 chemical species and 48 reaction-steps. Consider that a fast chemistry can not be used to describe the fine details (such as the instability) of chemically reacting flows while a reduced chemical kinetics can not be used confidently due to the uncertainty contained in the reaction mechanisms. However, the use of a detailed finite rate chemistry may make it difficult to obtain a fully converged solution due to the coupling between the large number of flow, turbulence, and chemical equations. The numerical results obtained in the present study are in good agreement with the measured data. The good comparison is attributed to the numerical method that can yield strongly converged results for the reacting flow and to the use of the multiple-time-scale turbulence equations that can accurately describe the mixing of the fuel and the oxidant.

Kim, S.-W.↗
NASA NTRSDate not supplied

X-57 “Maxwell” High-Lift Propeller Testing and Model Development

NASA’s X-57 “Maxwell” distributed electric propulsion flight demonstrator has a high-lift system that includes 12 fixed-pitch high-lift propellers located upstream of the wing leading edge for lift augmentation at low speeds. These high-lift propellers are not required at higher speeds and are folded conformally along the nacelles to reduce drag when not in operation. Aircraft performance models and flight simulations incorporate propeller performance to predict thrust, moments, and power consumption and propeller model accuracy is important in identifying safe operating regimes for this aircraft. Current high-lift propeller performance models have been verified and calibrated against numerous computational fluid dynamics analyses under a variety of flight conditions. We performed a series of full-scale wind tunnel tests at the NASA Langley Research Center Low Speed Aeroacoustic Wind Tunnel to further validate these models, to identify any adverse operating conditions for these propellers, and to assess the accuracy of the facility’s new Propeller Test Stand. The results indicate that the models accurately predicted performance and that the right-handed propeller showed lower torque and thrust for a given propeller speed compared to the left-handed version. Acoustic data support these measurements, showing higher tonal noise for the left versus right propeller. We suspect this is due to slight differences in the material properties between the early production (right) and more recent (left) blade sets leading to greater detwisting in the right-handed blades under load. Both propellers demonstrated very stable operation throughout the test including during deployment, stowing, and windmilling. Furthermore, high-speed imagery of the spinning propeller indicated no blade oscillations or instabilities for any test case. We also leveraged these images to analyze the effects of propeller loading on blade detwisting and tip position. Unfortunately, the thrust measurements varied significantly during the test due to cross-loading of the Propeller Test Stand load cell, though trends in the thrust data were still observable. The torque measurement, however, remained very precise for all test points and was used for performance comparisons. As a result of this test, we validated and improved the X-57 high-lift propeller performance models, obtained detailed propeller acoustic data and high-speed imagery, observed the propeller operating safely during a variety of dynamic events, and identified two potentially adverse, low-propeller-speed conditions which will be avoided during normal aircraft operation.

X-57↗
NASA NTRS1993

LIF measurements of scalar mixing in turbulent shear layers

The structure of shear layer flows at high Reynolds numbers remains a very interesting problem. Straight mixing layers have been studied and yielded information on the probability density function (pdf) of a passive scalar across the layer. Konrad and Koochesfahani & Dimotakis measured the pdf of the mixture fraction for mixing layers of moderate Reynolds numbers, each about 25,000 (Re based on velocity difference and visual thickness). Their measurements showed a 'non-marching' pdf (central hump which is invariant from edge to edge across the layer), a result which is linked to the visualizations of the spanwise Kelvin-Helmholtz (K-H) instability mode, which is the primary instability for plane shear layer flows. A secondary instability mode, the Taylor-Gortler (T-G) instability, which is associated with streamwise vortical structures, has also been observed in shear layers. Image reconstruction by Jimenez et al. and volume renderings by Karasso & Mungal at low Re numbers have demonstrated that the K-H and the T-G instability modes occur simultaneously in a non-mutually destructive way, evidence that supports the quasi two-dimensional aspect of these flows and the non-marching character of the pdf at low Reynolds numbers. At higher Re numbers though, the interaction of these two instability modes is still unclear and may affect the mixing process. In this study, we perform measurements of the concentration pdf of plane mixing layers for different operating conditions. At a speed ratio of r = U(sub 1)/U(sub 2) = 4:1, we examine three Reynolds number cases: Re = 14,000, Re = 31,000, and Re = 62,000. Some other Re number cases' results, not presented in detail, are invoked to explain the behavior of the pdf of the concentration field. A case of r = 2.6:1 at Re = 20,000 is also considered. The planar laser-induced fluorescence technique is used to yield quantitative measurements. The different Re are obtained by changing the velocity magnitudes of the two streams. The question of resolution of these measurements is addressed. In order to investigate the effects of the initial conditions on the development and the structure of the mixing layer, the boundary layer on the high-speed side of the splitter plate is tripped. The average concentration and the average mixed fluid concentration are also calculated to further understand the changes in the shear layer for the different cases examined.

Karasso, Paris S.↗
DOE OSTI2020

Characterizing magnetized plasmas with dynamic mode decomposition

Accurate and efficient plasma models are essential to understand and control experimental devices. Existing magnetohydrodynamic or kinetic models are nonlinear and computationally intensive and can be difficult to interpret, while often only approximating the true dynamics. In this work, data-driven techniques recently developed in the field of fluid dynamics are leveraged to develop interpretable reduced-order models of plasmas that strike a balance between accuracy and efficiency. In particular, dynamic mode decomposition (DMD) is used to extract spatio-temporal magnetic coherent structures from the experimental and simulation datasets of the helicity injected torus with steady inductive (HIT-SI) experiment. Three-dimensional magnetic surface probes from the HIT-SI experiment are analyzed, along with companion simulations with synthetic internal magnetic probes. A number of leading variants of the DMD algorithm are compared, including the sparsity-promoting and optimized DMD. Optimized DMD results in the highest overall prediction accuracy, while sparsity-promoting DMD yields physically interpretable models that avoid overfitting. These DMD algorithms uncover several coherent magnetic modes that provide new physical insights into the inner plasma structure. These modes were subsequently used to discover a previously unobserved three-dimensional structure in the simulation, rotating at the second injector harmonic. Finally, using data from probes at experimentally accessible locations, DMD identifies a resistive kink mode, a ubiquitous instability seen in magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

SRM Internal Flow Tests and Computational Fluid Dynamic Analysis: RSRM Full Scale Motor Analyses - Volume 2

The purpose of the RSRM Nozzle Slag Ejection Precursor Test is to investigate the effect that slag ejection from the RSRM nozzle has on the chamber pressure and trust of the SRB's. In past firings of the Reusable Solid Rocket Motor (RSRM) both static test and flight motors have shown small pressure perturbations occurring primarily between 65 and 80 seconds. A joint NASA/Thiokol team investigation concluded that the cause of the pressure perturbations was the periodic ingestion and ejection of molten aluminum oxide slag from the cavity around the submerged nozzle nose which tends to trap and collect individual aluminum oxide droplets from the approach flow. The conclusions of the team were supported by numerous data and observations from special tests including high speed photographic films, real time radiography, plume calorimeters, accelerometers, strain gauges, nozzle TVC system force gauges, and motor pressure and thrust data. A simplistic slag ballistics model was formulated to relate a given pressure perturbation to a required slag quantity. Also, a cold flow model using air and water was developed to provide data on the relationship between the slag flow rate and the chamber pressure increase. Both the motor and the cold flow model exhibited low frequency oscillations in conjunction with periods of slag ejection. Motor and model frequencies were related to scaling parameters. The data indicate that there is a periodicity to the slag entrainment and ejection phenomena which is possibly related to organized oscillations from instabilities in the dividing streamline shear layer which impinges on the underneath surface of the nozzle.

Source record↗
NASA NTRS1989

Time-dependent computational studies of flames in microgravity

The research performed at the Center for Reactive Flow and Dynamical Systems in the Laboratory for Computational Physics and Fluid Dynamics, at the Naval Research Laboratory, in support of the NASA Microgravity Science and Applications Program is described. The primary focus was on investigating fundamental questions concerning the propagation and extinction of premixed flames in Earth gravity and in microgravity environments. The approach was to use detailed time-dependent, multispecies, numerical models as tools to simulate flames in different gravity environments. The models include a detailed chemical kinetics mechanism consisting of elementary reactions among the eight reactive species involved in hydrogen combustion, coupled to algorithms for convection, thermal conduction, viscosity, molecular and thermal diffusion, and external forces. The external force, gravity, can be put in any direction relative to flame propagation and can have a range of values. A combination of one-dimensional and two-dimensional simulations was used to investigate the effects of curvature and dilution on ignition and propagation of flames, to help resolve fundamental questions on the existence of flammability limits when there are no external losses or buoyancy forces in the system, to understand the mechanism leading to cellular instability, and to study the effects of gravity on the transition to cellular structure. A flame in a microgravity environment can be extinguished without external losses, and the mechanism leading to cellular structure is not preferential diffusion but a thermo-diffusive instability. The simulations have also lead to a better understanding of the interactions between buoyancy forces and the processes leading to thermo-diffusive instability.

Oran, Elaine S.↗
NASA NTRS2005

The Dynamics of Miscible Interfaces: Simulations

This research project focuses on the dynamics of interfacial regions between miscible fluids. While much attention has focused on immiscible interfaces in the past, miscible interfaces have been explored to a much lesser degree, so that there are many open questions regarding their dynamics at this time. Among the more pressing issues is the role that nonconventional stresses can play in such interfacial regions. Such stresses are typically not accounted for in efforts to model the dynamics of miscible flows. Our research aims to clarify under which circumstances these stresses do have to be taken into account, and what quantitative approaches are most suitable in this regard. In order to address these issues, we have focused on conducting linear stability analyses and nonlinear simulations for capillary tube and Hele-Shaw flows, and to compare the results with corresponding experiments performed in the labs of our co-investigators Prof. Maxworthy at USC, and Dr. Balasubramaniam at NASA. Over the duration of the project we have, among other things, focused on the effects of variable diffusion coefficients in such flows, and specifically on their influence in the growth of instabilities. Furthermore, our three-dimensional spectral element simulations have made good progress, so that we have come to a point where we can conduct more detailed comparisons with experimental observations. We are currently focusing our efforts on reproducing the tip-splitting instability observed by Maxworthy. Finally, we have discovered a new core-annular flow instability in the Stokes flow regime during the last year. This represents a significant finding, as this instability does not have an immiscible counterpart.

Meiburg, Eckart↗
DOE OSTI2023

Staggered scheme for the compressible fluctuating hydrodynamics of multispecies fluid mixtures

Here, we present a numerical formulation for the solution of nonisothermal, compressible Navier-Stokes equations with thermal fluctuations to describe mesoscale transport phenomena in multispecies fluid mixtures. The novelty of our numerical method is the use of staggered grid momenta along with a finite volume discretization of the thermodynamic variables to solve the resulting stochastic partial differential equations. The key advantages of the numerical scheme are that it significantly simplifies the discretization of diffusive and stochastic momentum fluxes into a more compact form, and it provides an unambiguous prescription of boundary conditions involving pressure. The staggered grid scheme more accurately reproduces the equilibrium static structure factor of hydrodynamic fluctuations in gas mixtures compared to a collocated scheme described previously by Balakrishnan et al. [Phys. Rev. E 89, 013017 (2014)1539-375510.1103/PhysRevE.89.013017]. The numerical method is tested for ideal noble gases mixtures under various nonequilibrium conditions, such as applied thermal and concentration gradients, to assess the role of cross-diffusion effects, such as Soret and Dufour, on the long-ranged correlations of hydrodynamic fluctuations, which are also more accurately reproduced compared to the collocated scheme. We numerically study giant nonequilibrium fluctuations driven by concentration gradients and fluctuation-driven Rayleigh-Taylor instability in gas mixtures. Wherever applicable, excellent agreement is observed with theory and measurements from the direct simulation Monte Carlo method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Optimizing $\mathrm{ACRT}$ to reduce inclusion formation during the $\mathrm{VGF}$ growth of cadmium zinc telluride: I. Computational approach

Theoretical computations reveal insight into parameters that govern the accelerated crucible rotation technique (ACRT) applied to the gradient freeze (VGF) growth of cadmium zinc telluride (CZT). A metric based on the classic Mullins and Sekerka instability criterion is put forth that provides a quantitative means of assessing the impact of different ACRT rotation schedules on morphological stability and inclusion formation during CZT growth. Here this metric is employed in a 2-k factorial design of experiments that identifies the maximum rotation rate as the most important ACRT parameter, followed by the acceleration period and the rest period. Subsequent optimization of rotation schedule is performed for a specific VGF growth system, and mechanistic explanations are provided for the preferred ACRT schedule, which improves interface stability via slower rotation rates and longer acceleration periods than would be applied using classical ACRT schedules.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Direct measurements of the 3D plasma velocity in single-helical-axis RFP plasmas

The first local velocity measurements of helical equilibrium plasmas in the Reversed Field Pinch (RFP) Single Helical Axis (SHAx) state using a Charge Exchange Recombination Spectroscopy (CHERS) diagnostic are presented. Measurements show strong axisymmetric and non-axisymmetric flow, with n ¼ 5 components of flow related to the ðm; nÞ¼ð1; 5Þ dominant magnetic mode on the order of the axisymmetric flow in certain regions of the plasma, as well as significant n > 5 flow. Flow measurements are compared with NIMROD simulations of visco-resistive, single-fluid MHD in toroidal and cylindrical geometries with limited axial periodicity. Both measurements and the simulation with toroidal geometry show stronger inboard flows relative to the outboard flows, which is attributed to the toroidal geometry of the device. In the experiment, the n ¼ 5 component of flow is phase shifted from the reconnection-like flow pattern observed in the single-fluid simulations, possibly due to decoupling of the ion and electron fluids over much of the plasma. Finally, the strength of the helical angular flow shear relative to the critical shear necessary to disrupt nonlinear coupling between tearing modes is calculated around the helical magnetic axis. The shear in the measured flow is on the order of the theoretical critical threshold needed to nonlinearly decouple modes, but the measurement uncertainty in the gradient of the flow is large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Excitation of toroidally localized harmonics of global Alfvén eigenmodes

Abstract The spherical tokamak NSTX and the upgraded, higher toroidal field, NSTX-U, routinely create plasmas with large un-thermalized populations of super-Alfvénic fast ions. This population of non-thermal energetic ions excites a broad spectrum of Alfvénic waves, from the lower frequency (≈30 kHz–≈200 kHz) toroidal Alfvén eigenmodes (Cheng and Chance 1986 Phys. Fluids 29 3695), to the higher frequency (400 kHz–3 MHz) global Alfvén eigenmodes (GAEs) and compressional Alfvén eigenmodes (Goedbloed 1975 Phys. Fluids 18 1258). In this paper we present evidence that the GAE non-linearly excite modes, presumably GAE, at frequencies consistent with non-linear or 3-wave coupling. The observation of the excitation of 2 nd harmonic GAE through the intrinsic non-linearity of Global Alfvén modes demonstrates that the non-linear terms can act as an exciter-antenna inside the plasma, broadcasting at harmonics of the mode frequency and with concomitantly shorter wavelengths. As with experiments using external antenna to excite otherwise weakly stable Alfvénic modes, this data can provide information about the stability of modes at harmonics of the GAE. It may also provide information on the nature of the non-linearities in the wave dispersion equation or a direct measure of the mode amplitude (Smith et al 2006 Phys. Plasmas 13 042504). We also report that the short wavelength 2 nd harmonic GAE can be strongly toroidally localized. The observation that the shorter wavelength ( n ≈ 20) 2 nd harmonic modes can be toroidally localized potentially has implications for the impact of fast-ion driven instabilities on fast-ion confinement in ITER and future fusion reactors (Gorelenkov et al 2014 Nucl. Fusion 54 125001).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗