Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “stability 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 1,333 records · Page 74

Prediction of the energetic particle redistribution by an improved critical gradient model and analysis of the transport threshold

Based on the theory of critical gradient model (CGM) and following the simulation method proposed by Waltz et al., a combination of TGLFEP and EPtran code is employed to predict the energetic particle (EP) transport induced by Alfvén eigenmodes (AEs). To be consistent with the experiment, recent improvements to the simulation method include consideration of threshold evolution and orbit loss due to finite orbit width. The revised CGM is applied to simulate two DIII-D experimental discharges (#142111 and #153071). It well reproduces the experimental profiles with multiple unstable AEs and large-scale EP transport. Discharge #142111 had previously been simulated using a nonlinear MHD-kinetic code MEGA with a transport mechanism based on stochasticity induced by overlapping AE. By comparing the simulated EP profiles, we find that the AE transport threshold is approximated by both the MEGA nonlinear stability threshold and the proposed CGM threshold (error <5% for single n and <17% for multiple n simulation). Both of them are larger than the linear stability threshold of the most unstable AE mode by a quantity of the order of the flux needed to sustain EP transport by the background turbulence. Additionally, we have also applied the improved CGM to simulate the α particle redistribution for a China Fusion Engineering Test Reactor steady state scenario. Because of the clear separation between the AE unstable region and the loss cone, only a moderate α particle loss of ~9.6% is predicted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mathematical solutions in internal dose assessment: A comparison of Python-based differential equation solvers in biokinetic modeling

Abstract In biokinetic modeling systems employed for radiation protection, biological retention and excretion have been modeled as a series of discretized compartments representing the organs and tissues of the human body. Fractional retention and excretion in these organ and tissue systems have been mathematically governed by a series of coupled first-order ordinary differential equations (ODEs). The coupled ODE systems comprising the biokinetic models are usually stiff due to the severe difference between rapid and slow transfers between compartments. In this study, the capabilities of solving a complex coupled system of ODEs for biokinetic modeling were evaluated by comparing different Python programming language solvers and solving methods with the motivation of establishing a framework that enables multi-level analysis. The stability of the solvers was analyzed to select the best performers for solving the biokinetic problems. A Python-based linear algebraic method was also explored to examine how the numerical methods deviated from an analytical or semi-analytical method. Results demonstrated that customized implicit methods resulted in an enhanced stable solution for the inhaled 60 Co (Type M) and 131 I (Type F) exposure scenarios for the inhalation pathway of the International Commission on Radiological Protection (ICRP) Publication 130 Human Respiratory Tract Model (HRTM). The customized implementation of the Python-based implicit solvers resulted in approximately consistent solutions with the Python-based matrix exponential method ( expm ). The differences generally observed between the implicit solvers and expm are attributable to numerical precision and the order of numerical approximation of the numerical solvers. This study provides the first analysis of a list of Python ODE solvers and methods by comparing their usage for solving biokinetic models using the ICRP Publication 130 HRTM and provides a framework for the selection of the most appropriate ODE solvers and methods in Python language to implement for modeling the distribution of internal radioactivity.

61 RADIATION PROTECTION AND DOSIMETRY↗

Resistive wall tearing mode disruptions in DIII-D and ITER tokamaks

Disruptions are a serious problem in tokamaks, in which thermal and magnetic energy confinement is lost. This paper uses data from the DIII-D experiment, theory, and simulations to demonstrate that resistive wall tearing modes (RWTMs) produce the thermal quench (TQ) in a typical locked mode shot. Analysis of the linear RWTM dispersion relation shows the parameter dependence of the growth rate, particularly on the resistive wall time. Linear simulations of the locked mode equilibrium show that it is unstable with a resistive wall and stable with an ideally conducting wall. Nonlinear simulations demonstrate that the RWTM grows to sufficient amplitude to cause a complete thermal quench. The RWTM growth time is proportional to the thermal quench time. Here, the nonlinearly saturated RWTM magnetic perturbation amplitude agrees with experimental measurements. The onset condition is that the q = 2 rational surface is sufficiently close to the resistive wall. Collectively, this identifies the RWTM as the cause of the TQ. In ITER, RWTMs will produce long TQ times compared to present-day experiments. ITER disruptions may be significantly more benign than previously predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

New Safety Feedback Control Design to Guarantee Adequate Frequency Performance in Microgrids

Safety analysis of power systems is concerned with the system's ability to maintain critical variables within specified limits following a disturbance. Frequency control adequacy has become increasingly important as the system inertia decreases due to the increase in renewable energy penetration. Various controllers for inverters have been proposed to improve the system frequency response and few are capable to ensure the safety of the response. In this article, a diesel-wind energy system is considered and modeled as a switching system between normal, faulted, and post-fault modes. Further, a safety feedback controller is designed as a supplementary signal for a wind turbine generator such that the speed of the diesel generator stays within a permissible range in the presence of a finite energy disturbance. Numerical results on the modified 33-bus microgrid system obtained of the proposed novel approach indicate that the suggested control configuration can guarantee adequate frequency response without excessive conservativeness.

barrier function↗

Sufficient Conditions for the Stability of a Lur'e System With Two Nonzero Inputs Based on Lyapunovs Direct Method

The Lur'e problem is a famous problem in nonlinear control theory. The analysis of this system is usually restricted to Lur'e systems in which the external input is zero. One form of this problem involves the stability analysis of the Lur'e system with two nonzero inputs. The stability analysis of this type of system is traditionally performed based on L stability theory. The objective of this paper is to present sufficient conditions for the stability of this system based on Lyapunov's direct method. This is a conceptually simpler approach to the existing analysis based on L stability theory and it provides some additional information about the stability of the system.

Adetona, Olawale↗

Dynamics of Orbits near 3:1 Resonance in the Earth-Moon System

The Interstellar Boundary Explorer (IBEX) spacecraft is currently in a highly elliptical orbit around Earth with a period near 3:1 resonance with the Moon. Its orbit is oriented so that apogee does not approach the Moon. Simulations show this orbit to be remarkably stable over the next twenty years. This article examines the dynamics of such orbits in the Circular Restricted 3-Body Problem (CR3BP). We look at three types of periodic orbits, each exhibiting a type of symmetry of the CR3BP. For each of the orbit types, we assess the local stability using Floquet analysis. Although not all of the periodic solutions are stable in the mathematical sense, any divergence is so slow as to produce practical stability over several decades. We use Poincare maps with twenty-year propagations to assess the nonlinear stability of the orbits, where the perturbation magnitudes are related to the orbit uncertainty for the IBEX mission. Finally we show that these orbits belong to a family of orbits connected in a bifurcation diagram that exhibits exchange of stability. The analysis of these families of period orbits provides a valuable starting point for a mission orbit trade study.

Dichmann, Donald J.↗

Networks of electrochemical oxidation of common lithium-ion Battery solvents revealed by NMR spectroscopy

Raising the upper cutoff voltage of lithium-ion batteries (LIBs) to increase energy density often exceeds the electrolyte’s anodic stability limit, accelerating degradation and creating a major durability tradeoff. Designing electrolytes that can sustain long-term high-voltage cycling requires a clearer understanding of the fundamental mechanisms occurring when commercial carbonate solvents oxidize. Here, to this end, simplified single-salt, single-solvent formulations of LiClO 4 and LiPF 6 in dimethyl carbonate (DMC), ethylene carbonate (EC), or ethyl methyl carbonate (EMC) were anodically electrolyzed on inert electrodes and monitored for extended periods of time using 1 H, 13 C, 19 F, and 35 Cl nuclear magnetic resonance (NMR) spectroscopy. The controlled environment of the experiments, coupled to the unique sensitivity of NMR, unveiled novel metastable intermediates and the formation of branching networks of products with temporal evolution. Oxidation of the pristine solvent primarily proceeds through a radical pathway that also produces highly reactive protons but faces competition from a second pathway involving a radical carbocation intermediate. In all cases, the intermediates follow a variety of downstream pathways that can intersect with each other. The concomitant network of reactions represents a significant increase in complexity compared to common descriptions in the literature, yet, critically, it helps explain the wide range of products typically identified in electrolyte oxidation in complete cells. The results highlight the need for refocusing fundamental research on anodic stability to analysis of the hierarchy of reaction networks to better inform efforts to mitigate the detrimental effects on battery performance, including prevention and harvesting of proton and radical products.

Electrolytes↗

Properties and Autoignition Reactivity of Diesel Boiling Range Ethers Produced from Guerbet Alcohols

We examine the properties of diesel boiling range ethers made from coupling of alcohols produced by oligomerization of ethanol (Guerbet alcohols) for their utility as low-carbon liquid fuel blendstocks. Basic properties of boiling point, flash point, freezing point, density and viscosity are well suited for blending into diesel fuels. For the mixture of ethers the lightest component, di-n-butyl ether, can be present at up to 20 vol% while still having adequately high flashpoint for safe handling. Soot formation tendency (as yield sooting index) is well below that of conventional diesel. The ethers have similar compatibility with elastomers as conventional diesel, based on Hansen solubility parameter analysis. Oxidation stability was assessed for 30 vol% blends of individual ethers in a conventional diesel fuel using a long-term storage test. Over 6 weeks we observed no formation of peroxides or degradation. n-alkyl ethers with carbon number of 8 or higher have cetane number over 100, which is outside the defined range of cetane number, while branched ethers are over 70. The ethers also blend antagonistically into conventional diesel for cetane number, meaning that the blend cetane value is lower than predicted based on a linear by volume, mass, or mole model. We show that aromatics and naphthenes likely act as radical scavengers to slow or shut down autoignition of the highly reactive ethers at low to medium blend levels. Overall, diesel boiling range ethers show significant promise as high quality low-net carbon diesel blendstocks.

09 BIOMASS FUELS↗

Omega-Limit Sets and Input-To-State Stability in Power Grids with Switching Equilibria

This paper studies a power transmission system with both conventional generators (CGs) and distributed energy assets (DEAs) providing frequency control. We consider an operating condition with demand aggregating two dynamic components: one that switches between different values on a finite set, and one that varies smoothly over time. Such dynamic operating conditions may result from protection scheme activations, external cyber-attacks, or due to the integration of dynamic loads, such as data centers. Mathematically, the dynamics of the resulting system are captured by a system that switches between a finite number of vector fields - or modes -, with each mode having a distinct equilibrium point induced by the demand aggregation. To analyze the stability properties of the resulting switching system, we leverage tools from hybrid dynamic inclusions and the concept of ..omega.. -limit sets from sets. Specifically, we characterize a compact set that is semi-globally practically asymptotically stable under the assumption that the switching frequency and load variation rate are sufficiently slow. For arbitrarily fast variations of the load, we use a level-set argument with multiple Lyapunov functions to establish input-to-state stability of a larger set and with respect to the rate of change of the loads. The theoretical results are illustrated via numerical simulations on the IEEE 39-bus test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Stable Pentagonal Layered Palladium Diselenide Enables Rapid Electrosynthesis of Hydrogen Peroxide

Electrosynthesis of hydrogen peroxide (H 2 O 2 ) via the two-electron oxygen reduction reaction (2e – ORR) is promising for various practical applications, such as wastewater treatment. However, few electrocatalysts are active and selective for 2e – ORR yet are also resistant to catalyst leaching under realistic operating conditions. Here, a joint experimental and computational study reveals active and stable 2e – ORR catalysis in neutral media over layered PdSe 2 with a unique pentagonal puckered ring structure type. Computations predict active and selective 2e – ORR on the basal plane and edge of PdSe 2 , but with distinct kinetic behaviors. Further, electrochemical measurements of hydrothermally synthesized PdSe 2 nanoplates show a higher 2e – ORR activity than other Pd–Se compounds (Pd 4 Se and Pd 17 Se 15 ). PdSe 2 on a gas diffusion electrode can rapidly accumulate H 2 O 2 in buffered neutral solution under a high current density. The electrochemical stability of PdSe 2 is further confirmed by long device operational stability, elemental analysis of the catalyst and electrolyte, and synchrotron X-ray absorption spectroscopy. This work establishes a new efficient and stable 2e – ORR catalyst at practical current densities and opens catalyst designs utilizing the unique layered pentagonal structure motif.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nonlinear error field response in the presence of plasma rotation and real frequencies due to favorable curvature

This works describes the interaction of an phenomenon known to arise in tokamak fusion experiments with the magnetic perturbations that can arise at the edge of such experiments. We present nonlinear NIMROD resistive MHD simulations of the response of a rotating plasma to an error field when the plasma has weakly damped linear tearing modes (TMs), stabilized by a pressure gradient and favorable curvature. The favorable curvature leads to the Glasser effect: the occurrence of real frequencies and stabilization with positive stability index Δ'.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A comprehensive conductivity model for drift and micro-tearing modes

The parallel electrical conductivity is a crucial parameter in the study of the linear stability of drift-modes like the micro-tearing mode (MTM). The conductivity enters by closing the electromagnetic tearing layer equations. Recent progress in the understanding of the pedestal suggests that the MTM could play an important role in its structure and evolution. For this reason, we revisit and improve previous model conductivities. This parameter is generally derived from the linearized drift kinetic equation. In the past literature, it has been computed using either simplified collision operators or neglecting the spatial dependence away from the rational surface. Here, a fully consistent expression for the conductivity that would accurately model the pedestal has not been available. By applying a novel variational procedure and with the full Fokker Plank collision operator, including electron–electron collisions, we compute a closed expression for the parallel electrical conductivity in the form of a rational function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Consistent coupling algorithms for coupled core-edge simulations of plasma turbulence

Two schemes for coupling gyrokinetic simulations of microturbulence in tokamaks are proposed. The first scheme is based on an additive Schwarz domain decomposition. In this work, we show that, because the goal of turbulence is long-time averages of the dynamics rather than temporal accuracy, the iteration to self-consistency across domains, which is typically required by Schwarz schemes, can be avoided, thereby accelerating the computation. Second, we propose a coupling scheme that relies entirely on the addition of source terms, leaving the boundary conditions arbitrary. The practical motivations for such a scheme are discussed, and forms of the source terms that ensure consistency and stability are derived. The schemes are tested on a nonlinear, one-dimensional model problem, and the first scheme is further tested on the Hasegawa–Wakatani model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Microtearing instabilities and electron thermal transport in low and high collisionality NSTX discharges

Microtearing mode (MTM) real frequency, growth rate, magnetic fluctuation amplitude and resulting electron thermal transport are studied in systematic NSTX scans of relevant plasma parameters. Here, the dependency of the MTM real frequency and growth rate on plasma parameters, suitable for low and high collision NSTX discharges, is obtained by using the reduced MTM transport model. The plasma parameter dependencies are compared and found to be consistent with the results obtained from MTM using the Gyrokinetic GYRO code. The scaling trend of collision frequency and plasma beta is found to be consistent with the global energy confinement trend observed in the NSTX experiment. The strength of the magnetic fluctuation is found to be consistent with the gyrokinetic estimate. In earlier studies, it was found that the version of the Multi-Mode (MM) anomalous transport model, which did not contain the effect of MTMs, provided an appropriate description of the electron temperature profiles in standard tokamak discharges and not in spherical tokamaks. When the MM model, which involves transport associated with MTMs, is incorporated in the TRANSP code and is used in the study of electron thermal transport in NSTX discharges, it is observed that the agreement with the experimental electron temperature profile is substantially improved.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗