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 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

GIST (Grid Impedance Scan Tool) [SWR-22-73]

NREL has developed GIST (Grid Impedance Scan Tool) to scan the impedance response of inverter-based resources such as land-based and offshore wind power plants, solar power plants, etc., from their electromagnetic transient simulation models to evaluate their impact on grid stability. GIST can scan impedance response of inverter-based resources from vendor-supplied black-box simulation models without requiring internal proprietary details, enabling detailed stability analysis of power systems with high levels of inverter-based generation.

Shah, Shahil↗
DOE OSTI2024

Control Algorithms for Multi-Wavefront Sensor Adaptive Optics for Astronomical Exoplanet Imaging

This package consists of simulations of control algorithms being proposed for upcoming adaptive optics projects for exoplanet direct imaging, in which one wavefront affector (deformable mirror) is controlled by two wavefront sensors. Simulations include time-domain behavior under various assumed conditions, parameter optimization routines for the controllers, and stability analysis.

Sengupta, Aditya↗
DOE OSTI2023

Semi-Implicit Computation of Fast Modes in a Scheme Integrating Slow Modes by a Leapfrog Method Based on a Selective Implicit Time Filter

Abstract A scheme for integration of atmospheric equations containing terms with differing time scales is developed. The method employs a filtered leapfrog scheme utilizing a fourth-order implicit time filter with one function evaluation per time step to compute slow-propagating phenomena such as advection and rotation. The terms involving fast-propagating modes are handled implicitly with an unconditionally stable method that permits application of larger time steps and faster computations compared to fully explicit treatment. Implementation using explicit and recurrent formulation is provided. Stability analysis demonstrates that the method is conditionally stable for any combination of frequencies involved in the slow and fast terms as they approach the origin. The implicit filter used in the method damps the computational modes without noticeably sacrificing the accuracy of the physical mode. TheO[(Δt 4 )] accuracy for amplitude errors achieved by the implicitly filtered leapfrog is preserved in applications where terms responsible for fast propagation are integrated with a semi-implicit method. Detailed formulation of the method for soundproof nonhydrostatic anelastic equations is provided. Procedures for implementation in global spectral shallow-water models are also given. Examples comparing numerical and analytical solutions for linear gravity waves demonstrate the accuracy of the scheme. The performance is also shown in more practical nonlinear applications, where numerical solutions accomplished by the method are evaluated against those computed from a scheme where the slow terms are handled by the third-order Runge–Kutta scheme. It demonstrates that the method is able to accurately resolve fine-scale dynamics of Kelvin–Helmholtz shear instabilities, the evolution of density current, and nonlinear drifts of twin tropical cyclones.

Meteorology & Atmospheric Sciences↗
DOE OSTI2020

Numerical Experiments and Identification of Areas for Further Consideration (B639388 Subcontract Quarter II Report)

The second quarter of the project was spent carrying out numerical experiments on a variety of test matrices, using a combination of precisions, with the intent to study the numerical properties (namely, accuracy and convergence behavior) of the Conjugate Gradient (CG) method. We summarize our findings in the remainder of the document. Other activities include attending biweekly xSDK meetings. We are also currently collaborating with Steven Thomas (NREL) on the numerical stability analysis of “low-synch” Gram-Schmidt routines, for potential application within high-performance GMRES variants. This work is in progress. The subsequent quarter will be spent delving into the numerical analysis in order to provide theoretical explanation for the behavior observed in our experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Adaptive Recovery Model: Designing Systems for Testing Tracing and Vaccination to Support COVID-19 Recovery Planning.

This report documents a new approach to designing disease control policies that allocate scarce testing, contact tracing, and vaccination resources to better control community transmission of COVID19 or similar diseases. The Adaptive Recovery Model (ARM) combines a deterministic compartmental disease model with a stochastic network disease propagation model to enable us to simulate COVID-19 community spread through the lens of two complementary modeling motifs. ARM contact networks are derived from cell-phone location data that have been anonymized and interpreted as individual arrivals to specic public locations. Modeling disease spread over these networks allows us to identify locations within communities conducive to rapid disease spread. ARM applies this model- and data-derived abstractions of community transmission to evaluate the effectiveness of disease control measures including targeted social distancing, contact tracing, testing and vaccination. The architecture of ARM provides a unique capacity to help decision makers understand how best to deploy scarce testing, tracing and vaccination resources to minimize disease-spread potential in a community. This document details the novel mathematical formulations underlying ARM, presents a dynamical stability analysis of the deterministic model components, a sensitivity analysis of control parameters and network structure, and summarizes a process for deriving contact networks from cell-phone location data. An example use case steps through applying ARM to evaluate three targeted social distancing policies using Bernalillo County, New Mexico as an exemplar test locale. This step-by-step analysis demonstrates how ARM can be used to measure the relative performance of competing public health policies. Initial scenario tests of ARM shows that ARMs design focus on resource utilization rather than simple incidence prediction can provide decision makers with additional quantitative guidance for managing ongoing public health emergencies and planning future responses.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

The DPG Method for the Convection-Reaction Problem Revisited

We study both conforming and non-conforming versions of the practical DPG method for the convection-reaction problem. We determine that the most common approach for DPG stability analysis (construction of a local Fortin operator) is infeasible for the convection-reaction problem. We then develop a line of argument based on the direct construction of a global Fortin operator; we find that employing a polynomial enrichment for the test space does not suffice for this purpose, motivating the introduction of a (two-element) subgrid mesh. The argument combines mathematical analysis with numerical experiments

42 ENGINEERING↗
DOE OSTI2022

Accelerated Wind-Turbine Wake Recovery Through Actuation of the Tip-Vortex Instability

Advances in wind-plant control have often focused on more effectively balancing power between neighboring turbines. Wake steering is one such method that provides control-based improvements in a quasi-static way, but this does little to fundamentally change the wake recovery process, and thus, it has limited potential. This study investigates use of another control paradigm known as dynamic wake control (DWC) to excite the mutual inductance instability between adjacent tip-vortex structures, thereby accelerating the breakdown of the structures. The current work carries this approach beyond the hypothetical by applying the excitation via turbine control vectors that already exist on all modern wind turbines: blade pitch and rotor speed control. The investigation leverages a free-vortex wake method (FVWM) that allows a thorough exploration of relevant frequencies and amplitudes of harmonic forcing for each control vector (as well as the phase difference between the vectors for a tandem configuration) while still capturing the essential tip-vortex dynamics. The FVWM output feeds into a Fourier stability analysis working to pinpoint candidate DWC strategies suggesting fastest wake recovery. Near-wake length reductions of >80% are demonstrated, although without considering inflow turbulence. Analysis is provided to interpret these predictions considering the presence of turbulence in a real atmospheric inflow.

17 WIND ENERGY↗
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 OSTI2020

Invisible QCD as Dark Energy

We account for the late time acceleration of the Universe by extending the Quantum Chromodynamics (QCD) color to a SU(3) invisible sector (IQCD). If the Invisible Chiral symmetry is broken in the early universe, a condensate of dark pions (dpions) and dark gluons (dgluons) forms. The condensate naturally forms due to strong dynamics similar to the Nambu–Jona-Lasinio mechanism. As the Universe evolves from early times to present times the interaction energy between the dgluon and dpion condensate dominates with a negative pressure equation of state and causes late time acceleration. We conclude with a stability analysis of the coupled perturbations of the dark pions and dark gluons.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Nonlocal Magneto-curvature Instability in a Differentially Rotating Disk

A global mode is shown to be unstable to nonaxisymmetric perturbations in a differentially rotating Keplerian disk containing either vertical or azimuthal magnetic fields. In an unstratified cylindrical disk model, using both global eigenvalue stability analysis and linear global initial-value simulations, it is demonstrated that this instability dominates at strong magnetic fields where local standard magnetorotational instability (MRI) becomes stable. Unlike the standard MRI mode, which is concentrated in the high flow shear region, these distinct global modes (with low azimuthal mode numbers) are extended in the global domain and are Alfvén-continuum-driven unstable modes. As its mode structure and relative dominance over MRI are inherently determined by the global spatial curvature as well as the flow shear in the presence of a magnetic field, we call it the magneto-curvature (magneto-spatial-curvature) instability. Consistent with the linear analysis, as the field strength is increased in the nonlinear simulations, a transition from MRI-driven turbulence to a state dominated by global nonaxisymmetric modes is obtained. This global instability could therefore be a source of nonlinear transport in accretion disks at a higher magnetic field than predicted by local models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Kinetic Simulations of the Kruskal–Schwarzschild Instability in Accelerating Striped Outflows: Dynamics and Energy Dissipation

Astrophysical relativistic outflows are launched as Poynting-flux dominated, yet the mechanism governing efficient magnetic dissipation, which powers the observed emission, is still poorly understood. We study magnetic energy dissipation in relativistic "striped" jets, which host current sheets separating magnetically dominated regions with opposite field polarity. The effective gravity force g in the rest frame of accelerating jets drives the Kruskal–Schwarzschild instability (KSI), a magnetic analog of the Rayleigh–Taylor instability. By means of 2D and 3D particle-in-cell simulations, we study the linear and nonlinear evolution of the KSI. The linear stage is well described by linear stability analysis. The nonlinear stages of the KSI generate thin (skin-depth-thick) current layers, with length comparable to the dominant KSI wavelength. There, the relativistic drift-kink mode and the tearing mode drive efficient magnetic dissipation. The dissipation rate can be cast as an increase in the effective width Δ eff of the dissipative region, which follows dΔ eff /dt ≃ 0.05 $\sqrt{Δ_{eff} g}$. Our results have important implications for the location of the dissipation region in gamma-ray burst and active galactic nuclei jets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Bistable dynamics and Hopf bifurcation in a refined model of early stage HIV infection

Recent clinical studies have shown that HIV disease pathogenesis can depend strongly on many factors at the time of transmission, including the strength of the initial viral load and the local availability of CD4+ T-cells. In this article, a new within-host model of HIV infection that incorporates the homeostatic proliferation of T-cells is formulated and analyzed. Due to the effects of this biological process, the influence of initial conditions on the proliferation of HIV infection is further elucidated. The identifiability of parameters within the model is investigated and a local stability analysis, which displays additional complexity in comparison to previous models, is conducted. The current study extends previous theoretical and computational work on the early stages of the disease and leads to interesting nonlinear dynamics, including a parameter region featuring bistability of infectious and viral clearance equilibria and the appearance of a Hopf bifurcation within biologically relevant parameter regimes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Integrated Reliability and Economic Modeling for Transmission Across Large Regions: A Space Odyssey

Power flow modeling and stability analysis are needed to more-comprehensively assess system reliability but the development of the system portfolios and conditions require use of economic models (e.g., production cost). What are the state of art methods for efficiently linking economic and reliability models to enable examination of multiple snapshots and perform detailed nodal analyses?

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Anomalous transport in the H-mode pedestal of Alcator C-Mod discharges

Anomalous transport in the H-mode pedestal region of five Alcator C-Mod discharges, representing a collisionality scan is analyzed. The understanding of anomalous transport in the pedestal region is important for the development of a comprehensive model for the H-mode pedestal slope. In this research, a possible role of the drift resistive inertial ballooning modes (Rafiq et al 2010 Phys. Plasmas 17 082511) in the edge of Alcator C-Mod discharges is analyzed. The stability analysis, carried out using the TRANSP code, indicates that the DRIBM modes are strongly unstable in Alcator C-Mod discharges with large electron collisionality. An improved interpretive analysis of H-mode pedestal experimental data is carried out utilizing the additive flux minimization technique (Pankin et al 2013 Phys. Plasmas 20 102501) together with the guiding-center neoclassical kinetic XGC0 code. The neoclassical and neutral physics are simulated in the XGC0 code and the anomalous fluxes are computed using the additive flux minimization technique. The anomalous fluxes are reconstructed and compared with each other for the collisionality scan Alcator C-Mod discharges. It is found that the electron thermal anomalous diffusivities at the pedestal top increase with the electron collisionality. This dependence can also point to the drift resistive inertial ballooning modes as the modes that drive the anomalous transport in the plasma edge of highly collisional discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1963

INVESTIGATION OF PROBLEMS IN THERMAL CONVECTION

Numerical solution of thermal convection equations for velocity and temperature fluctuations as well as the stability analysis and heat transport for free boundary conditions

THERMAL CONVECTION↗