A MnO[subscript x] enhanced atomically dispersed iron?nitrogen?carbon catalyst for the oxygen reduct
Explore the source record for details and available documents.
Ding, Shichao↗
SEARCH · Search NASA
Results for “stabilization method”
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.
On the formulation of a minimal uncertainty model for robust control with structured uncertainty
In the design and analysis of robust control systems for uncertain plants, representing the system transfer matrix in the form of what has come to be termed an M-delta model has become widely accepted and applied in the robust control literature. The M represents a transfer function matrix M(s) of the nominal closed loop system, and the delta represents an uncertainty matrix acting on M(s). The nominal closed loop system M(s) results from closing the feedback control system, K(s), around a nominal plant interconnection structure P(s). The uncertainty can arise from various sources, such as structured uncertainty from parameter variations or multiple unsaturated uncertainties from unmodeled dynamics and other neglected phenomena. In general, delta is a block diagonal matrix, but for real parameter variations delta is a diagonal matrix of real elements. Conceptually, the M-delta structure can always be formed for any linear interconnection of inputs, outputs, transfer functions, parameter variations, and perturbations. However, very little of the currently available literature addresses computational methods for obtaining this structure, and none of this literature addresses a general methodology for obtaining a minimal M-delta model for a wide class of uncertainty, where the term minimal refers to the dimension of the delta matrix. Since having a minimally dimensioned delta matrix would improve the efficiency of structured singular value (or multivariable stability margin) computations, a method of obtaining a minimal M-delta would be useful. Hence, a method of obtaining the interconnection system P(s) is required. A generalized procedure for obtaining a minimal P-delta structure for systems with real parameter variations is presented. Using this model, the minimal M-delta model can then be easily obtained by closing the feedback loop. The procedure involves representing the system in a cascade-form state-space realization, determining the minimal uncertainty matrix, delta, and constructing the state-space representation of P(s). Three examples are presented to illustrate the procedure.
Belcastro, Christine M.↗
Extended virtual element method for the Laplace problem with singularities and discontinuities
In this paper, we propose the extended virtual element method (X-VEM) to treat singularities and crack discontinuities that arise in the Laplace problem. The virtual element method (VEM) is a stabilized Galerkin formulation on arbitrary polytopal meshes, wherein the basis functions are implicit (virtual)—they are not known explicitly nor do they need to be computed within the problem domain. In this work, suitable projection operators are used to decompose the bilinear form on each element into two parts: a consistent term that reproduces the first-order polynomial space and a correction term that ensures stability. A similar approach is pursued in the X-VEM with a few notable extensions. To capture singularities and discontinuities in the discrete space, we augment the standard virtual element space with an additional contribution that consists of the product of virtual nodal basis (partition-of-unity) functions with enrichment functions. For discontinuities, basis functions are discontinuous across the crack and for singularities a weakly singular enrichment function that satisfies the Laplace equation is chosen. For the Laplace problem with a singularity, we devise an extended projector that maps functions that lie in the extended virtual element space onto linear polynomials and the enrichment function, whereas for the discontinuous problem, the consistent element stiffness matrix has a block-structure that is readily computed. An adaptive homogeneous numerical integration method is used to accurately and efficiently (no element-partitioning is required) compute integrals with integrands that are weakly singular. Once the element projection matrix is computed, the same steps as in the standard VEM are followed to compute the element stabilization matrix. Numerical experiments are performed on quadrilateral and polygonal (convex and nonconvex elements) meshes for the problem of an $\textit{L}$-shaped domain with a corner singularity and the problem of a cracked membrane under mode III loading, and results are presented that affirm the sound accuracy and demonstrate the optimal rates of convergence in the $\textit{L}^2$ norm and energy of the proposed method.
42 ENGINEERING↗
CFL Optimized Forward–Backward Runge–Kutta Schemes for the Shallow-Water Equations
Abstract We present the formulation and optimization of a Runge–Kutta-type time-stepping scheme for solving the shallow-water equations, aimed at substantially increasing the effective allowable time step over that of comparable methods. This scheme, called FB-RK(3,2), uses weighted forward–backward averaging of thickness data to advance the momentum equation. The weights for this averaging are chosen with an optimization process that employs a von Neumann–type analysis, ensuring that the weights maximize the admittable Courant number. Through a simplified local truncation error analysis and numerical experiments, we show that the method is at least second-order in time for any choice of weights and exhibits low dispersion and dissipation errors for well-resolved waves. Further, we show that an optimized FB-RK(3,2) can take time steps up to 2.8 times as large as a popular three-stage, third-order strong stability-preserving Runge–Kutta method in a quasi-linear test case. In fully nonlinear shallow-water test cases relevant to oceanic and atmospheric flows, FB-RK(3,2) outperforms SSPRK3 in admittable time step by factors roughly between 1.6 and 2.2, making the scheme approximately twice as computationally efficient with little to no effect on solution quality. Significance Statement The purpose of this work is to develop and optimize time-stepping schemes for models relevant to oceanic and atmospheric flows. Specifically, for the shallow-water equations we optimize for schemes that can take time steps as large as possible while retaining solution quality. We find that our optimized schemes can take time steps between 1.6 and 2.2 times larger than schemes that cost the same number of floating point operations, translating directly to a corresponding speedup. Our ultimate goal is to use these schemes in climate-scale simulations.
54 ENVIRONMENTAL SCIENCES↗
Porosity evolution during heating of copper made from powder by friction extrusion
Here, friction extrusion was used to compact and extrude solid copper rod from feedstock powders. Following extrusion, considerable porosity was observed throughout the extrudate cross section due to the entrained porosity in the feedstock material and the extrusion process. The thermal stability of the extrudate was investigated via a series of heat treatments. Porosity evolution exhibited three distinct stages—an unchanged plateau (0°C–300°C) followed by an abrupt increase (400°C–500°C) and ultimately a reduction, as the temperature increases (>500°C). The peak porosity measured was ~25%. The underlying driving force for pore evolution is described as the competition between the internal pore pressure, material strength, and sintering kinetics, as a function of temperature. The observed porosity evolution and driving force are not expected to be limited to copper. Thus, this manuscript reveals an important consideration regarding elevated-temperature thermal stability as advanced manufacturing methods involving direct powder extrusion are explored.
36 MATERIALS SCIENCE↗
Measurement of Stability Diagram at IOTA at Fermilab
Nonlinear focusing elements can enhance the stability of particle beams in high-energy colliders by means of Landau Damping, through the tune spread which is introduced. We propose an experiment at Fermilab’s Integrable Optics Test Accelerator (IOTA) to investigate the influence of nonlinear focusing elements on the transverse stability of the beam. In this experiment, we employ an anti-damper, an active transverse feedback system, as a controlled mechanism to induce coherent beam instability. By utilizing the anti-damper, we can examine the impact of the nonlinear focusing element on the beam’s transverse stability. The stability diagram, a tool used to determine the system’s stability, will be measured using a recently demonstrated method at the LHC. This measurement is carried out experimentally by selecting specific threshold gains and measuring them for a range of phases. The stability diagram is represented by gei¿ on the complex plane. The experiment at IOTA adds insight towards the stability diagram measurement method by supplying a reduced machine impedance, to investigate the impedance’s effect on the stability diagram, as well as a larger range of phase measurements.
43 PARTICLE ACCELERATORS↗
Efficient IMEX Runge-Kutta Methods for Nonhydrostatic Dynamics
We develop new implicit-explicit Runge-Kutta (IMEX RK) methods for integrating horizontally explicit, vertically implicit (HEVI) partitionings of nonhydrostatic atmosphere models (HEVI models). These new methods, termed IMKG methods, are IMEX RK methods whose explicit part has optimal or near-optimal stability on the imaginary axis and whose implicit part is I-stable. A specialized stability region is presented for characterizing the stability of IMKG and other IMEX methods integrating HEVI models. Subsequently, we formulate two families of IMEX RK methods to enable deriving IMKG methods with a high explicit stage count for integrating HEVI models with large, stable time-steps. Here, we then derive a HEVI partitioning of the HOMME-NH nonhydrostatic model and use this model to compare the accuracy and efficiency of several IMKG methods with other IMEX RK methods from the literature.
58 GEOSCIENCES↗
A Cyber Threat Mitigation Approach For Wide Area Control of SVCs using Stability Monitoring
We propose a stability monitoring approach for the mitigation of cyber threats directed at the wide area control (WAC) system used for coordinated control of Flexible AC Transmission Systems (FACTS) used for power oscillation damping (POD) of active power flow on inter-area tie lines. The approach involves monitoring the modes of the active power oscillation on an inter-area tie line using the Matrix Pencil (MP) method. We use the stability characteristics of the observed modes as a proxy for the presence of destabilizing cyber threats. We monitor the system modes to determine whether any destabilizing modes appear after the WAC system engages to control the POD. If the WAC signal exacerbates the POD performance, the FACTS falls back to POD using local measurements. The proposed approach does not require an expansive system-wide view of the network. We simulate replay, control integrity, and timing attacks for a test system and present results that demonstrate the performance of the SM approach for mitigation.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Impact of Aerodynamics and Structures Technology on Heavy Lift Tiltrotors
Rotor performance and aeroelastic stability are presented for a 124,000-lb Large Civil Tilt Rotor (LCTR) design. It was designed to carry 120 passengers for 1200 nm, with performance of 350 knots at 30,000 ft altitude. Design features include a low-mounted wing and hingeless rotors, with a very low cruise tip speed of 350 ft/sec. The rotor and wing design processes are described, including rotor optimization methods and wing/rotor aeroelastic stability analyses. New rotor airfoils were designed specifically for the LCTR; the resulting performance improvements are compared to current technology airfoils. Twist, taper and precone optimization are presented, along with the effects of blade flexibility on performance. A new wing airfoil was designed and a composite structure was developed to meet the wing load requirements for certification. Predictions of aeroelastic stability are presented for the optimized rotor and wing, along with summaries of the effects of rotor design parameters on stability.
Acree, C. W., Jr.↗
Nanoparticle Assembly in High Polymer Concentration Solutions Increases Superlattice Stability
Explore the source record for details and available documents.
Lee, Margaret S.↗
Large-signal Stability Analysis of Grid-forming Inverters with Equivalent-circuit Models
Here, this paper proposes an energy function-based direct method for large-signal stability assessment of grid-forming (GFM) inverters leveraging an equivalent-circuit representation of all involved control- and physical-layer dynamics. Three different primary controls, a standard inner-current outer-voltage cascaded-control architecture, output LCL filter, and reference-current saturation limiting are featured in the modeling and analysis framework. A composite energy function for the GFM inverter is obtained by summing up individual energy contributions gleaned from the circuit representation. The approach can readily be generalized to different primary controls, output-filter arrangements, and current limiters since it is based on a circuit-theoretic foundation. Numerical simulations validate the efficacy of the approach in estimating the critical clearing time following a large-signal disturbance.
24 POWER TRANSMISSION AND DISTRIBUTION↗
A globally efficient means of distributing UTC time and frequency through GPS
Time and frequency outputs comparable in quality to the best laboratories have been demonstrated on an integrated system suitable for field application on a global basis. The system measures the time difference between 1 pulse-per-second (pps) signals derived from local primary frequency standards and from a multi-channel GPS C/A receiver. The measured data is processed through optimal SA Filter algorithms that enhance both the stability and accuracy of GPS timing signals. Experiments were run simultaneously at four different sites. Even with large distances between sites, the overall results show a high degree of cross-correlation of the SA noise. With sufficiently long simultaneous measurement sequences, the data shows that determination of the difference in local frequency from an accepted remote standard to better than 1 x 10(exp -14) is possible. This method yields frequency accuracy, stability, and timing stability comparable to that obtained with more conventional common-view experiments. In addition, this approach provides UTC(USNO MC) in real time to an accuracy better than 20 ns without the problems normally associated with conventional common-view techniques. An experimental tracking loop was also set up to demonstrate the use of enhanced GPS for dissemination of UTC(USNO MC) over a wide geographic area. Properly disciplining a cesium standard with a multi-channel GPS receiver, with additional input from USNO, has been found to permit maintaining a timing precision of better than 10 ns between Palo Alto, CA and Washington, DC.
Kusters, John A.↗
Diagonally Implicit Runge-Kutta Methods for Ordinary Differential Equations. A Review
A review of diagonally implicit Runge-Kutta (DIRK) methods applied to rst-order ordinary di erential equations (ODEs) is undertaken. The goal of this review is to summarize the characteristics, assess the potential, and then design several nearly optimal, general purpose, DIRK-type methods. Over 20 important aspects of DIRKtype methods are reviewed. A design study is then conducted on DIRK-type methods having from two to seven implicit stages. From this, 15 schemes are selected for general purpose application. Testing of the 15 chosen methods is done on three singular perturbation problems. Based on the review of method characteristics, these methods focus on having a stage order of two, sti accuracy, L-stability, high quality embedded and dense-output methods, small magnitudes of the algebraic stability matrix eigenvalues, small values of aii, and small or vanishing values of the internal stability function for large eigenvalues of the Jacobian. Among the 15 new methods, ESDIRK4(3)6L[2]SA is recommended as a good default method for solving sti problems at moderate error tolerances.
Kennedy, Christopher A.↗
On strategies for inverting remote sensing data
Attention is given to a number of methods for inverting remote sensing data obtained in a variety of astronomical applications. Applications include image restoration, inversion of helioseismological data to obtain the internal rotation rate of stars such as the sun, fitting of spectra (especially thermal line spectra) to grating or other dispersed observed spectra, differential emission measure analysis, and reconstruction of images derived from interferometric observations. The results consider the tradeoff between resolution and variance and the stability properties for each method and propose an inversion stragegy using the available techniques.
Jeffrey, W.↗
Fourth Order Finite Difference Methods for the Wave Equation with Mesh Refinement Interfaces
In this work, we analyze two types of summation-by-parts finite difference operators for approximating the second derivative with variable coefficient. The first type uses ghost points, while the second type does not use any ghost points. A previously unexplored relation between the two types of summation-by-parts operators is investigated. By combining them we develop a new fourth order accurate finite difference discretization with hanging nodes on the mesh refinement interface. We take the model problem as the two-dimensional acoustic wave equation in second order form in terms of acoustic pressure, and we prove energy stability for the proposed method. Compared to previous approaches using ghost points, the proposed method leads to a smaller system of linear equations that needs to be solved for the ghost point values. Another attractive feature of the proposed method is that the explicit time step does not need to be reduced relative to the corresponding periodic problem. Numerical experiments, both for smoothly varying and discontinuous material properties, demonstrate that the proposed method converges to fourth order accuracy. A detailed comparison of the accuracy and the time-step restriction with the simultaneous-approximation-term penalty method is also presented. (An erratum is attached.)
97 MATHEMATICS AND COMPUTING↗
A Virtual Element Method for the Full MHD system
In this manuscript we present a novel discretization for the incompressible MHD system. Our approach follows the framework of the Virtual Element Method and offers two main advantages. The method can be implemented in unstructured meshes making it highly versatile and capable of handling a wide array of problems involving interfaces, free-boundaries or adaptive refinements on the mesh. The second advantage involves the divergence of the magnetic field, our approach guarantees that it remains solenoidal. We include a theoretical proof of the condition on the magnetic field as well as energy estimates and a well-posedness study. The latter sheds light as to the stability properties of the method.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Anisotropic Goal-Based Mesh Adaptation Metric Clarification and Development
Adaptive unstructured mesh techniques have a limited, but growing impact on production analysis workflows to control discretization error for reliable simulation results. Multiple independent implementations of flow solvers, anisotropic metric construction methods, and anisotropic mesh adaptation mechanics have matured. Goal-based metrics target estimated error in output functions, such as lift and drag, through the guidance of an adjoint solution. A unification of goal-based anisotropic metrics is presented for steady viscous flows, which is an active area of research. These goal-based metrics drive robust and efficient anisotropic mesh adaptation for the calculation of output functions. The super-convergent functional output error behavior of stabilized finite-element methods is exploited without a formal proof, and evidence of super-convergence is shown in numerical experiments. Mesh adapted drag and lift outputs for two simple bodies in compressible viscous flow show convergence of error to less than a single drag count. Asymptotic behavior established for relatively coarse meshes shows the efficiency of this goal-based metric when compared to solution interpolation error control and expert-guided meshing. Anisotropic mesh adaptation techniques are applied to a transport aircraft in a high-lift configuration where variation between approaches decreases with mesh refinement, but asymptotic behavior is not observed with available resources.
goal-based↗