Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 685 records · Page 38

Implementing intensity modulated radiotherapy to the prostate bed: Dosimetric study and early clinical results

Salvage intensity modulated radiotherapy (IMRT) to the prostate bed has hardly been studied so far. We present here a feasibility study and early clinical results for 10 patients. These patients were selected on the basis of having either a biochemical relapse or high risk histology after prostatectomy. They were treated using “sliding-window” IMRT to 68 Gy in 34 fractions. Three-dimensional conformal radiotherapy (3D-CRT) plans were generated using the same planning computed tomography data set. Dose coverage of planning target volumes (PTVs) and of organs-at-risk (OAR, namely: rectum, bladder, and femoral heads) were compared. Acute toxicity and chronic toxicity were measured using the Common Toxicity Criteria for Adverse Events version 3.0 scale. IMRT significantly reduces the dose above the prescription dose given to the PTV1 (mean dose: IMRT 67.2 Gy vs 3D-CRT 67.7 Gy (p = 0.0137)), without altering dose coverage for PTV2 (mean dose: IMRT 68.1 Gy vs 3D-CRT 68.0 Gy (p = 0.3750)). Doses to OAR were lower with IMRT and differences were statistically significant (mean dose: IMRT 51.4 Gy vs 3D-CRT 56.6 Gy for rectum (p = 0.002), IMRT 45.1 Gy vs 3D-CRT 53.1 Gy for bladder (p = 0.002), and IMRT 26.1 Gy vs 3D-CRT 28.4 Gy for femoral heads (p = 0.0059)). There was no acute or chronic genitourinary or gastrointestinal toxicity >1 with a median follow-up of 38 months. IMRT to the prostatic fossa is feasible and reduces dose to OAR, with consequential limited toxicity.

61 RADIATION PROTECTION AND DOSIMETRY↗

Interior-Point Methods for Estimating Seasonal Parameters in Discrete-Time Infectious Disease Models

Infectious diseases remain a significant health concern around the world. Mathematical modeling of these diseases can help us understand their dynamics and develop more effective control strategies. In this work, we show the capabilities of interior-point methods and nonlinear programming (NLP) formulations to efficiently estimate parameters in multiple discrete-time disease models using measles case count data from three cities. These models include multiplicative measurement noise and incorporate seasonality into multiple model parameters. Our results show that nearly identical patterns are estimated even when assuming seasonality in different model parameters, and that these patterns show strong correlation to school term holidays across very different social settings and holiday schedules. We show that interior-point methods provide a fast and flexible approach to parameterizing models that can be an alternative to more computationally intensive methods.

60 APPLIED LIFE SCIENCES↗

Sequential Inference of Hospitalization Electronic Health Records Using Probabilistic Models

In the dynamic hospital setting, decision support can be a valuable tool for improving patient outcomes. Data-driven inference of future outcomes is challenging in this dynamic setting, where long sequences such as laboratory tests and medications are updated frequently. This is due in part to heterogeneity of data types and mixed-sequence types contained in variable length sequences. In this work we design a probabilistic unsupervised model for multiple arbitrary-length sequences contained in hospitalization Electronic Health Record (EHR) data. The model uses a latent variable structure and captures complex relationships between medications, diagnoses, laboratory tests, neurological assessments, and medications. It can be trained on original data, without requiring any lossy transformations or time binning. Inference algorithms are derived that use partial data to infer properties of the complete sequences, including their length and presence of specific values. We train this model on data from subjects receiving medical care in the Kaiser Permanente Northern California integrated healthcare delivery system. The results are evaluated against held-out data for predicting the length of sequences and presence of Intensive Care Unit (ICU) in hospitalization bed sequences. Our method outperforms a baseline approach, showing that in these experiments the trained model captures information in the sequences that is informative of their future values.

97 MATHEMATICS AND COMPUTING↗

Data Fusion for a Vision-Radiological System for Source Tracking and Discovery

A multidisciplinary approach to allow the tracking of the movement of radioactive sources by fusing data from multiple radiological and visual sensors is under development. The goal is to improve the ability to detect, locate, track and identify nuclear/radiological threats. The key concept is that such widely available visual and depth sensors can impact radiological detection, since the intensity fall-off in the count rate can be correlated to movement in three dimensions. To enable this, we pose an important question; what is the right combination of sensing modalities and vision algorithms that can best compliment a radiological sensor, for the purpose of detection and tracking of radioactive material? Similarly what is the best radiation detection methods and unfolding algorithms suited for data fusion with tracking data? Data fusion of multi-sensor data for radiation detection have seen some interesting developments lately. Significant examples include intelligent radiation sensor systems (IRSS), which are based on larger numbers of distributed similar or identical radiation sensors coupled with position data for network capable to detect and locate radiation source. Other developments are gamma-ray imaging systems based on Compton scatter in segmented detector arrays. Similar developments using coded apertures or scatter cameras for neutrons have recently occurred. The main limitation of such systems is not so much in their capability but rather in their complexity and cost which is prohibitive for large scale deployment. Presented here is a fusion system based on simple, low-cost computer vision and radiological sensors for tracking of multiple objects and identifying potential radiological materials being transported or shipped. The main focus of this work is the development on two separate calibration algorithms for characterizing the fused sensor system. The deviation from a simple inverse square-root fall-off of radiation intensity is explored and accounted for. In particular, the computer vision system enables a map of distance-dependence of the sources being tracked. Infrared, laser or stereoscopic vision sensors are all options for computer-vision implementation depending on interior vs exterior deployment, resolution desired and other factors. Similarly the radiation sensors will be focused on gamma-ray or neutron detection due to the long travel length and ability to penetrate even moderate shielding. There is a significant difference between the vision sensors and radiation sensors in the way the 'source' or signals are generated. A vision sensor needs an external light-source to illuminate the object and then detects the re-emitted illumination (or lack thereof). However, for a radiation detector, the radioactive material is the source itself. The only exception to this is the field of active interrogations where radiation is beamed into a material to entice new/additional radiation emission beyond what the material would emit spontaneously. The aspect of the nuclear material being the source itself means that all other objects in the environment are 'illuminated' or irradiated by the source. Most radiation will readily penetrate regular material, scatter in new directions or be absorbed. Thus if a radiation source is located near a larger object that object will in turn scatter some radiation that was initially emitted in a direction other than the direction of the radiation detector, this can add to the count rate that is observed. The effect of these scatter is a deviation from the traditional distance dependence of the radiation signal and is a key challenge that needs a combined system calibration solution and algorithms. Thus both an algebraic approach as well as a statistical approach have been developed and independently evaluated to investigate the sensitivity to this deviation from the simplified radiation fall-off as a function of distance. The resulting calibrated system algorithms are used and demonstrated in various laboratory scenarios, and later in realistic tracking scenarios. The selection and testing of radiological and computer-vision sensors for the additional specific scenarios will be the subject of ongoing and future work. (authors)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

SU-E-QI-14: Quantitative Variogram Detection of Mild, Unilateral Disease in Elastase-Treated Rats

Purpose: Determining the presence of mild or early disease in the lungs can be challenging and subjective. We present a rapid and objective method for evaluating lung damage in a rat model of unilateral mild emphysema based on a new approach to heterogeneity assessment. We combined octree decomposition (used in three-dimensional (3D) computer graphics) with variograms (used in geostatistics to assess spatial relationships) to evaluate 3D computed tomography (CT) lung images for disease. Methods: Male, Sprague-Dawley rats (232 ± 7 g) were intratracheally dosed with 50 U/kg of elastase dissolved in 200 μL of saline to a single lobe (n=6) or with saline only (n=5). After four weeks, 3D micro-CT images were acquired at end expiration on mechanically ventilated rats using prospective gating. Images were masked, and lungs were decomposed to homogeneous blocks of 2×2×2, 4×4×4, and 8×8×8 voxels using octree decomposition. The spatial variance – the square of the difference of signal intensity – between all pairs of the 8×8×8 blocks was calculated. Variograms – graphs of distance vs. variance - were made, and data were fit to a power law and the exponent determined. The mean HU values, coefficient of variation (CoV), and the emphysema index (EI) were calculated and compared to the variograms. Results: The variogram analysis showed that significant differences between groups existed (p<0.01), whereas the mean HU (p=0.07), CoV (p=0.24), and EI (p=0.08) did not. Calculation time for the variogram for a typical 1000 block decomposition was ∼6 seconds, and octree decomposition took ∼2 minutes. Decomposing the images prior to variogram calculation resulted in a ∼700x decrease in time as compared to other published approaches. Conclusions: Our results suggest that the approach combining octree decomposition and variogram analysis may be a rapid, non-subjective, and sensitive imaging-based biomarker for quantitative characterization of lung disease.

60 APPLIED LIFE SCIENCES↗

Advanced Manufacturing for Thermal and Environmental Control Systems: Achieving National Energy Goals

This project was part of a regional initiative in the five counties of Central New York (CNY) that received funding from the U.S. Department of Energy (DOE) and four other federal agencies through the 2012 Advanced Manufacturing Jobs and Innovation Accelerator Challenge (AMJIAC). The CNY initiative was focused on cultivating the emergent regional cluster in “Advanced Manufacturing for Thermal and Environmental Control (AM-TEC).” As one component of the CNY AM-TEC initiative, the DOE-funded project supported five research & development seed projects that strategically targeted: 1) needs and opportunities of CNY AM-TEC companies, and 2) the goal of DOE’s Advanced Manufacturing Office (AMO) to reduce energy consumption by 50% across product life-cycles over 10 years. The project also sought to fulfill the AMO mission of developing and demonstrating new, energy-efficient processing and materials technologies at a scale adequate to prove their value to manufacturers and spur investment. The five seed projects demonstrated technologies and processes that can reduce energy intensity and improve production as well as use less energy throughout their lifecycles. The project was conducted over three years in two 18-month budget periods. During the first budget period, two projects proposed in the original AMJAIC application were successfully completed: Seed Project 1 focused on saving energy in heat transfer processes via development of nano structured surfaces to significantly increase heat flux; Seed Project 2 addressed saving energy in data centers via subzero cooling of the computing processors. Also during the first budget period, a process was developed and executed to select a second round of seed projects via a competitive request for proposals from regional companies and university collaborators. Applicants were encouraged to form industry-academic partnerships to leverage experience and resources of public and private sectors in the CNY region. Proposals were evaluated by a national panel of experts. Three projects were selected for awards and were completed successfully during the second budget period: Seed Project 3 focused on enabling self-powered furnaces to permit residents to shelter in place during power outages; Seed Project 4 addressed development of a novel method of controlling air conditioning systems that could enable flexible load matching in market segments not possible with existing technologies; and Seed Project 5 focused on the creation of smarter occupancy sensors to enable effective highly localized demand based ventilation.

30 DIRECT ENERGY CONVERSION↗

Effectively Using Remote I/O For Work Composition in Distributed Workflows

Distributed scientific workflows are becoming more important with the interest in incorporating AI into their loops. A critical programming and performance question is how to compose workflow tasks when data is produced on one system but must be consumed on another. Since the dominant technique is composition with remote I/O, this paper explores its performance expectations. We describe BigFlowSim, a workflow I/O simulator that captures key implementation choices for remote I/O, including intensity, reuse, locality, access pattern, and data movement.With BigFlowSim, we generate a synthetic benchmark. We quantify the effects of each parameter with a performance sensitivity study. We explain trends in terms of data movement reduction and show that, under certain conditions, it is possible to establish a total order among most parameters. We apply these insights to a high energy physics workflow, Belle II Monte Carlo and simulate several I/O optimizations. Speedups range from 5% to 2×, without changing compute time.

Friese, Ryan D.↗

An Advanced Network and distributed Storage Laboratory (ANDSL) for Data Intensive Science (Final Report)

The original intent of this project was to build and operate an Advanced Network and Distributed Storage Laboratory (ANDSL) for Data Intensive Science that will prepare the Open Science Grid (OSG) community for a new generation of wide area communication capabilities operating at a 100Gb rate. The plan was to engage a team of four dedicated scientists that will design, develop and operate a state of the art laboratory that will be fully integrated into the organizational and operational infrastructure of the OSG. The ANDSL was planned to include host-based hardware and software required to profile, emulate and study real-life distributed storage workloads at rates of 100 Gbps. The laboratory was designed to utilize the Department Of Energy’s Advanced Network Initiative (ANI) 100GbTestbed as connections to the fabric become available so that OSG stakeholders can use it to profile real-life data intensive science applications as well as experiment with advanced storage and data management technologies. The goal of the original program of work was to close the adoption gap, decreasing the time between availability of production 100Gbps links and the end-to-end use of such a network to enable scientific discovery.

100Gb networks↗

Development, verification, and validation of comprehensive acoustic fluid-structure interaction capabilities in an open-source computational platform

The acoustic fluid-structure interaction (FSI) formulation is a practical numerical approach for the seismic analysis of fluid-filled tanks. However, there are no verification and validation studies reported in the literature that demonstrate the ability of an acoustic FSI numerical model to predict responses important to structural and mechanical design for intense translational and rotational earthquake inputs. Herein, an acoustic FSI formulation is implemented in the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE), and is formally verified and validated using analytical solutions and code-to-code verification, and experimental data, respectively. The analytical solutions are for small amplitude, unidirectional seismic inputs. The code-to-code verification utilizes a previously verified and validated Arbitrary Lagrangian-Eulerian (ALE) numerical model in the commercial finite element code LS-DYNA. The validation studies utilize a comprehensive data set assembled from results of 3D earthquake-simulator tests of a fluid-filled vessel. The acoustic numerical model in MOOSE is verified and validated for hydrodynamic pressures and support reactions except for cases that involve significant convective response. For small amplitude inputs, numerically predicted wave heights match those of the analytical solutions. The numerical model is not verified and validated for wave height calculations under intense 3D seismic inputs. The run times for the acoustic FSI simulations in MOOSE are an order of magnitude, or more, shorter than for the corresponding ALE simulations in LS-DYNA. The utility of the MOOSE acoustic FSI implementation is demonstrated by seismic analysis of a building equipped with a fluid-filled, advanced nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

ON A NEW THEORETICAL FRAMEWORK FOR RR LYRAE STARS. I. THE METALLICITY DEPENDENCE

We present new nonlinear, time-dependent convective hydrodynamical models of RR Lyrae stars computed assuming a constant helium-to-metal enrichment ratio and a broad range in metal abundances (Z = 0.0001–0.02). The stellar masses and luminosities adopted to construct the pulsation models were fixed according to detailed central He-burning horizontal-branch evolutionary models. The pulsation models cover a broad range in stellar luminosity and effective temperatures and the modal stability is investigated for both fundamental (FU) and first overtone polsators (FOs). We predict the topology of the instability strip (IS) as a function of the metal content and new analytical relations for the edges of the IS in the observational plane. Moreover, a new analytical relation to constrain the pulsation mass of double pulsators as a function of the period ratio and the metal content is provided. We derive new Period–Radius–Metallicity relations for FU and FO pulsators. They agree quite well with similar empirical and theoretical relations in the literature. From the predicted bolometric light curves, transformed into optical (UBVRI) and near-infrared (NIR; JHK) bands, we compute the intensity-averaged mean magnitudes along the entire pulsation cycle and in turn new and homogenous metal-dependent (RIJHK) Period–Luminosity relations. Moreover, we compute new dual and triple-band optical, optical–NIR, and NIR Period–Wesenheit–Metallicity relations. Interestingly, we find that the optical Period-W(V, B–V) is independent of the metal content and that the accuracy of individual distances is a balance between the adopted diagnostics and the precision of photometric and spectroscopic data sets.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulating Blade-Strike on Fish passing through Marine Hydrokinetic Turbines

The study reported here evaluated the occurrence, frequency, and intensity of blade strike of fish on an axial-flow marine hydrokinetic turbine by using two modeling approaches: a conventional kinematic formulation and a proposed Lagrangian particle- based scheme. The kinematic model included simplifying assumptions of fish trajectories such as distribution and velocity. The proposed method overcame the need for such simplifications by integrating the following components into a computational fluid dynamics (CFD) model: (i) advanced eddy-resolving flow simulation, (ii) generation of ambient turbulence based on field data, (iii) moving turbine blades in highly transient flows, and (iv) Lagrangian particles to mimic the potential fish pathways. The test conditions to evaluate the blade-strike probability and fish survival rate were: (i) the turbulent environment, (ii) the fish size, and (iii) the approaching flow velocity. The proposed method offered the ability to produce potential fish trajectories and their interaction with the rotating turbine. Depending upon the scenario, the percentile of particles that registered a collision event ranged from 6% to 19% of the released sample size. Next, by using a set of experimental correlations of the exposure-response of living fish colliding with moving blades, the simulated collision data were used as input variables to estimate the survival rate of fish passing through the operating turbine. The resulting survival rates were greater than 96% in all scenarios, which is comparable to or better than known survival rates for conventional hydropower turbines. The figures of strike probability and mortality rate were amplified by the kinematic model. The proposed method offered the advantage of expanding the evaluation of other mechanisms of stress and injury on fish derived from hydrokinetic turbines and related devices.

CFD↗

Automatic Lane-Level Road Network Extraction from Aerial Imagery for Transportation Digital Twins

Accurate road networks are essential for credible traffic microsimulation and transportation digital twins, yet high-definition maps are often difficult to obtain due to limited availability, high cost, or proprietary restrictions. Some build networks from crowdsourced data, such as OpenStreetMap, but these sources often contain geometric and semantic inconsistencies. Others create networks manually, a process that is labor-intensive and difficult to scale. To address these limitations, this work presents an end-to-end pipeline that automatically extracts georeferenced, lane-level road networks from publicly available high-resolution satellite imagery and converts them into simulation-ready assets. The developed end-to-end pipeline has three primary modules: (1) A computer-vision-based module first detects directed lane geometries and intersection layouts. (2) A heuristic-based topology construction module then identifies approach and exit legs and establishes conflict-free lane-to-lane connections. (3) Finally, an automatic simulation-building module converts the extracted network into standard formats, e.g., OpenDRIVE, and generates routable SUMO networks. The framework supports both complete network construction from scratch and local-scale refinement of existing networks through lane-count correction, transition recovery, and geometric regularization. The proposed pipeline provides a practical pathway to generate traffic simulation networks from satellite imagery, significantly reducing manual reconstruction effort and enabling scalable, continuously updated transportation digital twins.

Guo, Hetian [University of Georgia, Athens] (ORCID↗

Beam Based RF Voltage Measurements and Longitudinal Beam Tomography at the Fermilab Booster

Increasing proton beam power on neutrino production targets is one of the major goals of the Fermilab long term accelerator programs. In this effort, the Fermilab 8 GeV Booster synchrotron plays a critical role for at least the next two decades. Therefore, understanding the Booster in great detail is important as we continue to improve its performance. For example, it is important to know accurately the available RF power in the Booster by carrying out beam-based measurements in order to specify the needed upgrades to the Booster RF system. Since the Booster magnetic field is changing continuously measuring/calibrating the RF voltage is not a trivial task. Here, we present a beam based method for the RF voltage measurements. Data analysis is carried out using computer programs developed in Python and MATLAB. The method presented here is applicable to any RCS which do not have flat-bottom and flat-top in the acceleration magnetic ramps. We have also carried out longitudinal beam tomography at injection and extraction energies with the data used for RF voltage measurements. Beam based RF voltage measurements and beam tomography were never done before for the Fermilab Booster. The results from these investigations will be very useful in future intensity upgrades.

43 PARTICLE ACCELERATORS↗

Broadband (0–4 Hz) Ground Motions for a Magnitude 7.0 Hayward Fault Earthquake With Three‐Dimensional Structure and Topography

Abstract We performed fully deterministic broadband (0–4 Hz) high‐performance computing ground motion simulations of a magnitude 7.0 scenario earthquake on the Hayward Fault (HF) in the San Francisco Bay Area of Northern California. Simulations consider average one‐dimensional (1‐D) and three‐dimensional (3‐D) anelastic structure with flat and topographic free surfaces. Ground motion intensity measures (GMIMs) for the 3‐D model display dramatic differences across the HF due to geologic heterogeneity, with low wave speeds east of the HF amplifying motions. The median GMIMs agree well with Ground Motion Prediction Equations (GMPEs); however, the 3‐D model generates more scatter than the 1‐D model. Ratios of 3‐D/1‐D GMIMs from the same source allow isolation of path and site effects for the 3‐D model. These ratios show remarkably similar trends as site‐specific factors for the GMPE predictions, suggesting that wave propagation effects in our 3‐D simulations are on average consistent with empirical data.

58 GEOSCIENCES↗

Fiber optic computing using distributed feedback

Abstract The widespread adoption of machine learning and other matrix intensive computing algorithms has renewed interest in analog optical computing, which has the potential to perform large-scale matrix multiplications with superior energy scaling and lower latency than digital electronics. However, most optical techniques rely on spatial multiplexing, requiring a large number of modulators and detectors, and are typically restricted to performing a single kernel convolution operation per layer. Here, we introduce a fiber-optic computing architecture based on temporal multiplexing and distributed feedback that performs multiple convolutions on the input data in a single layer. Using Rayleigh backscattering in standard single mode fiber, we show that this technique can efficiently apply a series of random nonlinear projections to the input data, facilitating a variety of computing tasks. The approach enables efficient energy scaling with orders of magnitude lower power consumption than GPUs, while maintaining low latency and high data-throughput.

97 MATHEMATICS AND COMPUTING↗

Solvent effects on the absorption spectrum and first hyperpolarizability of keto-enol tautomeric forms of anil derivatives: A Monte Carlo/quantum mechanics study

Theoretical results for the absorption spectrum and electric properties of the enol and keto tautomeric forms of anil derivatives in the gas-phase and in solution are presented. The electronic properties in chloroform, acetonitrile, methanol, and water were determined by carrying out sequential Monte Carlo simulations and quantum mechanics calculations based on the time dependent density functional theory and on the second-order Møller–Plesset perturbation theory method. The results illustrate the role played by electrostatic interactions in the electronic properties of anil derivatives in a liquid environment. There is a significant increase of the dipole moment in solution (20%-100%) relative to the gas-phase value. Solvent effects are mild for the absorption spectrum and linear polarizability but they can be particularly important for first hyperpolarizability. A large first hyperpolarizability contrast between the enol and keto forms is observed when absorption spectra present intense lowest-energy absorption bands. Dynamic results for the first hyperpolarizability are in qualitative agreement with the available experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Crystal-momentum-resolved contributions to high-order harmonic generation in solids

We analytically and numerically investigate the emission of high-order harmonic radiation from model solids by intense few-cycle midinfrared laser pulses. In single-active-electron approximation, we expand the active electron's wave function in a basis of adiabatic Houston states and describe the solid's electronic band structure in terms of an adjustable Kronig-Penney model potential. For high-order harmonic generation (HHG) from MgO crystals, we examine spectra from two-band and converged multiband numerical calculations. We discuss the characteristics of intra- and interband contributions to the HHG spectrum for computations including initial crystal momenta either from the Γ point at the center of the first Brioullin zone (BZ) only or from the entire first BZ. For sufficiently high intensities of the driving laser field, we find relevant contributions to HHG from the entire first BZ. Here, based on numerically calculated spectra, we scrutinize the cutoff harmonic orders as a function of the laser peak intensity and find good qualitative agreement with our analytical saddle-point-approximation predictions and published theoretical data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Building Halo Merger Trees from the Q Continuum Simulation

Cosmological N-body simulations rank among the most computationally intensive efforts today. A key challenge is the analysis of structure, substructure, and the merger history for many billions of compact particle clusters, called halos. Effectively representing the merging history of halos is essential for many galaxy formation models used to generate synthetic sky catalogs, an important application of modern cosmological simulations. Generating realistic mock catalogs requires computing the halo formation history from simulations with large volumes and billions of halos over many time steps, taking hundreds of terabytes of analysis data. We present fast parallel algorithms for producing halo merger trees and tracking halo substructure from a single-level, density-based clustering algorithm. Merger trees are created from analyzing the halo-particle membership function in adjacent snapshots, and substructure is identified by tracking the "cores" of merging halos - sets of particles near the halo center. Core tracking is performed after creating merger trees and uses the relationships found during tree construction to associate substructures with hosts. The algorithms are implemented with MPI and evaluated on a Cray XK7 supercomputer using up to 16,384 processes on data from HACC, a modern cosmological simulation framework. We present results for creating merger trees from 101 analysis snapshots taken from the Q Continuum, a large volume, high mass resolution, cosmological simulation evolving half a trillion particles.

79 ASTRONOMY AND ASTROPHYSICS↗