Search NASA⌕ Search

SEARCH · Search NASA

Results for “Data Intensive Computing”

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 973 records · Page 54

8760-Based Method for Representing Variable Generation Capacity Value in Capacity Expansion Models: Preprint

Capacity expansion models (CEMs) are widely used to evaluate the least-cost portfolio of electricity generators, transmission, and storage needed to reliably serve demand over the evolution of many years or decades. Various CEM formulations are used to evaluate systems ranging in scale from states or utility service territories to national or multi-national systems. CEMs can be computationally complex, and to achieve acceptable solve times, key parameters are often estimated using simplified methods. In this paper, we focus on two of these key parameters associated with the integration of variable generation (VG) resources: capacity value and curtailment. We first discuss common modeling simplifications used in CEMs to estimate capacity value and curtailment, many of which are based on a representative subset of hours that can miss important tail events or which require assumptions about the load and resource distributions that may not match actual distributions. We then present an alternate approach that captures key elements of chronological operation over all hours of the year without the computationally intensive economic dispatch optimization typically employed within more detailed operational models. The updated methodology characterizes the (1) contribution of VG to system capacity during high load and net load hours, (2) the curtailment level of VG, and (3) the potential reductions in curtailments enabled through deployment of storage and more flexible operation of select thermal generators. We apply this alternate methodology to an existing CEM, the Regional Energy Deployment System (ReEDS). Results demonstrate that this alternate approach provides more accurate estimates of capacity value and curtailments by explicitly capturing system interactions across all hours of the year. This approach could be applied more broadly to CEMs at many different scales where hourly resource and load data is available, greatly improving the representation of challenges associate with integration of variable generation resources.

14 SOLAR ENERGY↗

Structural Modification of Graphene on Copper Substrates Irradiated by Nanosecond High-Intensity Ion Beams

Interaction of a nanosecond high-intensity pulsed ion beam with thin graphene films on copper substrates is analyzed. Methods of Raman spectroscopy are used to investigate the degree of graphene degradation depending on the integral implanted dose. The role of the substrate in the structural degradation of graphene irradiated by charged particle beams is demonstrated using the software package SRIM, intended for modeling radiative defect cascades under irradiation by charged particle beams, and the data on radiation resistance of graphene available in the literature.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Cosmological radio emission induced by WIMP Dark Matter

We present a detailed analysis of the radio synchrotron emission induced by WIMP dark matter annihilations and decays in extragalactic halos. We compute intensity, angular correlation, and source counts and discuss the impact on the expected signals of dark matter clustering, as well as of other astrophysical uncertainties as magnetic fields and spatial diffusion. Bounds on dark matter microscopic properties are then derived, and, depending on the specific set of assumptions, they are competitive with constraints from other indirect dark matter searches. At GHz frequencies, dark matter sources can become a significant fraction of the total number of sources with brightness below the microJansky level. We show that, at this level of fluxes (which are within the reach of the next-generation radio surveys), properties of the faint edge of differential source counts, as well as angular correlation data, can become an important probe for WIMPs.

79 ASTRONOMY AND ASTROPHYSICS↗

Threshold for Onset of Injury in Chinook Salmon from Exposure to Impulsive Pile Driving Sounds

The risk of effects to fishes and other aquatic life from impulsive sound produced by activities such as pile driving and seismic exploration is increasing throughout the world, particularly with the increased exploitation of oceans for energy production. At the same time, there are few data that provide insight into the effects of these sounds on fishes. The goal of this study was to provide quantitative data to define the levels of impulsive sound that could result in the onset of barotrauma to fish. A High Intensity Controlled Impedance Fluid filled wave Tube was developed that enabled laboratory simulation of high-energy impulsive sound that were characteristic of aquatic far-field, plane-wave acoustic conditions. The sounds used were based upon the impulsive sounds generated by an impact hammer striking a steel shell pile. Neutrally buoyant juvenile Chinook salmon (Oncorhynchus tshawytscha) were exposed to impulsive sounds and subsequently evaluated for barotrauma injuries. Observed injuries ranged from mild hematomas at the lowest sound exposure levels to organ hemorrhage at the highest sound exposure levels. Frequency of observed injuries were used to compute a biological response weighted index (RWI) to evaluate the physiological impact of injuries at the different exposure levels. As single strike and cumulative sound exposure levels (SELss, SELcum respectively) increased, RWI values increased. Based on the results, tissue damage associated with adverse physiological costs occurred when the RWI was greater than 2. In terms of sound exposure levels a RWI of 2 was achieved for 1920 strikes by 177 dB re 1 µPa 2 ∙s SEL ss yielding a SEL cum of 210 dB re 1 µPa 2 ∙s, and for 960 strikes by 180 dB re 1 µPa 2 ∙s SEL ss yielding a SEL cum of 210 dB re 1 µPa 2 ∙s. These metrics define thresholds for onset of injury in juvenile Chinook salmon.

59 BASIC BIOLOGICAL SCIENCES↗

On the minimum number of radiation field parameters to specify gas cooling and heating functions

Fast and accurate approximations of gas cooling and heating functions are needed for hydrodynamic galaxy simulations. We use machine learning to analyze atomic gas cooling and heating functions in the presence of a generalized incident local radiation field computed by Cloudy. We characterize the radiation field through binned radiation field intensities instead of the photoionization rates used in our previous work. We find a set of 6 energy bins whose intensities exhibit relatively low correlation. We use these bins as features to train machine learning models to predict Cloudy cooling and heating functions at fixed metallicity. We compare the relative SHapley Additive exPlanation (SHAP) value importance of the features. From the SHAP analysis, we identify a feature subset of 3 energy bins (0.5-1, 1-4, and 13-16Ry) with the largest importance and train additional models on this subset. We compare the mean squared errors and distribution of errors on both the entire training data table and a randomly selected 20% test set withheld from model training. The machine learning models trained with 3 and 6 bins, as well as 3 and 4 photoionization rates, have comparable accuracy everywhere, with errors ≳10 times smaller than for the interpolation table of Gnedin and Hollon (2012). We conclude that 3 energy bins (or 3 analogous photoionization rates: molecular hydrogen photodissociation, neutral hydrogen HI, and fully ionized carbon CVI) are sufficient to characterize the dependence of the gas cooling and heating functions on our assumed incident radiation field model.

79 ASTRONOMY AND ASTROPHYSICS↗

Tokamak-independent Software Analysis Suite for Multi-Spectral Line-Polarization MSE diagnostics

A tokamak-independent analysis suite has been developed to process data from Motional Stark Effect (MSE) diagnostics. The software supports multi-spectral line-polarization MSE diagnostics which simultaneously measure emission at the MSE sigma and pi lines as well as at two `background' wavelengths that are displaced from the MSE spectrum by a few nanometers. This analysis accurately estimates the amplitude of partially-polarized background light at the sigma and pi wavelengths even in situations where the background light changes rapidly in time and space. The signal amplitude at many frequencies is computed using a numerical-beat algorithm which allows the retardance of the MSE photo-elastic modulators (PEMs) to be monitored during routine operation. It also allows the use of summed intensities at multiple frequencies in the calculation of polarization direction, which increases the effective signal strength and reduces sensitivity to PEM retardance drift. The software allows the polarization angles to be corrected for calibration drift using a system that illuminates the MSE diagnostic with polarized light at four known polarization angles within ten seconds of a plasma discharge The software suite is modular, parallelized, and portable to other facilities. (2016-06-06)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tokamak-independent Software Analysis Suite for Multi-Spectral Line-Polarization MSE diagnostics

A tokamak-independent analysis suite has been developed to process data from Motional Stark Effect (MSE) diagnostics. The software supports multi-spectral line-polarization MSE diagnostics which simultaneously measure emission at the MSE sigma and pi lines as well as at two `background' wavelengths that are displaced from the MSE spectrum by a few nanometers. This analysis accurately estimates the amplitude of partially-polarized background light at the sigma and pi wavelengths even in situations where the background light changes rapidly in time and space. The signal amplitude at many frequencies is computed using a numerical-beat algorithm which allows the retardance of the MSE photo-elastic modulators (PEMs) to be monitored during routine operation. It also allows the use of summed intensities at multiple frequencies in the calculation of polarization direction, which increases the effective signal strength and reduces sensitivity to PEM retardance drift. The software allows the polarization angles to be corrected for calibration drift using a system that illuminates the MSE diagnostic with polarized light at four known polarization angles within ten seconds of a plasma discharge The software suite is modular, parallelized, and portable to other facilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Full Integration of Lipidomics Data into Multi-OMIC Functional Enrichment

Lipids have numerous roles critical to living systems. Many of these vital functions are tightly linked or regulated by proteins. Even though the function of many lipid-related enzymes is known, direct links between most individual lipid species and the enzymes that metabolize them are severely lacking. This missing connection impedes biological progress. A resource that allows users to quickly identify protein-lipid relationships will advance our understanding of lipid homeostasis and metabolic perturbed systems. Databases are readily available for mapping small molecule metabolites (i.e., polar metabolites) and specific oxidized fatty acid lipids (e.g., arachidonic acid derivatives) to proteins/enzymes known to be directly involved in metabolism of these molecules (e.g. SMPDB, Reactome, HMDB, BioCyc). This enables the integration of metabolite and protein (as well as transcript) data when performing functional enrichment and metabolic modeling allowing researchers to rapidly examine metabolic pathways from multiple omics sources. Currently, linking lipidomics data with proteomics data is conducted using a time-intensive, user defined, manual process. Here we propose creating a tool that will link lipidomics and proteomics data and also conduct functional enrichment analysis and visualization allowing for full statistical integration. Using text-mining approaches, we will extract lipid references from uniprot and create a database linking lipid species to enzymes directly related to their metabolism. We will also extract directional annotation, which will allow even greater specificity when linking data types and will help eliminate false positives from functional associations.

59 BASIC BIOLOGICAL SCIENCES↗

GlyQ-IQ: Glycomics Quintavariate-Informed Quantification with High-Performance Computing and GlycoGrid 4D Visualization

Dense LC-MS datasets have convoluted extracted ion chromatograms with multiple chromatographic peaks that cloud the differentiation between intact compounds with their overlapping isotopic distributions, peaks due to insource ion fragmentation, and noise. Making this differentiation is critical in glycomics datasets because chromatographic peaks correspond to different intact glycan structural isomers. The GlyQ-IQ software is targeted chromatography centric software designed for chromatogram and mass spectra data processing and subsequent glycan composition annotation. The targeted analysis approach offers several key advantages to LC-MS data processing and annotation over traditional algorithms. A priori information about the individual target’s elemental composition allows for exact isotope profile modeling for improved feature detection and increased sensitivity by focusing chromatogram generation and peak fitting on the isotopic species in the distribution having the highest intensity and data quality. Glycan target annotation is corroborated by glycan family relationships and in source fragmentation detection. The GlyQ-IQ software is developed in this work (Part 1) and was used to profile N-glycan compositions from human serum LC-MS Datasets. The companion manuscript GlyQ-IQ Part 2 discusses developments in human serum N-glycan sample preparation, glycan isomer separation, and glycan electrospray ionization. A case study is presented to demonstrate how GlyQ-IQ identifies and removes confounding chromatographic peaks from high mannose glycan isomers from human blood serum. In addition, GlyQ-IQ was used to generate a broad N-glycan profile from a high resolution (100K/60K) nESI-LS-MS/MS dataset including CID and HCD fragmentation acquired on a Velos Pro Mass spectrometer. 101 glycan compositions and 353 isomer peaks were detected from a single sample. 99% of the GlyQ-IQ glycan-feature assignments passed manual validation and are backed with high resolution mass spectra and mass accuracies less than 7 ppm.

Environmental Molecular Sciences Laboratory↗

Upgrade of the Minos+ Experiment Data Acquisition for the High Energy NuMI Beam Run

The Minos+ experiment is an extension of the Minos experiment at a higher energy and more intense neutrino beam, with the data collection having begun in the fall of 2013. The neutrino beam is provided by the Neutrinos from the Main Injector (NuMI) beam-line at Fermi National Accelerator Laboratory (Fermilab). The detector apparatus consists of two main detectors, one underground at Fermilab and the other in Soudan, Minnesota with the purpose of studying neutrino oscillations at a base line of 735 km. The original data acquisition system has been running for several years collecting data from NuMI, but with the extended run from 2013, parts of the system needed to be replaced due to obsolescence, reliability problems, and data throughput limitations. Specifically, we have replaced the front-end readout controllers, event builder, and data acquisition computing and trigger processing farms with modern, modular and reliable devices with few single points of failure. The new system is based on gigabit Ethernet TCP/IP communication to implement the event building and concatenation of data from many front-end VME readout crates. The simplicity and partitionability of the new system greatly eases the debugging and diagnosing process. As a result, the new system improves throughput by about a factor of three compared to the old system, up to 800 megabits per second, and has proven robust and reliable in the current run.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The feasibility of head motion tracking in helical CT: A step toward motion correction

Purpose: To establish a practical and accurate motion tracking method for the development of rigid motion correction methods in helical x-ray computed tomography (CT). Methods: A commercially available optical motion tracking system provided 6 degrees of freedom pose measurements at 60 Hz. A 4 Multiplication-Sign 4 calibration matrix was determined to convert raw pose data acquired in tracker coordinates to a fixed CT coordinate system with origin at the isocenter of the scanner. Two calibration methods, absolute orientation (AO), and a new method based on image registration (IR), were compared by means of landmark analysis and correlation coefficient in phantom images coregistered using the derived motion transformations. Results: Transformations calculated using the IR-derived calibration matrix were found to be more accurate, with positional errors less than 0.5 mm (mean RMS), and highly correlated image voxel intensities. The AO-derived calibration matrix yielded larger mean RMS positional errors ( Asymptotically-Equal-To 1.0 mm), and poorer correlation coefficients. Conclusions: The authors have demonstrated the feasibility of accurate motion tracking for retrospective motion correction in helical CT. Their new IR-based calibration method based on image registration and function minimization was simpler to perform and delivered more accurate calibration matrices. This technique is a useful tool for future work on rigid motion correction in helical CT and potentially also other imaging modalities.

60 APPLIED LIFE SCIENCES↗

Precursor Analysis Report: SQL Slammer Worm Infection of Davis-Besse Nuclear Power Plant 2003

The SQL Slammer Worm Infection of Davis-Besse Nuclear Power Plant 2003 Precursor Analysis Report leverages publicly available information about Davis-Besse’s 2003 cyber attack and catalogs anomalous observables for each technique employed in the attack. This analysis is based upon the methodology of the Cybersecurity for the Operational Technology Environment (CyOTE) program. On 25 January 2003, the SQL Slammer worm infected more than 90% of vulnerable hosts and crashed the internet in 10 to 15 minutes, making it one of the fastest spreading worms in history. SQL Slammer is a fileless, memory-resident worm that remotely exploits a stack-based buffer overflow vulnerability on local hosts to intensively scan and rapidly self-propagate across the internet. The worm infected approximately 300,000 unpatched hosts running Microsoft Structured Query Language (SQL) Server 2000 or Microsoft Desktop Engine (MSDE) 2000 with SQL Server Resolution Service. The SQL Slammer worm indirectly infected FirstEnergy’s Davis-Besse nuclear power plant by first infecting a consultant’s company network server and then propagating through an external misconfigured connection into Davis-Besse’s site network. The infection caused major network congestion, slow performance, data overloads, and the inability of local hosts to communicate with each other, which eventually caused a loss of availability and a loss of view when the Safety Parameter Display System (SPDS) and Plant Process Computer (PPC) crashed. At the time of the infection, the plant was already offline, the digital monitoring systems had redundant analog backups, and the plant control and safety functions were not affected, so there were no concerns of a safety breach. However, this incident resulted in many lessons learned and spawned important discussions about cybersecurity’s role in nuclear safety and electric power reliability regulation, policy, and guidance. Researchers and analysts identified 10 unique techniques utilized during the attack with a total of 640 observables using MITRE ATT&CK® for Industrial Control Systems. The CyOTE program assesses observables accompanying techniques used prior to the triggering event to identify opportunities to detect malicious activity. If observables accompanying the attack techniques are perceived and investigated prior to the triggering event, earlier comprehension of malicious activity can take place. Eight of the identified techniques used during Davis-Besse cyber attack were precursors to the triggering event. Analysis identified 596 observables associated with these precursor techniques, 428 of which were assessed to have an increased likelihood of being perceived in the 331 days preceding the triggering event. The response and comprehension time could have been reduced if the observables had been identified earlier. The information gathered in this report contributes to a library of observables tied to a repository of artifacts, data sources, and technique detection references for practitioners and developers to support the comprehension of indicators of attack. Asset owners and operators can use these products if they experience similar observables or to prepare for comparable scenarios.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Improvement in the Accuracy of Flux Measurement of Radio Sources by Exploiting an Arithmetic Pattern in Photon Bunching Noise

A hierarchy of statistics of increasing sophistication and accuracy is proposed to exploit an interesting and fundamental arithmetic structure in the photon bunching noise of incoherent light of large photon occupation number, with the purpose of suppressing the noise and rendering a more reliable and unbiased measurement of the light intensity. The method does not require any new hardware, rather it operates at the software level with the help of high-precision computers to reprocess the intensity time series of the incident light to create a new series with smaller bunching noise coherence length. The ultimate accuracy improvement of this method of flux measurement is limited by the timing resolution of the detector and the photon occupation number of the beam (the higher the photon number the better the performance). The principal application is accuracy improvement in the signal-limited bolometric flux measurement of a radio source.

79 ASTRONOMY AND ASTROPHYSICS↗

Star-forming galaxies as the origin of diffuse high-energy backgrounds: gamma-ray and neutrino connections, and implications for starburst history

Star-forming galaxies have been predicted to contribute considerably to the diffuse gamma-ray background as they are guaranteed reservoirs of cosmic rays. Assuming that the hadronic interactions responsible for high-energy gamma rays also produce high-energy neutrinos and that O(100) PeV cosmic rays can be produced and confined in starburst galaxies, we here discuss the possibility that star-forming galaxies are also the main sources of the high-energy neutrinos observed by the IceCube experiment. First, we compute the diffuse gamma-ray background from star-forming galaxies, adopting the latest Herschel PEP/HerMES luminosity function and relying on the correlation between the gamma-ray and infrared luminosities reported by Fermi observations. Then we derive the expected intensity of the diffuse high-energy neutrinos from star-forming galaxies including normal and starburst galaxies. Our results indicate that starbursts, including those with active galactic nuclei and galaxy mergers, could be the main sources of the high-energy neutrinos observed by the IceCube experiment. We find that assuming a cosmic-ray spectral index of 2.1–2.2 for all starburst-like galaxies, our predictions can be consistent with both the Fermi and IceCube data, but larger indices readily fail to explain the observed diffuse neutrino flux. Taking the starburst high-energy spectral index as free parameter, and extrapolating from GeV to PeV energies, we find that the spectra harder than E{sup -2.15} are likely to be excluded by the IceCube data, which can be more constraining than the Fermi data for this population.

79 ASTRONOMY AND ASTROPHYSICS↗

Fundamental Mechanisms Driving the Amorphous to Crystalline Phase Transformation

Phase transformations are ubiquitous, fundamental phenomena that lie at the heart of many structural, optical and electronic properties in condensed matter physics and materials science. Many transformations, especially those occurring under extreme conditions such as rapid changes in the thermodynamic state, are controlled by poorly understood processes involving the nucleation and quenching of metastable phases. Typically these processes occur on time and length scales invisible to most experimental techniques ({micro}s and faster, nm and smaller), so our understanding of the dynamics tends to be very limited and indirect, often relying on simulations combined with experimental study of the ''time infinity'' end state. Experimental techniques that can directly probe phase transformations on their proper time and length scales are therefore key to providing fundamental insights into the whole area of transformation physics and materials science. LLNL possesses a unique dynamic transmission electron microscope (DTEM) capable of taking images and diffraction patterns of laser-driven material processes with resolution measured in nanometers and nanoseconds. The DTEM has previously used time-resolved diffraction patterns to quantitatively study phase transformations that are orders of magnitude too fast for conventional in situ TEM. More recently the microscope has demonstrated the ability to directly image a reaction front moving at {approx}13 nm/ns and the nucleation of a new phase behind that front. Certain compound semiconductor phase change materials, such as Ge{sub 2}Sb{sub 2}Te{sub 5} (GST), Sb{sub 2}Te and GeSb, exhibit a technologically important series of transformations on scales that fall neatly into the performance specifications of the DTEM. If a small portion of such material is heated above its melting point and then rapidly cooled, it quenches into an amorphous state. Heating again with a less intense pulse leads to recrystallization into a vacancy-stabilized metastable rock salt structure. Each transformation takes {approx}10-100 ns, and the cycle can be driven repeatedly a very large number of times with a nanosecond laser such as the DTEM's sample drive laser. These materials are widely used in optical storage devices such as rewritable CDs and DVDs, and they are also applied in a novel solid state memory technology - phase change memory (PCM). PCM has the potential to produce nonvolatile memory systems with high speed, extreme density, and very low power requirements. For PCM applications several materials properties are of great importance: the resistivities of both phases, the crystallization temperature, the melting point, the crystallization speed, reversibility (number of phase-transformation cycles without degradation) and stability against crystallization at elevated temperature. For a viable technology, all these properties need to have good scaling behavior, as dimensions of the memory cells will shrink with every generation. In this LDRD project, we used the unique single-shot nanosecond in situ experimentation capabilities of the DTEM to watch these transformations in GST on the time and length scales most relevant for device applications. Interpretation of the results was performed in conjunction with atomistic and finite-element computations. Samples were provided by collaborators at IBM and Stanford University. We observed, and measured the kinetics of, the amorphous-crystalline and melting-solidification transitions in uniform thin-film samples. Above a certain threshold, the crystal nucleation rate was found to be enormously high (with many nuclei appearing per cubic {micro}m even after nanosecond-scale incubation times), in agreement with atomistic simulation and consistent with an extremely low nucleation barrier. We developed data reduction techniques based on principal component analysis (PCA), revealing the complex, multi-dimensional evolution of the material while suppressing noise and irrelevant information. Using a novel specimen geometry, we also achieved repeated switching between the amorphous and crystalline phases enabling in situ study of structural change after phase cycling, which is relevant to device performance. We also observed the coupling between the phase transformations and the evolution of morphology on the nanometer scale, revealing the gradual development of striations in uniform films and preferential melting at sharp edges in laser-heated nanopatterned GST. This nonuniform melting, interpreted through simulation as being a direct result of geometrical laser absorption effects, appears to be responsible for a marked loss in morphological stability even at moderate laser intensities and may be an important factor in the longevity of nanostructured phase change materials in memory applications.

36 MATERIALS SCIENCE↗

Role of excited state photoionization in the 852.1 nm Cs laser pumped by Cs-Ar photoassociation

Photoionization of Cs (6p {sup 2}P{sub 3/2}) atoms during the operation of a Cs D{sub 2} line (852.1 nm: 6p {sup 2}P{sub 3/2}{yields}6s {sup 2}S{sub 1/2}) laser, pumped by free{yields}free transitions of thermal Cs-Ar ground state pairs, has been investigated experimentally and computationally. Photoexcitation of Cs vapor/Ar mixtures through the blue satellite of the D{sub 2} transition (peaking at 836.7 nm) selectively populates the {sup 2}P{sub 3/2} upper laser level by the dissociation of the CsAr excited complex. Comparison of laser output energy data, for instantaneous pump powers up to 3 MW, with the predictions of a numerical model sets an upper bound of 8 Multiplication-Sign 10{sup -26} cm{sup 4} W{sup -1} on the Cs (6p {sup 2}P{sub 3/2}) two photon ionization cross-section at 836.7 nm which corresponds to a single photon cross-section of 2.4 Multiplication-Sign 10{sup -19} cm{sup 2} for a peak pump intensity of 3 MW cm{sup -2}.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High Energy Electromagnetic Particle Transportation on the GPU

We present massively parallel high energy electromagnetic particle transportation through a finely segmented detector on a Graphics Processing Unit (GPU). Simulating events of energetic particle decay in a general-purpose high energy physics (HEP) detector requires intensive computing resources, due to the complexity of the geometry as well as physics processes applied to particles copiously produced by primary collisions and secondary interactions. The recent advent of hardware architectures of many-core or accelerated processors provides the variety of concurrent programming models applicable not only for the high performance parallel computing, but also for the conventional computing intensive application such as the HEP detector simulation. The components of our prototype are a transportation process under a non-uniform magnetic field, geometry navigation with a set of solid shapes and materials, electromagnetic physics processes for electrons and photons, and an interface to a framework that dispatches bundles of tracks in a highly vectorized manner optimizing for spatial locality and throughput. Core algorithms and methods are excerpted from the Geant4 toolkit, and are modified and optimized for the GPU application. Program kernels written in C/C++ are designed to be compatible with CUDA and OpenCL and with the aim to be generic enough for easy porting to future programming models and hardware architectures. To improve throughput by overlapping data transfers with kernel execution, multiple CUDA streams are used. Issues with floating point accuracy, random numbers generation, data structure, kernel divergences and register spills are also considered. Performance evaluation for the relative speedup compared to the corresponding sequential execution on CPU is presented as well.

Canal, P.↗