Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Predictions of core plasma performance for the SPARC tokamak

SPARC is designed to be a high-field, medium-size tokamak aimed at achieving net energy gain with ion cyclotron range-of-frequencies (ICRF) as its primary auxiliary heating mechanism. Empirical predictions with conservative physics indicate that SPARC baseline plasmas would reach $\textit{Q}$ ≈ 11, which is well above its mission objective of $\textit{Q}$ > 2. To build confidence that SPARC will be successful, physics-based integrated modelling has also been performed. The TRANSP code coupled with the theory-based trapped gyro-Landau fluid (TGLF) turbulence model and EPED predictions for pedestal stability find that $\textit{Q}$ ≈ 9 is attainable in standard H-mode operation and confirms $\textit{Q}$ > 2 operation is feasible even with adverse assumptions. In this analysis, ion cyclotron waves are simulated with the full wave TORIC code and alpha heating is modelled with the Monte–Carlo fast ion NUBEAM module. Detailed analysis of expected turbulence regimes with linear and nonlinear CGYRO simulations is also presented, demonstrating that profile predictions with the TGLF reduced model are in reasonable agreement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Effect of polycyclosilane microstructure on thermal properties

Thermal characterization of polysilanes has focused on the influence of organic side chains, whereas little is understood about the influence of silane backbone microstructure on thermal stability, phase properties, and pyrolysis. To address this knowledge gap, we prepared three distinct polycyclosilanes: linear polymers synthesized from the cyclosilane building blocks 1,4Si 6 and 1,3Si 6 , as well as a cyclic polymer of 1,3Si 6 . Thermal properties across the temperature range 25 to 600 °C were investigated using differential scanning calromietry (DSC) and thermogravimetric analysis (TGA). We found differences between linear and cyclic materials, including a phase transition unique to the cyclic polymer and lower rates of mass loss during pyrolysis. Density functional theory (DFT) calculations provided insight into microstructure-dependent pyrolysis.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Estimating and Calibrating DER Model Parameters Using Levenberg–Marquardt Algorithm in Renewable Rich Power Grid

The proliferation of inverter-based distributed energy resources (IBDERs) has increased the number of control variables and dynamic interactions, leading to new grid control challenges. For stability analysis and designing appropriate protection controls, it is important that IBDER models are accurate. This paper focuses on the accurate estimation and parameter calibration of DER_A, a recently proposed aggregated IBDER model. In particular, we focus on the parameters of the reactive power–voltage regulation module. We formulate the problem of parameter tuning as a non-linear least square minimization problem and solve it using the Levenberg–Marquardt (LM) method. The LM method is primarily chosen due to its flexibility in adaptively selecting between the steepest descent and Gauss–Newton methods through a damping parameter. The LM approach is used to minimize the error between the actual measurements and the estimated response of the model. Further, the computational challenges posed by the numerical calculation of the Jacobian are tackled using a quasi-Newton root-finding approach. The proposed method is validated on a real feeder model in the northeastern part of the United States. The feeder is modeled in OpenDSS and the measurements thus obtained are fed to the DER_A model for calibration. The simulation results indicate that our approach is able to successfully calibrate the relevant model parameters quickly and with high accuracy, with a total sum of square error of 3.57 × 10 –7 .

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Polymorphism and Negative Linear Compressibility in Pyrazine-d4

Pyrazine (1,4-diazine) is a widely used linker in coordination polymers, molecular magnets, and metal–organic frameworks, yet its solid-state polymorphism is controlled by weak C–H···N interactions. Using variable-temperature and high-pressure neutron diffraction, a previously unrecognized low-temperature transition to phase IV is identified, showing that the form stabilized under pressure can also be accessed below ∼90 K at ambient pressure in the perdeuterated compound. The transition from phase III → IV involves a symmetry-breaking lattice distortion driven by the reorganization of C–H···N hydrogen bonds, accompanied by negative linear compressibility, which is primarily accommodated through the reduction in void-space within the structure. Structural analysis reveals that while the overall monoclinic framework of phase IV is preserved under both temperature and pressure routes, the evolution of hydrogen-bonded network is distinct depending on the pathway taken. These results establish pyrazine as a model molecular solid for studying pressure–temperature polymorphism and demonstrate how weak directional interactions can produce anomalous mechanical responses in organic crystals.

Ridley, Christopher J. [Oak Ridge National Laborat↗
DOE OSTI2025

Navigating thermal stability intricacies of high-nickel cathodes for high-energy lithium batteries

High-nickel oxide cathodes, LiNi x M 1−x O 2 (x ≥ 0.8), are preferred in automotive lithium batteries, but they face thermal instability challenges. Inconsistent literature reports and unstandardized testing protocols further complicate quantitative assessments of the thermal stability of these cathodes. Here, we present here a statistical thermal analysis based on the differential scanning calorimetry measurements of 15 representative cathode materials with different compositions, morphologies and states of charge. The findings reveal that each cathode has a critical state of charge that defines its safe operating limit, which is affected by the metal–oxygen bond strength and surface reactivity. The thermal runaway temperature is dictated by the layered Li 1−x NiO 2 to LiNi 2 O 4 spinel-like phase transition, which is thermodynamically determined by the metal–oxygen bond covalency and kinetically influenced by the cation mixing and particle size. Raman spectroscopy is used to predict the thermal runaway temperature on the basis of the linear relationship between them. Finally, we propose a thermal stability index to quantify cathode thermal stability as a guide for developing safer high-nickel cathodes.

batteries↗
DOE OSTI2021

Unveiling decaying mechanism through quantitative structure-activity relationship in electrolytes for lithium-ion batteries

Understanding the decaying mechanism in lithium-ion batteries (LIBs) is critical to establishing a stable electrolyte system. Despite the advent of various novel electrolyte solvents designated for high-voltage LIBs, their working principles are not fully understood. Currently, oxidative decomposition of electrolytes is believed to be the major cause of capacity fade, and tremendous effort has been devoted to discovering a new electrolyte with enhanced anodic stability. However, the oxidative decomposition process cannot solely explain the rapid decay of some electrolyte systems with intrinsic high anodic stability when used with a high-nickel layered oxide cathode such as LiNi 0.6 Mn 0.2 Co 0.2 O 2 (NMC622). In this report, a study of the quantitative structure-activity relationship was conducted to deepen the mechanistic understanding of the decay in high-voltage LIBs. Here, the results obtained from the newly introduced molecular pair analysis and linear free-energy relationship (LFER) studies were highly consistent with the solvation-involved decaying mechanism in a high-nickel layered oxide cathode cycling at high voltage (> 4.5 V vs. Li/Li + ). There was no evidence correlating the solvation ability of electrolyte solvents with the decay of a high-nickel layered oxide cathode cycling at a relatively low voltage (< 4.3 V vs. Li/Li + ), nor with the high-voltage spinel cathode LiNi 0.5 Mn 1.5 O 4 (LNMO). Undoubtedly, the unveiled mechanistic insight provides a critical guideline for the development of an appropriate electrolyte system targeting different high-voltage cathode materials.

25 ENERGY STORAGE↗
DOE OSTI2024

Stability of the first-order character of phase transition in HoCo 2

HoCo 2 exhibits a giant magnetocaloric (MC) effect at its first-order magnetostructural phase transition around 77 K, and understanding the thermodynamic nature of this transition in response to external magnetic fields is crucial for its MC applications. In this study, we present a comprehensive investigation of specific heat and magnetization measurements of HoCo 2 under varying magnetic fields. The specific heat measurements qualitatively indicate a transformation from first- to second-order behavior of this phase transition at higher magnetic fields. However, analysis of the power-law dependence of the magnetic entropy change (Δ S M ∝H n ) and the breakdown of universal behavior in the temperature dependence of Δ S M suggest that the first-order nature remains intact, even up to 7 T. This stability of the first-order nature is further manifested through the distinctive non-linear behavior of modified Arrott plots, with a negative slope in the 6–7 T range.

36 MATERIALS SCIENCE↗
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↗
DOE OSTI2020

H-mode pedestal improvements with neon injection in DIII-D

Improvements in the H-mode pedestal pressure and global energy confinement were determined on DIII-D with both 3 and 6 MW of neutral beam (NBI) heating with neon injection. The pedestal pressure increased, primarily associated with a significant increase in the pedestal density with a relatively smaller decrease in the pedestal temperature. At the same neon injection rate, a rapid decrease in edge localized mode (ELM) frequency and ~35% increase in plasma stored energy were observed in the 3 MW discharge, while in the 6 MW discharge, the ELM frequency was gradually decreased and the plasma stored energy gradually increased by up to 20%. The measured pedestal widths from both 3 and 6 MW discharges matched the EPED1.0 width scaling, ${\Delta _{{\psi _N}}} = 0.089\beta _{p,ped}^{1/2}$ ,data-, within 20% deviations before and after neon injection. The peeling-ballooning mode (PBM) stability analysis showed that after neon injection, the ballooning stability boundary was increased, while the peeling stability boundary only dropped once the neon level was relatively high. The ballooning stability boundary increase with neon injection is consistent with increase of the ion diamagnetic stabilization (i.e. including the contribution from carbon, neon as well as deuterium ions) which is used in the criterion to determine when the calculated linear growth rates of PBM are unstable. Consequently, the new ELITE computed PBM boundaries were more consistent with the operating points than when only deuterium was included. This in turn suggests that the impurity may help to improve pedestal pressure through affecting ion diamagnetic frequency in the pedestal region, which could positively lead to the ballooning stability boundary extension.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Visual Tool for Assessing Stability of DER Configurations on Three-Phase Radial Networks

Here, we present a method and tool for evaluating the placement of Distributed Energy Resources (DER) on distribution circuits in order to control voltages and power flows. Our previous work described Phasor-Based Control (PBC), a novel control framework where DERs inject real and reactive power to track voltage magnitude and phase angle targets. Here, we employ linearized power flow equations and integral controllers to develop a linear state space model for PBC acting on a three-phase unbalanced network. We use this model to evaluate whether a given inverter-based DER configuration admits a stable set of controller gains, which cannot be done by analyzing controllability nor by using the Lyapunov equation. Instead, we sample over a parameter space to identify a stable set of controller gains. Our stability analysis requires only a line impedance model and does not entail simulating the system or solving an optimization problem. We incorporate this assessment into a publicly available visualization tool and demonstrate three processes for evaluating many control configurations on the IEEE 123-node test feeder (123NF).

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Asymptotically compatible reproducing kernel collocation and meshfree integration for the peridynamic Navier equation

Here, we study reproducing kernel (RK) collocation method for peridynamic Navier equation. In the first part, we apply a linear RK approximation to both displacement and dilatation, and then back-substitute dilatation and solve the peridynamic Navier equation in a pure displacement form. The RK collocation scheme converges to the nonlocal limit for a fixed nonlocal interaction length and also to the local limit as nonlocal interactions vanish. The stability is shown by comparing the collocation scheme with the standard Galerkin scheme using Fourier analysis. In the second part, we apply the RK collocation to the quasi-discrete peridynamic Navier equation and show its convergence to the correct local limit when the ratio between the nonlocal length scale and the discretization parameter is fixed. The analysis is carried out on a special family of rectilinear Cartesian grids for the RK collocation method with a designated kernel with finite support. We assume the Lamé parameters satisfy λ ≥ μ to avoid extra assumptions on the nonlocal kernel. Finally, numerical experiments are conducted to validate the theoretical results.

42 ENGINEERING↗
DOE OSTI2020

Linear beam stability in periodic focusing systems: Krein signature and band structure

The general question of how a beam becomes unstable has been one of the fundamental research topics among beam and accelerator physicists for several decades. In this study, we revisited the general problem of linear beam stability in periodic focusing systems by applying the concepts of Krein signature and band structure. We numerically calculated the eigenvalues and other associated characteristics of one-period maps, and discussed the stability properties of single-particle motions with skew quadrupoles and envelope perturbations in high-intensity beams on an equal footing. In particular, an application of the Krein theory to envelope instability analysis was newly attempted in this study. Finally, the appearance of instabilities is interpreted as the result of the collision between eigenmodes of opposite Krein signatures and the formation of a band gap.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Computation of linear MHD instabilities with the multi-region relaxed MHD energy principle

Over the last decade, a variational principle based on a generalisation of Taylor's relaxation, referred to as multi-region relaxed magnetohydrodynamics (MRxMHDs) has been developed. The numerical solutions of the MRxMHD equilibria have been constructed using the Stepped Pressure Equilibrium Code (SPEC) (Hudson et al 2012 Phys. Plasmas 19 112502). In principle, SPEC could also be established to describe the MRxMHD stability of an equilibrium. In this work, a theoretical framework is developed to relate the second variation of the energy functional to the so-called Hessian matrix, enabling the prediction of MHD linear instabilities of cylindrical plasmas, and is implemented in SPEC. The negative and positive eigenvalues of the Hessian matrix predict the stability of an equilibrium. Verification studies of SPEC stability results with the M3D-C 1 code and the tearing mode $\Delta^{^{\prime}}$ criterion have been conducted for ideal and resistive MHD instabilities, respectively, in a pressureless cylindrical tokamak, and have shown good agreement. Overall, our stability analysis is a critical step towards understanding the MHD stability of three-dimensional MHDs where nested flux surfaces, magnetic islands and stochastic regions co-exist.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Including the parallel mass flow in calculating the steady-state solutions and stability of the momentum balance equations for a quasisymmetric stellarator

The Helically Symmetric Experiment (HSX) is a quasisymmetric stellarator with minimal parallel viscous damping in a helical direction. The parallel flow (Vǁ) along the magnetic field is similarly weakly damped by viscosity. In this paper, the self-consistent steady-state parallel and poloidal momentum balance equations are used to show that a large Vǁ on the order of the ion thermal velocity can increase the ion resonant radial electric field (Er) beyond the value calculated using the typical approximation that Vǁ is zero. By altering the damping of Vǁ, either by degrading the quasisymmetry or varying the neutral density, the ion resonant Er can shift in a controllable fashion. It is shown explicitly that there exist stable and unstable steady-state solutions in the two-dimensional space of Vǁ and Er. A stability analysis of each solution is performed by calculating the eigenvalues and eigenvectors of the Jacobian. The unstable solution corresponds to a saddle point in which the eigenvalues have opposite signs. The analysis leads to the conclusion that unstable solutions occur when the derivative of the total poloidal damping with respect to Er is positive. A hysteresis in Er and Vǁ is observed when the radial current density is linearly increased to a maximum and then decreased back to zero. Jumps in the radial electric field and the parallel flow are observed as the radial current density drives the evolution from one stable point to the next. This result is similar to experimental data observed on several devices.

Physics↗
DOE OSTI2023

Extended isogeometric analysis of multi-material and multi-physics problems using hierarchical B-splines

Here, this paper presents an immersed, isogeometric finite element framework to predict the response of multi-material, multi-physics problems with complex geometries using locally refined discretizations. To circumvent the need to generate conformal meshes, this work uses an extended finite element method (XFEM) to discretize the governing equations on non-conforming, embedding meshes. A flexible approach to create truncated hierarchical B-splines discretizations is presented. This approach enables the refinement of each state variable field individually to meet field-specific accuracy requirements. To obtain an immersed geometry representation that is consistent across all hierarchically refined B-spline discretizations, the geometry is immersed into a single mesh, the XFEM background mesh, which is constructed from the union of all hierarchical B-spline meshes. An extraction operator is introduced to represent the truncated hierarchical B-spline bases in terms of Lagrange shape functions on the XFEM background mesh without loss of accuracy. The truncated hierarchical B-spline bases are enriched using a generalized Heaviside enrichment strategy to accommodate small geometric features and multi-material problems. The governing equations are augmented by a formulation of the face-oriented ghost stabilization enhanced for locally refined B-spline bases. We present examples for two- and three-dimensional linear elastic and thermo-elastic problems. The numerical results validate the accuracy of our framework. The results also demonstrate the applicability of the proposed framework to large, geometrically complex problems.

42 ENGINEERING↗
DOE OSTI2021

Covalent Functionalization of Nickel Phosphide Nanocrystals with Aryl-Diazonium Salts

Covalent functionalization of Ni 2 P nanocrystals was demonstrated using aryl-diazonium salts. Spontaneous adsorption of aryl functional groups was observed, with surface coverages ranging from 20 to 96% depending on the native reactivity of the salt as determined by the aryl substitution pattern. Increased coverage was possible for low reactivity species using a sacrificial reductant. Functionalization was confirmed using thermogravimetric analysis, Fourier transform infrared spectroscopy, and X-ray photoelectron spectroscopy. The structure and energetics of this nanocrystal electrocatalyst system, as a function of ligand coverage, were explored with density functional theory calculations. The Hammett parameter of the surface functional group was found to linearly correlate with the change in Ni and P core-electron binding energies and the nanocrystal's experimentally and computationally determined work function. The electrocatalytic activity and stability of the functionalized nanocrystals for hydrogen evolution were also improved when compared to the unfunctionalized material, but a simple trend based on electrostatics was not evident. Density functional theory was used to understand this discrepancy, revealing that H adsorption energies on the covalently functionalized Ni 2 P also do not follow the electrostatic trend and are predictive descriptors of the experimental results.

14 SOLAR ENERGY↗
DOE OSTI2020

Uncovering the Relationship between Diameter and Height of Electrodeposited Lithium Protrusions in a Rigid Electrolyte

A promising approach for enabling rechargeable batteries with significantly higher energy densities than current lithium-ion batteries is by deploying lithium-metal anodes. However, the growth of lithium protrusions during charging presents significant challenges. Since these protrusions are often branched and filamentous in conventional liquid electrolytes, this problem is referred to in the literature as the “dendrite problem”. While solid electrolytes have the potential to solve this problem, protrusions grow in all electrolytes when the current density exceeds a critical value. Fundamentally understanding the formation is necessary to develop a rational approach for increasing the critical current density, but it is challenging due to the complex interplay between electrochemical and material properties. The diameters and heights of protrusions on lithium-metal anodes stabilized by a rigid block copolymer electrolyte were measured in situ by synchrotron hard X-ray microtomography. The diameter of the shorting protrusions increased linearly with increasing electrolyte thickness. Further, a universal linear relationship between protrusion height and diameter of both shorting and non-shorting protrusions was observed. A model based on the concentrated solution theory was used to establish the electrochemical and mechanical sources for our observations. Here, the computational analysis indicates that elastic and plastic deformation of both the lithium metal and the polymer are important to describe protrusion growth. Both stress-induced current density effects due to the deformation of the electrolyte near the protrusion and plastic deformation of lithium metal combine to give the counterintuitive result: the fastest-growing protrusions have the largest diameter.

25 ENERGY STORAGE↗
DOE OSTI2022

Analysis of cell-based diffusion acceleration for the slice balance approach

In this work, we perform analysis on the use of cell-based diffusion acceleration methodologies to accelerate the convergence of transport solutions discretized with the slice balance approach (SBA) on unstructured polygonal grids.We investigated both linear diffusion synthetic acceleration (DSA) and non linear diffusion acceleration (NDA), including its partial-current variant (pNDA). DSA and NDA were both shown to diverge for intermediate ranges of mesh optical thicknesses. However, pNDA and Krylov methods like GMRES and Broyden stabilized the acceleration schemes, including problems with degenerate cells formed by mesh refinement. (author)

42 ENGINEERING↗