Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “code optimization”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Fast multipurpose Monte Carlo simulation for proton therapy using multi- and many-core CPU architectures

Purpose: Accuracy in proton therapy treatment planning can be improved using Monte Carlo (MC) simulations. However the long computation time of such methods hinders their use in clinical routine. This work aims to develop a fast multipurpose Monte Carlo simulation tool for proton therapy using massively parallel central processing unit (CPU) architectures. Methods: A new Monte Carlo, called MCsquare (many-core Monte Carlo), has been designed and optimized for the last generation of Intel Xeon processors and Intel Xeon Phi coprocessors. These massively parallel architectures offer the flexibility and the computational power suitable to MC methods. The class-II condensed history algorithm of MCsquare provides a fast and yet accurate method of simulating heavy charged particles such as protons, deuterons, and alphas inside voxelized geometries. Hard ionizations, with energy losses above a user-specified threshold, are simulated individually while soft events are regrouped in a multiple scattering theory. Elastic and inelastic nuclear interactions are sampled from ICRU 63 differential cross sections, thereby allowing for the computation of prompt gamma emission profiles. MCsquare has been benchmarked with the GATE/GEANT4 Monte Carlo application for homogeneous and heterogeneous geometries. Results: Comparisons with GATE/GEANT4 for various geometries show deviations within 2%–1 mm. In spite of the limited memory bandwidth of the coprocessor simulation time is below 25 s for 10{sup 7} primary 200 MeV protons in average soft tissues using all Xeon Phi and CPU resources embedded in a single desktop unit. Conclusions: MCsquare exploits the flexibility of CPU architectures to provide a multipurpose MC simulation tool. Optimized code enables the use of accurate MC calculation within a reasonable computation time, adequate for clinical practice. MCsquare also simulates prompt gamma emission and can thus be used also for in vivo range verification.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

WE-DE-201-03: Combined Use of 192Ir, 60Co, and 169Yb Sources with the Novel Direction Modulated Brachytherapy Tandem Applicator for High Dose Rate Brachytherapy Planning of Cervical Cancer

Purpose: To evaluate the improvement in plan quality when various combinations of 192Ir, 60Co, and 169Yb sources are used in combination with a novel direction modulated brachytherapy (DMBT) tandem applicator for high dose rate brachytherapy of cervical cancer. Methods: The proposed DMBT tandem applicator is designed for image-guided adaptive brachytherapy (IGABT), especially MRI, of cervical cancer. It has 6 peripheral holes of 1.3-mm width, grooved along a 5.4-mm diameter nonmagnetic tungsten alloy rod of density 18.0 g/cc, capable of generating directional dose profiles - leading to enhanced dose sculpting capacity through inverse planning. Monte Carlo simulations of the three HDR sources individually inside the DMBT applicator were performed and imported into an in-house developed inverse optimization code. We then performed inverse planning with 14 cervical cancer patients enrolled in EMBRACE study. In all patients, 3D MRI-based planning was performed while utilizing 1) tandem-ring and needles attached-to-ring (7 patients) and 2) tandem-ring and needles both attached-to-ring and free-hand-loaded (7 patients), in accordance with the GEC-ESTRO recommendations. All plans were normalized to receive the same HRCTV D90 and DVH parameters were evaluated. Results: The DMBT tandem was used in all cases. Overall, the combined use of two sources (192Ir-60Co and 192Ir-169Yb, but not 60Co-169Yb) generally produced better quality plans than with the 192Ir source alone in terms of sparing OARs. For example, up to 3.5, 4.4, and 3.9% individual reductions in D2cc were observed for the bladder, rectum, and sigmoid, respectively, between 192Ir-60Co and 192Ir-only plans for patient cases in #1. While up to 5.5, 2.0, and 5.7% individual reductions were observed for patient cases in #2. Conclusion: We have demonstrated that, in addition to “directional modulation” of DMBT, use of multiple sources with sufficient differences in energy can be utilized to achieve additional improvement in plan quality for IGABT of cervical cancer.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Distributed Halide

Many image processing tasks are naturally expressed as a pipeline of small computational kernels known as stencils. Halide is a popular domain-specific language and compiler designed to implement image processing algorithms. Halide uses simple language constructs to express what to compute and a separate scheduling co-language for expressing when and where to perform the computation. This approach has demonstrated performance comparable to or better than hand-optimized code. Until now, however, Halide has been restricted to parallel shared memory execution, limiting its performance for memory-bandwidth-bound pipelines or large-scale image processing tasks. We present an extension to Halide to support distributed-memory parallel execution of complex stencil pipelines. These extensions compose with the existing scheduling constructs in Halide, allowing expression of complex computation and communication strategies. Existing Halide applications can be distributed with minimal changes, allowing programmers to explore the tradeoff between recomputation and communication with little effort. Approximately 10 new of lines code are needed even for a 200 line, 99 stage application. On nine image processing benchmarks, our extensions give up to a 1.4× speedup on a single node over regular multithreaded execution with the same number of cores, by mitigating the effects of non-uniform memory access. The distributed benchmarks achieve up to 18× speedup on a 16 node testing machine and up to 57× speedup on 64 nodes of the NERSC Cori supercomputer.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

BCL: A Cross-Platform Distributed Data Structures Library

One-sided communication is a useful paradigm for irregular parallel applications, but most one-sided programming environments, including MPI's one-sided interface and PGAS programming languages, lack application-level libraries to support these applications. We present the Berkeley Container Library, a set of generic, cross-platform, high-performance data structures for irregular applications, including queues, hash tables, Bloom filters and more. BCL is written in C++ using an internal DSL called the BCL Core that provides one-sided communication primitives such as remote get and remote put operations. The BCL Core has backends for MPI, OpenSHMEM, GASNet-EX, and UPC++, allowing BCL data structures to be used natively in programs written using any of these programming environments. Along with our internal DSL, we present the BCL ObjectContainer abstraction, which allows BCL data structures to transparently serialize complex data types while maintaining efficiency for primitive types. We also introduce the set of BCL data structures and evaluate their performance across a number of high-performance computing systems, demonstrating that BCL programs are competitive with hand-optimized code, even while hiding many of the underlying details of message aggregation, serialization, and synchronization.

Brock, Benjamin↗
DOE OSTI2011

FY2011 Status Report on SN2ND Neutronics Solver Development

A change in focus of NEAMS to a broader set of reactor types necessitated significant effort being diverted away from SN2ND development in FY2011. Nevertheless, significant gains were made with regard to building the new preconditioner and meeting the objectives of 1) parallelism in energy, 2) a spatial multigrid concept, 3) using SN vector spaces to accommodate the adjoint and time dependent functionalities needed for NEAMS work, and 4) complete manuals and documentation. While the desired goal of an exportable SN2ND solver usable within the NEAMS community for neutronics analysis was not completed, we can be confident that this year’s work puts the SN2ND solver well on its way to meeting that objective. It is important to note that SN2ND is just as capable as any other methodology for doing thermal reactor analysis, but the motivation to execute on smaller scale machines has initiated focused development on simplified treatments such as those proposed in the MOCARV tool. All of the cited coding tasks were completed but not all are fully working (i.e. debugged). The spatial multigrid technique was by far the most time consuming part of the SN2ND development process. It required multiple new data structure concepts along with identifying the necessary changes in the vectors and their associated functions. Combining the spatial multigrid as part of a GMRES preconditioner operating on the full vector space, necessary for parallelization in energy, increased the complexity of that development. Research into spatial multigrid techniques for unstructured meshes, void treatments, and code optimization and acceleration are key development needs for SN2ND. Similar to previous years, a considerable amount of effort was spent on validation work. Again we focused on fast reactor problems such as ZPR experiments and the MONJU reactor in Japan. Conventional drawer homogenized models and partially heterogeneous models were used on the ZPR calculations. The drawer homogenized models produced excellent results for both eigenvalue and foil measurements on four loadings of ZPR-6/7. The partially heterogeneous models of ZPR-6/7 also produced good results, but additional work is necessary to generate consistent effective cross section data. Calculations for MONJU were homogeneous and also gave excellent results. Finally, mock-up heterogeneous calculations were setup and attempted to scope out the performance limitations of SN2ND on whole core heterogeneous calculations. Specifically, heterogeneous models of a PWR, a VHTR, and the MONJU reactor were created using anywhere from 2 group to 33 group cross section data. Meshes with greater than 100,000,000 vertices are necessary for accuracy and a specific non-scaling memory issue in PETSc restricted the SN2ND modeling (VHTR and PWR could not be solved). All of the calculations performed this year indicate that significant computational resources beyond those available today are required to obtain a high fidelity solution possible with SN2ND. The future goals for SN2ND are to get the new version fully working, reliable in execution, and to have comparable performance to the old solver such that it can be released within the NEAMS community. After its export to RSICC, we will rely upon further preconditioner research to be done in collaboration with university partners and through the exascale hub activities. Some of that research will focus on making SN2ND applicable to thermal reactor problems, spatial multigrid, and void treatments.

Smith, M. A.↗
DOE OSTI2013

30-kW Maintenance-Free Stirling Engine for High-Performance Dish Concentrating Solar

This final report summarizes the final results of the 30-kW Maintenance-Free Stirling Engine for High-Performance Dish Concentrating Solar project performed by Infinia Corporation for the U.S. Department of Energy under contract DE-FC36-08GO18032 during the project period of March 1, 2008 – March 31, 2012. Over the course of the project, Infinia Corporation achieved substantial progress in the design, fabrication, assembly, testing and evaluation of a 30-kW multi-cylinder, free-piston Stirling engine (MCFPSE) prototype. In this effort, prior MCFPSE conceptual design approaches were adapted for the development of a 30-kW design, thereby producing a cohesive design package and detailed performance specifications. Numerical analysis using the Stirling Engine Analysis and Optimization Code (SAGE) was completed culminating in the selection of a fully balanced 6-cylinder configuration composed of two 3-cylinder engine modules. The model was refined throughout the program to include all wall conduction losses, piston and rod seal losses, appendix, shuttle and pumping losses, and piston radiation loss to closely represent the real system. Performance mapping of the engine was then completed using the detailed model. Based on knowledge and experience from its existing commercial products, and using the results from the extensive numerical modeling discussed above, Infinia successfully fabricated, assembled and tested two sets of 30-kW engine modules. One of these configurations was driven by an external power supply and operated as a cooler in order to verify the mechanical integrity of the system. In the engine configuration, induction heaters provided the energy to the heater head to operate the engine. Considerable effort was expended on correcting problems associated with contact and rubbing of internal components, identification and subsequent elimination of internal leaks, and undesirable amounts of piston bias. Many details concerning the fabrication and testing efforts are considered proprietary. As a result, those findings are not disclosed herein but rather can be found in proprietary reports supplied to the DOE on a periodic basis. The prototype did not fully demonstrate performance nor did it realize the projected cost targets. Furthermore, the basis of the mechanical and electrical architecture of the 30-kW multi-cylinder prototype design is now at least two generations old. Alternatively, Infinia’s 8-kW Bonneville generator, which was successfully demonstrated in April 2012, and is slated to enter commercial production in late 2013, offers immediate benefits of lower cost and significantly higher power density such that the generator will likely meet DOE’s Sunshot LCOE target goal of $\$$0.07-0.10/kWh by 2015 and $\$$0.05-0.07/kWh by 2020. Therefore, the design of a 30-kW Stirling engine as disclosed in the project will not be pursued further. However, Infinia remains keenly interested in developing higher power Stirling machines in the future based on the new Bonneville mechanical design.

14 SOLAR ENERGY↗
DOE OSTI2014

Final Report

The purpose of this project was to develop tools and techniques to improve the ability of computational scientists to investigate and correct problems (bugs) in their programs. Specifically, the University of Maryland component of this project focused on the problems associated with the finite number of bits available in a computer to represent numeric values. In large scale scientific computation, numbers are frequently added to and multiplied with each other billions of times. Thus even small errors due to the representation of numbers can accumulate into big errors. However, using too many bits to represent a number results in additional computation, memory, and energy costs. Thus it is critical to find the right size for numbers. This project focused on several aspects of this general problem. First, we developed a tool to look for cancelations, the catastrophic loss of precision in numbers due to the addition of two numbers whose actual values are close to each other, but whose representation in a computer is identical or nearly so. Second, we developed a suite of tools to allow programmers to identify exactly how much precision is required for each operation in their program. This tool allows programmers to both verify that enough precision is available, but more importantly find cases where extra precision could be eliminated to allow the program to use less memory, computer time, or energy. These tools use advanced binary modification techniques to allow the analysis of actual optimized code. The system, called Craft, has been applied to a number of benchmarks and real applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

FY17 Status Report on NEAMS Neutronics Activities

Under the U.S. DOE NEAMS program, the high-fidelity neutronics code system has been developed to support the multiphysics modeling and simulation capability named SHARP. The neutronics code system includes the high-fidelity neutronics code PROTEUS, the cross section library and preprocessing tools, the multigroup cross section generation code MC2-3, the in-house meshing generation tool, the perturbation and sensitivity analysis code PERSENT, and post-processing tools. The main objectives of the NEAMS neutronics activities in FY17 are to continue development of an advanced nodal solver in PROTEUS for use in nuclear reactor design and analysis projects, implement a simplified sub-channel based thermal-hydraulic (T/H) capability into PROTEUS to efficiently compute the thermal feedback, improve the performance of PROTEUS-MOCEX using numerical acceleration and code optimization, improve the cross section generation tools including MC2-3, and continue to perform verification and validation tests for PROTEUS.

47 OTHER INSTRUMENTATION↗
DOE OSTI2017

Advanced Electronic Structure Methods for Heterogeneous Catalysis and Separation of Heavy Metals

For code optimizations on Theta, graphene sheet systems were used for benchmarking purposes. Graphene sheets are of interest to researchers in the area of microlubricants. The proposed scientific applications for this project were not achievable with the support awarded for this project. The performance improvements on Theta reported in this work will directly support the scientific applications described in our recently submitted 2018 INCITE proposal.

42 ENGINEERING↗
DOE OSTI2019

Updated Groundwater Flow Simulations of the Savannah River Site General Separations Area

The groundwater flow model supporting Performance Assessments (PAs) and Composite Analyses (CAs) at the Savannah River Site was significantly revised in 2016 and 2017 using new hydrostratigraphic surfaces, updated well water level calibration targets, and semi-automated model calibration with the PEST optimization code. This model is referred to as “GSA_2016”. This report documents further refinement of the GSA_2016 model in 2018 to incorporate updates to model calibration targets, closure of the H-Area Ash Basin, construction of E-Area Slit Trench operational covers, and plume information from the Mixed Waste Management Facility and Low-Level Radioactive Waste Disposal Facility. Another objective was to lower hydraulic head residuals by adding another calibration zone. The resulting model is referred to as “GSA_2018”.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Usable Data Abstractions for Next-Generation Scientific Workflows

Data- and computationally-intensive scientific research, such as numerical simulations and inversions or the training of large neural networks in machine learning applications, that are well suited for HPC environments also often require expert insight and evaluation throughout the computation which can be greatly facilitated with the use of interactive computing tools, such as those in the Jupyter ecosystem. HPC workflows and interactive workflows are typically treated as orthogonal, however, the next generation of research will require both. The first challenge we face in this project is thus designing the right level of abstractions to allow interactive capabilities in the JupyterLab environment to allow the working scientist to flexibly explore and query their data at multiple levels, with a minimal amount of customization required of the underlying optimized codes. In addition to these questions regarding the high-level representation of data for interactive use in HPC, we tackled two additional issues that are part of the entire lifecycle of research and that become particularly acute in HPC contexts: how to improve the experience of interfacing with the HPC system's scheduling environment for a scientist focused on exploratory questions, and how can that scientist then best share the results of their work with others in a self-contained, reproducible manner.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

The Origins of Matter on Near-Exascale Supercomputing

This LDRD addressed one of the central questions in physics today: the lack of antimatter relative to matter in the universe. If antimatter and matter existed in equal amounts they would have annihilated into radiation and we would not exist today. DOE is funding experiments to investigate processes that could lead to lepton number violation and the resulting neutrino-less double beta decay (0nbb) signal that may be observed (lepton number is a property of elementary particles similar to electric charge and lepton number violation would explain the matter-antimatter imbalance). However, even if the 0nbb decay is detected in these experiments, the physics mechanism that violates lepton number would not be known, as there are several candidates. Determining the appropriate mechanism from the experimental data requires theoretical calculations. Theoretical effort in the field has focused on long-range 0nbb decay mechanisms associated with the exchange of light neutrinos between two nucleons. An equally important set of mechanisms are associated with the exchange of heavy particles and short-range physics, including the exchange of massive Majorana neutrinos. Calculating the event rate for 0nbb decays for candidate theories involving heavy particle exchange requires large-scale Lattice simulations of the underlying fundamental theory of Quantum Chromodynamics (QCD). We have performed these simulations on the LLNL supercomputers and have calculated the necessary matrix elements of short-range operators that arise from the exchange of heavy mediators that contribute to 0nbb decay via the leading pion exchange diagrams. Our published results can now be used as input to effective field theory calculations to calculate the 0nbb decay rate for the nuclei involved in the various experimental detectors. In the process of this work we entered the 2018 Gordon Bell Prize competition (the highest prize in scientific supercomputing) with our innovative and highly optimized code and methods and were selected as one of the six international finalists. During this effort we helped with the bring-up of the new LLNL Sierra and Lassen supercomputers by identifying hardware and software problems. Aiming towards the next step of this research, involving the sub-leading nucleon-nucleon direct vertex, we have developed algorithms and codes that will play a significant role in future Lattice QCD calculations involving nucleon-nucleon interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Integrated Optimization and Control of a Hybrid Gas Turbine/sCO 2 Power System

During phase-I, the project team led by Echogen Power Systems (EPS) had two primary objectives based on investigating the application of gas turbines with supercritical carbon dioxide (sCO 2 ) power cycles. The first objective was to improve the overall efficiency and performance of a hybrid gas turbine/sCO 2 power system through a joint optimization of the two subsystems (gas turbine and sCO 2 power cycle) using non-linear optimization techniques that simultaneously evaluate thermal performance of the combined cycle. The hybrid power system included several points of interaction, including (but not limited to) gas turbine exhaust, fuel heating, inlet chilling and turbine cooling. The second objective was to establish a baseline transient response model of the hybrid power system and a notional microgrid and begin steps to integrate the control systems of the three major elements (gas turbine, sCO 2 cycle and grid controller). The project team established a baseline performance for a combined cycle power plant using a production gas turbine and scaled sCO 2 power cycle only utilizing exhaust heat recovery. Echogen’s non-linear techno-economic optimization code was extended by adding gas turbine component models derived from a in-house developed gas turbine design code. With the two cycles coupled by the gas turbine exhaust, design parameters of both cycles were allowed to vary simultaneously to determine performance opportunity versus isolated designs. Returning to the baseline gas turbine/sCO 2 power cycle transient models: Echogen had in-house developed sCO 2 cycle transient model in GT-Suite system simulation software, and had partnered with Siemens Finspång for gas turbine transient model, and Siemens PTI group to provide micro-grid load profile as well as hybrid power cycle generated load (power and frequency) analysis. The transient model for the SGT-750 Siemens gas turbine was a “black-box” functional mock-up interface (FMI) model developed by Siemens Industrial Turbomachinery in Finspång, Sweden. The SGT-750 is a twin-shaft gas turbine that produces 40 MW electricity with an efficiency of about 40% at ISO conditions. At 100% gas turbine throttle (load), the SGT-750 has average exhaust conditions of 114.6 kg/s and 469.8°C. The transient model for sCO 2 power cycle was developed by Echogen in GT-SUITE 1D system simulation software platform. The basic CO 2 flow circuit has single-shaft turbomachinery with net 11.5 MW electrical power output at design conditions. The power turbine has a double-ended shaft with one end connected to synchronous generator through a fixed-ratio gearbox. The other end of power turbine is connected to the compressor through a continuously variable transmission. The major components of the sCO 2 power cycle modeled include air cooled condenser/cooler, CO 2 compressor, recuperator, two waste heat exchanger coils, power turbine, continuous variable transmission, gearbox and generator. Integration of SGT-750 transient model and sCO 2 power cycle transient model was done in Matlab Simulink. In the integrated model, the gas turbine and sCO 2 power cycle interacted at two points, first one being the gas turbine exhaust gas flow rate and temperature, which were inputs to sCO 2 power cycle model. The second point was the distribution of micro-grid load demand signal between the SGT-750 generator and sCO 2 cycle generator. For a given combined-cycle load demand, the gas turbine load demand was equal to the total demand minus the sCO 2 cycle power generated. In the present study the integrated model was simulated for two cases of grid load demand: (i) for a step change, both positive-step and negative-step, in grid load demand (ii) for a micro-grid load demand curve provided by Siemens PTI group. Finally, the time series plots representing load demand versus integrated system response were presented including the sCO 2 power cycle control system performance plots. The actual generated power and frequency of both the generators, gas turbine and sCO 2 power cycle, was supplied to Siemens PTI group for dynamic grid assessment, results of which are provided in appendices.

03 NATURAL GAS↗
DOE OSTI2021

Development of Explainable, Knowledge-Guided AI Models to Enhance the E3SM Land Model Development and Uncertainty Quantification

Focal Area(s): (2)Predictive modeling using AI techniques and AI-derived model components; use of AI and other tools to design a prediction system comprising of a hierarchy of models. (3) Insight gleaned from complex data (both observed and simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge- guided AI. Science Challenge: The Energy Exascale Earth System Model (E3SM) is a fully coupled, state-of-the-science Earth system model that uses code optimized for DOE's advanced computers to address the most critical scientific questions facing our nation and society (Golaz et al., 2019). The E3SM Land model (ELM) is designed to understand how the changes in terrestrial land surfaces will interact with other Earth system components and has been used to understand hydrologic cycles, biogeophysics, and ecosystem dynamics. In spite of great successes, the ELM has several known issues that restrain rapid improvements. For example, the ELM uses equilibrium models to simulate dynamic land-climate interactions and it requires long model spin-up time to identify suitable initial conditions for transient simulations. The ELM lacks built-in uncertainty mechanisms that can improve the robustness of model predictions. The ELM is a holistic, deterministic model system with a rigid design, and in many situations, it is hard to modify the ELM system to incorporate new theory/hypothesis and new data across scales to address emerging science problems (such as predicting the impacts of water cycle extremes). In addition, The ELM is technically optimized for traditional CPU-centric computers and it cannot fully utilize the current and incoming leadership computers for model simulations and uncertainty quantification (UQ). The success of artificial intelligence (AI) has inspired scientists to use AI models to discover intrinsic features from simulation data (Chattopadhyay et al., 2020) and observational data (Reichstein et al., 2019) to gain further process understanding of Earth science problems. However, autonomous AI model training through deep learning usually requires a huge amount of annotated data. To overcome the limitations from the data and computing resources, knowledge-guided AI models are necessary where human-knowledge is ingested in model construction (Banino et al., 2018) and training process (Silver et al., 2016) for efficient learning. Herein, we present a new way that leverages the process understanding from the ELM to guide AI model development for the ELM enhancement and UQ. We hope this study can inspire further Earth and environmental system model developments and transformations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

The US Naval Observatory VLBI Spectroscopic Catalog

Despite their importance for astrometry and navigation, the active galactic nuclei (AGNs) that comprise the International Celestial Reference Frame (ICRF) are relatively poorly understood, with key information such as their spectroscopic redshifts, AGN spectral type, and emission-/absorption-line properties generally missing from the literature. Using an updated, publicly available, state-of-the-art spectroscopic fitting code optimized for the spectra of AGNs from low to high redshift, we present a catalog of emission-line and spectral continuum parameters for 1014 unique ICRF3 objects with single-fiber spectra from the Sloan Digital Sky Survey DR16. We additionally present black hole virial mass scaling relationships that use Hα-, Hβ-, Mg II -, and C IV -based line widths, all consistent with each other, which can be used in studies of radio-loud objects across a wide range of redshifts, and we use these scaling relationships to provide derived properties such as black hole masses and bolometric luminosities for the catalog. We briefly comment on these properties for the ICRF objects, as well as their overall spectroscopic characteristics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Modular Expression Language for Ordinary Differential Equation Editing

MELODEEis a system for describing systems of initial value problem ordinary differential equations, and a compiler for the language that produces optimized code to integrate the differential equations. Features include rational polynomial approximation for expensive functions and automatic differentiation for symbolic jacobians

Blake, Robert C.↗
DOE OSTI2016

SPOC Benchmark Case: SNRE Model

The Small Nuclear Rocket Engine (SNRE) was modeled in the Center for Space Nuclear Research’s (CSNR) Space Propulsion Optimization Code (SPOC). SPOC aims to create nuclear thermal propulsion (NTP) geometries quickly to perform parametric studies on design spaces of historic and new NTP designs. The SNRE geometry was modeled in SPOC and a critical core with a reasonable amount of criticality margin was found. The fuel, tie-tubes, reflector, and control drum masses were predicted rather well. These are all very important for neutronics calculations so the active reactor geometries created with SPOC can continue to be trusted. Thermal calculations of the average and hot fuel channels agreed very well. The specific impulse calculations used historically and in SPOC disagree so mass flow rates and impulses differed. Modeling peripheral and power balance components that do not affect nuclear characteristics of the core is not a feature of SPOC and as such, these components should continue to be designed using other tools. A full paper detailing the available SNRE data and comparisons with SPOC outputs will be submitted as a follow-up to this abstract.

MCNP↗