Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “software behavior computation”

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 55 records · Page 3

Virtual Observation System for Earth System Model: An Application to ACME Land Model Simulations

Investigating and evaluating physical-chemical-biological processes within an Earth system model (EMS) can be very challenging due to the complexity of both model design and software implementation. A virtual observation system (VOS) is presented to enable interactive observation of these processes during system simulation. Based on advance computing technologies, such as compiler-based software analysis, automatic code instrumentation, and high-performance data transport, the VOS provides run-time observation capability, in-situ data analytics for Earth system model simulation, model behavior adjustment opportunities through simulation steering. A VOS for a terrestrial land model simulation within the Accelerated Climate Modeling for Energy model is also presented to demonstrate the implementation details and system innovations.

58 GEOSCIENCES↗

An Unsteady State Gas Transport Model in Matrix Laboratory (MATLAB) for PWR Gas Management

Reactor Coolant System (RCS) hydrogen gas management can be a challenge at Pressurized Water Reactors (PWRs) during shutdown, startup, and steady state operation. Poor hydrogen control can result in higher general corrosion rates for RCS components and increases in critical path time for shutdown or startup. Also, performance indicator penalties may be accrued by entry into Electric Power Research Institute (EPRI) Action Level conditions or by accrual of Chemistry Effectiveness Index (CEI) penalties established by the Institute of Nuclear Power Operations (INPO). Furthermore, poor hydrogen control may increase industrial safety risk due to the creation of an explosive environment. Part of the hydrogen management difficulty arises from the interaction of several phenomena that require complex simulations to fully understand. While hydrogen gas transport models have been developed, the application of these models is limited by the computational intensity required to accurately model the dynamism inherent in plant shutdown or startup. Subsequently, a gas transport model is developed in a suitable computational program to meet these needs and provide input to an optimized RCS hydrogen management strategy for all operational modes. The gas transport model is used to model RCS hydrogen behavior for an example PWR during an End of Cycle (EOC) Forced Outage (FO) startup, a refueling outage (RFO) shutdown, and during nominal steady state operation. The computational software selected is Matrix Laboratory (MATLAB{sup R})

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Coastal Ohio Wind Project

The Coastal Ohio Wind Project intends to address problems that impede deployment of wind turbines in the coastal and offshore regions of Northern Ohio. The project evaluates different wind turbine designs and the potential impact of offshore turbines on migratory and resident birds by developing multidisciplinary research, which involves wildlife biology, electrical and mechanical engineering, and geospatial science. Firstly, the project conducts cost and performance studies of two- and three-blade wind turbines using a turbine design suited for the Great Lakes. The numerical studies comprised an analysis and evaluation of the annual energy production of two- and three-blade wind turbines to determine the levelized cost of energy. This task also involved wind tunnel studies of model wind turbines to quantify the wake flow field of upwind and downwind wind turbine-tower arrangements. The experimental work included a study of a scaled model of an offshore wind turbine platform in a water tunnel. The levelized cost of energy work consisted of the development and application of a cost model to predict the cost of energy produced by a wind turbine system placed offshore. The analysis found that a floating two-blade wind turbine presents the most cost effective alternative for the Great Lakes. The load effects studies showed that the two-blade wind turbine model experiences less torque under all IEC Standard design load cases considered. Other load effects did not show this trend and depending on the design load cases, the two-bladed wind turbine showed higher or lower load effects. The experimental studies of the wake were conducted using smoke flow visualization and hot wire anemometry. Flow visualization studies showed that in the downwind turbine configuration the wake flow was insensitive to the presence of the blade and was very similar to that of the tower alone. On the other hand, in the upwind turbine configuration, increasing the rotor blade angle of attack reduced the wake size and enhanced the vortices in the flow downstream of the turbine-tower compared with the tower alone case. Mean and rms velocity distributions from hot wire anemometer data confirmed that in a downwind configuration, the wake of the tower dominates the flow, thus the flow fields of a tower alone and tower-turbine combinations are nearly the same. For the upwind configuration, the mean velocity shows a narrowing of the wake compared with the tower alone case. The downwind configuration wake persisted longer than that of an upwind configuration; however, it was not possible to quantify this difference because of the size limitation of the wind tunnel downstream of the test section. The water tunnel studies demonstrated that the scale model studies could be used to adequately produce accurate motions to model the motions of a wind turbine platform subject to large waves. It was found that the important factors that affect the platform is whether the platform is submerged or surface piercing. In the former, the loads on the platform will be relatively reduced whereas in the latter case, the structure pierces the wave free surface and gains stiffness and stability. The other important element that affects the movement of the platform is depth of the sea in which the wind turbine will be installed. Furthermore, the wildlife biology component evaluated migratory patterns by different monitoring systems consisting of marine radar, thermal IR camera and acoustic recorders. The types of radar used in the project are weather surveillance radar and marine radar. The weather surveillance radar (1988 Doppler), also known as Next Generation Radar (NEXRAD), provides a network of weather stations in the US. Data generated from this network were used to understand general migratory patterns, migratory stopover habitats, and other patterns caused by the effects of weather conditions. At a local scale our marine radar was used to complement the datasets from NEXRAD and to collect additional monitoring parameters such as passage rates, flight paths, flight directions, and flight altitudes of nocturnal migrating species. Our work focused on the design and development of custom built marine radar that used t-bar and parabolic dish antennas. The marine radar used in the project was Furuno (XANK250) which was coupled with a XIR3000B digitizing card from Russell Technologies for collection of the radar data. The radar data was processed by open source radR processing software using different computational techniques and methods. Additional data from thermal IR imaging cameras were collected to detect heat emitted from objects and provide information on movements of birds and bats, data which we used for different animal flight behavior analysis. Lastly, the data from the acoustic recorders were used to provide the number of bird calls for assessing patterns and peak passage rates during migration. The development of the geospatial database included collection of different data sources that are used to support offshore wind turbine development. Many different data sets were collected and organized using initial version of web-based repository software tools that can accommodate distribution of rectified pertinent data sets such as the lake depth, lake bottom engineering parameters, extent of ice, navigation pathways, wind speed, important bird habitats, fish efforts and other layers that are relevant for supporting robust offshore wind turbine developments. Additional geospatial products developed during the project included few different prototypes for offshore wind farm suitability which can involve different stakeholders and participants for solving complex planning problems and building consensus. Some of the prototypes include spatial decision support system (SDSS) for collaborative decision making, a web-based Participatory Geographic Information System (PGIS) framework for evaluating importance of different decision alternatives using different evaluation criteria, and an Android application for collection of field data using mobile and tablet devices . In summary, the simulations of two- and three-blade wind turbines suggested that two-bladed machines could produce comparable annual energy as the three-blade wind turbines but have a lighter tower top weight, which leads to lower cost of energy. In addition, the two-blade rotor configuration potentially costs 20% less than a three blade configuration that produces the same power at the same site. The cost model analysis predicted a potential cost savings of approximately 15% for offshore two-blade wind turbines. The foundation design for a wind turbine in Lake Erie is likely to be driven by ice loads based on the currently available ice data and ice mechanics models. Hence, for Lake Eire, the cost savings will be somewhat smaller than the other lakes in the Great Lakes. Considering the size of cranes and vessels currently available in the Great Lakes, the cost optimal wind turbine size should be 3 MW, not larger. The surveillance data from different monitoring systems suggested that bird and bat passage rates per hour were comparable during heavy migrations in both spring and fall seasons while passage rates were significantly correlated to wind directions and wind speeds. The altitude of migration was higher during heavy migrations and higher over water relative to over land. Notable portions of migration on some spring nights occurred parallel the shoreline, often moving perpendicular to southern winds. The birds approaching the Western basin have a higher propensity to cross than birds approaching the Central basin of Lake Erie and as such offshore turbine development might be a better option further east towards Cleveland than in the Western basin. The high stopover density was more strongly associated with migration volume the following night rather than the preceding night. The processed mean scalar wind speeds with temporal resolutions as fine as 10-minute intervals near turbine height showed that August is the month with the weakest winds while December is the month, which typically has the strongest winds. The ice data suggests that shallow western basin of Lake Erie has higher ice cover duration many times exceeding 90 days during some winters.

17 WIND ENERGY↗

Exploration of multifidelity UQ sampling strategies for computer network applications.

Network modeling is a powerful tool to enable rapid analysis of complex systems that can be challenging to study directly using physical testing. Here, two approaches are considered: emulation and simulation. The former runs real software on virtualized hardware, while the latter mimics the behavior of network components and their interactions in software. Although emulation provides an accurate representation of physical networks, this approach alone cannot guarantee the characterization of the system under realistic operative conditions. Operative conditions for physical networks are often characterized by intrinsic variability (payload size, packet latency, etc.) or a lack of precise knowledge regarding the network configuration (bandwidth, delays, etc.); therefore uncertainty quantification (UQ) strategies should be also employed. UQ strategies require multiple evaluations of the system with a number of evaluation instances that roughly increases with the problem dimensionality, i.e., the number of uncertain parameters. It follows that a typical UQ workflow for network modeling based on emulation can easily become unattainable due to its prohibitive computational cost. In this paper, a multifidelity sampling approach is discussed and applied to network modeling problems. The main idea is to optimally fuse information coming from simulations, which are a low-fidelity version of the emulation problem of interest, in order to decrease the estimator variance. By reducing the estimator variance in a sampling approach it is usually possible to obtain more reliable statistics and therefore a more reliable system characterization. Several network problems of increasing difficulty are presented. For each of them, the performance of the multifidelity estimator is compared with respect to the single fidelity counterpart, namely, Monte Carlo sampling. For all the test problems studied in this work, the multifidelity estimator demonstrated an increased efficiency with respect to MC.

97 MATHEMATICS AND COMPUTING↗

Computational Fluid Dynamics Modeling of the Fuel Reactor in NETL's 50 kW th Chemical Looping Facility

The National Energy Technology Laboratory (NETL) has explored chemical looping in its 50 kW th facility using a number of oxygen carriers. In this work, the results for methane conversion in the fuel reactor with a hematite iron ore as the oxygen carrier are analyzed. The experimental results are compared to predictions using CPFD's barracuda computational fluid dynamics (CFD) code with kinetics derived from the analysis of fixed bed data. It has been found through analytical techniques from thermal gravimetric analysis data as well as the same fixed bed data that the kinetics for the methane–hematite reaction follows a nucleation and growth or Johnson–Mehl–Avrami (JMA) reaction mechanism. barracuda does not accept nucleation and growth kinetics; however, there is enough sufficient variability of the solids dependence within the software such that the nucleation and growth behavior can be mimicked. This paper presents the method to develop the pseudo-JMA kinetics for barracuda extracted from the fixed bed data and then applies these values to the fuel reactor data to compare the computational results to experimental data obtained from 50 kW th unit for validation. Lastly, a fuel reactor design for near complete conversion is proposed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Determination of anharmonic free energy contributions: Low temperature phases of the Lennard-Jones system

We investigate a general method to calculate the free energy of crystalline solids by considering the harmonic approximation and quasistatically switching the anharmonic contribution. The advantage of this method is that the harmonic approximation provides an already very accurate estimate of the free energy, and therefore the anharmonic term is numerically very small and can be determined to high accuracy. We further show that the anharmonic contribution to the free energy satisfies a number of exact inequalities that place constraints on its magnitude and allows approximate but fast and accurate estimates. The method is implemented into a readily available general software by combining the code HOODLT (Highly Optimized Object Oriented Dynamic Lattice Theory) for the harmonic part and the molecular dynamics (MD) simulation package HOOMD-blue for the anharmonic part. We use the method to calculate the low temperature phase diagram for Lennard-Jones particles. We demonstrate that hcp is the equilibrium phase at low temperature and pressure and obtain the coexistence curve with the fcc phase, which exhibits reentrant behavior. Furthermore, several implications of the method are discussed.

36 MATERIALS SCIENCE↗

MDDC Multi-Length Scale Data Architecture Contribution Report – PNNL, INL, ANL, LANL and ORNL

This report offers a comprehensive view of data streams currently generated at Pacific Northwest National Laboratory, Idaho National Laboratory, Argonne National Laboratory, Los Alamos National Laboratory, and Oak Ridge National Laboratory set to integrate into the evolving Multi-Dimensional Data Correlation framework at Oak Ridge National Laboratory. Developed by the Advanced Materials and Manufacturing Technologies program, the Multi-Dimensional Data Correlation framework serves as a cutting-edge software to manage data relevant to advanced manufacturing and material behavior in advanced reactors. The report defines data streams, highlights their generation methods and visualization methods both for experimental and computational aspects relevant to the Advanced Materials and Manufacturing Technologies project. A logical next step for this work is to integrate the MDDC framework into PNNL’s, INL’s, ANL’s, LANL’s and ORNL’s fabrication, experimentation, and modelling workflows. This would require setting up the MDDC framework at PNNL, INL, ANL, and LANL and integrating it into the data collection and storage for these different activities.

36 MATERIALS SCIENCE↗

Assessing VQLS for Fluid Dynamics on a Hybrid Quantum-HPC Stack

Recent advances in quantum linear solvers offer a promising direction for accelerating extreme scientific computations such as fluid dynamics. However, the deep and complex circuits required by many quantum algorithms limit their practical use on current quantum hardware. The Variational Quantum Linear Solver (VQLS) presents a viable alternative for near-term quantum devices (NISQ), and initial efforts have explored its application to select fluid dynamics problems. In this work, we evaluate the use of VQLS for canonical fluid dynamics problems, aiming to identify pathways for generalizing its implementation across a broader class of systems. We analyze the impact of various circuit ansatz and classical optimizers on solution quality and convergence behavior. Furthermore, we assess the algorithm's feasibility within a hybrid quantum–high-performance computing (HPC) framework by porting it to QFw, a state-of-the-art quantum-HPC software stack. 11This manuscript has been authored by UT-Battelle, LLC, under contract DE-AC05-00OR22725 with the US Department of Energy (DOE). The US government retains and the publisher, by accepting the article for publication, acknowledges that the US government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for US government purposes. DOE will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan. This research used resources of the Oak Ridge Leadership Computing Facility at the Oak Ridge National Laboratory, which is supported by the Office of Science of the US DOE under Contract No. DE-AC05-00OR22725.

Gopalakrishnan Meena, Murali [ORNL] (ORCID:0000000↗

Preliminary Thermal-Hydraulic Analysis of the UTA–2 Subcritical Linear Accelerator–Driven System

The National Nuclear Security Administration’s mission includes establishing a reliable supply of 99 Mo without highly enriched uranium. Oak Ridge National Laboratory (ORNL) supports this objective through collaborative research and development with industrial partners. Niowave Inc., a current partner, is currently designing a subcritical linear accelerator-driven system (ADS) and preparing for the US Nuclear Regulatory Commission’s licensing process. Niowave’s system has the potential to efficiently supply medical radioisotopes. The technology includes a superconducting electron accelerator and a pile of both natural and low-enriched uranium (LEU) targets. The process fissions uranium and many valuable isotopes can then be extracted from the targets. Niowave is currently iterating through conceptual and detailed design processes for several system sizes. This report discusses UTA–2, which is at the demonstration stage. UTA–2 will validate numerical modeling results with experimental measurements before progressing to the detailed design of UTA–3, the commercial-sized ADS. The thermal-hydraulic behavior of the UTA–2 core design was numerically investigated using STAR-CCM+, as described in this report. STAR-CCM+, a state-of-the-art computational fluid dynamics (CFD) software that was commercially developed by Siemens, has an extensive user base and a set of validation studies. It is also compliant with the American Society of Mechanical Engineers’ Nuclear Quality Assurance 1 standard. Three cases are investigated in this report. Case 1 quantified the temperature field within the UTA–2 assembly and water tank using only conduction as the method of thermal energy transport. This simplified approach was overly conservative and yielded wetted cladding temperatures above the coolant saturation temperature. In this case, the maximum temperature of the wetted cladding surface of the highest power rod exceeded the saturation temperature of water by 63.8°C. Because of the overly conservative approach taken in Case 1 and its negative subcooled margin, buoyancy-driven natural circulation flow physics were implemented in Case 2. Adding coolant motion significantly distributed the thermal energy of the system through convective heat transfer. This relatively small amount of convective heat transfer significantly reduced system temperatures and increased the subcooled margin from –63.8 to 61.3°C. This margin confirmed that no boiling was expected during normal operation of UTA–2 at 230 W. Case 3 had no additional physics models but considered an overpower event during which the power of each LEU and natural uranium rod was at its respective peak values. This resulted in a study with the total assembly power equal to 176% of the nominal power of 230 W considered in Cases 1 and 2. The resulting natural circulation flows were slightly enhanced. The subcooled margin decreased slightly to 46.3°C, which is still a significant margin to local boiling of the water within the tank. This margin confirmed that no boiling was expected during an abnormal operation of UTA–2 at 406 W.

07 ISOTOPE AND RADIATION SOURCES↗

Numerical Investigation of Enhanced Dehumidification Processes By Using Dielectrophoresis Principles in Moist Airflows

Dispersed particle-laden flows are encountered in many building and industrial applications, such as flow in a fluidized bed, hydrocarbon transportation in pipelines, and the fouling of air-cooled heat exchangers (Kuruneru et al., 2016; Ray et al., 2019; Wang et al., 2019). Computational fluid dynamic (CFD) models have been developed in recent years to depict particle-fluid and particle-particle interactions in laminar or turbulent flows with increasing accuracy and stability. One particular particle-laden system of interest for moisture control is electrically-enhanced condensation in air and water droplet flows. Electrically-enhanced condensation consists of the use of highly charged water droplets injected in the moist air. The droplets become electric seeds that attract polar water vapor molecules to their surfaces and promote condensation. The nucleation and growth of the charged droplets deplete the vapor phase near a droplet, which is compensated for by the dielectrophoresis flow and diffusion. Dielectrophoresis flow involves surrounding vapor at a distance of about 10 to 100 nm for droplets charged by an electrospray compared to ~2 nm for a single electron charge in a droplet. As the vapor molecules collapse on the surface of the droplets, their initial electrical charge decreases with time due to the neutralization of the ions. While the physics of this phenomena is well known, engineering models for predicting the condensation rates are not available. This work computationally investigates dehumidification of moist airflow in a converging rectangular duct. The objective is to develop an engineering model that predicts water vapor condensation by employing dielectrophoresis principles. We construct a Computational Fluid Dynamics (CFD) model of the duct with electrically-enhanced condensation. The model is implemented in the open-source software OpenFOAM. We utilize the Multi-Phase Particle-In-Cell (MP-PIC) method coupled with a Population Balance Equation (PBE) approach to simulate the particle-laden system. This methodology is an Eulerian-Lagrangian approach used to simulate the droplets' behavior in the humid air. The MP-PIC approach (Andrews and O'Rourke, 1996) mitigates the computational cost by parceling several fundamental particles with similar properties (such as types, sizes, and temperature) into one computational particle. Thus, the billions of particles can be substituted by millions of computational particles without significant loss of information. The PBE was considered with the Lagrangian frame to combine the particle distribution function used in MP-PIC (Kim et al., 2020). This approach preserves mass and energy conservation between the phases in the Eulerian and Lagrangian structures. The PBE in this procedure was directly linked to the discrete parcels, making the simulation of the particle distribution computationally efficient and robust. The MP-PIC-PBE approach used in the present work was applied to the dehumidification of air. Water droplets were injected in the air stream and forced to grow according to experimentally derived correlation. The experiments were conducted on a converging duct with the same geometry and boundary conditions used to build the CFD model. This approach enabled us to approximate the effect of dielectrophoresis phenomena on the droplet and air interface. This presentation will discuss the details of the new CFD model built for the duct, the implementation of the model in OpenFOAM CFD programming language, and the experimental validation of the newly developed model. The results revealed a moderate yet measurable increase in droplet diameter due to water vapor condensation at the vapor-liquid interface of the electrically charged droplets' surface. The seed water droplet particles grew in size by capturing the water vapor in the surrounding air. The OpenFOAM model predicted reductions of humidity in the air from 5 to 10 percent.

Yel Mahi, Maliha↗

SU-F-I-43: A Software-Based Statistical Method to Compute Low Contrast Detectability in Computed Tomography Images

Purpose: The low contrast detectability (LCD) of a CT scanner is its ability to detect and display faint lesions. The current approach to quantify LCD is achieved using vendor-specific methods and phantoms, typically by subjectively observing the smallest size object at a contrast level above phantom background. However, this approach does not yield clinically applicable values for LCD. The current study proposes a statistical LCD metric using software tools to not only to assess scanner performance, but also to quantify the key factors affecting LCD. This approach was developed using uniform QC phantoms, and its applicability was then extended under simulated clinical conditions. Methods: MATLAB software was developed to compute LCD using a uniform image of a QC phantom. For a given virtual object size, the software randomly samples the image within a selected area, and uses statistical analysis based on Student’s t-distribution to compute the LCD as the minimal Hounsfield Unit’s that can be distinguished from the background at the 95% confidence level. Its validity was assessed by comparison with the behavior of a known QC phantom under various scan protocols and a tissue-mimicking phantom. The contributions of beam quality and scattered radiation upon the computed LCD were quantified by using various external beam-hardening filters and phantom lengths. Results: As expected, the LCD was inversely related to object size under all scan conditions. The type of image reconstruction kernel filter and tissue/organ type strongly influenced the background noise characteristics and therefore, the computed LCD for the associated image. Conclusion: The proposed metric and its associated software tools are vendor-independent and can be used to analyze any LCD scanner performance. Furthermore, the method employed can be used in conjunction with the relationships established in this study between LCD and tissue type to extend these concepts to patients’ clinical CT images.

60 APPLIED LIFE SCIENCES↗

DT-HYDRO

The software solves the time dependent, one dimensional (1D) coupled mass and momentum balance equations governing the elastic flow of water through the penstock, turbine and draft tube in a hydroelectric facility using a high order finite volume based method. The numerical method is based on the Kurganov-Tadmor central method paired with the Monotonic Upstream-centered Scheme for Conservation Laws (MUSCL). This solution method accurately resolves the fast transient behavior of the flow, including water hammer. Additionally, the software estimates the full 3D flow field within the turbine chamber in real time, a feat that is made possible by leveraging pre-computed CFD results by utilizing a reduced order modeling method based on an efficient randomized singular value decomposition (SVD) driven proper orthogonal decomposition (POD) with POD-mode weight regression. The reduced order model of the 3D flow is directly coupled to the 1D elastic flow model so the entire flow field through the penstock and turbine system is resolved quickly and with high fidelity.

Gurecky, William [Oak Ridge National Laboratory (O↗

Parallel time integration software

This package implements an optimal-scaling multigrid solver for the (non) linear systems that arise from the discretization of problems with evolutionary behavior. Typically, solution algorithms for evolution equations are based on a time-marching approach, solving sequentially for one time step after the other. Parallelism in these traditional time-integrarion techniques is limited to spatial parallelism. However, current trends in computer architectures are leading twards system with more, but not faster. processors. Therefore, faster compute speeds must come from greater parallelism. One approach to achieve parallelism in time is with multigrid, but extending classical multigrid methods for elliptic poerators to this setting is a significant achievement. In this software, we implement a non-intrusive, optimal-scaling time-parallel method based on multigrid reduction techniques. The examples in the package demonstrate optimality of our multigrid-reduction-in-time algorithm (MGRIT) for solving a variety of parabolic equations in two and three sparial dimensions. These examples can also be used to show that MGRIT can achieve significant speedup in comparison to sequential time marching on modern architectures.

Source record↗

Parallel time integration software

This package implements an optimal-scaling multigrid solver for the (non) linear systems that arise from the discretization of problems with evolutionary behavior. Typically, solution algorithms for evolution equations are based on a time-marching approach, solving sequentially for one time step after the other. Parallelism in these traditional time-integrarion techniques is limited to spatial parallelism. However, current trends in computer architectures are leading twards system with more, but not faster. processors. Therefore, faster compute speeds must come from greater parallelism. One approach to achieve parallelism in time is with multigrid, but extending classical multigrid methods for elliptic poerators to this setting is a significant achievement. In this software, we implement a non-intrusive, optimal-scaling time-parallel method based on multigrid reduction techniques. The examples in the package demonstrate optimality of our multigrid-reduction-in-time algorithm (MGRIT) for solving a variety of parabolic equations in two and three sparial dimensions. These examples can also be used to show that MGRIT can achieve significant speedup in comparison to sequential time marching on modern architectures.

Source record↗

GridOPTICS/GridPACK

GridPACK is a software framework consisting of a set of modules designed to simplify the development of programs that model the power grid and run on parallel, high performance computing platforms. It also contains several fully developed applications, including powerflow, dynamic simulation, state estimation, Kalman filter analysis (dynamic state estimation), contingency analysis and real time path rating. These applications can be used either standalone or as components in more complicated workflows that combine several different types of application together. The framework modules are available as a combination of libraries and software templates and consist of components for setting up and distributing power grid networks, support for modeling the behavior of individual buses and branches in the network, converting the network models to the corresponding algebraic equations, and parallel routines for manipulating and solving large algebraic systems. The framework also contains a module for distributing tasks evenly amongst computing resources, even if individual tasks vary widely in their execution times. Additional modules support input and output, basic statistical analysis of contingency based calculations, distributed data structures, as well as basic profiling and error management.

Palmer, Bruce↗

Streamlining the Design-to-Build Transition with Build-Optimization Software Tools

Scaling-up capabilities for the design, build, and test of synthetic biology constructs holds great promise for the development of new applications in fuels, chemical production, or cellular-behavior engineering. Construct design is an essential component in this process; however, not every designed DNA sequence can be readily manufactured, even using state-of-the-art DNA synthesis methods. Current biological computer-aided design and manufacture tools (bioCAD/CAM) do not adequately consider the limitations of DNA synthesis technologies when generating their outputs. Designed sequences that violate DNA synthesis constraints may require substantial sequence redesign or lead to price-premiums and temporal delays, which adversely impact the efficiency of the DNA manufacturing process. In this paper, we have developed a suite of build-optimization software tools (BOOST) to streamline the design-build transition in synthetic biology engineering workflows. BOOST incorporates knowledge of DNA synthesis success determinants into the design process to output ready-to-build sequences, preempting the need for sequence redesign. The BOOST web application is available at https://boost.jgi.doe.gov and its Application Program Interfaces (API) enable integration into automated, customized DNA design processes. Finally, the herein presented results highlight the effectiveness of BOOST in reducing DNA synthesis costs and timelines.

59 BASIC BIOLOGICAL SCIENCES↗

RELAP5-3D Developer Guidelines and Programming Practices

Our ultimate goal is to create and maintain RELAP5-3D as the best software tool available to analyze nuclear power plants. This begins with writing excellent programming and requires thorough testing. This document covers development of RELAP5-3D software, the behavior of the RELAP5-3D program that must be maintained, and code testing. RELAP5-3D must perform in a manner consistent with previous code versions with backward compatibility for the sake of the users. Thus file operations, code termination, input and output must remain consistent in form and content while adding appropriate new files, input and output as new features are developed. As computer hardware, operating systems, and other software change, RELAP5-3D must adapt and maintain performance. The code must be thoroughly tested to ensure that it continues to perform robustly on the supported platforms. The coding must be written in a consistent manner that makes the program easy to read to reduce the time and cost of development, maintenance and error resolution. The programming guidelines presented her are intended to institutionalize a consistent way of writing FORTRAN code for the RELAP5-3D computer program that will minimize errors and rework. A common format and organization of program units creates a unifying look and feel to the code. This in turn increases readability and reduces time required for maintenance, development and debugging. It also aids new programmers in reading and understanding the program. Therefore, when undertaking development of the RELAP5-3D computer program, the programmer must write computer code that follows these guidelines. This set of programming guidelines creates a framework of good programming practices, such as initialization, structured programming, and vector-friendly coding. It sets out formatting rules for lines of code, such as indentation, capitalization, spacing, etc. It creates limits on program units, such as subprograms, functions, and modules. It establishes documentation guidance on internal comments. The guidelines apply to both existing and new subprograms. They are written for both FORTRAN 77 and FORTRAN 95. The guidelines are not so rigorous as to inhibit a programmer’s unique style, but do restrict the variations in acceptable coding to create sufficient commonality that new readers will find the coding in each new subroutine familiar. It is recognized that this is a “living” document and must be updated as languages, compilers, and computer hardware and software evolve.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

INTEGRATED WORKFLOW MANAGEMENT FOR PARTICLE ACCELERATOR SIMULATION

Supercomputing systems are used for a wide range of computationally demanding tasks in many fields of science and engineering. They play a key role in numerical simulation, in which mathematical models are computed in order to simulate the behavior of physical systems. Scientists and engineers that use supercomputers for numerical simulation often have their productivity limited by the need to manually organize and manage extremely large amounts of data that are often produced and consumed by the software programs run on these systems. Recognizing these limitations, Kitware Inc. (Clifton Park, NY) and SLAC National Accelerator Laboratory (Menlo Park, CA) are developing an advanced software platform that can reduce the cognitive overhead required by knowledge workers when using supercomputers for numerical simulation. Phase I of the project is complete and includes the development of new capabilities for organizing simulation project files, improvements to the user interface and overall usability, and deployment of a “middle tier” server to sit between user desktop machines and supercomputers to offload much of the data management workload. The project also developed prototype software for executing sequences of numerical simulations, and a prototype for migrating supercomputing software to cloud-based computing systems to provide a potential alternative to supercomputers with different logistical and price-to-performance tradeoffs.

Tourtellott, John↗