Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “source 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,243 records · Page 69

UPC++ as_eager Working Group Draft, Revision 2020.6.2

This draft proposes an extension for a new future-based completion variant that can be more effectively streamlined for RMA and atomic access operations that happen to be satisfied at runtime using purely node-local resources. Many such operations are most efficiently performed synchronously using load/store instructions on shared-memory mappings, where the actual access may only require a few CPU instructions. In such cases we believe it’s critical to minimize the overheads imposed by the UPC++ runtime and completion queues, in order to enable efficient operation on hierarchical node hardware using shared-memory bypass. The new upcxx::{source,operation}_cx::as_eager_future() completion variant accomplishes this goal by relaxing the current restriction that future-returning access operations must return a non-ready future whose completion is deferred until a subsequent explicit invocation of user-level progress. This relaxation allows access operations that are completed synchronously to instead return a ready future, thereby avoiding most or all of the runtime costs associated with deferment of future completion and subsequent mandatory entry into the progress engine. We additionally propose to make this new as_eager_future() completion variant the new default completion for communication operations that currently default to returning a future. This should encourage use of the streamlined variant, and may provide performance improvements to some codes without source changes. A mechanism is proposed to restore the legacy behavior on-demand for codes that might happen to rely on deferred completion for correctness. Finally, we propose a new as_eager_promise() completion variant that extends analogous improvements to promise-based completion, and corresponding changes to the default behavior of as_promise().

97 MATHEMATICS AND COMPUTING↗

TU-AB-BRC-10: Modeling of Radiotherapy Linac Source Terms Using ARCHER Monte Carlo Code: Performance Comparison of GPU and MIC Computing Accelerators

Purpose: (1) To perform phase space (PS) based source modeling for Tomotherapy and Varian TrueBeam 6 MV Linacs, (2) to examine the accuracy and performance of the ARCHER Monte Carlo code on a heterogeneous computing platform with Many Integrated Core coprocessors (MIC, aka Xeon Phi) and GPUs, and (3) to explore the software micro-optimization methods. Methods: The patient-specific source of Tomotherapy and Varian TrueBeam Linacs was modeled using the PS approach. For the helical Tomotherapy case, the PS data were calculated in our previous study (Su et al. 2014 41(7) Medical Physics). For the single-view Varian TrueBeam case, we analytically derived them from the raw patient-independent PS data in IAEA’s database, partial geometry information of the jaw and MLC as well as the fluence map. The phantom was generated from DICOM images. The Monte Carlo simulation was performed by ARCHER-MIC and GPU codes, which were benchmarked against a modified parallel DPM code. Software micro-optimization was systematically conducted, and was focused on SIMD vectorization of tight for-loops and data prefetch, with the ultimate goal of increasing 512-bit register utilization and reducing memory access latency. Results: Dose calculation was performed for two clinical cases, a Tomotherapy-based prostate cancer treatment and a TrueBeam-based left breast treatment. ARCHER was verified against the DPM code. The statistical uncertainty of the dose to the PTV was less than 1%. Using double-precision, the total wall time of the multithreaded CPU code on a X5650 CPU was 339 seconds for the Tomotherapy case and 131 seconds for the TrueBeam, while on 3 5110P MICs it was reduced to 79 and 59 seconds, respectively. The single-precision GPU code on a K40 GPU took 45 seconds for the Tomotherapy dose calculation. Conclusion: We have extended ARCHER, the MIC and GPU-based Monte Carlo dose engine to Tomotherapy and Truebeam dose calculations.

60 APPLIED LIFE SCIENCES↗

Cloud Property Retrieval and 3D Radiative Transfer

Cloud thickness and photon mean-free-path together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale is observed as a change in the spatial spectrum of cloud radiances, and also as the "halo size" seen by off beam lidar such as THOR and WAIL. Such of beam lidar returns are now being used to retrieve cloud layer thickness and vertical scattering extinction profile. We illustrate with recent measurements taken at the Oklahoma ARM site, comparing these to the-dependent 3D simulations. These and other measurements sensitive to 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are providing a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval algorithms appropriate for inhomogeneous clouds. The international "Intercomparison of 3D Radiation Codes" or I3RC, program is coordinating and evaluating the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://i3rc.gsfc.nasa.gov/. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large-scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate. I3RC is currently implementing an "open source" 3d code capable of solving the baseline cases. Maintenance of this effort is one of the goals of a new 3DRT Working Group under the International Radiation Commission. It is hoped that the 3DRT WG will include active participation by land and ocean modelers as well, such as 3D vegetation modelers participating in RAMI.

Cahalan, Robert F.↗

Advanced capabilities for materials modelling with Quantum ESPRESSO

Quantum ESPRESSO is an integrated suite of open-source computer codes for quantum simulations of materials using state-of-the-art electronic-structure techniques, based on density-functional theory, density-functional perturbation theory, and many-body perturbation theory, within the plane-wave pseudopotential and projector-augmented-wave approaches. Quantum ESPRESSO owes its popularity to the wide variety of properties and processes it allows to simulate, to its performance on an increasingly broad array of hardware architectures, and to a community of researchers that rely on its capabilities as a core open-source development platform to implement their ideas. In this paper we describe recent extensions and improvements, covering new methodologies and property calculators, improved parallelization, code modularization, and extended interoperability both within the distribution and with external software.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling Z-Pinches with the FLASH code

Z-Pinches occur when large amounts of current are run axially through a cylindrical conducting material, and the magnetic pressure causes said material to implode radially to its z-axis. They are useful concepts in a large range of physics and fields like fusion. Simulating these experiments are important to not only national laboratories but also in the private sector as well. FLASH is an open source MHD code that is a useful option for modeling these experiments. However, the physics needed to be verified against an actual theory. Using a paper written by S.A Slutz, Z-Pinch experiments were run with FLASH based on the parameters described in the paper, and then results were validated against the theory presented. This verified that FLASH, although has some discrepancies in the theory due to a difference in physics calculations, remains a powerful tool for modeling Z-Pinch experiments not only on the scale of the Z-Machine from Sandia National Laboratory, but also at the scale of the proposed Next Generation Pulsed Power (NGPP) machine.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Homogeneous Cooling State as a Verification Test for Kinetic Theory-Based Continuum Models of Gas–Solid Flows

Granular and multiphase (gas–solids) kinetic theory-based models have emerged a leading modeling strategy for the simulation of particle flows. Similar to the Navier–Stokes equations of single-phase flow, although substantially more complex, kinetic theory-based continuum models are typically solved with computational fluid dynamic (CFD) codes. Under the assumptions of the so-called homogeneous cooling state (HCS), the governing equations simplify to an analytical solution describing the “cooling” of fluctuating particle velocity, or granular temperature. The HCS is used here to verify the implementation of a recent multiphase kinetic theory-based model in the open source mfix code. In this work, results from the partial verification test show that the available implicit (backward) Euler time integration scheme converges to the analytical solution with the expected first-order rate. A second-order accurate backward differentiation formula (BDF) is also implemented and observed to converge at a rate consistent with its formal accuracy.

Computational fluid dynamics↗

A code-independent generalized actuator line model for wind farm aerodynamics over simple and complex terrain

Actuator line model representations of wind turbines reduce the simulation cost of wind farm aerodynamics relative to blade geometry resolving simulations. However, most implementations are code specific and may not have proper load balancing when run in parallel. Here, a generalized actuator line model (GALM) is developed to overcome this implementation drawback of existing approaches. The GALM can be coupled to any existing microscale solver with minimal or no modifications. The coupling of our GALM model with the open-source microscale code CgWind is used to demonstrate both its ease of use and fidelity. The results obtained from the coupled-model simulations are validated for both isolated turbine in wind tunnel experiments and measurements in offshore wind farm. Finally, an operational wind farm over complex terrain is demonstrated. The simulations show that wake meandering and power production are strongly influenced by terrain impacts. Access information for obtaining the code is provided.

17 WIND ENERGY↗

Underwater Mapping Results for Seabotix vLBV300 Vehicle with Tritech Gemini 720i Imaging Sonar near Newport, OR

This document presents results from tests to demonstrate underwater mapping capabilities of an underwater vehicle in conditions typically found in marine renewable energy arrays. These tests were performed with a tethered Seabotix vLBV300 underwater vehicle. The vehicle is equipped with an inertial navigation system (INS) based on a Gladiator Landmark 40 IMU and Teledyne Explorer Doppler Velocity Log, as well as a Gemini 720i scanning sonar acquired from Tritech. The results presented include both indoor pool and offshore deployments. The indoor pool deployments were performed on October 7, 2016 and February 3, 2017 in Corvallis, OR. The offshore deployment was performed on April 20, 2016 off the coast of Newport, OR (44.678 degrees N, 124.109 degrees W). During the mission period, the sea state varied between 3 and 4, with an average significant wave height of 1.6 m. Data was recorded from both the INS and the sonar.

16 TIDAL AND WAVE POWER↗

Autocoding for Verifiability

A viewgraph presentation of an automatic code scheme for source verification issues is shown. The topics include: 1) Model-Based Control Law Development with Automatic Code Generation; 2) Source Verification Issues of Automatically Generated Code; 3) MXZ Code Generator; 4) Benefits and Limitations of MXZ; and 5) Current State of the Tool.

Smith, Tim↗

A Wideband Satcom Based Avionics Network with CDMA Uplink and TDM Downlink

The purpose of this paper is to describe some key technical ideas behind our vision of a future satcom based digital communication network for avionics applications The key features of our design are as follows: (a) Packetized transmission to permit efficient use of system resources for multimedia traffic; (b) A time division multiplexed (TDM) satellite downlink whose physical layer is designed to operate the satellite link at maximum power efficiency. We show how powerful turbo codes (invented originally for linear modulation) can be used with nonlinear constant envelope modulation, thus permitting the satellite amplifier to operate in a power efficient nonlinear regime; (c) A code division multiple access (CDMA) satellite uplink, which permits efficient access to the satellite from multiple asynchronous users. Closed loop power control is difficult for bursty packetized traffic, especially given the large round trip delay to the satellite. We show how adaptive interference suppression techniques can be used to deal with the ensuing near-far problem; (d) Joint source-channel coding techniques are required both at the physical and the data transport layer to optimize the end-to-end performance. We describe a novel approach to multiple description image encoding at the data transport layer in this paper.

Agrawal, D.↗

Virtual Flow Solver - Geophysics: A 3D Incompressible Navier-Stokes Solver

Virtual Flow Solver - Geophysics (VFS-Geophysics) is a three-dimensional (3D) incompressible Navier-Stokes solver based on the Curvilinear Immersed Boundary (CURVIB) method. The CURVIB is a sharp interface type of immersed boundary (IB) method that enables the simulation of fluid flows in the presence of geometrically complex moving bodies. The CURVIB method can be applied to wind/MHK turbine simulations and energy applications. VFS-Geophysics is the result of many years of research work by several graduate students, post-docs, and research associates that have been involved in the Computational Hydrodynamics and Biofluids Laboratory directed by Professor Fotis Sotiropoulos. The preparation of the present manual has been supported by the U.S. Department of Energy (DE-EE 0005482).

16 TIDAL AND WAVE POWER↗

A Comparison of GPU-Accelerated Multiphase CFD Solvers on the Polaris Supercomputer: Part 1

This report is in support of the Innovative and Novel Computational Impact on Theory and Experiment (INCITE) program sponsored by the U.S. Department of Energy (USDOE). With INCITE-level resources, one project, titled BubblyFlow, was granted computational resources for the 2025 calendar year on the Polaris supercomputer at the Argonne Leadership Computing Facility (ALCF). The project aims to conduct simulations to understand the fundamental characteristics of turbulent bubbly flow phenomena in nature. Staff at the ALCF and Argonne’s Computational Science division, along with collaborators at the City College of New York and University of Illinois at Chicago, helped a summer student to assess the accuracy and performance of two high performance computing (HPC) codes. Both codes, ImExLBM and FluTAS, are fundamentally different in their mathematical and numerical modeling. However, both may be used to solve the same physical problem. The collaboration sought to better understand the differences between both codes in terms of accuracy and efficiency. This would ultimately help the BubblyFlow project better utilize resources and establish a knowledge-base of code capabilities in future simulation campaigns. We compare ImExLBM and FluTAS, two high-performance multiphase computational fluid dynamics (CFD) solvers, in terms of physical fidelity, time-to-solution, and parallel efficiency. We validate ImExLBM (Implicit-Explicit Lattice Boltzmann Method) against a canonical benchmark and assess it’s performance relative to FluTAS (Fluid Transport Accelerated Solver), a well-established open-source CFD code.

97 MATHEMATICS AND COMPUTING↗

High performance jet-engine flight test data base for HSR

The primary acoustic priority of the flight test data base for HSR is the validation of the NASA Aircraft Noise Prediction Program (ANOPP) and other source noise codes. Also, the noise measurements are an important support function for the High Lift Program devoted to HSR. Another concern that will be addressed is a possible noise problem 7-20 miles from take-off during climbout. The attention arises from the higher speeds envisioned for the HSCT compared to conventional aircraft causing levels to increase because of Doppler amplification in conjunction with high source levels due to jet noise. An attempt may be made to measure airframe noise for the F-16XL test which would provide an assessment of this noise component for delta wing aircraft.

Kelly, Jeffrey↗

Hybrid-VPIC: An open-source kinetic/fluid hybrid particle-in-cell code

Hybrid-VPIC is an extension of the open-source high-performance particle-in-cell (PIC) code VPIC incorporating hybrid kinetic ion/fluid electron solvers. This paper describes the models that are available in the code and gives an overview of applications of the code to space and laboratory plasma physics problems. Particular choices in how the hybrid solvers were implemented are documented for reference by users. A few solutions for handling numerical complications particular to hybrid codes are also described. Finally, special emphasis is given to the computationally taxing problem of modeling mix in collisional high-energy-density regimes, for which more accurate electron fluid transport coefficients have been implemented for the first time in a hybrid PIC code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Practical Guide to Writing a Radiative Transfer Code

Using our decades-long experience in radiative transfer (RT) code development for Earth science, we endeavor to reduce the knowledge gap of bringing RT from theory to code quickly. Despite numerous classic and recent literature, it is still hard to develop anRT code from scratch within a few weeks. It is equally hard to understand, not to mention modify, an existing “monster” RT code, for which the developer is either located remotely or has retired. Following the format of “Numerical Recipes” by Press et al., we collocate in this paper small pieces of necessary theory with corresponding small pieces of RT code. These are arranged in an order that is natural for code development, which is often opposite of the natural order for laying out the theoretical basis. We focus on the transfer of unpolarized monochromatic solar radiation in a plane-parallel atmosphere over a reflecting surface. Both the surface and the atmosphere are homogeneous (uniform) at all directions. The multiple scattering is numerically solved using the deterministic method of Gauss-Seidel iterations. Except for the presented Python-Numba open-source RT code gsit, the paper does not report any new scientific results, but rather serves as an academic demonstration. If development time is an issue or the reader is familiar with basic concepts of RT theory, we recommend proceeding directly to Sec.3 “RT code development.

multiple light scattering↗

Development and Validation of Physics-Based Sub-Models of High-Pressure Supercritical Fuel Injection at Diesel Conditions Insert

The project sought to develop and validate accurate, physics-based, numerical submodels and implement them in computational fluid dynamics (CFD) software codes to enable accurate predictions of diesel sprays at pressures higher than the current systems. The combined experimental-computational effort has four main objectives: (1) acquire spatially and temporally resolved scalar measurements of fuel-oxidizer mixing in the near field of the fuel jet, (2) develop software modules to calculate thermo-physical properties of diesel surrogates and their mixtures with oxidizer of interest in real-gas operating regime, (3) integrate the real-fluid model into a commercial and an open-source CFD code to simulate fuel-air mixing in high-pressure diesel sprays, and (4) assess the robustness, accuracy, and uncertainty of the integrated CFD solvers in computing diesel sprays. Project accomplishments are consistent with the original goals of the project. Although we acquired spatially and temporally resolved scalar measurements, they were limited to refractive index difference and not to thermodynamic properties such as temperature and equivalence ratio. However, we successfully developed the theoretical background needed to acquire such measurements at diesel conditions and demonstrated the capability of the RSD for scalar measurements in a canonical turbulent jet flow replicating features of the diesel jet. Additional work will be needed to implement these procedures at high-pressure, high-temperature engine-like conditions of the CPFR. Mainly, we did not anticipate the need to find the liquid boundary and its statistical features or the state relationships between refractive index and thermodynamic properties, assuming the widely known Dale-Gladstone equation would suffice. Still, efforts devoted to these tasks were essential to fulfil the goals of the project and serves as the foundation on which many of the conclusions and key contributions of this project are based. Further assessment and application of the developed CFD codes is recommended to fully realize the benefits of this research.

33 ADVANCED PROPULSION SYSTEMS↗

GJ 229B: Solving the Puzzle of the First Known T Dwarf with the APOLLO Retrieval Code

GJ 229B was the first T dwarf to be discovered in 1995, and its spectrum has been more thoroughly observed than that of most other brown dwarfs. Yet a full spectroscopic analysis of its atmosphere has not been done with modern techniques. This spectrum has several peculiar features, and recent dynamical estimates of GJ 229B’s mass and orbit have disagreed widely, both of which warrant closer investigation. With a separation of tens of astronomical units from its host star, GJ 229B falls near the border of the planet and stellar population formation regimes, so its atmosphere could provide clues to formation processes for intermediate objects of this type. In an effort to resolve these questions, we performed retrievals on published spectra of GJ 229B over a wide range of wavelengths (0.5–5.1 μm) using the open-source APOLLO code. Based on these retrievals, we present a more precise mass estimate of 41.6 ± 3.3MJ and an effective temperature estimate of 869 +5 −7 K, which are more consistent with evolutionary models for brown dwarfs and suggest an older age for the system of >1.0 Gyr. We also present retrieved molecular abundances for the atmosphere, including replicating the previously observed high CO abundance, and discuss their implications for the formation and evolution of this object. This retrieval effort will give us insight into how to study other brown dwarfs and directly imaged planets, including those observed with JWST and other next-generation telescopes.

Alex R Howe↗