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

Geoid-to-topography ratios on Venus: A global perspective

Recently available spherical harmonic solutions for the geoid and topography of Venus are sufficiently high resolution that they can be used to address questions concerning the relationship between geoid and topography on a regional scale. We have approached this question by mapping the geoid-to-topography ratio (GTR) on a systematic global basis. For a given point on the surface, we consider the geoid and elevation values at all points on a gridded representation of those fields located within a specified distance of the reference point. From the set of paired values, we determine the correlation coefficient and the best-fitting straight line. The latter is the GTR at that position, and the former is a measure of the significance of the derived ratio. This procedure is then repeated for all points on the global grid, yielding maps of the GTR and the correlation coefficient. Unlike previous studies of the GRT on Venus, this apprach permits us to make an objective and systematic search for regions with anomalous GTR's as well as areas that do not demonstrate any strong correlation between geoid and topography. These maps can be updated regularly as new harmonic models of the Venus geoid are produced from new Magellan tracking data. This procedure permits the development of a global perspective on the relationship between GTR and venusian surface tectonics.

Simons, Mark↗

Climate Reanalysis: Progress and Future Prospects

Reanalysis is the process whereby an unchanging data assimilation system is used to provide a consistent reprocessing of observations, typically spanning an extended segment of the historical data record. The process relies on an underlying model to combine often-disparate observations in a physically consistent manner, enabling production of gridded data sets for a broad range of applications including the study of historical weather events, preparation of climatologies, business sector development and, more recently, climate monitoring. Over the last few decades, several generations of reanalyses of the global atmosphere have been produced by various operational and research centers, focusing more or less on the period of regular conventional and satellite observations beginning in the mid to late twentieth century. There have also been successful efforts to extend atmospheric reanalyses back to the late nineteenth and early twentieth centuries, using mostly surface observations. Much progress has resulted from (and contributed to) advancements in numerical weather prediction, especially improved models and data assimilation techniques, increased computing capacity, the availability of new observation types and efforts to recover and improve the quality of historical ones. The recent extension of forecast systems that allow integrated modeling of meteorological, oceanic, land surface, and chemical variables provide the basic elements for coupled data assimilation. This has opened the door to the development of a new generation of coupled reanalyses of the Earth system, or integrated Earth system analyses (IESA). Evidence so far suggests that this approach can improve the analysis of currently uncoupled components of the Earth system, especially at their interface, and lead to increased predictability. However, extensive analysis coupling as envisioned for IESA, while progressing, still presents significant challenges. These include model biases that can be exacerbated when coupled, component systems with different physical characteristics and different spatial and temporal scales, and component observations in different media with different spatial and temporal frequencies and different latencies. Quantification of uncertainty in reanalyses is also a critical challenge and is important for expanding their utility as a tool for climate change assessment. This talk provides a brief overview of the progress of reanalysis development during recent decades, and describes remaining challenges in the progression toward coupled Earth system reanalyses.

IES↗

The GISMO-2 Bolometer Camera

We present the concept for the GISMO-2 bolometer camera) which we build for background-limited operation at the IRAM 30 m telescope on Pico Veleta, Spain. GISM0-2 will operate Simultaneously in the 1 mm and 2 mm atmospherical windows. The 1 mm channel uses a 32 x 40 TES-based Backshort Under Grid (BUG) bolometer array, the 2 mm channel operates with a 16 x 16 BUG array. The camera utilizes almost the entire full field of view provided by the telescope. The optical design of GISM0-2 was strongly influenced by our experience with the GISMO 2 mm bolometer camera which is successfully operating at the 30m telescope. GISMO is accessible to the astronomical community through the regular IRAM call for proposals.

Staguhn, Johannes G.↗

Some Numerical Experiments on Detonation Wave Propagation

In this paper we present the results of a series of numerical experiments done on the propagation and initiation of a detonation wave. The calculations are performed in one-dimension, with considerable grid resolution. Of particular interest are the following questions: (1) the nature of periodic and chaotic instabilities generated by the wave; (2) the influence of the grid resolution on these instabilities; (3) the influence of the 'quality' of the numerical scheme; and (4) the influence of 'noise'. In the calculations, we use a second-order Total Variation Diminishing (TVD) scheme as the basic numerical method, with grid spacings as low as a fraction of a micron. Detonations waves are generated at the closed end of a tube, and allowed to propagate for approximately 20 cm. The required energy for successful initiation of the detonation will be measured for different cases of grid resolution and numerical schemes. A modified version of the TVD scheme has also been devised, which allows for much lower numerical diffusion of the radical species in the exponentially growing region behind the shock. The effect of this modification will be demonstrated. Oscillations in peak pressure and induction length are seen to develop in some cases: the oscillations can go through a sequence of modes, from a regular, high frequency mode to a low frequency mode with period doubling. A chaotic regime can also be obtained. General conclusions on the quality of algorithms will be presented. We will also discuss the performance of a version of the code developed on the IBM SP2 parallel computer.

Cambier, Jean-Luc↗

Ecosystem age-class dynamics and distribution in the LPJ-wsl v2.0 global ecosystem model

Forest ecosystem processes follow classic responses with age, peaking production around canopy closure and declining thereafter. Although age dynamics might be more dominant in certain regions over others, demographic effects on net primary production (NPP) and heterotrophic respiration (Rh) are bound to exist. Yet, explicit representation of ecosystem demography is notably absent in many global ecosystem models. This is concerning because the global community relies on these models to regularly update our collective understanding of the global carbon cycle. This paper aims to present the technical developments of a computationally efficient approach for representing age-class dynamics within a global ecosystem model, the Lund–Potsdam–Jena – Wald, Schnee, Landschaft version 2.0 (LPJ-wsl v2.0) dynamic global vegetation model and to determine if explicit representation of demography influenced ecosystem stocks and fluxes at global scales or at the level of a grid cell. The modeled age classes are initially created by simulated fire and prescribed wood harvesting or abandonment of managed land, otherwise aging naturally until an additional disturbance is simulated or prescribed. In this paper, we show that the age module can capture classic demographic patterns in stem density and tree height compared to inventory data, and that simulated patterns of ecosystem function follow classic responses with age. We also present two scientific applications of the model to assess the modeled age-class distribution over time and to determine the demographic effect on ecosystem fluxes relative to climate. Simulations show that, between 1860 and 2016, zonal age distribution on Earth was driven predominately by fire, causing a 45- to 60-year difference in ages between older boreal (50–90° N) and younger tropical (23° S–23° N) ecosystems. Between simulation years 1860 and 2016, land-use change and land management were responsible for a decrease in zonal age by −6 years in boreal and by −21 years in both temperate (23–50° N) and tropical latitudes, with the anthropogenic effect on zonal age distribution increasing over time. A statistical model helped to reduce LPJ-wsl v2.0 complexity by predicting per-grid-cell annual NPP and Rh fluxes by three terms: precipitation, temperature, and age class; at global scales, R2 was between 0.95 and 0.98. As determined by the statistical model, the demographic effect on ecosystem function was often less than 0.10 kg C/sq. myr but as high as 0.60 kg C/sq. myr where the effect was greatest. In the eastern forests of North America, the simulated demographic effect was of similar magnitude, or greater than, the effects of climate; simulated demographic effects were similarly important in large regions of every vegetated continent. Simulated spatial datasets are provided for global ecosystem ages and the estimated coefficients for effects of precipitation, temperature and demography on ecosystem function. The discussion focuses on our finding of an increasing role of demography in the global carbon cycle, the effect of demography on relaxation times (resilience) following a disturbance event and its implications at global scales, and a finding of a 40 Pg C increase in biomass turnover when including age dynamics at global scales. Whereas time is the only mechanism that increases ecosystem age, any additional disturbance not explicitly modeled will decrease age. The LPJ-wsl v2.0 age module represents another step forward towards understanding the role of demography in global ecosystems.

LPJ-wsl↗

Computations of Unsteady Viscous Compressible Flows Using Adaptive Mesh Refinement in Curvilinear Body-fitted Grid Systems

A methodology for accurate and efficient simulation of unsteady, compressible flows is presented. The cornerstones of the methodology are a special discretization of the Navier-Stokes equations on structured body-fitted grid systems and an efficient solution-adaptive mesh refinement technique for structured grids. The discretization employs an explicit multidimensional upwind scheme for the inviscid fluxes and an implicit treatment of the viscous terms. The mesh refinement technique is based on the AMR algorithm of Berger and Colella. In this approach, cells on each level of refinement are organized into a small number of topologically rectangular blocks, each containing several thousand cells. The small number of blocks leads to small overhead in managing data, while their size and regular topology means that a high degree of optimization can be achieved on computers with vector processors.

Steinthorsson, E.↗

In-flight Far-Infrared Performance of the CIRS Instrument on Cassini

The Composite Infrared Spectrometer (CIRS) on-board Cassini consists of two interferometers: a conventional Michelson for the mid-infrared; and a Martin-Puplett type in the far-infrared employing wire grid polarizers to split, recombine and analyze the radiation. The far-IR focal plane (FP1) assembly uses two thermopile detectors to measure the final transmitted and reflected beams at the polarizer-analyzer: if one fails, the interferometer can still operate, albeit with a lower efficiency. The combined effect is for good response from 10 to 300/cm, and declining response to 600/cm. This paper will examine in-flight performance of the far-IR interferometer, including NESR and response. Regular noise spikes, resulting from pickup from other electrical sub-systems has been found on the CIRS interferograms, and the removal of these effects is discussed. The radiometric calibration is described, and then we show how the calibration was applied to science data taken during the Jupiter flyby of December 2000. Finally, we discuss signal-to-noise on the calibrated spectra, emphasizing limitations of the current instrument and the potential for improvement in future missions.

Nixon, Conor A.↗

Recent developments in FEM-CFD

The current status of CFD with regard to unstructured grids employing finite element methods and Eulerian frames is reviewed. Algorithms suitable for the computation of large three-dimensional problems involving flow past arbitrary geometries are developed. Adaptive mesh refinement strategy is reviewed, and domain splitting or local time-stepping are briefly addressed. The development of search algorithms of optimal order, variable time-stepping Jacobi smoothers for elliptic problems, and transport concepts for hyperbolics to help achieve good performance for unstructured multigrid processes is discussed. As examples, transient supersonic flow in a channel, regular shock reflection of a wall, viscous flow past a protruberance, potential flow past a cylinder, and Burgers equation are considered.

Loehner, R.↗

Interpolating Fields of Carbon Monoxide Data Using a Hybrid Statistical-Physical Model

Atmospheric Carbon Monoxide (CO) is a pollutant gas of which the US congress has mandated regular monitoring, and satellite sensors can be used to retrieve regional concentrations of CO over several vertical layers. However, CO at cloudy locations cannot be observed and have to be estimated from the observed data set, resulting in an interpolation problem. The current state-of-the-art solution is to combine prior information, computed by a deterministic physical model, with observations. However, the deterministic model may introduce uncertainties that do not derive from the data. While sharing certain features with the physical model, this paper presents a Bayesian hierarchical model for interpolating CO on a 3-dimensional spatial grid, across time. To our knowledge such a model has not been considered before. The model is applied to a hypothetical air-quality monitoring scenario, and is compared to existing interpolation methods. The results provide motivation for the use of the statistical model for regional to local applications.

Arellano, A. A.↗

Multigrid methods for bifurcation problems: The self adjoint case

This paper deals with multigrid methods for computational problems that arise in the theory of bifurcation and is restricted to the self adjoint case. The basic problem is to solve for arcs of solutions, a task that is done successfully with an arc length continuation method. Other important issues are, for example, detecting and locating singular points as part of the continuation process, switching branches at bifurcation points, etc. Multigrid methods have been applied to continuation problems. These methods work well at regular points and at limit points, while they may encounter difficulties in the vicinity of bifurcation points. A new continuation method that is very efficient also near bifurcation points is presented here. The other issues mentioned above are also treated very efficiently with appropriate multigrid algorithms. For example, it is shown that limit points and bifurcation points can be solved for directly by a multigrid algorithm. Moreover, the algorithms presented here solve the corresponding problems in just a few work units (about 10 or less), where a work unit is the work involved in one local relaxation on the finest grid.

Taasan, Shlomo↗

Total Variation Majorization Minimization (TV-MM) Approach to Radiometer Brightness Temperature Gridding and Reconstruction

This paper presents the implementation of an algorithm to enhance the image resolution of the Earth's surface brightness temperature (T B ) data measured by radiometers such as the one onboard of the Soil Moisture Active Passive (SMAP) mission. A key step in radiometer T B processing is the conversion of the swath-based calibrated antenna temperature (T A ) measurements to the Level 3 Earth-centered grid. The simplest algorithm to transform this data from swath to gridded format is called drop-in-the-bucket which simply averages surrounding noisy T A samples to form a T B value at the gridded location. This method reduces noise, however produces low resolution products. To obtain a higher resolution product, SMAP uses other techniques such the Backus-Gilbert (BG) algorithm, which is the conventional method used in microwave radiometry. Although this method performs the required interpolation, it is not effective in denoising and removing blurring effects due to antenna filtering of the radiometer image data. Our motivation for this development is to further improve the resolution through post-processing of the radiometer T B image, a highly cost-effective method of image enhancement. The approach adapted in this work is based on the minimization of the Total Variation (TV) regularized objective function that is used extensively in solving general ill-posed linear inverse problems in image processing. Since the TV-based objective function is convex but not everywhere differentiable, there exists many numerical algorithms that can estimate the solution and the one selected for this work is called Majorization- Minimization (MM). By applying this algorithm, simulation experiments were performed based on synthetic data from the Geophysical model as well as real SMAP data to demonstrate the effectiveness of the technique. Results were then compared against the BG method.

Wing Lee↗

A NASA GISTEMP Observational Uncertainty Ensemble: Regional and Monthly Uncertainty

The historical global temperature record is an essential data product for quantifying the variability and change of the Earth system. In recent years, better characterization of observational uncertainty in global and hemispheric trends has become available, but the methodologies are not necessarily applicable to analyses at smaller regional areas, or monthly means, where station sparsity and other systematic issues contribute to greater uncertainty. This work details a gridded uncertainty ensemble of historical temperature anomalies from the Goddard Institute for Space Studies (GISS) Surface Temperature product (GISTEMP) product. This ensemble characterizes the complex spatial and temporal correlation structure of uncertainty in gridded historical temperature, enabling proper uncertainty propagation for climate and social science at regional and monthly scales. This work details the methodology for generating the uncertainty ensemble, key statistics of the uncertainty evolution over space and time, and provides best practices for using the uncertainty ensemble in future studies. Summary statistics from the uncertainty ensemble are in good agreement with production GISTEMP. Two applications of the uncertainty ensemble are also presented. First, the warmest year on record is shown to most likely be 2016 with a 53.2% chance and 2020 as the second most likely with a 44.4% chance. Second, it is shown that the arctic is warming 2.5 - 5 times faster than the globe, significantly faster than the regularly quoted twice as fast.

GISTEMP↗

Stratospheric circulation in the southern summer/northern winter 1980-81 - Behavior of zonal waves 1-10

Daily global grids of high quality, satellite-borne microwave measurements have been used to study the behavior of stratospheric waves during November 1980-March 1981. The primary motivation was to investigate, by comparing with earlier analyses, the interannual variation of recently reported medium-scale wave domination of Southern Hemisphere (SH) summer circulation. Zonal means, as well as time-mean and transient zonal waves 1-10, are studied. In addition, space-time spectra are determined from the 138-day record. Prominent medium-scale features are found to occur in 1980-81, although wave 4 is strongest, in contrast to the Global Weather Experiment period (1978-79), when wave 5 dominated the SH summer circulation. The same highly continuous, regular eastward phase movement of SH medium-scale waves is found in 1980-81 as in 1978-79. A second noteworthy observation is that of strong, time-mean waves 3-4 with symmetry about the equator, with positive temperature perturbations over the western Pacific and Atlantic equatorial regions, and negative anomalies over South America, Africa, and Indonesia. Confidence in the SH spectral results is enhanced by good agreement between our Northern Hemisphere (NH) medium-scale analyses and recent geopotential height studies by other investigators. These results serve to bring into sharper focus the question of why, in contrast with the more wellknown NH, the SH circulation evidences such robust and highly continuous medium-scale eddies in summer.

Yu, W.-B.↗

Virtual and Experimental Visualization of Flows in Packed Beds of Spheres Simulating Porous Media Flows

A videotape presentation of flow in a packed bed of spheres is provided. The flow experiment consisted of three principal elements: (1) an oil tunnel 76.2 mm by 76.2 mm in cross section, (2) a packed bed of spheres in regular and irregular arrays, and (3) a flow characterization methodology, either (a) full flow field tracking (FFFT) or (b) computational fluid dynamic (CFD) simulation. The refraction indices of the oil and the test array of spheres were closely matched, and the flow was seeded with aluminum oxide particles. Planar laser light provided a two-dimensional projection of the flow field, and a traverse simulated a three-dimensional image of the entire flow field. Light focusing and reflection rendered the spheres black, permitting visualization of the planar circular interfaces in both the axial and transverse directions. Flows were observed near the wall-sphere interface and within the set of spheres. The CFD model required that a representative section of a packed bed be formed and gridded, enclosing and cutting six spheres so that symmetry conditions could be imposed at all cross-boundaries. Simulations had to be made with the flow direction at right angles to that used in the experiments, however, to take advantage of flow symmetry. Careful attention to detail was required for proper gridding. The flow field was three-dimensional and complex to describe, yet the most prominent finding was flow threads, as computed in the representative 'cube' of spheres with face symmetry and conclusively demonstrated experimentally herein. Random packing and bed voids tended to disrupt the laminar flow, creating vortices.

Hendricks, R. C.↗

The Computational Complexity, Parallel Scalability, and Performance of Atmospheric Data Assimilation Algorithms

The computational complexity of algorithms for Four Dimensional Data Assimilation (4DDA) at NASA's Data Assimilation Office (DAO) is discussed. In 4DDA, observations are assimilated with the output of a dynamical model to generate best-estimates of the states of the system. It is thus a mapping problem, whereby scattered observations are converted into regular accurate maps of wind, temperature, moisture and other variables. The DAO is developing and using 4DDA algorithms that provide these datasets, or analyses, in support of Earth System Science research. Two large-scale algorithms are discussed. The first approach, the Goddard Earth Observing System Data Assimilation System (GEOS DAS), uses an atmospheric general circulation model (GCM) and an observation-space based analysis system, the Physical-space Statistical Analysis System (PSAS). GEOS DAS is very similar to global meteorological weather forecasting data assimilation systems, but is used at NASA for climate research. Systems of this size typically run at between 1 and 20 gigaflop/s. The second approach, the Kalman filter, uses a more consistent algorithm to determine the forecast error covariance matrix than does GEOS DAS. For atmospheric assimilation, the gridded dynamical fields typically have More than 10(exp 6) variables, therefore the full error covariance matrix may be in excess of a teraword. For the Kalman filter this problem can easily scale to petaflop/s proportions. We discuss the computational complexity of GEOS DAS and our implementation of the Kalman filter. We also discuss and quantify some of the technical issues and limitations in developing efficient, in terms of wall clock time, and scalable parallel implementations of the algorithms.

Lyster, Peter M.↗

Modifications to Axially Symmetric Simulations Using New DSMC (2007) Algorithms

Several modifications aimed at improving physical accuracy are proposed for solving axially symmetric problems building on the DSMC (2007) algorithms introduced by Bird. Originally developed to solve nonequilibrium, rarefied flows, the DSMC method is now regularly used to solve complex problems over a wide range of Knudsen numbers. These new algorithms include features such as nearest neighbor collisions excluding the previous collision partners, separate collision and sampling cells, automatically adaptive variable time steps, a modified no-time counter procedure for collisions, and discontinuous and event-driven physical processes. Axially symmetric solutions require radial weighting for the simulated molecules since the molecules near the axis represent fewer real molecules than those farther away from the axis due to the difference in volume of the cells. In the present methodology, these radial weighting factors are continuous, linear functions that vary with the radial position of each simulated molecule. It is shown that how one defines the number of tentative collisions greatly influences the mean collision time near the axis. The method by which the grid is treated for axially symmetric problems also plays an important role near the axis, especially for scalar pressure. A new method to treat how the molecules are traced through the grid is proposed to alleviate the decrease in scalar pressure at the axis near the surface. Also, a modification to the duplication buffer is proposed to vary the duplicated molecular velocities while retaining the molecular kinetic energy and axially symmetric nature of the problem.

Liechty, Derek S.↗

Visualization of Atmospheric Water Vapor Data for SAGE

The goal of this project was to develop visualization tools to study the water vapor dynamics using the Stratospheric Aerosol and Gas Experiment 11 (SAGE 11) water vapor data. During the past years, we completed the development of a visualization tool called EZSAGE, and various Gridded Water Vapor plots, tools deployed on the web to provide users with new insight into the water vapor dynamics. Results and experiences from this project, including papers, tutorials and reviews were published on the main Web page. Additional publishing effort has been initiated to package EZSAGE software for CD production and distribution. There have been some major personnel changes since Fall, 1998. Dr. Mou-Liang Kung, a Professor of Computer Science assumed the PI position vacated by Dr. Waldo Rodriguez who was on leave. However, former PI, Dr. Rodriguez continued to serve as a research adviser to this project to assure smooth transition and project completion. Typically in each semester, five student research assistants were hired and trained. Weekly group meetings were held to discuss problems, progress, new research direction, and activity planning. Other small group meetings were also held regularly for different objectives of this project. All student research assistants were required to submit reports for conference submission.

Kung, Mou-Liang↗

Global-scale distribution of ozone in the remote troposphere from ATom and HIPPO airborne field missions

Ozone is a key constituent of the troposphere, where it drives photochemical processes, impacts air quality, and acts as a climate forcer. Large-scale in situ observations of ozone commensurate with the grid resolution of current Earth system models are necessary to validate model outputs and satellite retrievals. In this paper, we examine measurements from the Atmospheric Tomography (ATom; four deployments in 2016–2018) and the HIAPER Pole-to-Pole Observations (HIPPO; five deployments in 2009–2011) experiments, two global-scale airborne campaigns covering the Pacific and Atlantic basins. ATom and HIPPO represent the first global-scale, vertically resolved measurements of O3 distributions throughout the troposphere, with HIPPO sampling the atmosphere over the Pacific and ATom sampling both the Pacific and Atlantic. Given the relatively limited temporal resolution of these two campaigns, we first compare ATom and HIPPO ozone data to longer-term observational records to establish the representativeness of our dataset. We show that these two airborne campaigns captured on average 53 %, 54 %, and 38 % of the ozone variability in the marine boundary layer, free troposphere, and upper troposphere–lower stratosphere (UTLS), respectively, at nine well-established ozonesonde sites. Additionally, ATom captured the most frequent ozone concentrations measured by regular commercial aircraft flights in the northern Atlantic UTLS. We then use the repeated vertical profiles from these two campaigns to confirm and extend the existing knowledge of tropospheric ozone spatial and vertical distributions throughout the remote troposphere. We highlight a clear hemispheric gradient, with greater ozone in the Northern Hemisphere, consistent with greater precursor emissions and consistent with previous modeling and satellite studies. We also show that the ozone distribution below 8 km was similar in the extra-tropics of the Atlantic and Pacific basins, likely due to zonal circulation patterns. However, twice as much ozone was found in the tropical Atlantic as in the tropical Pacific, due to well-documented dynamical patterns transporting continental air masses over the Atlantic. Finally, we show that the seasonal variability of tropospheric ozone over the Pacific and the Atlantic basins is driven year-round by transported continental plumes and photochemistry, and the vertical distribution is driven by photochemistry and mixing with stratospheric air. This new dataset provides additional constraints for global climate and chemistry models to improve our understanding of both ozone production and loss processes in remote regions, as well as the influence of anthropogenic emissions on baseline ozone.

ATom, aircraft profiles, Atlantic pollution, ozone↗