Search NASA⌕ Search

SEARCH · Search NASA

Results for “Regular 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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Analysis models for the estimation of oceanic fields

A general model for statistically optimal estimates is presented for dealing with scalar, vector and multivariate datasets. The method deals with anisotropic fields and treats space and time dependence equivalently. Problems addressed include the analysis, or the production of synoptic time series of regularly gridded fields from irregular and gappy datasets, and the estimate of fields by compositing observations from several different instruments and sampling schemes. Technical issues are discussed, including the convergence of statistical estimates, the choice of representation of the correlations, the influential domain of an observation, and the efficiency of numerical computations.

Carter, E. F.↗
NASA NTRS1989

Data traffic reduction schemes for Cholesky factorization on asynchronous multiprocessor systems

Communication requirements of Cholesky factorization of dense and sparse symmetric, positive definite matrices are analyzed. The communication requirement is characterized by the data traffic generated on multiprocessor systems with local and shared memory. Lower bound proofs are given to show that when the load is uniformly distributed the data traffic associated with factoring an n x n dense matrix using n to the alpha power (alpha less than or equal 2) processors is omega(n to the 2 + alpha/2 power). For n x n sparse matrices representing a square root of n x square root of n regular grid graph the data traffic is shown to be omega(n to the 1 + alpha/2 power), alpha less than or equal 1. Partitioning schemes that are variations of block assignment scheme are described and it is shown that the data traffic generated by these schemes are asymptotically optimal. The schemes allow efficient use of up to O(n to the 2nd power) processors in the dense case and up to O(n) processors in the sparse case before the total data traffic reaches the maximum value of O(n to the 3rd power) and O(n to the 3/2 power), respectively. It is shown that the block based partitioning schemes allow a better utilization of the data accessed from shared memory and thus reduce the data traffic than those based on column-wise wrap around assignment schemes.

Naik, Vijay K.↗
NASA NTRS1989

Methods for local gravity field approximation

The most widely known modern method for estimating gravity field values from observed data is least-squares collocation. Its advantages are that it can make estimates at arbitrary locations based on irregularly spaced observations, and that it makes use of statistical information about errors in the input data while providing corresponding information about the quality of the output estimates. Disadvantages of collocation include the necessity of inverting square matrices of dimension equal to the number of data values and the need to assume covariance models for the gravity field and the data errors. Fourier methods are an important alternative to collocation; having the advantage of greater computational efficiency, but requiring data estimates to be on a regular grid and not using or providing statistical accuracy information. The GEOFAST algorithm is an implementation of collocation that achieves high computational efficiency by transforming the estimation equations into the frequency domain where an accurate approximation may be made to reduce the workload. The forward and inverse Fast Fourier Transforms (FFTs) are utilized. The accuracy and computational efficiency of the GEOFAST algorithm is demonstrated using two sets of synthetic gravity data: marine gravity for an ocean trench region including wavelengths longer than 200 km; and local land gravity containing wavelengths as short as 5 km. These results are discussed along with issues such as the advantages of first removing reference field models before carrying out the estimation algorithm.

Sailor, R. V.↗
NASA NTRS1989

Fourier transform methods in local gravity modeling

New algorithms were derived for computing terrain corrections, all components of the attraction of the topography at the topographic surface and the gradients of these attractions. These algoriithms utilize fast Fourier transforms, but, in contrast to methods currently in use, all divergences of the integrals are removed during the analysis. Sequential methods employing a smooth intermediate reference surface were developed to avoid the very large transforms necessary when making computations at high resolution over a wide area. A new method for the numerical solution of Molodensky's problem was developed to mitigate the convergence difficulties that occur at short wavelengths with methods based on a Taylor series expansion. A trial field on a level surface is continued analytically to the topographic surface, and compared with that predicted from gravity observations. The difference is used to compute a correction to the trial field and the process iterated. Special techniques are employed to speed convergence and prevent oscillations. Three different spectral methods for fitting a point-mass set to a gravity field given on a regular grid at constant elevation are described. Two of the methods differ in the way that the spectrum of the point-mass set, which extends to infinite wave number, is matched to that of the gravity field which is band-limited. The third method is essentially a space-domain technique in which Fourier methods are used to solve a set of simultaneous equations.

Harrison, J. C.↗
NASA NTRS1990

Continuous assimilation of Geosat altimetric sea level observations into a numerical synoptic ocean model of the California Current

The Geosat altimetric sea level observations for the period from January to December 1987 were continuously assimilated into a realistic wind-driven numerical synoptic ocean model of the California Current in order to evaluate the effectiveness of using a realistic synoptic ocean model for interpolating (dynamically) real altimetric sea level observations onto a regular grid. The accuracy of dynamical interpolation was tested by comparing the gridded sea level residuals to ones estimated from in situ observations (by expendable bathythermographs) collected in the California Current region during the same period. The comparison yielded nearly exact agreement at low frequency (i.e., semiannual cycle), but less agreement on month-to-month time scales of variability, possibly due to the unfiltered nature of the in situ estimates.

White, Warren B.↗
NASA NTRS1990

The role of shocks in NGC 6240

The case of the infrared-bright system NGC 6240 encapsulates in many ways the major issues in understanding processes in infrared-bright systems in general. Its morphology strongly suggests that this is a merger in progress, with energetic phenomena being triggered as evidenced by strong radio emission and optical emission lines as well as the strong far-infrared output. Evidence has been produced supporting various schemes for producing the energy being radiated by this object. Early work on the optical spectrum suggested a Seyfert 2 classification, but the ionization level and spatial extent of the emission implied a more distributed source of energy. IR line studies found evidence for young supergiants, and for immense amounts of H2 excited by low-velocity shocks, which might contribute to a large enhancement in star-formation rate. Finally, the radio-continuum morphology shows compact sources possibly associated with active nuclei in both remnants of the merger. This is one of the nearest IR-bright galaxies objects at high luminosity (Cz = 7500 km s(exp -1), so that spatially resolved studies can yield further insight into which of these processes are most important - star formation, nuclear activity, or shock excitation. The velocity field, linewidth distribution, and ionization structure of the gas in NGC 6240 was mapped with the Dense-Pak fiber-optic array on the Kitt Peak 4-m telescope. In all, 135 spectra were obtained through an array of 2 inch apertures, with considerable overlap between many of the adjacent apertures produced by offsetting the telescope between the exposures. As the strongest features with lambda less than 1.1 ma, H alpha and the adjacent (NII) lines were observed, along with the neighboring (O I) and (S II) features. Maps of velocity and linewidth were constructed by interpolating measurements from these spectra onto a regular grid, and show irregular patterns with large velocity gradients and, particularly, an extended region of large local (shear?) velocity width.

Keel, William C.↗
NASA NTRS1991

Atmospheric water vapor transport: Estimation of continental precipitation recycling and parameterization of a simple climate model

The advective transport of atmospheric water vapor and its role in global hydrology and the water balance of continental regions are discussed and explored. The data set consists of ten years of global wind and humidity observations interpolated onto a regular grid by objective analysis. Atmospheric water vapor fluxes across the boundaries of selected continental regions are displayed graphically. The water vapor flux data are used to investigate the sources of continental precipitation. The total amount of water that precipitates on large continental regions is supplied by two mechanisms: (1) advection from surrounding areas external to the region; and (2) evaporation and transpiration from the land surface recycling of precipitation over the continental area. The degree to which regional precipitation is supplied by recycled moisture is a potentially significant climate feedback mechanism and land surface-atmosphere interaction, which may contribute to the persistence and intensification of droughts. A simplified model of the atmospheric moisture over continents and simultaneous estimates of regional precipitation are employed to estimate, for several large continental regions, the fraction of precipitation that is locally derived. In a separate, but related, study estimates of ocean to land water vapor transport are used to parameterize an existing simple climate model, containing both land and ocean surfaces, that is intended to mimic the dynamics of continental climates.

Brubaker, Kaye L.↗
NASA NTRS1995

Test particle study of minor ions in solar wind turbulence

We perform a parameter study of the temporal evolution of a test particle distribution function in MHD turbulence. The turbulent fields are calculated using a pseudo-spectral method and periodic boundary conditions on a regular grid of 180(exp 3) points, appropriate for incompressible, homogeneous and isotropic turbulence. Initially, the kinetic and the magnetic energy are equal on the average. Both, deterministic and random initial conditions are used, in the former case with zeros of the magnetic field located at grid points, in the latter case located by interpolation between grid points. The evolution of the minor ion distribution function is studied in detail as these turbulent fields evolve, developing strong current and vorticity sheets. Using the full collisionless equation of motion for the test particles, the efficiency of nonlinear interactions can be studied. The results are compared to theoretical predictions and are then discussed in connection with the observations of the dynamical properties of solar wind minor ions derived from in situ observations.

Zurbuchen, Th.↗
NASA NTRS1996

Multilevel Algorithm for Atmospheric Data Assimilation

A multiscale algorithm for the problem of optimal statistical interpolation of observed data has been developed. This problem includes the calculation of the vector of the 'analyzed' (best estimated) atmosphere flow field w(sup a) by the formula: w(sup a) = w(sup f) + P(sup f) H(sup T) y, where the quantity y is defined by the equation (H P(sup f) H(sup T) + R)y = w(sup o) - H w(sup f), using the given model forecast first guess w(sup f) and the vector of observations w(sup o); H is an interpolation operator from the regular grid to the observation network, P(sup f) is the forecast error covariance matrix, and R is the observation error covariance matrix. At this initial stage the case of univariate analysis of single level radiosonde height data is considered. The matrix R is assumed to be diagonal, and the matrix P(sup f) is assumed to be given by the formula P(sub ij)(sup f) = sigma(sub i)(sup f) mu(sub ij) sigma(sub j)(sub f), where mu(sub ij) is a smooth, decreasing function of the distance between the i-th and the j-th points. In this paper we describe a multiscale iterative process based on a multiresolution, simultaneous displacement technique and a localized variational calculation of iteration parameters.

Brandt, Achi↗
NASA NTRS1998

Incorporating Parallel Computing into the Goddard Earth Observing System Data Assimilation System (GEOS DAS)

Atmospheric data assimilation is a method of combining actual observations with model forecasts to produce a more accurate description of the earth system than the observations or forecast alone can provide. The output of data assimilation, sometimes called the analysis, are regular, gridded datasets of observed and unobserved variables. Analysis plays a key role in numerical weather prediction and is becoming increasingly important for climate research. These applications, and the need for timely validation of scientific enhancements to the data assimilation system pose computational demands that are best met by distributed parallel software. The mission of the NASA Data Assimilation Office (DAO) is to provide datasets for climate research and to support NASA satellite and aircraft missions. The system used to create these datasets is the Goddard Earth Observing System Data Assimilation System (GEOS DAS). The core components of the the GEOS DAS are: the GEOS General Circulation Model (GCM), the Physical-space Statistical Analysis System (PSAS), the Observer, the on-line Quality Control (QC) system, the Coupler (which feeds analysis increments back to the GCM), and an I/O package for processing the large amounts of data the system produces (which will be described in another presentation in this session). The discussion will center on the following issues: the computational complexity for the whole GEOS DAS, assessment of the performance of the individual elements of GEOS DAS, and parallelization strategy for some of the components of the system.

Larson, Jay W.↗
NASA NTRS1998

I/O Parallelization for the Goddard Earth Observing System Data Assimilation System (GEOS DAS)

The National Aeronautics and Space Administration (NASA) Data Assimilation Office (DAO) at the Goddard Space Flight Center (GSFC) has developed the GEOS DAS, a data assimilation system that provides production support for NASA missions and will support NASA's Earth Observing System (EOS) in the coming years. The GEOS DAS will be used to provide background fields of meteorological quantities to EOS satellite instrument teams for use in their data algorithms as well as providing assimilated data sets for climate studies on decadal time scales. The DAO has been involved in prototyping parallel implementations of the GEOS DAS for a number of years and is now embarking on an effort to convert the production version from shared-memory parallelism to distributed-memory parallelism using the portable Message-Passing Interface (MPI). The GEOS DAS consists of two main components, an atmospheric General Circulation Model (GCM) and a Physical-space Statistical Analysis System (PSAS). The GCM operates on data that are stored on a regular grid while PSAS works with observational data that are scattered irregularly throughout the atmosphere. As a result, the two components have different data decompositions. The GCM is decomposed horizontally as a checkerboard with all vertical levels of each box existing on the same processing element(PE). The dynamical core of the GCM can also operate on a rotated grid, which requires communication-intensive grid transformations during GCM integration. PSAS groups observations on PEs in a more irregular and dynamic fashion.

Lucchesi, Rob↗
NASA NTRS1998

A Catchment-Based Land Surface Model for GCMs and the Framework for its Evaluation

A new GCM-scale land surface modeling strategy that explicitly accounts for subgrid soil moisture variability and its effects on evaporation and runoff is now being explored. In a break from traditional modeling strategies, the continental surface is disaggregated into a mosaic of hydrological catchments, with boundaries that are not dictated by a regular grid but by topography. Within each catchment, the variability of soil moisture is deduced from TOP-MODEL equations with a special treatment of the unsaturated zone. This paper gives an overview of this new approach and presents the general framework for its off-line evaluation over North-America.

Ducharen, A.↗
NASA NTRS1999

The Kalman Filter and High Performance Computing at NASA's Data Assimilation Office (DAO)

Atmospheric data assimilation is a method of combining actual observations with model simulations to produce a more accurate description of the earth system than the observations alone provide. The output of data assimilation, sometimes called "the analysis", are accurate regular, gridded datasets of observed and unobserved variables. This is used not only for weather forecasting but is becoming increasingly important for climate research. For example, these datasets may be used to assess retrospectively energy budgets or the effects of trace gases such as ozone. This allows researchers to understand processes driving weather and climate, which have important scientific and policy implications. The primary goal of the NASA's Data Assimilation Office (DAO) is to provide datasets for climate research and to support NASA satellite and aircraft missions. This presentation will: (1) describe ongoing work on the advanced Kalman/Lagrangian filter parallel algorithm for the assimilation of trace gases in the stratosphere; and (2) discuss the Kalman filter in relation to other presentations from the DAO on Four Dimensional Data Assimilation at this meeting. Although the designation "Kalman filter" is often used to describe the overarching work, the series of talks will show that the scientific software and the kind of parallelization techniques that are being developed at the DAO are very different depending on the type of problem being considered, the extent to which the problem is mission critical, and the degree of Software Engineering that has to be applied.

Lyster, Peter M.↗
NASA NTRS2001

Global Discrete Artificial Boundary Conditions for Time-Dependent Wave Propagation

We construct global artificial boundary conditions (ABCs) for the numerical simulation of wave processes on unbounded domains using a special non-deteriorating algorithm that has been developed previously for the long-term computation of wave-radiation solutions. The ABCs are obtained directly for the discrete formulation of the problem; in so doing, neither a rational approximation of 'non-reflecting kernels,' nor discretization of the continuous boundary conditions is required. The extent of temporal nonlocality of the new ABCs appears fixed and limited; in addition, the ABCs can handle artificial boundaries of irregular shape on regular grids with no fitting/adaptation needed and no accuracy loss induced. The non-deteriorating algorithm, which is the core of the new ABCs is inherently three-dimensional, it guarantees temporally uniform grid convergence of the solution driven by a continuously operating source on arbitrarily long time intervals, and provides unimprovable linear computational complexity with respect to the grid dimension. The algorithm is based on the presence of lacunae, i.e., aft fronts of the waves, in wave-type solutions in odd-dimension spaces, It can, in fact, be built as a modification on top of any consistent and stable finite-difference scheme, making its grid convergence uniform in time and at the same time keeping the rate of convergence the same as that of the non-modified scheme. In the paper, we delineate the construction of the global lacunae-based ABCs in the framework of a discretized wave equation. The ABCs are obtained for the most general formulation of the problem that involves radiation of waves by moving sources (e.g., radiation of acoustic waves by a maneuvering aircraft). We also present systematic numerical results that corroborate the theoretical design properties of the ABCs' algorithm.

Ryabenkii, V. S.↗
NASA NTRS2002

Surface Temperature Assimilation in the Global Land Data Assimilation System (GLDAS)

The Global Land Data Assimilation System (GLDAS) is a global land parameterization that uses prescribed meteorology as forcing in order to determine regular gridded land surface states (temperature and moisture) and other properties (e.g. water and heat fluxes). In the present experiment, the assimilation of surface skin temperature is incorporated into the land parameterizations. The meteorological forcing was derived from the Goddard Earth Observing System (GEOS-3) Data Assimilation System (DAS) for the full year of 1998 GLDAS can use several land parameterizations, but here we use the Mosaic land surface model and the Common Land Model (CLM). TOVS surface temperature observations are assimilated into GLDAS. The TOVS observations are less frequent that observations used in previous experiments (ISCCP). The purpose of this presentation is to evaluate the impact of the TOVS assimilation on both Mosaic and CLM. We will especially consider the impact of coarse temporal observations on the assimilation and bias correction.

Bosilovich, Michael G.↗
NASA NTRS2017

Climate Data Guide - Modern Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2)

The Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) is a global atmospheric reanalysis produced by the NASA Global Modeling and Assimilation Office (GMAO). It spans the satellite observing era from 1980 to the present. The goals of MERRA-2 are to provide a regularly-gridded, homogeneous record of the global atmosphere, and to incorporate additional aspects of the climate system including trace gas constituents (stratospheric ozone), and improved land surface representation, and cryospheric processes. MERRA-2 is also the first satellite-era global reanalysis to assimilate space-based observations of aerosols and represent their interactions with other physical processes in the climate system. The inclusion of these additional components are consistent with the overall objectives of an Integrated Earth System Analysis (IESA). MERRA-2 is intended to replace the original MERRA product, and reflects recent advances in atmospheric modeling and data assimilation. Modern hyperspectral radiance and microwave observations, along with GPS-Radio Occultation and NASA ozone datasets are now assimilated in MERRA-2. Much of the structure of the data files remains the same in MERRA-2. While the original MERRA data format was HDF-EOS, the MERRA-2 supplied binary data format is now NetCDF4 (with lossy compression to save space).

MERRA-2↗
NASA NTRS2019

Numerical Reanalyses as a Gateway to Arctic Synthesis

Reanalyses are regularly gridded, retrospective depictions of the physical earth system, which are produced through the correction of a short-term forecast to available observations. In the Arctic, reanalyses are particularly well suited to marshal the sparse observing network to provide a plausible, multivariate representation of conditions. Atmospheric reanalyses such as MERRA-2 (NASA Modern-Era Retrospective analysis for Research and Applications, version 2) and ocean reanalyses such as SODA3 (Univ. Maryland Simple Ocean Data Assimilation version 3) are widely used in Arctic research for diagnostic studies of circulation, model evaluation, and as boundary conditions for a variety of process models. Here, we provide examples that illustrate the utility of reanalyses for providing information on the spatial and temporal scales of recent, rapid changes in the Arctic. Recent trends in Arctic surface temperatures, surface melt over Greenland and Arctic glaciers, and evolving freshwater conditions in the Arctic Ocean are examples where reanalyses can provide information that cannot easily be obtained via other means. These examples provide information on the scale, magnitude, and the uncertainty of recent Arctic change and provide a context for future scenarios. We further quantify uncertainties in key reanalyses variables and approaches for addressing these issues.

Arctic↗
NASA NTRSDate not supplied

Expanding the Time Dimension of Hyperspectral Infrared Sounding Observations: Designing the NUCAPS-Forecast System

Low Earth orbit (LEO) hyper-spectral infrared (IR) sounders have significant yet untapped potential for characterizing thermodynamic environments of convective initiation and ongoing convection. While LEO soundings are of enormous value to weather forecasters, they cannot provide the temporal resolution needed to resolve the rapidly evolving thermodynamics of the convective environment. A novel methodology was developed to extend snapshots of LEO soundings forward in time up to six hours into the future and has been made available within the NOAA National Weather Service Advanced Weather Interactive Processing System (AWIPS) for testing in a quasi-operational setting. This methodology is based on parcel forward-trajectory calculations from the satellite observing time to recreate future soundings of temperature and moisture at regularly gridded intervals in space and time. The soundings are based on NOAA-Unique Combined Atmospheric Processing System (NUCAPS) retrievals from the Suomi NPP and NOAA-20 satellite platforms to allow for rapid near-real time product processing. Science code was transitioned to operational production in 2019 and the product, NUCAPS-Forecast, was tested at the NOAA Hazardous Weather Testbed Spring Experiment in 2019 and 2021. Forecaster feedback has led to improving the product and processing system to address product limitations and improve its applicability for pre-convective forecasting. This presentation highlights the design of the NUCAPS-Forecast system. Solutions were tested and implemented to remove excess moisture and conjoin forward-trajected parcels into vertical profiles for calculating grids of stability indices. Additionally, use cases that explore the tendencies of derived convective available potential energy (CAPE) and convective inhibition (CIN) were investigated. When evaluated against gridded, hourly accumulated rainfall obtained from the Multi-Radar Multi-Sensor (MRMS) observational data set areas with elevated CAPE (reduced CIN) are shown to be associated with areas of precipitation and increase confidence the LEO soundings advected forward in time physically represent atmospheric processes. These results suggest that adiabatic parcel advection of LEO satellite sounding snapshots forward in time exhibits skill in discriminating locations and times that are more likely to lead to convective initiation. Leveraging trajectory modeling to increase the spatial and temporal resolution of LEO IR sounding observations today demonstrates the value of these observations from the future constellation of JPSS satellites and GeoXO mission.

weather↗