Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “production simulation codes”

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 613 records · Page 34

Learning the mechanisms of chemical disequilibria

When at equilibrium, large-scale systems obey thermodynamics because they have microscopic configurations that are typical. “Typical” states are a fraction of those possible with the majority of the probability. A more precise definition of typical states underlies the transmission, coding, and compression of information. However, this definition does not apply to natural systems that are transiently away from equilibrium. Here, we introduce a variational measure of typicality and apply it to atomistic simulations of a model for hydrogen oxidation. While a gaseous mixture of hydrogen and oxygen combusts, reactant molecules transform through a variety of ephemeral species en route to the product, water. Out of the exponentially growing number of possible sequences of chemical species, we find that greater than 95% of the probability concentrates in less than 1% of the possible sequences. Overall, these results extend the notion of typicality across the nonequilibrium regime and suggest that typical sequences are a route to learning mechanisms from experimental measurements. They also open up the possibility of constructing ensembles for computing the macroscopic observables of systems out of equilibrium.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Application of Hydraulic and Thermal Stimulation Techniques to Create Enhanced Geothermal Systems

Geothermal energy production requires temperature at depth, a working fluid to transport heat to the surface, and a network of fluid pathways connecting injection and production wells. These requirements constrain the availability of geothermal energy resources. Enhanced Geothermal System (EGS) development expands geothermal resource availability by artificially generating networks of fluid pathways in the reservoir. This is accomplished through hydraulic and thermal stimulation of the geothermal reservoir. Hydraulic and thermal techniques were applied at the Raft River geothermal field as part of a Department of Energy (DOE) EGS development project. The test well RRG-9 ST1 was hydraulically stimulated three times between February 2012 and April 2014. Since June, 2013 continuous injection of sub-reservoir temperature plant water has been used to thermally stimulate the well. The stimulation program has resulted in a large increase in the injectivity of the well, an increase from less than 20 gpm in June 2013 to nearly 1,000 gpm in April 2016. Wellbore imaging and temperature analysis identified a northeast striking fracture zone intersecting the well between 5,640 and 5,660 ft. measured depth (MD) has nominally accepted all of the injected fluid. Microseismic activity, related to stimulation program and plant activity was used to track possible fluid pathways within the reservoir. The microseismic data indicates that the injected fluid moves primarily to the northeast through the Narrows Zone, a northeast striking fault structure that bisects the field. Tracers injected into RRG-9 ST1 and recovered at the production wells RRG-2 and RRG-4, located to the northeast of RRG-9 ST1, support this theory. It is hypothesized that water injected into RRG-9 ST1 moves through the intersecting fracture zone which connects into the Narrows Zone, and then moves along the Narrows Zone to the northeast. This conceptual model was numerical simulated using FALCON a finite element reservoir simulation code developed by Idaho National Laboratory. The modeling suggests that both thermal and hydraulic fracturing mechanisms played an important role in improving the injectivity at RRG-9 ST1. The stimulation program at Raft River has successfully demonstrated the effectiveness of hydraulic and thermal stimulation techniques. RRG-9 ST1 is now in commercial use.

15 GEOTHERMAL ENERGY↗

Error-Bounded Learned Scientific Data Compression with Preservation of Derived Quantities

Scientific applications continue to grow and produce extremely large amounts of data, which require efficient compression algorithms for long-term storage. Compression errors in scientific applications can have a deleterious impact on downstream processing. Thus, it is crucial to preserve all the “known” Quantities of Interest (QoI) during compression. To address this issue, most existing approaches guarantee the reconstruction error of the original data or primary data (PD), but cannot directly control the problem of preserving the QoI. In this work, we propose a physics-informed compression technique that is composed of two parts: (i) reduction of the PD with bounded errors and (ii) preservation of the QoI. In the first step, we combine tensor decompositions, autoencoders, product quantizers, and error-bounded lossy compressors to bound the reconstruction error at high levels of compression. In the second step, we use constraint satisfaction post-processing followed by quantization to preserve the QoI. To illustrate the challenges of reducing the reconstruction errors of the PD and QoI, we focus on simulation data generated by a large-scale fusion code, XGC, which can produce tens of petabytes in a single day. The results show that our approach can achieve a high compression amount while accurately preserving the QoI within scientifically acceptable bounds.

97 MATHEMATICS AND COMPUTING↗

3D TRISO particle-explicit compact meshing

The TRI-structural ISOtropic (TRISO) layered fuel particle is a robust nuclear fuel form offering enhanced safety and performance for advanced reactor concepts, including high-temperature gas-cooled reactors and other Generation IV designs. These poppy-seed-sized particles are embedded in a graphite matrix to form fuel elements that must withstand elevated temperatures and high burn-up levels. The heterogeneous nature of these fuel elements — comprising thousands of randomly distributed TRISO particles — produces complex stress fields and thermal gradients that one- and two-dimensional models cannot accurately capture. While three-dimensional modeling has improved predictions of dimensional changes, internal pressure buildup, and fission product transport under irradiation, current approaches rely on homogenized material properties that are known to have considerable divergence from experimental observations. This work presents a methodology for optimized random packing of TRISO fuel compacts and full three-dimensional mesh generation within the BISON fuel performance code, with each particle coating layer individually discretized. The resulting mesh was demonstrated through heat conduction simulations under representative in-reactor operating conditions, showing strong agreement with expected behavior. This capability enables detailed analysis of particle-to-particle interactions, matrix cracking mechanisms, and the statistical distribution of coating layer failures — all of which directly govern fuel performance and safety margins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The state of Monte Carlo neutron transport on GPU accelerators - A review - 25609

A review of the current state of Monte Carlo neutron transport of Graphical Processing Units (GPUs) has been performed. Up-to-date data has been solicited from researchers in the field and combined with data available in literature to provide the current state of Monte Carlo neutron transport on GPU hardware. Nine implementations of neutron transport on GPUs have been summarized and the performance results have been presented to allow the reader to compare performance. A comparison of the results has shown that researchers obtained similar performance when scaled to the performance of NVIDIA Tesla K40 GPUs. A strong scaling equation will be presented that can be used to describe the GPU performance against different numbers of CPU cores. One implementation of a production Monte Carlo code on GPU hardware, the code Shift, showed that newer NVIDIA Pascal and Volta based GPUs provide a significant performance improvement. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

ORIGEN Integration into MPACT

This report documents the work performed to meet the objectives of the RTM milestone on ORIGEN Integration into MPACT. The main objective of the milestone was to demonstrate the capability to perform depletion problems with MPACT where the point depletion calculation, which solves the nuclide transmutation equation with ORIGEN. While MPACT already has a depletion capability based on ORIGEN, it is desired to transition to using ORIGEN because ORIGEN has undergone far more validation and software quality assurance than the implementation in MPACT. The integration was completed by compiling and linking the two codes together at build-time. Within the TriBITS build structure used for VERA, the MPACT and SCALE code repositories are already present. An optional dependency was created between MPACT to SCALE. The API available within ORIGEN was modified slightly to be able to be used by MPACT and two new software classes were implemented in MPACT to call the necessary ORIGEN subroutines. To test the ORIGEN integration two benchmark problems were analyzed. The first being the spent fuel isotopics from the Takahama-3 PWR reactor. This involved comparisons to experimentally measured nuclide concentrations for select actinides and fission products, and comparisons to other numerical results. The benchmark defines a 1-D problem limiting most of the effects due to the accuracy of a multi-dimensional transport solver, and thus differences can be assumed to be more likely caused by the point depletion solver methodology and associated nuclear data. The other benchmark analyzed two problems from the reactor physics benchmarks for LWR next generation fuel. The two configurations simulated were a PWR UO 2 pin cell and PWR UO 2 assembly. For all the benchmark problems the results were generally in good agreement with the other numerical results which were computed by several codes, including MPACT with its original depletion implementation. Thus it is concluded that ORIGEN was successfully integrated into MPACT and is capable of providing reasonable results. An additional cursory analysis of the performance of the ORIGEN point depletion compared to native MPACT point depletion solver was also performed. In this analysis it was observed that ORIGEN, on average, needs about 0.28 to 0.30 seconds to perform the point depletion calculation. This is greater than the target performance of 0.1 seconds which means future work to investigate smaller depletion chains will be needed. The MPACT internal depletion solver completes the point depletion calculation in 0.02 seconds, so it may be used for routine analysis until the use of a simplified burnup chain and performance improvements are implemented. Other future work should focus on implementing the post-correction fix-up and quadratic depletion and computation of decay heat and gamma sources.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Undergraduate education at the University of Kentucky accelerator laboratory

The University of Kentucky accelerator laboratory provides undergraduate students an opportunity to experience scientific research in nuclear physics. They are introduced to our research program, which focuses on the examination of basic nuclear properties using the 7 MV CN Van de Graaff accelerator and the time-of-flight facility to perform inelastic neutron scattering measurements with g-ray and neutron detection. Furthermore, students are given a series of projects to support the measurements including experimental equipment setup, kinematic calculations, data collection, experimental diagnostics, MCNP simulations, analysis code development, and preliminary analysis of the data. Undergraduate physics majors from the University of Dallas took part in (n,n'g) measurements for LiF. Neutron scattering cross sections for these nuclei, which are of significant interest in coolant designs for molten salt reactors, were determined for incident neutron energies between 0.8 and 4.5 MeV. By measuring the g-ray production cross sections and feeding contributions from higher-lying levels, students extracted the neutron level cross sections for Li and F. Additional (n,n'g) measurements for Li metal will be performed in the near future to supplement the g-ray data.

43 PARTICLE ACCELERATORS↗

Theory Manual for the Fuel Cycle Analysis Code REBUS

REBUS is a system of programs designed for the fuel cycle analysis of fast reactors and accelerator-driven systems. The first version of REBUS was developed in 1970 for calculating the equilibrium operating conditions of liquid metal fast breeder reactors. During the past five decades, various extensions and refinements have been made in the simulation of material motions in fuel cycle facilities and the solution of the neutronics and isotopic transmutation equations that govern the reactor behavior. The latest version complies with the standard code practices and interface dataset specifications of the committee on computer code coordination (CCCC). This report provides a detailed description of the physical processes and their mathematical models and solution methods of the current production version of REBUS. The models of the facilities and physical processes involved in fast reactor fuel cycle are first discussed. Then the mathematical formulation of these physical models is presented, followed by detailed discussions of solution methods of the resulting set of equations. Finally the output edit quantities are discussed along with their definitions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Developing an Automated Uncertainty Quantification Tool to Improve Watershed-Scale Predictions of Water and Nutrient Cycling

Managing the flow of water, nutrients, and contaminants in watersheds is vital to addressing pressing issues related to water scarcity, access to clean drinking water, energy production, resilience to natural and anthropogenic perturbations, and ecological restoration. Decisions about the management of watersheds critically depend on the accuracy with which the flow of water and chemicals through the watershed can be predicted by computer models. Prediction uncertainty can be reduced by matching the model to data, which are collected in the field at great expense. The contribution of watershed characterization data to reducing uncertainty of relevant model predictions can be evaluated in a so-called data-worth analysis, which provides transparent, quantitative metrics about a data set’s value for the support of relevant watershed management objectives. To achieve this goal, we developed a software package that implements the data-worth analysis approach for use with state-of-the-art watershed models. The purpose of the proposed data-worth analysis is to help decision-makers allocate resources for watershed characterization such that the uncertainty in model predictions can be significantly reduced, which leads to better, more effective management decisions. At the same time, watershed characterization costs can be reduced. The specific technical objectives of this SBIR/STTR Phase II project were to develop a framework and associated software toolsets that implement the uncertainty quantification and data-worth analysis approach for use with state-of-the-art watershed models. This goal was achieved by (A) developing a user-friendly, robust software package that is accessible to a wide audience, including watershed managers, policy-makers, and public stakeholders; (B) by demonstrating application of the prototype on several use cases that are representative of complex watershed management challenges spanning a range of scales and that consider different open-source, DOE-based codes and other modeling platforms; and (C) by gathering information about the needs and requirements from potential users to help guide future developments, ensuring that the final product will be commercially viable. The developed software consists of a graphical user interface that guides the user through a sequence of analysis steps, supported by toolsets that leverage state-of-the-art computational simulation-optimization capabilities. A prototype of the software runs on multiple platforms (PC, Mac, multi-processor Linux environment), is linked to diverse watershed simulators (e.g., ECOSYS, TOUGH2, TOUGHREACT, Amanzi-ATS), performs multiple analysis tasks (predictive simulations, sensitivity analysis, uncertainty analysis, automatic parameter estimation, and data-worth analysis, multicomponent geothermometry), and is readily extensible to include external simulators and analysis tools. The software is being commercialized and will be continually updated to address user needs.

58 GEOSCIENCES↗

Wind-US User's Guide, Version 2.0

Wind-US is a computational platform which may be used to numerically solve various sets of equations governing physical phenomena. Currently, the code supports the solution of the Euler and Navier-Stokes equations of fluid mechanics, along with supporting equation sets governing turbulent and chemically reacting flows. Wind-US is a product of the NPARC Alliance, a partnership between the NASA Glenn Research Center (GRC) and the Arnold Engineering Development Center (AEDC) dedicated to the establishment of a national, applications-oriented flow simulation capability. The Boeing Company has also been closely associated with the Alliance since its inception, and represents the interests of the NPARC User's Association. The "Wind-US User's Guide" describes the operation and use of Wind-US, including: a basic tutorial; the physical and numerical models that are used; the boundary conditions; monitoring convergence; the files that are read and/or written; parallel execution; and a complete list of input keywords and test options.

Towne, Charles E.↗

Wind-US Users Guide Version 4.0

Wind-US is a computational platform which may be used to numerically solve various sets of equations governing physical phenomena. Currently, the code supports the solution of the Euler and Navier-Stokes equations of fluid mechanics, along with supporting equation sets governing turbulent and chemically reacting flows. Wind-US is a product of the NPARC Alliance, a partnership between the NASA Glenn Research Center (GRC) and the Arnold Engineering Development Complex (AEDC) dedicated to the establishment of a national, applications-oriented flow simulation capability. The Boeing Company has also been closely associated with the Alliance since its inception, and represents the interests of the NPARC User's Association. The "Wind-US User's Guide" describes the operation and use of Wind-US, including: a basic tutorial; the physical and numerical models that are used; the boundary conditions; monitoring convergence; the files that are read and/or written; parallel execution; and a complete list of input keywords and test options. For current information about Wind-US and the NPARC Alliance, please see the Wind-US home page at http://www.grc.nasa.gov/WWW/winddocs/ and the NPARC Alliance home page at http://www.grc.nasa.gov/WWW/wind/.

fluid dynamics↗

CFD Simulations of the Space Launch System Ascent Aerodynamics and Booster Separation

This paper presents details of Computational Fluid Dynamic modeling of the Space Launch System during ascent. The primary focus of the paper is the flow simulation of the vehicle during ascent using the Overflow Navier-Stokes code. Computations of 739 first-stage flight conditions covering a range of Mach numbers, angles of attack, and roll angles were computed. The overset grid system contained 375 million grid points, and over 28 mil- lion CPU hours were used in the simulations. The simulations were run on the Pleiades supercomputer at the NASA Advanced Supercomputer Center at Ames Research Center. The data products from this work include integrated line-loads, surface pressure coefficients, venting pressures, and protuberance air-loads. Detailed comparisons were made of the aerodynamic performance predicted by Overflow and the wind-tunnel derived aero- dynamic database. A small number of the cases were run with two different turbulence models and with two differencing schemes. These results were used to quantify the sensitivity to the choice of the turbulence model and to the differencing scheme. The paper also introduces an effort to use the inviscid, unstructured Cartesian solver Cart3D to compute the aerodynamics during booster separation. Adaptive mesh refinement is being used to enable accurate simulations of sixteen booster-separation-motor plumes. The use of this tool is explored in preparation for building a booster-separation aerodynamic database.

Ascent Aerodynamics↗

Real-Time Digital Simulation of Inertial Response with Hardware-in-the-Loop Implementation on the CART3 Wind Turbine at the National Wind Technology Center

With increasing penetrations of wind power on electric grids, the stability and reliability of interconnected power systems may be impacted. In some countries that have developed renewable energy sources and systems, grid codes have been revised to require wind power plants (WPPs) to provide ancillary services to support the power system frequency in case of severe grid events. To do this, wind turbine generators (WTGs) should be deloaded to reserve a certain amount of active power for primary frequency response; however, deloading curtails annual energy production, and the market for this type of service needs to be further developed. In this report, we focus on the temporary frequency support provided by WTGs through inertial response. WTGs have potential to provide inertial response, but appropriate control methods should be implemented. With the implemented inertial control methods, wind turbines are capable of increasing their active power output by releasing some of their stored kinetic energy when a frequency excursion occurs. Active power can be temporarily boosted above the maximum power points, after which the rotor speed decelerates, and subsequently an active power output reduction restores the kinetic energy. In this report, we develop two types of models for wind power systems: the first is common, based on the wind power aerodynamic equation, and the power coefficient can be regressed using nonlinear functions; the second is much more complicated, wherein the wind turbine system is modeled using the Fatigue, Aerodynamics, Structures, and Turbulence Modeling (FAST) tool with several degrees of freedoms. A nine-bus test power system is built in Simulink and the Real-Time Digital Simulator, respectively, which are used to evaluate the frequency support performance of the WPPs. We implement two distinct types of inertial control methods in the modeled wind turbines: frequency-based inertial control (FBIC) and stepwise inertial control (SIC). We compare the performances of the two methods in terms of their frequency nadirs, rates of change of frequency, and recovery times. We conclude the results under various wind speeds and penetration cases, which provide insight into designing the inertial response of WTGs. Further, we discuss the impact of the parameters on the performance of the inertial control methods. We evaluate both the scaling factors for the FBIC method and the slope values for the TLIC methods. The simulation work shows the characteristics of different inertial responses compared to conventional synchronous generators. Based on the simulation results, we modify, improve, and test the inertial control methods under a more realistic wind turbine model based on FAST. We then validate the inertial responses under highly turbulent wind conditions generated by TurbSim, and we examine their influences on the turbine mechanical components. The extensive simulation proves the effectiveness of the proposed inertial control methods as well as the nine-bus test power system. We then reconsider the parameters. We rebuild the same test power system using Real time Simulator Computer Aided Design (RSCAD), and we implement the inertial control methods in the real Controls Advanced Research Turbine (CART3), which is prepared for the hardware-in-the-loop field-test simulation. After the setups for the hardware and software hybrid simulation platform are complete, the inertial response is further tested on a real wind turbine for the first time, in which CART3 release the controlled inertial response against the emulated frequency excursion, provided by the real-time simulated power system test bed in RTDS.

17 WIND ENERGY↗

LLNL Response to the DOE ASCR RFI, "Stewardship of Software for Scientific and High-Performance Computing"

For decades, Lawrence Livermore National Laboratory (LLNL) has been engaged in significant research, development, and support for software to enable scientific computing and, particularly, the use of high performance computing (HPC) in the NNSA mission space. In particular, the move in the mid-1990’s to simulation as a leading component of stockpile stewardship through the ASCI and the successor ASC programs, as well as the need for reliable data acquisition and control software for the National Ignition Facility, have been important drivers in building expertise in production-quality software development at LLNL. LLNL has also been a leader in the DOE SciDAC FASTMath Institute and the DOE Exascale Computing Project (ECP), both of which have striven to make scientific computing software – in particular, the enabling technologies underpinning simulation capabilities – more widely adopted and sustainable. As such, we believe that our experience can inform the broader goal of software stewardship for scientific and high-performance computing. LLNL strongly supports the formation of a new DOE ASCR program element in software stewardship and sustainment. Historically, DOE ASCR has funded applied mathematics and computer science research that has led to the development of important new capabilities and algorithms that are expressed as artifacts in research software. Such frameworks, libraries, and tools have seldom been directly funded to address the important issues of code maintenance, documentation, robustness, and community building. Software engineering and support have typically been done on the side in support of the ASCR-driven research products. DOE funding priorities have been slow to recognize that good software engineering, the kind that ensures research investments have more adoption and longevity, requires significant resources. Based upon our experiences, we have prepared this response to highlight the concerns and issues we believe to be important as DOE ASCR considers its role in scientific software stewardship. We believe that role is important and will require a significant investment of new funding to legitimately support the technologies past and future DOE ASCR investments have and will produce to facilitate their uptake and adoption in the broader scientific computing community. Following a summary of our involvement in scientific software development, the remainder our response is organized around the nine topics specifically identified in the RFI.

97 MATHEMATICS AND COMPUTING↗

Monte Carlo simulation of initial breakdown phase for magnetised toroidal ICRF discharges

The radio-frequency (RF) plasma production technique in the ion cyclotron range of frequency (ICRF) attracts growing attention among fusion experts because of its high potential for solving several basic problems of reactor-oriented superconducting fusion machines, such as ICRF wall conditioning in tokamaks and stellarators (T{sub e} = 3−5eV, n{sub e}<10{sup 12}cm{sup −3}), ICRF-assisted tokamak start-up and target plasma production (n{sub e} = 10{sup 13}cm{sup −3}) in stellarators. Plasma initiation by ICRF has been studied intensively using single particle descriptions and basic analytic models. To further improve the present understanding on plasma production employing the vacuum RF field of ICRF antennas in toroidal devices in presence of the toroidal magnetic field, and its parametric dependencies a Monte Carlo code has been developed. The 1D code RFdinity1D describes the motion of electrons, accelerated by the RF field in front of the ICRF antenna, along one toroidal magnetic field line. Dependent on their individual energies and the related electron collision cross sections (ionisation, excitation and dissociation) weighted by a Monte Carlo procedure, an electron avalanche may occur. Breakdown conditions are discussed as function of RF discharge parameters (i) RF vacuum electric field strength, (ii) RF frequency and (iii) neutral pressure (H2). The slope of the exponential density increase, taken as measure for the breakdown speed, shows qualitative agreement to experimental breakdown times as found in literature and experimental data of the ASDEX upgrade and TEXTOR tokamak, and is interpreted by studying the characteristic electron velocity distribution functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Adoption of High Performance Computational (HPC) Modeling Software for Widespread Use in the Manufacture of Welded Structures

This report summarizes the final product developed for the US DOE Small Business Innovation Research (SBIR) Phase II grant made to Engineering Mechanics Corporation of Columbus (Emc 2 ) between April 16, 2014 and August 31, 2016 titled ‘Adoption of High Performance Computational (HPC) Modeling Software for Widespread Use in the Manufacture of Welded Structures’. Many US companies have moved fabrication and production facilities off shore because of cheaper labor costs. A key aspect in bringing these jobs back to the US is the use of technology to render US-made fabrications more cost-efficient overall with higher quality. One significant advantage that has emerged in the US over the last two decades is the use of virtual design for fabrication of small and large structures in weld fabrication industries. Industries that use virtual design and analysis tools have reduced material part size, developed environmentally-friendly fabrication processes, improved product quality and performance, and reduced manufacturing costs. Indeed, Caterpillar Inc. (CAT), one of the partners in this effort, continues to have a large fabrication presence in the US because of the use of weld fabrication modeling to optimize fabrications by controlling weld residual stresses and distortions and improving fatigue, corrosion, and fracture performance. This report describes Emc 2 ’s DOE SBIR Phase II final results to extend an existing, state-of-the-art software code, Virtual Fabrication Technology (VFT®), currently used to design and model large welded structures prior to fabrication - to a broader range of products with widespread applications for small and medium-sized enterprises (SMEs). VFT® helps control distortion, can minimize and/or control residual stresses, control welding microstructure, and pre-determine welding parameters such as weld-sequencing, pre-bending, thermal-tensioning, etc. VFT® uses material properties, consumable properties, etc. as inputs. Through VFT®, manufacturing companies can avoid costly design changes after fabrication. This leads to the concept of joint design/fabrication where these important disciplines are intimately linked to minimize fabrication costs. Finally service performance (such as fatigue, corrosion, and fracture/damage) can be improved using this product. Emc 2 ’s DOE SBIR Phase II effort successfully adapted VFT® to perform efficiently in an HPC environment independent of commercial software on a platform to permit easy and cost effective access to the code. This provides the key for SMEs to access this sophisticated and proven methodology that is quick, accurate, cost effective and available “on-demand” to address weld-simulation and fabrication problems prior to manufacture. In addition, other organizations, such as Government agencies and large companies, may have a need for spot use of such a tool. The open source code, WARP3D, a high performance finite element code used in fracture and damage assessment of structures, was significantly modified so computational weld problems can be solved efficiently on multiple processors and threads with VFT®. The thermal solver for VFT®, based on a series of closed form solution approximations, was extensively enhanced for solution on multiple processors greatly increasing overall speed. In addition, the graphical user interface (GUI) was re-written to permit SMEs access to an HPC environment at the Ohio Super Computer Center (OSC) to integrate these solutions with WARP3D. The GUI is used to define all weld pass descriptions, number of passes, material properties, consumable properties, weld speed, etc. for the structure to be modeled. The GUI was enhanced to make it more user-friendly so that non-experts can perform weld modeling. Finally, an extensive outreach program to market this capability to fabrication companies was performed. This access will permit SMEs to perform weld modeling to improve their competitiveness at a reasonable cost.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Infrastructure for Multiphysics Software Integration in High Performance Computing-Aided Science and Engineering

The project described in this report constructed and exercised an innovative multiphysics coupling toolkit called the Illinois Rocstar MultiPhysics Application Coupling Toolkit (IMPACT). IMPACT is an open source, flexible, natively parallel infrastructure for coupling multiple uniphysics simulation codes into multiphysics computational systems. IMPACT works with codes written in several high-performance-computing (HPC) programming languages, and is designed from the beginning for HPC multiphysics code development. It is designed to be minimally invasive to the individual physics codes being integrated, and has few requirements on those physics codes for integration. The goal of IMPACT is to provide the support needed to enable coupling existing tools together in unique and innovative ways to produce powerful new multiphysics technologies without extensive modification and rewrite of the physics packages being integrated. There are three major outcomes from this project: 1) construction, testing, application, and open-source release of the IMPACT infrastructure, 2) production of example open-source multiphysics tools using IMPACT, and 3) identification and engagement of interested organizations in the tools and applications resulting from the project. This last outcome represents the incipient development of a user community and application echosystem being built using IMPACT. Multiphysics coupling standardization can only come from organizations working together to define needs and processes that span the space of necessary multiphysics outcomes, which Illinois Rocstar plans to continue driving toward. The IMPACT system, including source code, documentation, and test problems are all now available through the public gitHUB.org system to anyone interested in multiphysics code coupling. Many of the basic documents explaining use and architecture of IMPACT are also attached as appendices to this document. Online HTML documentation is available through the gitHUB site. There are over 100 unit tests provided that run through the Illinois Rocstar Application Development (IRAD) lightweight testing infrastructure that is also supplied along with IMPACT. The package as a whole provides an excellent base for developing high-quality multiphysics applications using modern software development practices. To facilitate understanding how to utilize IMPACT effectively, two multiphysics systems have been developed and are available open-source through gitHUB. The simpler of the two systems, named ElmerFoamFSI in the repository, is a multiphysics, fluid-structure-interaction (FSI) coupling of the solid mechanics package Elmer with a fluid dynamics module from OpenFOAM. This coupling illustrates how to combine software packages that are unrelated by either author or architecture and combine them into a robust, parallel multiphysics system. A more complex multiphysics tool is the Illinois Rocstar Rocstar Multiphysics code that was rebuilt during the project around IMPACT. Rocstar Multiphysics was already an HPC multiphysics tool, but now that it has been rearchitected around IMPACT, it can be readily expanded to capture new and different physics in the future. In fact, during this project, the Elmer and OpenFOAM tools were also coupled into Rocstar Multiphysics and demonstrated. The full Rocstar Multiphysics codebase is also available on gitHUB, and licensed for any organization to use as they wish. Finally, the new IMPACT product is already being used in several multiphysics code coupling projects for the Air Force, NASA and the Missile Defense Agency, and initial work on expansion of the IMPACT-enabled Rocstar Multiphysics has begun in support of a commercial company. These initiatives promise to expand the interest and reach of IMPACT and Rocstar Multiphysics, ultimately leading to the envisioned standardization and consortium of users that was one of the goals of this project.

42 ENGINEERING↗

Accelerating $\mathrm{LHC}$ event generation with simplified pilot runs and fast $\mathrm{PDF}$s

Poor computing efficiency of precision event generators for LHC physics has become a bottleneck for Monte-Carlo event simulation campaigns. We provide solutions to this problem by focusing on two major components of general-purpose event generators: The PDF evaluator and the matrix-element generator. For a typical production setup in the ATLAS experiment, we show that the two can consume about 80% of the total runtime. Using NLO simulations of pp→ℓ + ℓ - +jets and pp→$t\bar{t}$+jets as an example, we demonstrate that the computing footprint of LHAPDF and SHERPA can be reduced by factors of order 10, while maintaining the formal accuracy of the event sample. The improved codes are made publicly available.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗