Impact of Ada in the Flight Dynamics Division: Excitement and frustration

In 1985, NASA Goddard's Flight Dynamics Division (FDD) began investigating how the Ada language might apply to their software development projects. Although they began cautiously using Ada on only a few pilot projects, they expected that, if the Ada pilots showed promising results, they would fully transition their entire development organization from FORTRAN to Ada within 10 years. However, nearly 9 years later, the FDD still produces 80 percent of its software in FORTRAN, despite positive results on Ada projects. This paper reports preliminary results of an ongoing study, commissioned by the FDD, to quantify the impact of Ada in the FDD, to determine why Ada has not flourished, and to recommend future directions regarding Ada. Project trends in both languages are examined as are external factors and cultural issues that affected the infusion of this technology. This paper is the first public report on the Ada assessment study, which will conclude with a comprehensive final report in mid 1994

Possible climate change impacts on water resources availability in a large semi-arid catchment in Northeast Brazil.

The semiarid region of Northeast Brazil is characterized by water scarcity, vulnerability of natural resources, and pronounced climatic variability. An integrated model has been developed to simulate this complex situation with an emphasis on a large-scale representation of hydrological processes and on the sensitivity to climate change. Regional climate change scenarios were obtained by empirical downscaling with large-scale climate information from different GCMs which differ strongly in their projections for future precipitation. The results show that due to these differences, it is still impossible to give quantitative values of the water availability in a forecast sense, i.e. to assign probabilities to the simulated results. However, it becomes clear that efficient and ecologically sound water management is a key question for further development. The results show that, independent of the climate change, agriculture is more vulnerable to drought impacts in the case of rainfed compared to irrigated farming. However, the capacity of irrigation and water infrastructure to enhance resilience with respect to climatic fluctuations is significantly constrained in the case of a negative precipitation trend

Modeling sea level changes and geodetic variations by glacial isostasy: the improved SELEN code

We describe the basic features of SELEN, an open source Fortran 90 program for the numerical solution of the so-called "Sea Level Equation" for a spherical, layered, non-rotating Earth with Maxwell viscoelastic rheology. The Sea Level Equation was introduced in the 70s to model the sea level variations in response to the melting of late-Pleistocene ice-sheets, but it can be also employed for predictions of geodetic quantities such as vertical and horizontal surface displacements and gravity variations on a global and a regional scale. SELEN (acronym of SEa Level EquatioN solver) is particularly oriented to scientists at their first approach to the glacial isostatic adjustment problem and, according to our experience, it can be successfully used in teaching. The current release (2.9) considerably improves the previous versions of the code in terms of computational efficiency, portability and versatility. In this paper we describe the essentials of the theory behind the Sea Level Equation, the purposes of SELEN and its implementation, and we provide practical guidelines for the use of the program. Various examples showing how SELEN can be configured to solve geodynamical problems involving past and present sea level changes and current geodetic variations are also presented and discussed

Software development tools: A bibliography, appendix C.

A bibliography containing approximately 200 citations on tools which help software developers perform some development task (such as text manipulation, testing, etc.), and which would not necessarily be found as part of a computing facility is given. The bibliography comes from a relatively random sampling of the literature and is not complete. But it is indicative of the nature and range of tools currently being prepared or currently available

Data Driven Surrogate Based Optimization in the Problem Solving Environment WBCSim

Large scale, multidisciplinary, engineering designs are always difficult due to the complexity and dimensionality of these problems. Direct coupling between the analysis codes and the optimization routines can be prohibitively time consuming due to the complexity of the underlying simulation codes. One way of tackling this problem is by constructing computationally cheap(er) approximations of the expensive simulations, that mimic the behavior of the simulation model as closely as possible. This paper presents a data driven, surrogate based optimization algorithm that uses a trust region based sequential approximate optimization (SAO) framework and a statistical sampling approach based on design of experiment (DOE) arrays. The algorithm is implemented using techniques from two packages—SURFPACK and SHEPPACK that provide a collection of approximation algorithms to build the surrogates and three different DOE techniques—full factorial (FF), Latin hypercube sampling (LHS), and central composite design (CCD)—are used to train the surrogates. The results are compared with the optimization results obtained by directly coupling an optimizer with the simulation code. The biggest concern in using the SAO framework based on statistical sampling is the generation of the required database. As the number of design variables grows, the computational cost of generating the required database grows rapidly. A data driven approach is proposed to tackle this situation, where the trick is to run the expensive simulation if and only if a nearby data point does not exist in the cumulatively growing database. Over time the database matures and is enriched as more and more optimizations are performed. Results show that the proposed methodology dramatically reduces the total number of calls to the expensive simulation runs during the optimization process

A Simple Multiprocessor Management System for Event-Parallel Computing

Offline software using TCP/IP sockets to distribute particle physics events to multiple UNIX/RISC workstations is described. A modular, building block approach was taken, which allowed tailoring to solve specific tasks efficiently and simply as they arose. The modest, initial cost was having to learn about sockets for interprocess communication. This multiprocessor management software has been used to control the reconstruction of eight billion raw data events from Fermilab Experiment E791.Comment: 10 pages, 3 figures, compressed Postscript, LaTeX. Submitted to NI

Geothermal systems simulation: A case study

Geothermal reservoir simulation is a key step for developing sustainable and efficient strategies for the exploitation of geothermal resources. It is applied in the assessment of several areas of reservoir engineering, such as reservoir performance and re-injection programs, pressure decline in depletion, phase transition conditions, and natural evolution of hydrothermal convection systems. Fluid flow and heat transfer in rock masses, fluid-rock chemical interaction and rock mass deformation are some of the processes addressed in reservoir modelling. The case study of the Las Tres Virgenes (LTV) geothermal field (10 MWe), Baja California Sur, Mexico is presented. Three dimensional (3D) natural state simulations were carried out from emplacement and cooling of two spherical magma chambers using a conductive approach. A conceptual model of the volcanic system was developed on a lithostratigraphic and geochronological basis. Magma chamber volumes were established from eruptive volumes estimations. The thermophysical properties of the medium were assumed to correspond to the dominant rock in each lithological unit as an initial value, and further calibration was made considering histograms of experimentally obtained thermophysical properties of rocks. As the boundaries of the model lie far from the thermal anomaly, we assumed specified temperature boundaries. A Finite Volume (FV) numerical scheme was implemented in a Fortran 90 code to solve the heat equation. Static formation temperatures from well logs were used for validation of the numerical results. Good agreement was observed in those geothermal wells dominated by conductive heat transfer. For other wells, however, it is clear that conduction alone cannot explain observed behaviour, three-dimensional convective models are being implemented for future multiphysics simulations