Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple time stepping”

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 181 records · Page 10

Identifying Basic Factors for Communal Prosperity - Space Technologies are Bridging this Gap

There are many aspects, which are important for maintaining environmentally clean and safe conditions for a healthy and economically self-sufficient community. This problem was somewhat of a lesser concern in earlier days because many communities were small, isolated and solely dependent upon their owners or landlords. Due to an astronomical growth in human population within the last century, extensive use of combustion technologies, and changing environmental conditions has resulted in scarcity of natural resources. In reality, the societal sustainability issues are becoming much more acute and complex. Therefore, the researchers and social scientists are joining forces to address these topics and find solutions to many contentious areas such as public health and diseases, water resources, agriculture production, survivability during and after the natural disasters, energy needs and many others. Forthrightly speaking, there is no canned solution or a methodology to go about solving these issues since the magnitude and complexity of these issues are multi-dimensional and are further inter-locked with other areas. A common sense tells us that we need data, resources and technologies to begin addressing these problems. This is where space observations have provided us with tremendous information and opportunities, which are of great assets to the science, economist, and social scientists. This paper specifically addresses what are critical areas for a successful societal sustainability and growth; and how can we take advantage of multiple sensors and models already in existence. Increasing our knowledge of the home planet, via amplified set of observations, is certainly a right step in a right direction. Furthermore, this is a pre-requisite in understanding multiple hazard phenomena's. This paper further examines various space sensors and observing architectures that can be useful specifically in addressing some of these complex issues. The ultimate goal is to serve the society by providing such valuable information so the decision makers can take full advantage in making timely societal decisions.

Habib, Shahid↗

Affordability Approaches for Human Space Exploration

The design and development of historical NASA Programs (Apollo, Shuttle and International Space Station), have been based on pre-agreed missions which included specific pre-defined destinations (e.g., the Moon and low Earth orbit). Due to more constrained budget profiles, and the desire to have a more flexible architecture for Mission capture as it is affordable, NASA is working toward a set of Programs that are capability based, rather than mission and/or destination specific. This means designing for a performance capability that can be applied to a specific human exploration mission/destination later (sometime years later). This approach does support developing systems to flatter budgets over time, however, it also poses the challenge of how to accomplish this effectively while maintaining a trained workforce, extensive manufacturing, test and launch facilities, and ensuring mission success ranging from Low Earth Orbit to asteroid destinations. NASA Marshall Space Flight Center (MSFC) in support of Exploration Systems Directorate (ESD) in Washington, DC has been developing approaches to track affordability across multiple Programs. The first step is to ensure a common definition of affordability: the discipline to bear cost in meeting a budget with margin over the life of the program. The second step is to infuse responsibility and accountability for affordability into all levels of the implementing organization since affordability is no single person s job; it is everyone s job. The third step is to use existing data to identify common affordability elements organized by configuration (vehicle/facility), cost, schedule, and risk. The fourth step is to analyze and trend this affordability data using an affordability dashboard to provide status, measures, and trends for ESD and Program level of affordability tracking. This paper will provide examples of how regular application of this approach supports affordable and therefore sustainable human space exploration architecture.

Holladay, Jon↗

Genesis Solar Wind Collector Cleaning Assessment: Update on 60336 Sample Case Study

To maximize the scientific return of Genesis Solar Wind return mission it is necessary to characterize and remove a crash-derived particle and thin film surface contamination. A small subset of Genesis mission collector fragments are being subjected to extensive study via various techniques. Here we present an update on the sample 60336, a Czochralski silicon (Si-CZ) based wafer from the bulk array (B/C). This sample has undergone multiple cleaning steps (see the table below): UPW spin wash, aggressive chemical cleanings (including aqua regia, hot xylene and RCA1), as well as optical and chemical (EDS, ToF-SIMS) imaging. Contamination appeared on the surface of 60336 after the initial 2007 UPW cleaning. Aqua regia and hot xylene treatment (8/13/2013) did little to remove contaminants. The sample was UPW cleaned for the third time and imaged (9/16/13). The UPW removed the dark stains that were visible on the sample. However, some features, like "the Flounder" (a large, 100 micron feature in Fig. 1b) appeared largely intact, resisting all previous cleaning efforts. These features were likely from mobilized adhesive, derived from the Post-It notes used to stabilize samples for transport from Utah after the hard landing. To remove this contamination, an RCA step 1 organic cleaning (RCA1) was employed. Although we are still uncertain on the nature of the Flounder and why it is resistant to UPW and aqua regia/hot xylene treatment, we have found RCA1 to be suitable for its removal. It is likely that the glue from sticky pads used during collector recovery may have been a source for resistant organic contamination [9]; however [8] shows that UPW reaction with crash-derived organic contamination does not make particle removal more difficult.

Goreva, Y. S.↗

Time-History Statistics of Soot Formation in A Model Gas Turbine Combustor

Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.

LES↗

A Framework for Simulation of Aircraft Flyover Noise Through a Non-Standard Atmosphere

This paper describes a new framework for the simulation of aircraft flyover noise through a non-standard atmosphere. Central to the framework is a ray-tracing algorithm which defines multiple curved propagation paths, if the atmosphere allows, between the moving source and listener. Because each path has a different emission angle, synthesis of the sound at the source must be performed independently for each path. The time delay, spreading loss and absorption (ground and atmosphere) are integrated along each path, and applied to each synthesized aircraft noise source to simulate a flyover. A final step assigns each resulting signal to its corresponding receiver angle for the simulation of a flyover in a virtual reality environment. Spectrograms of the results from a straight path and a curved path modeling assumption are shown. When the aircraft is at close range, the straight path results are valid. Differences appear especially when the source is relatively far away at shallow elevation angles. These differences, however, are not significant in common sound metrics. While the framework used in this work performs off-line processing, it is conducive to real-time implementation.

Arntzen, Michael↗

Cooperative Three-Robot System for Traversing Steep Slopes

Teamed Robots for Exploration and Science in Steep Areas (TRESSA) is a system of three autonomous mobile robots that cooperate with each other to enable scientific exploration of steep terrain (slope angles up to 90 ). Originally intended for use in exploring steep slopes on Mars that are not accessible to lone wheeled robots (Mars Exploration Rovers), TRESSA and systems like TRESSA could also be used on Earth for performing rescues on steep slopes and for exploring steep slopes that are too remote or too dangerous to be explored by humans. TRESSA is modeled on safe human climbing of steep slopes, two key features of which are teamwork and safety tethers. Two of the autonomous robots, denoted Anchorbots, remain at the top of a slope; the third robot, denoted the Cliffbot, traverses the slope. The Cliffbot drives over the cliff edge supported by tethers, which are payed out from the Anchorbots (see figure). The Anchorbots autonomously control the tension in the tethers to counter the gravitational force on the Cliffbot. The tethers are payed out and reeled in as needed, keeping the body of the Cliffbot oriented approximately parallel to the local terrain surface and preventing wheel slip by controlling the speed of descent or ascent, thereby enabling the Cliffbot to drive freely up, down, or across the slope. Due to the interactive nature of the three-robot system, the robots must be very tightly coupled. To provide for this tight coupling, the TRESSA software architecture is built on a combination of (1) the multi-robot layered behavior-coordination architecture reported in "An Architecture for Controlling Multiple Robots" (NPO-30345), NASA Tech Briefs, Vol. 28, No. 10 (October 2004), page 65, and (2) the real-time control architecture reported in "Robot Electronics Architecture" (NPO-41784), NASA Tech Briefs, Vol. 32, No. 1 (January 2008), page 28. The combination architecture makes it possible to keep the three robots synchronized and coordinated, to use data from all three robots for decision- making at each step, and to control the physical connections among the robots. In addition, TRESSA (as in prior systems that have utilized this architecture) , incorporates a capability for deterministic response to unanticipated situations from yet another architecture reported in Control Architecture for Robotic Agent Command and Sensing (NPO-43635), NASA Tech Briefs, Vol. 32, No. 10 (October 2008), page 40. Tether tension control is a major consideration in the design and operation of TRESSA. Tension is measured by force sensors connected to each tether at the Cliffbot. The direction of the tension (both azimuth and elevation) is also measured. The tension controller combines a controller to counter gravitational force and an optional velocity controller that anticipates the motion of the Cliffbot. The gravity controller estimates the slope angle from the inclination of the tethers. This angle and the weight of the Cliffbot determine the total tension needed to counteract the weight of the Cliffbot. The total needed tension is broken into components for each Anchorbot. The difference between this needed tension and the tension measured at the Cliffbot constitutes an error signal that is provided to the gravity controller. The velocity controller computes the tether speed needed to produce the desired motion of the Cliffbot. Another major consideration in the design and operation of TRESSA is detection of faults. Each robot in the TRESSA system monitors its own performance and the performance of its teammates in order to detect any system faults and prevent unsafe conditions. At startup, communication links are tested and if any robot is not communicating, the system refuses to execute any motion commands. Prior to motion, the Anchorbots attempt to set tensions in the tethers at optimal levels for counteracting the weight of the Cliffbot; if either Anchorbot fails to reach its optimal tension level within a specified time, it sends message to the other robots and the commanded motion is not executed. If any mechanical error (e.g., stalling of a motor) is detected, the affected robot sends a message triggering stoppage of the current motion. Lastly, messages are passed among the robots at each time step (10 Hz) to share sensor information during operations. If messages from any robot cease for more than an allowable time interval, the other robots detect the communication loss and initiate stoppage.

Stroupe, Ashley↗

Terminal Sequencing and Spacing (TSS)

The Federal Aviation Administration's (FAA) Next Generation Air Transportation System (or NextGen) is being designed to support the predicted increases in traffic volume and to increase the capacity, efficiency and safety of the National Airspace System (NAS). The Federal Aviation Administration (FAA) identifies Performance-Based Navigation (PBN) as a key enabling capability of NextGen and is actively publishing PBN procedures at major airports throughout the United States. Standard Terminal Arrival Routes (STARs), procedures, and approaches are designed to facilitate fuel-efficient continuous descent operations. However, their use is limited during periods of high traffic demand due to the complexity of merging multiple streams of aircraft to the same airport. As a result, most arrivals in the Terminal Radar Approach Control (TRACON) area continue to be controlled using radar vectoring and step-down descents, resulting in high workload for controllers and diverting aircraft from efficient PBN trajectories. To address this issue, NASA developed the Terminal Sequencing and Spacing (TSS) system, an advanced arrival management technology that combines time-based scheduling and controller-based precision spacing tools. TSS is a ground-based controller automation tool that facilitates sequencing and merging arrivals on Performance-Based Navigation (PBN) routes, especially during highly congested demand periods. The two main components of TSS are: 1) a scheduler that de-conflicts merging arrivals in the terminal area by computing appropriate arrival times to the runway threshold and upstream terminal merge points, and 2) a set of Controller-Managed Spacing (CMS) decision support tools to efficiently assist schedule conformance. Sixteen high-fidelity human-in-the-loop simulations involving more than five hundred hours of evaluation time, were conducted to mature TSS from proof-of-concept design to a fully functional prototype. Results indicate high controller use and acceptability of the CMS tools as well as improved PBN route conformance (Figure 2). The TSS technology was transferred to the FAA in 2014, and it is targeted for deployment to several busy airports in the U.S. starting in 2018. Potential enhancements to TSS using DataComm will also be presented.

controller-managed spacing↗

An Overview of the Patch Integral Method (PIM), a New Heat Transfer Analysis Tool for Hypersonic Wind Tunnel Facilities at NASA Langley

NASA Langley’s hypersonic wind tunnels are heavily leveraged for planetary missions. The data collection method in these tunnels is thermography, and surface temperature measurements of the model surface are collected and reduced to produce surface heating data, as seen in Fig 1. However, during model injection, no temperature data are collected, and thus conventional, integral heat transfer methods cannot be used to solve for surface heating. A method was developed in the 1990’s to reduce this data despite the data gap, known as the step approximation method. The method assumes that the film coefficient behaves as a step function, the model is semi-infinite, and thermal properties are constant. With these simplifying assumptions, a Laplace transform can be performed to result in an equation that takes the initial temperature of the model and a temperature at some point in time to back out the film coefficient at that time. This is the method that is used in the current thermographic data reduction software, IHEAT. While computationally light-weight, the step approximation has several issues associated with it. The time-history of temperature is not accounted for, which is vital as heat transfer is an integral process. Additionally, the required semi-infinite assumption is unnecessary and might be violated during runtime. Thermal variation of material properties can have a sizeable impact on heating results and are not modeled by the method. This method also takes multiple seconds to “collapse” to a steady state, which is undesirable from both a facility and data reduction standpoint. The method is also very sensitive to the “effective time” approximation, an approximation of when heating instantaneously starts (which is a nonphysical simplification), and a small variation in this value can result in an error in heating results.

J. S. Cheatwood↗

Digital Earth Initiative: A Joint Interagency Program

The Digital Earth is a virtual representation of our planet that enables a person to explore and interact with the vast amounts of natural and cultural information gathered about the Earth. The Digital Earth comprises data interfaces and standards enabling access to geo-referenced data from remote sensing, cartographic, demographic, medical, and other sources to respond to questions posed by the user. In a recent address at the California Science Center in Los Angeles, Vice President Al Gore articulated a Digital Earth Vision. That vision spoke to developing a multi-resolution, three-dimensional representation of the planet, into which we can roam and zoom into vast quantities of embedded geo-referenced data. The vision was not limited to moving through space but also allowing travel over a time-line, which can be set for days, years, centuries, or even geological epochs. As prototypes become available, it would also be possible to interact with the Digital Earth in multiple places around the country with access to high-speed networks and at a more limited level of access over the Internet. NASA was asked by the Vice President to lead an interagency initiative that would take steps to bring this vision to the public. This talk describes the start-up and plans of the Digital Earth Interagency Working Group in the formulation of its charter, an architecture reference model for Digital Earth, public/private partnerships, cooperative agreement notices, Digital Earth prototypes, and testbeds. Animations employing technologies for virtual roaming and zooming through multi-resolution satellite data set as prototype systems will be presented along with examples of potential user scenarios. Plans for engaging academia and industry in implementing the Digital Earth initiative will be discussed.

Halem, Milton↗

Air Traffic Management Technology Demostration-1 (ATD-1) Avionics Phase 2 Flight Test and Results

NASA’s first Air Traffic Management Technology Demonstration (ATD-1) subproject successfully completed 19 days of flight test validation in January and February 2017 of an Interval Management (IM) avionics prototype and the procedures used to conduct IM arrival and approach operations. IM is one of the three elements integrated into NASA’s ATD-1 concept of operations with the subproject goal of improving aircraft efficiency and airport throughput during high-density arrival operations. The ATD-1 concept of operations combines advanced arrival scheduling, controller decision support tools, and interval management (IM) avionics to enable merging of multiple, time-based, efficient arrival streams. IM contributes to the operation by calculating speeds that enable an aircraft to precisely achieve a specific time or distance behind another aircraft. When precise spacing intervals can be calculated, achieved, and then maintained during high-density operations, aircraft efficiency should be improved by enabling the aircraft to remain closer to the optimum descent trajectory instead of using vectors and step-down altitudes, and airport throughput should be maintained or improved by each aircraft arriving at the runway threshold closer to the assigned spacing interval. This avionics development and flight test was conducted under a NASA contract by Boeing Research and Technology, with Boeing Commercial Aircraft, Honeywell, United Airlines, and Jeppesen as sub-contractors. The Honeywell built IM avionics were the first ever prototype built based on NASA requirements as well as developing and non-flight tested international IM standards, integrated into two test aircraft, and then flown in real-world conditions at the Grant County International Airport (KMWH). The IM prototype flown in the flight test used data from the Ownship and the assigned lead, or Target, aircraft to calculate the airspeed necessary for the Ownship to achieve the desired spacing. The flight test demonstrated that the IM avionics prototype generally met the IM requirement for spacing accuracy. However, the control laws implemented require further development to reduce the high IM speed command rate and the number of speed reversals observed during the test. Pilots assessed the IM procedure as acceptable, and issues requiring further attention were identified. In summary, the IM avionics prototype showed significant promise in contributing to the goals of improving aircraft efficiency and airport throughput. The flight test results also provided important data to the FAA and the working group developing the follow-on version of the international IM standards.

Baxley, Brian T.↗

Development of a three dimensional numerical water quality model for continental shelf applications

A model to predict the distribution of water quality parameters in three dimensions was developed. The mass transport equation was solved using a non-dimensional vertical axis and an alternating-direction-implicit finite difference technique. The reaction kinetics of the constituents were incorporated into a matrix method which permits computation of the interactions of multiple constituents. Methods for the computation of dispersion coefficients and coliform bacteria decay rates were determined. Numerical investigations of dispersive and dissipative effects showed that the three-dimensional model performs as predicted by analysis of simpler cases. The model was then applied to a two dimensional vertically averaged tidal dynamics model for the Providence River. It was also extended to a steady state application by replacing the time step with an iteration sequence. This modification was verified by comparison to analytical solutions and applied to a river confluence situation.

Spaulding, M.↗

A simple hydrologically based model of land surface water and energy fluxes for general circulation models

A generalization of the single soil layer variable infiltration capacity (VIC) land surface hydrological model previously implemented in the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation model (GCM) is described. The new model is comprised of a two-layer characterization of the soil column, and uses an aerodynamic representation of the latent and sensible heat fluxes at the land surface. The infiltration algorithm for the upper layer is essentially the same as for the single layer VIC model, while the lower layer drainage formulation is of the form previously implemented in the Max-Planck-Institut GCM. The model partitions the area of interest (e.g., grid cell) into multiple land surface cover types; for each land cover type the fraction of roots in the upper and lower zone is specified. Evapotranspiration consists of three components: canopy evaporation, evaporation from bare soils, and transpiration, which is represented using a canopy and architectural resistance formulation. Once the latent heat flux has been computed, the surface energy balance is iterated to solve for the land surface temperature at each time step. The model was tested using long-term hydrologic and climatological data for Kings Creek, Kansas to estimate and validate the hydrological parameters, and surface flux data from three First International Satellite Land Surface Climatology Project Field Experiment (FIFE) intensive field campaigns in the summer-fall of 1987 to validate the surface energy fluxes.

Liang, XU↗

Implementation and Characterization of Three-Dimensional Particle-in-Cell Codes on Multiple-Instruction-Multiple-Data Massively Parallel Supercomputers

A three-dimensional electrostatic particle-in-cell (PIC) plasma simulation code has been developed on coarse-grain distributed-memory massively parallel computers with message passing communications. Our implementation is the generalization to three-dimensions of the general concurrent particle-in-cell (GCPIC) algorithm. In the GCPIC algorithm, the particle computation is divided among the processors using a domain decomposition of the simulation domain. In a three-dimensional simulation, the domain can be partitioned into one-, two-, or three-dimensional subdomains ("slabs," "rods," or "cubes") and we investigate the efficiency of the parallel implementation of the push for all three choices. The present implementation runs on the Intel Touchstone Delta machine at Caltech; a multiple-instruction-multiple-data (MIMD) parallel computer with 512 nodes. We find that the parallel efficiency of the push is very high, with the ratio of communication to computation time in the range 0.3%-10.0%. The highest efficiency (> 99%) occurs for a large, scaled problem with 64(sup 3) particles per processing node (approximately 134 million particles of 512 nodes) which has a push time of about 250 ns per particle per time step. We have also developed expressions for the timing of the code which are a function of both code parameters (number of grid points, particles, etc.) and machine-dependent parameters (effective FLOP rate, and the effective interprocessor bandwidths for the communication of particles and grid points). These expressions can be used to estimate the performance of scaled problems--including those with inhomogeneous plasmas--to other parallel machines once the machine-dependent parameters are known.

Lyster, P. M.↗

Cyclotron Line Measurements with INTEGRAL

Due to its broadband energy coverage, INTEGRAL has made important contributions to observing and interpreting cyclotron lines, which are present in the 10-100 keV range of a sample of accreting pulsars. In these systems photons with energies fulfilling the resonance condition inelastically Compton scatter off electrons quantized in the accretion column above the neutron star's magnetic pole(s). This process gives rise to the broad, absorption-like lines or 'cyclotron resonant scattering features' (CRSF). The observed lines allow to directly measure the B-fields of these sources, resulting in values of a few times 1E12G. In this overview I will present recent highlights regarding CRSF observations as well as discuss current ideas and models for the physical conditions in the accretion column. Among the former are the stability of the spectrum of Vela X-1 during giant flares in 2003, the observation of three cyclotron lines during the 2004 outburst of V0332+53, the confirmation of the fundamental line at approximately 45 keV during a 2005 normal outburst of A0535-26, and the simultaneous detection of the two lines in the dipping source 4U 1907+09 (for which also a torque reversal was detected for the first time). Through these and other observations it has become increasingly apparent that two types of observations can potentially be used to constrain the accretion column geometry: the determination of energy ratios for multiple harmonic lines (only two sources with greater than 2 lines are known), was well as the evolution of the fundamental line centroid, which, for different sources, may or may not be correlated with flux. Furthermore, first steps have been taken away from the usual phenomenological description of the lines, towards a physical approach based on self-consistent CRSF modeling. Initial applications are presented.

Pottschmidt, K.↗

Recent Enhancements To The FUN3D Flow Solver For Moving-Mesh Applications

An unsteady Reynolds-averaged Navier-Stokes solver for unstructured grids has been extended to handle general mesh movement involving rigid, deforming, and overset meshes. Mesh deformation is achieved through analogy to elastic media by solving the linear elasticity equations. A general method for specifying the motion of moving bodies within the mesh has been implemented that allows for inherited motion through parent-child relationships, enabling simulations involving multiple moving bodies. Several example calculations are shown to illustrate the range of potential applications. For problems in which an isolated body is rotating with a fixed rate, a noninertial reference-frame formulation is available. An example calculation for a tilt-wing rotor is used to demonstrate that the time-dependent moving grid and noninertial formulations produce the same results in the limit of zero time-step size.

Biedron, Robert T,↗

Process for Patterning Indium for Bump Bonding

An innovation was created for the Cosmology Large Angular Scale Surveyor for integration of low-temperature detector chips with a silicon backshort and a silicon photonic choke through flipchip bonding. Indium bumps are typically patterned using liftoff processes, which require thick resist. In some applications, it is necessary to locate the bumps close to high-aspect-ratio structures such as wafer through-holes. In those cases, liftoff processes are challenging, and require complicated and time-consuming spray coating technology if the high-aspect-ratio structures are delineated prior to the indium bump process. Alternatively, processing the indium bumps first is limited by compatibility of the indium with subsequent processing. The present invention allows for locating bumps arbitrarily close to multiple-level high-aspect-ratio structures, and for indium bumps to be formed without liftoff resist. The process uses the poor step coverage of indium deposited on a silicon wafer that has been previously etched to delineate the location of the indium bumps. The silicon pattern can be processed through standard lithography prior to adding the high-aspect-ratio structures. Typically, high-aspectratio structures require a thick resist layer so this layer can easily cover the silicon topography. For multiple levels of topography, the silicon can be easily conformally coated through standard processes. A blanket layer of indium is then deposited onto the full wafer; bump bonding only occurs at the high points of the topography.

Denis, Kevin↗

Multiple steady states for characteristic initial value problems

The time dependent, isentropic, quasi-one-dimensional equations of gas dynamics and other model equations are considered under the constraint of characteristic boundary conditions. Analysis of the time evolution shows how different initial data may lead to different steady states and how seemingly anamolous behavior of the solution may be resolved. Numerical experimentation using time consistent explicit algorithms verifies the conclusions of the analysis. The use of implicit schemes with very large time steps leads to erroneous results.

Salas, M. D.↗

Performance Characteristics of the Multi-Zone NAS Parallel Benchmarks

We describe a new suite of computational benchmarks that models applications featuring multiple levels of parallelism. Such parallelism is often available in realistic flow computations on systems of grids, but had not previously been captured in bench-marks. The new suite, named NPB Multi-Zone, is extended from the NAS Parallel Benchmarks suite, and involves solving the application benchmarks LU, BT and SP on collections of loosely coupled discretization meshes. The solutions on the meshes are updated independently, but after each time step they exchange boundary value information. This strategy provides relatively easily exploitable coarse-grain parallelism between meshes. Three reference implementations are available: one serial, one hybrid using the Message Passing Interface (MPI) and OpenMP, and another hybrid using a shared memory multi-level programming model (SMP+OpenMP). We examine the effectiveness of hybrid parallelization paradigms in these implementations on three different parallel computers. We also use an empirical formula to investigate the performance characteristics of the multi-zone benchmarks.

Jin, Haoqiang↗