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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Metrics of Balance Control for Use in Screening Tests of Vestibular Function

Decrements in balance control have been documented in astronauts after space flight. Reliable measures of balance control are needed for use in postflight field tests at remote landing sites. Diffusion analysis (DA) is a statistical mechanical tool that shows the average difference of the dependent variable on varying time scales. These techniques have been shown to measure differences in open-loop and closed-loop postural control in astronauts and elderly subjects. The goal of this study was to investigate the reliability of these measures of balance control. Eleven subjects were tested using the Clinical Test of Sensory Interaction on Balance: the subject stood with feet together and arms crossed on a stable or compliant surface, with eyes open or closed and with or without head movements in the pitch or yaw plane. Subjects were instrumented with inertial motion sensors attached to their trunk segment. The DA curves for linear acceleration measures were characterized by linear fits measuring open- (Ds) and closed-loop (Dl) control, and their intersection point (X-int, Y-int). Ds and Y-int showed significant differences between the test conditions. Additionally, Ds was correlated with the root mean square (RMS) of the signal, indicating that RMS was dominated by open-loop events (< 0.5 seconds). The Y-int was found to be correlated with the average linear velocity of trunk movements. Thus DA measures could be applied to derive reliable metrics of balance stability during field tests.

Fiedler, Matthew↗
DOE OSTI2019

Automated thermal extraction-desorption gas chromatography mass spectrometry: A multifunctional tool for comprehensive characterization of polymers and their degradation products

The TED-GC–MS analysis is a two-step method. A sample is first decomposed in a thermogravimetric analyzer (TGA) and the gaseous decomposition products are then trapped on a solid-phase adsorber. Subsequently, the solid-phase adsorber is analyzed with thermal desorption gas chromatography mass spectrometry (TDU-GC–MS). This method is ideally suited for the analysis of polymers and their degradation processes. A new entirely automated system is introduced which enables high sample throughput and reproducible automated fractioned collection of decomposition products. The fractionated collection together with low temperatures reduces the risk of contamination, improves instrumental stability and minimizes maintenance efforts. Through variation of the two main parameters (purge gas flow and heating rate) it is shown how the extraction process can be optimized. By measuring the decomposition products of polyethylene it is demonstrated that compounds with masses of up to 434 Da can be detected. This is achieved despite the low temperature (˜40 °C) of the solid-phase adsorber and the low thermal desorption temperature of 200 °C in the TDU unit. It is now shown that automated TED-GC–MS represents a new flexible multi-functional method for comprehensive polymer analyses. Comparable polymer characterization was previously only achievable through a combination of multiple independent analytical methods. This is demonstrated by three examples focused on practical challenges in materials analysis and identification: The first one is the analysis of wood plastic composites for which the decomposition processes of the polymer and the bio polymer (wood) could be clearly distinguished by fractionated collection using sequential adsorbers. Secondly, a fast quantitative application is shown by determining the weight concentrations of an unknown polyolefin blend through comparison with a reference material. Additionally, the determination of microplastic concentrations in environmental samples is becoming an increasingly important analytical necessity. It is demonstrated that with TED-GC–MS calibration curves showing good linearity for the most important precursors for microplastic, even complex matrix materials (suspended particulate matter) can be successfully analyzed.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Asymmetric charge separation and recombination in symmetrically functionalized σ–π hybrid oligosilanes

The flexibility of σ-conjugated silanes presents new opportunities for manipulating charge generation, transport, and non-linear optical properties of materials. Recently we synthesized a series of acceptor–donor–acceptor (ADA) compounds in which a methylated oligosilane core (D) is flanked by electron-deficient cyanovinyl-substituted arenes (A). Based on a detailed characterization of the photophysics of ADA and donor–acceptor (DA) architectures using both steady state and ultrafast spectroscopic measurements we illustrate that asymmetric charge separation occurs directly following light absorption. Lippert analysis of solvatochromic emission indicates large changes in dipole moments on excitation consistent with the formation of dipolar emissive states. Time resolved absorption measurements reveal common excited-state relaxation behavior across molecular structures: spectral dynamics associated with the relaxation of nascent excited states occur on a common timescale for all structures within the same solvent environment, whereas charge recombination via excited-state decay consistently follows a common energy gap law. Ultrafast time-resolved Raman measurements reveal that reduction of the cyanovinyl moieties is instantaneous with excitation, with only minor shifts in vibrational features over the course of excited-state relaxation. We conclude that excited-state symmetry breaking that gives rise to asymmetric intramolecular charge transfer (ICT) is associated with the conformation of the central Si chain. In contrast, ultrafast solvent reorganization or solvent-controlled intramolecular dynamics only serve to stabilize nascent dipolar excited states, rather than induce charge separation from an initial quadrupolar state.

Zhou, Jiawang↗
DOE OSTI2024

Exploration of the polarization angle variability of the Crab Nebula with POLARBEAR and its application to the search for axionlike particles

The Crab Nebula, also known as Tau A, is a polarized astronomical source at millimeter wavelengths. It has been used as a stable light source for polarization angle calibration in millimeter-wave astronomy. However, it is known that its intensity and polarization vary as a function of time at a variety of wavelengths. Thus, it is of interest to verify the stability of the millimeter-wave polarization. If detected, polarization variability may be used to better understand the dynamics of Tau A, and for understanding the validity of Tau A as a calibrator. One intriguing application of such observation is to use it for the search of axionlike particles (ALPs). Ultralight ALPs couple to photons through a Chern-Simons term, and induce a temporal oscillation in the polarization angle of linearly polarized sources. After assessing a number of systematic errors and testing for internal consistency, we evaluate the variability of the polarization angle of the Crab Nebula using 2015 and 2016 observations with the 150 GHz P instrument. We place a median 95% upper bound of polarization oscillation amplitude A < 0.06 5 ° over the oscillation frequencies from 0.75 year − 1 to 0.66 hour − 1 . Assuming that no sources other than ALP are causing Tau A’s polarization angle variation, that the ALP constitutes all the dark matter, and that the ALP field is a stochastic Gaussian field, this bound translates into a median 95% upper bound of ALP-photon coupling g a γ γ < 2.16 × 10 − 12 GeV − 1 × ( m a / 10 − 21 eV ) in the mass range from 9.9 × 10 − 23 eV to 7.7 × 10 − 19 eV . This demonstrates that this type of analysis using bright polarized sources is as competitive as those using the polarization of cosmic microwave background in constraining ALPs. Published by the American Physical Society 2024

Adachi, Shunsuke (ORCID:0000000204007555)↗
DOE OSTI2014

Influence of starting precursors and synthesis methods on the physiochemical properties of zirconia

Graphical abstract: Crystallite size of tetragonal phase of the zirconia samples prepared using different synthesis parameters and precursors as a function of calcination temperature. Surface area values of the zirconia samples calcined at 500 and 700 °C are in given brackets. - Highlights: • Zirconia prepared with modified sol–gel method is less stable compared with zirconia prepared by precipitation method. • Optimized synthesis conditions shifted the glow exotherm to higher temperature range indicating better thermal stability. • Tetragonal-zirconia could be synthesized in cost-effective manner using zirconium oxy-nitrate. • In our studies no co-relation between the surface area and crystallite size was observed. - Abstract: Under identical and judiciously pre-optimized synthesis conditions, the influence of different combinations of zirconium sources and/or post treatment conditions on structural properties, thermal stability, phase composition and morphology of zirconia has been investigated. High surface area tetragonal zirconia could be synthesized in a cost-effective manner from 1 M solution of zirconium oxy-nitrate at pH 11 using aqueous ammonia solution as a precipitant when calcined at 400 °C for 3 h. Irrespective of the preparation method, pH and starting precursor, zirconia samples prepared without digestion contained dominant monoclinic phase with some traces of tetragonal phase when calcined at 700 °C. Even though there is linear decrease in surface area with increase in the crystallite size for each sample as a function of calcination temperature, no co-relation between the surface area and crystallite size could be achieved. SEM images show agglomerated and irregular shape particles between 10 to 20 μm.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Comprehensive study of carbon dioxide adsorption in the metal–organic frameworks M 2 (dobdc) (M = Mg, Mn, Fe, Co, Ni, Cu, Zn)

Analysis of the CO 2 adsorption properties of a well-known series of metal–organic frameworks M 2 (dobdc) (dobdc 4– = 2,5-dioxido-1,4-benzenedicarboxylate; M = Mg, Mn, Fe, Co, Ni, Cu, and Zn) is carried out in tandem with in situ structural studies to identify the host–guest interactions that lead to significant differences in isosteric heats of CO 2 adsorption. Neutron and X-ray powder diffraction and single crystal X-ray diffraction experiments are used to unveil the site-specific binding properties of CO 2 within many of these materials while systematically varying both the amount of CO 2 and the temperature. Unlike previous studies, we show that CO 2 adsorbed at the metal cations exhibits intramolecular angles with minimal deviations from 180°, a finding that indicates a strongly electrostatic and physisorptive interaction with the framework surface and sheds more light on the ongoing discussion regarding whether CO 2 adsorbs in a linear or nonlinear geometry. This has important implications for proposals that have been made to utilize these materials for the activation and chemical conversion of CO 2 . For the weaker CO 2 adsorbents, significant elongation of the metal–O(CO 2 ) distances are observed and diffraction experiments additionally reveal that secondary CO 2 adsorption sites, while likely stabilized by the population of the primary adsorption sites, significantly contribute to adsorption behavior at ambient temperature. In conclusion, density functional theory calculations including van der Waals dispersion quantitatively corroborate and rationalize observations regarding intramolecular CO 2 angles and trends in relative geometric properties and heats of adsorption in the M 2 (dobdc)–CO 2 adducts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

A Compact, Single Stage, >1 kV Medium-Voltage Line Impedance Stabilization Network

This paper presents the design of a single-phase, single-stage line impedance stabilization network (LISN) for medium-voltage (MV) applications. More than 1 kV rated widebandgap (WBG) power semiconductor switches are increasingly utilized in new, emerging power electronics energy systems to improve power density and efficiency. However, due to inherently fast-switching speeds, WBG switch modules emit considerable electromagnetic interference (EMI) (e.g., common mode (CM) or differential mode (DM)). State-of-the-art offers standardized LISN solutions to validate and certify the new energy system for electromagnetic compatibility (EMC). However, most are for low voltage applications (i.e., < 1 kV). MV LISNs (i.e., > 1 kV) are rare until recently and literature does not provide sufficient guidelines for designing and characterizing such devices. It has been a critical challenge for many scientists and engineers to reliably certify the emerging MV energy systems (e.g., electric ships and aircraft). This paper addresses such a technology gap. Specifically, a CISPR 16-1-2 compliant 50Ω/50μH LISN with 1.5kV, 75A and 30MHz measurement capability has been proposed. Detailed performance study versus non-linear parasitic parameters variations in MV inductors and capacitors have been done. Based on new understandings, novel techniques to intuitively mitigate unwanted parasitic have been proposed to develop the proposed LISN successfully. Thorough characterization of important LISN parameters are presented and factors influencing them are analyzed. A rigorous analysis, experimental tests, and in-depth comparisons over state-of-the-art have been made to validate the effectiveness. This is done through the state-of-the-art 300kVA MV EMI testbed.

42 ENGINEERING↗
DOE OSTI2020

Rapid Open-Air Fabrication of Perovskite Solar Modules

We report on the open-air fabrication of perovskite solar modules with key advances including scalable large-area spray deposition, new monolithic integration scribing techniques, advanced photoluminescence characterization, and reproducible high-throughput manufacturability. Perovskite deposition with linear speeds of 12 m/min without a post-anneal are demonstrated with improved device performance, luminescent yield and >10X carrier lifetimes. Manufacturability using monolithic integration of series-connected modules is accomplished with a new indirect fiber laser ablation scribing method. A stable cell and module power output of 18.0% and 15.5%, respectively, was achieved with a subcell V oc > 1.06V. A comprehensive supporting technoeconomic analysis details the entire in-line manufacturing process from the glass substrate to the junction box of the encapsulated module. Here, the module manufacturing cost, balance of system costs and levelized cost of energy for a range of module efficiencies and lifetimes provides insights for the necessary tool speeds, efficiencies, and lifetimes for utility-scale energy generation.

14 SOLAR ENERGY↗
DOE OSTI2026

Spectrograph stabilization using a single-delay interferometer on the Hale Telescope

We describe a technique for spectrograph stabilization useful when conventional mitigation techniques of vacuum tanks, thermal insulation, and laser frequency comb may be impractical, expensive, heavy, or bulky. This includes spectrographs on airborne platforms or mounted on telescopes where they suffer a changing gravity vector or other drifts. Placing a fixed-delay interferometer in series with a spectrograph forms an externally dispersed interferometer (EDI). This produces a uniform sinusoidal comb multiplying input spectrum, creating (through heterodyning) beats (moiré patterns). In Fourier space for low frequencies up to the comb frequency, the moiré generated signal counter-rotates to ordinary spectra under an unknown disperser wavenumber drift Δx. This generates a large negative feedback signal useful in a conceptual control loop, to converge rapidly to a stable spectrum and yield Δx. A modified EDI data analysis algorithm (“crossfading”) combines frequency-weighted moiré with conventional spectrum to cancel net output spectrum reaction to Δx. Needing only a single-delay, this is a practical improvement over prior crossfading analyses requiring multiple delays. We test crossfading on ThAr data near 4850 cm−1 taken on Hale telescope in an earlier project. In a single pass, we reduce drift 20 times. Using seven iterations, we reduce 0.5 cm−1 (31 km/s Doppler equivalent) drift to 4×10−7 cm−1 (2.5 cm/s). The interferometer delay can wander, because linearity of phase versus wavenumber interpolates science features between bracketing calibrating spectral references. Second, mathematically reversing the heterodyning effect doubles effective spectral resolution without changing disperser slit.

Erskine, David J [Lawrence Livermore National Labo↗
NASA NTRS1994

A Feasibility Study for Perioperative Ventricular Tachycardia Prognosis and Detection and Noise Detection Using a Neural Network and Predictive Linear Operators

To locate the accessory pathway(s) in preexicitation syndromes, epicardial and endocardial ventricular mapping is performed during anterograde ventricular activation via accessory pathway(s) from data originally received in signal form. As the number of channels increases, it is pertinent that more automated detection of coherent/incoherent signals is achieved as well as the prediction and prognosis of ventricular tachywardia (VT). Today's computers and computer program algorithms are not good in simple perceptual tasks such as recognizing a pattern or identifying a sound. This discrepancy, among other things, has been a major motivating factor in developing brain-based, massively parallel computing architectures. Neural net paradigms have proven to be effective at pattern recognition tasks. In signal processing, the picking of coherent/incoherent signals represents a pattern recognition task for computer systems. The picking of signals representing the onset ot VT also represents such a computer task. We attacked this problem by defining four signal attributes for each potential first maximal arrival peak and one signal attribute over the entire signal as input to a back propagation neural network. One attribute was the predicted amplitude value after the maximum amplitude over a data window. Then, by using a set of known (user selected) coherent/incoherent signals, and signals representing the onset of VT, we trained the back propagation network to recognize coherent/incoherent signals, and signals indicating the onset of VT. Since our output scheme involves a true or false decision, and since the output unit computes values between 0 and 1, we used a Fuzzy Arithmetic approach to classify data as coherent/incoherent signals. Furthermore, a Mean-Square Error Analysis was used to determine system stability. The neural net based picking coherent/incoherent signal system achieved high accuracy on picking coherent/incoherent signals on different patients. The system also achieved a high accuracy of picking signals which represent the onset of VT, that is, VT immediately followed these signals. A special binary representation of the input and output data allowed the neural network to train very rapidly as compared to another standard decimal or normalized representations of the data.

Moebes, T. A.↗
NASA NTRS2023

JWST NIRCam Defocused Imaging Photometric Stability Performance and How It Can Sense Mirror Tilts

We use JWST NIRCam short-wavelength photometry to capture a transit lightcurve of the exoplanet HAT-P-14 b to assess performance as part of instrument commissioning. The short-wavelength precision is 152 ppm per 27 s integration as measured over the full time series compared to a theoretical limit of 107 ppm, after corrections to spatially correlated 1/f noise. Persistence effects from charge trapping are well fit by an exponential function with short characteristic timescales, settling on the order of 5–15 minutes. The short-wavelength defocused photometry is also uniquely well suited to measure the real-time wave-front error of JWST. Analysis of the images and reconstructed wave-front maps indicates that two different hexagonal primary mirror segments exhibited “tilt events,” where they changed orientation rapidly in less than ∼1.4 s. In some cases, the magnitude and timing of the flux jumps caused by tilt events can be accurately predicted with a telescope model. These tilt events can be sensed by simultaneous longer-wavelength NIRCam grism spectral images alone in the form of changes to the point-spread function, diagnosed from the full width at half maximum. They can also be sensed with the fine guidance sensor instrument from difference images. Tilt events possibly from sudden releases of stress in the backplane structure behind the mirrors were expected during the commissioning period because they were found in ground-based testing. Tilt events have shown signs of decreasing in frequency but have not disappeared completely. The detectors exhibit some minor (less than 1%) deviations from linear behavior in the first few groups of each integration, potentially impacting absolute fluxes and transit depths on bright targets, where only a handful of groups are possible. Overall, the noise is within 50% of the theoretical photon noise and read noise. This bodes well for high-precision measurements of transiting exoplanets and other time variable targets.

Everett Schlawin↗
NASA NTRS1993

Jupiter's winds and Arnol'd's second stability theorem: Slowly moving waves and neutral stability

Since the Voyager encounters in 1979, it has been known that Jupiter's cloud-top zonal winds violate the barotropic stability criterion. A vortex-tube stretching analysis of the Voyager wind data indicates that the more general Charney-Stern stability criterion is also violated. On the other hand, the zonal winds determined by tracking cloud features in Hubble Space Telescope images taken in 1991 precisely match the zonal winds determined by tracking cloud features in Voyager images, and it is hard to understand how a complicated zonal wind profile like Jupiter's could be unstable and yet not change at all in 12 years. In fact, there are at least two unknown ways to violate the Charney-Stern stability criterion and still have a stable flow. The better known of these is called Fjortoft's theorem, or Arnol'd's 1st theorem for the case of large-amplitude perturbations. Although the Fjortoft-Arnol'd theorem has been extended from the quasi-geostrophic equations to the primitive equations, the basic requirement that the potential vorticity be an increasing function of streamfunction is opposite to the case found in Jupiter, where the Voyager data indicate that the potential vorticity is a decreasing function of streamfunction. But this second case is precisely that which is covered by Arnol'd's 2nd stability theorem. In fact, the Voyager data suggest that Jupiter's zonal winds are neutrally stable with respect to Arnol'd's 2nd stability theorem. Here, we analyze the linear stability problem of a one-parameter family of sinusoidal zonal wind profiles that are close to neutral stability with respect to Arnol'd's 2nd stability theorem. We find numerically that the most unstable mode is always stationary, which may help to explain the slowly moving mode 10 waves observed on Jupiter. We find that violation of Arnol'd's 2nd stability theorem is both necessary and sufficient for instability of sinusoidal profiles. However, there appears to be no simple extension of Arnol'd's 2nd stability theorem to the primitive equations. Nevertheless, the primitive growth rates are small, and the primitive system is still governed by the quasi-geostrophic neutral-stability configuration.

Stamp, Andrew P.↗
NASA NTRS1991

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.↗
NASA NTRS2015

In-Flight Aeroelastic Stability of the Thermal Protection System on the NASA HIAD, Part II: Nonlinear Theory and Extended Aerodynamics

Conical shell theory and a supersonic potential flow aerodynamic theory are used to study the nonlinear pressure buckling and aeroelastic limit cycle behavior of the thermal protection system for NASA's Hypersonic Inflatable Aerodynamic Decelerator. The structural model of the thermal protection system consists of an orthotropic conical shell of the Donnell type, resting on several circumferential elastic supports. Classical Piston Theory is used initially for the aerodynamic pressure, but was found to be insufficient at low supersonic Mach numbers. Transform methods are applied to the convected wave equation for potential flow, and a time-dependent aerodynamic pressure correction factor is obtained. The Lagrangian of the shell system is formulated in terms of the generalized coordinates for all displacements and the Rayleigh-Ritz method is used to derive the governing differential-algebraic equations of motion. Aeroelastic limit cycle oscillations and buckling deformations are calculated in the time domain using a Runge-Kutta method in MATLAB. Three conical shell geometries were considered in the present analysis: a 3-meter diameter 70 deg. cone, a 3.7-meter 70 deg. cone, and a 6-meter diameter 70 deg. cone. The 6-meter configuration was loaded statically and the results were compared with an experimental load test of a 6-meter HIAD. Though agreement between theoretical and experimental strains was poor, the circumferential wrinkling phenomena observed during the experiments was captured by the theory and axial deformations were qualitatively similar in shape. With Piston Theory aerodynamics, the nonlinear flutter dynamic pressures of the 3-meter configuration were in agreement with the values calculated using linear theory, and the limit cycle amplitudes were generally on the order of the shell thickness. The effect of axial tension was studied for this configuration, and increasing tension was found to decrease the limit cycle amplitudes when the circumferential elastic supports were neglected, but resulted in more complex behavior when the supports were included. The nominal flutter dynamic pressure of the 3.7-meter configuration was significantly lower than that of the 3-meter, and it was found that two sets of natural modes coalesce to flutter modes near the same dynamic pressure. This resulted in a significant drop in the limit cycle frequencies at higher dynamic pressures, where the flutter mode with the lower frequency becomes more critical. Pre-buckling pressure loads and the aerodynamic pressure correction factor were studied for all geometries, and these effects resulted in significantly lower flutter boundaries compared with Piston Theory alone. The maximum dynamic pressure predicted by aerodynamic simulations of a proposed 3.7-meter HIAD vehicle was still lower than any of the calculated flutter dynamic pressures, suggesting that aeroelastic effects for this vehicle are of little concern.

Goldman, Benjamin D.↗
NASA NTRS2002

Stability and Heat Transfer Characteristics of Condensing Films

The overall objective of this research is to investigate the fundamental physics of film condensation in reduced gravity. The condensation of vapor on a cool surface is important in many engineering problems,including spacecraft thermal control and also the behavior of condensate films that may form on the interior surfaces of spacecraft. To examine the effects of body force on condensing films, two different geometries have been tested in the laboratory: (1) a stabilizing gravitational body force (+1g, or condensing surface facing 'upwards') and (2) de-stabilizing gravitational body force (-1g, or 'downwards'). For each geometry, different fluid configurations are employed to help isolate the fluid mechanical and thermal mechanisms operative in condensing films. The fluid configurations are (a) a condensing film, and (b) a non-condensing film with film growth by mass addition by through the plate surface. Condensation experiments are conducted in a test cell containing a cooled copper or brass plate with an exposed diameter of 12.7 cm. The metal surface is polished to allow for double-pass shadowgraph imaging, and the test surface is instrumented with imbedded heat transfer gauges and thermocouples. Representative shadowgraph images of a condensing, unstable (-1g) n-pentane film are shown. The interfacial disturbances associated with the de-stabilizing body force leading to droplet formation and break-off can be clearly seen. The heat transfer coefficient associated with the condensing film is shown. The heat transfer coefficient is seen to initially decrease, consistent with the increased thermal resistance due to layer growth. For sufficiently long time, a steady value of heat transfer is observed, accompanied by continuous droplet formation and break-off. The non-condensing cell consists of a stack of thin stainless steel disks 10 cm in diameter mounted in a brass enclosure. The disks are perforated with a regular pattern of 361 holes each 0.25 mm in diameter. Non-condensing experiments in -1g have employed 50 cSt and 125 cSt silicone oil pumped through the perforated disks at a specified rate by a syringe micropump. The time to droplet break-off and the disturbance wavelengths appear to decrease with increasing pumping rate. The ability to reliably perform multi-point, ultrasonic measurements of the film thickness has been demonstrated. A linear array of eight transducers of 6 mm diameter (with a beam footprint of comparable size) are pulsed with a square-wave signal at a frequency of 5 MHz and a pulse duration of approximately 0.3 s. For thin films (60 m to 2-3 mm in thickness) the layer thickness is determined by frequency analysis, where the received ultrasound pulse is Fourier transformed and the spacing between the peaks in the frequency spectrum is analyzed. For thicker layers (up to at least 1 cm in thickness), time-domain analysis is performed of the received ultrasound pulses to generate directly the layer thickness. A time-trace of the film thickness at a point using a single transducer in the linear array is shown for the case of an unstable (-1g) n-pentane film. The oscillations in film thickness are evidently due to the passage and/or shedding of droplets from the cooled plate surface. The entire transducer array was used to measure the changes in film thickness resulting from the passage of gravity waves generated either by an oscillating wall or the impact of a single droplet on the free surface of a film. The enclosure in both cases was 14 cm square and the transducer spacing was 12 mm. Best results were obtained using as test fluid a mixture of 50% glycerol and 50% water with a fluid layer thickness of 3-5 mm. In both cases the measured wavelengths and wave propagation speeds using the ultrasound technique compared reasonably well with those observed by optical imaging. Additional information can be found in the original extended abstract.

Hermanson, J. C.↗
DOE OSTI2017

Copper Sensing in Alkaline Electrolyte Using Anodic Stripping Voltammetry by Means of a Lead Mediator

Anodic stripping voltammetry (ASV) is an analysis technique that permits the selective and quantitative analysis of metal ion species in solution. It is most commonly applied in neutral to acidic electrolyte largely due to inherent metal ion solubility. Bismuth (Bi) is a common film used for ASV due to its good sensitivity, overall stability and insensitivity to O2. ASV, utilizing a Bi film, along with cadmium (Cd) and lead (Pb) as the plating mediators, has recently been adapted to determine zinc (Zn) concentrations in highly alkaline environments (30 % NaOH or 35 % M KOH). Successful analysis of Zn in alkaline relies on the ability of the hydroxide to form soluble metal anion species, such as Bi(OH) 4 – and Zn(OH) 4 2– . Here, we look to extend this technique to detect and quantify copper (Cu) ions in these highly basic electrolytes. However, in general, the use of ASV to detect and quantify Cu ion concentrations is notoriously difficult as the Cu stripping peak potential overlays with that of Bi from the common Bi film electrode. Here, an ASV method for determining Cu concentration in alkaline solutions is developed utilizing Pb as a deposition mediator. As such, it was found that when analyzing Cu solutions in the presence of Pb, the stripping voltammetry curves present separate and defined Cu stripping peaks. Different analyzes were made to find the best stripping voltammetry performance conditions. As such, an accumulation time of 5 minutes, an accumulation potential of ≤–1.45 V vs. Hg/HgO, and a concentration of 35 wt% KOH were determined to be the conditions that presented the best ASV results. Utilizing these conditions, calibration curves in the presence of 5.0 ppm Pb showed the best linear stripping signal correlation with an r-squared value of 0.991 and a limit of detection (LOD) of 0.67 ppm. Lastly, these results give way to evaluating Cu concentrations using ASV in aqueous alkaline solutions.

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

MASCOT - MATLAB Stability and Control Toolbox

MASCOT software was created to provide the conceptual aircraft designer accurate predictions of air vehicle stability and control characteristics. The code takes as input mass property data in the form of an inertia tensor, aerodynamic loading data, and propulsion (i.e. thrust) loading data. Using fundamental non-linear equations of motion, MASCOT then calculates vehicle trim and static stability data for any desired flight condition. Common predefined flight conditions are included. The predefined flight conditions include six horizontal and six landing rotation conditions with varying options for engine out, crosswind and sideslip, plus three takeoff rotation conditions. Results are displayed through a unique graphical interface developed to provide stability and control information to the conceptual design engineers using a qualitative scale indicating whether the vehicle has acceptable, marginal, or unacceptable static stability characteristics. This software allows the user to prescribe the vehicle s CG location, mass, and inertia tensor so that any loading configuration between empty weight and maximum take-off weight can be analyzed. The required geometric and aerodynamic data as well as mass and inertia properties may be entered directly, passed through data files, or come from external programs such as Vehicle Sketch Pad (VSP). The current version of MASCOT has been tested with VSP used to compute the required data, which is then passed directly into the program. In VSP, the vehicle geometry is created and manipulated. The aerodynamic coefficients, stability and control derivatives, are calculated using VorLax, which is now available directly within VSP. MASCOT has been written exclusively using the technical computing language MATLAB . This innovation is able to bridge the gap between low-fidelity conceptual design and higher-fidelity stability and control analysis. This new tool enables the conceptual design engineer to include detailed static stability and trim constraints in the conceptual design loop. The unique graphical interface developed for this tool presents the stability data in a format that is understandable by the conceptual designer, yet also provides the detailed quantitative results if desired.

Kenny, Sean↗
DOE OSTI2024

On the Stability of Power Transmission Systems Under Persistent Inverter Attacks: A Bi-Linear Matrix Approach

We investigate the stability and robustness properties of a power transmission system under persistent deceiving attacks on inverter-interfaced energy resources. The attacks can corrupt the damping coefficients in the inverters' controllers and measurements of the frequency at the points of coupling. Leveraging tools from hybrid dynamical systems theory, we characterize a broad family of persistent (and not necessarily periodic) attacks acting on the inverters, under which the stability properties of the transmission system can be shown to not be compromised. To address potentially conservative conditions identified through conventional bounding techniques, sufficient conditions on the average activation time of the attacks are identified via Lyapunov theory, as well as the formulation and solution of a class of bilinear matrix inequalities (BMI). The results are obtained for constant and slowly time-varying loads via input-to-state stability (ISS) tools. Numerical simulations on the IEEE 39-bus test system are also presented.

24 POWER TRANSMISSION AND DISTRIBUTION↗