Search NASA⌕ Search

SEARCH · Search NASA

Results for “Time parallelization”

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,621 records · Page 90

Highly Efficient Cooperative Catalysis by Co III (Porphyrin) Pairs in Interpenetrating Metal-Organic Frameworks

A series of porous twofold interpenetrated In-Co III (porphyrin) metal–organic frameworks (MOFs) were constructed by in situ metalation of porphyrin bridging ligands and used as efficient cooperative catalysts for the hydration of terminal alkynes. The twofold interpenetrating structure brings adjacent Co III (porphyrins) in the two networks parallel to each other with a distance of about 8.8 Å, an ideal distance for the simultaneous activation of both substrates in alkyne hydration reactions. As a result, the In-Co III (porphyrin) MOFs exhibit much higher (up to 38 times) catalytic activity than either homogeneous catalysts or MOF controls with isolated Co III (porphyrin) centers, thus highlighting the potential application of MOFs in cooperative catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Additive manufactured bipolar plate for high-efficiency hydrogen production in proton exchange membrane electrolyzer cells

Additive manufacturing (AM) technology is capable of fast and low-cost prototyping from complex 3D digital models. To take advantage of this technology, a stainless steel (SS) plate with parallel flow field served as a combination of a cathode bipolar plate and a current distributor; it was fabricated using selective laser melting (SLM) techniques and investigated in a proton exchange membrane electrolyzer cell (PEMEC) in-situ for the first time. The experimental results show that the PEMEC with an AM SS cathode bipolar plate can achieve an excellent performance for hydrogen production for a voltage of 1.779 V and a current density of 2.0 A/cm 2 . The AM SS cathode bipolar plate was also characterized by SEM and EDS, and the results show a uniform elemental distribution across the plate with very limited oxidization. As a result, this research demonstrates that AM method could be a route to aid cost-effective and rapid development of PEMECs.

08 HYDROGEN↗

Efficient algorithms and implementations of entropy-based moment closures for rarefied gases

We present efficient algorithms and implementations of the 35-moment system equipped with the maximum-entropy closure in the context of rarefied gases. While closures based on the principle of entropy maximization have been shown to yield very promising results for moderately rarefied gas flows, the computational cost of these closures is in general much higher than for closure theories with explicit closed-form expressions of the closing fluxes, such as Grad's classical closure. Following a similar approach as Garrett et al. (2015) , we investigate efficient implementations of the computationally expensive numerical quadrature method used for the moment evaluations of the maximum-entropy distribution by exploiting its inherent fine-grained parallelism with the parallelism offered by multi-core processors and graphics cards. We show that using a single graphics card as an accelerator allows speed-ups of two orders of magnitude when compared to a serial CPU implementation. To accelerate the time-to-solution for steady-state problems, we propose a new semi-implicit time discretization scheme. The resulting nonlinear system of equations is solved with a Newton type method in the Lagrange multipliers of the dual optimization problem in order to reduce the computational cost. Additionally, fully explicit time-stepping schemes of first and second order accuracy are presented. We investigate the accuracy and efficiency of the numerical schemes for several numerical test cases, including a steady-state shock-structure problem.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Capacitively coupled pickup in MCP-based photodetectors using a conductive metallic anode

We have designed and tested a robust 20×20 cm 2 thin metal film internal anode capacitively coupled to an external array of signal pads or micro-strips for use in fast microchannel plate photodetectors. The internal anode, in this case a 10 nm-thick NiCr film deposited on a 96% pure Al 2 O 3 3 mm-thick ceramic plate and connected to HV ground, provides the return path for the electron cascade charge. The multi-channel pickup array consists of a printed-circuit card or glass plate with metal signal pickups on one side and the signal ground plane on the other. The pickup can be put in close proximity to the bottom outer surface of the sealed photodetector, with no electrical connections through the photodetector hermetic vacuum package other than a single ground connection to the internal anode. Two pickup patterns were tested using a small commercial MCP-PMT as the signal source: 1) parallel 50 Ω 25-cm-long micro-strips with an analog bandwidth of 1.5 GHz, and 2) a 20×20 cm 2 array of 2-dimensional square ‘pads’ with sides of 1.27 cm or 2.54 cm. The rise-time of the fast input pulse is maintained for both pickup patterns. For the pad pattern, we observe 80% of the directly coupled amplitude. For the strip pattern we measure 34% of the directly coupled amplitude on the central strip of a broadened signal. The physical decoupling of the photodetector from the pickup pattern allows easy customization for different applications while maintaining high analog bandwidth.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Impact of detector simulation in particle physics collider experiments

Through the last three decades, precise simulation of the interactions of particles with matter and modeling of detector geometries has proven to be of critical importance to the success of the international high-energy physics experimental programs. For example, the detailed detector modeling and accurate physics of the Geant4-based simulation software of the CMS and ATLAS particle physics experiments at the European Center of Nuclear Research (CERN) Large Hadron Collider (LHC) was a determinant factor for these collaborations to deliver physics results of outstanding quality faster than any hadron collider experiment ever before. This review article highlights the impact of detector simulation on particle physics collider experiments. It presents numerous examples of the use of simulation, from detector design and optimization, through software and computing development and testing, to cases where the use of simulation samples made a difference in the accuracy of the physics results and publication turnaround, from data-taking to submission. It also presents the economic impact and cost of simulation in the CMS experiment. Future experiments will collect orders of magnitude more data, taxing heavily the performance of simulation and reconstruction software for increasingly complex detectors. Consequently, it becomes urgent to find solutions to speed up simulation software in order to cope with the increased demand in a time of flat budgets. The study ends with a short discussion on the potential solutions that are being explored, by leveraging core count growth in multicore machines, using new generation coprocessors, and re-engineering of HEP code for concurrency and parallel computing.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Cryogenic light detectors with thermal signal amplification for 0 νββ search experiments

As a step towards the realization of cryogenic-detector experiments to search for neutrinoless double-beta decay (such as CROSS, BINGO, and CUPID), we investigated a batch of 10 Ge light detectors (LDs) assisted by Neganov-Trofimov-Luke (NTL) signal amplification. Each LD was assembled with a large cubic light-emitting crystal (45 mm side) using the recently developed CROSS mechanical structure. The detector array was operated at milli-Kelvin temperatures in a pulse-tube cryostat at the Canfranc underground laboratory in Spain. We achieved good performance with scintillating bolometers from CROSS, made of Li 2 100 MoO 4 crystals and used as reference detectors of the setup, and with all LDs tested (except for a single device that encountered an electronics issue). No leakage current was observed for 8 LDs with an electrode bias up to 100 V. Operating the LDs at an 80 V electrode bias applied in parallel, we obtained a gain of around 9 in the signal-to-noise ratio of these devices, allowing us to achieve a baseline noise RMS of O(10 eV). Thanks to the strong current polarization of the temperature sensors, the time response of the devices was reduced to around half a millisecond in rise time. The achieved performance of the LDs was extrapolated via simulations of pile-up rejection capability for several configurations of the CUPID detector structure. Despite the sub-optimal noise conditions of the LDs (particularly at high frequencies), we demonstrated that the NTL technology provides a viable solution for background reduction in CUPID.

47 OTHER INSTRUMENTATION↗

Quantum Dynamics at Scale: Ultrafast Control of Emergent Functional Materials

Confluence of extreme-scale quantum dynamics simulations (i.e. quantum@scale) and cutting-edge x-ray free-electron laser experiments are revolutionizing materials science. An archetypal example is the exciting concept of using picosecond light pulses to control emergent material properties on demand in atomically-thin layered materials. This paper describes efforts to scale our quantum molecular dynamics engine toward the United States' first exaflop/s computer, under an Aurora Early Science Program project named "Metascalable layered material genome". Key algorithmic and computing techniques incorporated are: (1) globally-scalable and locally-fast solvers within a linear-scaling divide-conquer-recombine algorithmic framework; (2) algebraic `BLASification' of computational kernels; and (3) data alignment and loop restructuring, along with register and cache blocking, for enhanced vectorization and efficient memory access. The resulting weak-scaling parallel efficiency was 0.93 on 131,072 Intel Xeon Phi cores for a 56.6 million atom (or 169 million valence-electron) system, whereas the various code transformations achieved 5-fold speedup. The optimized simulation engine allowed us for the first time to establish a significant effect of substrate on the dynamics of layered material upon electronic excitation.

Algebraic BLASification↗

GridPACK™ : A Framework for Developing Power Grid Simulations on High-Performance Computing Platforms

This paper describes the GridPACK™ framework, which is designed to help power grid engineers develop modeling software capable of running on high performance computers. The framework makes extensive use of software templates to provide high level functionality while at the same time allowing developers the freedom to express whatever models and algorithms they are using. GridPACK™ contains modules for setting up distributed power grid networks, assigning buses and branches with arbitrary behaviors to the network, creating distributed matrices and vectors and using parallel linear and non-linear solvers to solve algebraic equations. It also provides mappers to create matrices and vectors based on properties of the network and functionality to support IO and to mana

Palmer, Bruce J.↗

Polarized 3He gas circulating technologies for neutron analyzers

We outline our project to develop a circulating polarized helium-3 system for developing of large, quasi-continuously operating neutron analyzers. The project consisted of four areas: 1) Development of robust external cavity narrowed diode laser output with spectral line width < 0.17 nm and power of 2000 W. 2) Development of large glass polarizing cells using cell surface treatments to obtain long relaxation lifetimes. 3) Refinements of the circulation system with an emphasis on gas purification and materials testing. 4) Design/fabrication of a new polarizer system. 5) Preliminary testing of the new polarizer. 1. Developed Robust High-Power Narrowed Laser The optical configuration of the laser was discussed in the proposal and will be reviewed in the body of this report. The external cavity is configured to mutually lock the wavelength of five 10-bar laser stacks. All the logistical milestones were been met and critical subsystems- laser stack manifold and power divider, external laser cavity, and output telescope- were assembled and tested at low power. Each individual bar is narrowed to ~0.05 nm; when combined the laser has a cumulative spectral width of 0.17 nm across the entire beam due to variations of the bars central wavelength by +/- 0.1 nm, which is similar to that of Volume Bragg Grating narrowed laser bars. This configuration eliminates the free-running “pedestal” that occurs in other external cavity diode lasers. The full-scale laser was completed in 2016 and was used in both the older and newer helium polarizers. This laser was operated at 75% power for periods of up to 8 hours. Once installed, the spectrum became slightly broader (~.25 nm) at full power; this is likely due to very slight misalignments that occurred during handling. 2. Developed the processes to create uniform sintered sol-gel coatings. Our work on cell development comprised: 1) Production of large GE180 cells and explore different means of cell preparation, and 2) Development of apply sol-gel coatings to the interior of both borosilicate and aluminosilicate cells. We applied six sol-gel coatings. By modifying the mixture and developing procedures to drain and dry the cell, we produced visually uniform coatings on the interior of the cells. We now have perfected that process as described below in our report. We were able to accelerate the testing of cells using an ex situ method that avoids installing each cell into a polarizer. In the project’s last year, we conducted 38 external tests of 8 different cells. We also installed two sol-gel coated cells in our polarizers. We created cell with long ex situ relaxation lifetimes, one of which exceeded 40 hours. However, when installed in the polarizer the measured lifetime is 8 hours or less. 3. Demonstrated cycling of polarized gas and ex situ cell testing We are now cycling polarized gas from the polarizer to glass vessels and back. This has allowed us, for the first time, to make ex situ T1 measurements of polarizing cells without installing them into the polarizer itself. This has greatly improved our productivity in producing cells and evaluating our cell preparation processes. We continued development of the gas handling system in parallel with fabricating new polarizer. The integrated system was tested by the end of 2016. We now regularly cycle gas into and out of the polarizer. 4. Completed new polarizer infrastructure and control systems. We completed the new polarizer infrastructure in November 2016. The polarizer subsystems are 1) the frame, 2) the oil flow system, 3) the gas handling system, 4) the pressure vessel, with embedded solenoid, 5) cell mounting hardware with heat spreaders, and 6) electrical power and instrumentation. 5. Carried out initial tests of polarizer. We completed initial testing of the polarizer in April and May of 2017. These tests were carried out for periods up to 6 hours with laser power between 750 and 1300 Watts. The laser performed well and the polarization with asymptotic to 45 percent, which was below expectations. This low value resulted from a stationary thermal inversion in the cell that caused most of the laser power to be absorbed near the laser inlet window and deprived the lower portions of the cell of pumping laser light. Possible solutions to this problem include enhanced cooling of the cell near the laser entry and slight detuning of the laser. 6. Ongoing work. Our polarizer development efforts are ongoing to pursue our interest in neutron analyzers, nuclear targets, and providing helium for medical imaging. Current tests in the pipeline include: 1. Testing cooling enhancements to improve laser penetration of spectrally narrow lasers; 2. Testing of a cell with isolation valves that minimizes diffusive contact with gas handling hardware during polarization; 3. Testing of smaller hybrid cells with reduced alkali loading; 4. Producing polarized helium-3 for MRI imaging at the University of Missouri.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

GPU Acceleration of Volume Fraction and Centroid Computation from General Shapes on Unstructured Mesh

One important requirement for coupling multi-physics models is a means to remap fields between meshes in an accurate and conservative manner. Since the remapping algorithm is embarrassingly parallel, this motivates the use of Graphical Processing Units (GPU) for such routine. However, since our remapping code, portage, is fairly complex, we explore the use of GPUs through a simpler problem of generating material’s volume fraction and centroid inside each cell of the mesh. We use the same programming structure that is used in portage, namely, on-node parallelism via NVIDIA’s thrust library. Using a similar programming structure enables us to draw lessons from this study that can be applied to the more complex problem. We have found that while the speed-up is pronounced (10-100 times), significant effort is required to re-factor the code so that thrust can be used with a CUDA backend. Moreover, thrust takes away the fine grain control that vanilla CUDA programming gives us.

97 MATHEMATICS AND COMPUTING↗

Challenges and Recommended Policies for Simultaneous Global Implementation of Low-GWP Refrigerants and High Efficiency in Room Air Conditioners

Increasing incomes, electrification, and urbanization—as well as a warming world—are driving up the global stock of air conditioners (ACs), particularly in emerging economies with hot climates. AC energy consumption is expected to increase substantially as the global stock of room ACs rises to 1.5 billion in 2030 and 2.5 billion in 2050. Hence, improving AC energy efficiency will be critical to reducing AC energy, cost (consumer lifecycle cost, electricity generation cost, etc.), peak load, and emissions impacts. The 2016 Kigali Amendment to the Montreal Protocol offers an opportunity to improve AC energy efficiency in tandem with the phasedown of high global warming potential (GWP) hydrofluorocarbon (HFC) refrigerants. Based on the most recent information, a literature review, and interviews with manufacturers and industry experts, we find the main barriers to deploying high-efficiency ACs include concerns about market demand and cost, which could be mitigated by appropriately improved design of market-transformation programs such as standards and labeling, incentive, and procurement programs. The main barriers to the low-GWP refrigerant transition include the need for timely revision of safety standards and associated costs for capacity-building activities allowing safe use of low-GWP refrigerants in ACs. Policy action and the market transformation can be accelerated by advancing the refrigerant transition and efficiency improvements in parallel.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Grizzly: A Multi-scale and Multi-Physics Tool to Model Aging of Nuclear Power Plant Components

Grizzly is a multi-scale multi-physics tool being developed at Idaho National Laboratory (INL) as part of the US Department of Energy's Light Water Reactor Sustainability program to provide improved safety assessments of systems, components, and structures in nuclear power plants subjected to age-related degradation. Its goal is to provide an improved scientific basis for decisions surrounding license renewal, which would permit operation of commercial nuclear power plants beyond 60 years. Grizzly is based on INL's MOOSE framework, which enables multi-physics simulations in a massively parallel computing environment. The reactor pressure vessel (RPV) has been chosen as the initial application for Grizzly. The primary concern for the RPV is that it becomes increasingly susceptible to fracture initiation at the sites of pre-existing flaws over time. This is because exposure to irradiation and high temperatures embrittles the steel. Assessing the safety of an RPV under transient conditions involves simulations at multiple scales. At the engineering scale, strongly coupled equations of heat conduction and solid mechanics are solved to simulate the global response of the RPV to accident conditions. Sub-models are used to represent regions with pre-existing flaws and the stress intensity factors are evaluated from the domain integrals. Crack propagation is then determined using a model that provides the temperature and neutron fluence dependent fracture toughness of the RPV steel. The ability of Grizzly to provide accurate engineering scale estimates of RPV safety beyond 60 years of operation strongly relies on the predictability of the embrittlement model. Multiple models for embrittlement have been proposed, and the physically based empirical model is implemented in Grizzly. However, the formation of late blooming phases greatly reduces the confidence on the extrapolated values obtained from embrittlement model, which is further exacerbated by the absence of real field data for validation. Experimentally validated lower length scale models can significantly improve the predictability of the embrittlement model beyond the calibrated range and the development of such models is another ongoing effort in Grizzly. The modeling activity in Grizzly spans across different length scales, and can be divided into: (i) Irradiation induced microstructure evolution at the atomistic and single crystal length scale using molecular dynamics, atomistic kinetic Monte-Carlo, phase-field and rate theory; (ii) Effect of the irradiation induced defects on the flow stress and fracture processes at the polycrystalline scale using crystal plasticity. Subsequently, the prediction from these models can be used to obtain irradiation induced ductile-to-brittle transition temperature (DBTT) shifts from specimen-scale simulations; (iii) Global thermo-mechanical response of the RPV under normal and accident scenarios, and probabilistic risk analyses at the engineering scale using the DBTT shifts obtained from the lower length scale models. The item (i) has been discussed in a previous report. In this extended abstract an overview of items (ii) and (iii) is provided. The extended abstract summarizes the current status of the modeling effort under Grizzly. At the engineering scale, reduced order model development is underway using the statistical tool RAVEN. This will allow probabilistic analysis considering a wide range of flaw topologies in a computationally efficient manner. Coupling the irradiation induced defect growth predictions to crystal plasticity is also in progress. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Mechanisms of material removal and mass transport in focused ion beam nanopore formation

Despite the widespread use of focused ion beam (FIB) processing as a material removal method for applications ranging from electron microscope sample preparation to nanopore processing for DNA sequencing, the basic material removal mechanisms of FIB processing are not well understood. We present the first complete atomistic simulation of high-flux FIB using large-scale parallel molecular dynamics (MD) simulations of nanopore fabrication in freestanding thin films. We focus on the root mechanisms of material removal and rearrangement and describe the role of explosive boiling in forming nanopores. FIB nanopore fabrication is typically understood to occur via sputter erosion. This can be shown to be the case in low flux systems, where individual ion impacts are sufficiently separated in time that they may be considered as independent events. But our detailed MD simulations show that in high flux FIB processing, above a threshold level at which thermal effects become significant, the primary mechanism of material removal changes to a significantly accelerated, thermally dominated process. Under these conditions, the target is heated by the ion beam faster than heat is conducted away by the material, leading quickly to melting, and then continued heating to nearly the material critical temperature. This leads to explosive boiling of the target material with spontaneous bubble formation and coalescence. Mass is rapidly rearranged at the atomistic scale, and material removal occurs orders of magnitude faster than would occur by simple sputtering. While the phenomenology is demonstrated computationally in silicon, it can be expected to occur at lower beam fluxes in other cases where thermal conduction is suppressed due to material properties, geometry, or ambient thermal conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Testing the contribution of dispersal to microbial succession following a wildfire

ABSTRACT Given increased wildfire activity, there is growing interest in understanding the drivers of microbial succession after fire. Dispersal may be especially important to post-fire succession as biotic communities can be more susceptible to invasion following a disturbance. Here, we experimentally manipulated dispersal into disturbed leaf litter communities collected following a wildfire and tracked bacterial and fungal dispersal assemblages over time. We show that the identity and source of microbes immigrating into the soil surface post-fire change across time with seasonal shifts and the reemergence of aboveground vegetation. Further, dispersal significantly contributed to the reassembly of leaf litter microbial communities after the fire, producing an increasingly distinct assembly trajectory. The effect of dispersal on α-diversity and β-diversity was ecosystem dependent but, unexpectedly, influenced bacterial and fungal communities in a similar manner within ecosystems. Collectively, these results demonstrate that dispersal explicitly alters the course of microbial community succession following a wildfire and may impact bacteria and fungi in parallel ways, despite differing in traits expected to alter dispersal patterns. IMPORTANCE Identifying the mechanisms underlying microbial community succession is necessary for predicting how microbial communities, and their functioning, will respond to future environmental change. Dispersal is one mechanism expected to affect microbial succession, yet the difficult nature of manipulating microorganisms in the environment has limited our understanding of its contribution. Using a dispersal exclusion experiment, this study isolates the specific effect of environmental dispersal on bacterial and fungal community assembly over time following a wildfire. The work demonstrates the potential to quantify dispersal impacts on soil microbial communities over time and test how dispersal might further interact with other assembly processes in response to environmental change.

59 BASIC BIOLOGICAL SCIENCES↗

Parallel Fokker–Planck-DSMC algorithm for rarefied gas flow simulation in complex domains at all Knudsen numbers

A major challenge for the conventional Direct Simulation Monte Carlo (DSMC) technique lies in the fact that its computational cost becomes prohibitive in the near continuum regime, where the Knudsen number (Kn)—characterizing the degree of rarefaction—becomes small. In contrast, the Fokker–Planck (FP) based particle Monte Carlo scheme allows for computationally efficient simulations of rarefied gas flows in the low and intermediate Kn regime. The Fokker–Planck collision operator—instead of performing binary collisions employed by the DSMC method—integrates continuous stochastic processes for the phase space evolution in time. This allows for time step and grid cell sizes larger than the respective collisional scales required by DSMC. Dynamically switching between the FP and the DSMC collision operators in each computational cell is the basis of the combined FP-DSMC method, which has been proven successful in simulating flows covering the whole Kn range. Until recently, this algorithm had only been applied to two-dimensional test cases. In this contribution, we present the first general purpose implementation of the combined FP-DSMC method. Utilizing both shared- and distributed-memory parallelization, this implementation provides the capability for simulations involving many particles and complex geometries by exploiting state of the art computer cluster technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Nonlinear electronic excitations in crystalline solids using meta-generalized gradient approximation and hybrid functional in time-dependent density functional theory

We develop methods to calculate electron dynamics in crystalline solids in real-time time-dependent density functional theory employing exchange-correlation potentials which reproduce band gap energies of dielectrics; a meta-generalized gradient approximation was proposed by Tran and Blaha [Phys. Rev. Lett. 102, 226401 (2009)] (TBm-BJ) and a hybrid functional was proposed by Heyd, Scuseria, and Ernzerhof [J. Chem. Phys. 118, 8207 (2003)] (HSE). In time evolution calculations employing the TB-mBJ potential, we have found it necessary to adopt the predictor-corrector step for a stable time evolution. We have developed a method to evaluate electronic excitation energy without referring to the energy functional which is unknown for the TB-mBJ potential. For the HSE functional, we have developed a method for the operation of the Fock-like term in Fourier space to facilitate efficient use of massive parallel computers equipped with graphic processing units. We compare electronic excitations in silicon and germanium induced by femtosecond laser pulses using the TB-mBJ, HSE, and a simple local density approximation (LDA). At low laser intensities, electronic excitations are found to be sensitive to the band gap energy: they are close to each other using TB-mBJ and HSE and are much smaller in LDA. At high laser intensities close to the damage threshold, electronic excitation energies do not differ much among the three cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Development and applications of retro-reflective surfaces for ultrasound in LBE

The Belgian Nuclear Research Centre SCK-CEN is in the process of developing MYRRHA - a fast flux research reactor to replace the aging BR2. MYRRHA is conceptualized as an accelerator driven system cooled by lead bismuth eutectic mixture (LBE). As LBE is opaque to visual light, ultrasonic measurement techniques are employed as the main technology to provide feedback to submerged operations when needed. Conceptually, MYRRHA is a pool type reactor divided in a cold lower plenum and hot upper plenum separated by a diaphragm that forces the main flow through the core. The main flow is cooled by four heat exchangers and driven by two liquid metal pumps. One of the tasks tackled using ultrasound is locating a potentially lost fuel assembly to assist a recovery operation. As all fuel manipulations in MYRRHA are performed in the lower plenum, a potentially lost fuel assembly is located in the lower plenum. Buoyancy will force the lost fuel assembly to float against the diaphragm unless it is still partially inserted in the core. Because of the latter situation, an ultrasonic scan localizing the fuel assembly should be performed from a large distance to avoid a collision with such a partially inserted fuel assembly. Unfortunately, standard machined stainless steel objects, such as a fuel assembly, reflect an ultrasonic pulse in a specular way which induces stringent requirements on the alignment of the ultrasonic sensor with respect to the fuel assembly as we cannot rely on diffuse reflections and/or scattering of the ultrasonic pulse. Moreover, increasing the distance also increases geometric spreading and absorption of the pulse weakening the signal amplitude even faster to noise levels when deviating from perfect alignment. An alternative approach consists in relying on reflections from the known surroundings: a lost fuel assembly will block the line-of-sight to the diaphragm resulting in an anomaly in the reflection - either a shorter than expected time-of-flight of the pulse or a complete absence of a reflection like a shadow. In that way, it suffices to align the sensor with the diaphragm instead of the fuel assembly which is much easier to achieve as the robotics on which the sensor is mounted move parallel with the diaphragm. The alignment requirement in the latter approach can be further relaxed by using a tiling of retro-reflectors on the lower surface of the diaphragm. In that way, alignment becomes less vital and the main source of acoustic energy loss - geometric spread of the beam - is almost completely removed, leaving only absorption losses. In this paper, we present the first results in developing a retro reflectance surface for ultrasound in LBE. We present experimental results for different designs of retro-reflectors in both water and LBE. We discuss both linear and array retro-reflectors of different sizes and investigate the influence of the main relevant ultrasonic parameters such as wavelength and spot size on the strength of the received reflection under different alignment angles. We also demonstrate how retro-reflective surfaces can be exploited when localizing objects using linear and rotating scanning methods. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Portable Data-Parallel Surface Reconstruction on a Uniform Rectilinear Grid

With the increasing heterogeneity and on-node parallelism of high-performance computing hardware, a major challenge is to develop portable and efficient algorithms and software. In this work, we present our implementation of a portable code to perform surface reconstruction using NVIDIA's Thrust library. Surface reconstruction is a technique commonly used in volume tracking methods for simulations of multimaterial flow with interfaces. Here, we have designed a 3D mesh data structure that is easily mapped to the 1D vectors used by Thrust and at the same time is simple to use and uses familiar data structure terminology (such as cells, faces, vertices, and edges). With this new data structure in place, we have implemented a piecewise linear interface reconstruction algorithm in 3 dimensions that effectively exploits the symmetry present in a uniform rectilinear computational cell. Finally, we report performance results, which show that a single implementation of these algorithms can be compiled to multiple backends (specifically, multi-core CPUs, NVIDIA GPUs, and Intel Xeon Phi processors), making efficient use of the available parallelism on each. We also compare performance of our implementation to a legacy FORTRAN implementation in Message Passing Interface (MPI) and show performance parity on single and multi-core CPU and achieved good parallel speed-ups on GPU. In conclusion, our research demonstrates the advantage of performance portability of the underlying data-parallel programming model.

97 MATHEMATICS AND COMPUTING↗