Search NASA⌕ Search

SEARCH · Search NASA

Results for “Grid Interactions 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 253 records · Page 14

Comparison of Observed Longwave, Shortwave Irradiance and Surface Temperature from “MOSAiC” to CERES Radiative Transfer Calculations and Inputs

As atmospheric temperatures rise due to increased anthropogenic forcing, the effect is expected to be larger over the arctic than midlatitude and tropics, known as polar amplification. A multi-national program, the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) was ran between Sept 2019 and Oct 2020. The project deployed scientific instruments on board the German research vessel Polarstern with the ship remaining across a year to observe all aspects of polar climate. The US DOE deployed the ARM AMF2 aboard the ship along with several off-ship sites for extended spatial observations. Of particular importance was the measurement of the energetics of the ice/ocean/atmosphere interactions through observations of surface irradiance and surface temperatures. NASA’s Clouds and the Earth’s Radiant Energy System (CERES) project produces the SYN1deg Edition 4 data product. This product provides high quality, hourly, globally gridded and temporally complete maps of top of atmosphere (TOA), in atmosphere, and surface irradiances. TOA fluxes are derived from CERES instruments and geostationary satellites. In addition, TOA, in atmosphere, and surface irradiances are computed using the Langley Fu and Liou radiative transfer model. The radiative transfer model is run hourly at 1 degree ~equal area spatial resolution. Meteorological profiles are provided by Global Modeling and Assimilation Office’s GEOS-541 reanalysis product. Cloud properties are derived solely from Terra and Aqua MODIS imagers northward of 60°. Here we compare SYN1deg hourly calculations of surface irradiance to observations from the AMF2 to validate the products estimates of surface irradiance in this challenging area. Along with the irradiance comparisons we take a close look at the surface temperature record in the GEOS-541 product and compare it over both time and space to several surface observations provided by MOSAiC. Along with the comparison to the re-analysis record we will compare these observed surface temperatures to those derived from the AIRS product, which is known to have some difficulty in extreme high latitude areas. The goals of this study are 1) understand computed surface downward irradiance and temperature error separated by surface type (e.g. sea ice or open water) and by season, and 2) error covariance in spatial and temporal space. We seek to use this information to extrapolate the error from the MOSAiC domain to a larger arctic region.

David A Rutan↗

GEOS S2S-2_1 File Specification: GMAO Seasonal and Sub-Seasonal Forecast Output

The NASA GMAO seasonal (9 months) and subseasonal (45 days) forecasts are produced with the Goddard Earth Observing System (GEOS) Atmosphere-Ocean General Circulation Model and Data Assimilation System Version S2S-2_1. The new system replaces version S2S-1.0 described in Borovikov et al (2017), and includes upgrades to many components of the system. The atmospheric model includes an upgrade from a pre-MERRA-2 version running on a latitude-longitude grid at approx. 1 degree resolution to a current version running on a cubed sphere grid at approximately 1/2 degree resolution. The important developments are related to the dynamical core (Putman et al., 2011), the moist physics (''two-moment microphysics'' of Barahona et al., 2014) and the cryosphere (Cullather et al., 2014). As in the previous GMAO S2S system, the land model is that of Koster et al (2000). GMAO S2S-2_1 now includes the Goddard Chemistry Aerosol Radiation and Transport (GOCART, Colarco et al., 2010) single moment interactive aerosol model that includes predictive aerosols including dust, sea salt and several species of carbon and sulfate. The previous version of GMAO S2S specified aerosol amounts from climatology, which were used to inform the atmospheric radiation only. The ocean model includes an upgrade from MOM4 to MOM5 (Griffies 2012), and continues to be run on the tripolar grid at approximately 1/2 degree resolution in the tropics with 40 vertical levels. As in S2S-1.0, the sea ice model is from the Los Alamos Sea Ice model (CICE4, Hunke and Lipscomb 2010). The Ocean Data Assimilation System (ODAS) has been upgraded from the one described in Borovikov et al., 2017 to one that uses a modified version of the Penny, 2014 Local Ensemble Transform Kalman Filter (LETKF), and now assimilates along-track altimetry. The ODAS also does a nudging to MERRA-2 SST and sea ice boundary conditions. The atmospheric data assimilation fields used to constrain the atmosphere in the ODAS have been upgraded from MERRA to a MERRA-2 like system. The system is initialized using a MERRA-2-like atmospheric reanalysis (Gelaro et al. 2017) and the GMAO S2S-2_1 ocean analysis. Additional ensemble members for forecasts are produced with initial states at 5-day intervals, with additional members based on perturbations of the atmospheric and ocean states. Both subseasonal and seasonal forecasts are submitted to the National MultiModel Ensemble (NMME) project, and are part of the US/Canada multimodel seasonal forecasts (http://www.cpc.ncep.noaa.gov/products/NMME/). A large suite of retrospective forecasts (''hindcasts'') have been completed, and contribute to the calculation of the model's baseline climatology and drift, anomalies from which are the basis of the seasonal forecasts.

Seasonal Forecast↗

The Impacts of an Observationally-Based Cloud Fraction and Condensate Overlap Parameterization on a GCM's Cloud Radiative Effect

It has been shown that the details of how cloud fraction overlap is treated in GCMs has substantial impact on shortwave and longwave fluxes. Because cloud condensate is also horizontally heterogeneous at GCM grid scales, another aspect of cloud overlap should in principle also be assessed, namely the vertical overlap of hydrometeor distributions. This type of overlap is usually examined in terms of rank correlations, i.e., linear correlations between hydrometeor amount ranks of the overlapping parts of cloud layers at specific separation distances. The cloud fraction overlap parameter and the rank correlation of hydrometeor amounts can be both expressed as inverse exponential functions of separation distance characterized by their respective decorrelation lengths (e-folding distances). Larger decorrelation lengths mean that hydrometeor fractions and probability distribution functions have high levels of vertical alignment. An analysis of CloudSat and CALIPSO data reveals that the two aspects of cloud overlap are related and their respective decorrelation lengths have a distinct dependence on latitude that can be parameterized and included in a GCM. In our presentation we will contrast the Cloud Radiative Effect (CRE) of the GEOS-5 atmospheric GCM (AGCM) when the observationally-based parameterization of decorrelation lengths is used to represent overlap versus the simpler cases of maximum-random overlap and globally constant decorrelation lengths. The effects of specific overlap representations will be examined for both diagnostic and interactive radiation runs in GEOS-5 and comparisons will be made with observed CREs from CERES and CloudSat (2B-FLXHR product). Since the radiative effects of overlap depend on the cloud property distributions of the AGCM, the availability of two different cloud schemes in GEOS-5 will give us the opportunity to assess a wide range of potential cloud overlap consequences on the model's climate.

Oreopoulos, Lazaros↗

An analysis of supersonic flows with low-Reynolds number compressible two-equation turbulence models using LU finite volume implicit numerical techniques

A generalized flow solver using an implicit Lower-upper (LU) diagonal decomposition based numerical technique has been coupled with three low-Reynolds number kappa-epsilon models for analysis of problems with engineering applications. The feasibility of using the LU technique to obtain efficient solutions to supersonic problems using the kappa-epsilon model has been demonstrated. The flow solver is then used to explore limitations and convergence characteristics of several popular two equation turbulence models. Several changes to the LU solver have been made to improve the efficiency of turbulent flow predictions. In general, the low-Reynolds number kappa-epsilon models are easier to implement than the models with wall-functions, but require much finer near-wall grid to accurately resolve the physics. The three kappa-epsilon models use different approaches to characterize the near wall regions of the flow. Therefore, the limitations imposed by the near wall characteristics have been carefully resolved. The convergence characteristics of a particular model using a given numerical technique are also an important, but most often overlooked, aspect of turbulence model predictions. It is found that some convergence characteristics could be sacrificed for more accurate near-wall prediction. However, even this gain in accuracy is not sufficient to model the effects of an external pressure gradient imposed by a shock-wave/ boundary-layer interaction. Additional work on turbulence models, especially for compressibility, is required since the solutions obtained with base line turbulence are in only reasonable agreement with the experimental data for the viscous interaction problems.

Lee, J.↗

MFLOP to GFLOP: The Impact on High Fidelity Based Computational Aeroelasticity

Aeroelasticity which involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of modern aircraft which can experience complex flow/structure interactions. Even at moderate angles of attack supersonic aircraft can experience vortex induced aeroelastic oscillations. Near transonic speeds buffet associated structural oscillations are possible. Aircraft flying in transonic regime may experience a dip in the flutter speed. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Computations using these high fidelity equations require large computational resources both in memory and speed. Current conventional supercomputers have reached their limitations both in memory and speed. As a result, parallel computers have evolved to overco me the limitations of conventional computers. This paper will address the transition that is taking place in computational aeroelasticity from conventional computers to parallel computers. The paper will address special techniques needed to take advantage of the architecture of new parallel computers. Results will be illustrated from computations made on iPSC/860 and IBM SP2 computer by using ENSAERO code that directly couples the Euler/Navier-Stokes flow equations with high resolution finite-element structural equations. Modifications required in both fluids and structural solvers in order to run efficiently on parallel computers will be discussed. Implementation of moving grids and fluid/structural interface on parallel computers will be discussed.

Guruswamy, Guru P.↗

The GGCMI Phase 2 emulators: global gridded crop model responses to changes in CO2, temperature, water, and nitrogen (version 1.0)

Statistical emulation allows combining advantageous features of statistical and process-based crop models for understanding the effects of future climate changes on crop yields. We describe here the development of emulators for nine process-based crop models and five crops using output from the Global Gridded Model Intercomparison Project (GGCMI) Phase 2. The GGCMI Phase 2 experiment is designed with the explicit goal of producing a structured training dataset for emulator development that samples across four dimensions relevant to crop yields: atmospheric carbon dioxide (CO2) concentrations, temperature, water supply, and nitrogen inputs (CTWN). Simulations are run under two different adaptation assumptions: that growing seasons shorten in warmer climates, and that cultivar choice allows growing seasons to remain fixed. The dataset allows emulating the climatological-mean yield response of all models with a simple polynomial in mean growing-season values. Climatological-mean yields are a central metric in climate change impact analysis; we show here that they can be captured without relying on interannual variations. In general, emulation errors are negligible relative to differences across crop models or even across climate model scenarios; errors become significant only in some marginal lands where crops are not currently grown. We demonstrate that the resulting GGCMI emulators can reproduce yields under realistic future climate simulations, even though the GGCMI Phase 2 dataset is constructed with uniform CTWN offsets, suggesting that the effects of changes in temperature and precipitation distributions are small relative to those of changing means. The resulting emulators therefore capture relevant crop model responses in a lightweight, computationally tractable form, providing a tool that can facilitate model comparison, diagnosis of interacting factors affecting yields, and integrated assessment of climate impacts.

James A Franke↗

Interannual and Decadal Variability of Ocean Surface Latent Heat Flux as Seen from Passive Microwave Satellite Algorithms

Ocean surface turbulent fluxes are critical links in the climate system since they mediate energy exchange between the two fluid systems (ocean and atmosphere) whose combined heat transport determines the basic character of Earth's climate. Deriving physically-based latent and sensible heat fluxes from satellite is dependent on inferences of near surface moisture and temperature from coarser layer retrievals or satellite radiances. Uncertainties in these "retrievals" propagate through bulk aerodynamic algorithms, interacting as well with error properties of surface wind speed, also provided by satellite. By systematically evaluating an array of passive microwave satellite algorithms, the SEAFLUX project is providing improved understanding of these errors and finding pathways for reducing or eliminating them. In this study we focus on evaluating the interannual variability of several passive microwave-based estimates of latent heat flux starting from monthly mean gridded data. The algorithms considered range from those based essentially on SSM/I (e.g. HOAPS) to newer approaches that consider additional moisture information from SSM/T-2 or AMSU-B and lower tropospheric temperature data from AMSU-A. On interannual scales, variability arising from ENSO events and time-lagged responses of ocean turbulent and radiative fluxes in other ocean basins (as well as the extratropical Pacific) is widely recognized, but still not well quantified. Locally, these flux anomalies are of order 10-20 W/sq m and present a relevant "target" with which to verify algorithm performance in a climate context. On decadal time scales there is some evidence from reanalyses and remotely-sensed fluxes alike that tropical ocean-averaged latent heat fluxes have increased 5-10 W/sq m since the early 1990s. However, significant uncertainty surrounds this estimate. Our work addresses the origin of these uncertainties and provides statistics on time series of tropical ocean averages, regional space / time correlation analysis, and separation of contributions by variations in wind and near surface humidity deficit. Comparison to variations in reanalysis data sets is also provided for reference.

Robertson, Franklin R.↗

Updraft Dynamics and Microphysics: On the Added Value of The Cumulus Thermal Reference Frame in Simulations of Aerosol–Deep Convection Interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether, or under which conditions, a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy-permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using the following two different reference frames: static cloudy updraft grid cells versus tracked cumulus thermals. In both frameworks, we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

aerosols↗

Interaction of the terrestrial and atmospheric hydrological cycles in the context of the North American southwest summer monsoon

Work under this grant has used information on precipitation and water vapor fluxes in the area of the Mexican Monsoon to analyze the regional precipitation climatology, to understand the nature of water vapor transport during the monsoon using model and observational data, and to analyze the ability of the TRMM remote sensing algorithm to characterize precipitation. An algorithm for estimating daily surface rain volumes from hourly GOES infrared images was developed and compared to radar data. Estimates were usually within a factor of two, but different linear relations between satellite reflectances and rainfall rate were obtained for each day, storm type and storm development stage. This result suggests that using TRMM sensors to calibrate other satellite IR will need to be a complex process taking into account all three of the above factors. Another study, this one of the space-time variability of the Mexican Monsoon, indicate that TRMM will have a difficult time, over the course of its expected three year lifetime, identifying the diurnal cycle of precipitation over monsoon region. Even when considering monthly rainfalls, projected satellite estimates of August rainfall show a root mean square error of 38 percent. A related examination of spatial variability of mean monthly rainfall using a novel method for removing the effects of elevation from gridded gauge data, show wide variation from a satellite-based rainfall estimates for the same time and space resolution. One issue addressed by our research, relating to the basic character of the monsoon circulation, is the determination of the source region for moisture. The monthly maps produced from our study of monsoon variability show the presence of two rainfall maxima in the analysis normalized to sea level, one in south-central Arizona associated with the Mexican monsoon maximum and one in southeastern New Mexico associated with the Gulf of Mexico. From the point of view of vertically-integrated fluxes and flux divergence of water vapor from ECMWF data, most moisture at upper levels arrives from the Gulf of Mexico, while low level moisture comes from the northern Gulf of California. Composites of ECMWF analyses for wet and dry periods (classified by rain gauge data) show that both regimes show low level moisture arriving from northern and central Gulf of California. Above 700 MB, moisture comes from both source regions and the Sierra Madre Occidental. During wet periods a longer fetch through the moist air mass above western Mexico results in a greater moisture flux into the Sonoran Desert region, while there is less moisture from the Gulf of Mexico both above and below 700 mb. Work on the grant subcontract at the University of Colorado concentrated on the development of a technique useful to TRMM combining visible, infrared and passive microwave data for measuring precipitation. Two established techniques using either visible or infrared data applied over the US Southwest correlated with gauges at the 0.58 to 0.70 level. The application of some established passive microwave techniques were less successful for a variety of reason, including problems in both the gauge and satellite data quality, sampling problems and weaknesses inherent in the algorithms themselves. A more promising solution for accurate rainfall estimation was explored using visible and infrared data to perform a cloud classification, which when combined with information about the background (e.g. Iand/ocean), was used to select the most appropriate microwave algorithm from a suite of possibilities.

Dickinson, Robert E.↗

Updraft Dynamics and Microphysics: On the Added Value of the Cumulus Thermal Reference Frame in Simulations of Aerosol–Deep Convection Interactions

One fundamental question about atmospheric moist convection processes that remains debated is whether or under what conditions a relevant variability in background aerosol concentrations may have a significant dynamical impact on convective clouds and their associated precipitation. Furthermore, current climate models must parameterize both the microphysical and the cumulus convection processes, but this is usually implemented separately, whereas in nature there is a strong coupling between them. As a first step to improve our understanding of these two problems, we investigate how aerosol concentrations modify key properties of updrafts in eight large-eddy permitting regional simulations of a case study of scattered convection over Houston, Texas, in which convection is explicitly simulated and microphysical processes are parameterized. Dynamical and liquid-phase microphysical responses are investigated using two different reference frames: static cloudy-updraft grid cells versus tracked cumulus thermals. In both frameworks we observe the expected microphysical responses to higher aerosol concentrations, such as higher cloud number concentrations and lower rain number concentrations. In terms of the dynamical responses, both frameworks indicate weak impacts of varying aerosol concentrations relative to the noise between simulations over the observationally derived range of aerosol variability for this case study. On the other hand, results suggest that thermals are more selective than cloudy-updraft grid cells in terms of sampling the most active convective air masses. For instance, vertical velocity from thermals is significantly higher at upper levels than when sampled from cloudy-updraft grid points, and several microphysical variables have higher average values in the cumulus thermal framework than in the cloudy-updraft framework. In addition, the thermal analysis is seen to add rich quantitative information about the rates and covariability of microphysical processes spatially and throughout tracked thermal lifecycles, which can serve as a stronger foundation for improving subgrid-scale parameterizations.

Daniel Hernandez-Deckers↗

Development of an Innovative Algorithm for Aerodynamics-Structure Interaction Using Lattice Boltzmann Method

The lattice Boltzmann equation (LBE) is a kinetic formulation which offers an alternative computational method capable of solving fluid dynamics for various systems. Major advantages of the method are owing to the fact that the solution for the particle distribution functions is explicit, easy to implement, and the algorithm is natural to parallelize. In this final report, we summarize the works accomplished in the past three years. Since most works have been published, the technical details can be found in the literature. Brief summary will be provided in this report. In this project, a second-order accurate treatment of boundary condition in the LBE method is developed for a curved boundary and tested successfully in various 2-D and 3-D configurations. To evaluate the aerodynamic force on a body in the context of LBE method, several force evaluation schemes have been investigated. A simple momentum exchange method is shown to give reliable and accurate values for the force on a body in both 2-D and 3-D cases. Various 3-D LBE models have been assessed in terms of efficiency, accuracy, and robustness. In general, accurate 3-D results can be obtained using LBE methods. The 3-D 19-bit model is found to be the best one among the 15-bit, 19-bit, and 27-bit LBE models. To achieve desired grid resolution and to accommodate the far field boundary conditions in aerodynamics computations, a multi-block LBE method is developed by dividing the flow field into various blocks each having constant lattice spacing. Substantial contribution to the LBE method is also made through the development of a new, generalized lattice Boltzmann equation constructed in the moment space in order to improve the computational stability, detailed theoretical analysis on the stability, dispersion, and dissipation characteristics of the LBE method, and computational studies of high Reynolds number flows with singular gradients. Finally, a finite difference-based lattice Boltzmann method is developed for inviscid compressible flows.

Mei, Ren-Wei↗

Community-Based Services that Facilitate Interoperability and Intercomparison of Precipitation Datasets from Multiple Sources

Over the past 12 years, large volumes of precipitation data have been generated from space-based observatories (e.g., TRMM), merging of data products (e.g., gridded 3B42), models (e.g., GMAO), climatologies (e.g., Chang SSM/I derived rain indices), field campaigns, and ground-based measuring stations. The science research, applications, and education communities have greatly benefited from the unrestricted availability of these data from the Goddard Earth Sciences Data and Information Services Center (GES DISC) and, in particular, the services tailored toward precipitation data access and usability. In addition, tools and services that are responsive to the expressed evolving needs of the precipitation data user communities have been developed at the Precipitation Data and Information Services Center (PDISC) (http://disc.gsfc.nasa.gov/precipitation or google NASA PDISC), located at the GES DISC, to provide users with quick data exploration and access capabilities. In recent years, data management and access services have become increasingly sophisticated, such that they now afford researchers, particularly those interested in multi-data set science analysis and/or data validation, the ability to homogenize data sets, in order to apply multi-variant, comparison, and evaluation functions. Included in these services is the ability to capture data quality and data provenance. These interoperability services can be directly applied to future data sets, such as those from the Global Precipitation Measurement (GPM) mission. This presentation describes the data sets and services at the PDISC that are currently used by precipitation science and applications researchers, and which will be enhanced in preparation for GPM and associated multi-sensor data research. Specifically, the GES-DISC Interactive Online Visualization ANd aNalysis Infrastructure (Giovanni) will be illustrated. Giovanni enables scientific exploration of Earth science data without researchers having to perform the complicated data access and match-up processes. In addition, PDISC tool and service capabilities being adapted for GPM data will be described, including the Google-like Mirador data search and access engine; semantic technology to help manage large amounts of multi-sensor data and their relationships; data access through various Web services (e.g., OPeNDAP, GDS, WMS, WCS); conversion to various formats (e.g., netCDF, HDF, KML (for Google Earth)); visualization and analysis of Level 2 data profiles and maps; parameter and spatial subsetting; time and temporal aggregation; regridding; data version control and provenance; continuous archive verification; and expertise in data-related standards and interoperability. The goal of providing these services is to further the progress towards a common framework by which data analysis/validation can be more easily accomplished.

Liu, Zhong↗

An Operational Configuration of the ARPS Data Analysis System to Initialize WRF in the NM'S Environmental Modeling System

The Weather Research and Forecasting (WRF) model is the next generation community mesoscale model designed to enhance collaboration between the research and operational sectors. The NM'S as a whole has begun a transition toward WRF as the mesoscale model of choice to use as a tool in making local forecasts. Currently, both the National Weather Service in Melbourne, FL (NWS MLB) and the Spaceflight Meteorology Group (SMG) are running the Advanced Regional Prediction System (AIRPS) Data Analysis System (ADAS) every 15 minutes over the Florida peninsula to produce high-resolution diagnostics supporting their daily operations. In addition, the NWS MLB and SMG have used ADAS to provide initial conditions for short-range forecasts from the ARPS numerical weather prediction (NWP) model. Both NM'S MLB and SMG have derived great benefit from the maturity of ADAS, and would like to use ADAS for providing initial conditions to WRF. In order to assist in this WRF transition effort, the Applied Meteorology Unit (AMU) was tasked to configure and implement an operational version of WRF that uses output from ADAS for the model initial conditions. Both agencies asked the AMU to develop a framework that allows the ADAS initial conditions to be incorporated into the WRF Environmental Modeling System (EMS) software. Developed by the NM'S Science Operations Officer (S00) Science and Training Resource Center (STRC), the EMS is a complete, full physics, NWP package that incorporates dynamical cores from both the National Center for Atmospheric Research's Advanced Research WRF (ARW) and the National Centers for Environmental Prediction's Non-Hydrostatic Mesoscale Model (NMM) into a single end-to-end forecasting system. The EMS performs nearly all pre- and postprocessing and can be run automatically to obtain external grid data for WRF boundary conditions, run the model, and convert the data into a format that can be readily viewed within the Advanced Weather Interactive Processing System. The EMS has also incorporated the WRF Standard Initialization (SI) graphical user interface (GUT), which allows the user to set up the domain, dynamical core, resolution, etc., with ease. In addition to the SI GUT, the EMS contains a number of configuration files with extensive documentation to help the user select the appropriate input parameters for model physics schemes, integration timesteps, etc. Therefore, because of its streamlined capability, it is quite advantageous to configure ADAS to provide initial condition data to the EMS software. One of the biggest potential benefits of configuring ADAS for ingest into the EMS is that the analyses could be used to initialize either the ARW or NMM. Currently, the ARPS/ADAS software has a conversion routine only for the ARW dynamical core. However, since the NIvIM runs about 2.5 times faster than the ARW, it is quite advantageous to be able to run an ADAS/NMM configuration operationally due to the increased efficiency.

Case, Jonathan↗

NASA Center for Climate Simulation (NCCS) Presentation

The NASA Center for Climate Simulation (NCCS) offers integrated supercomputing, visualization, and data interaction technologies to enhance NASA's weather and climate prediction capabilities. It serves hundreds of users at NASA Goddard Space Flight Center, as well as other NASA centers, laboratories, and universities across the US. Over the past year, NCCS has continued expanding its data-centric computing environment to meet the increasingly data-intensive challenges of climate science. We doubled our Discover supercomputer's peak performance to more than 800 teraflops by adding 7,680 Intel Xeon Sandy Bridge processor-cores and most recently 240 Intel Xeon Phi Many Integrated Core (MIG) co-processors. A supercomputing-class analysis system named Dali gives users rapid access to their data on Discover and high-performance software including the Ultra-scale Visualization Climate Data Analysis Tools (UV-CDAT), with interfaces from user desktops and a 17- by 6-foot visualization wall. NCCS also is exploring highly efficient climate data services and management with a new MapReduce/Hadoop cluster while augmenting its data distribution to the science community. Using NCCS resources, NASA completed its modeling contributions to the Intergovernmental Panel on Climate Change (IPCG) Fifth Assessment Report this summer as part of the ongoing Coupled Modellntercomparison Project Phase 5 (CMIP5). Ensembles of simulations run on Discover reached back to the year 1000 to test model accuracy and projected climate change through the year 2300 based on four different scenarios of greenhouse gases, aerosols, and land use. The data resulting from several thousand IPCC/CMIP5 simulations, as well as a variety of other simulation, reanalysis, and observationdatasets, are available to scientists and decision makers through an enhanced NCCS Earth System Grid Federation Gateway. Worldwide downloads have totaled over 110 terabytes of data.

Webster, William P.↗

Vortex Generators in a Streamline-Traced, External-Compression Supersonic Inlet

Vortex generators within a streamline-traced, external-compression supersonic inlet for Mach 1.66 were investigated to determine their ability to increase total pressure recovery and reduce total pressure distortion. The vortex generators studied were rectangular vanes arranged in counter-rotating and co-rotating arrays. The vane geometric factors of interest included height, length, spacing, angle-of-incidence, and positions upstream and downstream of the inlet terminal shock. The flow through the inlet was simulated numerically through the solution of the steady-state, Reynolds-averaged Navier-Stokes equations on multi-block, structured grids using the Wind-US flow solver. The vanes were simulated using a vortex generator model. The inlet performance was characterized by the inlet total pressure recovery and the radial and circumferential total pressure distortion indices at the engine face. Design of experiments and statistical analysis methods were applied to quantify the effect of the geometric factors of the vanes and search for optimal vane arrays. Co-rotating vane arrays with negative angles-of-incidence positioned on the supersonic diffuser were effective in sweeping low-momentum flow from the top toward the sides of the subsonic diffuser. This distributed the low-momentum flow more evenly about the circumference of the subsonic diffuser and reduced distortion. Co-rotating vane arrays with negative angles-of-incidence or counter-rotating vane arrays positioned downstream of the terminal shock were effective in mixing higher-momentum flow with lower-momentum flow to increase recovery and decrease distortion. A strategy of combining a co-rotating vane array on the supersonic diffuser with a counter-rotating vane array on the subsonic diffuser was effective in increasing recovery and reducing distortion.

computational fluid dynamics↗

CFD Analysis of Modular Thrusters Performance

The effective performance of modular thrusters in an aerospike configuration is difficult to determine. Standard analytical tools are applicable to conventional nozzle shapes, but are limited when applied to an aerospike nozzle (An aerospike nozzle is an altitude compensating external nozzle). Three baseline nozzle shapes are derived using standard analytical procedures. The baseline nozzle sizes are restricted to fill a volume envelope. The three shapes are an axi-symmetric round nozzle, a two dimensional planar square exit nozzle, and a super elliptic round to nearly square nozzle. The integrated (thruster/aerospike) performance of the three nozzles is determined through the use of three dimensional viscous computational fluid dynamic (CFD) calculation where complex features of the flow field can be accurately captured. The resulting installed performance is then used to evaluate the efficiency of these nozzle shapes for aerospike applications. The determination of effective performance of a thruster nozzle integrated into an aerospike nozzle require the solution of the three dimensional turbulent Navier-Stokes equations. The model used in this study consisted of two zones; one of the upstream thruster cowl surface so freestream conditions can be accurately predicted, and two, the aerospike surface beginning with the thruster outflow and extending to the end of the aerospike surface. The numerical grid consisted of over 120,000 nodes and used symmetry on the thruster centerline and edge. A two species non-reacting chemistry model was used to capture the variation of fluid properties between the hot plume base and freestream air. From the results of the three baseline nozzle aerospike calculations, the effictive performance of the nozzle was determined. The flow fields of these calculations do show some variation between the cases. Recirculation zones on the cowl surface is predicted for the two dimensional planar nozzle and a smaller one for the super elliptic nozzle. The recirculation is caused by the strong pressure gradient between the plume and freestream flows. The axi-symmetric nozzle results indicate recirculation zones on the thruster face. These recirculation zones smooth the pressure gradient between the plume and freestream flow limiting the formation of recirculation on the cowl surface. Thruster to thruster interaction is evident for the axi-symmetric and supper elliptic calculation while the two dimensional planar nozzle did not have any lateral expansion in the nozzle, so thruster to thruster interaction is limited. The integrated performance results, at the altitude choosen, show very little variation between the three thruster shapes. This result allows for nozzle shape determination based on additional considerations (thermal, structural, weight) besides performance.

Ungewitter, Ronald J.↗

A General Interface Method for Aeroelastic Analysis of Aircraft

The aeroelastic analysis of an aircraft requires an accurate and efficient procedure to couple aerodynamics and structures. The procedure needs an interface method to bridge the gap between the aerodynamic and structural models in order to transform loads and displacements. Such an interface method is described in this report. This interface method transforms loads computed by any aerodynamic code to a structural finite element (FE) model and converts the displacements from the FE model to the aerodynamic model. The approach is based on FE technology in which virtual work is employed to transform the aerodynamic pressures into FE nodal forces. The displacements at the FE nodes are then converted back to aerodynamic grid points on the aircraft surface through the reciprocal theorem in structural engineering. The method allows both high and crude fidelities of both models and does not require an intermediate modeling. In addition, the method performs the conversion of loads and displacements directly between individual aerodynamic grid point and its corresponding structural finite element and, hence, is very efficient for large aircraft models. This report also describes the application of this aero-structure interface method to a simple wing and an MD-90 wing. The results show that the aeroelastic effect is very important. For the simple wing, both linear and nonlinear approaches are used. In the linear approach, the deformation of the structural model is considered small, and the loads from the deformed aerodynamic model are applied to the original geometry of the structure. In the nonlinear approach, the geometry of the structure and its stiffness matrix are updated in every iteration and the increments of loads from the previous iteration are applied to the new structural geometry in order to compute the displacement increments. Additional studies to apply the aero-structure interaction procedure to more complicated geometry will be conducted in the second phase of the present contract.

Tzong, T.↗

Observations on the long-period variability of the Gulf Stream downstream of Cape Hatteras

To examine the long-period variability of the Gulf Stream, sea level residuals relative to a 2-year mean sea level in the Gulf Stream downstream of Cape Hatteras (between 75 deg W and 60 deg W longitude) are used. Residuals, as derived from Geosat altimetry between November 1986 and December 1988, were gridded in space and time at a temporal resolution of 10 days and spatial resolution of 1/4 deg. Complex empirical orthogonal function (CEOF) analysis was applied to the data set to extract the spatially correlated signal with the original data subsampled to 1/2 deg. In addition to determining the space-time scales and propagation characterisitics of the different modes, wavenumber-frequency spectral techniques were used to separate the variability into propagating and stationary components. The CEOF technique applied to the data set indicated that the first four CEOF modes accounted for 60% of the variability and were found to be above the noise leve 99% of the time. CEOF 1 was associated with westward propagation at 5 km/d at a wavelength of 2000 km and eastward propagation at 1-2 km/d centered at a 500-km wavelength. This first CEOF is in good agreement with thin-jet equivalent barotropic models which predict westward propagation for wavelengths greater than 1130 km. A deflection of the wavelike pattern at 65 deg W also indicates a possible topographic effect. A simple scaling of the effect of topography indicates that for length scales longer than the internal Rossby radius of deformation, the topographic term is at least of the same order of magnitude as the beta effect. The second CEOF was more broadbanded in wavenumber space, with eastward propagation occurring in a wavenumber-frequency band between 300 and 1400 km and 0.5 and 2.0 cycles/yr. The third CEOF is similar in structure to the first, but with less energy. CEOF 4 was clearly identifiable with higher frequencies than the first three with westward propagation at 4 km/d. The spatial location of this mode along with the westward propagation indicates possible influences from eddy-stream interactions. Thus topography, Rossby wave dynamics and eddy-stream interactions all appear to have a significant role in determining the space-time scales and propagation properties of the long-period response of sea level in the Gulf Stream.

Vazquez, Jorge↗