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 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Amplitude Analysis of ωπ0 Photoproduction at GlueX

spectrum of light mesons produced from a linearly polarized photon beam. The production and decays of a light meson resonance X such as γp → Xp′ →ωπ0p′ can be modeled with polarized vector-pseudoscalar ampli- tudes, which can describe the contribution of individual amplitudes to the total measured intensity. The status of mass-independent fits to the ωπ0 mass spectrum over a wide range of mass and momentum transfer−twill be presented, with an emphasis on interactions with the b1(1235) meson. We will also present in parallel an analysis of moments of the angular dis- tributions for the same process, as a means of verifying the stability of the amplitude-based results. These results will help to improve the broader knowledge of light meson states and how they are produced.

Scheuer, Kevin [College of William and Mary, Willi↗
NASA NTRS1997

SiC Recession Due to SiO2 Scale Volatility Under Combustor Conditions

One of today's most important and challenging technological problems is the development of advanced materials and processes required to design and build a fleet of supersonic High Speed Civil Transport (HSCT) airliners, a follow-up to the Concorde SST. The innovative combustor designs required for HSCT engines will need high-temperature materials with long-term environmental stability. Higher combustor liner temperatures than today's engines and the need for lightweight materials will require the use of advanced ceramic-matrix composites (CMC's) in hot-section components. The HSCT is just one example being used to demonstrate the need for such materials. This thesis evaluates silicon carbide (SiC) as a potential base material for HSCT and other similar applications. Key issues are the environmental durability for the materials of interest. One of the leading combustor design schemes leads to an environment which will contain both oxidizing and reducing gas mixtures. The concern is that these environments may affect the stability of the silica (SiO2) scale on which SiC depends for environmental protection. A unique High Pressure Burner Rig (HPBR) was developed to simulate the combustor conditions of future gas turbine engines, and a series of tests were conducted on commercially available SiC material. These tests are intended as a feasibility study for the use of these materials in applications such as the HSCT. Linear weight loss and surface recession of the SiC is observed as a result of SiO2 volatility for both fuel-lean and fuel-rich gas mixtures. These observations are compared and agree well with thermogravimetric analysis (TGA) experiments. A strong Arrhenius-type temperature dependence exists. In addition, the secondary dependencies of pressure and gas velocity are defined. As a result, a model is developed to enable extrapolation to points outside the experimental space of the burner rig, and in particular, to potential gas turbine engine conditions.

Robinson, Raymond Craig↗
DOE OSTI2019

Experimental and Theoretical Study of Cyclopentanone as a Catalyst for Low Temperature Alkene Oxidation: Preprint

Cyclopentanone was evaluated for spark-ignition fuel properties to determine potential as a bio-blendstock. Although key properties such as octane numbers and energy density are promising, when blended with commercial gasoline the oxidation stability was significantly reduced. Oxidation stability determination exposes fuel to mild oxidative conditions of 700 kPa oxygen followed by temperature increase to 100 degrees C, resulting in a final pressure of ~ 900 kPa. Oxygen consumption is determined by monitoring for decrease in pressure over time - the results evaluate potential for autoxidation during storage. An increase in oxidation rate from cyclopentanone was not anticipated given the relatively low reactivity of this compound. Detailed hydrocarbon analysis of blends before and after oxidation revealed cyclopentanone was not consumed, but alkenes were oxidized to a significantly higher degree with this compound present. Experiments with surrogate blends of isooctane, cyclopentanone, and linear isomers of hexene demonstrated that cyclopentanone catalyzes alkene oxidation to form epoxides under mild oxidative conditions. This unexpected behavior observed may have implications in low temperature combustion of alkenes when ketone fuels are present, as well as epoxidation chemistry. This study proposes a detailed mechanism for a catalytic cycle evaluated with quantum chemical calculations performed with density functional theory.

09 BIOMASS FUELS↗
DOE OSTI2021

Perspectives on the Calibration of CNN Energy Reconstruction in Highly Granular Calorimeters

We present a study which shows encouraging stability of the response linearity for a simulated high granularity calorimeter module reconstructed by a CNN model to miscalibration, bias, and noise effects. Our results also show an intuitive, quantifiable relationship between these factors and the calibration parameters. We trained a CNN model to reconstruct energy in the calorimeter module using simulated single-pion events; we then observed the response of the model under various miscalibration, bias, and noise conditions that affected the model input. From these data, we estimated linear response models to calibrate the CNN. We also quantified the relationship between these factors and the calibration parameters by regression analysis.

Akchurin, Nural↗
DOE OSTI2019

Quench-Back Management for Fast Decaying Currents in SHMS Superconducting Magnets at Jefferson Lab

The Super High Momentum Spectrometer (SHMS) of Hall C, part of the 12-GeV upgrade at Jefferson Lab, was successfully commissioned in 2017. Early operation demonstrated that fast current dumps decays of the SHMS Q2/Q3 and dipole superconducting magnets triggered quenches, causing some level of operational difficulty. Tests and detailed assessments suggest that nonquench induced fast current decay could result in substantial ac loss in the conductor and subsequently trigger a quench-back effect. The OPERA/ELECTRA software package was used to calculate the amount of heat deposited in the copper stabilizer from a fast current decay. The resistances of external energy dump resistors were lowered to slow the fast dumping of the magnet's current to reduce or eliminate the quench-back effect. A worst-case adiabatic quench scenario was analyzed, assuming no external dump resistor and no liquid helium surrounding the coil, to ensure the safety of the magnets. The stress levels in the coil imposed by winding, collaring preload, Lorentz force, and temperature gradient during a quench were also examined. The Tsai-Wu material failure criterion was used to determine the acceptable combined stress level. Linear orthotropic analysis of the coil indicates that the magnets can be operated safely with appropriately sized dump resistors. Fast dump tests with the modified dump resistors have been planned to verify the performance and suitability.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1992

Flutter analysis of supersonic axial flow cascades using a high resolution Euler solver. Part 1: Formulation and validation

This report presents, in two parts, a dynamic aeroelastic stability (flutter) analysis of a cascade of blades in supersonic axial flow. Each blade of the cascade is modeled as a typical section having pitching and plunging degrees of freedom. Aerodynamic forces are obtained from a time accurate, unsteady, two-dimensional cascade solver based on the Euler equations. The solver uses a time marching flux-difference splitting (FDS) scheme. Flutter stability is analyzed in the frequency domain. The unsteady force coefficients required in the analysis are obtained by harmonically oscillating (HO) the blades for a given flow condition, oscillation frequency, and interblade phase angle. The calculated time history of the forces is then Fourier decomposed to give the required unsteady force coefficients. An influence coefficient (IC) method and a pulse response (PR) method are also implemented to reduce the computational time for the calculation of the unsteady force coefficients for any phase angle and oscillation frequency. Part 1, this report, presents these analysis methods and their validation by comparison with results obtained from linear theory for a selected flat plate cascade geometry. A typical calculation for a rotor airfoil is also included to show the applicability of the present solver for airfoil configurations. The predicted unsteady aerodynamic forces for a selected flat plate cascade geometry and flow conditions correlated well with those obtained from linear theory for different interblade phase angles and oscillation frequencies. All the three methods of predicting unsteady force coefficients, namely, HO, IC, and PR, showed good correlations with each other. It was established that only a single calculation with four blade passages is required to calculate the aerodynamic forces for any phase angle for a cascade consisting of any number of blades, for any value of the oscillation frequency. Flutter results, including mistuning effects, for a cascade of stator airfoils are presented in Part 2 of the report.

Reddy, T. S. R.↗
NASA NTRSDate not supplied

On A Stabilization of the Ingard-Myers Impedance Boundary Condition

It has been well-known that the Ingard-Myers impedance condition, while simple to apply, is subject to the hydrodynamic Kelvin-Helmholtz instability due to its use of a vortex sheet in modeling the flow at the liner boundary. Recently, in the development of a time domain boundary element method for acoustic scattering by treated surfaces, it was found that by neglecting a certain second-order spatial derivative term in the Ingard-Myers formulation, the hydrodynamic instability can be avoided. The present paper aims to provide further analysis of this modified condition, hereby referred to as Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC). It will be shown, based on the dispersion relations of linear waves, that the instability intrinsic to the Ingard-Myers condition is eliminated in the proposed new formulation. Quantitative assessments on the accuracy of TIMIBC for scattering of acoustic waves by lined surfaces will be carried out, and its effectiveness will be demonstrated by numerical examples. Specifically, the accuracy is assessed by comparing solutions obtained by the Ingard-Myers condition with that by the proposed TIMIBC where theoretical reflection coefficients at a lined surface are derived for cases of plane and spherical incident waves. It is found that the TIMIBC provides a good approximation to the original Ingard-Myers condition for flows of low to mid subsonic Mach numbers. Time domain implementations of TIMIBC are also discussed and illustrated with a numerical example using a finite difference scheme. As many studies have shown the Ingard-Myers condition to be the correct limit of boundary layer thickness going to zero, the proposed TIMIBC can offer a practical solution for overcoming the intrinsic instability associated with the Ingard-Myers condition.

Time Domain BEM↗
NASA NTRSDate not supplied

On A Stabilization of the Ingard-Myers Impedance Boundary Condition

It has been well-known that the Ingard-Myers impedance condition, while simple to apply, is subject to the hydrodynamic Kelvin-Helmholtz instability due to its use of a vortex sheet in modeling the flow at the liner boundary. Recently, in the development of a time domain boundary element method for acoustic scattering by treated surfaces, it was found that by neglecting a certain second-order spatial derivative term in the Ingard-Myers formulation, the hydrodynamic instability can be avoided. The present paper aims to provide further analysis of this modified condition, hereby referred to as Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC). It will be shown, based on the dispersion relations of linear waves, that the instability intrinsic to the Ingard-Myers condition is eliminated in the proposed new formulation. Quantitative assessments on the accuracy of TIMIBC for scattering of acoustic waves by lined surfaces will be carried out, and its effectiveness will be demonstrated by numerical examples. Specifically, the accuracy is assessed by comparing solutions obtained by the Ingard-Myers condition with that by the proposed TIMIBC where theoretical reflection coefficients at a lined surface are derived for cases of plane and spherical incident waves. It is found that the TIMIBC provides a good approximation to the original Ingard-Myers condition for flows of low to mid subsonic Mach numbers. Time domain implementations of TIMIBC are also discussed and illustrated with a numerical example using a finite difference scheme. As many studies have shown the Ingard-Myers condition to be the correct limit of boundary layer thickness going to zero, the proposed TIMIBC can offer a practical solution for overcoming the intrinsic instability associated with the Ingard-Myers condition.

Time Domain BEM↗
DOE OSTI2014

Phosphine-stabilised Au{sub 9} clusters interacting with titania and silica surfaces: The first evidence for the density of states signature of the support-immobilised cluster

Chemically made, atomically precise phosphine-stabilized clusters Au{sub 9}(PPh{sub 3}){sub 8}(NO{sub 3}){sub 3} were deposited on titania and silica from solutions at various concentrations and the samples heated under vacuum to remove the ligands. Metastable induced electron spectroscopy was used to determine the density of states at the surface, and X-ray photoelectron spectroscopy for analysing the composition of the surface. It was found for the Au{sub 9} cluster deposited on titania that the ligands react with the titania substrate. Based on analysis using the singular value decomposition algorithm, the series of MIE spectra can be described as a linear combination of 3 base spectra that are assigned to the spectra of the substrate, the phosphine ligands on the substrate, and the Au clusters anchored to titania after removal of the ligands. On silica, the Au clusters show significant agglomeration after heat treatment and no interaction of the ligands with the substrate can be identified.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1999

Aeroelastic Analysis of a Trimmed Generic Hypersonic Vehicle

The aeroelastic equations of motion governing a hypersonic vehicle in free flight are derived. The equations of motion for a translating and rotating flexible body using Lagrange's equations in terms of quasi-coordinates are presented. These equations are simplified for the case of a vehicle with pitch and plunge rigid body degrees of freedom and small elastic displacements. The displacements are approximated by a truncated series of the unrestrained mode shapes, which are obtained using equivalent plate theory. Subsequently, the nonlinear equations of motion are linearized about the trim state, which is obtained using a rigid body trim model and steady hypersonic aerodynamics. The appropriate flutter derivatives are calculated from piston theory. Results describing mode shapes, trim behavior, and aeroelastic stability of a generic hypersonic vehicle are presented.

Nydick, I.↗
DOE OSTI2023

Evolution of coupled weakly driven waves in a dissipative plasma

The nonlinear collisional dynamics of coupled driven plasma waves in the presence of background dissipation is studied analytically within kinetic theory. Sufficiently near marginal stability, phase space correlations are poorly preserved and time delays become unimportant. The system is then shown to be governed by two first-order coupled autonomous differential equations of cubic order for the wave amplitudes and two complementary first-order equations for the evolution of their phases. That system of equations can be decoupled and further simplified to a single second-order differential equation of Liénard's type for each amplitude. Numerical solutions for this equation are obtained in the general case, while analytic solutions are obtained for special cases in terms of parameters related to the spacing of the resonances of the two waves in frequency space, e.g., wave lengths and oscillation frequencies. These parameters are further analyzed to find classes of quasi-steady saturation and pulsating scenarios. To classify equilibrium points, local stability analysis is applied, and bifurcation conditions are determined. When the two waves saturate at similar amplitude levels, their combined signal is shown to invariably exhibit amplitude beating and phase jumps of nearly π. In conclusion, the obtained analytical results can be used to benchmark simulations and to interpret eigenmode amplitude measurements in fusion experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A Lyapunov-Based Generalized Dc-Side Controller Design for PV-Connected Systems: Preprint

The objective of this paper is to realize a universal dc-side controller for photovoltaic (PV) systems where the control is agnostic to the downstream converter configuration. To achieve this, the downstream power converter and its controls are manipulated into an effective power control loop that is then cast into a generalized multi-loop design framework. On the dc side, a nonlinear small-signal model of the PV input is realized exclusively in terms of PV datasheet parameters (i.e., open-circuit voltage, short-circuit current, and maximum power point). Finally, a linear controller is used to modulate the dc-side PV system with the generic downstream power controller. A Lyapunov candidate is proposed to analyze the stability of the interconnected system and provide a streamlined approach for the controller design. The proposed design is validated on a 1-kVA experimental setup that interfaces a PV module to the grid.

dc-dc control↗
NASA NTRS1997

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized discretization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin discretization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobian linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. These variants have been implemented in the "ELF" library for which example calculations will be shown. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Some prevalent limiting strategies will be reviewed. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc. will be addressed as needed.

Barth, Timothy↗
NASA NTRS1998

Discretization and Preconditioning Algorithms for the Euler and Navier-Stokes Equations on Unstructured Meshes

Several stabilized demoralization procedures for conservation law equations on triangulated domains will be considered. Specifically, numerical schemes based on upwind finite volume, fluctuation splitting, Galerkin least-squares, and space discontinuous Galerkin demoralization will be considered in detail. A standard energy analysis for several of these methods will be given via entropy symmetrization. Next, we will present some relatively new theoretical results concerning congruence relationships for left or right symmetrized equations. These results suggest new variants of existing FV, DG, GLS, and FS methods which are computationally more efficient while retaining the pleasant theoretical properties achieved by entropy symmetrization. In addition, the task of Jacobean linearization of these schemes for use in Newton's method is greatly simplified owing to exploitation of exact symmetries which exist in the system. The FV, FS and DG schemes also permit discrete maximum principle analysis and enforcement which greatly adds to the robustness of the methods. Discrete maximum principle theory will be presented for general finite volume approximations on unstructured meshes. Next, we consider embedding these nonlinear space discretizations into exact and inexact Newton solvers which are preconditioned using a nonoverlapping (Schur complement) domain decomposition technique. Elements of nonoverlapping domain decomposition for elliptic problems will be reviewed followed by the present extension to hyperbolic and elliptic-hyperbolic problems. Other issues of practical relevance such the meshing of geometries, code implementation, turbulence modeling, global convergence, etc, will. be addressed as needed.

Barth, Timothy J.↗
NASA NTRS2002

Hubble Space Telescope Pointing Performance Due to Micro-Dynamic Disturbances from the NICMOS Cryogenic Cooler

The Vibration Admittance Test (VET) was performed to measure the emitted disturbances of the Near Infrared Camera and Multi-Object Spectrometer (NICMOS) Cryogenic Cooler (NCC) in preparation for NCC installation onboard the Hubble Space Telescope (HST) during Servicing Mission 3B (SM3B). Details of the VET ground-test are described, including facility characteristics, sensor complement and configuration, NCC suspension, and background noise measurements. Kinematic equations used to compute NCC mass center displacements and accelerations from raw measurements are presented, and dynamic equations of motion for the NCC VET system are developed and verified using modal test data. A MIMO linear frequency-domain analysis method is used to compute NCC-induced loads and HST boresight jitter from VET measurements. These results are verified by a nonlinear time-domain analysis approach using a high-fidelity structural dynamics and pointing control simulation for HST. NCC emitted acceleration levels not exceeding 35 micro-g rms were measured in the VET and analysis methods herein predict 3.1 milli-areseconds rms jitter for HST on-orbit. Because the NCC is predicted to become the predominant disturbance source for HST, VET results indicate that HST will continue to meet the 7 milli-arcsecond pointing stability mission requirement in the post-SM3B era.

Clapp, Brian R.↗
NASA NTRS1950

Theoretical Investigation of a Proportional-Plus-Flicker Automatic Pilot

The proportional-plus-flicker automatic pilot operates by a nonlinear principle whereby a fast-acting flicker servomotor response is combined with a low-speed proportional servomotor response for the purpose of obtaining supersonic stability and control. Essentially, the autopilot maintains a zero reference about which the output is proportional to the input. However, a flicker response overrides this proportional response at a fixed angle of gimbal displacement on either side of the zero gyroscope reference. Therefore, in contrast to other high speed control systems, the design requirements are simplified because the two components of the proportional-flicker control system are easy to build separately and they can be combined in a relatively simple manner. By application of the proportional-flicker principle, satisfactory stability can be obtained by the proper adjustment of the variable factors in the autopilot mechanism; namely, the proportional gain, the amplitude of flicker control deflection, the autopilot time-lag factor (the time-lag between flicker and proportional operation), and the point in the range that the autopilot switches from a flicker to a proportional system. There is a possibility that these factors can be adjusted so that a more rapid response time (the time to reach steady state) is obtained with the non-linear proportional-flicker autopilot than with a purely linear proportional autopilot. For the main part of this analysis, the proportional part of the system is approximated by a zero-phase-lag proportional autopilot with the assumption that the control surface moves instantaneously at the point where the system switches from flicker to proportional. Good correlation is shown between the results obtained by this method and results obtained by using a close approximation of an actual autopilot transfer function for proportional autopilot operation.

Seaberg, Ernest C.↗
DOE OSTI2016

Fourier Analysis of the Multilevel QD Method with {alpha}-Approximation For Time-Dependent Radiative Transfer Problems

Thermal radiation transfer (TRT) problems describe interaction of photon radiation with matter. They are defined by the time-dependent radiative transfer (RT) equation for the specific intensity I coupled with the energy balance (EB) equation. This class of problems is characterised by high dimensionality, multiple scales and strong nonlinearity. In this paper we analyze the multilevel quasi-diffusion (QD) method for solving TRT problems. This is a nonlinear method of moments with exact closures. The QD method enables one to reduce dimensionality of the original radiative transfer problem. It is formulated by a nonlinear system of equations consisting of (i) the high-order RT (HORT) equation for the intensity of photon radiation and (ii) low-order equations for the angular and energy moments of the intensity. The nonlinear closures are formulated by means of stable factors. The low-order problem brings the radiative transfer model to the level of multi-physics equations. The multilevel QD (MLQD) method possesses attractive features for efficient multi-physics coupling. To calculate the high-order RT solution we apply the {alpha}-approximation which assumes that the intensity changes exponentially in time over a time step. The {alpha}-approximation leads to the steady-state RT equation with a modified opacity and enables one to avoid storing the high-dimensional solution from the previous time step. The rate of change is evaluated by the solution of the time-dependent low-order QD (LOQD) equations for moments. We note that the HORT equation in the {alpha}-approximation contains an extra nonlinear term. It affects the iterative behaviour of the MLQD method. The obtained computational results indicated that the rate of convergence is slower compared to the rate when the time-dependent HORT is used to calculate the intensity. In this paper we perform an iterative stability analysis of the MLQD method in a discretized form and analyze convergence of outer (transport) iterations on a time step. We consider multigroup TRT problems in 1D slab geometry. The fully implicit scheme is used to discretize equations in time. The HORT equation is approximated with the method of step characteristics. A finite volume method is applied for spatial discretization of multigroup LOQD equations. The grey (one-group) LOQD equations are algebraically consistent with the multigroup LOQD equations. To perform Fourier analysis a model TRT problem in an infinite homogeneous medium is considered. The system of nonlinear equations of the MLQD method is linearized near an equilibrium solution at some given temperature. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Day-ahead photovoltaic power production forecasting methodology based on machine learning and statistical post-processing

A main challenge towards ensuring large-scale and seamless integration of photovoltaic systems is to improve the accuracy of energy yield forecasts, especially in grid areas of high photovoltaic shares. The scope of this paper is to address this issue by presenting a unified methodology for hourly-averaged day-ahead photovoltaic power forecasts with improved accuracy, based on data-driven machine learning techniques and statistical post-processing. More specifically, the proposed forecasting methodology framework comprised of a data quality stage, data-driven power output machine learning model development (artificial neural networks), weather clustering assessment (K-means clustering), post-processing output optimisation (linear regressive correction method) and the final performance accuracy evaluation. The results showed that the application of linear regression coefficients to the forecasted outputs of the developed day-ahead photovoltaic power production neural network improved the performance accuracy by further correcting solar irradiance forecasting biases. The resulting optimised model provided a mean absolute percentage error of 4.7% when applied to historical system datasets. Finally, the model was validated both, at a hot as well as a cold semi-arid climatic location, and the obtained results demonstrated close agreement by yielding forecasting accuracies of mean absolute percentage error of 4.7% and 6.3%, respectively. Finally, the validation analysis provides evidence that the proposed model exhibits high performance in both forecasting accuracy and stability.

42 ENGINEERING↗