Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “global refinement”

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 415 records · Page 23

Incremental triangulation by way of edge swapping and local optimization

This document is intended to serve as an installation, usage, and basic theory guide for the two dimensional triangulation software 'HARLEY' written for the Silicon Graphics IRIS workstation. This code consists of an incremental triangulation algorithm based on point insertion and local edge swapping. Using this basic strategy, several types of triangulations can be produced depending on user selected options. For example, local edge swapping criteria can be chosen which minimizes the maximum interior angle (a MinMax triangulation) or which maximizes the minimum interior angle (a MaxMin or Delaunay triangulation). It should be noted that the MinMax triangulation is generally only locally optical (not globally optimal) in this measure. The MaxMin triangulation, however, is both locally and globally optical. In addition, Steiner triangulations can be constructed by inserting new sites at triangle circumcenters followed by edge swapping based on the MaxMin criteria. Incremental insertion of sites also provides flexibility in choosing cell refinement criteria. A dynamic heap structure has been implemented in the code so that once a refinement measure is specified (i.e., maximum aspect ratio or some measure of a solution gradient for the solution adaptive grid generation) the cell with the largest value of this measure is continually removed from the top of the heap and refined. The heap refinement strategy allows the user to specify either the number of cells desired or refine the mesh until all cell refinement measures satisfy a user specified tolerance level. Since the dynamic heap structure is constantly updated, the algorithm always refines the particular cell in the mesh with the largest refinement criteria value. The code allows the user to: triangulate a cloud of prespecified points (sites), triangulate a set of prespecified interior points constrained by prespecified boundary curve(s), Steiner triangulate the interior/exterior of prespecified boundary curve(s), refine existing triangulations based on solution error measures, and partition meshes based on the Cuthill-McKee, spectral, and coordinate bisection strategies.

Wiltberger, N. Lyn↗

Global Radiative Flux Profile Dataset: Revised and Extended

The third generation of the radiative flux profile data product, called ISCCP-FH, is described. The revisions over the previous generation (called ISCCP-FD) include improvements in the radiative model representation of gaseous and aerosol effects, as well as a refined statistical model of cloud vertical layer variations with cloud types, and increased spatial resolution. The new product benefits from the changes in the new H-version of the ISCCP cloud products (called ISCCP-H): higher spatial resolution, revised radiance calibration and treatment of ice clouds, treatment of aerosol effects, and revision of all the ancillary atmosphere and surface property products. The ISCCP-FH product is evaluated against more direct measurements from the Clouds and the Earth’s Radiant Energy System and the Baseline Surface Radiation Network products, showing some small, overall reductions in average flux uncertainties; but the main results are similar to ISCCP-FD: the ISCCP-FH uncertainties remain ≲10 Wm −2 at the top-of-atmosphere (TOA) and ≲15 Wm −2 at surface for monthly, regional averages. The long-term variations of TOA, surface and in-atmosphere net fluxes are documented and the possible transient cloud feedback implications of a long-term change of clouds are investigated. The cloud and flux variations from 1998 to 2012 suggest a positive cloud-radiative feedback on the oceanic circulation and a negative feedback on the atmospheric circulation. This example demonstrates that the ISCCP-FH product can provide useful diagnostic information about weather-to-interannual scale variations of radiation induced by changes in cloudiness as well as atmospheric and surface properties.

ISCCP-FH↗

U.S. GPM Status

The joint U.S.-Japan Global Precipitation Measurement (GPM) mission, now in its tenth year of operations, continues to pursue research on scientific and operational shortcomings in algorithms for retrieving global precipitation from satellite observations. A summary of operations, including the orbit boost, processing, science, and applications is provided. This includes key improvements to the GPM products for which NASA has lead responsibility. This includes refinements to ensure continuity for each product across the boundary in 2014 between the predecessor Tropical Rainfall Measuring Mission (TRMM) and the GPM Core Observatory. Early evaluations of the U.S. Science Team’s Integrated Multi-satellitE Retrievals for GPM (IMERG) in V07 show interesting behavior across the TRMM orbit boost that is highly relevant to the planned orbit boost for the GPM Core observatory. We will present summary statistics that illustrate the tremendous demand that exists in the user community for GPM precipitation products, which on the U.S. side is primarily focused on IMERG, and end by considering major issues that might be addressed in up-coming versions. These include the use of machine learning algorithms, the operational challenge of swarms of “small”, perhaps short-lived satellites, and continuation of precessing calibration satellites after GPM.

GPM↗

Pathways of bio-jet adoption in the US aviation industry with implications for the overall transportation and energy sectors: an integrated, multi-sectoral analysis of future scenarios

Bio-jet adoption has emerged as an attractive option to complement and supplement the use of refined fossil liquid fuels in the aviation industry in the US. However, there are significant uncertainties surrounding the costs of bio-jet including but not limited to costs of feedstock, transformation costs and the competition with co-products of bio-jet that may be demanded elsewhere in the transportation or energy sectors. This study models alternative trajectories of bio-jet adoption in the US aviation industry by 2050 through the use of a global integrated multi sector dynamics model. Three bio-jet production and consumption pathways are presented- soybean oil to jet, corn ethanol to jet (ETJ) and Fisher–Tropsch-based bio-jet, with each pathway explicitly considering the co-production of renewable diesel and renewable gasoline alongside the bio-jet. Without explicit actions or technology changes to offset the technology cost of bio-jet, scenarios where bio-jet displaces refined liquids result in higher aviation fuel prices (ranging from a 25% increase to 120% increase by mid-century) and lower demand (ranging from −14% to −43%). Corn ethanol will play an important role in the US if large scale amounts of bio-jet are to be produced with smaller effects on demand and prices. While scenarios with high levels of bio-jet availability without the availability of ETJ in the US can significantly reduce emissions in the aviation sector, these reductions are achieved more through the reduction in overall aviation fuel demand rather than technology adoption.

09 BIOMASS FUELS↗

Synchronization-Avoiding Graph Algorithms

Because they were developed for optimal sequential complexity, classical graph algorithms as found in textbooks have strictly-defined orders of operations. Enforcing a prescribed order of operations, or even an approximate order, in a distributed memory setting requires significant amounts of synchronization, which in turn can severely limit scalability. As a result, new algorithms are typically required to achieve scalable performance, even for solving well-known graph problems. Yet, even in these cases, parallel graph algorithms are written according to parallel programming models that evolved for, e.g., scientific computing, and that still have inherent, and scalability-limiting, amounts of synchronization. In this paper we present a new approach to parallel graph algorithms: synchronization-avoiding algorithms. To eliminate synchronization and its associated overhead, synchronization-avoiding algorithms perform work in an unordered and fully asynchronous fashion in such a way that the result is constantly refined toward its final state. “Wasted” work is minimized by locally prioritizing tasks using problem-dependent task utility metrics. We classify algorithms for graph applications into two broad categories: algorithms with monotonic updates (which evince global synchronization) and algorithms with nonmonotonic updates (which evince vertex-centric synchronization). We apply our approach to both classes and develop novel, synchronization-avoiding algorithms for solving exemplar problems: SSSP and connected components for the former, graph coloring for the latter. We demonstrate that eliminating synchronization in conjunction with effective scheduling policies and optimizations in the runtime results in improved scalability for both classes of algorithms.

Firoz, Jesun S.↗

Large eddy simulation of a turbulent non-reacting spray jet

We performed Large Eddy Simulation (LES) of a turbulent non-reacting n-Heptane spray jet, referred to as Spray H in the Engine Combustion Network (ECN), and executed a data analysis focused on key LES metrics such as fraction of resolved turbulent kinetic energy and similarity index. In the simulation, we used the dynamic structure model for the sub-grid stress, and the Lagrangian-based spray-parcel models coupled with the blob-injection model. The finest mesh-cell size used was characterized by an Adaptive Mesh Refinement (AMR) cell size of 0.0625 mm. To obtain ensemble statistics, we performed 28 numerical realizations of the simulation. Demonstrated by the comparison with experimental data in a previous study [7], this LES has accurately predicted global quantities, such as liquid and vapor penetrations. The analysis in this work shows that 14 realizations of LES are sufficient to provide a reasonable representation of the average flow behavior that is benchmarked against the 28-realization ensemble. With the current mesh, numerical schemes, and sub-grid scale turbulence model, more than 95% of the turbulent kinetic energy is directly resolved in the flow regions of interest. The large-scale flow structures inferred from a statistical analysis reveal a region of disorganized flow around the peripheral region of the spray jet, which appears to be linked to the entrainment process.

Hu, Bing↗

Chooz-A steam generators characterization - 15205

EDF nuclear waste management requires a deep understanding of characterization, classification and waste sorting operations. In fact, French nuclear waste management defines several classes with specific management, treatment and storage facilities. Based on particular criteria, the more the radiological risk of the nuclear waste is important, the more its management will be complex and expensive. During the dismantling of the first French pressurized reactor Chooz-A, decontamination of the primary water circuits (not including the reactor vessel), the steam generators and the pressurizer have been carried out in order to reduce their activity levels. Thanks to these decontamination operations, and a specific characterization methodology, EDF was able to de-classify the 4 steam generators and store them in one piece at the ANDRA Very Low Level Activity disposal facilities instead of the Low and Intermediate Level activity one. This de-classification allowed EDF to avoid important cutting and packaging processes. Therefore, significant dose reductions for the operators and financial gain have been made. To characterize and declare Chooz-A SG activity, EDF-CIDEN uses a methodology defined by the French institute of atomic energy, CEA. The method is based on external gamma spectrometry measurements performed with NaI collimated detectors, associated with MERCURAD simulations providing the transfer functions for the detectors and activity sources. Internal measurements are carried out with a CZT (CdZnTe) probe inside the SG tubes to refine the 3D model. In fact, the primary side represents the main source of activity, and understanding its contamination distribution is important to reduce the model and calculation uncertainties. Measurements eventually provide SG {sup 60}Co global activity, from which the activities of other radionuclides of the spectrum are determined using scaling factors. The final activity declaration takes into account the standard deviation of the measurements in order to cover the uncertainties of the methodology. Thereby, the declaration is 'slightly over-estimative', enough to meet the waste storage facility acceptance criteria. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The June 2012 North American Derecho: A Testbed for Evaluating Regional and Global Climate Modeling Systems at Cloud‐Resolving Scales

Abstract In this paper, we introduce a testbed for evaluating and comparing climate modeling systems at cloud resolving scales using hindcasts of the June 2012 North American derecho. To demonstrate its utility for model intercomparison, the testbed is applied to two models: the regionally‐refined Simple Cloud‐Resolving E3SM Atmosphere Model (SCREAM) at 6.5, 3.25 and 1.625 km grid spacing and the Weather Research and Forecasting (WRF) model with 3.2 and 1.6 km grid spacing. We find the simulation results to be highly sensitive to the initial conditions (ICs), initialization time, and model configurations, with ICs from the Rapid Refresh producing the best simulation. Significant improvement is identified in both models as horizontal grid spacing is refined. While a propagation delay of approximately 2 hr is found in both models, SCREAM at 1.625 km simulates the observed bow echo structure of the derecho well and predicts strong surface gusts that exceed 30 m/s. In comparison, WRF has difficulty producing surface wind over 25 m/s, with wind gusts in WRF 42%–46% lower than in SCREAM. However, WRF has a lower bias in simulating cloud top temperature and extent, but overestimates precipitation intensity. Both models reproduce the observed outgoing longwave radiation spatial patterns well (Pearson correlation >0.88), but, compared with NEXRAD observations, simulate generally larger areas of composite radar reflectivity >40 dBZ and underestimate the precipitating area by ∼47%.

2012 North American derecho↗

Numerical Analysis of Combustion Dynamics in a Full-Scale Rotating Detonation Rocket Engine using Large Eddy Simulations

Large eddy simulations (LES) using detailed chemistry and leveraging adaptive mesh refinement (AMR) are performed to gain insights into the combustion dynamics within a full-scale methane-oxygen non-premixed rotating detonation rocket engine (RDRE) employing impinging discrete injection schemes. In particular, a comparative analysis of two operating conditions corresponding to the same global equivalence ratio but different mass flow rates is carried out to investigate the resultant impact on detonation wave characteristics and RDRE global performance. Multiple co-rotating detonation waves with spatially-distributed wave structure and preferential alignment with the inner wall of the annulus (due to asymmetry in fuel distribution) are encountered under both conditions. Both cases exhibit pre-detonation deflagrative burning in the fill region, while one of the cases shows higher susceptibility to backflow into the feed plenums due to lower plenum pressures. Furthermore, heat release analysis shows that the thrust obtained from the RDRE is closely linked to the distribution of total heat release between detonative and deflagrative combustion. On the other hand, combustion efficiency is associated with the fraction of heat release occurring in fuel-rich versus fuel-lean regions within the RDRE.

33 ADVANCED PROPULSION SYSTEMS↗

LiDAR-Aided Interior Orientation Parameters Refinement Strategy for Consumer-Grade Cameras Onboard UAV Remote Sensing Systems

Unmanned aerial vehicles (UAVs) are quickly emerging as a popular platform for 3D reconstruction/modeling in various applications such as precision agriculture, coastal monitoring, and emergency management. For such applications, LiDAR and frame cameras are the two most commonly used sensors for 3D mapping of the object space. For example, point clouds for the area of interest can be directly derived from LiDAR sensors onboard UAVs equipped with integrated global navigation satellite systems and inertial navigation systems (GNSS/INS). Imagery-based mapping, on the other hand, is considered to be a cost-effective and practical option and is often conducted by generating point clouds and orthophotos using structure from motion (SfM) techniques. Mapping with photogrammetric approaches requires accurate camera interior orientation parameters (IOPs), especially when direct georeferencing is utilized. Most state-of-the-art approaches for determining/refining camera IOPs depend on ground control points (GCPs). However, establishing GCPs is expensive and labor-intensive, and more importantly, the distribution and number of GCPs are usually less than optimal to provide adequate control for determining and/or refining camera IOPs. Moreover, consumer-grade cameras with unstable IOPs have been widely used for mapping applications. Therefore, in such scenarios, where frequent camera calibration or IOP refinement is required, GCP-based approaches are impractical. To eliminate the need for GCPs, this study uses LiDAR data as a reference surface to perform in situ refinement of camera IOPs. The proposed refinement strategy is conducted in three main steps. An image-based sparse point cloud is first generated via a GNSS/INS-assisted SfM strategy. Then, LiDAR points corresponding to the resultant image-based sparse point cloud are identified through an iterative plane fitting approach and are referred to as LiDAR control points (LCPs). Finally, IOPs of the utilized camera are refined through a GNSS/INS-assisted bundle adjustment procedure using LCPs. Seven datasets over two study sites with a variety of geomorphic features are used to evaluate the performance of the developed strategy. The results illustrate the ability of the proposed approach to achieve an object space absolute accuracy of 3–5 cm (i.e., 5–10 times the ground sampling distance) at a 41 m flying height.

47 OTHER INSTRUMENTATION↗

Studies of QCD structure in high-energy collisions

”Studies of QCD structure in high-energy collisions” is a research project in theoretical particle physics at Southern Methodist University funded by US DOE Award DE-SC0013681. The award furnished bridge funding for one year (2015/04/15-2016/03/31) between the periods funded by Nadolsky’s DOE Early Career Research Award DE-SC0003870 (in 2010-2015) and a DOE grant DE-SC0010129 for SMU Department of Physics (starting in April 2016). The primary objective of the research is to provide theoretical predictions for Run-2 of the CERN Large Hadron Collider (LHC). The LHC physics program relies on state-of-the-art predictions in the field of quantum chromodynamics. The main effort of our group went into the global analysis of parton distribution functions (PDFs) employed by the bulk of LHC computations. Parton distributions describe internal structure of protons during ultrarelivistic collisions. A new generation of CTEQ parton distribution functions (PDFs), CT14, was released in summer 2015 and quickly adopted by the HEP community. The new CT14 parametrizations of PDFs were obtained using benchmarked NNLO calculations and latest data from LHC and Tevatron experiments. The group developed advanced methods for the PDF analysis and estimation of uncertainties in LHC predictions associated with the PDFs. We invented and refined a new ’meta-parametrization’ technique that streamlines usage of PDFs in Higgs boson production and other numerous LHC processes, by combining PDFs from various groups using multivariate stochastic sampling. In 2015, the PDF4LHC working group recommended to LHC experimental collaborations to use ’meta-parametrizations’ as a standard technique for computing PDF uncertainties. Finally, to include new QCD processes into the global fits, our group worked on several (N)NNLO calculations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

One Structure, Two Elements—LuGe 2 Superconductor vs Ordinary Metallic Conductor LuSn 2 . A Case Study on How Site-Selective Germanium for Tin Atom Substitution Leads to Modulating of the Charge Distribution

The substitution of chemically similar elements in a given crystal structure is an effective way to enhance physical properties, but the understanding on such improvements is usually impeded because the substitutions are random, and the roles of the different atoms cannot be distinguished by crystallographic symmetry. Herein, we provide a detailed crystallographic analysis and property measurements for the continuous solid solutions LuGe x Sn 2–x (0 < x < 2). The results show that there is no apparent change of the global symmetry, with the end-members LuGe 2 and LuSn 2 , as well as the intermediate LuGe x Sn 2–x compositions adopting the ZrSi 2 type structure (space group Cmcm, Pearson index oC12). Yet, the refinements of the crystal structures from single-crystal X-ray diffraction data show that Ge–Sn atom substitutions are not random, but occur preferentially at the zig-zag chain. The patterned distribution of two group 14 elements leads to a significant variation in chemical bonding and charge ordering within the other structural fragment, the 2D square nets, thereby resulting in tuned electron transport. The enhancement is greater than the typical Bloch-Gruneisen model, and more akin to the parallel-resistor model. Magnetization measurements on single crystals show bulk superconductivity in all LuGe x Sn 2–x samples with shielding fractions as high as 90%. Specific heat data confirm the effect to originate from residual metallic tin in the material, indicating that Sn atom substitutions in the 2D square nets cause disruptions of the hypervalent bonding and local anisotropy, which ultimately leads to vanishing of the superconducting state in the end-member LuGe 2 . This work sheds light on how the complexity in chemical interactions by two different carbon congeners leads to changes in the physical properties and how they can be correlated with the induced charge distribution. Furthermore, these studies also provide a general approach to modulation of charge density, and thus, of emerging physical properties in other classes of intermetallic systems based on the main-group elements of groups 13 to 15.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The GFDL‐CM4X Climate Model Hierarchy, Part I: Model Description and Thermal Properties

We present the GFDL‐CM4X (Geophysical Fluid Dynamics Laboratory Climate Model version 4X) coupled climate model hierarchy. The primary application for CM4X is to investigate ocean and sea ice physics as part of a realistic coupled Earth climate model. CM4X utilizes an updated MOM6 (Modular Ocean Model version 6) ocean physics package relative to CM4.0, and there are two members of the hierarchy: one that uses a horizontal grid spacing of 0.25° (referred to as CM4X‐p25) and the other that uses a grid 0.125° (CM4X‐p125). CM4X also refines its atmospheric grid from the nominally 100 km (cubed sphere C96) of CM4.0–50 km (C192). Finally, CM4X simplifies the land model to allow for a more focused study of the role of ocean changes to global mean climate. CM4X‐p125 reaches a global ocean area mean heat flux imbalance of -0.02 W m -2 within $\mathcal{O}$ (150) years in a pre‐industrial simulation, and retains that thermally equilibrated state over the subsequent centuries. This 1850 thermal equilibrium is characterized by roughly 400 ZJ less ocean heat than present‐day, which corresponds to estimates for anthropogenic ocean heat uptake between 1870 and present‐day. CM4X‐p25 approaches its thermal equilibrium only after more than 1000 years, at which time its ocean has roughly 1100 ZJ more heat than its early 21st century ocean initial state. Furthermore, the root‐mean‐square sea surface temperature bias for historical simulations is roughly 20% smaller in CM4X‐p125 relative to CM4X‐p25 (and CM4.0). We offer the mesoscale dominance hypothesis for why CM4X‐p125 shows such favorable thermal equilibration properties.

54 ENVIRONMENTAL SCIENCES↗

Report of the panel on lithospheric structure and evolution, section 3

The panel concluded that NASA can contribute to developing a refined understanding of the compositional, structural, and thermal differences between continental and oceanic lithosphere through a vigorous program in solid Earth science with the following objectives: determine the most fundamental geophysical property of the planet; determine the global gravity field to an accuracy of a few milliGals at wavelengths of 100 km or less; determine the global lithospheric magnetic field to a few nanoTeslas at a wavelength of 100 km; determine how the lithosphere has evolved to its present state via acquiring geologic remote sensing data over all the continents.

Chase, Clement G.↗

Troubleshooting deep-learner training data problems using an evolutionary algorithm on Summit

Architectural and hyper-parameter design choices can influence deep-learner (DL) model fidelity but can also be affected by malformed training and validation data. However, practitioners may spend significant time refining layers and hyper-parameters before discovering that distorted training data was impeding training progress. We found that an evolutionary algorithm (EA) can be used to diagnose this kind of DL problem. An EA evaluated thousands of DL configurations on Summit that yielded no overall improvement in DL performance, which suggested problems with the training and validation data. We suspected that Contrast Limited Adaptive Histogram Equalization (CLAHE) enhancement that was applied to previously generated digital surface models (DSMs), for which we were training DLs to find errors, had damaged the training data. Subsequent runs with an alternative global normalization yielded significantly improved DL performance. However, the DL Intersection Over Union (IOU) still exhibited consistent sub-par performance, which suggested further problems with the training data and DL approach. Nonetheless, we were able to diagnose this problem within a 12-hour span via Summit runs, which prevented several weeks of unproductive trial-and-error DL configuration refinement, and allowed for a more timely convergence on an ultimately viable solution.

97 MATHEMATICS AND COMPUTING↗

Structure-guided design of a broadly cross-reactive multivalent group a streptococcal vaccine

The M protein of group A streptococci (Strep A) is a major virulence determinant and protective antigen. The N-terminal region of the M protein is variable in sequence, defines the M/emm type, and contains epitopes that elicit opsonic antibodies that protect animals from challenge infections. Although there are >200 M types of Strep A, there is now evidence that structurally related M proteins can be grouped into clusters and that immunity may be cluster-specific in addition to M type-specific. This observation has led to recent studies of structure-based design of multivalent M peptide vaccines to select peptides predicted to cross-react with heterologous M types to improve vaccine coverage. In the current study, we have applied a refined series of peptide structural algorithms to predict immunological cross-reactivity among 117N-terminal M peptides representing the most prevalent M types of Strep A. Based on the results of the structural analyses, in combination with global M type prevalence data, we constructed a 32-valent vaccine containing 19 cross-reactive vaccine candidates predicted to cross-react with 37 heterologous M peptides to which were added 13 type-specific M peptides. Further, the 4-protein recombinant vaccine was immunogenic in rabbits and elicited significant levels of antibodies against 31/32 (97%) vaccine peptides and 28/37 (76%) peptides predicted to cross-react. The vaccine antisera also promoted opsonophagocytic killing of vaccine and cross-reactive M types of Strep A. Based on a recent analysis of M type prevalence of Strep A, the potential global coverage of the 32-valent vaccine is ~90%, ranging from 68% in Africa to 95% in North America. Our results indicate the utility of structure-based design that may be applied to future studies of broadly protective M peptide vaccines.

60 APPLIED LIFE SCIENCES↗

Ensemble-Based Parameter Estimation in a Coupled GCM Using the Adaptive Spatial Average Method

Ensemble-based parameter estimation for a climate model is emerging as an important topic in climate research. And for a complex system such as a coupled ocean–atmosphere general circulation model, the sensitivity and response of a model variable to a model parameter could vary spatially and temporally. An adaptive spatial average (ASA) algorithm is proposed to increase the efficiency of parameter estimation. Refined from a previous spatial average method, the ASA uses the ensemble spread as the criterion for selecting “good” values from the spatially varying posterior estimated parameter values; these good values are then averaged to give the final global uniform posterior parameter. In comparison with existing methods, the ASA parameter estimation has a superior performance: faster convergence and enhanced signal-to-noise ratio.

54 ENVIRONMENTAL SCIENCES↗

A Methodology for the Design and Verification of Globally Asynchronous/Locally Synchronous Architectures

Recent advanced in model-checking have made it practical to formally verify the correctness of many complex synchronous systems (i.e., systems driven by a single clock). However, many computer systems are implemented by asynchronously composing several synchronous components, where each component has its own clock and these clocks are not synchronized. Formal verification of such Globally Asynchronous/Locally Synchronous (GA/LS) architectures is a much more difficult task. In this report, we describe a methodology for developing and reasoning about such systems. This approach allows a developer to start from an ideal system specification and refine it along two axes. Along one axis, the system can be refined one component at a time towards an implementation. Along the other axis, the behavior of the system can be relaxed to produce a more cost effective but still acceptable solution. We illustrate this process by applying it to the synchronization logic of a Dual Fight Guidance System, evolving the system from an ideal case in which the components do not fail and communicate synchronously to one in which the components can fail and communicate asynchronously. For each step, we show how the system requirements have to change if the system is to be implemented and prove that each implementation meets the revised system requirements through modelchecking.

Miller, Steven P.↗