Search NASA⌕ Search

SEARCH · Search NASA

Results for “System Code”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,783 records · Page 99

THE MULTI-COLOR LIGHT CURVES OF THE W UMa TYPE CONTACT BINARY EP ANDROMEDAE

New multi-color CCD photometric light curves of the eclipsing binary EP And were obtained over six nights in 2006, 2011, and 2012. Using the Wilson-Devinney code, we computed the photometric elements of this system. It was discovered that EP And is a W-type W UMa contact binary system with a mass ratio of q = 2.685 and a degree of contact factor f = 24.9%, rather than an A-type system. Combining 28 newly determined times of minimum light derived from 1999 to 2012 with others collected from the literature, a long-term increase (dP/dt = +5.22 Multiplication-Sign 10{sup -8} days yr{sup -1}) with a sinusoidal variation (A = 0.0109 days; T = 40.89 yr) in the orbital period was found. The orbital period secular increases may be interpreted as conservative mass transfer from the less massive component to the more massive one, and cyclic variations of the orbital period may be caused by the light-travel time effect through the presence of a third body. The evolutionary status and the age of the system are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Mixed Stochastic-Deterministic Approach for Many-Body Perturbation Theory Calculations

Here, we present an approach for GW calculations of quasiparticle energies with quasi-quadratic scaling by approximating high-energy contributions to the Green's function in its Lehmann representation with effective stochastic vectors. The method is easy to implement without altering the GW code, converges rapidly with stochastic parameters, and treats systems of various dimensionality and screening response. Our calculations on a 5.75° twisted MoS 2 bilayer show how large-scale GW methods include geometry relaxations and electronic correlations on an equal basis in structurally nontrivial materials.

36 MATERIALS SCIENCE↗

PyOMP: Multithreaded Parallel Programming in Python

We know that Python is a widely used language in scientific computing. When the goal is high performance, however, Python lags far behind low-level languages such as C and Fortran. To support applications that stress performance, Python needs to access the full capabilities of modern CPUs. That means support for parallel multithreading. In this paper, we describe PyOMP, a system that enables OpenMP in Python. Programmers write code in Python with OpenMP, Numba generates code that compiles to LLVM, and the resulting programs run with performance that approaches that from code written with C and OpenMP. In this paper we provide an update on the PyOMP project and explain how to install it and use it to write parallel multithreaded code in Python.

97 MATHEMATICS AND COMPUTING↗

GPU Lifetimes on Titan Supercomputer: Survival Analysis and Reliability

The Cray XK7 Titan was the top supercomputer system in the world for a long time and remained critically important throughout its nearly seven year life. It was an interesting machine from a reliability viewpoint as most of its power came from 18,688 GPUs whose operation was forced to execute three rework cycles, two on the GPU mechanical assembly and one on the GPU circuitboards. We write about the last rework cycle and a reliability analysis of over 100,000 years of GPU lifetimes during Titan’s 6-year-long productive period. Using time between failures analysis and statistical survival analysis techniques, we find that GPU reliability is dependent on heat dissipation to an extent that strongly correlates with detailed nuances of the cooling architecture and job scheduling. We describe the history, data collection, cleaning, and analysis and give recommendations for future supercomputing systems. We make the data and our analysis codes publicly available.

Ostrouchov, George↗

STEPS: A Portable Numerical Simulation Toolkit for Electrical Power System Dynamic Studies

Numerical simulation is the key technique for large scale power system analysis. Redistribution of global renewable power via international interconnections requires new simulation tools to study the interconnected systems with different nominal frequencies as a whole. In this paper we introduce an open source simulation toolkit for electrical power systems (STEPS) which is hosted at Github. Its kernel is coded in C++ with major functions of power flow and electro-mechanical dynamic simulation. Flexible options are provided and configurable to improve power flow solution and dynamic simulation. Common devices and models are supported in STEPS for AC/DC hybrid system studies. Studies of interconnected systems with different nominal frequencies is supported in STEPS for research of international interconnection. Application program interfaces are provided and wrapped with Python to enable high-level interfaces for general applications. STEPS is thread safe and parallel computation is supported in both kernel and script levels to accelerate simulation. It is portable and works on Windows and GNU/Linux platforms. Cases from small to large scale systems are thoroughly tested to validate the toolkit with commercial packages as benchmarks.

42 ENGINEERING↗

Endoluminal ultrasound applicators for MR-guided thermal ablation of pancreatic tumors: Preliminary design and evaluation in a porcine pancreas model

Purpose: Endoluminal ultrasound may serve as a minimally invasive option for delivering thermal ablation to pancreatic tumors adjacent to the stomach or duodenum. The objective of this study was to explore the basic feasibility of this treatment strategy through the design, characterization, and evaluation of proof-of-concept endoluminal ultrasound applicators capable of placement in the gastrointestinal (GI) lumen for volumetric pancreas ablation under MR guidance. Methods: Two variants of the endoluminal applicator, each containing a distinct array of two independently powered transducers (10 × 10 mm 3.2 MHz planar; or 8 × 10 × 20 mm radius of curvature 3.3 MHz curvilinear geometries) at the distal end of a meter long flexible catheter assembly, were designed and fabricated. Transducers and circulatory water flow for acoustic coupling and luminal cooling were contained by a low-profile polyester balloon covering the transducer assembly fixture. Each applicator incorporated miniature spiral MR coils and mechanical features (guiding tips and hinges) to facilitate tracking and insertion through the GI tract under MRI guidance. Acoustic characterization of each device was performed using radiation force balance and hydrophone measurements. Device delivery into the upper GI tract, adjacent to the pancreas, and heating characteristics for treatment of pancreatic tissue were evaluated in MR-guided ex vivo and in vivo porcine experiments. MR guidance was utilized for anatomical target identification, tracking/positioning of the applicator, and MR temperature imaging (MRTI) for PRF-based multislice thermometry, implemented in the real-time RTHawk software environment. Results: Force balance and hydrophone measurements indicated efficiencies of 48.8% and 47.8% and −3 dB intensity beam-widths of 3.2 and 1.2 mm for the planar and curvilinear transducers, respectively. Ex vivo studies on whole-porcine carcasses revealed capabilities of producing ablative temperature rise (ΔT > 15 °C) contours in pancreatic tissue 4–40 mm long and 4–28 mm wide for the planar transducer applicator (1–13 min sonication duration, ∼4 W/cm{sup 2} applied acoustic intensity). Curvilinear transducers produced more selective heating, with a narrower ΔT > 15 °C contour length and width of up to 1–24 mm and 2–7 mm, respectively (1–7 min sonication duration, ∼4 W/cm{sup 2} applied acoustic intensity). Active tracking of the miniature spiral coils was achieved using a Hadamard encoding tracking sequence, enabling real-time determination of each coil’s coordinates and automated prescription of imaging planes for thermometry. In vivo MRTI-guided heating trials in three pigs demonstrated capability of ∼20 °C temperature elevation in pancreatic tissue at 2 cm depths from the applicator, with 5–7 W/cm{sup 2} applied intensity and 6–16 min sonication duration. Dimensions of thermal lesions in the pancreas ranged from 12 to 28 mm, 3 to 10 mm, and 5 to 10 mm in length, width, and depth, respectively, as verified through histological analysis of tissue sections. Multiple-baseline reconstruction and respiratory-gated acquisition were demonstrated to be effective strategies in suppressing motion artifacts for clear evolution of temperature profiles during MRTI in the in vivo studies. Conclusions: This study demonstrates the technical feasibility of generating volumetric ablation in pancreatic tissue using endoluminal ultrasound applicators positioned in the stomach lumen. MR guidance facilitates target identification, device tracking/positioning, and treatment monitoring through real-time multislice PRF-based thermometry.

60 APPLIED LIFE SCIENCES↗

An Efficient Multicore Implementation of a Novel HSS-Structured Multifrontal Solver Using Randomized Sampling

Here, we present a sparse linear system solver that is based on a multifrontal variant of Gaussian elimination and exploits low-rank approximation of the resulting dense frontal matrices. We use hierarchically semiseparable (HSS) matrices, which have low-rank off-diagonal blocks, to approximate the frontal matrices. For HSS matrix construction, a randomized sampling algorithm is used together with interpolative decompositions. The combination of the randomized compression with a fast ULV HSS factoriz ation leads to a solver with lower computational complexity than the standard multifrontal method for many applications, resulting in speedups up to 7 fold for problems in our test suite. The implementation targets many-core systems by using task parallelism with dynamic runtime scheduling. Numerical experiments show performance improvements over state-of-the-art sparse direct solvers. The implementation achieves high performance and good scalability on a range of modern shared memory parallel systems, including the Intel Xeon Phi (MIC). The code is part of a software package called STRUMPACK - STRUctured Matrices PACKage, which also has a distributed memory component for dense rank-structured matrices.

97 MATHEMATICS AND COMPUTING↗

High-Fidelity Numerical Investigation on Elucidating Sodium Heat Transfer Characteristics for 37-Pin Wire-Wrapped Fuel Bundle in the PLANDTL Facility

This study involved a Reynolds-averaged Navier-Stokes- (RANS-) based computational fluid dynamics (CFD) analysis of the 37-pin wire-wrapped fuel bundle of the PNC Plant dynamics test loop (PLANDTL) facility. Previously, mainly the hydrodynamic phenomena of the wire-wrapped fuel bundle were analyzed, but the present study additionally included heat transfer analysis through conjugate heat transfer. The main purpose of the study was to benchmark the experimental data of the PLANDTL 37-pin wire-wrapped fuel bundle to investigate the heat transfer phenomena. In addition, the aim was to verify the accuracy of the RANS-based CFD analysis method using the STAR-CCM+ simulation software in comparison with the experimental data. The grid used for verification was an innovative grid system consisting of hexahedra using Fortran-based code. The development of the RANS-based CFD methodology included grid sensitivity analysis, turbulence model sensitivity analysis, and turbulent Prandtl number sensitivity analysis. Information on the temperature, mass flow rate, and area of the CFD results for each subchannel was provided for the top of the heated section and is expected to serve as a reference for future studies aiming to perform the validation and verification of a PLANDTL facility. In addition, the dependence of the peak temperature on the azimuth angle of each pin was analyzed.

97 MATHEMATICS AND COMPUTING↗

Stochastic Hybrid Models and Algorithms for Fluids (FHDeX) v1

This software contains numerical modules and application codes for studying the behavior of fluid dynamical systems where fluctuations play an important role in the dynamics. There is support for finite-volume/finite-difference solvers for partial differential equations, particle methods, and hybrid fluid-particle approaches.

Bell, JohnB.↗

Screamer version 4.1

Screamer is a special purpose circuit code developed for the design of Pulsed Power systems. It models electrical circuits which have a restricted topology in order to provide a fast-running tool while still allowing configurations general enough for most Pulsed Power system designs

Spielman, Rick↗

Icypaw

The Intermode Communication Protocol Wrapper, or Icypaw, structures and facilitates communication between endpoints in a network. It is used to connect various independent software and hardware components into one loosely organized system controllable from one or more locations. The code uses the MQTT protocol to communicate and enhances it with a useful framework. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525. SAND2021-8293 O

Van Der Wall, JayW.↗

Structural Code Considerations for Solar Rooftop Installations

Residential rooftop solar panel installations are limited in part by the high cost of structural related code requirements for field installation. Permitting solar installations is difficult because there is a belief among residential permitting authorities that typical residential rooftops may be structurally inadequate to support the additional load associated with a photovoltaic (PV) solar installation. Typical engineering methods utilized to calculate stresses on a roof structure involve simplifying assumptions that render a complex non-linear structure to a basic determinate beam. This method of analysis neglects the composite action of the entire roof structure, yielding a conservative analysis based on a rafter or top chord of a truss. Consequently, the analysis can result in an overly conservative structural analysis. A literature review was conducted to gain a better understanding of the conservative nature of the regulations and codes governing residential construction and the associated structural system calculations.

Dwyer, Stephen F.↗

Numerical Estimation of the Spent Fuel Ratio

Sabotage of spent nuclear fuel casks remains a concern nearly forty years after attacks against shipment casks were first analyzed and has a renewed relevance in the post-9/11 environment. A limited number of full-scale tests and supporting efforts using surrogate materials, typically depleted uranium dioxide (DUO 2 ), have been conducted in the interim to more definitively determine the source term from these postulated events. However, the validity of these large- scale results remain in question due to the lack of a defensible spent fuel ratio (SFR), defined as the amount of respirable aerosol generated by an attack on a mass of spent fuel compared to that of an otherwise identical surrogate. Previous attempts to define the SFR in the 1980's have resulted in estimates ranging from 0.42 to 12 and include suboptimal experimental techniques and data comparisons. Because of the large uncertainty surrounding the SFR, estimates of releases from security-related events may be unnecessarily conservative. Credible arguments exist that the SFR does not exceed a value of unity. A defensible determination of the SFR in this lower range would greatly reduce the calculated risk associated with the transport and storage of spent nuclear fuel in dry cask systems. In the present work, the shock physics codes CTH and ALE3D were used to simulate spent nuclear fuel (SNF) and DUO 2 targets impacted by a high-velocity jet at an ambient temperature condition. These preliminary results are used to illustrate an approach to estimate the respirable release fraction for each type of material and ultimately, an estimate of the SFR. This page intentionally blank

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Verification of the DIF3D Software to Support Fast Reactor Analysis

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the polynomial order used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Terry Turbopump Expanded Operating Band Modeling and Simulation Efforts in Fiscal Year 2020 - Progress Report

The Terry Turbine Expanded Operating Band Project is currently conducting testing at Texas A&M University as part of a revised experimental program meant to supplant previous full-scale testing plans under the headings of Milestone 5 and Milestone 6. In consultation with Sandia National Laboratories technical staff and with modeling and simulation support from the same, the hybrid Milestone 5&6 plan is moving forward with experiments aimed at addressing knowledge gaps regarding scale, working fluid, and turbopump self-regulation. Modeling and simulation efforts at Sandia National Laboratories in FY20 fell under the broad umbrella of Milestone 7 and consisted exclusively of MELCOR-related tasks aimed at: 1) Constructing/improving input models of Texas A&M University experiments, 2) Constructing a generic boiling water reactor input model according to best practices with systems-level Teny turbine capabilities, and 3) Adding code capability in order to leverage experimental data/findings, address bugs, and improve general code robustness Project impacts of the Covid-19 pandemic have fortunately been minimal thus far but are mentioned as necessary when discussing the hybrid Milestone 5&6 progress as well as the corresponding Milestone 7 modeling and simulation progress.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Data-Driven Exploration of Climate Attractor Manifolds For Long-Term Predictability

Focal Area: This white paper responds to Focal Area 3. We seek to gain insight into decadal-scale climate predictability by applying novel manifold-finding probabilistic AI techniques to the complex data produced by Earth System models (ESMs) such as E3SM. The associated portfolio of research activities leverages DOE’s asset mix of HPC platforms, climate expertise, climate simulation codes, and AI expertise. Science Challenge: Climate and climate models are dynamical systems exhibiting properties that are interpretable through chaos theory. The theory contains an important concept that is relevant to multi-decade-scale climate prediction: a chaotic attractor. While the space containing all the possible states of the Earth’s atmosphere and ocean, the possible weather, is large, the realized states tend to stay near the smaller-dimensioned attractor. This behavior is responsible for the “order behind the irregularity” [1] of climate phenomena. Climate change can be thought of as a change in the properties of the attractor, and predicting the climate over years to decades is equivalent to predicting how those properties will change. To date, the attractor has been a useful conceptual tool, but has not been amenable to direct characterization. A new development is the advent of efficient high-dimensional manifold-finding probabilistic AI techniques, which permit a data-driven characterization of the ESM attractor and its probability distribution over weather states. Such a characterization would result in a natural dimensional reduction — a “non-linear Principal Components Analysis (PCA) adapted to climate simulation data” — leading to important advances in scenario-based long-term climate prediction, long-term prediction of water cycle extremes, ESM verification, inter-model comparison, and process model development.

54 ENVIRONMENTAL SCIENCES↗

Instabilities driven by the fundamental crabbing mode

The Electron Ion Collider (EIC) requires crab cavities to correct for crossing angle effects. These cavities produce a time varying horizontal kick that makes the collision head on in the center of mass frame of the bunches. This results in a luminosity that is very close to the the luminosity obtained without a crossing angle. Producing the kick requires a transverse RF field. This field requires a superconducting resonator and reasonable power levels require loaded quality factors of order a million. This results is a very large transverse impedance for the crabbing mode. The horizontal beta function is large and this, coupled with the large transverse impedance of the resonator can lead to strong instabilities. In the first section of the paper a simple model is used to describe the situation without feedback. There are large ranges of tune leading to strong instabilities in both the electron and hadron storage rings of the EIC. These estimates are confirmed with simulations. In the second section, RF feedback is used to reduce the apparent impedance of the crabbing mode. Simulating the actual feedback would require both a significant rewrite of the tracking code and a significant reduction in computational speed. Instead, the apparent impedance is fitted using one pole filters, which allows use of the existing simulation code. In the third section a transverse damping system is outlined.

43 PARTICLE ACCELERATORS↗