Search NASA⌕ Search

SEARCH · Search NASA

Results for “gridded data”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

The Principle of Energetic Consistency: Application to the Shallow-Water Equations

If the complete state of the earth's atmosphere (e.g., pressure, temperature, winds and humidity, everywhere throughout the atmosphere) were known at any particular initial time, then solving the equations that govern the dynamical behavior of the atmosphere would give the complete state at all subsequent times. Part of the difficulty of weather prediction is that the governing equations can only be solved approximately, which is what weather prediction models do. But weather forecasts would still be far from perfect even if the equations could be solved exactly, because the atmospheric state is not and cannot be known completely at any initial forecast time. Rather, the initial state for a weather forecast can only be estimated from incomplete observations taken near the initial time, through a process known as data assimilation. Weather prediction models carry out their computations on a grid of points covering the earth's atmosphere. The formulation of these models is guided by a mathematical convergence theory which guarantees that, given the exact initial state, the model solution approaches the exact solution of the governing equations as the computational grid is made more fine. For the data assimilation process, however, there does not yet exist a convergence theory. This book chapter represents an effort to begin establishing a convergence theory for data assimilation methods. The main result, which is called the principle of energetic consistency, provides a necessary condition that a convergent method must satisfy. Current methods violate this principle, as shown in earlier work of the author, and therefore are not convergent. The principle is illustrated by showing how to apply it as a simple test of convergence for proposed methods.

Cohn, Stephen E.↗
NASA NTRSDate not supplied

Advancing NASA SatCORPS Global Data Products with Cloud Computing and Machine Learning

Operational satellite imager radiances are valuable for deriving many different physical parameters that can be used for a variety of weather, aviation, and energy applications. The NASA Satellite ClOud and Radiation Property retrieval System (SatCORPS) applies a suite of algorithms to meteorological satellite data to provide cloud properties, radiative fluxes and other parameters on a global scale. The use of cloud computing has enabled recent enhancements to process a constellation of geostationary satellites at higher spatiotemporal resolutions than previously possible that meet low latency and near real-time needs. Data taken from Meteosat-8 and -11, Himawari, GOES-16, and -17 are processed and combined with operational polar orbiting satellite data and composited on a 3-km grid to provide global coverage. To improve the utility of the data products, machine learning and other innovative methods are applied in various ways to help minimize data product uncertainties under the most challenging conditions and to improve their consistency at all times of day. An update on recent SatCORPS enhancements is presented, highlighting the community benefits achieved with the use of cloud computing and machine learning.

William L Smith↗
NASA NTRS1992

The NASA-IGES geometry data exchange standard

Described here are the data exchange efforts and plans supported by the NASA Steering Committee for Surface Modeling and Grid Generation. Current methods for geometry data exchange between computer aided design (CAD) systems and NASA computational fluid dynamics (CFD) analysis systems are tedious and induce errors. A Geometry Data Exchange Standard is proposed, utilizing a subset of an existing national standard titled Initial Graphic Exchange Standard (IGES). Future plans for data exchange standardization include all aspects of CFD data. Software systems to utilize this NASA-IGES Geometry Data Exchange Specification are under development.

Blake, Matthew W.↗
NASA NTRS2015

Rapid, High-Resolution Detection of Environmental Change over Continental Scales from Satellite Data - the Earth Observation Data Cube

The effort and cost required to convert satellite Earth Observation (EO) data into meaningful geophysical variables has prevented the systematic analysis of all available observations. To overcome these problems, we utilise an integrated High Performance Computing and Data environment to rapidly process, restructure and analyse the Australian Landsat data archive. In this approach, the EO data are assigned to a common grid framework that spans the full geospatial and temporal extent of the observations - the EO Data Cube. This approach is pixel-based and incorporates geometric and spectral calibration and quality assurance of each Earth surface reflectance measurement. We demonstrate the utility of the approach with rapid time-series mapping of surface water across the entire Australian continent using 27 years of continuous, 25 m resolution observations. Our preliminary analysis of the Landsat archive shows how the EO Data Cube can effectively liberate high-resolution EO data from their complex sensor-specific data structures and revolutionise our ability to measure environmental change.

Lewis, Adam↗
NASA NTRS2012

Globally Gridded Satellite (GridSat) Observations for Climate Studies

Geostationary satellites have provided routine, high temporal resolution Earth observations since the 1970s. Despite the long period of record, use of these data in climate studies has been limited for numerous reasons, among them: there is no central archive of geostationary data for all international satellites, full temporal and spatial resolution data are voluminous, and diverse calibration and navigation formats encumber the uniform processing needed for multi-satellite climate studies. The International Satellite Cloud Climatology Project set the stage for overcoming these issues by archiving a subset of the full resolution geostationary data at approx.10 km resolution at 3 hourly intervals since 1983. Recent efforts at NOAA s National Climatic Data Center to provide convenient access to these data include remapping the data to a standard map projection, recalibrating the data to optimize temporal homogeneity, extending the record of observations back to 1980, and reformatting the data for broad public distribution. The Gridded Satellite (GridSat) dataset includes observations from the visible, infrared window, and infrared water vapor channels. Data are stored in the netCDF format using standards that permit a wide variety of tools and libraries to quickly and easily process the data. A novel data layering approach, together with appropriate satellite and file metadata, allows users to access GridSat data at varying levels of complexity based on their needs. The result is a climate data record already in use by the meteorological community. Examples include reanalysis of tropical cyclones, studies of global precipitation, and detection and tracking of the intertropical convergence zone.

Knapp, Kenneth R.↗
NASA NTRSDate not supplied

WMLES for the Fifth High-Lift PredictionWorkshop Cases Using FUN3D

This paper presents solution assessments and grid convergence studies for the test cases outlined in the Fifth High-Lift Prediction Workshop (HLPW-5), focusing on the high-lift Common Research Models (CRM-HL). The study utilizes a wall-modeled large-eddy simulation (WMLES) methodology developed in the unstructured-grid, node-centered flow solver FUN3D. The second-order accurate simulations conducted in this study utilize a finite-volume spatial discretization and an implicit temporal scheme. Large-scale turbulent features are resolved away from the wall, with small-scale effects captured by the Vreman subgrid-scale model. An equilibrium wall function uses the first grid point off the wall serving as the critical interface between the wall model and the large-eddy simulation region, thus requiring careful placement in grid design. WMLES solutions are assessed for HLPW-5 cases, including a clean wing-body configuration and geometry-buildup configurations corresponding to the 5.1\% ONERA CRM-HL model. Grid-convergence studies are systematically conducted using uniformly refined grids. Moreover, simulation results and grid sensitivity are presented for the NASA 5.2\% CRM-HL configuration at both moderate and flight-scale Reynolds numbers. Overall, WMLES results are satisfactory and agree well with available experimental data, especially on sufficiently fine grids.

high-lift aerodynamics↗
NASA NTRS1993

Incorporation of a hierarchical grid component structure into GRIDGEN

The underlying framework of the GRIDGEN multiple block grid generation system has been refined so that grid components are now stored within a hierarchical data structure. This restructuring has enhanced the usability of the software by allowing grids to be generated on a more intuitive level. This new framework also provides a means by which the multiple block system can be edited at most any level in the grid generation process. Editing tools are currently being added to GRIDGEN so that a change to the grid can be propagated backward and forward in the data hierarchy. The new data structure, the editing tools, and other recent GRIDGEN improvements are described in this paper.

Steinbrenner, John P.↗
NASA NTRS1995

Managing Complexity in Multidisciplinary Visualization

As high performance computing technology progresses, computational simulations are becoming more advanced in their capabilities. In the computational aerosciences domain, single discipline steady-state simulations computed on a single grid are far from the state-of-the-art. In their place are complex, time-dependent multidisciplinary simulations that attempt to model a given geometry more realistically. The product of these multidisciplinary simulations is a massive amount of data stored in different formats, grid topologies, units of measure, etc., as a result of the differences in the simulated physical domains. In addition to the challenges posed by setting up and performing the simulation, additional challenges exist in analyzing computational results. Visualization plays an important role in the advancement of multidisciplinary simulations. To date, visualization has been used to aid in the interpretation of large amounts of simulation data. Because the human visual system is effective in digesting a large amount of information presented graphically, visualization has helped simulation scientists to understand complex simulation results. As these simulations become even more complex, integrating several different physical domains, visualization will be critical to digest the massive amount of information. Another important role for visualization is to provide a common communication medium from which the domain scientists can use to develop, debug, and analyze their work. Multidisciplinary analyses are the next step in simulation technology, not only in computational aerosciences, but in many other areas such as global climate modeling. Visualization researchers must understand and work towards the challenges posed by multidisciplinary simulation scenarios. This paper addresses some of these challenges, describing technologies that must be investigated to create a useful visualization analysis tool for domain scientists.

Miceli, Kristina D.↗
NASA NTRS2004

Grid Enabled Geospatial Catalogue Web Service

Geospatial Catalogue Web Service is a vital service for sharing and interoperating volumes of distributed heterogeneous geospatial resources, such as data, services, applications, and their replicas over the web. Based on the Grid technology and the Open Geospatial Consortium (0GC) s Catalogue Service - Web Information Model, this paper proposes a new information model for Geospatial Catalogue Web Service, named as GCWS which can securely provides Grid-based publishing, managing and querying geospatial data and services, and the transparent access to the replica data and related services under the Grid environment. This information model integrates the information model of the Grid Replica Location Service (RLS)/Monitoring & Discovery Service (MDS) with the information model of OGC Catalogue Service (CSW), and refers to the geospatial data metadata standards from IS0 19115, FGDC and NASA EOS Core System and service metadata standards from IS0 191 19 to extend itself for expressing geospatial resources. Using GCWS, any valid geospatial user, who belongs to an authorized Virtual Organization (VO), can securely publish and manage geospatial resources, especially query on-demand data in the virtual community and get back it through the data-related services which provide functions such as subsetting, reformatting, reprojection etc. This work facilitates the geospatial resources sharing and interoperating under the Grid environment, and implements geospatial resources Grid enabled and Grid technologies geospatial enabled. It 2!so makes researcher to focus on science, 2nd not cn issues with computing ability, data locztic~, processir,g and management. GCWS also is a key component for workflow-based virtual geospatial data producing.

Chen, Ai-Jun↗
NASA NTRS1999

Smagglce: Surface Modeling and Grid Generation for Iced Airfoils: Phase 1 Results

SmaggIce (Surface Modeling and Grid Generation for Iced Airfoils) is a software toolkit used in the process of aerodynamic performance prediction of iced airfoils with grid-based Computational Fluid Dynamics (CFD). It includes tools for data probing, boundary smoothing, domain decomposition, and structured grid generation and refinement. SmaggIce provides the underlying computations to perform these functions, a GUI (Graphical User Interface) to control and interact with those functions, and graphical displays of results, it is being developed at NASA Glenn Research Center. This paper discusses the overall design of SmaggIce as well as what has been implemented in Phase 1. Phase 1 results provide two types of software tools: interactive ice shape probing and interactive ice shape control. The ice shape probing tools will provide aircraft icing engineers and scientists with an interactive means to measure the physical characteristics of ice shapes. On the other hand, the ice shape control features of SmaggIce will allow engineers to examine input geometry data, correct or modify any deficiencies in the geometry, and perform controlled systematic smoothing to a level that will make the CFD process manageable.

Vickerman, Mary B.↗
NASA NTRS1992

Monolithic millimeter-wave diode array beam controllers: Theory and experiment

In the current work, multi-function beam control arrays have been fabricated and have successfully demonstrated amplitude control of transmitted beams in the W and D bands (75-170 GHz). While these arrays are designed to provide beam control under DC bias operation, new designs for high-speed electronic and optical control are under development. These arrays will fill a need for high-speed watt-level beam switches in pulsed reflectometer systems under development for magnetic fusion plasma diagnostics. A second experimental accomplishment of the current work is the demonstration in the 100-170 GHz (D band) frequency range of a new technique for the measurement of the transmission phase as well as amplitude. Transmission data can serve as a means to extract ('de-embed') the grid parameters; phase information provides more complete data to assist in this process. Additional functions of the array beam controller yet to be tested include electronically controlled steering and focusing of a reflected beam. These have application in the areas of millimeter-wave electronic scanning radar and reflectometry, respectively.

Sjogren, L. B.↗
NASA NTRS2006

Convective Weather Forecast Quality Metrics for Air Traffic Management Decision-Making

Since numerical weather prediction models are unable to accurately forecast the severity and the location of the storm cells several hours into the future when compared with observation data, there has been a growing interest in probabilistic description of convective weather. The classical approach for generating uncertainty bounds consists of integrating the state equations and covariance propagation equations forward in time. This step is readily recognized as the process update step of the Kalman Filter algorithm. The second well known method, known as the Monte Carlo method, consists of generating output samples by driving the forecast algorithm with input samples selected from distributions. The statistical properties of the distributions of the output samples are then used for defining the uncertainty bounds of the output variables. This method is computationally expensive for a complex model compared to the covariance propagation method. The main advantage of the Monte Carlo method is that a complex non-linear model can be easily handled. Recently, a few different methods for probabilistic forecasting have appeared in the literature. A method for computing probability of convection in a region using forecast data is described in Ref. 5. Probability at a grid location is computed as the fraction of grid points, within a box of specified dimensions around the grid location, with forecast convection precipitation exceeding a specified threshold. The main limitation of this method is that the results are dependent on the chosen dimensions of the box. The examples presented Ref. 5 show that this process is equivalent to low-pass filtering of the forecast data with a finite support spatial filter. References 6 and 7 describe the technique for computing percentage coverage within a 92 x 92 square-kilometer box and assigning the value to the center 4 x 4 square-kilometer box. This technique is same as that described in Ref. 5. Characterizing the forecast, following the process described in Refs. 5 through 7, in terms of percentage coverage or confidence level is notionally sound compared to characterizing in terms of probabilities because the probability of the forecast being correct can only be determined using actual observations. References 5 through 7 only use the forecast data and not the observations. The method for computing the probability of detection, false alarm ratio and several forecast quality metrics (Skill Scores) using both the forecast and observation data are given in Ref. 2. This paper extends the statistical verification method in Ref. 2 to determine co-occurrence probabilities. The method consists of computing the probability that a severe weather cell (grid location) is detected in the observation data in the neighborhood of the severe weather cell in the forecast data. Probabilities of occurrence at the grid location and in its neighborhood with higher severity, and with lower severity in the observation data compared to that in the forecast data are examined. The method proposed in Refs. 5 through 7 is used for computing the probability that a certain number of cells in the neighborhood of severe weather cells in the forecast data are seen as severe weather cells in the observation data. Finally, the probability of existence of gaps in the observation data in the neighborhood of severe weather cells in forecast data is computed. Gaps are defined as openings between severe weather cells through which an aircraft can safely fly to its intended destination. The rest of the paper is organized as follows. Section II summarizes the statistical verification method described in Ref. 2. The extension of this method for computing the co-occurrence probabilities in discussed in Section HI. Numerical examples using NCWF forecast data and NCWD observation data are presented in Section III to elucidate the characteristics of the co-occurrence probabilities. This section also discusses the procedure for computing throbabilities that the severity of convection in the observation data will be higher or lower in the neighborhood of grid locations compared to that indicated at the grid locations in the forecast data. The probability of coverage of neighborhood grid cells is also described via examples in this section. Section IV discusses the gap detection algorithm and presents a numerical example to illustrate the method. The locations of the detected gaps in the observation data are used along with the locations of convective weather cells in the forecast data to determine the probability of existence of gaps in the neighborhood of these cells. Finally, the paper is concluded in Section V.

Chatterji, Gano B.↗
NASA NTRS2002

Summary of Data from the First AIAA CFD Drag Prediction Workshop

The results from the first AIAA CFD Drag Prediction Workshop are summarized. The workshop was designed specifically to assess the state-of-the-art of computational fluid dynamics methods for force and moment prediction. An impartial forum was provided to evaluate the effectiveness of existing computer codes and modeling techniques, and to identify areas needing additional research and development. The subject of the study was the DLR-F4 wing-body configuration, which is representative of transport aircraft designed for transonic flight. Specific test cases were required so that valid comparisons could be made. Optional test cases included constant-C(sub L) drag-rise predictions typically used in airplane design by industry. Results are compared to experimental data from three wind tunnel tests. A total of 18 international participants using 14 different codes submitted data to the workshop. No particular grid type or turbulence model was more accurate, when compared to each other, or to wind tunnel data. Most of the results overpredicted C(sub Lo) and C(sub Do), but induced drag (dC(sub D)/dC(sub L)(exp 2)) agreed fairly well. Drag rise at high Mach number was underpredicted, however, especially at high C(sub L). On average, the drag data were fairly accurate, but the scatter was greater than desired. The results show that well-validated Reynolds-Averaged Navier-Stokes CFD methods are sufficiently accurate to make design decisions based on predicted drag.

Levy, David W.↗
NASA NTRS2005

Data Format Standardization of Space Weather Model Output at the Community Coordinated Modeling Center

The disparate nature of space weather model output provides many challenges with regards to the portability and reuse of not only the data itself, but also any tools that are developed for analysis and visualization. We are developing and implementing a comprehensive data format standardization methodology that allows heterogeneous model output data to be stored uniformly in any common science data format. We will discuss our approach to identifying core meta-data elements that can be used to supplement raw model output data, thus creating self-descriptive files. The meta-data should also contain information describing the simulation grid. This will ultimately assists in the development of efficient data access tools capable of extracting data at any given point and time. We will also discuss our experiences standardizing the output of two global magnetospheric models, and how we plan to apply similar procedures when standardizing the output of the solar, heliospheric, and ionospheric models that are also currently hosted at the Community Coordinated Modeling Center.

Maddox, M.↗
NASA NTRS1999

Local Data Integration in East Central Florida

The Applied Meteorology Unit has configured a Local Data Integration System (LDIS) for east central Florida which assimilates in-situ and remotely-sensed observational data into a series of high-resolution gridded analyses. The ultimate goal for running LDIS is to generate products that may enhance weather nowcasts and short-range (less than 6 h) forecasts issued in support of the 45th Weather Squadron (45 WS), Spaceflight Meteorology Group (SMG), and the Melbourne National Weather Service (NWS MLB) operational requirements. LDIS has the potential to provide added value for nowcasts and short-ten-n forecasts for two reasons. First, it incorporates all data operationally available in east central Florida. Second, it is run at finer spatial and temporal resolutions than current national-scale operational models such as the Rapid Update Cycle and Eta models. LDIS combines all available data to produce grid analyses of primary variables (wind, temperature, etc.) at specified temporal and spatial resolutions. These analyses of primary variables can be used to compute diagnostic quantities such as vorticity and divergence. This paper demonstrates the utility of LDIS over east central Florida for a warm season case study. The evolution of a significant thunderstorm outflow boundary is depicted through horizontal and vertical cross section plots of wind speed, divergence, and circulation. In combination with a suitable visualization too], LDIS may provide users with a more complete and comprehensive understanding of evolving mesoscale weather than could be developed by individually examining the disparate data sets over the same area and time.

Case, Jonathan L.↗
NASA NTRS2018

Computational Investigation of Nominally-Orthogonal Pneumatic Active Flow Control for High-Lift Systems

We explore the feasibility of using nominally-orthogonal jets as active aerodynamic load control for multi-element high-lift systems, and whether the nominally-orthogonal jets can offer a variety of performance improvements. These nominally-orthogonal jets inject momentum normal to the airfoil surface near the flap trailing edge, where they create a vortex that entrains flow from the opposing side and change the airfoil circulation. Lift-enhancement opportunities of trailing edge nominally-orthogonal jets have previously been studied by Malavard et al. and Blaylock et al. on single-element airfoils; however, their effect on drag was not thoroughly investigated. In this study, we investigate two- dimensional nominally-orthogonal jet effects on both lift and drag on a two-element airfoil, NLR7301. We utilize Chimera Grid Tools to generate structured curvilinear overset grids, and the Reynolds-averaged Navier-Stokes solver OVERFLOW-2 to solve for the flow field around the airfoil. We perform various computational sensitivity studies on the baseline airfoil without a jet to validate computational results against benchmark experimental data. Using a Chimera overset grid topology, we demonstrate a similar lift-enhancement effect between a nominally-orthogonal jet and a nominally-orthogonal physical tab employed at the same location on the studied airfoil. After we introduce the nominally-orthogonal jet concept, we investigate nominally-orthogonal jets with various momentum coefficient settings, C(sub µ) = 0.00 − 0.04 and present a lift-enhancement relationship ∆C(sub l) ≃ 3.59 (C(sub µ)) for this airfoil. We discuss that utilizing a nominally-orthogonal jet with C(sub µ) = 0.01 can shift the linear region of the lift curve by a ∆C(sub l) = 0.36 for the pressure side jet and by a ∆C(sub l) = −0.27 for the suction side jet. Employing a nominally-orthogonal jet is also shown effective in altering the drag. To study the impact, we carry out a drag decomposition study in the form of drag polars. We show for a given C(sub l) = 2.50, a nominally-orthogonal jet with C(sub µ) = 0.01 on the pressure and suction side of the airfoil results in 113 drag count decrements and 41 drag count increments, respectively, compared to the baseline airfoil with no jet. These results show that large and controllable changes in aerodynamic performance can be achieved by relatively small active flow control inputs using the nominally-orthogonal jets presented in this study.

Nominally-Orthogonal↗
NASA NTRS2022

Coupling Carbon Oxidation and Surface Recession in Direct-Simulation Monte Carlo Code, SPARTA

Ablative thermal protection system (TPS) materials for spacecraft are composites that are often made out of carbon-based reinforcement and a polymeric matrix. They endure high-temperature oxidation and surface recession when re-entering Earth’s atmosphere. Ablation is the result of many coupled and competing thermal, mechanical, and chemical phenomena, and it is difficult to isolate the role of each on the overall degradation of the TPS. Here we develop an ablation model for material recession coupled explicitly to finite rate carbon oxidation in complex microstructures. In this work, Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al, was adapted to tally surface reactions and other surface data on a per-grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also develop robust algorithms which handle the evolution of the flow cells and solid material regions, including split cells (flow cell divided in two by a solid surface). Finally, we demonstrate our implicit chemistry model for 2D and 3D geometries by producing reaction statistics and detailed visualization of oxidation-induced material recession at the microscale.

V Arias↗
NASA NTRS2000

Incremental Parallelization of Non-Data-Parallel Programs Using the Charon Message-Passing Library

Message passing is among the most popular techniques for parallelizing scientific programs on distributed-memory architectures. The reasons for its success are wide availability (MPI), efficiency, and full tuning control provided to the programmer. A major drawback, however, is that incremental parallelization, as offered by compiler directives, is not generally possible, because all data structures have to be changed throughout the program simultaneously. Charon remedies this situation through mappings between distributed and non-distributed data. It allows breaking up the parallelization into small steps, guaranteeing correctness at every stage. Several tools are available to help convert legacy codes into high-performance message-passing programs. They usually target data-parallel applications, whose loops carrying most of the work can be distributed among all processors without much dependency analysis. Others do a full dependency analysis and then convert the code virtually automatically. Even more toolkits are available that aid construction from scratch of message passing programs. None, however, allows piecemeal translation of codes with complex data dependencies (i.e. non-data-parallel programs) into message passing codes. The Charon library (available in both C and Fortran) provides incremental parallelization capabilities by linking legacy code arrays with distributed arrays. During the conversion process, non-distributed and distributed arrays exist side by side, and simple mapping functions allow the programmer to switch between the two in any location in the program. Charon also provides wrapper functions that leave the structure of the legacy code intact, but that allow execution on truly distributed data. Finally, the library provides a rich set of communication functions that support virtually all patterns of remote data demands in realistic structured grid scientific programs, including transposition, nearest-neighbor communication, pipelining, gather/scatter, and redistribution. At the end of the conversion process most intermediate Charon function calls will have been removed, the non-distributed arrays will have been deleted, and virtually the only remaining Charon functions calls are the high-level, highly optimized communications. Distribution of the data is under complete control of the programmer, although a wide range of useful distributions is easily available through predefined functions. A crucial aspect of the library is that it does not allocate space for distributed arrays, but accepts programmer-specified memory. This has two major consequences. First, codes parallelized using Charon do not suffer from encapsulation; user data is always directly accessible. This provides high efficiency, and also retains the possibility of using message passing directly for highly irregular communications. Second, non-distributed arrays can be interpreted as (trivial) distributions in the Charon sense, which allows them to be mapped to truly distributed arrays, and vice versa. This is the mechanism that enables incremental parallelization. In this paper we provide a brief introduction of the library and then focus on the actual steps in the parallelization process, using some representative examples from, among others, the NAS Parallel Benchmarks. We show how a complicated two-dimensional pipeline-the prototypical non-data-parallel algorithm- can be constructed with ease. To demonstrate the flexibility of the library, we give examples of the stepwise, efficient parallel implementation of nonlocal boundary conditions common in aircraft simulations, as well as the construction of the sequence of grids required for multigrid.

VanderWijngaart, Rob F.↗