Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed computation”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

SU-D-16A-04: Accuracy of Treatment Plan TCP and NTCP Values as Determined Via Treatment Course Delivery Simulations

Purpose: To to determine if tumor control probability (TCP) and normal tissue control probability (NTCP) values computed on the treatment planning image are representative of TCP/NTCP distributions resulting from probable positioning variations encountered during external-beam radiotherapy. Methods: We compare TCP/NTCP as typically computed on the planning PTV/OARs with distributions of those parameters computed for CTV/OARs via treatment delivery simulations which include the effect of patient organ deformations for a group of 19 prostate IMRT pseudocases. Planning objectives specified 78 Gy to PTV1=prostate CTV+5 mm margin, 66 Gy to PTV2=seminal vesicles+8 mm margin, and multiple bladder/rectum OAR objectives to achieve typical clinical OAR sparing. TCP were computed using the Poisson Model while NTCPs used the Lyman-Kutcher-Bruman model. For each patient, 1000 30-fraction virtual treatment courses were simulated with each fractional pseudo- time-oftreatment anatomy sampled from a principle component analysis patient deformation model. Dose for each virtual treatment-course was determined via deformable summation of dose from the individual fractions. CTVTCP/ OAR-NTCP values were computed for each treatment course, statistically analyzed, and compared with the planning PTV-TCP/OARNTCP values. Results: Mean TCP from the simulations differed by <1% from planned TCP for 18/19 patients; 1/19 differed by 1.7%. Mean bladder NTCP differed from the planned NTCP by >5% for 12/19 patients and >10% for 4/19 patients. Similarly, mean rectum NTCP differed by >5% for 12/19 patients, >10% for 4/19 patients. Both mean bladder and mean rectum NTCP differed by >5% for 10/19 patients and by >10% for 2/19 patients. For several patients, planned NTCP was less than the minimum or more than the maximum from the treatment course simulations. Conclusion: Treatment course simulations yield TCP values that are similar to planned values, while OAR NTCPs differ significantly, indicating the need for probabilistic methods or PRVs for OAR risk assessment. Presenting author receives support from Philips Medical Systems.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Variational Dirichlet Blur Kernel Estimation

Blind image deconvolution involves two key objectives: 1) latent image and 2) blur estimation. For latent image estimation, we propose a fast deconvolution algorithm, which uses an image prior of nondimensional Gaussianity measure to enforce sparsity and an undetermined boundary condition methodology to reduce boundary artifacts. For blur estimation, a linear inverse problem with normalization and nonnegative constraints must be solved. However, the normalization constraint is ignored in many blind image deblurring methods, mainly because it makes the problem less tractable. In this paper, we show that the normalization constraint can be very naturally incorporated into the estimation process by using a Dirichlet distribution to approximate the posterior distribution of the blur. Making use of variational Dirichlet approximation, we provide a blur posterior approximation that considers the uncertainty of the estimate and removes noise in the estimated kernel. Furthermore, experiments with synthetic and real data demonstrate that the proposed method is very competitive to the state-of-the-art blind image restoration methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Calculation of the Phenix end-of-life test 'Control Rod Withdrawal' with the ERANOS code

The Inst. of Radiological Protection and Nuclear Safety (IRSN) acts as technical support to French public authorities. As such, IRSN is in charge of safety assessment of operating and under construction reactors, as well as future projects. In this framework, one current objective of IRSN is to evaluate the ability and accuracy of numerical tools to foresee consequences of accidents. Neutronic studies step in the safety assessment from different points of view among which the core design and its protection system. They are necessary to evaluate the core behavior in case of accident in order to assess the integrity of the first barrier and the absence of a prompt criticality risk. To reach this objective one main physical quantity has to be evaluated accurately: the neutronic power distribution in core during whole reactor lifetime. Phenix end of life tests, carried out in 2009, aim at increasing the experience feedback on sodium cooled fast reactors. These experiments have been done in the framework of the development of the 4. generation of nuclear reactors. Ten tests have been carried out: 6 on neutronic and fuel aspects, 2 on thermal hydraulics and 2 for the emergency shutdown. Two of them have been chosen for an international exercise on thermal hydraulics and neutronics in the frame of an IAEA Coordinated Research Project. Concerning neutronics, the Control Rod Withdrawal test is relevant for safety because it allows evaluating the capability of calculation tools to compute the radial power distribution on fast reactors core configurations in which the flux field is very deformed. IRSN participated to this benchmark with the ERANOS code developed by CEA for fast reactors studies. This paper presents the results obtained in the framework of the benchmark activity. A relatively good agreement was found with available measures considering the approximations done in the modeling. The work underlines the importance of burn-up calculations in order to have a fine core concentrations mesh for the calculation of the power distribution. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

An iterative reconstruction algorithm for pole figure inversion using total variation regularization

This paper describes a new discrete method for inverting X-ray pole figures to estimate the orientation distribution function (ODF). The method employs the equal-volume `cubochoric' representation for uniform discretization of orientation space, SO(3). The forward-projection model is combined with an anisotropic total variation term to iteratively determine the ODF. Here, the efficacy of the new method is evaluated with both model and experimental data and compared with existing discrete and series expansion methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Distributed Energy Resources Cybersecurity Framework: Applying the NIST Risk Management Process

In an effort to strengthen the cybersecurity posture for federal agencies and reduce the time and complexities of following the Risk Management Framework (RMF) six-step process, the National Renewable Energy Laboratory has dedicated research into expanding the existing Distributed Energy Resources Cybersecurity Framework to provide functionality that aids in the RMF steps. Users will have the opportunity to learn, document, and review the framework, saving time and resources. Existing functionality within the web application will continue to be available with added features for usability.

14 SOLAR ENERGY↗
DOE OSTI2016

Testing of Novel Application of GOTHIC{sup TM} to Modelling of Hydrogen Distribution in Intermediate Level Waste Storage Facilities - 16499

One of the strategies for treatment of intermediate level waste in the United Kingdom is to retrieve the waste from the legacy storage area, condition in a suitable storage container and then hold on-site in an intermediate level waste store until permanent storage in a geological disposal facility becomes available. During storage, hydrogen may be generated inside the containers at a slow rate by a number of mechanisms. The containers can be equipped with vents that release hydrogen and other gasses into the intermediate level waste store, a building which would be constructed to high standards to minimise ingress of moisture. With the low number of air changes associated with the building design there is the potential for gas to accumulate over time. There is, therefore, a requirement to be able to predict future gas concentrations inside the building over a range of operating conditions to ensure that a safe operating envelope is maintained. The usual approach to this problem is to model the evolution of the hydrogen distribution by use of computational fluid dynamics tools. In a situation where weeks of real time need to be simulated in order to understand the slow build-up of a very small volumetric concentration of hydrogen, the use of computational fluid dynamics is not practicable within the time scales demanded by industry. This paper investigates a new approach to tackling this problem using the thermal hydraulic code GOTHIC{sup TM} (GOTHIC{sup TM} incorporates technology developed for the electric power industry under the sponsorship of EPRI, the Electric Power Research Institute). (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Effective oscillator strength distributions of spherically symmetric atoms for calculating polarizabilities and long-range atom–atom interactions

Effective oscillator strength distributions are systematically generated and tabulated for the alkali atoms, the alkaline-earth atoms, the alkaline-earth ions, the rare gases and some miscellaneous atoms. These effective distributions are used to compute the dipole, quadrupole and octupole static polarizabilities, and are then applied to the calculation of the dynamic polarizabilities at imaginary frequencies. These polarizabilities can be used to determine the long-range C{sub 6}, C{sub 8} and C{sub 10} atom–atom interactions for the dimers formed from any of these atoms and ions, and we present tables covering all of these combinations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Soot and Spectral Radiation Modeling for a High-Pressure Turbulent Spray Flame

Simulations are performed of a transient high-pressure turbulent n-dodecane spray flame under engine-relevant conditions. An unsteady RANS formulation is used, with detailed chemistry, a semi-empirical two-equation soot model, and a particle-based transported composition probability density function (PDF) method to account for unresolved turbulent fluctuations in composition and temperature. Results from the PDF model are compared with those from a locally well-stirred reactor (WSR) model to quantify the effects of turbulence-chemistry-soot interactions. Computed liquid and vapor penetration versus time, ignition delay, and flame lift-off height are in good agreement with experiment, and relatively small differences are seen between the WSR and PDF models for these global quantities. Computed soot levels and spatial soot distributions from the WSR and PDF models show large differences, with PDF results being in better agreement with experimental measurements. An uncoupled photon Monte Carlo method with line-by-line spectral resolution is used to compute the spectral intensity distribution of the radiation leaving the flame. This provides new insight into the relative importance of molecular gas radiation versus soot radiation, and the importance of turbulent fluctuations on radiative heat transfer.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2018

Uncertainties in secondary angular distributions and their impact on shielding, criticality and kinetics analysis - 25556

The availability of the covariance data improved over the last decades. Use of covariance matrices of different origin (such as JEFF, JENDL, ENDF, TENDL, SCALE, etc.) and types (files MF31, MF33, MF34, MF35) in sensitivity and uncertainty analysis of vast list of cases ranging from fission to fusion and from criticality, kinetics and shielding to adjustment applications was already studied in the past. The present work focuses on the status and the availability of the covariance of the angular distributions of scattered neutrons (secondary angular distributions - SAD covariances). SAD uncertainties are usually neglected in classical uncertainty propagation codes, even though the methods and codes for calculating the sensitivities are available for decades and excellent agreement was demonstrated between different deterministic and Monte Carlo codes. Of particular concern is still the availability, quality and completeness of covariance matrices. Indeed, in most evaluations the SAD covariances are limited to P1 term, if available at all. Furthermore, some differences were observed between different evaluations, e.g. in the case of iron. (author)

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2024

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

DFTTK: Density Functional Theory ToolKit for high-throughput lattice dynamics calculations

In this work, we present a software package in Python for high-throughput first-principles calculations of thermodynamic properties at finite temperatures, which we refer to as DFTTK (Density Functional Theory ToolKit). DFTTK is based on the atomate package and integrates our experiences in the last decades on the development of theoretical methods and computational softwares. It includes task submissions on all major operating systems and task executions on high-performance computing environments. Furthermore, the distribution of the DFTTK package comes with examples of calculations of phonon density of states, heat capacity, entropy, enthalpy, and free energy under the quasi-harmonic phonon scheme for the stoichiometric phases of Al, Ni, Al 3 Ni, AlNi, AlNi 3 , Al 3 Ni 4 , and Al 3 Ni 5 , and the fcc solution phases treated using the special quasirandom structures at the compositions of Al 3 Ni, AlNi, and AlNi 3 .

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Communication: Dynamical and structural analyses of solid hydrogen under vapor pressure

Nuclear quantum effects play a dominant role in determining the phase diagram of H{sub 2}. With a recently developed quantum molecular dynamics simulation method, we examine dynamical and structural characters of solid H{sub 2} under vapor pressure, demonstrating the difference from liquid and high-pressure solid H{sub 2}. While stable hexagonal close-packed lattice structures are reproduced with reasonable lattice phonon frequencies, the most stable adjacent configuration exhibits a zigzag structure, in contrast with the T-shape liquid configuration. The periodic angular distributions of H{sub 2} molecules indicate that molecules are not a completely free rotor in the vapor-pressure solid reflecting asymmetric potentials from surrounding molecules on adjacent lattice sites. Discrete jumps of librational and H–H vibrational frequencies as well as H–H bond length caused by structural rearrangements under vapor pressure effectively discriminate the liquid and solid phases. The obtained dynamical and structural information of the vapor-pressure H{sub 2} solid will be useful in monitoring thermodynamic states of condensed hydrogens.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Numerical study of the temperature dependence of the NMR relaxation rate across the superfluid–Bose glass transition in one dimension

Here, we study the nuclear magnetic resonance (NMR) spin-lattice relaxation rate 1/ T 1 1 in random one-dimensional spin chains as a function of the temperature and disorder strength. In the zero temperature limit, the system displays a disorder-induced quantum phase transition between a critical Tomonaga-Luttinger liquid (TLL) phase and a localized Bose glass phase. The 1/ T 1 is investigated across this transition using large-scale simulations based on matrix product state techniques. We find that this quantity can detect the transition and probe the value of the dimensionless TLL parameter K . We also compute the NMR relaxation rate distributions for each temperature and disorder strength considered. In particular, we discuss the applicability of the stretched exponential fit to the return-to-equilibrium function to extract the 1/ T 1 experimentally. The results presented here should provide valuable insights in regards of future NMR experiments in realistic disordered spin compounds.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

A Monte Carlo approach to validation of FFF VMAT treatment plans for the TrueBeam linac

Purpose: To commission and benchmark a vendor-supplied (Varian Medical Systems) Monte Carlo phase-space data for the 6 MV flattening filter free (FFF) energy mode on a TrueBeam linear accelerator for the purpose of quality assurance of clinical volumetric modulated arc therapy (VMAT) treatment plans. A method for rendering the phase-space data compatible with BEAMnrc/DOSXYZnrc simulation software package is presented. Methods: Monte Carlo (MC) simulations were performed to benchmark the TrueBeam 6 MV FFF phase space data that have been released by the Varian MC Research team. The simulations to benchmark the phase space data were done in three steps. First, the original phase space which was created on a cylindrical surface was converted into a format that was compatible with BEAMnrc. Second, BEAMnrc was used to create field size specific phase spaces located underneath the jaws. Third, doses were calculated with DOSXYZnrc in a water phantom for fields ranging from 1 Multiplication-Sign 1 to 40 Multiplication-Sign 40 cm{sup 2}. Calculated percent depth doses (PDD), transverse profiles, and output factors were compared with measurements for all the fields simulated. After completing the benchmarking study, three stereotactic body radiotherapy (SBRT) VMAT plans created with the Eclipse treatment planning system (TPS) were calculated with Monte Carlo. Ion chamber and film measurements were also performed on these plans. 3D gamma analysis was used to compare Monte Carlo calculation with TPS calculations and with film measurement. Results: For the benchmarking study, MC calculated and measured values agreed within 1% and 1.5% for PDDs and in-field transverse profiles, respectively, for field sizes >1 Multiplication-Sign 1 cm{sup 2}. Agreements in the 80%-20% penumbra widths were better than 2 mm for all the fields that were compared. With the exception of the 1 Multiplication-Sign 1 cm{sup 2} field, the agreement between measured and calculated output factors was within 1%. It is of note that excellent agreement in output factors for all field sizes including highly asymmetric fields was achieved without accounting for backscatter into the beam monitor chamber. For the SBRT VMAT plans, the agreement between Monte Carlo and ion chamber point dose measurements was within 1%. Excellent agreement between Monte Carlo, treatment planning system and Gafchromic film dose distribution was observed with over 99% of the points in the high dose volume passing the 3%, 3 mm gamma test. Conclusions: The authors have presented a method for making the Varian IAEA compliant 6 MV FFF phase space file of the TrueBeam linac compatible with BEAMnrc/DOSXYZnrc. After benchmarking the modified phase space against measurement, they have demonstrated its potential for use in MC based quality assurance of complex delivery techniques.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

On the Computational Viability of Quantum Optimization for PMU Placement

Using optimal phasor measurement unit placement as a prototypical problem, we assess the computational viability of the current generation D-Wave Systems 2000Q quantum annealer for power systems design problems. We reformulate minimum dominating set for the annealer hardware, solve the reformulation for a standard set of IEEE test systems, and benchmark solution quality and time to solution against the CPLEX optimizer and simulated annealing. For some problem instances the 2000Q outpaces CPLEX. For instances where the 2000Q underperforms with respect to CPLEX and simulated annealing, we suggest hardware improvements for the next generation of quantum annealers.

hardware↗
DOE OSTI2023

Building the I (Interoperability) of FAIR for performance reproducibility of large-scale composable workflows in RECUP

Abstract-Scientific computing communities increasingly run their experiments using complex data- and compute-intensive workflows that utilize distributed and heterogeneous architectures targeting numerical simulations and machine learning, often executed on the Department of Energy Leadership Computing Facilities (LCFs). We argue that a principled, systematic approach to implementing FAIR principles at scale, including fine-grained metadata extraction and organization, can help with the numerous challenges to performance reproducibility posed by such workflows. We extract workflow patterns, propose a set of tools to manage the entire life cycle of performance metadata, and aggregate them in an HPC-ready framework for reproducibility (RECUP). We describe the challenges in making these tools interoperable, preliminary work, and lessons learned from this experiment.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Methods and apparatus for multi-resolution replication of files in a parallel computing system using semantic information

Techniques are provided for storing files in a parallel computing system using different resolutions. A method is provided for storing at least one file generated by a distributed application in a parallel computing system. The file comprises one or more of a complete file and a sub-file. The method comprises the steps of obtaining semantic information related to the file; generating a plurality of replicas of the file with different resolutions based on the semantic information; and storing the file and the plurality of replicas of the file in one or more storage nodes of the parallel computing system. The different resolutions comprise, for example, a variable number of bits and/or a different sub-set of data elements from the file. A plurality of the sub-files can be merged to reproduce the file.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Decentralized Distributed Proximal Policy Optimization (DD-PPO) for High Performance Computing Scheduling on Multi-User Systems

Resource allocation in High Performance Computing (HPC) environments presents a complex and multifaceted challenge for job scheduling algorithms. Beyond the efficient allocation of system resources, schedulers must account for and optimize multiple performance metrics, including job wait time and system throughput. Traditional heuristic-based scheduling algorithms increasingly struggle and lack the efficiency needed to meet the demands and address the complexity and scale of modern HPC systems. Consequently, recent research efforts have focused on leveraging advancements in Artificial Intelligence (AI) and Deep Learning (DL), particularly Reinforcement Learning (RL), to develop more adaptable and intelligent scheduling strategies. Previous RL-based scheduling approaches have explored a range of algorithms, from Deep Q-Networks (DQN) to Proximal Policy Optimization (PPO), and more recently, hybrid methods that integrate Graph Neural Networks (GNNs) with RL techniques. However, a common limitation across these methods is their reliance on relatively small datasets, with few methods being evaluated using large-scale, multi-million-job trace datasets representative of real-world HPC workloads. Moreover, existing RL schedulers face scalability issues due to centralized policy updates, which hinder training efficiency and performance when applied to large datasets. This study introduces a novel RL-based scheduler utilizing Decentralized Distributed Proximal Policy Optimization (DD-PPO) algorithm, which supports large-scale distributed training across multiple workers without requiring parameter synchronization at every step. By eliminating reliance on centralized updates to a shared policy, the DD-PPO scheduler enhances scalability, training efficiency, and sample utilization. Experimental validation using a large real-world dataset containing over 11.5 million job traces collected from petascale HPC systems over six years assesses the influence of dataset scale on training effectiveness and compares DD-PPO performance to traditional and advanced scheduling approaches. The experimental results demonstrate improved scheduling performance in comparison to both heuristic-based schedulers and existing RL-based scheduling algorithms.

AI↗