Search NASA⌕ Search

SEARCH · Search NASA

Results for “STABILITY CALCULATIONS”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1952

Summary of methods for calculating dynamic lateral stability and response and for estimating aerodynamic stability derivatives

A summary of methods for making dynamic lateral stability and response calculations and for estimating the aerodynamic stability derivatives required for use in these calculations is presented. The processes of performing calculations of the time histories of lateral motions, of the period and damping of these motions, and of the lateral stability boundaries are presented as a series of simple straightforward steps. Existing methods for estimating the stability derivatives are summarized and, in some cases, simple new empirical formulas are presented. Detailed estimation methods are presented for low-subsonic-speed conditions but only a brief discussion and a list of references are given for transonic and supersonic speed conditions.

Campbell, John P↗
DOE OSTI2025

Low- n stability and plasma response to RMP in various STEP scenarios

The low-n (n is the toroidal mode number) magnetohydrodynamic (MHD) stability and plasma response are numerically investigated for various scenarios designed for STEP, that are relevant for the H-mode pedestal analysis. Control of the edge-localized modes (ELMs) with externally applied resonant magnetic perturbations (RMPs) is considered. Optimization of the ELM control coil current configuration, based on the computed plasma MHD response and well-established figures of merit validated on present-day experiments, finds reasonable robustness of a fixed coil phasing (for a given n-number) to control ELMs in all five STEP plasmas considered. Based on certain semi-empirical criteria, the required coil current to achieve ELM suppression is estimated to be about 10–20 kAt with the n = 1 or 2 RMP configuration and about 100–200 kAt for the n = 4 RMP. Systematic linear stability calculations are used to map out stability windows for the low-n kink-peeling modes, in terms of the ideal-wall location and variation of the edge safety factor q 95 with respect to the target design. The kink-peeling stability boundary is found to be generally sensitive to the q 95 variation, which has implications for achieving the quiescent H-mode regime in STEP. Full toroidal quasilinear initial-value simulations for these STEP plasmas find that generation of the edge-harmonic oscillations (EHOs) depends sensitively on the plasma scenario, the initial linear stability of the kink-peeling modes, the initial plasma toroidal flow and q 95 . In general, it is easier (more robust) to access the EHO-regime for two of the cases considered with smaller plasma volume and higher on-axis safety factor. Finally, quasilinear simulations find robust density pumpout due to applied RMPs in these STEP plasmas, but the effect on the plasma toroidal flow varies among different cases.

EHO↗
NASA NTRS1987

Stability/Instability Analysis of Rotating Machinery

Numerical index of stability calculated for nonlinear system. Technique for determining rotor stability or instability from analysis of measurements adapted for use with computer simulations of rotor motion. Involves calculation of log decrement or increment of vibration amplitude. Applicable to rotors mounted in loose bearings and to similar problems in which load-versus-deflection characteristics nonlinear. Developed for assessments of vibrational characteristics of turbopump rotors, technique also usable with such mechanisms as high-speed ball bearings.

Powers, Richard W.↗
DOE OSTI2026

Stability Constants of Lanthanide-nitrate Complexes in Aqueous Solutions: A Theoretical Study

Calculating stability constants for lanthanide--nitrate complexes in aqueous solution is challenging due to the complex free-energy landscapes of the participating species. In this work, we compare cluster-continuum solvation and condensed-phase approaches using universal machine learning interatomic potentials (MLIPs) for determining the stability constants of lanthanide--nitrate complexes in aqueous solutions. Within the cluster--continuum solvation framework at the B3LYP level of theory, reactions involving lanthanide coordination numbers of both 8 and 9 are found to be relevant. After an empirical linear free-energy correction, the cluster--continuum results fall on the same order-of-magnitude scale as the experimental stability constants. By contrast, MACE-MP0 MLIP underestimates lanthanide hydration numbers and gives PMF-derived stability constants with large deviations from experiment, whereas MACE-MATPES-R2SCAN improves both hydration structure and the stability-constant scale but still does not quantitatively reproduce the detailed lanthanide trend. Across both approaches, nitrate binding is best viewed as a labile coordination motif rather than a fixed mono- or bidentate structure, with hydration-shell structure influencing which configurations are favored. Overall, the cluster--continuum calculations provide a practical semi-quantitative baseline for the experimental stability-constant scale, while the explicit-solvent MLIP benchmarks show clear progress from MACE-MP0 to MACE-MATPES-r2SCAN but also highlight the need for lanthanide-targeted training, fine-tuning, or improved long-range and polarization treatments to obtain predictive thermodynamics for complex aqueous lanthanide chemistry.

Dinpajooh, Mohammadhasan↗
DOE OSTI2025

A review of collaborative studies between the NSTX/-U and MAST/-U spherical tokamaks

The National Spherical Torus Experiment (NSTX) at the Princeton Plasma Physics Laboratory in the United States, and the mega ampere spherical tokamak (MAST) at the United Kingdom Atomic Energy Authority in the United Kingdom, and their respective upgrades (NSTX-U and MAST-U) are two MAST fusion devices that have operated roughly over the past two decades. Both devices have made significant contributions to understanding spherical tokamak (ST) plasma physics, and fusion plasmas in general, and both have contributed data to multi-machine database studies. Several diagnostics have been physically moved from one machine to the other by diagnostic teams working on both devices. Collaboration has benefited both research teams in the areas of operational expertise, scenario development, and equilibrium reconstruction techniques. More focused comparative studies between the two devices have been pursued over the years in many areas as well, including stability calculations, disruption characterization, pedestal and edge localized mode stability, confinement and transport, energetic particles, and heating and current drive modelling. Together NSTX/-U and MAST/-U set the stage for the future of STs, which is entering the phase of design of demonstration power plant devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Type-I ELM mitigation by continuous lithium granule gravitational injection into the upper tungsten divertor in EAST

Large edge-localized modes (ELMs) were mitigated by gravitational injection of lithium granules into the upper X-point region of the experimental advanced superconducting tokamak (EAST) device with tungsten plasma-facing components. The maximum ELM size was reduced by ~70% in high β N H-mode plasmas. Large ELM stabilization was sustained for up to about 40 energy confinement times, with constant core radiated power and no evidence of high-Z or low-Z impurity accumulation. The lithium granules injection reduced the edge plasma pedestal density and temperature and their gradients, due to increased edge radiation and reduced recycling from the plasma-facing components. Ideal stability calculations using the ELITE code indicate that the stabilization of large ELMs correlates with improved stability of intermediate-n peeling-ballooning modes, due to reduced edge current resulting from the profile changes. The pedestal pressure reduction was partially offset by a core density increase, which resulted in a modest ~7% drop in core stored energy and normalized energy confinement time. We surmise that the remnant small ELMs are triggered by the penetration of multiple Li granules just past the separatrix, similar to small ELMs triggered by deuterium pellet Futatani et al (2014 Nucl. Fusion 54 073008). This study extends previous ELM elimination with Li powder injection Maingi et al (2018 Nucl. Fusion 58 024003) in EAST because (1) use of small, dust-like powder and the related potential health hazards were eliminated, and (2) use of macroscopic granules should be more applicable to future devices, due to deeper penetration than dust particles, e.g. inside the separatrix with velocities ~10 m s -1 in EAST.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Type-I ELM mitigation by continuous lithium granule gravitational injection into the upper tungsten divertor in EAST

Large edge-localized modes (ELMs) were mitigated by gravitational injection of lithium granules into the upper X-point region of the EAST device with tungsten plasma-facing components. The maximum ELM size was reduced by ~ 70% in high βN H-mode plasmas. Large ELM stabilization was sustained for up to about 40 energy confinement times, with constant core radiated power and no evidence of high-Z or low-Z impurity accumulation. The lithium granules injection reduced the edge plasma pedestal density and temperature and their gradients, due to increased edge radiation and reduced recycling from the plasma-facing components. Ideal stability calculations using the ELITE code indicate that the stabilization of large ELMs correlates with improved stability of intermediate-n peeling-ballooning modes, due to reduced edge current resulting from the profile changes. The pedestal pressure reduction was partially offset by a core density increase, which resulted in a modest ~ 7% drop in core stored energy and normalized energy confinement time. We surmise that the remnant small ELMs are triggered by the penetration of multiple Li granules just past the separatrix, similar to small ELMs triggered by deuterium pellet [S. Futatani et al., Nucl. Fusion 54 (2014) 073008]. This study extends previous ELM elimination with Li powder injection [R. Maingi et al., Nucl. Fusion 58 (2018) 024003] in EAST because 1) use of small, dust-like powder and the related potential health hazards were eliminated, and 2) use of macroscopic granules should be more applicable to future devices, due to deeper penetration than dust particles, e.g. inside the separatrix with velocities ~ 10 m/s in EAST.

EAST↗
NASA NTRS1983

Formulation and solution of rotary-wing aeroelastic stability and response problems

The state of the art in the formulation and solution of rotary-wing aeroelastic stability and response problems is reviewed in detail. The approximations used in the structural, inertia and aerodynamic operators are discussed. The important role of geometric nonlinearities, due to moderate deflections, and aerodynamic stall in the aeroelastic stability and response problem are identified. It is also shown that geometric nonlinearities are of primary importance in aeroelastic stability calculations, and have a more limited, though important, role in response calculations. Next, formulation of coupled rotor/fuselage problems is described, for both air and ground resonance type problems. Both topics, the isolated blade problem and the coupled rotor/fuselage problem, are treated for both hover and forward flight. Solution of aeroelastic stability and response problems proceeds in two stages. First, the spatial dependence is eliminated by using Galerkin's method, or by using the finite element method. Next the nonlinear, or linear, ordinary differential equation with periodic coefficients have to be solved for stability or response. Efficient numerical methods for accomplishing these objectives are presented in a comprehensive manner. The paper contains a number of illustrative numerical results which are intended to clarify various aspects of the modeling process and serve as representative results for both aeroelastic stability and response calculations for a variety of blade and rotor configurations.

Friedmann, P. P.↗
NASA NTRS2000

Interpreting Assemblages with Titanite (Sphene): It Does not have to be Greek to You.

Assemblages with titanite, pyroxene(s), olivine, ilmenite, magnetite, and quartz can be used to constrain the intensive and compositional variables that operate during crystallization. Such assemblages are relatively rare in metamorphic rocks, but they are more common in igneous rocks and more frequently reported in plutonic than volcanic rocks. We used the program QUILF, enhanced with thermodynamic data for titanite, to compute stable reactions among titanite (CaTiOSiO4), Fe-Mg-Ti ilmenite and magnetite (hereafter ilmenite and magnetite), Ca-Mg-Fe pyroxenes and olivine, and quartz, and to evaluate some of the factors that control titanite stability. Calculations at 1, 3, and 6 Kbar and 650, 850, 1100 0 C, in the system CaO - MgO - FeO Fe2O3 - TiO2 - SiO2, suggest that the reactions: Augitc + Ilmenite = Titanite + Magnetite + Quartz and Augite + Ilmenite + Quartz = Titanite + Orthopyroxene, impose well defined fugacity of O2, alpha(sub SiO2), , and compositional restrictions to the assemblages: (1) Titanite + Magnetite + Quartz, (2) Titanite + Orthopyroxene, (3) Augite + Ilmenite, and consequently titanite stability. From our calculations in this system we can draw the following general conclusions: (1) The assemblage Titanite + Magnetite + Quartz is always a good indicator of relatively high fugacity of O2, and it is likely more common in relatively iron-rich bulk compositions and for decreasing temperature and pressure conditions. (2) At high temperatures (>= 650 C) titanite is not stable in quartz-saturated rocks that contain the assemblage Orthopyroxene + Augite + Ilmenite + Magnetite. (3) In quartz-saturated rocks the coexistence of titanite and magnetite with either orthopyroxene or olivine requires a confluence of conditions relating bulk composition, fugacity of O2, and slow cooling. Thus, such assemblages must be rare. (4) Regardless of T and fugacity of O2 conditions, and bulk-composition, titanite is not stable in quartz-absent rocks that contain Olivine + Orthopyroxene + Augite + Ilmenite + Magnetite. Decreasing temperature and pressure conditions appear to favor titanite crystallization, thus, it is not unsurprising that titanite is frequently observed in slowly cooled rocks, albeit, in association with amphibole. We argue that the titanite + amphibole association is likely favored by high water activity, regardless of oxygen fugacity. Because water activity increases during crystallization of a pluton, the association titanite + amphibole, and consequently titanite, is likely to be more common in plutonic rocks than in volcanic rocks.

Xirouchakis, Dimitrios M.↗
DOE OSTI2025

Doubling fusion power with volumetric optimization in magnetic confinement fusion devices

A technique, volumetric power optimization, is presented for enhancing the power output of magnetic confinement fusion devices. Applied to a tokamak, this approach involves shifting the burning plasma region to a larger plasma volume while introducing minimal perturbations to the plasma boundary shape. This edge perturbation—squareness—is analogous to pinching and stretching the edge boundary. Stability calculations confirm that this edge alteration is compatible with maintaining plasma stability. This optimization method for optimizing fusion power output could improve the performance of magnetic confinement fusion power plants.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Thresholds for the onset of fluid and magnetofluid turbulence

Linear stability calculations conducted in plasma physics are based on the theory of hydrodynamic stability of neutral fluids. The present investigation is concerned with the validity of procedures based on linear stability analysis in six much-studied situations. All have the property that the fluid goes from laminar to turbulent at critical values of some dimensionless number, such as the Reynolds number or the Rayleigh number. The question is, at what threshold values of the dimensionless number do the unstable motions set in, and can these thresholds be predicted by a linear analysis of the stability of the laminar state as the threshold is approached from the stable side. In three cases linear stability analysis clearly seems to fail. These cases include Plane Poiseuille Flow, Plane Couette Flow, and Cylindrical Pipe Flow (Hagen-Poiseuille Flow). Situations in which predictions provided by linear stability analysis are correct are related to Rotating Couette Flow, Thermally-Driven Convection, and Instability of Laminar Boundary Layers.

Montgomery, D.↗
NASA NTRS1986

Calculating Aerodynamic-Stability Derivatives

VORSTAB program developed to calculate lateral-directional characteristics of nonplanar wing/body combinations in subsonic flow. Mathematically determines effects of edge-separated vortex flow, including augmented vortex lift, strake-induced downwash, and vortex breakdown. VORSTAB written in FORTRAN IV.

Lan, C. E.↗
NASA NTRS1974

Measurements of the short-term stability of quartz crystal resonators: A window on future developments in crystal oscillators

Recent measurements of the inherent short-term stability of quartz crystal resonators are presented. These measurements show that quartz resonators are much more stable for times less than 1s than the best available commercial quartz oscillators. A simple model appears to explain the noise mechanism in crystal controlled oscillators and points the way to design changes which should permit more than 2 orders of magnitude improvement in their short-term stability. Calculations show that a reference signal at 1 THz, derived from frequency multiplying a 5 MHz source with the above measured crystal stability, should have an instantaneous or fast linewidth of order 1 Hz. These calculations explicitly include the noise contribution of our present multiplier chains and are briefly outlined.

Walls, F. L.↗
NASA NTRS1980

Stability analysis for laminar flow control, part 2

Topics covered include: (1) optimization of the numerics of the SALLY stability analysis code; (2) relation between temporal and spatial stability theory; (3) compressible flow stability calculations; (4) spectral methods for the boundary layer equations; and (5) numerical study of nonlinear, nonparallel stability of incompressible flows.

Orszag, S. A.↗
DOE OSTI2026

Density limit in peeling-limited pedestals at and above the Greenwald value in DIII-D high poloidal beta plasmas

High pressure, peeling limited pedestals with pedestal normalized beta β N,ped >2 and pedestal top density n ped at or above the Greenwald density nG have been achieved in DIII-D high poloidal beta plasmas, with high global normalized beta β N >3 and energy confinement H 98 ~1.2-1.7. Higher β N allows higher pedestal density above the Greenwald value and higher pedestal pressure, even with a low injected torque. MHD modeling confirms that the experimental profiles lie near the peeling-mode unstable boundary with high normalized pressure gradient βMHD and high edge current density. Experimental analysis and stability calculations indicate that the high poloidal beta with strong Shafranov shift, high βMHD and weak/negative magnetic shear improves the pedestal stability by decoupling the peeling and ballooning modes and stabilizing the ballooning modes, thus facilitating access to the second stable region of peeling-ballooning mode. The access to the second stable peeling-ballooning stability region opens Super-H-like channels without extremely strong shaping or strong torque injection. The high-pressure peeling pedestal allows the pedestal density to go beyond the Greenwald limit with strong ExB shear maintained: pedestal pressure increases with pedestal density even when n ped >n G , until reaching the ideal MHD instability boundary, where giant ELMs occur. The giant ELMs are dominated by a strong n=1 component and cause a large reduction of the edge pressure, but a negligible change of the core pressure, consistent with kink/peeling-mode induced instability. The pedestal recovers from the collapse and typically sustains a high baseline density, around the Greenwald limit, during the whole discharge duration. Experiments also found that internal transport barriers and n ped ~0.9nG, peeling limited pedestals could be simultaneously achieved in high β N plasmas, while an internal feedback between ITB strength and pedestal performance is found.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1947

Systematic investigations of the influence of the shape of the profile upon the position of the transition point.

The position of the beginning of transition laminar/turbulent as a function of the thickness and the camber of the profile at various Reynolds numbers and lift coefficients was investigated for a series of Joukowsky profiles. The calculation of the boundary layer was carried out according to the Pohlhausen method which may be continued by a simplified stability calculation according to H. Schlichting. A list of tables is given which permits the reading off of the position of the transition point on suction and pressure side for each Joukowsky profile.

Bussmann, K↗