Search NASA⌕ Search

SEARCH · Search NASA

Results for “Overlapping grid”

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 145 records · Page 8

Numerical simulation of axisymmetric turbulent flow in combustors and diffusors

A method for predicting turbulent flow in combustors and diffusers is developed. The Navier-Stokes equations, incorporating a turbulence kappa-epsilon model equation, were solved in a nonorthogonal curvilinear coordinate system. The solution applied the finite volume method to discretize the differential equations and utilized the SIMPLE algorithm iteratively to solve the differenced equations. A zonal grid method, wherein the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. However, treatment of the zonal boundaries required special handling. Boundary overlap and interpolating techniques were used and an adjustment of the flow variables was required to assure conservation of mass, momentum and energy fluxes. The numerical accuracy was assessed using different finite differencing methods, i.e., hybrid, quadratic upwind and skew upwind, to represent the convection terms. Flows in different geometries of combustors and diffusers were simulated and results compared with experimental data and good agreement was obtained.

Yung, Chain Nan↗

The correlated-k method for radiation calculations in nonhomogeneous atmospheres

The accuracy of the correlated-k method, which is a technique for radiation calculations with spectrally averaged data in nonhomogeneous atmospheres, is investigated. Comparisons are made for scattering and absorbing atmospheres containing CO2, H2O, and O3, and it is concluded that: (1) the errors in correlated-k are generally of order of magnitude 1 percent, (2) much larger errors occur only when a radiative quantity is very much smaller than its average value, (3) errors do not depend systematically on the gas molecule, the distributions of gases and aerosols, or on the aerosol optical properties, and (4) errors do not systematically increase with the order of differencing. It is shown that the multiplication property for transmission by overlapping bands can be incorporated into correlated-k, that temperature effects can be interpolated on a coarse grid, and that 10 quadrature points are often sufficient to average over complex spectral intervals containing thousands of lines.

Goody, Richard↗

Evaluating the design of an earth radiation budget instrument with system simulations. Part 2: Minimization of instantaneous sampling errors for CERES-I

Much of the new record of broadband earth radiation budget satellite measurements to be obtained during the late 1990s and early twenty-first century will come from the dual-radiometer Clouds and Earth's Radiant Energy System Instrument (CERES-I) flown aboard sun-synchronous polar orbiters. Simulation studies conducted in this work for an early afternoon satellite orbit indicate that spatial root-mean-square (rms) sampling errors of instantaneous CERES-I shortwave flux estimates will range from about 8.5 to 14.0 W/m on a 2.5 deg latitude and longitude grid resolution. Rms errors in longwave flux estimates are only about 20% as large and range from 1.5 to 3.5 W/sq m. These results are based on an optimal cross-track scanner design that includes 50% footprint overlap to eliminate gaps in the top-of-the-atmosphere coverage, and a 'smallest' footprint size to increase the ratio in the number of observations lying within to the number of observations lying on grid area boundaries. Total instantaneous measurement error also depends on the variability of anisotropic reflectance and emission patterns and on retrieval methods used to generate target area fluxes. Three retrieval procedures from both CERES-I scanners (cross-track and rotating azimuth plane) are used. (1) The baseline Earth Radiaton Budget Experiment (ERBE) procedure, which assumes that errors due to the use of mean angular dependence models (ADMs) in the radiance-to-flux inversion process nearly cancel when averaged over grid areas. (2) To estimate N, instantaneous ADMs are estimated from the multiangular, collocated observations of the two scanners. These observed models replace the mean models in computation of satellite flux estimates. (3) The scene flux approach, conducts separate target-area retrievals for each ERBE scene category and combines their results using area weighting by scene type. The ERBE retrieval performs best when the simulated radiance field departs from the ERBE mean models by less than 10%. For larger perturbations, both the scene flux and collocation methods produce less error than the ERBE retrieval. The scene flux technique is preferable, however, because it involves fewer restrictive assumptions.

Stowe, Larry↗

Development and Assessment of the SMAP Enhanced Passive Soil Moisture Product

Launched in January 2015, the National Aeronautics and Space Administration (NASA) Soil Moisture Active Passive (SMAP) observatory was designed to provide frequent global mapping of high-resolution soil moisture and freeze-thaw state every two to three days using a radar and a radiometer operating at L-band frequencies. Despite a hardware mishap that rendered the radar inoperable shortly after launch, the radiometer continues to operate nominally, returning more than two years of science data that have helped to improve existing hydrological applications and foster new ones.Beginning in late 2016 the SMAP project launched a suite of new data products with the objective of recovering some high-resolution observation capability loss resulting from the radar malfunction. Among these new data products are the SMAP Enhanced Passive Soil Moisture Product that was released in December 2016, followed by the SMAPSentinel-1 Active-Passive Soil Moisture Product in April 2017.This article covers the development and assessment of the SMAP Level 2 Enhanced Passive Soil Moisture Product (L2_SM_P_E). The product distinguishes itself from the current SMAP Level 2 Passive Soil Moisture Product (L2_SM_P) in that the soil moisture retrieval is posted on a 9 km grid instead of a 36 km grid. This is made possible by first applying the Backus-Gilbert optimal interpolation technique to the antenna temperature (TA) data in the original SMAP Level 1B Brightness Temperature Product to take advantage of the overlapped radiometer footprints on orbit. The resulting interpolated TA data then go through various correctioncalibration procedures to become the SMAP Level 1C Enhanced Brightness Temperature Product (L1C_TB_E). The L1C_TB_E product, posted on a 9 km grid, is then used as the primary input to the current operational SMAP baseline soil moisture retrieval algorithm to produce L2_SM_P_E as the final output. Images of the new product reveal enhanced visual features that are not apparent in the standard product. Based on in situ data from core validation sites and sparse networks representing different seasons and biomes all over the world, comparisons between L2_SM_P_E and in situ data were performed for the duration of April 1, 2015 October 30, 2016. It was found that the performance of the enhanced 9 km L2_SM_P_E is equivalent to that of the standard 36 km L2_SM_P, attaining a retrieval uncertainty below 0.040 m(exp 3)/m(exp 3) unbiased root-mean-square error (ubRMSE) and a correlation coefficient above 0.800. This assessment also affirmed that the Single Channel Algorithm using the V-polarized TB channel (SCA-V) delivered the best retrieval performance among the various algorithms implemented for L2_SM_P_E, a result similar to a previous assessment for L2_SM_P.

Enhanced↗

Kinematic domain decomposition to simulate flows past moving objects

A new method is developed to solve the unsteady Navier-Stokes equations on a composite grid which consists of subdomain grids moving with respect to each other. These subdomains are structured grids with different topologies. However, their interfaces are unstructured in the sense that they may or may not arbitrarily overlap with each other. This method captures the boundary-motion-induced flow component. The method is demonstrated for a transonic flow past an airfoil which experiences a combined motion of pitching and plunging.

Baysal, Oktay↗

A Daily Surface Melt Product From Passive Microwave Swath Observations for Arctic Sea Ice Studies

Most passive microwave sea ice data products are produced with daily gridded brightness temperatures. However, a single polar-orbiting passive microwave instrument makes many measurements per day where swaths overlap in the polar regions. These level 1 swath brightness temperature observations can be used to derive information on sub-daily sea ice processes. One such process is the rapid transition of snow and ice between freezing and melting due to diurnal freeze-thaw cycling or brief warming events in the pre-melt season. Here, a new method utilizing the sub-daily variability in swath-to-swath observations is used to produce daily, binary maps of the melting status on Arctic Sea ice using Level 1C brightness temperatures from the Special Sensor Microwave Imager and Sounder (SSM/I-SSMIS) data record. The resulting data set will enable investigations of the frequency, inter-annual variability, and spatial extent of brief surface melting events occurring in the winter and shoulder seasons outside of the traditional summer melt season.

Angela C. Bliss↗

Investigation of turbulent flow in highly curved ducts with application to turbomachinery components

Numerical predictions were performed using a semi-elliptic calculation procedure for the case of turbulent flow in passage through a 90 deg bend of square cross section. Two versions of the isotropic turbulent viscosity two equation k-epsilon model were used. The first, the wall function model (WFM), employs the logarithmic law-of-the-wall relation and the notion of equilibrium flow to set all the necessary boundary conditions at the first grid point adjacent to a solid wall. The second model, the Van Driest model (VDM), employs Prandtl's original mixing length formulation, in conjunction with Van Driest's semi-empirical relation for the mixing length, to calculate the turbulent viscosity in the near wall regions of the flow. In this case, boundary conditions for k and epsilon, required to calculate these quantities in the core of the flow, are obtained by matching the mixing length and Reynolds number model formulations in an overlapping region of the flow near the walls. In both cases the results obtained show an improvement over earlier calculations using an elliptic numerical procedure. This is attributed to the finer grids possible in the present work. Of the two models, the VDM formulation shows better overall conformity with the mean flow measurements. Neither model reproduces well the details of the stress distribution as a result of the implied isotropic turbulent viscosity.

Humphrey, J. A. C.↗

Western Rainier Seismic Zone Airborne Laser Swath Mapping

Airborne laser swath mapping (ALSM) of the Puget Lowland conducted by TerraPoint LLC for the Purget Sound Lidar Concortium (PSLC), has been successful in revealing Holocene fault scarps and lendsliders hidden beneath the dense, temperate rain forest cover and in quantifying shoreline terrace uplift. Expanding the PSLC efforts, NASA-USGS collaboration is now focusing on topographic mapping of seismogenic zones adjacent to volcanois in the western Cascades range in order to assess the presence of active faulting and tectonic deformation, better define the extend of lahars and understand their flow processes, and characterize landslide occurrence. Mapping of the western Rainier zone (WRZ) was conducted by TerraPoint in late 2002, after leaf fall and before snow accumulation. The WRZ is a NNW-trending, approx. 30 km-long zone of seismicity west of Mount Rainier National Park. The Puget Lowland ALSM methods were modified to accommodate challenges posed by the steep, high relief terrian. The laser data, acquired with a density of approx. 2 pulses /sq m, was filtered to identify returns from the ground from which a bare Earth digital elevation model (DEM) was produced with a grid size of 1.8 m. The RMS elevation accuracy of the DEM in flat, unvegetated areas is approx. 10cm based on consistency between overlapping flight swaths and comparisons to ground control points. The resulting DEM substantially improves upon Shuttle Radar Topography Mission and USGS photogrammetric mapping. For example, the DEM defines the size and spatial distribution of flood erratics left by the Electron lahar and of megaclasts within the Round Pass lahar, important for characterizing the lahar hydraulics. A previously unknown lateral levee on the Round Pass lahar is also revealed. In addition, to illustrating geomorfic feature within the WRZ, future plans for laser mapping of the Saint Helens and Darrington seismic zones will be described.

Harding, David J.↗

A Novel Method for Estimating Shortwave Direct Radiative Effect of Above-Cloud Aerosols Using CALIOP and MODIS Data

This paper describes an efficient and unique method for computing the shortwave direct radiative effect (DRE) of aerosol residing above low-level liquid-phase clouds using CALIOP and MODIS data. It accounts for the overlapping of aerosol and cloud rigorously by utilizing the joint histogram of cloud optical depth and cloud top pressure. Effects of sub-grid scale cloud and aerosol variations on DRE are accounted for. It is computationally efficient through using grid-level cloud and aerosol statistics, instead of pixel-level products, and a pre-computed look-up table in radiative transfer calculations. We verified that for smoke over the southeast Atlantic Ocean the method yields a seasonal mean instantaneous shortwave DRE that generally agrees with more rigorous pixel-level computation within 4. We have also computed the annual mean instantaneous shortwave DRE of light-absorbing aerosols (i.e., smoke and polluted dust) over global ocean based on 4 yr of CALIOP and MODIS data. We found that the variability of the annual mean shortwave DRE of above-cloud light-absorbing aerosol is mainly driven by the optical depth of the underlying clouds.

radiative forcing↗

A Novel Method for Estimating Shortwave Direct Radiative Effect of Above-cloud Aerosols over Ocean Using CALIOP and MODIS Data

This paper describes an efficient and unique method for computing the shortwave direct radiative effect (DRE) of aerosol residing above low-level liquid-phase clouds using CALIOP and MODIS data. It accounts for the overlapping of aerosol and cloud rigorously by utilizing the joint histogram of cloud optical depth and cloud top pressure. Effects of sub-grid scale cloud and aerosol variations on DRE are accounted for. It is computationally efficient through using grid-level cloud and aerosol statistics, instead of pixel-level products, and a pre-computed look-up table in radiative transfer calculations. We verified that for smoke over the southeast Atlantic Ocean the method yields a seasonal mean instantaneous shortwave DRE that generally agrees with more rigorous pixel-level computation within 4%. We have also computed the annual mean instantaneous shortwave DRE of light-absorbing aerosols (i.e., smoke and polluted dust) over global ocean based on 4 yr of CALIOP and MODIS data. We found that the variability of the annual mean shortwave DRE of above-cloud light-absorbing aerosol is mainly driven by the optical depth of the underlying clouds.

CALIOP and MODIS↗

Impact of Satellite Sea Surface Salinity Observations on ENSO Predictions from the GMAO Seasonal Forecast System

El Nino/Southern Oscillation (ENSO) has far reaching global climatic impacts and so extending useful ENSO forecasts would be of great benefit for society. However, one key variable that has yet to be fully exploited within coupled forecast systems is accurate estimation of near-surface ocean density. Satellite sea surface salinity (SSS), combined with temperature, help to identify ocean density changes and associated mixing near the ocean surface. We assess the impact of satellite SSS observations for improving near-surface dynamics within ocean analyses and how these impact dynamical ENSO forecasts using the NASA GMAO (Global Modeling and Assimilation Office) Sub-seasonal to Seasonal (S2S_v2.1) coupled forecast system (Molod et al. 2018 - i.e. NASA's contribution to the NMME (North American Multi-Model Ensemble) project). For all initialization experiments, all available along-track absolute dynamic topography and in situ observations are assimilated using the LETKF (Local Ensemble Transform Kalman Filter) scheme similar to Penny et al., 2013. A separate reanalysis additionally assimilates Aquarius V5 (September 2011 to June 2015) and SMAP (Soil Moisture Active Passive satellite) V4.1 (March 2015 to present) along-track data.We highlight the impact of satellite SSS on ocean reanalyses by comparing validation statistics of experiments that assimilate SSS versus our current prediction system that withholds SSS. We find that near-surface validation versus observed statistics for salinity are slightly degraded when assimilating SSS. This is an expected result due to known biases between SSS (measured by satellite at approximately 1-centimeter depth) and in situ measurements (typically measured by Argo floats at 3-meters depth). On the other hand, a very encouraging result is that both temperature, absolute dynamic topography, and mixed layer statistics are improved with SSS assimilation. Previous work has shown that correcting near-surface density structure via gridded SSS assimilation can improve coupled forecasts. Here we present results of coupled forecasts that are initialized from GMAO S2S spring reanalyses that assimilate/withhold along-track (L2) SSS. In particular, we contrast forecasts of the big 2015 El Nino, the 2017 La Nina and the 2018 weak El Nino. For each of these ENSO scenarios, assimilation of satellite SSS improves the forecast validation. Improved SSS and density upgrade the mixed layer depth leading to more accurate coupled air/sea interaction. From March to June 2015, the availability of two overlapping satellite SSS instruments, Aquarius and SMAP, allows a unique opportunity to compare and contrast forecasts initialized with the benefit of these two satellite SSS observation types. We assess the impact of gridded satellite sea surface salinity observations on dynamical ENSO forecasts for the big 2015 El Nino.

Hackert, E.↗

Airplane Mesh Development with Grid Density Studies

Automatic Grid Generation Wish List Geometry handling, including CAD clean up and mesh generation, remains a major bottleneck in the application of CFD methods. There is a pressing need for greater automation in several aspects of the geometry preparation in order to reduce set up time and eliminate user intervention as much as possible. Starting from the CAD representation of a configuration, there may be holes or overlapping surfaces which require an intensive effort to establish cleanly abutting surface patches, and collections of many patches may need to be combined for more efficient use of the geometrical representation. Obtaining an accurate and suitable body conforming grid with an adequate distribution of points throughout the flow-field, for the flow conditions of interest, is often the most time consuming task for complex CFD applications. There is a need for a clean unambiguous definition of the CAD geometry. Ideally this would be carried out automatically by smart CAD clean up software. One could also define a standard piece-wise smooth surface representation suitable for use by computational methods and then create software to translate between the various CAD descriptions and the standard representation. Surface meshing remains a time consuming, user intensive procedure. There is a need for automated surface meshing, requiring only minimal user intervention to define the overall density of mesh points. The surface mesher should produce well shaped elements (triangles or quadrilaterals) whose size is determined initially according to the surface curvature with a minimum size for flat pieces, and later refined by the user in other regions if necessary. Present techniques for volume meshing all require some degree of user intervention. There is a need for fully automated and reliable volume mesh generation. In addition, it should be possible to create both surface and volume meshes that meet guaranteed measures of mesh quality (e.g. minimum and maximum angle, stretching ratios, etc.).

Cliff, Susan E.↗

Development and preliminary tests of a new long-wave radiation parameterization

An efficient broad band longwave radiation code for CO2 and H2O and for O3 was developed. There are two bands each in the CO2 and H2O absorption regions, one for the band center and one for the band wings. One band covers O3 absorption and the overlapping H2O continuum. Overlap is also considered in the CO2 region, and there is H2O continuum absorption where applicable. Clouds are considered nonreflecting in the longwave. Therefore partial cover or partial transmission can be allowed for, by considering a cloud fraction at each atmospheric level. A special subroutine was written to allow for maximum or random overlap of clouds that may be used in the future. All algorithms were written with vectorization in mind with identical operations made for all horizontal grid points in a latitude circle. Where possible, operations are carried out covering the vertical grid points as well, yielding long vectors for efficient computations.

HARSHVARDHAN↗

Performance and Application of Parallel OVERFLOW Codes on Distributed and Shared Memory Platforms

The presentation discusses recent studies on the performance of the two parallel versions of the aerodynamics CFD code, OVERFLOW_MPI and _MLP. Developed at NASA Ames, the serial version, OVERFLOW, is a multidimensional Navier-Stokes flow solver based on overset (Chimera) grid technology. The code has recently been parallelized in two ways. One is based on the explicit message-passing interface (MPI) across processors and uses the _MPI communication package. This approach is primarily suited for distributed memory systems and workstation clusters. The second, termed the multi-level parallel (MLP) method, is simple and uses shared memory for all communications. The _MLP code is suitable on distributed-shared memory systems. For both methods, the message passing takes place across the processors or processes at the advancement of each time step. This procedure is, in effect, the Chimera boundary conditions update, which is done in an explicit "Jacobi" style. In contrast, the update in the serial code is done in more of the "Gauss-Sidel" fashion. The programming efforts for the _MPI code is more complicated than for the _MLP code; the former requires modification of the outer and some inner shells of the serial code, whereas the latter focuses only on the outer shell of the code. The _MPI version offers a great deal of flexibility in distributing grid zones across a specified number of processors in order to achieve load balancing. The approach is capable of partitioning zones across multiple processors or sending each zone and/or cluster of several zones into a single processor. The message passing across the processors consists of Chimera boundary and/or an overlap of "halo" boundary points for each partitioned zone. The MLP version is a new coarse-grain parallel concept at the zonal and intra-zonal levels. A grouping strategy is used to distribute zones into several groups forming sub-processes which will run in parallel. The total volume of grid points in each group are approximately balanced. A proper number of threads are initially allocated to each group, and in subsequent iterations during the run-time, the number of threads are adjusted to achieve load balancing across the processes. Each process exploits the multitasking directives already established in Overflow.

Djomehri, M. Jahed↗

Precision Lower Stratospheric Temperature Monitoring with the MSU: Technique, Validation, and Results, 1979-1991

Microwave Sounding Unit channel 4 data from the TIROS-N (Television and Infrared Operational Satellite)-N series of NOAA (National Oceanic and Atmospheric Administration) satellites are intercalibrated to provide a continuous global record of deep-layer averaged lower stratospheric temperatures during 1979-1991. A 13-year record of temperature anomalies is time averaged into pentads and months on a 2.5 deg. grid. The monthly gridpoint anomalies are validated with ten years of radiosonde data during 1979-88. The calibration stability of each satellite's measurements is evaluated during satellite overlap periods, the longest of which reveal no measurable instrumental drift at the level of 0.01 C yr (exp -1). Intercomparisons between NOAA-6 and NOAA- 7 anomalies indicate monthly gridpoint precision of 0.05 C in the tropics to around 0.10 C in the extratropics, and signal-to-noise ratios precision and stability statistics are much better than have been previously reported by other investigators for MSU channel 4. Pentad precision is about 0.10 C in the tropics to around 0.25 C at high latitudes and signal-to- noise ratios generally over 250 in the tropics and high latitudes, but 100-200 in the middle latitudes. Radiosonde comparisons to the monthly gridpoint anomalies have correlations ranging from 0.90 in the tropics (where the interannual variability is smallest) to as high as 0.99 at high-latitude stations. The corresponding standard error of estimate is generally around 0.3 C. A significant difference in decadal trends is found between the satellite and radiosonde systems, with a step change of 0.21 C (sondes cooler) compared to the satellite measurements. Investigations of the possible sources of the discrepancy lead us to suspect that the gradual transition from on-site calibration of sondes with thermometers to factory calibration of sondes around 1982 might have caused a change in the calibration, although this conclusion must be viewed as tentative. The largest globally averaged temperature variations during 1979-91 occur after the El Chichon (1982) and Pinatubo (1991) volcanic eruptions. These warm events are superimposed upon a net downward trend in temperatures during the period. This cooling trend has more of a step function than linear character, with the step occurring during the El Chichon warm event. It is strongest in polar regions and the Northern Hemisphere middle latitudes. These characteristics are qualitatively consistent with radiative adjustments expected to occur with observed ozone depletions.

Spencer, Roy W.↗

Planning for the utilization of the PCDS in studying the interaction of clouds (ISCCP-C data) and the Earth radiation budget (ERBE data)

The Pilot Climate Data System (PCDS) affords an opportunity to analyze data from different but highly complementary data sets. Two of these highly complementary data sets supported by the PCDS are the International Satellite Cloud Climatology Project (ISCCP) and the Earth Radiation Budget Experiment (ERBE). Both data set sponsors are aware of the utility of one data set to the other, and both projects utilize gridded data on a 2.5 deg by 2.5 deg grid. The ISCCP data have been collected since July 1983, and the NOAA-9 data for ERBE have been collected for more than a year. Therefore, there is a good chance to use these temporally overlapping data sets to investigate hypothesized relationships. Changes in cloudiness affect both cloud albedo feedback (shortwave) and the greenhouse effect (longwave). The relative importance of the effects of clouds on albedo versus outgoing longwave radiation (OLR) in determining the radiation balance has long been a matter of controversy. Now, however, changes in cloud amount as observed by the ISCCPO can be correlated to corresponding changes in the albedo and changes in the OLR from ERBE. Monthly means can be utilized in all instances.

Jacobowitz, H.↗

Embedded adaptation for time dependent real fluid/gas flows

An embedded adaptation algorithm is used to illustrate applications and requirements for Navier-Stokes and nonequilibrium unsteady two-dimensional flow fields. Both grid refinement and equation modification are included for unstructured distributions of quadrilateral cells. Multiple refinement parameters for distinct physical phenomena serve to monitor the evolving nonuniform field and to recognize flow features in isolated and overlapping subdomains. Examples from subsonic, hypersonic, internal, and external fields show substantial advantages for the adaptive approach.

Baron, Judson R.↗

A De-centralized Scheduling and Load Balancing Algorithm for Heterogeneous Grid Environments

In the past two decades, numerous scheduling and load balancing techniques have been proposed for locally distributed multiprocessor systems. However, they all suffer from significant deficiencies when extended to a Grid environment: some use a centralized approach that renders the algorithm unscalable, while others assume the overhead involved in searching for appropriate resources to be negligible. Furthermore, classical scheduling algorithms do not consider a Grid node to be N-resource rich and merely work towards maximizing the utilization of one of the resources. In this paper, we propose a new scheduling and load balancing algorithm for a generalized Grid model of N-resource nodes that not only takes into account the node and network heterogeneity, but also considers the overhead involved in coordinating among the nodes. Our algorithm is decentralized, scalable, and overlaps the node coordination time with that of the actual processing of ready jobs, thus saving valuable clock cycles needed for making decisions. The proposed algorithm is studied by conducting simulations using the Message Passing Interface (MPI) paradigm.

Arora, Manish↗