Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parallel methods”

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,531 records · Page 85
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Portable COTS RGB Wavefront System for Real-time Turbulence Mitigation

Atmospheric effects reduce the ability to detect or identify objects of interest due to spatial and spectral scintillation. Algorithms alone cannot correct for turbulence over long-range horizontal paths. While adaptive optics methods work well for astronomical viewing, they are not readily adapted for complex imagery over long-range horizontal paths. At MSS Parallel in 2016, we presented our work for the red/green/blue (RGB) wavefront sensor. Since then, we have completed our system to combine an off-the-shelf high-frame-rate USB3 color camera with a compact 90 mm Maksutov-Cassegrain telescope and an OKO deformable mirror (DM) to create a portable, real-time RGB wavefront sensor and turbulence mitigation system. The telescope’s aperture is covered with an aluminum plate holding three separate narrowband color filters. These filters, specifically matched to the relative spectral response of the sensor’s Bayer mask, provide nearly complete separation of the three color images. This device enables portable imaging through atmospheric turbulence and analysis thereof, where relative motion of the three-color Bayer-pixel video provides real-time wavefront information on the second-order statistics of Kolmogorov turbulence. We measured relative RGB pixel shifts and controlled this mirror using an ODROID single-board computer. We developed a bench setup to calibrate the relationship between pixel shifts and the deformable mirror actuators. We demonstrated the feasibility of this system on a benchtop and on scenes at 700 meters showing real time correction.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Linear differential

Apparatuses and methods of operating a linear differential (100, 600) are described herein. The linear differential (100, 600) contains a slide portion (102) with parallel right-hand and left-hand threaded rods (112, 114). Threaded onto the right-hand and left-hand threaded rods (112, 114) and attached to the slide portion (102) are right-hand and left-hand gears (116, 118). Meshed between the right-hand and left-hand gears (116, 118) and also attached to the slide portion (102) is a driven gear (200). An end effector (104) is attached to the driven gear (200) and is configured to translate along a translation axis (110) and rotate around a rotation axis (120).

Crawford, Anthony↗
DOE OSTI2023

Magnetic levitation

A magnetic levitation system is described, including a first cylinder-shaped magnet; a second cylinder-shaped magnet coaxially aligned with the first cylinder-shaped magnet; and a first cavity coaxially aligned with the first cylinder-shaped magnet; wherein the surfaces of the like-poles of the first and second cylinder-shaped magnets are parallel to each other and face each other to result in a linear magnetic field between the first and the second magnets. Methods of using a magnetic levitation system for analyzing a diamagnetic or paramagnetic sample are also described.

Ge, Shencheng↗
DOE OSTI2024

Data-Driven Unit Commitment Refinement - a Scalable Approach for Complex Modern Power Grids

Integration of renewable generation, which is often intermittent and decentralized, substantially increases the stochasticity and complexity of power grid operations. Future power systems planning will require significant computational capability to evaluate balance between demand and supply under varying conditions, both temporally and spatially. The standard approach for generation unit commitment is to use mixed-integer linear programming to find the optimal generation schedule considering ramping and generator constraints. In the future grid this poses computational scalability challenges because generation and demand are not known with certainty due to stochasticity in weather and complexity of the grid. To address this challenge, we present a data-driven unit commitment approach that can efficiently include stochastic weather impacts and contingency considerations to improve unit commitment. Our approach uses graph-based data analytics techniques on solutions to the security constrained (and possibly stochastic) economic dispatch problem to identify potential improvements to a given unit commitment. Recent breakthroughs in fully-parallel stochastic economic dispatch software allow this approach to be scalably deployed. Simulations on synthetic South Carolina and Texas grids show this method can improve grid reliability with security constraints over a set of contingencies, while also meaningfully lowering total generation cost.

Holt, Timothy↗
DOE OSTI2016

A domain-specific compiler for a parallel multiresolution adaptive numerical simulation environment

This paper describes the design and implementation of a layered domain-specific compiler to support MADNESS---Multiresolution ADaptive Numerical Environment for Scientific Simulation. MADNESS is a high-level software environment for the solution of integral and differential equations in many dimensions, using adaptive and fast harmonic analysis methods with guaranteed precision. MADNESS uses k-d trees to represent spatial functions and implements operators like addition, multiplication, differentiation, and integration on the numerical representation of functions. The MADNESS runtime system provides global namespace support and a task-based execution model including futures. MADNESS is currently deployed on massively parallel supercomputers and has enabled many science advances. Due to the highly irregular and statically unpredictable structure of the k-d trees representing the spatial functions encountered in MADNESS applications, only purely runtime approaches to optimization have previously been implemented in the MADNESS framework. This paper describes a layered domain-specific compiler developed to address some performance bottlenecks in MADNESS. The newly developed static compile-time optimizations, in conjunction with the MADNESS runtime support, enable significant performance improvement for the MADNESS framework.

Rajbhandari, Samyam↗
DOE OSTI2023

The Need for Nb3Sn Coated Cu Cavities for Future Accelerators

Based on current efforts in the U.S. on the novel concept of parallel-feed RF accelerator structures, and in the U.S. and abroad in producing Nb3Sn films on either Cu or bronze, we rec-ommend that the Particle Physics community foster R&D in Superconducting Nb3Sn coated Cu RF Cavities instead of costly bulk Nb. The paper includes methods to process the coated cavi-ties at temperatures consistent with Cu retaining its shape. A devoted global effort in develop-ing Cu cavity structures coated with Nb3Sn would make the ILC or Higgs factories more afforda-ble and more likely to be built. Not only do parallel-feed RF structures enable both higher ac-celerating gradients and higher efficiencies, but they would be applicable to both Cu and Nb3Sn coated Cu cells. Increased effort on these two techniques would synergize expenditures to-wards progress, which will converge on the choice of technology for the RF of an ILC or any fu-ture accelerator. The current methods of Nb3Sn coatings on Cu or bronze can be geared also towards standard cavity cells. In conclusion, the use of distributed coupling structure topology within improved performance parameters together with Nb3Sn coating technology can lead to a paradigm shift for superconducting linacs, with higher gradient, higher temperature of opera-tion, and reduced overall costs for any future collider.

43 PARTICLE ACCELERATORS↗
DOE OSTI2018

A conservative scheme for electromagnetic simulation of magnetized plasmas with kinetic electrons

A conservative scheme has been formulated and confirmed for gyrokinetic particle simulations of electromagnetic waves and instabilities in magnetized plasmas. An electron continuity equation derived from the drift kinetic equation is used to time advance the electron density perturbation by using the perturbed mechanical flow calculated from the parallel vector potential, and the parallel vector potential is solved by using the perturbed canonical flow from the perturbed distribution function. In gyrokinetic particle simulations using this new scheme, the shear Alfvén wave dispersion relation in the shearless slab and continuum damping in the sheared cylinder have been recovered. The new scheme overcomes the stringent requirement in the conventional perturbative simulation method that perpendicular grid size needs to be as small as electron collisionless skin depth even for the long wavelength Alfvén waves. The new scheme also circumvents the problem in the typical approach that an unphysically large parallel electric field arises due to the inconsistency between electrostatic potential calculated from the perturbed density and vector potential calculated from the perturbed canonical flow. Lastly, the gyrokinetic particle simulations of the Alfvén waves in sheared cylinder have superior numerical properties compared with the fluid simulations, which suffer from numerical difficulties associated with singular mode structures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

System Modeling and Robust Stability Region Analysis for Multi-inverters Based on VSG

As the number of inverter-interfaced distributed generators (IIDGs) increases, the decrease in the inertia of the microgrid affects the stability of microgrid voltage and current. The virtual synchronous generator (VSG) control method simulates the output characteristics of the synchronous machine to improve the inertia and damping of the microgrid. However, low-frequency oscillations can occur in systems with multi-VSGs connected in parallel, primarily when operating in island mode. Therefore, this article develops the small-signal state-space model for multi-inverters based on VSG and analyses the robust stability region (RSR) of the control parameters. Firstly, an extensible small-signal state-space model is established and simplified in the dq axis, which is a necessary preprocessing for the subsequent stability analysis. Secondly, the guardian map method based on bialternate product operation is constructed, and the RSR is established through translation mapping and rotation mapping. Therefore, the solution of RSR is transformed into the problem of identifying the Hurwitz matrix. Finally, the simulations and experiments verify that the calculated region can ensure stable operation during startup and meet the requirement to suppress low-frequency oscillations. In conclusion, the method of solving RSR proposed in this paper provides a basis for selecting controller parameters in practical engineering.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Modeling the spectral modification of lower hybrid wave in the presence of drift-wave type density fluctuation in the scrape-off-layer of the EAST tokamak

In this work, the spectrum change of the lower hybrid (LH) waves caused by a low-frequency density fluctuation in the scrape-off-layer is studied by applying the wave scattering model developed by Bonoli and Ott [Phys. Fluids 25, 359 (1982)] via a Monte Carlo method. A ray-tracing model specific to this purpose is developed to evaluate the probability distribution of both the poloidal refractive index (N θ ) and the parallel refractive index (N ∥ ) of the LH wave at the last closed flux surface (LCFS).The probability distributions of N ∥ and N θ at the LCFS are studied using the EAST parameters as a function of wave frequency, the initial N ∥ , and the polar injection position, which may influence the lower hybrid current drive efficiency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Fast synchronization-free algorithms for parallel sparse triangular solves with multiple right-hand sides

The sparse triangular solve kernels, SpTRSV and SpTRSM, are important building blocks for a number of numerical linear algebra routines. Parallelizing SpTRSV and SpTRSM on today's many-core platforms, such as GPUs, is not an easy task since computing a component of the solution may depend on previously computed components, enforcing a degree of sequential processing. As a consequence, most existing work introduces a preprocessing stage to partition the components into a group of level-sets or colour-sets so that components within a set are independent and can be processed simultaneously during the subsequent solution stage. However, this class of methods requires a long preprocessing time as well as significant runtime synchronization over-heads between the sets. To address this, we

HPC↗
DOE OSTI2016

A computationally efficient method for full-core conjugate heat transfer modeling of sodium fast reactors

For efficient and accurate temperature predictions of sodium fast reactor structures, a 3-D full-core conjugate heat transfer modeling capability is developed for an advanced system analysis tool, SAM. The hexagon lattice core is modeled with 1-D parallel channels representing the subassembly flow, and 2-D duct walls and inter-assembly gaps. The six sides of the hexagon duct wall and near-wall coolant region are modeled separately to account for different temperatures and heat transfer between coolant flow and each side of the duct wall. The Jacobian Free Newton Krylov (JFNK) solution method is applied to solve the fluid and solid field simultaneously in a fully coupled fashion. The 3-D full-core conjugate heat transfer modeling capability in SAM has been demonstrated by a verification test problem with 7 fuel assemblies in a hexagon lattice layout. In addition, the SAM simulation results are compared with RANS-based CFD simulations. Very good agreements have been achieved between the results of the two approaches.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Closed-Loop Error Learning Control for Uncertain Nonlinear Systems With Experimental Validation on a Mobile Robot

This article develops a closed-loop error learning control (CLELC) algorithm for feedback linearizable systems with experimental validation on a mobile robot. Traditional feedback and feedforward controllers are designed based on the nominal model by using the feedback linearization control (FLC) method. Then, an intelligent controller is designed based on the sliding mode learning algorithm that utilizes closed-loop error dynamics to learn the system behavior. The controllers are working in parallel, and the intelligent controller can gradually replace the feedback controller from the control of the system. In addition to the stability of the sliding mode learning algorithm, the closed-loop stability of an nth-order feedback linearizable system is proven. The simulation results demonstrate that the CLELC algorithm can improve control performance (e.g., smaller rise time, settling time, and overshoot) in the absence of uncertainties, and also provides robust control performance in the presence of uncertainties as compared to traditional FLC method. To test the efficiency and efficacy of the CLELC algorithm, the trajectory tracking problem of a tracked mobile robot is studied in real time. The experimental results demonstrate that the CLELC algorithm ensures highly accurate trajectory tracking performance than traditional FLC method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Pass-efficient methods for compression of high-dimensional turbulent flow data

The future of high-performance computing, specifically on future Exascale computers, will presumably see memory capacity and bandwidth fail to keep pace with data generated, for instance, from massively parallel partial differential equation (PDE) systems. Current strategies proposed to address this bottleneck entail the omission of large fractions of data, as well as the incorporation of in situ compression algorithms to avoid overuse of memory. To ensure that post-processing operations are successful, this must be done in a way that a sufficiently accurate representation of the solution is stored. Moreover, in situations where the input/output system becomes a bottleneck in analysis, visualization, etc., or the execution of the PDE solver is expensive, the number of passes made over the data must be minimized. In the interest of addressing this problem, this work focuses on the utility of pass-efficient, parallelizable, low-rank, matrix decomposition methods in compressing high-dimensional simulation data from turbulent flows. Additionally, a particular emphasis is placed on using coarse representation of the data – compatible with the PDE discretization grid – to accelerate the construction of the low-rank factorization. This includes the presentation of a novel single-pass matrix decomposition algorithm for computing the so-called interpolative decomposition. The methods are described extensively and numerical experiments on two turbulent channel flow data are performed. In the first (unladen) channel flow case, compression factors exceeding 400 are achieved while maintaining accuracy with respect to first- and second-order flow statistics. In the particle-laden case, compression factors of 100 are achieved and the compressed data is used to recover particle velocities. These results show that these compression methods can enable efficient computation of various quantities of interest in both the carrier and disperse phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

WE-D-BRF-04: Experimental Investigations On Ion Radiography with Beam Scanning Using a Range Telescope

Purpose: Ion beams exhibit a finite range and an inverted depth-dose profile, the Bragg peak. These favorable properties allow superior tumordose conformality, but introduce sensitivity to range uncertainties. Hence, imaging techniques play an increasingly important role to support the treatment planning and the in-vivo monitoring of the actual ion beam treatment. Methods: This work presents the experimental investigations carried out to address the feasibility of ion transmission imaging at the Heidelberg Ion Therapy center using an active raster scanning beam delivery system and a prototype range telescope set-up based on a stack of 61 parallel-plate ionization chambers (PPIC) interleaved with 3 mm absorber plates of PMMA. Results: An extensive characterization of the set-up in terms of beam parameters and settings of the read-out electronics was performed and results will be presented. A data processing method to increase the range resolution (MIRR) of the PPIC stack was developed. In this approach, the position of the maximum of the Bragg curve is deduced from the ratio of measured signals in adjacent PPIC channels. MIRR evaluation is based on Bragg curves obtained from Monte Carlo simulations and validated with experimental data acquired with the PPIC stack using ion beams. MIRR was applied to the carbon ion radiography of an anthropomorphic Alderson head phantom yielding a resolution of 0.8 mm water equivalent thickness (WET) compared to the nominal value of 3.495 mm WET given by the thickness of the absorber slabs in the PPIC stack. An absolute comparison of the Alderson phantom carbon ion transmitted image with an X-ray digitally reconstructed radiography, both converted into WET, will also be shown. Conclusion: The obtained results are very promising and motivate further developments of the system towards an eventual clinical use.This work is supported by the German Research Foundation and the German Academic Exchange Service. This work is supported by the German Research Foundation (DFG) and the German Academic Exchange Service (DAAD)

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2025

Multi-GPU porting of a phase-change cascaded lattice Boltzmann method for three-dimensional pool boiling simulations

The Lattice Boltzmann method (LBM) has proven effective in simulating phase-change phenomena, such as melting, solidification, evaporation, and boiling. In this work, we develop a highly parallelized multi-GPU implementation of LBM for three-dimensional pool boiling simulations. The code is based on the OpenACC programming model, which enables the code to be deployed efficiently on multi-core CPUs, GPUs, and potentially other accelerators, without the need for architecture-specific rewrites. To support large-scale simulations, the domain is decomposed and distributed across multiple compute nodes using MPI. We demonstrate that the code exhibits excellent scaling properties, with ideal strong-scaling running with up to 256 GPUs on the MareNostrum5 cluster.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

SU-E-T-462: Impact of the Radiochromic Film Energy Response On Dose Measurements of Low Energy Electronic Brachytherapy Sources

Purpose: We investigated the effect of the EBT3 GafChromicTM film model absorbed dose energy response when used for percent depth dose (PDD) measurements in low-energy photon beams. Methods: We measured PDDs in water from a Xoft 50 kVp source using EBT3 film, and compared them to PDD measurements acquired with a PTW-TN34013 parallel-plate ionization chamber. For the x-ray source, we simulated spectra using the EGSnrc (BEAMnrc) Monte Carlo code, and calculated Half Value Layer (HVL) at different distances from the source in water. Pieces of EBT3 film were irradiated in air and calibration curves were created in terms of air-kerma in air ((Kair)air) for different beam qualities. Pieces of EBT3 film were positioned at distances of 2–6 cm from the Xoft source in a water phantom using a custom-made holder, and irradiated at the same time. As scatter is incorporated in the measured film signal in water, measured (Kair)wat was subsequently converted into absorbed dose to water by the ratio of mass energy absorption coefficients following the AAPM TG-61 dosimetry protocol. Results: Our results show that film calibration curves obtained at beam qualities near the effective energy of the Xoft 50 kVp source in water lead to variation in absorbed dose energy dependence of the response of around 3%. However, if the calibration curve was established at MV beam quality, the error in absorbed dose could be as large as 15%. We observed agreement within 1% between PDD measurements using EBT3 film model (using a calibration curve obtained at 80 kVp, HVL=2.18 mm Al, Eeff=29.5 keV) and the parallel-plate ionization chamber. Conclusion: Accurate dose measurements using radiochromic films at low photon energies require that the radiochromic film dosimetry system be calibrated at corresponding low energies, as large absorbed dose errors are expected for calibrations performed at MV beam qualities.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Massively Parallel Bayesian Model Calibration and Uncertainty Quantification with Applications to Nuclear Fuels and Materials

The U.S. Department of Energy (DOE)’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program aims to develop predictive capabilities by applying computational methods to the analysis and design of advanced reactor and fuel cycle systems. This program has been providing engineering-scale support for the development of BISON, a high-fidelity and high-resolution fuel performance tool. Fuel behavior in a nuclear reactor is governed by a complex network of mechanisms interacting with various other physics aspects in the reactor system. Any model developed to represent the fuel behavior will likely be idealized resulting in uncertainties in their predictions compared to the observed data. As such, this report was motivated by the need to identify the sources of uncertainties and quantify and propagate them through the fuel model outputs. Such quantification of uncertainties will establish a level of model trustworthiness, identify approaches to improve the model trustworthiness, and even guide optimal experiment design for maximal information gain. To accomplish the uncertainty quantification for computational models, this report has relied on the Bayesian framework which provides probabilistic treatment of models their inputs and outputs. The current state-of-the-art on performing Bayesian Uncertainty Quantification (UQ) for nuclear engineering models using High Performance Computing (HPC) resources have been reviewed. Implementation of capabilities for massively parallel Bayesian UQ in Multiphysics Object-Oriented Simulation Environment (MOOSE) is discussed. Several verification cases are discussed to verify the accuracy of the quantified uncertainties using the developed computational capabilities in MOOSE. Then, the problem of quantifying the uncertainties in TRI-Structural isOtropic (TRISO) fuel silver release is addressed. For the first time, the uncertainties arising from the TRISO Fission Gas Release (FGR) model due to model inadequacy and experimental noise are quantified. Also, the Bayesian capabilities are applied to the calibration of the MATPRO creep model, a widely used model in several fuel assessment cases. The impact of the prediction uncertainties in the MATPRO model on the fuel cladding behavior as part of the TRIBULATION assessment case (which is an integral effects case) is investigated. This report concludes with a discussion on the future work for the UQ for computational models.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2018

Optimization of groundwater artificial recharge systems using a genetic algorithm: a case study in Beijing, China

An optimization approach is used for the operation of groundwater artificial recharge systems in an alluvial fan in Beijing, China. The optimization model incorporates a transient groundwater flow model, which allows for simulation of the groundwater response to artificial recharge. The facilities’ operation with regard to recharge rates is formulated as a nonlinear programming problem to maximize the volume of surface water recharged into the aquifers under specific constraints. This optimization problem is solved by the parallel genetic algorithm (PGA) based on OpenMP, which could substantially reduce the computation time. To solve the PGA with constraints, the multiplicative penalty method is applied. In addition, the facilities’ locations are implicitly determined on the basis of the results of the recharge-rate optimizations. Two scenarios are optimized and the optimal results indicate that the amount of water recharged into the aquifers will increase without exceeding the upper limits of the groundwater levels. Optimal operation of this artificial recharge system can also contribute to the more effective recovery of the groundwater storage capacity.

58 GEOSCIENCES↗