Search NASA⌕ Search

SEARCH · Search NASA

Results for “grid simulator testing”

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 1,243 records · Page 69

NREL's Advanced Distribution Management System (ADMS) Test Bed

NREL's advanced distribution management system (ADMS) test bed is a vendor-neutral resource for evaluating existing and future ADMS functionalities in a realistic laboratory setting. The test bed uses Opal-RT's ePHASORsim and eMEGAsim products interfaced with an OpenDSS power system simulation through the HELICS cosimulation platform developed by multiple DOE laboratories; and with power and controller hardware. NREL will present an overview of the test bed and the technologies included within it, and also some examples of how the test bed has been used to address industry challenges.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Robustness Against Disturbances in Power Systems Under Frequency Constraints

The wide deployment of renewable generation and the gradual decrease in the overall system inertia make modern power grids more vulnerable to transient instabilities and unacceptable frequency fluctuations. Time-domain simulation based assessment of the system robustness against uncertain and stochastic disturbances is extremely time-consuming. In this paper, we develop an alternative approach, which has its roots in the input-output stability analysis for Lur'e systems. Our approach consists of a mathematically rigorous characterization of the external disturbances that the power system is transiently stable and the frequency constraints are not violated. As a result, the derived certificate is efficiently constructed via convex optimization and is shown to be non-conservative for different IEEE test cases.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Aggregation and Grid Security Workshop Report

The Aggregation and Grid Security Workshop - held on June 17-18, 2025, National Laboratory of the Rockies (NLR) in Golden, Colorado - brought together approximately 40 external stakeholders from the energy sector, including VPP owner/operators, aggregators, OEMs, utilities, testing & certification labs, trade associations, and cybersecurity vendors. Led by key facilitators, the workshop focused on addressing cybersecurity challenges and enhancing grid resilience for aggregated Distributed Energy Resources (DERs) and Virtual Power Plants (VPPs). The workshop was catalyzed by recognition that traditional, rearward-looking regulatory frameworks are insufficient to keep pace with technological change. There is a "missing understanding" of risk, an "absent security basis" for managing it, and an "untenable responsibility" due to unclear ownership and requirements. The workshop aimed to shift the mindset from reacting to past crises to proactively preparing for emerging threats, fostering forward resilience through risk simulation and collaborative action. This report summarizes the outcomes of the workshop, marking it a significant step toward a secure, reliable and affordable energy future.

14 SOLAR ENERGY↗

Flow prediction for propfan engine installation effects on transport aircraft at transonic speeds

An Euler-based method for aerodynamic analysis of turboprop transport aircraft at transonic speeds has been developed. In this method, inviscid Euler equations are solved over surface-fitted grids constructed about aircraft configurations. Propeller effects are simulated by specifying sources of momentum and energy on an actuator disc located in place of the propeller. A stripwise boundary layer procedure is included to account for the viscous effects. A preliminary version of an approach to embed the exhaust plume within the global Euler solution has also been developed for more accurate treatment of the exhaust flow. The resulting system of programs is capable of handling wing-body-nacelle-propeller configurations. The propeller disks may be tractors or pushers and may represent single or counterrotation propellers. Results from analyses of three test cases of interest (a wing alone, a wing-body-nacelle model, and a wing-nacelle-endplate model) are presented. A user's manual for executing the system of computer programs with formats of various input files, sample job decks, and sample input files is provided in appendices.

Samant, S. S.↗

The Diurnal Cycle in TOGA-COARE: Regional Scale Model Simulations

The diurnal variation of precipitation processes over the tropics is a well-known phenomenon and has been studied using surface rainfall data, radar reflectivity data, and satellite-derived cloudiness and precipitation. Recently, analyzed observations from Tropical Oceans and Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE) in the tropical western Pacific ocean to study the relevant mechanisms producing diurnal variation of precipitation. They found that the diurnal Sea surface temperature (SST) cycle is important for afternoon showers in the undisturbed periods and diurnal radiative processes for nocturnal rainfall. Cloud resolving models (CRMS) have been used to determine the mechanisms associated with diurnal variation of precipitating processes. CRMs allow explicit cloud-radiation and air-sea interactive processes. However, CRMs can be only used for idealized simulations (i.e., no feedback between clouds and their embedded large-scale environments; cyclic lateral boundary conditions and idealized initial conditions). In this study, the Penn State/NCAR Mesoscale Model (MM5) with improved physics (i.e., cloud microphysics, radiation, land-soil-vegetation-surface processes, and TOGA COARE flux scheme) and a multiple level nesting technique (covers the TOGA COARE LSA/IFA with a 54 km grid and can nest down to 18, 6 and possibly even 2 km) will be adopted for studying the diurnal variations of rainfall. We will examine precipitation processes over open ocean and over land. We will also perform sensitivity tests to determine how the radiative forcing and diurnal SST cycle affects the development of convection.

Tao, W.-K.↗

Flexible Service Contracting for Risk Management within Integrated Transmission and Distribution Systems

The general objective of our project has been to investigate the ability of Independent Distribution System Operators (IDSOs), functioning as linkage agents for Integrated Transmission and Distribution (ITD) systems, to facilitate the flexible availability and usage of reserve in support of ITD system operations. This general objective is in accordance with Order 2222 of the U.S. Federal Energy Regulatory Commission (FERC), titled “Participation of Distributed Energy Resource Aggregations in Markets Operated by Regional Transmission Organizations and Independent System Operators”. The Final Rule for FERC Order 2222 was issued on September 17, 2020. The primary contribution of our project is that we have formulated an innovative energy management Transactive Energy Design (TES) approach for ITD systems that provides promising support for our project objective. Specifically, we have developed a new type of IDSO-managed TES design for distribution system operations, as well as new types of swing contracts permitting IDSOs to participate in transmission system operations as providers of reserve harnessed from distribution system resources in return for appropriate compensation. Together, these design elements constitute a scalable market-based approach facilitating efficient reserve procurement for ITD system operations from a fuller range of power resources. The efficacy of our approach has been demonstrated by means of detailed conceptual analyses as well as test case simulations. To implement the latter, we have developed the ITD TES Platform V2.0, a computational platform that permits the modeling and software implementation of ITD system operations over time. The Electric Reliability Council of Texas (ERCOT) energy region has been used as the empirical anchor for the development of this platform. Our conceptual and test-case work has been reported in a Wiley/IEEE Press book, two refereed book chapters, and seven refereed journal articles. The key components comprising the ITD TES Platform V2.0 have been released as documented open-source software at online GitHub repositories.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Towards a unified nonlocal, peridynamics framework for the coarse-graining of molecular dynamics data with fractures

Molecular dynamics (MD) has served as a powerful tool for designing materials with reduced reliance on laboratory testing. However, the use of MD directly to treat the deformation and failure of materials at the mesoscale is still largely beyond reach. In this work, we propose a learning framework to extract a peridynamics model as a mesoscale continuum surrogate from MD simulated material fracture data sets. Firstly, we develop a novel coarse-graining method, to automatically handle the material fracture and its corresponding discontinuities in the MD displacement data sets. Inspired by the weighted essentially non-oscillatory (WENO) scheme, the key idea lies at an adaptive procedure to automatically choose the locally smoothest stencil, then reconstruct the coarse-grained material displacement field as the piecewise smooth solutions containing discontinuities. Then, based on the coarse-grained MD data, a two-phase optimization-based learning approach is proposed to infer the optimal peridynamics model with damage criterion. In the first phase, we identify the optimal nonlocal kernel function from the data sets without material damage to capture the material stiffness properties. Then, in the second phase, the material damage criterion is learnt as a smoothed step function from the data with fractures. As a result, a peridynamics surrogate is obtained. As a continuum model, our peridynamics surrogate model can be employed in further prediction tasks with different grid resolutions from training, and hence allows for substantial reductions in computational cost compared with MD. We illustrate the efficacy of the proposed approach with several numerical tests for the dynamic crack propagation problem in a single-layer graphene. Our tests show that the proposed data-driven model is robust and generalizable, in the sense that it is capable of modeling the initialization and growth of fractures under discretization and loading settings that are different from the ones used during training.

97 MATHEMATICS AND COMPUTING↗

Simulations of the Modular Axisymmetric Scramjet Test Rig Under Reacting Flow Conditions

Simulations of the Air Force Research Laboratory (AFRL) Modular Axisymmetric Scramjet Test Rig (MASTeR) are presented. MASTeR is a parametric test article capable of investigating various scramjet cavity flameholder designs and fueling strategies with the goal to characterize and optimize flameholding capability. In the current work, three cavity aspect ratios, and two fueling strategies (upstream and in-cavity) with ethylene at a nominal facility pressure and temperature condition are evaluated. The reacting simulations are performed, and the resulting flow characteristics are discussed. For each configuration and fueling strategy, cavity residence time, entrainment rate, and fuel-air equivalence ratio are computed. The MASTeR geometry is defined in Engineering Sketch Pad (ESP) and the simulations use a sketch-to-solution (S2S) automated unstructured grid adaptation tool in VULCAN-CFD. This tool automatically generates a simulation grid from the ESP geometry and systematically adapts it to the numerical solution based on the Hessian error estimate of a specified flow field parameter. Reynolds averaged simulations (RAS) are used with a two-equation linear eddy viscosity and diffusivity model. The resulting database can be compared with the experimental data as those become available, and explored to develop models for cavity performance for scramjet propulsion design applications.

hypersonics↗

Stem Inc. SunShot Incubator Program Final Technical Report

In this Energy Storage Control Algorithms project, Stem sought to develop tools and control algorithms to increase the value and reduce balance-of-system and grid integration costs associated with adding distributed solar generation to the grid. These advances fell under the headings SolarScope and SolarController. Stem sought to create initial market traction with a fully commercialized product for the solar industry to size storage systems (SolarScope) as well as a solar intermittency-mitigation framework for utilities (SolarController) in the course of the project. The company sought to align strategic growth plans and enable the rollout of the products to broader audiences in multiple geographic regions by leveraging the major solar companies in the national market as partners. Both final products were both intended to be commercialized. They are: SolarScope: Analysis tool to identify viable PV + storage projects and thereby expedite the sales and interconnection processes. SolarScope combines customer load data, PV production estimates, utility rate tariff, and simulated storage into a simple user interface for PV developers. Developers can easily identify viable solar + storage sites without the need for complex and time consuming, site-by-site spreadsheet modeling. SolarContoller: Tool to autonomously dispatch distributed storage in order to mitigate voltage fluctuation and reduce curtailment. SolarController co-optimizes, in real time, storage dispatch for circuit stability and curtailment reduction, enabling higher penetrations of PV. SolarController is automated, not requiring utility dispatch or management, as Stem hardware senses grid voltage, frequency, customer load, PV production, and power factor. In the end the two products met with different outcomes. SolarScope was tested by potential users, and continues to be used as a foundational platform for partnership with key solar industry partners. SolarController, on the other hand, was successful in lab testing but was not commercialized due to a lack of marketability and lack of interested customer base. Together the development of these two products marked a material step forward for Stem; and a new milestone along the pathway of integration for the solar and storage industries. SolarScope is leading to real, out-of-the-lab project development in storage + solar for the commercial customer sector. Meanwhile SolarController has opened the eyes of regulators and utility executives alike to the potential of distributed solar and by doing so, has moved the conversation forward for the integration of distributed energy resources more broadly on the grid.

25 ENERGY STORAGE↗

VOLTTRON-Based System for Providing Ancillary Services with Residential Building Loads

Ancillary services entail controlled modulation of building equipment to maintain a stable balance of generation and load in the power system. Ancillary services include frequency regulation and contingency reserves, whose acting time ranges from several seconds to several minutes. Many pilot studies have been implemented to use industrial loads to provide ancillary services, and some have explored services from commercial building loads or electric vehicle charging loads. Residential loads, such as space conditioning and water heating, represent a largely untapped resource for providing ancillary services. The residential building sector accounts for a significant fraction of the total electricity use in the United States. Many loads in residential buildings are flexible and could potentially be curtailed or shifted at the request of the grid. However, there are many barriers that prevent residential loads being widely used for ancillary services. One of the major technical barriers is the lack of communication capabilities between end-use devices and the grid. End-use devices need to be able to receive the automatic generation control (AGC) signal from the grid operator and supply certain types of telemetry to verify response. With the advance of consumer electronics, communication-enabled, or 'connected,' residential equipment has emerged to overcome the communication barrier. However, these end-use devices have introduced a new interoperability challenge due to the existence of numerous standards and communication protocols among different end devices. In this paper, we present a VOLTTRON-based system that overcomes these technical challenges and provides ancillary services with residential loads. VOLTTRON is an open-source control and sensing platform for building energy management, facilitating interoperability solutions for end devices. We have developed drivers to communicate and control different types of end devices through standard-based interfaces, manufacturer-provided application programming interfaces, and proprietary communication interfaces. We document the ability to manage nine appliances, using four different standards or proprietary communication methods. A hardware-in-the-loop test was performed in a laboratory environment where the loads of a laboratory home and a large number of simulated homes are controlled by an aggregator. Upon receipt of an AGC signal, the VOLTTRON home energy management system (HEMS) of the laboratory home adjusts the end-device controls based on the comfort criteria set by the end users and sends telemetry to the aggregator to verify response. The aggregator then sends the AGC signal to other simulated homes in attempts to match the utility request as closely as possible. Frequency regulation is generally considered a higher value service than other ancillary services but it is also more challenging due to the constraint of short response time. A frequency regulation use case has been implemented with the regulation signals sent every 10 seconds. Experimental results indicate that the VOLTTRON-controlled residential loads are able to be controlled with sufficient fidelity to enable an aggregator to meet frequency regulation requirements. Future work is warranted, such as understanding the impact of this type of control on equipment life, and market requirements needed to open up residential loads to ancillary service aggregators.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

PRIME: a real-time cyber-physical systems testbed: from wide-area monitoring, protection, and control prototyping to operator training and beyond

As the power grid continues to evolve with advanced wide-area monitoring, protection, and control (WAMPAC) algorithms, there is an increasing need for realistic testbed environments with industry-grade software and hardware-in-the-loop (HIL) to perform verification and validation studies. Such testbed environments serve as ideal platforms to perform WAMPAC prototyping, operator training, and also to study the impacts of different types of cyberattack scenarios on the operation of the grid. In this paper, we introduce PRIME (PNNL cybeR physIcal systeMs tEstbed): the testbed that integrates real-time transmission system simulator with commercial industry grade energy management system (EMS) software and remote hardware-in-the-loop (RHIL). PRIME is an end-to-end, modular testbed that allows high-fidelity RHIL experimentation of a power system. We present two detailed case studies (fault location and clearing in transmission system, and operator training) to showcase the capabilities of our PRIME testbed. Finally, we briefly discuss some of the potential limitations of our testbed in terms of scalability and flexibility to setup larger test systems and identify directions for future work to address those limitations.

42 ENGINEERING↗

Models and methods for assessing the value of HVDC and MVDC technologies in modern power grids

This report reflects the results of U.S. Department of Energy’s (DOE) Grid Modernization project 0074 “Models and methods for assessing the value of HVDC [high-voltage direct current] and MTDC [multi-terminal direct current] technologies in modern power grids.” The work was done by the Pacific Northwest National Laboratory (PNNL) and Oak Ridge National Laboratory (ORNL) in cooperation with Mid-Continent Independent System Operator (MISO) and Siemens. The main motivation of this study was to show the benefit of using direct current (DC) systems larger than those in existence today as they overlap with the alternating current (AC) systems. Proper use of their flexibility in terms of active/reactive power control and fast response can provide much-needed services to the grid at the same time as moving large blocks of energy to take advantage of cost diversity. Ultimately, the project’s success will enable decision-makers and investors to make well-informed decisions regarding this use of DC systems. This project showed the technical feasibility of HVDC macrogrid for frequency control and congestion relief in addition to bulk power transfers. Industry-established models for commonly used technologies were employed, along with high-fidelity models for recently developed HVDC converter technologies; like the modular multilevel converters (MMCs), a voltage source converters (VSC). Models for General Electric Positive Sequence Load Flow (GE PSLF) and Siemens Power System Simulator (PSS/E), widely used analysis programs, were for the first time adapted to include at the same time both Western Electricity Coordinating Council (WECC) and Eastern Interconnection (EI), the two largest North American interconnections. The high-fidelity models and their control were developed in detail for MMC system and extended to HVDC systems in point-to-point and in three-node multi-terminal configurations. Using a continental-level mixed AC-DC grid model, and using a HVDC macrogrid power flow and transient stability model, the results showed that the HVDC macrogrid relieved congestion and mitigated loop flows in AC networks, and provided up to 24% improvement in frequency responses. These are realistic studies, based on the 2025 heavy summer and EI multi-regional modeling working group (MMWG) 2026 summer peak cases. This work developed high-fidelity models and simulation algorithms to understand the dynamics of MMC. The developed models and simulation algorithms are up to 25 times faster than the existing algorithms. Models and control algorithms for high-fidelity models were designed and tested for point-to-point and multi-terminal configurations. The multi-terminal configuration was tested connecting simplified models of EI, WI, and Electric Reliability Council of Texas (ERCOT). The developed models showed up to 45% improvement in frequency response with the connection of all the three asynchronous interconnections in the United States using fast and advanced DC technologies like the multi-terminal MMC-DC system. Future work will look into developing high-fidelity models of other advanced DC technologies, combining high-fidelity models with the continental-level model, incorporating additional services. More scenarios involving large-scale HVDC and MTDC will be evaluated.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A parallel-in-time algorithm for variable step multistep methods

This paper presents a multigrid reduction in time (MGRIT) algorithm for achieving time parallelism using multistep backward difference formula (BDF) methods on variably-spaced temporal grids. This MGRIT approach transforms the linear multistep methods into single step methods applied to groups of time steps. Stability considerations are addressed through lowering of the order on coarse grids. The methods are presented for fixed and variable time step formulations. Furthermore, numerical results show moderate speedups for the heat equation as well as two IEEE power grid test problems characterized by nonlinear differential algebraic equations (DAE) of index 1. The performance of MGRIT with BDF is compared to MGRIT with Runge-Kutta for both the fixed and variable step methods of the same order.

97 MATHEMATICS AND COMPUTING↗

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗

A Posteriori Study of a DNS Database Describing Super critical Binary-Species Mixing

Currently, the modeling of supercritical-pressure flows through Large Eddy Simulation (LES) uses models derived for atmospheric-pressure flows. Those atmospheric-pressure flows do not exhibit the particularities of high densitygradient magnitude features observed both in experiments and simulations of supercritical-pressure flows in the case of two species mixing. To assess whether the current LES modeling is appropriate and if found not appropriate to propose higher-fidelity models, a LES a posteriori study has been conducted for a mixing layer that initially contains different species in the lower and upper streams, and where the initial pressure is larger than the critical pressure of either species. An initially-imposed vorticity perturbation promotes roll-up and a double pairing of four initial span-wise vortices into an ultimate vortex that reaches a transitional state. The LES equations consist of the differential conservation equations coupled with a real-gas equation of state, and the equation set uses transport properties depending on the thermodynamic variables. Unlike all LES models to date, the differential equations contain, additional to the subgrid scale (SGS) fluxes, a new SGS term that is a pressure correction in the momentum equation. This additional term results from filtering of Direct Numerical Simulation (DNS) equations, and represents the gradient of the difference between the filtered pressure and the pressure computed from the filtered flow field. A previous a priori analysis, using a DNS database for the same configuration, found this term to be of leading order in the momentum equation, a fact traced to the existence of high-densitygradient magnitude regions that populated the entire flow; in the study, models were proposed for the SGS fluxes as well as this new term. In the present study, the previously proposed constantcoefficient SGS-flux models of the a priori investigation are tested a posteriori in LES, devoid of or including, the SGS pressure correction term. The present pressure-correction model is different from, and more accurate as well as less computationally intensive than that of the a priori study. The constant-coefficient SGS-flux models encompass the Smagorinsky (SMC), in conjunction with the Yoshizawa (YO) model for the trace, the Gradient (GRC) and the Scale Similarity (SSC) models, all exercised with the a priori study constant coefficients calibrated at the transitional state. The LES comparison is performed with the filtered- and-coarsened (FC) DNS, which represents an ideal LES solution. Expectably, an LES model devoid of SGS terms is shown to be considerably inferior to models containing SGS effects. Among models containing SGS effects, those including the pressure-correction term are substantially superior to those devoid of it. The sensitivity of the predictions to the initial conditions and grid size are also investigated. Thus, it has been discovered that, additional to the atmospheric-pressure models currently used, a new model is necessary to simulate supercritical-pressure flows. This model depends on the thermodynamic characteristics of the chemical species involved.

Bellan, Josette↗

Sub-grid scale model classification and blending through deep learning

In this article we detail the use of machine learning for spatio-temporally dynamic turbulence model classification and hybridization for large eddy simulations (LES) of turbulence. Here, our predictive framework is devised around the determination of local conditional probabilities for turbulence models that have varying underlying hypotheses. As a first deployment of this learning, we classify a point on our computational grid as that which requires the functional hypothesis, the structural hypothesis or no modelling at all. This ensures that the appropriate model is specified from a priori knowledge and an efficient balance of model characteristics is obtained in a particular flow computation. In addition, we also utilize the conditional-probability predictions of the same machine learning to blend turbulence models for another hybrid closure. Our test case for the demonstration of this concept is given by Kraichnan turbulence, which exhibits a strong interplay of enstrophy and energy cascades in the wavenumber domain. Our results indicate that the proposed methods lead to robust and stable closure and may potentially be used to combine the strengths of various models for complex flow phenomena prediction.

42 ENGINEERING↗

Experimental and analytical study of close-coupled ventral nozzles for ASTOVL aircraft

Flow in a generic ventral nozzle system was studied experimentally and analytically with a block version of the PARC3D computational fluid dynamics program (a full Navier-Stokes equation solver) in order to evaluate the program's ability to predict system performance and internal flow patterns. For the experimental work a one-third-size model tailpipe with a single large rectangular ventral nozzle mounted normal to the tailpipe axis was tested with unheated air at steady-state pressure ratios up to 4.0. The end of the tailpipe was closed to simulate a blocked exhaust nozzle. Measurements showed about 5 1/2 percent flow-turning loss, reasonable nozzle performance coefficients, and a significant aftward axial component of thrust due to flow turning loss, reasonable nozzle performance coefficients, and a significant aftward axial component of thrust due to flow turning more than 90 deg. Flow behavior into and through the ventral duct is discussed and illustrated with paint streak flow visualization photographs. For the analytical work the same ventral system configuration was modeled with two computational grids to evaluate the effect of grid density. Both grids gave good results. The finer-grid solution produced more detailed flow patterns and predicted performance parameters, such as thrust and discharge coefficient, within 1 percent of the measured values. PARC3D flow visualization images are shown for comparison with the paint streak photographs. Modeling and computational issues encountered in the analytical work are discussed.

Mcardle, Jack G.↗