SEARCH · Engineering Papers
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.
A method of forming a nanocomposite thermoelectric material having microstructural stability at temperatures greater than 1000 C. The method includes creating nanocrystalline powder by cryomilling. The method is particularly useful in forming SiGe alloy powder.
Active control of aeroelastic response is a complex problem in which the designer usually tries to satisfy many design criteria which are often conflicting in nature. To further complicate the design problem, the state space equations describing this type of control problem are usually of high order, involving a large number of states to represent the flexible structure and unsteady aerodynamics. Control laws based on the standard Linear - Quadratic - Gaussian method are of the same high order as the aeroelastic plant and may be difficult to implement in the flight computer. To overcome this disadvantage a new approach was developed for designing low-order optimized robust control laws. In this approach, a nonlinear programming algorithm is used to search for the values of control law design variables that minimize a performance index while satisfying several inequality constraints that describe the design criteria on the stability robustness and responses. The method is applied to a gust load alleviation problem and a stability robustness improvement problem of a drone aircraft.
Throughout history, two similar grid-forming droop controls are widely reported in literature—the single-and multi-loop droop controls. Although being very similar, the authors find that the dynamic performance and stability characteristics of each control method are very different in a microgrid. Compared with single-loop droop control, the multi-loop droop control is prone to be less damped and loses stability more easily under some circumstances. This paper provides a novel insight into the different dynamic responses of the two basic controls. It points out that the two similar controls adjust the angular frequency and voltage magnitude at different locations within the inverter, resulting in different coupling reactances that impact the dynamic response and stability of microgrids differently. The use of single-loop droop control results in a larger coupling reactance, which helps to improve dynamic response and stability. This novel insight is verified through full-order small signal analysis, offline electromagnetic transient simulation, and real-time hardware-in-the-loop simulation experiments. Results show that the microgrid has a larger small signal stability boundary when using single-loop droop control, and this difference increases as the value of an inverter’s inner filter inductance increases.
A stability analysis of a thin horizontal rotating fluid layer which is subjected to arbitrary horizontal and vertical temperature gradients is presented. The basic state is a nonlinear Hadley cell which contains both Ekman and thermal boundary layers; it is given in closed form. The stability analysis is based on the linearized Navier-Stokes equations, and zonally symmetric perturbations in the form of waves propagating in the meridional direction are considered. Numerical methods were used for the stability problem. It was found that the instability sets in when the Richardson number is close to unity and that the critical Richardson number is a non-monotonic function of the Prandtl number. Further, it was found that the critical Richardson number decreases with increasing Ekman number until a critical value of the Ekman number is reached beyond which the fluid is stable.
Abstract. Phase correlation (PC) is a well-known method for estimating cloud motion vectors (CMVs) from infrared and visible spectrum images. Commonly, phase shift is computed in the small blocks of the images using the fast Fourier transform. In this study, we investigate the performance and the stability of the blockwise PC method by changing the block size, the frame interval, and combinations of red, green, and blue (RGB) channels from the total sky imager (TSI) at the United States Atmospheric Radiation Measurement user facility's Southern Great Plains site. We find that shorter frame intervals, followed by larger block sizes, are responsible for stable estimates of the CMV, as suggested by the higher autocorrelations. The choice of RGB channels has a limited effect on the quality of CMVs, and the red and the grayscale images are marginally more reliable than the other combinations during rapidly evolving low-level clouds. The stability of CMVs was tested at different image resolutions with an implementation of the optimized algorithm on the Sage cyberinfrastructure test bed. We find that doubling the frame rate outperforms quadrupling the image resolution in achieving CMV stability. The correlations of CMVs with the wind data are significant in the range of 0.38–0.59 with a 95 % confidence interval, despite the uncertainties and limitations of both datasets. A comparison of the PC method with constructed data and the optical flow method suggests that the post-processing of the vector field has a significant effect on the quality of the CMV. The raindrop-contaminated images can be identified by the rotation of the TSI mirror in the motion field. The results of this study are critical to optimizing algorithms for edge-computing sensor systems.
An experimental investigation was conducted at Mach 6 to determine the hypersonic aerodynamic characteristics of an all-body, delta-planform, hypersonic research aircraft (HYFAC configuration). The aerodynamic characteristics were obtained at Reynolds numbers based on model length of 2.84 million and 10.5 million and over an angle-of-attack range from minus 4 deg to 20 deg. The experimental results show that the HYFAC configuration is longitudinally stable and can be trimmed over the range of test conditions. The configuration had a small degree of directional stability over the angle-of-attack range and positive effective dihedral at angles of attack greater than 2 deg. Addition of canards caused a decrease in longitudinal stability and an increase in directional stability. Oil-flow studies revealed extensive areas of separated and vortex flow on the fuselage lee surface. A limited comparison of wind-tunnel data with several hypersonic approximations indicated that, except for the directional stability, the tangent-cone method gave adequate agreement at control settings between 5 deg and minus 5 deg and positive lift coefficient. A limited comparison indicated that the HYFAC configuration had greater longitudinal stability than an elliptical-cross-section configuration, but a lower maximum lift-drag ratio.
This paper discusses the application of robust control theory to evaluate robustness of the Ares-I control systems. Three techniques for estimating upper and lower bounds of uncertain parameters which yield stable closed-loop response are used here: (1) Monte Carlo analysis, (2) mu analysis, and (3) characteristic frequency response analysis. All three methods are used to evaluate stability envelopes of the Ares-I control systems with uncertain aerodynamics and flex dynamics. The results show that characteristic frequency response analysis is the most effective of these methods for assessing robustness.
The increasing penetration of converter-based renewable generation has resulted in faster frequency dynamics, and low and variable inertia. As a result, there is a need for frequency control methods that are able to stabilize a disturbance in the power system at timescales comparable to the fast converter dynamics. This paper proposes a combined linear and neural network controller for inverter-based primary frequency control that is stable at time-varying levels of inertia. We model the time-variance in inertia via a switched affine hybrid system model. We derive stability certificates for the proposed controller via a quadratic candidate Lyapunov function. We test the proposed control on a 12-bus 3-area test network, and compare its performance with a base case linear controller, optimized linear controller, and finite-horizon Linear Quadratic Regulator (LQR). Our proposed controller achieves faster mean settling time and over 50% reduction in average control cost across 100 inertia scenarios compared to the optimized linear controller. Unlike LQR which requires complete knowledge of the inertia trajectories and system dynamics over the entire control time horizon, our proposed controller is real-time tractable, and achieves comparable performance to LQR.
In this work, we propose a new method that extends conservative explicit multirate methods to implicit explicit-multirate methods. We develop extensions of order one and two with different stability properties on the implicit side. The method is suitable for time-stepping adaptive mesh refinement PDE discretizations with different degrees of stiffness. A numerical example with an advection-diffusion problem illustrates the new method's properties.
The properties of the flat spin mode of a general aviation configuration have been studied through analysis of rotary balance data, numerical simulation, and analytical study of the equilibrium state. The equilibrium state is predicted well from rotary balance data. The variations of yawing moment and pitching moment as functions of sideslip have been shown to be of great importance in obtaining accurate modeling. These dependencies are not presently available with sufficient accuracy from previous tests or theories. The stability of the flat spin mode has been examined extensively using numerical linearization, classical perturbation methods, and reduced order modeling. The stability exhibited by the time histories and the eigenvalue analyses is shown to be strongly dependent on certain static cross derivatives and more so on the dynamic derivatives. Explicit stability criteria are obtained from the reduced order models.
Doppler global velocimetry uses the absorption characteristics of iodine vapor to provide instantaneous three-component measurements of flow velocity within a plane defined by a laser light sheet. Although the technology is straightforward, its utilization as a flow diagnostics tool requires hardening of the optical system and careful attention to detail during data acquisition and processing if routine use in wind tunnel applications is to be achieved. A development program that reaches these goals is presented. Theoretical and experimental investigations were conducted on each technology element to determine methods that increase measurement accuracy and repeatability. Enhancements resulting from these investigations included methods to ensure iodine vapor calibration stability, single frequency operation of the laser and image alignment to sub-pixel accuracies. Methods were also developed to improve system calibration, and eliminate spatial variations of optical frequency in the laser output, spatial variations in optical transmissivity and perspective and optical distortions in the data images. Each of these enhancements is described and experimental examples given to illustrate the improved measurement performance obtained by the enhancement. The culmination of this investigation was the measured velocity profile of a rotating wheel resulting in a 1.75% error in the mean with a standard deviation of 0.5 m/s. Comparing measurements of a jet flow with corresponding Pitot measurements validated the use of these methods for flow field applications.
A recent NASA application of a remotely-piloted drop model to studies of the high angle-of-attack and spinning characteristics of a fighter configuration has provided an opportunity to evaluate and develop parameter estimation methods for the complex aerodynamic environment associated with high angles of attack. The paper discusses the overall drop model operation including descriptions of the model, instrumentation, launch and recovery operations, piloting concept, and parameter identification methods used. Static and dynamic stability derivatives were obtained for an angle-of-attack range from -20 deg to 53 deg. The results of the study indicated that the variations of the estimates with angle of attack were consistent for most of the static derivatives, and the effects of configuration modifications to the model (such as nose strakes) were apparent in the static derivative estimates. The dynamic derivatives exhibited greater uncertainty levels than the static derivatives, possibly due to nonlinear aerodynamics, model response characteristics, or additional derivatives.
Efficient and accurate poroelasticity models are critical in modeling geophysical problems such as oil exploration, gas-hydrate detection, and hydrogeology. We propose an efficient operator splitting method for Biot’s model of linear poroelasticity based on fixed-point iteration and constrained stress. In this method, we eliminate the constraint on time step via combining our fixed-point approach with a physics-based restraint between iterations. Three different cases are considered to demonstrate the stability and consistency of the method for constant and variable parameters. The results are validated against the results from the fully coupled approach. In case I, a single iteration is used for continuous coefficients. The relative error decreases with an increase in time. In case II, material coefficients are assumed to be linear. In the single iteration approach, the relative error grows significantly to 40% before rapidly decaying to zero. This is an artifact of the approximate solutions approaching the asymptotic solution. The error in the multiple iterations oscillates within 10 − 6 before decaying to the asymptotic solution. Nine iterations per time step are enough to achieve the relative error close to 10 − 7 . In the last case, the hybrid method with multiple iterations requires approximately 16 iterations to make the relative error 5 × 10 − 6 .
Compressible stability equations are solved using the spectral collocation method in an attempt to study the effects of temperature difference and compressibility on the stability of Taylor-Couette flow. It is found that the Chebyshev collocation spectral method yields highly accurate results using fewer grid points for solving stability problems. Comparisons are made between the result obtained by assuming small Mach number with a uniform temperature distribution and that based on fully incompressible analysis.
Stabilizer codes are a powerful method for implementing fault-tolerant quantum memory and in the case of topological codes, they form useful models for topological phases of matter. In this paper, we discuss the theory of stabilizer codes as a discrete version of a linear gauge structure, a concept we introduce here. A linear gauge structure captures all the familiar ingredients of U(1) gauge theory including a generalization of charge conservation, Maxwell equations and topological terms. Using this connection, we prove some important results for stabilizer codes which can be used to characterize their error-correction properties. However, this perspective does not depend on any particular Hamiltonian or Lagrangian, which as a consequence is agnostic to the source of the gauge redundancy. Based upon the connection to stabilizer codes, we view the source of the gauge redundancy as an ambiguity in the tensor product structure of the Hilbert space. That is, the gauge provides a set of equivalent but arbitrary ways to factorize the Hilbert space. We apply this formalism to the d=2 and 3 toric code as well as the paradigm fracton models, X-Cube and Haah’s cubic code. From this perspective, we are able to map all these models to some continuum version which captures all or some of the fractonic and topological features of their parent discrete models.
The Coarse-Mesh Finite Difference (CMFD) method has been used for over a decade to accelerate the convergence of the Method of Characteristics (MOC) solution to the two- dimensional particle transport equation in CASMO5. Numerical testing, along with widespread use in production-level calculations, have shown that the current CMFD implementation provides stability and robustness for a wide range of realistic reactor physics problems. However, the recent development of linear prolongation has attracted attention from the community as a way to further improve the performance and stability of CMFD. Two interpolation methods for linear prolongation are presented in this work and implemented into a test version of CASMO5. The performance of the proposed interpolations, referred to as the bilinear and linear directional schemes, is evaluated in terms of runtime relative to the default constant or uniform scaling update. Numerical results indicate that the use of linear prolongation can reduce the transport solver runtime on average by approximately 10% when tested with two hundred randomly selected cases. The new directional linear interpolation, combined with default constant boundary updates, is found to provide the highest reduction in runtime for the cases analyzed. (authors)
The principles of the cross method of moment distribution are used to check the stability of structural members under axial load. A brief theoretical treatment of the subject, together with an illustrative problem, is included as well as a discussion of the reduced modulus at high stresses and a set of tables to aid in the solution of practical problems.
A method of interpolating atmospheric soundings while reducing the errors associated with simple time interpolation was developed. The purpose of this was to provide a means to determine atmospheric stability at times between standard soundings and to relate changes in stability to intensity changes in an MCC. Four MCC cases were chosen for study with this method with four stability indices being included. The discussion centers on three aspects for each stability parameter examined: the stability field in the vicinity of the storm and its changes in structure and magnitude during the lifetime of the storm, the average stability within the storm boundary as a function of time and its relation to storm intensity, and the apparent flux of stability parameter into the storm as a consequence of low-level storm relative flow. It was found that the results differed among the four stability parameters, sometimes in a conflicting fashion. Thus, an interpolation of how the storm intensity is related to the changing environmental stability depends upon the particular index utilized. Some explanation for this problem is offered.