Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spectral modeling”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1981

Bifurcations from stationary to periodic solutions in a low-order model of forced, dissipative barotropic flow

The considered investigation is concerned with periodic solutions in the context of a forced, dissipative, barotropic spectral model truncated to three complex coefficients with constant forcing on only the intermediate scale. It is found that determining a periodic solution of this three-coefficient model also reduces to finding the algebraic roots of a real polynomial. In the derivation of this polynomial, a class of hydrodynamic spectral systems is described for which a periodic solution might be similarly specified. The existence of periodic solutions of the three-coefficient model is controlled by the roots of the stability polynomial of the basic stationary solution, which represents the simplest response to the constant forcing. When the forcing exceeds a critical value, the basic solution becomes unstable. Owing to the nature of the roots of the stability polynomial at critical forcing, bifurcation theory guarantees the existence of a periodic solution.

Mitchell, K. E.↗
NASA NTRS1991

Empirical models of the latitudinal variations of Te and Ne in the ionosphere at solar maximum

A global spectral model of electron temperature and density is developed in which only the latitudinal variations observed in narrow altitude and local time bands are used. The measurements for this purpose are obtained during intervals that are short as compared to a season. By suppressing all variables except latitude the model can use 17th-order polynomials and generates a latitudinal structure that is similar to that described by the observational data. The treatment is designed to correct for the eccentricity and slow evolution of the data-collection satellite orbit so that small-scale features can be resolved. By comparing the results of a global model and the present latitudinal model to observational measurements of electron temperature and density it is shown that auroral-zone and midlatitude trough structures are not resolved in latitudinal models. The observational data can be used to study the response of electron density and temperature to solar EUV flux variations and geomagnetic activity.

Brace, L. H.↗
DOE OSTI2024

The CXSFIT spectral fitting code: Past, present and future

Magnetically confined plasma experiments generate a wealth of spectroscopic data. The first step toward extracting physical parameters is to fit a spectral model to the often complex spectra. The CXSFIT (Charge eXchange Spectroscopy FITting) spectral fitting code was originally developed for fitting charge exchange spectra on JET from the late 1980s onward and has been further developed over decades to keep up with the needs of the users. The primary use is to efficiently fit a large number of spectra with many constrained Gaussian spectral lines of which the physical parameters can be coupled in a user-friendly manner. More recent additions to the code include time-dependent couplings between parameters, flexible background subtraction, and a non-linear coupling scheme between fit parameters. The latter was a pre-requisite for implementing Zeeman and motional Stark effect multiplets in the library of spectral features. The ability to save and replay “fit recipes,” even when multiple iterations are required, has ensured the traceability of the results and is one of the keys to the longevity and success of the code. The code is also in use on other tokamaks (AUG, ST-40) and to fit data from other spectroscopic diagnostics on JET. In this paper, we document the current capabilities and philosophy behind the structure of the code, including some of the algorithms used to calculate spectral features numerically efficiently. We also provide an outline of how CXSFIT could be transferred into a framework that would be able to meet the spectral fitting requirements of future devices, such as ITER.

Delabie, Ephrem G.↗
NASA NTRS1985

Atmospheric planetary wave response to external forcing

The tools of observational analysis, complex general circulation modeling, and simpler modeling approaches were combined in order to attack problems on the largest spatial scales of the earth's atmosphere. Two different models were developed and applied. The first is a two level, global spectral model which was designed primarily to test the effects of north-south sea surface temperature anomaly (SSTA) gradients between the equatorial and midlatitude north Pacific. The model is nonlinear, contains both radiation and a moisture budget with associated precipitation and surface evaporation, and utilizes a linear balance dynamical framework. Supporting observational analysis of atmospheric planetary waves is briefly summarized. More extensive general circulation models have also been used to consider the problem of the atmosphere's response, especially in the horizontal propagation of planetary scale waves, to SSTA.

Stevens, D. E.↗
NASA NTRS2008

Testing the Paradigm that Ultra-Luminous X-Ray Sources as a Class Represent Accreting Intermediate

To test the idea that ultraluminous X-ray sources (ULXs) in external galaxies represent a class of accreting Intermediate-Mass Black Holes (IMBHs), we have undertaken a program to identify ULXs and a lower luminosity X-ray comparison sample with the highest quality data in the Chandra archive. We establish a general property of ULXs that the most X-ray luminous objects possess the fattest X-ray spectra (in the Chandra band pass). No prior sample studies have established the general hardening of ULX spectra with luminosity. This hardening occurs at the highest luminosities (absorbed luminosity > or equals 5x10(exp 39) ergs/s) and is in line with recent models arguing that ULXs are actually stellar-mass black holes. From spectral modeling, we show that the evidence originally taken to mean that ULXs are IMBHs - i.e., the "simple IMBH model" - is nowhere near as compelling when a large sample of ULXs is looked at properly. During the last couple of years, XMM-Newton spectroscopy of ULXs has to some large extent begun to negate the simple IMBH model based on fewer objects. We confirm and expand these results, which validates the XMM-Newton work in a broader sense with independent X-ray data. We find (1) that cool disk components are present with roughly equal probability and total flux fraction for any given ULX, regardless of luminosity, and (2) that cool disk components extend below the standard ULX luminosity cutoff of 10(exp 39) ergs/s, down to our sample limit of 10(exp 38:3) ergs/s. The fact that cool disk components are not correlated with luminosity damages the argument that cool disks indicate IMBHs in ULXs, for which a strong statistical support was never made.

Berghea, C. T.↗
DOE OSTI2020

Laser Plasma Interaction in the MG Magnetic Field (Final Report)

Strong magnetic fields change properties of plasma. Employment of the magnetic fields can improve conditions for inertial confinement fusion. Compression and heating of the magnetized plasma is a base of the MagLIF pulsed power approach to the controlled fusion. Study of expansion and heating of plasma in the external magnetic fields, development of plasma instabilities, parametric and other effects are important for basic physics and applications. The main goal of the research for the grant was experimental investigation of laser produced plasma in the well characterized and controlled MG magnetic fields. We developed a platform for investigation of high intensity laser-plasma interaction (LPI) in the MG fields using a university-scale pulsed power machine. A Zebra machine at UNR generated transverse magnetic fields of 3 MG on the surface of the rod load and longitudinal fields of 1.4 MG in the coil loads. The magnetic fields were characterized with two-color Faraday rotation diagnostics. The pulsed power machine provided almost static magnetic fields for LPI compared to the much faster laser driven magnetic fields. We found that eddy currents in the rising magnetic field generate plasma on the metal targets. Eddy currents did not impact Si and CH dielectric targets. A Leopard laser with intensity in the focus up to 10 19 W/cm 2 produced plasma in the MG magnetic fields. The modified reflected Faraday diagnostics was used to study laser driven B-field at the MTW laser at Laboratory for Laser Energetics, UR. We demonstrated a slow dependence of the magnetic field in the coil on the applied laser intensity. The magnetic field increased by a factor of 2 if laser intensity increased by a factor of 30. This scaling is important for comparison of experiments at different laser facilities. Expansion of the plasma in the 2-3 MG external B-field was studied. Laser-produced plasma in the azimuthal magnetic field took the unique form of a thin disc expanding radially with a velocity of 250 km/s and confined in the vertical direction. A HYDRA MHD program at the University of Rochester was applied for simulations of LPI in the magnetic field B >3 MG. Simulations for the plasma disc parameters and the expansion dynamics were in a good agreement with experiments. Generation of narrow dense plasma jets in the longitudinal magnetic field of 0.6-0.8 MG was studied. Narrow plasma jets reached a length of 3-4 mm with the electron plasma density of (0.2-1.2)x10 20 cm -3 . A jet tip propagated with the velocity of 160-200 km/s. MHD simulations showed a good agreement of the dynamics of the formation of plasma jets with experiments. These jets are relevant to astrophysical jets. Simulations of K-shell spectra of Si plasma was performed for LPI in the B-field. The laser prepulse generated plasma near the laser target. A PrismSPECT model with a MeV electron beam produced by the laser was used. The spectral modelling showed the increased plasma density by a factor by 2-3 in the magnetic field. This was in agreement with the observed dynamics of the plasma plume confined by the axial magnetic field. The two-plasmon decay was studied in the 2.5-3 MG transversal magnetic fields. Strong 2-3.5 nm spectral widening and a 2-4 nm shifts of “red” and “blue” 3/2ω 0 spectral components were observed. The large shift and widening exceed the expected temperature and magnetic shift. PIC simulations are performing to clarify physics of this effect. The experimental research program for the grant is completed. New plasma effects in the MG fields were studied at the university-scale pulsed power generator. Two graduate and two undergraduate UNR and UR students carried experiments at the Zebra generator, supported and developed plasma diagnostics, provided data processing, and performed MHD simulations of LPI in the strong magnetic field. Research results are published in 9 referred papers and presented in 12 conferences. Program manager: Kramer U. Akli; Collaborators: R. C. Mancini and H. Sawada, University of Nevada, Reno (UNR); R. Betti and A. V. Maximov, University of Rochester (UR).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Modeling the Dynamics of the Middle Atmosphere and Lower Thermosphere Under the Influence of Gravity Waves

Our Numerical Spectral Model (NSM), which extends from the ground up into the thermosphere, is non-linear, time-dependent and has been employed for 2D and 3D applications. The standard version of the NSM incorporates Hines' Doppler Spread Parameterization for small scale gravity waves (GW), but planetary waves generated in the troposphere have also been incorporated. The NSM has been applied to describe: (1) the anomalous seasonal variations of the zonal circulation and temperature in the upper mesosphere, (2) the equatorial oscillations (quasi-biennial and semi-annual oscillations (QBO and SAO)) extending from the stratosphere into the upper mesosphere, (3) the diurnal and semi-diurnal tides, and (4) the planetary waves that are excited in the mesosphere. With the emphasis to provide understanding, we present here results from numerical experiments with the NSM that shed light on the GW processes that are of central importance in the mesosphere and lower thermosphere. These are our conclusions: (1) The large semiannual variations in the diurnal tide (DT), with peak amplitudes observed around equinox, are produced primarily by GW interactions that involve, in part, planetary waves. The DT, like planetary waves, tends to be amplified by GW momentum deposition, which reduces also the vertical wavelength, but variations in eddy viscosity associated with GW interactions are also important. (2) The semidiurnal tide (SDT) and its phase in particular, is strongly influenced by the mean zonal circulation. The SDT, individually, is also amplified by GW. But the DT filters out GW such that the GW interaction effectively reduces the amplitude of the SDT, producing a strong nonlinear interaction between the DT and SDT. (3) Without external time dependent energy or momentum sources, planetary waves (PW) are generated in the model for zonal wavenumbers 1 to 4, which have amplitudes in the mesosphere above 50 km as large as 40 m/s and periods between 50 and 2 days. The waves are generated primarily during solstice conditions, which indicates that the baroclinic instability (associated with the GW induced reversal in the latitudinal temperature gradient) is playing an important role. Numerical experiment show that GW, directly, also greatly amplify the PW. A common feature of the PW generated in summer and winter is that their vertical wavelengths throughout the mesosphere are large, which indicates that the waves are not propagating freely but are generated throughout the region. Another common feature is that the PW propagate preferentially westward in summer and eastward in winter, being launched from the westward and eastward zonal winds that prevail respectively in summer and winter at altitudes below 80 km. (4) Planetary waves generated internally by baroclinic instability and GW interaction produce large amplitude modulations of the DT and SDT. In summary we conclude that GW play major roles in generating and amplifying the dynamical components in the MLT region and, acting principally through wave filtering, produce important non-linear interactions between the components.

Mayr, H. G.↗
DOE OSTI2023

Data and scripts associated with a manuscript investigating impacts of solid phase extraction on freshwater organic matter optical signatures and mass spectrometry pairing

This data package is associated with the publication “Investigating the impacts of solid phase extraction on dissolved organic matter optical signatures and the pairing with high-resolution mass spectrometry data in a freshwater system” submitted to “Limnology and Oceanography: Methods.” This data is an extension of the River Corridor and Watershed Biogeochemistry SFA’s Spatial Study 2021 (https://doi.org/10.15485/1898914). Other associated data and field metadata can be found at the link provided. The goal of this manuscript is to assess the impact of solid phase extraction (SPE) on the ability to pair ultra-high resolution mass spectrometry data collected from SPE extracts with optical properties collected on ambient stream samples. Forty-seven samples collected from within the Yakima River Basin, Washington were analyzed dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC), absorbance, and fluorescence. Samples were subsequently concentrated with SPE and reanalyzed for each measurement. The extraction efficiency for the DOC and common optical indices were calculated. In addition, SPE samples were subject to ultra-high resolution mass spectrometry and compared with the ambient and SPE generated optical data. Finally, in addition to this cross-platform inter-comparison, we further performed and intra-comparison among the high-resolution mass spectrometry data to determine the impact of sample preparation on the interpretability of results. Here, the SPE samples were prepared at 40 milligrams per liter (mg/L) based on the known DOC extraction efficiency of the samples (ranging from ~30 to ~75%) compared to the common practice of assuming the DOC extraction efficiency of freshwater samples at 60%. This data package folder consists of one main data folder with one subfolder (Data_Input). The main data folder contains (1) readme; (2) data dictionary (dd); (3) file-level metadata (flmd); (4) final data summary output from processing script; and (5) the processing script. The R-markdown processing script (SPE_Manuscript_Rmarkdown_Data_Package.rmd) contains all code needed to reproduce manuscript statistics and figures (with the exception of that stated below). The Data_Input folder has two subfolders: (1) FTICR and (2) Optics. Additionally, the Data_Input folder contains dissolved organic carbon (DOC, measured as non-purgeable organic carbon, NPOC) data (SPS_NPOC_Summary.csv) and relevant supporting Solid Phase Extraction Volume information (SPS_SPE_Volumes.csv). Methods information for the optical and FTICR data is embedded in the header rows of SPS_EEMs_Methods.csv and SPS_FTICR_Methods.csv, respectively. In addition, the data dictionary (SPS_SPE_dd.csv), file level metadata (SPS_SPE_flmd.csv), and methods codes (SPS_SPE_Methods_codes.csv) are provided. The FTICR subfolder contains all raw FTICR data as well as instructions for processing. In addition, post processed FTICR molecular information (Processed_FTICRMS_Mol.csv) and sample data (Processed_FTICRMS_Data.csv) is provided that can be directly read into R with the associated R-markdown file. The Optics subfolder contains all Absorbance and Fluorescence Spectra. Fluorescence spectra have been blank corrected, inner filter corrected, and undergone scatter removal. In addition, this folder contains Matlab code used to make a portion of Figure 1 within the manuscript, derive various spectral parameters used within the manuscript, and used for parallel factor analysis (PARAFAC) modeling. Spectral indices (SPS_SpectralIndices.csv) and PARAFAC outputs (SPS_PARAFAC_Model_Loadings.csv and SPS_PARAFAC_Sample_Scores.csv) are directly read into the associated R-markdown file. We acknowledge the Yakama Nation as owners and caretakers of the lands where we collected some of these data. We thank the Confederated Tribes and Bands of the Yakama Nation Tribal Council and Yakama Nation Fisheries for working with us to facilitate sample collection and optimization of data usage according to their values and worldview.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1989

Land-surface temperature measurement from space - Physical principles and inverse modeling

To apply the multiple-wavelength (split-window) method used for satellite measurement of sea-surface temperature from thermal-infrared data to land-surface temperatures, the authors statistically analyze simulations using an atmospheric radiative transfer model. The range of atmospheric conditions and surface temperatures simulated is wide enough to cover variations in clear atmospheric properties and surface temperatures, both of which are larger over land than over sea. Surface elevation is also included in the simulation as the most important topographic effect. Land covers characterized by measured or modeled spectral emissivities include snow, clay, sands, and tree leaf samples. The empirical inverse model can estimate the surface temperature with a standard deviation less than 0.3 K and a maximum error less than 1 K, for viewing angles up to 40 degrees from nadir under cloud-free conditions, given satellite measurements in three infrared channels. A band in the region from 10.2 to 11.0 microns will usually give the most reliable single-band estimate of surface temperature. In addition, a band in either the 3.5-4.0-micron region or in the 11.5-12.6-micron region must be included for accurate atmospheric correction, and a band below the ozone absorption feature at 9.6 microns (e.g., 8.2-8.8 microns) will increase the accuracy of the estimate of surface temperature.

Wan, Zhengming↗
NASA NTRS1989

Explosive east coast cyclogenesis - Numerical experimentation and model-based diagnostics

Numerical experimentation of explosive east-coast cyclogenesis is performed using the Florida State University Global Spectral Model (FSUGSM). The three cases examined here are the Presidents' Day storm of February 18-19, 1979 and the North Atlantic and Pacific bombs of January 18-20, 1979 which formed off the east coasts of the United States and Japan, respectively. The use of a global model provides a framework for studying the phenomena on the 3-5 day time scale. The forecast verifications of the numerical experiments indicate that the FSUGSM was able to adequately predict the phase, intensity, and synoptic-scale structure. These results justify the use of model data for diagnostic studies of the bomb. The model data are used to quantify the role of the adiabatic and diabatic forcing in the explosive cyclogenetic process, using surface pressure tendency to gage development.

Manobianco, John↗
DOE OSTI2024

SCEC/USGS Community Stress-Drop Validation Study: How Spectral Fitting Approaches Influence Measured Source Parameters

Spectral source parameters used to estimate an earthquake’s stress drop (⁠Δσ⁠) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress‐Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low‐frequency displacement (⁠Ω 0 ⁠) and corner frequency (⁠ƒ c ⁠) by fitting a smooth function to benchmark displacement spectra. Participants used linear‐ or log‐sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate ⁠Ω 0 ⁠, ƒ c ⁠, and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in and can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade‐off between Ω 0 and ƒ c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω 0 and somewhat higher ƒ c compared to those assuming a Brune model, although resulting Δσ estimates are similar. Finally, when compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ⁠.

58 GEOSCIENCES↗
DOE OSTI2022

The MASSIVE Survey. XVI. The Stellar Initial Mass Function in the Center of MASSIVE Early-type Galaxies

The stellar initial mass function (IMF) is a fundamental property in the measurement of stellar masses and galaxy star formation histories. In this work, we focus on the most massive galaxies in the nearby universe $\mathrm{log}({M}_{\star }/\,{M}_{\odot })\gt 11.2$. We obtain high-quality Magellan/LDSS-3 long-slit spectroscopy with a wide wavelength coverage of 0.4–1.01 μm for 41 early-type galaxies (ETGs) in the MASSIVE survey and derive high signal-to-noise spectra within an aperture of R e /8. Using detailed stellar synthesis models, we constrain the elemental abundances and stellar IMF of each galaxy through full spectral modeling. All the ETGs in our sample have an IMF that is steeper than a Milky Way (Kroupa) IMF. The best-fit IMF mismatch parameter, α IMF = (M/L)/(M/L) MW , ranges from 1.1 to 3.1, with an average of $\langle$α IMF $\rangle$ = 1.84, suggesting that on average, the IMF is more bottom heavy than Salpeter. Comparing the estimated stellar masses with the dynamical masses, we find that most galaxies have stellar masses that are smaller than their dynamical masses within the 1σ uncertainty. We complement our sample with lower-mass galaxies from the literature and confirm that $\mathrm{log}({\alpha }_{\mathrm{IMF}})$ is positively correlated with $\mathrm{log}(\sigma )$, $\mathrm{log}({M}_{\star })$, and $\mathrm{log}({M}_{\mathrm{dyn}})$. From the combined sample, we show that the IMF in the centers of more massive ETGs is more bottom heavy. In addition, we find that $\mathrm{log}({\alpha }_{\mathrm{IMF}})$ is positively correlated with both [Mg/Fe] and the estimated total metallicity [Z/H]. We find suggestive evidence that the effective stellar surface density Σ Kroupa might be responsible for the variation of α IMF . We conclude that σ, [Mg/Fe], and [Z/H] are the primary drivers of the global stellar IMF variation.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1986

A global three-dimensional model of the circulation and chemistry of CFCl3, CF2Cl2, CH3CCl3, CCl4, and N2O

The use of a three-dimensional spectral model to study the tropospheric and stratospheric circulation, chemistry, and photochemistry of the CFCl3, CF2Cl2, CH3CCl3, CCl4, and N2O atmospheric species is examined. The components of the model are described. Lifetime, regional, and global trends and budgets for the species are evaluated. Calculated horizontal, vertical, and temporal distributions of the atmospheric species are compared with observations; good correlation is detected. The differences observed between calculated and observed surface distributions of CH3CCl3 and the vertical distribution of CCl4 are analyzed. The calculated global atmospheric lifetimes of CFCl3, CF2Cl2, CCl4, and N2O are 73, 210, 48, and 182 years, respectively.

Golombek, A.↗
NASA NTRS2001

Interannual Variability of Tropical Precipitation: How Well Do Climate Models Agree With Current Satellite Estimates?

Since the beginning of the World Climate Research Program's Global Precipitation Climatology Project (GPCP) satellite remote sensing of precipitation has made dramatic improvements, particularly for tropical regions. Data from microwave and infrared sensors now form the most critical input to precipitation data sets and can be calibrated with surface gauges to so that the strengths of each data source can be maximized in some statistically optimal sense. Recent availability of the TRMM (Tropical Rainfall Measuring Mission) has further aided in narrowing uncertainties in rainfall over die tropics and subtropics. Although climate modeling efforts have long relied on space-based precipitation estimates for validation, we now are in a position to make more quantitative assessments of model performance, particularly in tropical regions. An integration of the CCM3 using observed SSTs as a lower boundary condition is used to examine how well this model responds to ENSO forcing in terms of anomalous precipitation. An integration of the NCEP spectral model used for the Reanalysis-H effort is also examined. This integration is run with specified SSTs, but with no data assimilation. Our analysis focuses on two aspects of inter-annual variability. First are the spatial anomalies that are indicative of dislocations in Hadley and Walker circulations. Second, we consider the ability of models to replicate observed increases in oceanic precipitation that are noted in satellite observations for large ENSO events. Finally, we consider a slab ocean version of the CCM3 model with prescribed ocean beat transports that mimic upwelling anomalies, but which still allows the surface energy balance to be predicted. This less restrictive experiment is used to understand why model experiments with specified SSTs seem to have noticeably less interannual variability in precipitation than do the satellite observations.

Robertson, Franklin R.↗
NASA NTRS2000

Interannual Variability of the Tropical Energy Balance: Reconciling Observations and Models

Since the beginning of the World Climate Research Program's Global Precipitation Climatology Project (GPCP) satellite remote sensing of precipitation has made dramatic improvements, particularly for tropical regions. Data from microwave and infrared sensors now form the most critical input to precipitation data sets and can be calibrated with surface gauges to so that the strengths of each data source can be maximized in some statistically optimal sense. Recent availability of the TRMM (Tropical Rainfall Measuring Mission) has further aided in narrowing uncertainties in rainfall over the tropics and subtropics. Although climate modeling efforts have long relied on space-based precipitation estimates for validation, we now are in a position to make more quantitative assessments of model performance, particularly in tropical regions. An integration of the CCM3 using observed SSTs as a lower boundary condition is used to examine how well this model responds to ENSO forcing in terms of anomalous precipitation. An integration of the NCEP spectral model used for the Reanalysis-11 effort is also examined. This integration is run with specified SSTs, but no data assimilation. Our analysis focuses on two aspects. First are the spatial anomalies that are indicative of dislocations in Hadley and Walker circulations. Second, we consider the ability of models to replicate observed increases in oceanic precipitation that are noted in satellite observations for large ENSO events. Finally, we consider a slab ocean version of the CCM3 model with prescribed ocean heat transports that mimic upwelling anomalies, but which still allows the surface energy balance to be predicted. This less restrictive experiment is used to understand why model experiments with specified SSTs seem to have noticeably less interannual variability than do the satellite precipitation observations.

Robertson, Franklin R.↗
NASA NTRS1995

Enhancements to the SSME transfer function modeling code

This report details the results of a one year effort by Ohio University to apply the transfer function modeling and analysis tools developed under NASA Grant NAG8-167 (Irwin, 1992), (Bartholomew, 1992) to attempt the generation of Space Shuttle Main Engine High Pressure Turbopump transfer functions from time domain data. In addition, new enhancements to the transfer function modeling codes which enhance the code functionality are presented, along with some ideas for improved modeling methods and future work. Section 2 contains a review of the analytical background used to generate transfer functions with the SSME transfer function modeling software. Section 2.1 presents the 'ratio method' developed for obtaining models of systems that are subject to single unmeasured excitation sources and have two or more measured output signals. Since most of the models developed during the investigation use the Eigensystem Realization Algorithm (ERA) for model generation, Section 2.2 presents an introduction of ERA, and Section 2.3 describes how it can be used to model spectral quantities. Section 2.4 details the Residue Identification Algorithm (RID) including the use of Constrained Least Squares (CLS) and Total Least Squares (TLS). Most of this information can be found in the report (and is repeated for convenience). Section 3 chronicles the effort of applying the SSME transfer function modeling codes to the a51p394.dat and a51p1294.dat time data files to generate transfer functions from the unmeasured input to the 129.4 degree sensor output. Included are transfer function modeling attempts using five methods. The first method is a direct application of the SSME codes to the data files and the second method uses the underlying trends in the spectral density estimates to form transfer function models with less clustering of poles and zeros than the models obtained by the direct method. In the third approach, the time data is low pass filtered prior to the modeling process in an effort to filter out high frequency characteristics. The fourth method removes the presumed system excitation and its harmonics in order to investigate the effects of the excitation on the modeling process. The fifth method is an attempt to apply constrained RID to obtain better transfer functions through more accurate modeling over certain frequency ranges. Section 4 presents some new C main files which were created to round out the functionality of the existing SSME transfer function modeling code. It is now possible to go from time data to transfer function models using only the C codes; it is not necessary to rely on external software. The new C main files and instructions for their use are included. Section 5 presents current and future enhancements to the XPLOT graphics program which was delivered with the initial software. Several new features which have been added to the program are detailed in the first part of this section. The remainder of Section 5 then lists some possible features which may be added in the future. Section 6 contains the conclusion section of this report. Section 6.1 is an overview of the work including a summary and observations relating to finding transfer functions with the SSME code. Section 6.2 contains information relating to future work on the project.

Irwin, R. Dennis↗
NASA NTRS2009

Intra-Seasonal Monthly Oscillations in Stratospheric NCEP Data and Model Results

Intra-seasonal oscillations (ISO) are observed in the zonal-mean of mesospheric wind and temperature measurements-and the numerical spectral model (NSM) generates such oscillations. Relatively large temperature ISO are evident also in stratospheric CPC (NCEP) data at high latitudes, where the NSM produces amplitudes around 3 K at 30 km. Analyzing the NCEP data for the years 1996-2006, we find in Fourier spectra signatures of oscillations with periods between 1.7 and 3 months. With statistical confidence levels exceeding 70%, the spectral features are induced by nonlinear interactions involving the annual and semi-annual variations. The synthesized data show for the 10-year average that the temperature ISO peak in winter, having amplitudes close to 4 K. The synthesized complete spectrum for periods around 2 months produces oscillations, varying from year to year, which can reach peak amplitudes of 15 and 5 K respectively at northern and southern polar latitudes.

Mayr, H. G.↗