Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stability boundary”

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 253 records · Page 14

Mechanical and Tribological Performance of Additively Manufactured Nanocrystalline Aluminum via Cryomilling and Cold Spray

In this study, nanocrystalline (NC) aluminum (Al) and magnesium (Mg)-doped Al bulk components were fabricated using a hybrid manufacturing process that combines cryomilling and high-pressure cold spray (HPCS) additive deposition techniques. Yttria-stabilized zirconia (YSZ) was also added during the HPCS process to improve deposition efficiency and build-up thickness via peening. The evolution of morphology, crystallite size, and elemental composition of both cryomilled powders and cold-sprayed (CS’ed) components was systematically characterized using X-ray Diffraction (XRD), Scanning Electron Microscopy (SEM), and Transmission Electron Microscopy (TEM). Mechanical characterization was performed using Vickers microhardness and uniaxial tensile testing, while the tribological behavior was assessed using sliding wear tests under dry/lubricated conditions. XRD analysis revealed that increased cryomilling duration led to significant crystallite refinement, which directly correlated with enhanced hardness and strength. This mechanical strengthening was accompanied by an increase in coefficient of friction (COF) and lower wear rates. The results also showed that the Mg-doped Al exhibited superior hardness, tensile strength, and tribological performance compared to pure Al. The study further explores the underlying mechanisms responsible for these enhancements, highlighting the potential of solute-assisted grain boundary stabilization in tailoring high-performance NC Al alloys.

36 MATERIALS SCIENCE↗

Achieving a robust grassy-ELM operation regime in CFETR

We have identified a robust grassy-edge localized mode (ELM) operation regime for future tokamak reactors. The regime exists within a pedestal top electron collisionality (ν*) window at high global poloidal beta (β p ). The existence of an upper ν* limit for grassy-ELMs is consistent with results previously reported in experiments (Oyama et al 2010 Nucl. Fusion 50 064014), while the existence of a lower ν* limit has not been reported previously. Using EPED and BOUT++, a theoretical model that quantitatively explains the physics of the grassy-ELMs within the window, which distinguishes them from the small mixed-ELMs at lower ν*, is presented for the first time. A peeling-ballooning stability boundary is obtained by scanning the operating density space. The change in density corresponds to a change in ν* that affects the pedestal bootstrap current. High β p leads to a strong Shafranov shift, which affects the flux surface averaged pressure drive. The two effects combine to create a peeling-dominated window in intermediate ν* buffered by ballooning-dominated regimes. Only the peeling-dominated regime shows a cyclic behavior in the perturbed pressure during the nonlinear simulation of an ELM crash, reminiscent of grassy-ELM dynamics. Similarly, the energy released across the separatrix is demonstrated to be significantly smaller. The quick recovery of the ELM crash is explainable by the rapid rise of a low n kink-peeling instability when the pedestal current I ped exceeds a threshold at high β p . It minimizes the excursion beyond marginal stability and is absent in the ballooning-dominated regime. Comparison with recent experiments over a range of β p and ν* strongly supports the physical picture proposed by the modeling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aeroservoelastic Uncertainty Model Identification from Flight Data

Uncertainty modeling is a critical element in the estimation of robust stability margins for stability boundary prediction and robust flight control system development. There has been a serious deficiency to date in aeroservoelastic data analysis with attention to uncertainty modeling. Uncertainty can be estimated from flight data using both parametric and nonparametric identification techniques. The model validation problem addressed in this paper is to identify aeroservoelastic models with associated uncertainty structures from a limited amount of controlled excitation inputs over an extensive flight envelope. The challenge to this problem is to update analytical models from flight data estimates while also deriving non-conservative uncertainty descriptions consistent with the flight data. Multisine control surface command inputs and control system feedbacks are used as signals in a wavelet-based modal parameter estimation procedure for model updates. Transfer function estimates are incorporated in a robust minimax estimation scheme to get input-output parameters and error bounds consistent with the data and model structure. Uncertainty estimates derived from the data in this manner provide an appropriate and relevant representation for model development and robust stability analysis. This model-plus-uncertainty identification procedure is applied to aeroservoelastic flight data from the NASA Dryden Flight Research Center F-18 Systems Research Aircraft.

Brenner, Martin J.↗

Method for Selecting Rotor Suspension Design Criteria

A method for selecting rotor suspension design criteria that allow rotors with internal friction to negotiate critical speeds and avoid subsynchronous internal friction instabilities is described. A rotor model of a test rotor with internal friction is created and used to calculate the test rotor’s critical speeds, mode shapes, Campbell plots, and the rotational speed corresponding to the stability boundary for a nominal set of suspension parameters. The test rotor model is then used to demonstrate the suspension design criteria selection method.

Damiano, Brian↗

In-situ radiation response of additively manufactured modified Inconel 718 alloys

Here in this study, a novel alloy of modified Inconel 718 produced by laser powder bed fusion is studied before and after in-situ Kr irradiation up to 3 dpa at 200 and 450 °C. Before irradiation, the microstructure consists of dislocation cells having a misorientation angle less than 5° and with an average size of ~500 nm. There are also second phase particles of MC type carbides, Laves phase and oxides such as Y-O, Y-(Ti)-Al-O. While the microstructure consists of stacking fault tetrahedra, faulted and perfect loops after irradiation at 200 °C, dislocation loops are the primary defects at 450 °C. With increasing dose, the size of the defects remains similar at 200 °C while it increases at 450 °C. This has been attributed to the existence of vacancy type defects at 200 °C and the different defect transport mechanisms at different temperatures. Moreover, matrix and second phase particle compositions seem to be similar after irradiation. The sink strengths of the structures have been calculated and superior radiation resistance of this alloy has been attributed to the existence of fine cell boundaries stabilized by the second phase particles produced by additive manufacturing.

36 MATERIALS SCIENCE↗

Geochemical-hydromechanical couplings during water-serpentinized harzburgite interactions at 20° C and 30 bars

Interaction of groundwater with the serpentinized harzburgites of the upper mantle plays an important role in the geologic carbon cycle and serpentinization processes. Here, in this study, an intact serpentinized harzburgite with more than 60% lizardite from the Semail ophiolite in Oman was reacted with water at near neutral pH and low PCO 2 to observe its dissolution behavior and reaction path within the CaO-MgO-SiO 2 -CO 2 -H 2 O system at 20 °C and 30 bars. In the experiment, mass (i.e. H 2 O, ions) transfer is by diffusion and the selected temperature and pressure conditions mimic water-rock interactions at the shallow portions of most ultramafic rock hosted aquifers. Equilibrium activity-activity diagrams and thermodynamic data were used to compare the experimentally determined reaction path with the theoretical stability boundaries for minerals involved to investigate the evolution of bulk chemical composition, mineralogy, and water composition. The experiments demonstrated that (i) the water-rock interactions increased the pH of the aqueous phase from 5.9 to 7.5 over a period of 18 weeks; (ii) the dissolution of lizardite, orthopyroxene (enstatite), and clinopyroxene (diopside and augite) increased the concentration of Mg, Ca, and Si in the aqueous phase; (iii) the dissolution was incongruent with respect to Mg and Si, favoring Si release at pH > 6 due to (a) breaking of more reactive cation‑oxygen bonds that is consistent with the stoichiometry of magnesium for proton exchange reaction which favors a surface that is enriched in Mg at basic conditions, and (b) preferential dissolution of clinopyroxene Mg 1.30 Ca 0.52 Fe 0.06 Si 2 O 6 ; (iv) the aqueous phase was undersaturated with respect to carbonate (e.g. calcite, magnesite, and hydromagnesite) and hydrous (e.g. lizardite, chrysotile, brucite, and talc) minerals; (v) the bulk chemical composition and mineralogy of the intact serpentinized harzburgite matrix did not change; (vi) the thermodynamic data can successfully predict water-intact serpentinized harzburgite behavior if the water chemical composition can be constrained; and (vii) no detectable macroscopic form of damage (e.g. microcracks forming as a result of differential stress fields due to changes in volume during chemical reaction) was observed.

Carbonation↗

An experimental system and procedure of unsteady-state relative permeability test for gas hydrate-bearing sediments

Reliable estimations of the relative permeability of gas and water in hydrate-bearing sediments (HBS) and the dependency of the relative permeability on hydrate saturation are critical to predict the productivity of a hydrate reservoir. Yet, this remains poorly estimated owing to lack of experimental data associated with difficulties in conducting multiphase flow experiments in HBS. Recognizing the experimental challenges, this study intends to develop and validate a new experimental system and procedure of unsteady-state relative permeability test that can generate reliable and reproducible flow measurements in HBS. Furthermore, gas hydrate is considered as a part of solid matrix in the sediment, so one of the challenges is to maintain a constant hydrate saturation, which is achieved in this experimental study using tight pressure-temperature (P-T) control near the hydrate stability boundary. The measured differential pressure across the specimen, methane injection flow rate, and volume of displaced brine are used to calculate the relative permeability by adopting a conventional Buckley-Leverett theory-based interpretation method. Residual brine saturation calculated for the hydrate-bearing specimen is higher than that of hydrate-free specimen, presumably due to decrease in pore size, increase in heterogeneity of solid matrix, and increase in size distribution of solid matrix and pore in the presence of hydrates. Further studies are necessary to represent the results of the unsteady-state flow experiment in HBS with a gas hydrate-dependent relative permeability model.

03 NATURAL GAS↗

Multimode theory of electron hole transverse instability

We present Vlasov–Poisson three-dimensional linear stability analysis of an initially planar electron hole structure, solving for the distribution function by integration along unperturbed orbits. The non-sinusoidal potential perturbation shape (parallel to $B$ ) is expanded in eigenfunctions of the adiabatic Poisson operator, generalizing the prior assumption of a rigid shift of the equilibrium. We show that the shiftmode is then modified by a second discrete mode plus an integral over a continuum of wave-like modes. A rigorous treatment shows that the continuum can be approximated effectively by a single mode that satisfies the external wave dispersion relation, thus making the perturbation a weighted sum of three modes. We find numerically the solution for the complex instability frequency, and the corresponding three mode amplitudes determining the perturbation eigenmode. This multimode analysis refines the accuracy of the prior single-mode results, giving slightly higher growth rates at most parameters, as expected from the extra mode shape freedom. Oscillating modes near stability boundaries have larger mode distortions which help explain particle-in-cell simulations that observe instability up to ${\sim }20$ % beyond the prior shiftmode thresholds, and narrowing of the perturbation. At high magnetic field, the multimode analysis predicts a reduction of the already small growth rate.

Physics↗

How does ion temperature gradient turbulence depend on magnetic geometry? Insights from data and machine learning

Magnetic geometry has a significant effect on the level of turbulent transport in fusion plasmas. Here, we model and analyse this dependence using multiple machine learning methods and a dataset of >200 000 nonlinear gyrokinetic simulations of ion-temperature-gradient turbulence in diverse non-axisymmetric geometries. The dataset is generated using a large collection of both optimised and randomly generated stellarator equilibria. At fixed gradients and other input parameters, the turbulent heat flux varies between geometries by several orders of magnitude. Trends are apparent among the configurations with particularly high or particularly low heat flux. Regression and classification techniques from machine learning are then applied to extract patterns in the dataset. Due to a symmetry of the gyrokinetic equation, the heat flux and regressions thereof should be invariant to translations of the raw features in the parallel coordinate, similar to translation invariance in computer vision applications. Multiple regression models including convolutional neural networks (CNNs) and decision trees can achieve reasonable predictive power for the heat flux in held-out test configurations, with highest accuracy for the CNNs. Using Spearman correlation, sequential feature selection and Shapley values to measure feature importance, it is consistently found that the most important geometric lever on the heat flux is the flux surface compression in regions of bad curvature. The second most important geometric feature relates to the magnitude of geodesic curvature. These two features align remarkably with surrogates that have been proposed based on theory, while the methods here allow a natural extension to more features for increased accuracy. The dataset, released with this publication, may also be used to test other proposed surrogates, and we find that many previously published proxies do correlate well with both the heat flux and stability boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tandem Mass-Selective Cryogenic Digital Ion Traps for Enhanced Cluster Formation

We present the implementation of tandem mass-selective cryogenic ion traps, designed to enhance the range of ion processing capabilities that can be performed prior to spectroscopic interrogation. We show that both the formation of ion clusters and mass filtering steps can be combined in a single cryogenic linear quadrupole ion trap driven by RF square waves. Mass filtering and mass isolation can be achieved by manipulation of the RF frequency and duty cycle. Very importantly, this scheme circumvents the need for high-amplitude RF voltages that can be incompatible with typical cryogenic ion processing conditions. In addition, proper adjustment of the stability boundaries during the clustering process allows for the preferential formation of a specific cluster size rather than a broad distribution of sizes. Lastly, we show that a specific cluster size can be formed, mass-selected, and then transferred to another ion trap for a second completely separate ion processing step. Furthermore, the instrumentation and modular design developed here expand the scope of ionic species and clusters that can be accessed by processing electrosprayed ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the new universality class in structurally disordered n -vector model with long-range interactions

We study a stability boundary of a region where nontrivial critical behavior of an n-vector model with long-range power-law decaying interactions is induced by the presence of a structural disorder (e.g., weak quenched dilution). This boundary is given by the marginal dimension of the order parameter nc dependent on space dimension, d, and a control parameter of the interaction decay, σ, below which the model belongs to the new dilution-induced universality class. Exploiting the Harris criterion and recent field theoretical renormalization group results for the pure model with long-range interactions, we get nc as a three loop ɛ = 2σ – d-expansion. We provide numerical values for nc applying series resummation methods. Our results show that not only the Ising systems (n = 1) can belong to the new disorder-induced long-range universality class at d = 2 and 3.

Physics↗

Nonlinear bubble competition of the multimode ablative Rayleigh–Taylor instability and applications to inertial confinement fusion

The self-similar nonlinear evolution of the multimode ablative Rayleigh–Taylor instability (RTI) and the ablation-generated vorticity effect are studied for a range of initial conditions. We show that, unlike classical RTI, the nonlinear multimode bubble-front evolution remains in the bubble competition regime due to ablation-generated vorticity, which accelerates the bubbles, thereby preventing a transition into the bubble-merger regime. We develop an analytical bubble competition model to describe the linear and nonlinear stages of ablative RTI. Here, we show that vorticity inside the multimode bubbles is most significant at small scales with large initial perturbation. Since these small scales persist in the bubble competition regime, the self-similar growth coefficient ab can be enhanced by up to 30% relative to ablative bubble competition without vorticity effects. We use the ablative bubble competition model to explain the hydrodynamic stability boundary observed in OMEGA low-adiabat implosion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ideal MHD induced temperature flattening in spherical tokamaks

This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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

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

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On hybrid scenarios in KSTAR

We report the status of hybrid scenario experiments in Korea Superconducting Tokamak Advanced Research (KSTAR). The hybrid scenario is defined as stationary discharges with ${\beta _{\rm{N}}}$ $ \geqslant $ 2.4 and ${H_{89}}$ $ \geqslant $ 2.0 at ${q_{95}}$ < 6.5 without or with very mild sawtooth activities in KSTAR. It is being developed towards reactor-relevant conditions. High performance of ${\beta _{\rm{N}}} \lesssim 3.0$, ${H_{89}} \lesssim 2.4$ and G-factor ($ \equiv {\beta _{\rm{N}}}{H_{89}}/q_{95}^2$) $ \lesssim 0.4$6 has been achieved and sustained for $ \gtrsim 40{{\tau }_{\text{E}}}$ at ne/nGW ~ 0.7 with heating power of $ \lesssim $5 MW. Some KSTAR hybrid discharges exhibit a unique feature of a slow transition from conventional H-mode to hybrid mode after the third neutral beam injection. The reason for the confinement enhancement is extensively studied in this transition period of a representative discharge exhibiting a common feature of KSTAR hybrid scenarios. 0D performance analysis with magnetohydrodynamic activities, 1D kinetic profile dynamics, power balance analysis, linear gyro-kinetic analysis and edge pedestal stability analysis were conducted. The enhancement is thought to be from both the core and the pedestal. The improvement in the core region of the ion energy channel is observed from the linear gyro-kinetic analysis considering the electromagnetic, the fast ion, the Shafranov shift, ${\it{\omega }_{E \times B}}$, and the magnetic shear effect. The electromagnetic finite $\beta $ stabilisation plays a role in the inner core region at ${\rho _{\rm{tor}}}\sim0.35$ together with the fast ion effect. The alpha stabilisation effect is also found at ${\rho _{\rm{tor}}}\sim0.5$. ${\omega _{E \times B}}$, which could reduce the linear growth of the ion temperature gradient mode in the outer core region at ${\rho _{\rm{tor}}}\sim0.5 - 0.7$ with the highest contribution from the toroidal rotation. Furthermore, regarding the improvement in the pedestal, Shafranov shift broadens the stability boundary of the pedestal in support of the diamagnetic effect. The pedestal height and width could be reproduced by the EPED model, while a realistic current profile is used to calculate the internal inductance for Shafranov shift. Based on these findings, a comprehensive confinement enhancement mechanism has been proposed by considering the core-edge interplay.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Isotope dependence of the type I ELMy H-mode pedestal in JET-ILW hydrogen and deuterium plasmas

The pedestal structure, edge transport and linear MHD stability have been analyzed in a series of JET with the ITER-like wall hydrogen (H) and deuterium (D) type I ELMy H-mode plasmas. The pedestal pressure is typically higher in D than in H at the same input power and gas rate, with the difference mainly due to lower density in H than in D (Maggi et al (JET Contributors) 2018 Plasma Phys. Control. Fusion 60 014045). A power balance analysis of the pedestal has shown that higher inter-ELM separatrix loss power is required in H than in D to maintain a similar pedestal top pressure. This is qualitatively consistent with a set of interpretative EDGE2D-EIRENE simulations for H and D plasmas, showing that higher edge particle and heat transport coefficients are needed in H than in D to match the experimental profiles. It has also been concluded that the difference in neutral penetration between H and D leads only to minor changes in the upstream density profiles and with trends opposite to experimental observations. This implies that neutral penetration has a minor role in setting the difference between H and D pedestals, but higher ELM and/or inter-ELM transport are likely to be the main players. The interpretative EDGE2D-EIRENE simulations, with simultaneous upstream and outer divertor target profile constraints, have indicated higher separatrix electron temperature in H than in D for a pair of discharges at low fueling gas rate and similar stored energy (which required higher input power in H than in D at the same gas rate). The isotope dependence of linear MHD pedestal stability has been found to be small, but if a higher separatrix temperature is considered in H than in D, this could lead to destabilization of peeling-ballooning modes and shrinking of the stability boundary, qualitatively consistent with the reduced pedestal confinement in H.

08 HYDROGEN↗

A novel path to runaway electron mitigation via deuterium injection and current-driven MHD instability

Relativistic electron (RE) beams at high current density (low safety factor, qa) yet very low free-electron density accessed with D2 secondary injection in the DIII-D and JET tokamak are found to exhibit large-scale MHD instabilities that benignly terminate the RE beam. In JET, this technique has enabled termination of MA-level RE currents without measurable first-wall heating. This scenario thus offers an unexpected alternate pathway to achieve RE mitigation without collisional dissipation. Benign termination is explained by two synergistic effects. First, during the MHD-driven RE loss events both experiment and MHD orbit-loss modeling supports a significant increase in the wetted area of the RE loss. Second, as previously identified at JET and DIII-D, the fast kink loss timescale precludes RE beam regeneration and the resulting dangerous conversion of magnetic to RE kinetic energy. During the termination, the RE kinetic energy is lost to the wall, but the current fully transfers to the cold bulk thus enabling benign Ohmic dissipation of the magnetic energy on longer timescales via a conventional current quench. Hydrogenic (D2) secondary injection is found to be the only injected species that enables access to the benign termination. D2 injection: 1) facilitates access to low qa in existing devices (via reduced collisionality & resistivity), 2) minimizes the RE avalanche by ‘purging’ the high-Z atoms from the RE beam, 3) drives recombination of the background plasma, reducing the density and Alfven time, thus accelerating the MHD growth. Furthermore, this phenomenon is found to be accessible when crossing the low qa stability boundary with rising current, falling toroidal field, or contracting minor radius - the latter being the expected scenario for vertically unstable RE beams in ITER. While unexpected, this path scales favorably to fusion-grade tokamaks and offers a novel RE mitigation scenario in principle accessible with the day-one disruption mitigation system (DMS) of ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A dimensionality reduction algorithm for mapping tokamak operational regimes using a variational autoencoder (VAE) neural network

Abstract A variational autoencoder (VAE) is a type of unsupervised neural network which is able to learn meaningful data representations in a reduced dimensional space. We present an application of VAE in identifying the operational stability boundary of tokamak plasma discharges. This model was implemented using a dataset of over 3000 discharges from the high beta tokamak-extended pulse (HBT-EP) device. We found the VAE model to be capable of forming a continuous low-dimensional operational space map and identifying the operational boundaries using a specified warning time window. By projecting the operational parameters onto the same reduced space, this provides an intuitive way for the machine operator or an automated control system to perform disruption avoidance using a relevant control actuator as a discharge approaches a boundary. Pre-programmed GPU control experiments were conducted to demonstrate this control technique using HBT-EP’s saddle control coils as a horizontal position actuator, showing the ability to avoid the oncoming disruptive event and extend the duration of the discharge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗