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.

SourceResearch connectionFrom yearThrough year

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

At least 991 records · Page 55

All sources NASA NTRS DOE OSTI

DOE OSTI2016

SU-G-JeP3-10: Update On a Real-Time Treatment Guidance System Using An IR Navigation System for Pleural PDT

Purpose: Photodynamic therapy (PDT) is used in conjunction with surgical debulking of tumorous tissue during treatment for pleural mesothelioma. One of the key components of effective PDT is uniform light distribution. Currently, light is monitored with 8 isotropic light detectors that are placed at specific locations inside the pleural cavity. A tracking system with real-time feedback software can be utilized to improve the uniformity of light in addition to the existing detectors. Methods: An infrared (IR) tracking camera is used to monitor the movement of the light source. The same system determines the pleural geometry of the treatment area. Software upgrades allow visualization of the pleural cavity as a two-dimensional volume. The treatment delivery wand was upgraded for ease of light delivery while incorporating the IR system. Isotropic detector locations are also displayed. Data from the tracking system is used to calculate the light fluence rate delivered. This data is also compared with in vivo data collected via the isotropic detectors. Furthermore, treatment volume information will be used to form light dose volume histograms of the pleural cavity. Results: In a phantom study, the light distribution was improved by using real-time guidance compared to the distribution when using detectors without guidance. With the tracking system, 2D data can be collected regarding light fluence rather than just the 8 discrete locations inside the pleural cavity. Light fluence distribution on the entire cavity can be calculated at every time in the treatment. Conclusion: The IR camera has been used successfully during pleural PDT patient treatment to track the motion of the light source and provide real-time display of 2D light fluence. It is possible to use the feedback system to deliver a more uniform dose of light throughout the pleural cavity.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2023

PeleLMeX: an AMR Low Mach Number Reactive Flow Simulation Code without level sub-cycling

PeleLMeX simulates chemically reacting low Mach number flows with block-structured adaptive mesh refinement (AMR). The code is built upon the AMReX library, which provides the underlying data structures and tools to manage and operate on them across massively parallel computing architectures. PeleLMeX algorithmic features are inherited from its predecessor PeleLM but key improvements allow representation of more complex physical processes. Together with its compressible flow counterpart PeleC, the thermo-chemistry library PelePhysics and the multi-physics library PeleMP, it forms the Pele suite of open-source reactive flow simulation codes.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Covariance Structure of Pu Isotopic Compositions Obtained Using FRAM 5.2

In the nondestructive assay of special nuclear materials, isotopic determination based on the analysis of high-purity Ge high-resolution gamma-ray spectra plays a pivotal role. The isotopic results, which are based on a relative measurement of peak intensities, are used to interpret data acquired using quantitative (mass) nondestructive assay methods. In the case of plutonium measurements, this includes both passive and active neutron techniques, and also calorimetry (heat). Historically, precision estimates on the Pu- isotopic composition have been generated, on a case-by-case basis, by least squares fitting algorithms embedded in the analysis codes, but all other sources of uncertainty have generally been estimated using a top-down approach by item category. A largely neglected consideration is the covariance between the reported abundances and derived quantities – notably Pu-239-effective mass, Pu-240-effective mass and specific thermal power. In this work we study the covariance structure by statistically evaluating repeat spectrum measurements of two known low burn-up heavily filtered Pu items reduced using the spectrum analysis code FRAM 5.2 and alternate standard parameter sets. For routine measurements, however, only a single spectrum from an item will typically be acquired. We therefore discuss how to generate the sought after covariance information from an individual spectrum by spectrum re-sampling, an approach which could be added to provide item specific reports, and compare the two approaches to test its reliability.

Croft, Stephen↗

DOE OSTI2016

DNS Study of Supercritical CO{sub 2} Flow Using Spectral Element Code Nek5000

The study of supercritical CO{sub 2} (s-CO{sub 2}) heat transfer has become popular in recent decades. The driving factor is the growing interest of using s-CO{sub 2} in future power conversion cycles. The supercritical CO{sub 2} Brayton cycle has gained the most focus for its high efficiency, compact design, economics, etc, . In this cycle, heat exchangers are used to heat up or cool down CO{sub 2}. Thus, the heat transfer characteristics of s-CO{sub 2} need to be well understood for heat exchanger design. In this research, we use the open-source spectral-element code Nek5000, developed by Argonne National Laboratory, to perform a direct numerical simulation (DNS) for s-CO{sub 2} flow through a straight heated tube. The pioneering paper by Bae et al. is used to validate the results from the Nek5000 simulation. (authors)

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2013

THE SUPERNOVA THAT DESTROYED A PROTOGALAXY: PROMPT CHEMICAL ENRICHMENT AND SUPERMASSIVE BLACK HOLE GROWTH

The first primitive galaxies formed from accretion and mergers by z {approx} 15, and were primarily responsible for cosmological reionization and the chemical enrichment of the early cosmos. But a few of these galaxies may have formed in the presence of strong Lyman-Werner UV fluxes that sterilized them of H{sub 2}, preventing them from forming stars or expelling heavy elements into the intergalactic medium prior to assembly. At masses of 10{sup 8} M{sub Sun} and virial temperatures of 10{sup 4} K, these halos began to rapidly cool by atomic lines, perhaps forming 10{sup 4}-10{sup 6} M{sub Sun} Pop III stars and, later, the seeds of supermassive black holes. We have modeled the explosion of a supermassive Pop III star in the dense core of a line-cooled protogalaxy with the ZEUS-MP code. We find that the supernova (SN) expands to a radius of {approx}1 kpc, briefly engulfing the entire galaxy, but then collapses back into the potential well of the dark matter. Fallback fully mixes the interior of the protogalaxy with metals, igniting a violent starburst and fueling the rapid growth of a massive black hole at its center. The starburst would populate the protogalaxy with stars in greater numbers and at higher metallicities than in more slowly evolving, nearby halos. The SN remnant becomes a strong synchrotron source that can be observed with eVLA and eMERLIN and has a unique signature that easily distinguishes it from less energetic SN remnants. Such explosions, and their attendant starbursts, may well have marked the birthplaces of supermassive black holes on the sky.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

A Novel Algorithm for Solving the Multidimensional Neutron Transport Equation on Massively Parallel Architectures

We employ the Integral Transport Matrix Method (ITMM) as the kernel of new parallel solution methods for the discrete ordinates approximation of the within-group neutron transport equation. The ITMM abandons the repetitive mesh sweeps of the traditional source iterations (SI) scheme in favor of constructing stored operators that account for the direct coupling factors among all the cells' fluxes and between the cells' and boundary surfaces' fluxes. The main goals of this work are to develop the algorithms that construct these operators and employ them in the solution process, determine the most suitable way to parallelize the entire procedure, and evaluate the behavior and parallel performance of the developed methods with increasing number of processes, P. The fastest observed parallel solution method, Parallel Gauss-Seidel (PGS), was used in a weak scaling comparison with the PARTISN transport code, which uses the source iteration (SI) scheme parallelized with the Koch-baker-Alcouffe (KBA) method. Compared to the state-of-the-art SI-KBA with diffusion synthetic acceleration (DSA), this new method- even without acceleration/preconditioning-is completitive for optically thick problems as P is increased to the tens of thousands range. For the most optically thick cells tested, PGS reduced execution time by an approximate factor of three for problems with more than 130 million computational cells on P = 32,768. Moreover, the SI-DSA execution times's trend rises generally more steeply with increasing P than the PGS trend. Furthermore, the PGS method outperforms SI for the periodic heterogeneous layers (PHL) configuration problems. The PGS method outperforms SI and SI-DSA on as few as P = 16 for PHL problems and reduces execution time by a factor of ten or more for all problems considered with more than 2 million computational cells on P = 4.096.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2014

Physics Integration KErnels (PIKE)

Pike is a software library for coupling and solving multiphysics applications. It provides basic interfaces and utilities for performing code-to-code coupling. It provides simple "black-box" Picard iteration methods for solving the coupled system of equations including Jacobi and Gauss-Seidel solvers. Pike was developed originally to couple neutronics and thermal fluids codes to simulate a light water nuclear reactor for the Consortium for Simulation of Light-water Reactors (CASL) DOE Energy Innovation Hub. The Pike library contains no physics and just provides interfaces and utilities for coupling codes. It will be released open source under a BSD license as part of the Trilinos solver framework (trilinos.org) which is also BSD. This code provides capabilities similar to other open source multiphysics coupling libraries such as LIME, AMP, and MOOSE.

Pawlowski, Roger↗

DOE OSTI2014

Physics Integration KErnels (PIKE)

Pawlowski, Roger↗

DOE OSTI2021

Introducing PRIMRE’s MRE Software Knowledge Hub

This paper focuses on the role of the Marine Renewable Energy (MRE) Software Knowledge Hub on the Portal and Repository for Information on Marine Renewable Energy (PRIMRE). The MRE Software Knowledge Hub provides online services for MRE software users and developers, and seeks to develop assessments and recommendations for improving MRE software in the future. Online software discovery platforms, known as the Code Hub and the Code Catalog, are provided. The Code Hub is a collection of open-source MRE software that includes a landing page with search functionality, linked to files hosted on the MRE Code Hub GitHub organization. The Code Catalog is a searchable online platform for discovery of useful (open-source or commercial) software packages, tools, codes, and other software products. To gather information about the existing MRE software landscape, a software survey is being performed, the preliminary results of which are presented herein. Initially, the data collected in the MRE software survey will be used to populate the MRE Software knowledge hub on PRIMRE, and future work will use data from the survey to perform a gap analysis and develop a vision for future software development. Additionally, as one of PRIMRE’s roles is to support development of MRE software within project partners, a silo of knowledge relating to best practices has been gathered.

marine, energy, open source software, MHKiT, MRE C↗

DOE OSTI2015

A novel approach to correct the coded aperture misalignment for fast neutron imaging

Aperture alignment is crucial for the diagnosis of neutron imaging because it has significant impact on the coding imaging and the understanding of the neutron source. In our previous studies on the neutron imaging system with coded aperture for large field of view, “residual watermark,” certain extra information that overlies reconstructed image and has nothing to do with the source is discovered if the peak normalization is employed in genetic algorithms (GA) to reconstruct the source image. Some studies on basic properties of residual watermark indicate that the residual watermark can characterize coded aperture and can thus be used to determine the location of coded aperture relative to the system axis. In this paper, we have further analyzed the essential conditions for the existence of residual watermark and the requirements of the reconstruction algorithm for the emergence of residual watermark. A gamma coded imaging experiment has been performed to verify the existence of residual watermark. Based on the residual watermark, a correction method for the aperture misalignment has been studied. A multiple linear regression model of the position of coded aperture axis, the position of residual watermark center, and the gray barycenter of neutron source with twenty training samples has been set up. Using the regression model and verification samples, we have found the position of the coded aperture axis relative to the system axis with an accuracy of approximately 20 μm. Conclusively, a novel approach has been established to correct the coded aperture misalignment for fast neutron coded imaging.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2014

Benchmarking kinetic calculations of resistive wall mode stability

Validating the calculations of kinetic resistive wall mode (RWM) stability is important for confidently predicting RWM stable operating regions in ITER and other high performance tokamaks for disruption avoidance. Benchmarking the calculations of the Magnetohydrodynamic Resistive Spectrum—Kinetic (MARS-K) [Y. Liu et al., Phys. Plasmas 15, 112503 (2008)], Modification to Ideal Stability by Kinetic effects (MISK) [B. Hu et al., Phys. Plasmas 12, 057301 (2005)], and Perturbed Equilibrium Nonambipolar Transport (PENT) [N. Logan et al., Phys. Plasmas 20, 122507 (2013)] codes for two Solov'ev analytical equilibria and a projected ITER equilibrium has demonstrated good agreement between the codes. The important particle frequencies, the frequency resonance energy integral in which they are used, the marginally stable eigenfunctions, perturbed Lagrangians, and fluid growth rates are all generally consistent between the codes. The most important kinetic effect at low rotation is the resonance between the mode rotation and the trapped thermal particle's precession drift, and MARS-K, MISK, and PENT show good agreement in this term. The different ways the rational surface contribution was treated historically in the codes is identified as a source of disagreement in the bounce and transit resonance terms at higher plasma rotation. Calculations from all of the codes support the present understanding that RWM stability can be increased by kinetic effects at low rotation through precession drift resonance and at high rotation by bounce and transit resonances, while intermediate rotation can remain susceptible to instability. The applicability of benchmarked kinetic stability calculations to experimental results is demonstrated by the prediction of MISK calculations of near marginal growth rates for experimental marginal stability points from the National Spherical Torus Experiment (NSTX) [M. Ono et al., Nucl. Fusion 40, 557 (2000)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

The MolSSI QCArchive project: An open-source platform to compute, organize, and share quantum chemistry data

The Molecular Sciences Software Institute's (MolSSI) Quantum Chemistry Archive (QCArchive) project is an umbrella name that covers both a central server hosted by MolSSI for community data and the Python-based software infrastructure that powers automated computation and storage of quantum chemistry (QC) results. The MolSSI-hosted central server provides the computational molecular sciences community a location to freely access tens of millions of QC computations for machine learning, methodology assessment, force-field fitting, and more through a Python interface. Facile, user-friendly mining of the centrally archived quantum chemical data also can be achieved through web applications found at the website. The software infrastructure can be used as a standalone platform to compute, structure, and distribute hundreds of millions of QC computations for individuals or groups of researchers at any scale. The QCArchiveInfrastructure is open-source (BSD-3C), code repositories can be found at github, and releases can be downloaded via PyPI and Conda. This article is categorized under: Electronic Structure Theory > Ab Initio Electronic Structure Methods Software > Quantum Chemistry Data Science > Computer Algorithms and Programming

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Short-term electricity load forecasting: Application-driven evaluation of machine learning models across spatial and temporal scales

As we transition towards a decarbonized economy, the integration of variable renewable energy resources and new demands (e.g., electric vehicles, heat pumps) into the electricity grid places unprecedented pressure on grid operators to effectively anticipate and manage peak load. In this context, machine learning algorithms are proving to be indispensable for accurate short-term load forecasting, a crucial task to address these challenges. This study benchmarks 6 machine learning algorithms, including three neural networks and three tree-based algorithms, across various levels of spatial aggregation and time horizons (1, 4, 8, 24, and 48 h). The central contribution of this work is the comparison and analysis of load forecasting models not only based on statistical metrics, but also based on a novel error metric, which evaluates the cost implications of forecast errors for power system stakeholders. Results show that tree-based models outperform neural networks, based on statistical metrics, and yield less skewed error distributions for most spatial scales. However, through the lens of the novel error metric, neural networks are the more competitive choice, especially for forecast horizons that exceed 8 h. The study concludes with actionable recommendations to grid operators and highlights the need for the development of error metrics that link forecasting accuracy to operational costs. To promote transparency and open science, the datasets and Python code are open-sourced via a supplementary repository.

Houben, Nikolaus↗

DOE OSTI2023

Computer-controlled finishing via dynamically constraint position-velocity-time scheduler

In a Computer Numerical Controlled (CNC) finishing process, the target material removal from an optical surface is guided by the convolution between the influence function of a machine tool and its dwell time at certain points over the surface. To reduce dynamics stressing and increase machining efficiency, the dwell time must be converted to varying velocities, which are the actual inputs to the machine tool controller. Conventionally, the conversion assumed constant acceleration and relied on linear motion interpolation, which caused discontinuities in velocities. This unsmooth motion affects the material removal distribution, and, thus, the accuracy of the finished surface shape. Many modern CNC machines support the smoother, cubic-polynomial interpolated Position-Velocity-Time (PVT) motion mode; however, the conventional scheduler may fail to provide suitable velocities for the PVT. Here in this study answers this challenge by proposing a novel PVT-based velocity scheduler that achieves smooth motion while considering CNC dynamic limits. Firstly, the principle of the PVT is explained, and the PVT-based velocity scheduler is formulated. Secondly, a quadratic programming is used to optimize the velocities by imposing the CNC dynamic constraints and the C 1 continuities (zeroth and first derivatives are continuous) simultaneously. Thirdly, the smoothness and accuracy of the scheduled velocities are studied on different kinds of tool paths via simulation. Finally, a sub-0.3 nm level surface finishing experiment using ion beam figuring is demonstrated to verify the feasibility of the proposed method. The PVT-based scheduler and simulator code is open-sourced.

36 MATERIALS SCIENCE↗

DOE OSTI2018

Active learning for accelerated design of layered materials

Hetero-structures made from vertically stacked monolayers of transition metal dichalcogenides hold great potential for optoelectronic and thermoelectric devices. Discovery of the optimal layered material for specific applications necessitates the estimation of key material properties, such as electronic band structure and thermal transport coefficients. However, screening of material properties via brute force ab initio calculations of the entire material structure space exceeds the limits of current computing resources. Moreover, the functional dependence of material properties on the structures is often complicated, making simplistic statistical procedures for prediction difficult to employ without large amounts of data collection. Here, we present a Gaussian process regression model, which predicts material properties of an input hetero-structure, as well as an active learning model based on Bayesian optimization, which can efficiently discover the optimal hetero-structure using a minimal number of ab initio calculations. The electronic band gap, conduction/valence band dispersions, and thermoelectric performance are used as representative material properties for prediction and optimization. The Materials Project platform is used for electronic structure computation, while the BoltzTraP code is used to compute thermoelectric properties. Bayesian optimization is shown to significantly reduce the computational cost of discovering the optimal structure when compared with finding an optimal structure by building a regression model to predict material properties. The models can be used for predictions with respect to any material property and our software, including data preparation code based on the Python Materials Genomics (PyMatGen) library as well as python-based machine learning code, is available open source.

36 MATERIALS SCIENCE↗

DOE OSTI2013

Edible oil structures at low and intermediate concentrations. II. Ultra-small angle X-ray scattering of in situ tristearin solids in triolein

Ultra-small angle X-ray scattering has been used for the first time to elucidate, in situ, the aggregation structure of a model edible oil system. The three-dimensional nano- to micro-structure of tristearin solid particles in triolein solvent was investigated using 5, 10, 15, and 20% solids. Three different sample preparation procedures were investigated: two slow cooling rates of 0.5°/min, case 1 (22 days of storage at room temperature) and case 2 (no storage), and one fast cooling of 30°/min, case 3 (no storage). The length scale investigated, by using the Bonse-Hart camera at beamline ID-15D at the Advanced Photon Source, Argonne National Laboratory, covered the range from 300 Å to 10 μm. The unified fit and the Guinier-Porod models in the Irena software were used to fit the data. The former was used to fit 3 structural levels. Level 1 structures showed that the primary scatterers were essentially 2-dimensional objects for the three cases. The scatterers possessed lateral dimensions between 1000 and 4300 Å. This is consistent with the sizes of crystalline nanoplatelets present which were observed using cryo-TEM. Level 2 structures were aggregates possessing radii of gyration, R{sub g2} between 1800 Å and 12000 Å and fractal dimensions of either D{sub 2}=1 for case 3 or 1.8≤D{sub 2}≤2.1 for case 1 and case 2. D{sub 2} = 1 is consistent with unaggregated 1-dimensional objects. 1.8 ≤ D{sub 2} ≤ 2.1 is consistent with these 1-dimensional objects (below) forming structures characteristic of diffusion or reaction limited cluster-cluster aggregation. Level 3 structures showed that the spatial distribution of the level 2 structures was uniform, on the average, for case 1, with fractal dimension D{sub 3}≈3 while for case 2 and case 3 the fractal dimension was D{sub 3}≈2.2, which suggested that the large-scale distribution had not come to equilibrium. The Guinier-Porod model showed that the structures giving rise to the aggregates with a fractal dimension given by D{sub 2} in the unified fit level 2 model were cylinders described by the parameter s≈1 in the Guinier-Porod model. The size of the base of these cylinders was in agreement with the cryo-TEM observations as well as with the results of the level 1 unified fit model. By estimating the size of the nanoplatelets and understanding the structures formed via their aggregation, it will be possible to engineer novel lipids systems that embody desired functional characteristics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2015

The CARIBU EBIS control and synchronization system

The Californium Rare Isotope Breeder Upgrade (CARIBU) Electron Beam Ion Source (EBIS) charge breeder has been built and tested. The bases of the CARIBU EBIS electrical system are four voltage platforms on which both DC and pulsed high voltage outputs are controlled. The high voltage output pulses are created with either a combination of a function generator and a high voltage amplifier, or two high voltage DC power supplies and a high voltage solid state switch. Proper synchronization of the pulsed voltages, fundamental to optimizing the charge breeding performance, is achieved with triggering from a digital delay pulse generator. The control system is based on National Instruments realtime controllers and LabVIEW software implementing Functional Global Variables (FGV) to store and access instrument parameters. Fiber optic converters enable network communication and triggering across the platforms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Synthesis of polycaprolactone-titanium oxide multilayer films by nanosecond laser pulses and electrospinning technique for better implant fabrication

In this study, a combination of electrospinning and laser texturing is introduced as a novel method for increasing the biocompatibility of metal implants. Besides having a rough laser treated surface, the implant benefits from the high porosity and better wettability of an electrospun fibrous structure, which is a more favorable environment for cell proliferation. Titanium samples were patterned using a nanosecond laser beam and were placed as collectors in an electrospinning machine. They were then soaked in simulated body fluid for four weeks. Energy Dispersive X-ray and X-Ray Diffraction results indicate significantly more hydroxyapatite formation on laser treated samples with nanoscale fibers deposited on their surface. This shows that having a laser treated surface underneath the fibrous layer can improve short-term biocompatibility even before degradation of fibers. The thermal conductivity of the electrospun layer, measured using a Hot Disk Transient Plane Source instrument and computer code, was shown to be considerably lower than that of titanium and very close to bone. The presence of this layer can therefore be beneficial in making the implant more compatible to a biological medium. In case of dental implants, it was shown that this layer can act as a thermal barrier while a hot beverage is consumed and it can decrease the temperature rise by about 60%, which avoids any possible damage to newly formed cells during the healing period.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗