Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “AMPLITUDES”

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,783 records · Page 99

High Resolution and Large Dynamic Range Resonant Pressure Sensor Based on Q-Factor Measurement

A pressure sensor has a high degree of accuracy over a wide range of pressures. Using a pressure sensor relying upon resonant oscillations to determine pressure, a driving circuit drives such a pressure sensor at resonance and tracks resonant frequency and amplitude shifts with changes in pressure. Pressure changes affect the Q-factor of the resonating portion of the pressure sensor. Such Q-factor changes are detected by the driving/sensing circuit which in turn tracks the changes in resonant frequency to maintain the pressure sensor at resonance. Changes in the Q-factor are reflected in changes of amplitude of the resonating pressure sensor. In response, upon sensing the changes in the amplitude, the driving circuit changes the force or strength of the electrostatic driving signal to maintain the resonator at constant amplitude. The amplitude of the driving signals become a direct measure of the changes in pressure as the operating characteristics of the resonator give rise to a linear response curve for the amplitude of the driving signal. Pressure change resolution is on the order of 10(exp -6) torr over a range spanning from 7,600 torr to 10(exp -6) torr. No temperature compensation for the pressure sensor of the present invention is foreseen. Power requirements for the pressure sensor are generally minimal due to the low-loss mechanical design of the resonating pressure sensor and the simple control electronics.

Gutierrez, Roman C.↗

An Experimental Study of Nonlinear Standing Waves in Resonators with Numerical Comparison

Lawrenson et. al. [Journal of the Acoustic Society of America, Nov. 1998] described the generation of shock-free high-amplitude pressure waves in closed cavities using large equipment and resonators to produce the reported effects. An attempt is made to generate shock-free high-amplitude pressure waves using relatively small resonators. Ambient air is used as the working fluid. A small cylindrical resonator is tested resulting in the lack of a shocked waveform while a larger model of the same shape produces shock waves. A small conical resonator produces shock-free pressure waves at resonance, but the amplitude of these waves is small. A larger cone resonator model produces shock-free pressure waves of higher amplitude. A large horn-cone resonator also produces shock-free high amplitude pressure waves, A numerical model is used to compare the experimental results to theoretical results. The effects of structural resonances on the production of shock-free high-amplitude pressure waves are discussed, especially concerning difficulties encountered when these resonances were in the frequency ranges of interest. Identifying features of a structural resonance are presented.

Finkbeiner, Joshua R.↗

Computation of Free-Free Transitions in Atomic Physics: Foundations

The amplitude T for "free-free" processes, such as bremsstrahlung or photo- absorption by an electron in the continuum in the presence of an external field, is usually written as the matrix element of the radiation operator taken between two continuum states. However, unlike the case when at least one of the states is bound, as in radiative transitions, electron capture, or the photo-effect, this expression contains an unphysical term, proportional to a delta-function and is not really the physical amplitude Tphys. This continues to be true for both the velocity and length form of the dipole approximation to the amplitude T. We first give an a priori definition of Tphys in terms of the scattering parts of the continuum functions, which has an obvious interpretation in terms of time-ordered diagrams. We then show that when the formal amplitude is modified by a long- distance cutoff, the modified form approaches Tphys as the cutoff is removed. The modified form then serves as a basis for the definition of a physical velocity dipole amplitude and this in turn leads to an equivalent length form of the dipole amplitude. This exercise provides a clear theoretical basis for many extant calculations in which cutoff factors are introduces somewhat ad hoc, as needed.

Bhatia, A. K.↗

Nonlinear Longitudinal Mode Instability in Liquid Propellant Rocket Engine Preburners

Nonlinear pressure oscillations have been observed in liquid propellant rocket instability preburner devices. Unlike the familiar transverse mode instabilities that characterize primary combustion chambers, these oscillations appear as longitudinal gas motions with frequencies that are typical of the chamber axial acoustic modes. In several respects, the phenomenon is similar to longitudinal mode combustion instability appearing in low-smoke solid propellant motors. An important feature is evidence of steep-fronted wave motions with very high amplitude. Clearly, gas motions of this type threaten the mechanical integrity of associated engine components and create unacceptably high vibration levels. This paper focuses on development of the analytical tools needed to predict, diagnose, and correct instabilities of this type. For this purpose, mechanisms that lead to steep-fronted, high-amplitude pressure waves are described in detail. It is shown that such gas motions are the outcome of the natural steepening process in which initially low amplitude standing acoustic waves grow into shock-like disturbances. The energy source that promotes this behavior is a combination of unsteady combustion energy release and interactions with the quasi-steady mean chamber flow. Since shock waves characterize the gas motions, detonation-like mechanisms may well control the unsteady combustion processes. When the energy gains exceed the losses (represented mainly by nozzle and viscous damping), the waves can rapidly grow to a finite amplitude limit cycle. Analytical tools are described that allow the prediction of the limit cycle amplitude and show the dependence of this wave amplitude on the system geometry and other design parameters. This information can be used to guide corrective procedures that mitigate or eliminate the oscillations.

Sims, J. D.↗

Sensorimotor integration in human postural control

It is generally accepted that human bipedal upright stance is achieved by feedback mechanisms that generate an appropriate corrective torque based on body-sway motion detected primarily by visual, vestibular, and proprioceptive sensory systems. Because orientation information from the various senses is not always available (eyes closed) or accurate (compliant support surface), the postural control system must somehow adjust to maintain stance in a wide variety of environmental conditions. This is the sensorimotor integration problem that we investigated by evoking anterior-posterior (AP) body sway using pseudorandom rotation of the visual surround and/or support surface (amplitudes 0.5-8 degrees ) in both normal subjects and subjects with severe bilateral vestibular loss (VL). AP rotation of body center-of-mass (COM) was measured in response to six conditions offering different combinations of available sensory information. Stimulus-response data were analyzed using spectral analysis to compute transfer functions and coherence functions over a frequency range from 0.017 to 2.23 Hz. Stimulus-response data were quite linear for any given condition and amplitude. However, overall behavior in normal subjects was nonlinear because gain decreased and phase functions sometimes changed with increasing stimulus amplitude. "Sensory channel reweighting" could account for this nonlinear behavior with subjects showing increasing reliance on vestibular cues as stimulus amplitudes increased. VL subjects could not perform this reweighting, and their stimulus-response behavior remained quite linear. Transfer function curve fits based on a simple feedback control model provided estimates of postural stiffness, damping, and feedback time delay. There were only small changes in these parameters with increasing visual stimulus amplitude. However, stiffness increased as much as 60% with increasing support surface amplitude. To maintain postural stability and avoid resonant behavior, an increase in stiffness should be accompanied by a corresponding increase in damping. Increased damping was achieved primarily by decreasing the apparent time delay of feedback control rather than by changing the damping coefficient (i.e., corrective torque related to body-sway velocity). In normal subjects, stiffness and damping were highly correlated with body mass and moment of inertia, with stiffness always about 1/3 larger than necessary to resist the destabilizing torque due to gravity. The stiffness parameter in some VL subjects was larger compared with normal subjects, suggesting that they may use increased stiffness to help compensate for their loss. Overall results show that the simple act of standing quietly depends on a remarkably complex sensorimotor control system.

NASA Discipline Neuroscience↗

Soleus H-reflex gain in humans walking and running under simulated reduced gravity

The Hoffmann (H-) reflex is an electrical analogue of the monosynaptic stretch reflex, elicited by bypassing the muscle spindle and directly stimulating the afferent nerve. Studying H-reflex modulation provides insight into how the nervous system centrally modulates stretch reflex responses.A common measure of H-reflex gain is the slope of the relationship between H-reflex amplitude and EMG amplitude. To examine soleus H-reflex gain across a range of EMG levels during human locomotion, we used simulated reduced gravity to reduce muscle activity. We hypothesised that H-reflex gain would be independent of gravity level.We recorded EMG from eight subjects walking (1.25 m s-1) and running (3.0 m s-1) at four gravity levels (1.0, 0.75, 0.5 and 0.25 G (Earth gravity)). We normalised the stimulus M-wave and resulting H-reflex to the maximal M-wave amplitude (Mmax) elicited throughout the stride to correct for movement of stimulus and recording electrodes relative to nerve and muscle fibres. Peak soleus EMG amplitude decreased by ~30% for walking and for running over the fourfold change in gravity. As hypothesised, slopes of linear regressions fitted to H-reflex versus EMG data were independent of gravity for walking and running (ANOVA, P > 0.8). The slopes were also independent of gait (P > 0.6), contrary to previous studies. Walking had a greater y-intercept (19.9% Mmax) than running (-2.5% Mmax; P < 0.001). At all levels of EMG, walking H-reflex amplitudes were higher than running H-reflex amplitudes by a constant amount. We conclude that the nervous system adjusts H-reflex threshold but not H-reflex gain between walking and running. These findings provide insight into potential neural mechanisms responsible for spinal modulation of the stretch reflex during human locomotion.

Non-NASA Center↗

What the Sunspot Record Tells Us About Space Climate

The records concerning the number, sizes, and positions of sunspots provide a direct means of characterizing solar activity over nearly 400 years. Sunspot numbers are strongly correlated with modem measures of solar activity including: 10.7-cm radio flux, total irradiance, x-ray flares, sunspot area, the baseline level of geomagnetic activity, and the flux of galactic cosmic rays. The Group Sunspot Number provides information on 27 sunspot cycles, far more than any of the modem measures of solar activity, and enough to provide important details about long-term variations in solar activity or Space Climate. The sunspot record shows: 1) sunspot cycles have periods of 131 plus or minus 14 months with a normal distribution; 2) sunspot cycles are asymmetric with a fast rise and slow decline; 3) the rise time from minimum to maximum decreases with cycle amplitude; 4) large amplitude cycles are preceded by short period cycles; 5 ) large amplitude cycles are preceded by high minima; 6) although the two hemispheres remain linked in phase, there are significant asymmetries in the activity in each hemisphere; 7) the rate at which the active latitudes drift toward the equator is anti-correlated with the cycle period, 8) the rate at which the active latitudes drift toward the equator is positively correlated with the amplitude of the cycle after the next; 9) there has been a significant secular increase in the amplitudes of the sunspot cycles since the end of the Maunder Minimum (1715); and 10) there is weak evidence for a quasi-periodic variation in the sunspot cycle amplitudes with a period of about 90 years. These characteristics indicate that the next solar cycle should have a maximum smoothed sunspot number of about 1.45 plus or minus 30 in 2010 while the following cycle should have a maximum of about 70 plus or minus 30 in 2023.

Hathaway, David H.↗

Effects of spaceflight on rhesus quadrupedal locomotion after return to 1G

Effects of spaceflight on Rhesus quadrupedal locomotion after return to 1G. Locomotor performance, activation patterns of the soleus (Sol), medial gastrocnemius (MG), vastus lateralis (VL), and tibialis anterior (TA) and MG tendon force during quadrupedal stepping were studied in adult Rhesus before and after 14 days of either spaceflight (n = 2) or flight simulation at 1G (n = 3). Flight simulation involved duplication of the spaceflight conditions and experimental protocol in a 1G environment. Postflight, but not postsimulation, electromyographic (EMG) recordings revealed clonus-like activity in all muscles. Compared with preflight, the cycle period and burst durations of the primary extensors (Sol, MG, and VL) tended to decrease postflight. These decreases were associated with shorter steps. The flexor (TA) EMG burst duration postflight was similar to preflight, whereas the burst amplitude was elevated. Consequently, the Sol:TA and MG:TA EMG amplitude ratios were lower following flight, reflecting a "flexor bias." Together, these alterations in mean EMG amplitudes reflect differential adaptations in motor-unit recruitment patterns of flexors and extensors as well as fast and slow motor pools. Shorter cycle period and burst durations persisted throughout the 20-day postflight testing period, whereas mean EMG returned to preflight levels by 17 days postflight. Compared with presimulation, the simulation group showed slight increases in the cycle period and burst durations of all muscles. Mean EMG amplitude decreased in the Sol, increased in the MG and VL, and was unchanged in the TA. Thus adaptations observed postsimulation were different from those observed postflight, indicating that there was a response unique to the microgravity environment, i.e., the modulations in the nervous system controlling locomotion cannot merely be attributed to restriction of movement but appear to be the result of changes in the interpretation of load-related proprioceptive feedback to the nervous system. Peak MG tendon force amplitudes were approximately two times greater post- compared with preflight or presimulation. Adaptations in tendon force and EMG amplitude ratios indicate that the nervous system undergoes a reorganization of the recruitment patterns biased toward an increased recruitment of fast versus slow motor units and flexor versus extensor muscles. Combined, these data indicate that some details of the control of motor pools during locomotion are dependent on the persistence of Earth's gravitational environment.

Bion 11 Project↗

Enhanced vasomotion of cerebral arterioles in spontaneously hypertensive rats

Intrinsic rhythmic changes in the diameter of pial cerebral arterioles (30-70 microns) in anesthetized normotensive and hypertensive rats were assessed in vivo to determine if any significant differences exist between the two strains. All diameter measurements were analyzed using a traditional graphic analysis technique and a new frequency spectrum analysis technique known as the Prony Spectral Line Estimator. Graphic analysis of the data revealed that spontaneously hypertensive rats (SHR) possess a significantly greater fundamental frequency (5.57 +/- 0.28 cycles/min) of vasomotion compared to the control Wistar-Kyoto normotensive rats (WKY) (1.95 +/- 0.37 cycles/min). Furthermore, the SHR cerebral arterioles exhibited a significantly greater amplitude of vasomotion (10.07 +/- 0.70 microns) when compared to the WKY cerebral arterioles of the same diameter (8.10 +/- 0.70 microns). Diameter measurements processed with the Prony technique revealed that the fundamental frequency of vasomotion in SHR cerebral arterioles (6.14 +/- 0.39 cycles/min) was also significantly greater than that of the WKY cerebral arterioles (2.99 +/- 0.42 cycles/min). The mean amplitudes of vasomotion in the SHR and WKY strains obtained by the Prony analysis were found not to be statistically significant in contrast to the graphic analysis of the vasomotion amplitude of the arterioles. In addition, the Prony system was able to consistently uncover a very low frequency of vasomotion in both strains of rats that was typically less than 1 cycle/min and was not significantly different between the two strains. The amplitude of this slow frequency was also not significantly different between the two strains. The amplitude of the slow frequency of vasomotion (less than 1 cycle/min) was not different from the amplitude of the higher frequency (2-6 cycles/min) vasomotion by Prony or graphic analysis. These data suggest that a fundamental intrinsic defect exists in the spontaneously hypertensive rat that may contribute to the pathogenesis of hypertension in these animals.

NASA Discipline Cardiopulmonary↗

Analysis of variability in the burst oscillations of the accreting millisecond pulsar XTE J1814-338

The accreting millisecond pulsar XTE J1814-338 exhibits oscillations at the known spin frequency during Type I X-ray bursts. The properties of the burst oscillations reflect the nature of the thermal asymmetry on the stellar surface. We present an analysis of the variability of the burst oscillations of this source, focusing on three characteristics: fractional amplitude, harmonic content and frequency. Fractional amplitude and harmonic content constrain the size, shape and position of the emitting region, whilst variations in frequency indicate motion of the emitting region on the neutron star surface. We examine both long-term variability over the course of the outburst, and short-term variability during the bursts. For most of the bursts, fractional amplitude is consistent with that of the accretion pulsations, implying a low degree of fuel spread. There is however a population of bursts whose fractional amplitudes are substantially lower, implying a higher degree of fuel spread, possibly forced by the explosive burning front of a precursor burst. For the first harmonic, substantial differences between the burst and accretion pulsations suggest that hotspot geometry is not the only mechanism giving rise to harmonic content in the latter. Fractional amplitude variability during the bursts is low; we can only rule out the hypothesis that the fractional amplitude remains constant at the l(sigma) level for bursts that do not exhibit photospheric radius expansion (PRE). There are no significant variations in frequency in any of the bursts except for the one burst that exhibits PRE. This burst exhibits a highly significant but small (= 0.1Hz) drop in frequency in the burst rise. The timescale of the frequency shift is slower than simple burning layer expansion models predict, suggesting that other mechanisms may be at work.

Watts, Anna L.↗

Noninvasive Diagnosis of Coronary Artery Disease Using 12-Lead High-Frequency Electrocardiograms

A noninvasive, sensitive method of diagnosing certain pathological conditions of the human heart involves computational processing of digitized electrocardiographic (ECG) signals acquired from a patient at all 12 conventional ECG electrode positions. In the processing, attention is focused on low-amplitude, high-frequency components of those portions of the ECG signals known in the art as QRS complexes. The unique contribution of this method lies in the utilization of signal features and combinations of signal features from various combinations of electrode positions, not reported previously, that have been found to be helpful in diagnosing coronary artery disease and such related pathological conditions as myocardial ischemia, myocardial infarction, and congestive heart failure. The electronic hardware and software used to acquire the QRS complexes and perform some preliminary analyses of their high-frequency components were summarized in Real-Time, High-Frequency QRS Electrocardiograph (MSC- 23154), NASA Tech Briefs, Vol. 27, No. 7 (July 2003), pp. 26-28. To recapitulate, signals from standard electrocardiograph electrodes are preamplified, then digitized at a sampling rate of 1,000 Hz, then analyzed by the software that detects R waves and QRS complexes and analyzes them from several perspectives. The software includes provisions for averaging signals over multiple beats and for special-purpose nonrecursive digital filters with specific low- and high-frequency cutoffs. These filters, applied to the averaged signal, effect a band-pass operation in the frequency range from 150 to 250 Hz. The output of the bandpass filter is the desired high-frequency QRS signal. Further processing is then performed in real time to obtain the beat-to-beat root mean square (RMS) voltage amplitude of the filtered signal, certain variations of the RMS voltage, and such standard measures as the heart rate and R-R interval at any given time. A key signal feature analyzed in the present method is the presence versus the absence of reduced-amplitude zones (RAZs). In terms that must be simplified for the sake of brevity, an RAZ comprises several cycles of a high-frequency QRS signal during which the amplitude of the high-frequency oscillation in a portion of the signal is abnormally low (see figure). A given signal sample exhibiting an interval of reduced amplitude may or may not be classified as an RAZ, depending on quantitative criteria regarding peaks and troughs within the reduced-amplitude portion of the high-frequency QRS signal. This analysis is performed in all 12 leads in real time.

Schlegel, Todd T.↗

An Assessment of Cumulative Axial and Torsional Fatigue in a Cobalt-Base Superalloy

Cumulative fatigue under axial and torsional loading conditions can include both load-order (higMow and low/high) as well as load-type sequence (axial/torsional and torsional/axial) effects. Previously reported experimental studies on a cobalt-base superalloy, Haynes 188 at 538 C, addressed these effects. These studies characterized the cumulative axial and torsional fatigue behavior under high amplitude followed by low amplitude (Kalluri, S. and Bonacuse, P. J., "Cumulative Axial and Torsional Fatigue: An Investigation of Load-Type Sequance Effects," in Multiaxial Fatigue and Deformation: Testing and Prediction, ASTM STP 1387, S. Kalluri, and P. J. Bonacuse, Eds., American Society for Testing and Materials, West Conshohocken, PA, 2000, pp. 281-301) and low amplitude followed by high amplitude (Bonacuse, P. and Kalluri, S. "Sequenced Axial and Torsional Cumulative Fatigue: Low Amplitude Followed by High Amplitude Loading," Biaxial/Multiaxial Fatigue and Fracture, ESIS Publication 31, A. Carpinteri, M. De Freitas, and A. Spagnoli, Eds., Elsevier, New York, 2003, pp. 165-182) conditions. In both studies, experiments with the following four load-type sequences were performed: (a) axial/axial, (b) torsional/torsional, (c) axial/torsional, and (d) torsional/axial. In this paper, the cumulative axial and torsional fatigue data generated in the two previous studies are combined to generate a comprehensive cumulative fatigue database on both the load-order and load-type sequence effects. This comprehensive database is used to examine applicability of the Palmgren-langer-Miner linear damage rule and a nonlinear damage curve approach for Haynes 188 subjected to the load-order and load-type sequencing described above. Summations of life fractions from the experiments are compared to the predictions from both the linear and nonlinear cumulative fatigue damage approaches. The significance of load-order versus load-type sequence effects for axial and torsional loading conditions is discussed. Possible reasons for the observed differences between the computed and observed summations of cycle fractions are rationalized in terms of the observed ever lutions of cyclic axial and shear stress ranges in the cumulative fatigue tests.

Kalluri, Sreeramesh↗

Systems, Methods and Apparatus for Position Sensor Digital Conditioning Electronics

Systems, methods and apparatus are provided through which in some implementations determine the amplitude of an amplitude modulated signal, modulated by the position of an object being sensed. In some aspects, the apparatus accepts an excitation signal and the amplitude modulated signal and divides the amplitude modulated by the excitation signal to produce an output signal that is proportional to the position of the object being sensed. In other aspects, the division is performed only when the excitation signal is non-zero, such as close to the peaks in the excitation signal. In other aspects, the excitation signal and amplitude modulated signal are degraded due to an air gap and the degraded signals are used to correct for amplitude fluctuations due to the air gap, and produce an output signal, tolerant of the air gaps, that is proportional to the position of the object being sensed.

Howard, David E.↗

Geoelectric Hazard Maps for the Continental United States

In support of a multiagency project for assessing induction hazards, we present maps of extreme-value geoelectric amplitudes over about half of the continental United States. These maps are constructed using a parameterization of induction: estimates of Earth surface impedance, obtained at discrete geographic sites from magnetotelluric survey data, are convolved with latitude-dependent statistical maps of extreme-value geomagnetic activity, obtained from decades of magnetic observatory data. Geoelectric amplitudes are estimated for geomagnetic waveforms having 240 s sinusoidal period and amplitudes over 10 min that exceed a once-per-century threshold. As a result of the combination of geographic differences in geomagnetic activity and Earth surface impedance, once-per-century geoelectric amplitudes span more than 2 orders of magnitude and are an intricate function of location. For north-south induction, once-per-century geoelectric amplitudes across large parts of the United States have a median value of 0.26 Vkm; for east-west geomagnetic variation the median value is 0.23 Vkm. At some locations,once-per-century geoelectric amplitudes exceed 3 Vkm.

Love, Jeffrey J.↗

A Segmentation Algorithm for Characterizing Rise and Fall Segments in Seasonal Cycles: An Application to XCO2 to Estimate Benchmarks and Assess Model Bias

There is more useful information in the time series of satellite-derived column-averaged carbon dioxide (XCO2) than is typically characterized. Often, the entire time series is treated at once without considering detailed features at shorter timescales, such as nonstationary changes in signal characteristics – amplitude, period and phase. In many instances, signals are visually and analytically differentiable from other portions in a time series. Each rise (increasing) and fall (decreasing) segment in the seasonal cycle is visually discernable in a graph of the time series. The rise and fall segments largely result from seasonal differences in terrestrial ecosystem production, which means that the segment's signal characteristics can be used to establish observational benchmarks because the signal characteristics are driven by similar underlying processes. We developed an analytical segmentation algorithm to characterize the rise and fall segments in XCO2 seasonal cycles. We present the algorithm for general application of the segmentation analysis and emphasize here that the segmentation analysis is more generally applicable to cyclic time series. We demonstrate the utility of the algorithm with specific results related to the comparison between satellite- and model-derived XCO2 seasonal cycles (2009–2012) for large bioregions across the globe. We found a seasonal amplitude gradient of 0.74–0.77 ppm for every 10∘ of latitude in the satellite data, with similar gradients for rise and fall segments. This translates to a south–north seasonal amplitude gradient of 8 ppm for XCO2, about half the gradient in seasonal amplitude based on surface site in situ CO2 data (∼19 ppm). The latitudinal gradients in the period of the satellite-derived seasonal cycles were of opposing sign and magnitude (−9 d per 10∘ latitude for fall segments and 10 d per 10∘ latitude for rise segments) and suggest that a specific latitude (∼2∘ N) exists that defines an inversion point for the period asymmetry. Before (after) the point of asymmetry inversion, the periods of rise segments are lesser (greater) than the periods of fall segments; only a single model could reproduce this emergent pattern. The asymmetry in amplitude and the period between rise and fall segments introduces a novel pattern in seasonal cycle analyses, but, while we show these emergent patterns exist in the data, we are still breaking ground in applying the information for science applications. Maybe the most useful application is that the segmentation analysis allowed us to decompose the model biases into their correlated parts of biases in amplitude, period and phase independently for rise and fall segments. We offer an extended discussion on how such information about model biases and the emergent patterns in satellite-derived seasonal cycles can be used to guide future inquiry and model development.

Calle, Leonardo↗

A Segmentation Algorithm for Characterizing Rise and Fall Segments in Seasonal Cycles: an Application to Xco2 to Estimate Benchmarks and Assess Model Bias

There is more useful information in the time series of satellite-derived column-averaged carbon dioxide (XCO2) than is typically characterized. Often, the entire time series is treated at once without considering detailed features at shorter timescales, such as nonstationary changes in signal characteristics – amplitude, period and phase. In many instances, signals are visually and analytically differentiable from other portions in a time series. Each rise (increasing) and fall (decreasing) segment in the seasonal cycle is visually discernable in a graph of the time series. The rise and fall segments largely result from seasonal differences in terrestrial ecosystem production, which means that the segment's signal characteristics can be used to establish observational benchmarks because the signal characteristics are driven by similar underlying processes. We developed an analytical segmentation algorithm to characterize the rise and fall segments in XCO2 seasonal cycles. We present the algorithm for general application of the segmentation analysis and emphasize here that the segmentation analysis is more generally applicable to cyclic time series. We demonstrate the utility of the algorithm with specific results related to the comparison between satellite- and model-derived XCO2 seasonal cycles (2009–2012) for large bioregions across the globe. We found a seasonal amplitude gradient of 0.74–0.77 ppm for every 10∘ of latitude in the satellite data, with similar gradients for rise and fall segments. This translates to a south–north seasonal amplitude gradient of 8 ppm for XCO2, about half the gradient in seasonal amplitude based on surface site in situ CO2 data (∼19 ppm). The latitudinal gradients in the period of the satellite-derived seasonal cycles were of opposing sign and magnitude (−9 d per 10∘ latitude for fall segments and 10 d per 10∘ latitude for rise segments) and suggest that a specific latitude (∼2∘ N) exists that defines an inversion point for the period asymmetry. Before (after) the point of asymmetry inversion, the periods of rise segments are lesser (greater) than the periods of fall segments; only a single model could reproduce this emergent pattern. The asymmetry in amplitude and the period between rise and fall segments introduces a novel pattern in seasonal cycle analyses, but, while we show these emergent patterns exist in the data, we are still breaking ground in applying the information for science applications. Maybe the most useful application is that the segmentation analysis allowed us to decompose the model biases into their correlated parts of biases in amplitude, period and phase independently for rise and fall segments. We offer an extended discussion on how such information about model biases and the emergent patterns in satellite-derived seasonal cycles can be used to guide future inquiry and model development.

segmentation algorithm↗

Australian tidal currents – assessment of a barotropic model (COMPAS v1.3.0 rev6631) with an unstructured grid

While the variations of tidal range are large and fairly well known across Australia (less than 1 m near Perth but more than 14 m in King Sound), the properties of the tidal currents are not. We describe a new regional model of Australian tides and assess it against a validation dataset comprising tidal height and velocity constituents at 615 tide gauge sites and 95 current meter sites. The model is a barotropic implementation of COMPAS, an unstructured-grid primitive-equation model that is forced at the open boundaries by TPXO9v1. The mean absolute error (MAE) of the modelled M2 height amplitude is 8.8 cm, or 12 % of the 73 cm mean observed amplitude. The MAE of phase (10°), however, is significant, so the M2 mean magnitude of vector error (MMVE, 18.2 cm) is significantly greater. The root sum square over the eight major constituents is 26 % of the observed amplitude. We conclude that while the model has skill at height in all regions, there is definitely room for improvement (especially at some specific locations). For the M2 major axis velocity amplitude, the MAE across the 95 current meter sites, where the observed amplitude ranges from 0.1 to 156 cm s −1 , is 6.9 cm s −1 , or 22 % of the 31.7 cm s −1 observed mean. This nationwide average result is encouraging, but it conceals a very large regional variation. Relative errors of the tidal current amplitudes on the narrow shelves of New South Wales (NSW) and Western Australia exceed 100 %, but tidal currents are weak and negligible there compared to non-tidal currents, so the tidal errors are of little practical significance. Looking nationwide, we show that the model has predictive value for much of the 79 % of Australia's shelf seas where tides are a major component of the total velocity variability. In descending order this includes the Bass Strait, the Kimberley to Arnhem Land, and southern Great Barrier Reef regions. There is limited observational evidence to confirm that the model is also valuable for currents in other regions across northern Australia. We plan to commence publishing “unofficial” tidal current predictions for chosen regions in the near future based on both our COMPAS model and the validation dataset we have assembled.

Tidal currents↗

An Experimental Investigation of Fuel Staging In A Model Gas Turbine Combustor

Modern gas turbine engines are typically operated in a lean-premixed configuration in order to meet emissions requirements. Lean operation makes gas turbines more susceptible to combustion instability, which can cause catastrophic hardware failure if left unchecked. The spatially uneven distribution of fuel in a gas turbine combustor, known as fuel staging, is an effective way of suppressing instabilities. This dissertation examines the effect of both steady-state and transient fuel-staging in a model gas turbine combustor. First, the steady-state effect of fuel staging on nominally unstable operating conditions is examined. At 200 degrees Celsius and a bulk flow velocity of 26 m/s, the combustor undergoes self-excited instabilities when the equivalence ratio is between φ = 0.70 and φ = 0.74 when all nozzles are fueled equally. These instabilities are suppressed when the center nozzle equivalence ratio is increased to φ Staging = 0.85 and the outer nozzle equivalence ratio maintained at φ = 0.70 or when the global equivalence ratio is maintained at φ = 0.70 by increasing the center nozzle equivalence ratio to φ Staging = 0.74 and decreasing the outer nozzle equivalence ratio to φ Staging = 0.69. The staged p' RMS amplitudes are found to decrease with increasing equivalence ratio and the damping rates found to increase with increasing equivalence ratio. Next, axisymmetric fuel staging transients are conducted valve opening/closing timescales of 1, 16, 58, 4,000, and 10,000 ms, and equivalence ratios of φ staging = 0.75, 0.80, 0.85. It should be noted that these valve opening/closing timescales, referred to as τ Transient , differ from the actual fuel convection timescale that are summarized in Table 5.2. It is found that the combustor pressure fluctuation at the end of the transient strongly depends on staging equivalence ratio but not on τ Transient . Increasing the staging equivalence ratio decreased the staged- stable combustor pressure amplitude, similarly to the steady-state results. The staged-stable damping rate trends are insensitive to to τ Transient and increased with increasing staging equivalence ratio. The timescale of the transition, τ, is found to be sensitive to both τ Transient and the staging equivalence ratio. For a Modern gas turbine engines are typically operated in a lean-premixed configuration in order to meet emissions requirements. Lean operation makes gas turbines more susceptible to combustion instability, which can cause catastrophic hardware failure if left unchecked. The spatially uneven distribution of fuel in a gas turbine combustor, known as fuel staging, is an effective way of suppressing instabilities. This dissertation examines the effect of both steady-state and transient fuel-staging in a model gas turbine combustor. First, the steady-state effect of fuel staging on nominally unstable operating conditions is examined. At 200 degrees Celsius and a bulk flow velocity of 26 m/s, the combustor undergoes self-excited instabilities when the equivalence ratio is between φ = 0.70 and φ = 0.74 when all nozzles are fueled equally. These instabilities are suppressed when the center nozzle equivalence ratio is increased to φ Staging = 0.85 and the outer nozzle equivalence ratio maintained at φ = 0.70 or when the global equivalence ratio is maintained at φ = 0.70 by increasing the center nozzle equivalence ratio to φ Staging = 0.74 and decreasing the outer nozzle equivalence ratio to φ Staging = 0.69. The staged p' RMS amplitudes are found to decrease with increasing equivalence ratio and the damping rates found to increase with increasing equivalence ratio. Next, axisymmetric fuel staging transients are conducted valve opening/closing timescales of 1, 16, 58, 4,000, and 10,000 ms, and equivalence ratios of φ staging = 0.75, 0.80, 0.85. It should be noted that these valve opening/closing timescales, referred to as τ Transient , differ from the actual fuel convection timescale that are summarized in Table 5.2. It is found that the combustor pressure fluctuation at the end of the transient strongly depends on staging equivalence ratio but not on τ Transient . Increasing the staging equivalence ratio decreased the staged-stable combustor pressure amplitude, similarly to the steady-state results. The staged-stable damping rate trends are insensitive to to τ Transient and increased with increasing staging equivalence ratio. The timescale of the transition, τ, is found to be sensitive to both τ Transient and the staging equivalence ratio. For a Modern gas turbine engines are typically operated in a lean-premixed configuration in order to meet emissions requirements. Lean operation makes gas turbines more susceptible to combustion instability, which can cause catastrophic hardware failure if left unchecked. The spatially uneven distribution of fuel in a gas turbine combustor, known as fuel staging, is an effective way of suppressing instabilities. This dissertation examines the effect of both steady-state and transient fuel-staging in a model gas turbine combustor. First, the steady-state effect of fuel staging on nominally unstable operating conditions is examined. At 200 degrees Celsius and a bulk flow velocity of 26 m/s, the combustor undergoes self-excited instabilities when the equivalence ratio is between φ = 0.70 and φ = 0.74 when all nozzles are fueled equally. These instabilities are suppressed when the center nozzle equivalence ratio is increased to φ Staging = 0.85 and the outer nozzle equivalence ratio maintained at φ = 0.70 or when the global equivalence ratio is maintained at φ = 0.70 by increasing the center nozzle equivalence ratio to φ Staging = 0.74 and decreasing the outer nozzle equivalence ratio to φ Staging = 0.69. The staged p' RMS amplitudes are found to decrease with increasing equivalence ratio and the damping rates found to increase with increasing equivalence ratio. Next, axisymmetric fuel staging transients are conducted valve opening/closing timescales of 1, 16, 58, 4,000, and 10,000 ms, and equivalence ratios of φ staging = 0.75, 0.80, 0.85. It should be noted that these valve opening/closing timescales, referred to as τ Transient , differ from the actual fuel convection timescale that are summarized in Table 5.2. It is found that the combustor pressure fluctuation at the end of the transient strongly depends on staging equivalence ratio but not on τ Transient . Increasing the staging equivalence ratio decreased the staged- stable combustor pressure amplitude, similarly to the steady-state results. The staged-stable damping rate trends are insensitive to to τT Transient and increased with increasing staging equivalence ratio. The timescale of the transition, τ, is found to be sensitive to both τ Transient and the staging equivalence ratio. For a given τ Transient , the decay timescale decreased with increasing staging equivalence ratio. For 1, 16, and 58 ms transients, τ does not depend on τ Transient , though it does strongly depend on direction, where instability onset cases have longer and more variable τ than instability decay cases. The same dependence of τ on direction is seen for the 4,000 and 10,000 ms transients. However, unlike the short duration transients, τ does depend on τ Transient as would be expected. Phase space reconstruction shows the instability is a harmonic limit cycle with a single frequency. The instantaneous phase images shows the center region of the combustor goes out-of-phase first followed by the rest of the combustor. Finally, asymmetric fuel staging tests are conducted. It is found that non-axisymmetric staging can be as effective as axisymmetric staging, although small hardware variations between nozzles can affect the the minimum staging equivalence ratio that consistently suppresses the instability. These changes in the bifurcation equivalence ratio can make staging configurations appear less effective; however, both axisymmetric and non-axisymmetric staging configurations exhibit similar staged-stable p' RMS , damping rates, and instability transition times for equivalence ratios greater than the bifurcation equivalence ratio. Instantaneous phase images show similar instantaneous phase relationship structures in the staged nozzle (center or outer), suggesting that the instability suppression mechanism is similar in each case. This result is surprising given that the time-averaged flame shapes are very different between axisymmetric and non-axisymmetric staging configurations. The similar effectiveness of axisymmetric and non-axisymmetric staging in this geometry is attributed to the fact that the instability mode is a nominally axisymmetric plane wave. It is hypothesized that non-axisymmetric staging configurations will likely be more effective than axisymmetric staging configurations when the instability mode shape is asymmetric (such as a circumferential wave).

Culler, Wyatt↗