Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large scale 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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An automated workflow that generates atom mappings for large‐scale metabolic models and its application to Arabidopsis thaliana

SUMMARY Quantification of reaction fluxes of metabolic networks can help us understand how the integration of different metabolic pathways determines cellular functions. Yet, intracellular fluxes cannot be measured directly but are estimated with metabolic flux analysis (MFA), which relies on the patterns of isotope labeling of metabolites in the network. The application of MFA also requires a stoichiometric model with atom mappings that are currently not available for the majority of large‐scale metabolic network models, particularly of plants. While automated approaches such as the Reaction Decoder Toolkit (RDT) can produce atom mappings for individual reactions, tracing the flow of individual atoms of the entire reactions across a metabolic model remains challenging. Here we establish an automated workflow to obtain reliable atom mappings for large‐scale metabolic models by refining the outcome of RDT, and apply the workflow to metabolic models of Arabidopsis thaliana . We demonstrate the accuracy of RDT through a comparative analysis with atom mappings from a large database of biochemical reactions, MetaCyc. We further show the utility of our automated workflow by simulating 15 N isotope enrichment and identifying nitrogen (N)‐containing metabolites which show enrichment patterns that are informative for flux estimation in future 15 N‐MFA studies of A. thaliana . The automated workflow established in this study can be readily expanded to other species for which metabolic models have been established and the resulting atom mappings will facilitate MFA and graph‐theoretic structural analyses with large‐scale metabolic networks.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Conquering Data Chaos: Research Data Management with Kubernetes

Managing massive volumes of data and effectively making it accessible to researchers poses significant challenges and is a barrier to scientific discovery. In many cases, critical data is locked up in unwieldy file formats or one-off databases and is too large to effectively process on a single machine. This talk explores the role of Kubernetes, an open-source container orchestration platform, in addressing research data management challenges. I will discuss how we are using a set of publicly available open-source and home-grown tools in the National Renewable Energy Lab (NREL) Data, Analysis, and Visualization (DAV) group to help researchers overcome data-related bottlenecks. The talk will begin by providing an overview of the data challenges faced in research data management, including data storage, processing, and analysis. I will highlight Kubernetes' ability to handle large-scale data by leveraging containerization and distributed computing, including distributed storage. Kubernetes allows researchers to encapsulate data processing infrastructure and workflows into portable containers, enabling reproducibility and ease of deployment. Kubernetes can then schedule and manage the resource allocation of these containers to enable efficient utilization of limited computing resources, leading to more efficient data processing and analysis. I will discuss some limitations of traditional, siloed approaches to dealing with data and emphasize the need for solutions which foster collaboration. I will highlight how we are using Kubernetes at NREL to facilitate data sharing and cooperation among research teams. Kubernetes' flexible architecture enables the deployment of shared computing environments, such as Apache Superset, where researchers can seamlessly access and analyze shared datasets. Providing the ability to have one research team easily consume data generated by another, utilizing Kubernetes' as a central data platform, is one of the major wins we've encountered by adopting the platform. Finally, I will showcase real-world use cases from NREL where we have used Kubernetes to solve some persistent data challenges involving large volumes of sensor and monitoring data. I will discuss the challenges we encountered when creating our cluster and making it available as a production-ready resource. I will also discuss the specific suite of tools, including Postgres and Apache Druid for columnar and timeseries data, and Redpanda Kafka for streaming data we have deployed in our infrastructure, and the process that went into the selection of these tools.

collaborative environment↗
DOE OSTI2023

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗
NASA NTRS2004

Scale-Free Networks and Commercial Air Carrier Transportation in the United States

Network science, or the art of describing system structure, may be useful for the analysis and control of large, complex systems. For example, networks exhibiting scale-free structure have been found to be particularly well suited to deal with environmental uncertainty and large demand growth. The National Airspace System may be, at least in part, a scalable network. In fact, the hub-and-spoke structure of the commercial segment of the NAS is an often-cited example of an existing scale-free network After reviewing the nature and attributes of scale-free networks, this assertion is put to the test: is commercial air carrier transportation in the United States well explained by this model? If so, are the positive attributes of these networks, e.g. those of efficiency, flexibility and robustness, fully realized, or could we effect substantial improvement? This paper first outlines attributes of various network types, then looks more closely at the common carrier air transportation network from perspectives of the traveler, the airlines, and Air Traffic Control (ATC). Network models are applied within each paradigm, including discussion of implied strengths and weaknesses of each model. Finally, known limitations of scalable networks are discussed. With an eye towards NAS operations, utilizing the strengths and avoiding the weaknesses of scale-free networks are addressed.

Conway, Sheila R.↗
NASA NTRSDate not supplied

Investigation of the Glen Torridon Clay-Bearing Unit and Overlying Greenheugh Pediment by the Sample Analysis at Mars Instrument Suite

Exploring the Glen Torridon (GT) clay-bearing unit in Gale crater has long been a prime goal for the Mars Science Laboratory (MSL) because of spectral signatures of smectite identified from orbit. Smectite clay minerals can indicate habitable environments [1,2] and possibly facilitate the preservation of organic compounds [e.g., 3-7].Also, sulfate-bearing layers overlie GT and this sequence may record a large-scale change in aqueous conditions overtime [8].The Sample Analysis at Mars (SAM) instrument suitehas been essential in understanding volatile-bearing phases in Gale crater materials. SAM EGA has detected H2O, CO2, O2, H2, SO2, H2S, HCl, NO, and other trace gases, including organic fragments, evolved from many samples on heating. The identity and evolution temperature of evolved gases can support mineral detection by CheMin instrument X-ray diffraction (XRD)and place constraints on trace volatile-bearing phases, clay mineral compositions, and X-ray amorphous materials. In GT, SAM analyzed the sample Kilmarie (KM) from the Jura member, Glen Etive (GE), Mary Anning (MA) and Groken (GR)from the Knockfarril Hill member (KHm), Glasgow (GG) and Hutton (HU) from the fractured Intermediate Unit (fIU), and Edinburgh (EB) from the Stimson formation rocks of the Greenheugh Pediment (GP). The Jura, KHm and fIU are all members of the Murray formation. HU was sampled near the basal Siccar Point group (SPg) unconformity between the fIU and GP, in an area that exhibits a brighter tone and differences in geochemistry compared to Murray materials down-section [9,10]. CheMin XRD of most samples (except for HU and EB)showed large abundances of smectite clay minerals (~25-35%) and all samples showed X-ray amorphous materials and relatively low amounts of well-crystalline hematite[11]. SAM evaluated mineralogy through evolved gas analysis mass spectrometry (EGA) while detailed organic chemistry of several samples was examined by pyrolysis gas chromatography mass spectrometry and wet chemistry experiments [12,13].Here we discuss data and interpretations from SAM EGA analyses.

Mars↗
NASA NTRS2000

Fretting Stresses in Single Crystal Superalloy Turbine Blade Attachments

Single crystal nickel base superalloy turbine blades are being utilized in rocket engine turbopumps and turbine engines because of their superior creep, stress rupture, melt resistance and thermomechanical fatigue capabilities over polycrystalline alloys. Currently the most widely used single crystal nickel base turbine blade superalloys are PWA 1480/1493 and PWA 1484. These alloys play an important role in commercial, military and space propulsion systems. High Cycle Fatigue (HCF) induced failures in aircraft gas turbine and rocket engine turbopump blades is a pervasive problem. Blade attachment regions are prone to fretting fatigue failures. Single crystal nickel base superalloy turbine blades are especially prone to fretting damage because the subsurface shear stresses induced by fretting action at the attachment regions can result in crystallographic initiation and crack growth along octahedral planes. Furthermore, crystallographic crack growth on octahedral planes under fretting induced mixed mode loading can be an order of magnitude faster than under pure mode I loading. This paper presents contact stress evaluation in the attachment region for single crystal turbine blades used in the NASA alternate Advanced High Pressure Fuel Turbo Pump (HPFTP/AT) for the Space Shuttle Main Engine (SSME). Single crystal materials have highly orthotropic properties making the position of the crystal lattice relative to the part geometry a significant factor in the overall analysis. Blades and the attachment region are modeled using a large-scale 3D finite element (FE) model capable of accounting for contact friction, material orthotrophy, and variation in primary and secondary crystal orientation. Contact stress analysis in the blade attachment regions is presented as a function of coefficient of friction and primary and secondary crystal orientation, Stress results are used to discuss fretting fatigue failure analysis of SSME blades. Attachment stresses are seen to reach peak values at locations where fretting cracks have been observed. Fretting stresses at the attachment region are seen to vary significantly as a function of crystal orientation. Attempts to adapt techniques used for estimating fatigue life in the airfoil region, for life calculations in the attachment region, are presented. An effective model for predicting crystallographic crack initiation under mixed mode loading is required for life prediction under fretting action.

Arakere, Nagaraj K.↗
DOE OSTI2022

Influence of anomalous perpendicular transport on linear tearing mode dynamics in tokamak plasmas

The analysis of a previous paper [A. Cole and R. Fitzpatrick, Phys. Plasmas 13, 032503 (2006)] that maps out all of the two-fluid response regimes of a linear tearing layer interacting with an externally generated resonant magnetic perturbation, in a large aspect-ratio tokamak plasma, is generalized to take into account realistic levels of perpendicular particle transport. A new response-regime map is obtained that differs substantially, in many respects, from the old one. The improved analysis is first employed to find all of the two-fluid growth-rate regimes of a non-interacting low-mode-number tearing mode in a large aspect-ratio tokamak plasma. The analysis is then used to determine the scaling of the error-field penetration threshold with machine parameters in large aspect-ratio tokamak plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Boundary layer energization by means of optimized vortex generators

A three-dimensional, multi-block, multi-zone, Euler analysis has been developed and applied to analyze the flow processes induced by a lateral array of low profile vortex generators (VG). These vortex generators have been shown to alleviate boundary layer separation through the generation of streamwise vorticity. The analysis has been applied to help develop improved VG configurations in an efficient manner. Special attention has been paid to determining the accuracy requirements of the solver for calculations in which vortical mechanisms are dominant. The analysis has been used to assess the effectiveness or boundary layer energization capacity of different VG's, including the effect of scale and shape variation. Finally, the analysis has been validated through comparisons with experimental data obtained in a large-scale low-speed wind tunnel.

Barber, T. J.↗
DOE OSTI2022

How galaxies populate halos in extreme density environments: An analysis of the halo occupation distribution in SDSS

Recent works have shown that the properties of galaxy populations in simulated dark matter halos vary with large-scale environments. These results suggest a variation in the halo occupation distribution (HOD) in extreme density environments since the dynamical and astrophysical conditions prevailing in these regions may significantly affect the formation and evolution of their halos and residing galaxies, influencing the mean number of satellite galaxies. To analyse these effects, we identify cosmic voids and future virialised structures (FVS) in the Sloan Digital Sky Server Data Release 12 (SDSS-DR12) and estimate the HOD within these super-structures using group catalogues as dark matter halo proxies. Our goal is to use observational galaxy data to characterise the HOD within cosmic voids and FVS, explore the different properties of these galaxies populations and compare them with the general results outside these super-structures. We use a publicly available observational galaxy catalogue with information on redshifts, positions, magnitudes and other astrophysical features to build a volume complete galaxy sample and identify cosmic voids and FVS. Using a publicly available galaxy group catalogue as a proxy to dark matter halos, we compute the HOD within both types of super-structures for different absolute magnitude thresholds. We also study the dependence on the results on the main void and FVS properties, density and volume. We also analysed the main characteristics of the stellar content of galaxies inside these extreme-density regions such as the mean stellar age and the galaxy light concentration index. In all cases, we compare the results with those derived from the Field sample, defined by objects outside both types of environments. Inside cosmic voids, we find a strong decrease in HOD concerning the Field results. In the most extreme cases, the mean number of satellites fall to ~ 50%. Inside FVS, the HOD shows a significant increase to the Field, with a ~ 40% excess in the mean number of satellites. These results are present for the different galaxy luminosity ranges explored. In both environments, the differences with respect to the Field increases for the extreme values of the density environments. However, we obtain no signs of variations related to intrinsic characteristics of the super-structures, indicating that the effects depend mainly on the density of the large-scale environment. In addition, we find that the cumulative distribution of the mean age of stars of the central galaxy also varies in the different regions, this suggests that the history of the formation of the dark matter halos may be different. Finally, we explore the HOD for the 25% youngest (oldest) galaxies, based on the mean age of their stars. We find that for the low-mass groups the youngest galaxies are only present inside voids and are generally central galaxies. On the other hand, for the high-mass groups the FVS environments show the same increase in the HOD concerning the Field as previously mentioned. We find that cosmic voids lack a significant fraction of galaxies with the oldest stellar population.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1980

Improved subgrid-scale models for large-eddy simulation

The paper analyzes models for subgrid-scale turbulence. The analysis indicates that there is sufficient information in the resolved scales to determine some characteristics of the complete flow field. The kinetic energy of the small-scale motions can be decomposed into two parts: one results from the large scales and is correlated with them, and the other part is uncorrelated which leads to a two-component eddy-viscosity model. The 'production equals dissipation' argument does not apply to the small scales in the decay of turbulence because it does not account for the uncorrelated component. The exchange between the large and small scales takes place mainly between the smallest scales of the former and the largest scales of the latter; this argument is the basis of a new model shown to be superior to the Smagorinsky model (1963).

Bardina, J.↗
NASA NTRS2004

Modeling and Design Analysis Methodology for Tailoring of Aircraft Structures with Composites

Composite materials provide design flexibility in that fiber placement and orientation can be specified and a variety of material forms and manufacturing processes are available. It is possible, therefore, to 'tailor' the structure to a high degree in order to meet specific design requirements in an optimum manner. Common industrial practices, however, have limited the choices designers make. One of the reasons for this is that there is a dearth of conceptual/preliminary design analysis tools specifically devoted to identifying structural concepts for composite airframe structures. Large scale finite element simulations are not suitable for such purposes. The present project has been devoted to creating modeling and design analysis methodology for use in the tailoring process of aircraft structures. Emphasis has been given to creating bend-twist elastic coupling in high aspect ratio wings or other lifting surfaces. The direction of our work was in concert with the overall NASA effort Twenty- First Century Aircraft Technology (TCAT). A multi-disciplinary team was assembled by Dr. Damodar Ambur to work on wing technology, which included our project.

Rehfield, Lawrence W.↗
DOE OSTI2025

A detailed study of pre-heating effects in electron beam melting powder bed fusion process

Metal-based additive manufacturing processes, such as powder bed fusion with electron beam (PBF-EB) process, also referred to as electron beam melting (EBM), can produce high-density parts with minimal residual stresses due to the uniform and coherent preheating of the powder bed. However, understanding and controlling the multiple stages of preheating is required to enable the production of high-quality, consistent parts of various materials. This work presents a large-scale, multi-layer, three-dimensional numerical analysis focused on studying the preheating stages for predicting thermal history during the PBF-EB process. The model follows a continuous multi-stage cyclic process, that incorporates all the main stages of the PBF-EB process for 316 L stainless steel. This includes the gradual deposition of a new powder layer, the first and second preheating levels of the powder bed, and the energy deposition during melting (excluding the actual melt-pool behavior simulation). The model employs an adaptive time-scaling approach that automatically adjusts the energy deposition for each solution time-increment. This allows for localized changes in time-resolution over an otherwise computationally expensive multi-layer procedure. The material property variations are also taken into account, with an emphasis on the subtle irreversible changes in powder effective thermal conductivity after the two requisite preheating stages of the powder bed. This effect is studied using simplified conductivity models from the literature for partially sintered powder, validated by a dedicated experiment and numerical simulation. The large-scale model is then used to estimate the actual temperatures during first and second preheating levels for 316 L steel, which is not yet fully supported commercially for PBF-EB. Model predictions are corroborated by experiments, using and analyzing IR images, taken at the completion of each layer by the machine’s built-in infrared camera. The current model also incorporates a qualitative assessment for the effects of conductivity change during pre-heating, as well as evaluates the applicability of the time-scaling approach.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The Role of the Toroidal Vortex in Cumulus Clouds' Entrainment and Mixing

Abstract Shallow convective clouds play a crucial role in Earth's energy budget, as they modulate the radiative transfer in the atmosphere and participate in the vertical transport of aerosols, energy, and humidity. The parameterizations representing these complex, vital players in weather and climate models are mostly based on a description of steady‐state plumes and are a source of major uncertainty. Recently, several studies have shown that buoyant thermals are inherent in atmospheric convection and contain a toroidal (ring) vortex. This work studies those vortices in growing shallow cumulus (Cu) clouds using high‐resolution (10 m) Large Eddy Simulations that resolve these vortices in much detail. Recent analysis of such data showed that small‐scale turbulent diffusion is unable to explain the large diluted portion of the cloud. Here we advocate for the important role of the Cu toroidal vortex (TV) in cloud dilution and present the complex dynamics and structure of a Cu TV. Nevertheless, since the vortex dominates the cloud's dilution, simplicity emerges when considering the cloud's lateral mass flux profile. The cloud mixing is quantified using direct flux calculations and Eulerian tracers. In addition, Lagrangian tracers are used to identify the origin of the entrained air and its thermodynamic properties. It shows that most of the air entrained by the vortex is not recycled by the vortex, yet is significantly more humid than the environment. We suggest that the development of new models describing thermals, together with their toroidal vortices, might improve cloud parameterizations in weather and climate models.

Meteorology & Atmospheric Sciences↗
NASA NTRS2006

Mechanisms for Diurnal Variability of Global Tropical Rainfall Observed from TRMM

The behavior and various controls of diurnal variability in tropical-subtropical rainfall are investigated using Tropical Rainfall Measuring Mission (TRMM) precipitation measurements retrieved from the three level-2 TRMM standard profile algorithms for the 1998 annual cycle. Results show that diurnal variability characteristics of precipitation are consistent for all three algorithms, providing assurance that TRMM retrievals are producing consistent estimates of rainfall variability. As anticipated, most ocean areas exhibit more rainfall at night, while over most land areas, rainfall peaks during daytime; however, important exceptions are noted. The dominant feature of the oceanic diurnal cycle is a rainfall maximum in late-evening-early-morning (LE-EM) hours, while over land the dominant maximum occurs in the mid- to late afternoon (MLA). In conjunction with these maxima are pronounced seasonal variations of the diurnal amplitudes. Amplitude analysis shows that the diurnal pattern and its seasonal evolution are closely related to the rainfall accumulation pattern and its seasonal evolution. In addition, the horizontal distribution of diurnal variability indicates that for oceanic rainfall, there is a secondary MLA maximum coexisting with the LE-EM maximum at latitudes dominated by large-scale convergence and deep convection. Analogously, there is a preponderancy for an LE EM maximum over land coexisting with the stronger MLA maximum, although it is not evident that this secondary continental feature is closely associated with the large-scale circulation. Neither of the secondary maxima exhibit phase behavior that can be considered semidiurnal in nature. Diurnal rainfall variability over the ocean associated with large-scale convection is clearly an integral component of the general circulation. Phase analysis reveals differences in regional and seasonal features of the diurnal cycle, indicating that underlying forcing mechanisms differ from place to place. This is underscored by the appearance of secondary ocean maxima in the presence of large-scale convection, along with other important features. Among these, there are clear-cut differences between the diurnal variability of seasonal rainfall over the mid-Pacific and Indian Ocean Basins. The mid-Pacific exhibits double maxima in spring and winter but only LE-EM maxima in summer and autumn, while the Indian Ocean exhibits double maxima in spring and summer and only an LE-EM maximum in autumn and winter. There are also evident daytime maxima within the major large-scale marine stratocumulus regions off the west coasts of continents. The study concludes with a discussion concerning how the observational evidence either supports or repudiates possible forcing mechanisms that have been suggested to explain diurnal rainfall variability.

Yang, Song↗
NASA NTRSDate not supplied

Updates of MERRA-2 Data and Services at NASA GES DISC

Over40 years of NASA climate reanalysis datasets from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) are available at the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC). In addition to being used in traditional weather and climate research, MERRA-2 is also widely used in application studies of, e.g., wind and solar energy, air quality and health, food and drought, and heat waves. Two new MERRA-2 datasets were recently added at the GES DISC: (1) climate statistics derived fromMERRA-2 daily data to assist in the analysis of extreme temperature and precipitation events and of large-scale meteorological patterns from 1980 to the present and (2) gridded satellite and conventional observations processed in the MERRA-2 system, along with key statistics derived from the data assimilation, to help better understand how the quality of observations directly affect there analysis data. The GES DISC focuses its efforts on continually improving existing data services and to develop new data tools to satisfy various user communities. The newly added features include the following: Time series service: This is a new service for MERRA-2 data, which enables the easy and fast access of long-term hourly or daily time series at a location for popular parameters. The data is saved in a single le in Ascii format with a user- friendly structure. New analytic functions in the subsetter interface: Options for downloading daily minimum and maximum values have been added into the subsetter interface, in addition to the existing daily mean option, for all MERRA-2 and MERRA sub-daily products. Data format conversion to GeoTIFF has been implemented. New variables in Giovanni: Most monthly variables have been integrated into Giovanni, GESDISC’s online visualization and analysis tool. Due to the large data volume, hourly variables were selected based on user requests. More online information: New MERRA-2 documentation has been added: Data How-To, Data in Action, and FAQ. This presentation overviews two new MERRA-2 datasets and illustrates the new features of data services through a number of case studies. MERRA-2 data and services can be found at: https://disc.gsfc.gov/datasets?

Data management↗
NASA NTRS1985

Coherent structures in a turbulent mixing layer - A comparison between direct numerical simulations and experiments

An eduction scheme has been developed in an attempt to determine the characteristics of large-scale vortical structures in a turbulent mixing layer. This analysis scheme has been applied to a set of experimental data taken in a new, larger mixing layer facility designed to minimize boundary and resonance effects. A similar scheme has been developed to apply to the results of a direct numerical simulation of a temporally growing mixing layer. A comparison of the two approaches shows important similarities in the coherent structures. The numerical simulations indicate that low levels of coherent forcing can dramatically change the evolution of the mixing layer. In the absence of such forcing, the numerical simulations and experiments show a lack of regularity in the transverse position, spacing, amplitude, shape and spanwise coherence of the large-scale vortical structures.

Metcalfe, R. W.↗
NASA NTRS1999

Satellite Data Assimilation

Satellite altimetry provides the only routine observation of a dynamic variable of the global ocean. For instance, measurements from the TOPEX/POSEIDON altimeter have provided dramatic basin-wide images of the space-time evolution of sea level associated with the 1997-98 El Nino event. The nature of large-scale ocean circulation is being studied that underlies such sea level changes measured by satellite altimetry. To first approximation, properties of large-scale (order 1000 km and larger) sea level variability is latitudinally dependent. An analysis using a general circulation model shows that sea level changes are largely due to wind-driven baroclinic (depth dependent) circulation in the tropics (within 20 degrees of the equator), but are primarily due to the expansion and contraction of near surface water in temperate latitudes (between 20 deg. and 40 deg.) forced by seasonal heating and cooling. In contrast, wind-driven barotropic (depth independent) circulation dominates sea level variability in high latitudes (40 deg. and higher) characterized by periods that are as short as a few days. The presence of such inhomogeneity and the significance of high-frequency, large-scale sea level changes had not been fully recognized prior to this study, and are summarized in Fukumori et al. Numerical models provide theoretical relationships among properties that can be inverted using observations so as to estimate the entire state of the ocean, including properties that are otherwise difficult to measure remotely. The process is data assimilation. An approximate Kalman filter and smoother have been devised to assimilate three years of TOPEX/POSEIDON sea level data into a global ocean general circulation model. The figure below demonstrates the skill of the assimilation, and shows the altimeter assimilated estimate being in closer agreement than the simulation is with independent in situ measurements of subsur ace temperature and velocity, consistent with formal uncertainty estimates. Results from this study have been summarized and submitted for publication. The study demonstrates the feasibility of global ocean data assimilation and illustrates applications in monitoring and understanding of processes controlling the evolution of the ocean. Additional information is contained in the original.

Fukumori, Ichiro↗
NASA NTRS2008

Five-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Likelihoods and Parameters from the WMAP Data

The Wilkinson Microwave Anisotropy Probe (WMAP), launched in 2001, has mapped out the Cosmic Microwave Background with unprecedented accuracy over the whole sky. Its observations have led to the establishment of a simple concordance cosmological model for the contents and evolution of the universe, consistent with virtually all other astronomical measurements. The WMAP first-year and three-year data have allowed us to place strong constraints on the parameters describing the ACDM model. a flat universe filled with baryons, cold dark matter, neutrinos. and a cosmological constant. with initial fluctuations described by nearly scale-invariant power law fluctuations, as well as placing limits on extensions to this simple model (Spergel et al. 2003. 2007). With all-sky measurements of the polarization anisotropy (Kogut et al. 2003; Page et al. 2007), two orders of magnitude smaller than the intensity fluctuations. WMAP has not only given us an additional picture of the universe as it transitioned from ionized to neutral at redshift z approx.1100. but also an observation of the later reionization of the universe by the first stars. In this paper we present cosmological constraints from WMAP alone. for both the ACDM model and a set of possible extensions. We also consider tlle consistency of WMAP constraints with other recent astronomical observations. This is one of seven five-year WMAP papers. Hinshaw et al. (2008) describe the data processing and basic results. Hill et al. (2008) present new beam models arid window functions, Gold et al. (2008) describe the emission from Galactic foregrounds, and Wright et al. (2008) the emission from extra-Galactic point sources. The angular power spectra are described in Nolta et al. (2008), and Komatsu et al. (2008) present and interpret cosmological constraints based on combining WMAP with other data. WMAP observations are used to produce full-sky maps of the CMB in five frequency bands centered at 23, 33, 41, 61, and 94 GHz (Hinshaw et al. 2008). With five years of data, we are now able to place better limits on the ACDM model. as well as to move beyond it to test the composition of the universe. details of reionization. sub-dominant components, characteristics of inflation, and primordial fluctuations. We have more than doubled the amount of polarized data used for cosmological analysis. allowing a better measure of the large-scale E-mode signal (Nolta et al. 2008). To this end we describe an alternative way to remove Galactic foregrounds from low resolution polarization maps in which Galactic emission is marginalized over, providing a cross-check of our results. With longer integration we also better probe the second and third acoustic peaks in the temperature angular power spectrum, and have many more year-to-year difference maps available for cross-checking systematic effects (Hinshaw et al. 2008).

Dunkey, J.↗