Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NUMA”

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.

53 records · Page 3

An Integrated Framework for Memory-Centric Analysis: From Trace Collection to Co-Design

The memory wall phenomenon—where advances in processor performance significantly outpace those in memory subsystems—poses a fundamental challenge for contemporary computing systems. In memory-bound applications, memory subsystem behavior dominates performance, yet existing analysis approaches present significant limitations: detailed microarchitectural simulators require days to weeks to simulate modest workloads; hardware performance counters provide only aggregate statistics that obscure temporal and spatial access patterns; and scaled simulation approaches face challenges in capturing certain behaviors that emerge at larger scales. These limitations reflect a processor-centric design philosophy increasingly misaligned with memory-bound workloads where detailed understanding of memory access patterns, cache hierarchy interactions, and contention is critical for effective optimization. This paper presents an integrated framework for memory-centric analysis that enables effective hardware-software co-design. We describe practical trace collection techniques, including hardware-assisted processor tracing with minimal overhead and portable software-based instrumentation with statistical sampling. We present multi-perspective analysis methods that examine memory behavior from temporal, sequential, spatial, and relational viewpoints, revealing distinct optimization opportunities invisible in aggregate metrics. We detail an architectural modeling framework that uses sampled traces with temporal interpolation and confidence-based filtering to evaluate cache and memory configurations. Evaluation on representative benchmarks demonstrates that this framework achieves practical accuracy (L2 cache errors of 2.64\%, confidence-filtered L3 errors of 9.92\%, bandwidth errors of 7.33\%) while providing substantial speedup (26.8×) over cycle-accurate simulation, enabling rapid design space exploration. We demonstrate how this integrated framework enables systematic identification of both hardware optimizations (memory controller tuning, bank partitioning, NUMA configuration) and software optimizations (data layout restructuring, prefetching strategies, memory-aware scheduling). Through this comprehensive treatment of the memory-centric analysis pipeline—from trace collection through architectural modeling to co-design application—we provide researchers and practitioners with practical techniques for addressing memory bottlenecks in contemporary computing systems.

Gajaria, Dhruv Mayur↗

Efficient Parallelization of a Dynamic Unstructured Application on the Tera MTA

The success of parallel computing in solving real-life computationally-intensive problems relies on their efficient mapping and execution on large-scale multiprocessor architectures. Many important applications are both unstructured and dynamic in nature, making their efficient parallel implementation a daunting task. This paper presents the parallelization of a dynamic unstructured mesh adaptation algorithm using three popular programming paradigms on three leading supercomputers. We examine an MPI message-passing implementation on the Cray T3E and the SGI Origin2OOO, a shared-memory implementation using cache coherent nonuniform memory access (CC-NUMA) of the Origin2OOO, and a multi-threaded version on the newly-released Tera Multi-threaded Architecture (MTA). We compare several critical factors of this parallel code development, including runtime, scalability, programmability, and memory overhead. Our overall results demonstrate that multi-threaded systems offer tremendous potential for quickly and efficiently solving some of the most challenging real-life problems on parallel computers.

Oliker, Leonid↗

MLP: A Parallel Programming Alternative to MPI for New Shared Memory Parallel Systems

Recent developments at the NASA AMES Research Center's NAS Division have demonstrated that the new generation of NUMA based Symmetric Multi-Processing systems (SMPs), such as the Silicon Graphics Origin 2000, can successfully execute legacy vector oriented CFD production codes at sustained rates far exceeding processing rates possible on dedicated 16 CPU Cray C90 systems. This high level of performance is achieved via shared memory based Multi-Level Parallelism (MLP). This programming approach, developed at NAS and outlined below, is distinct from the message passing paradigm of MPI. It offers parallelism at both the fine and coarse grained level, with communication latencies that are approximately 50-100 times lower than typical MPI implementations on the same platform. Such latency reductions offer the promise of performance scaling to very large CPU counts. The method draws on, but is also distinct from, the newly defined OpenMP specification, which uses compiler directives to support a limited subset of multi-level parallel operations. The NAS MLP method is general, and applicable to a large class of NASA CFD codes.

Taft, James R.↗

Parallelization of Program to Optimize Simulated Trajectories (POST3D)

This paper describes the parallelization of the Program to Optimize Simulated Trajectories (POST3D). POST3D uses a gradient-based optimization algorithm that reaches an optimum design point by moving from one design point to the next. The gradient calculations required to complete the optimization process, dominate the computational time and have been parallelized using a Single Program Multiple Data (SPMD) on a distributed memory NUMA (non-uniform memory access) architecture. The Origin2000 was used for the tests presented.

Hammond, Dana P.↗

Roofline Analysis in the Intel® Advisor to Deliver Optimized Performance for applications on Intel® Xeon Phi™ Processor

In this session we show, in two case studies, how the roofline feature of Intel Advisor has been utilized to optimize the performance of kernels of the XGC1 and PICSAR codes in preparation for Intel Knights Landing architecture. The impact of the implemented optimizations and the benefits of using the automatic roofline feature of Intel Advisor to study performance of large applications will be presented. This demonstrates an effective optimization strategy that has enabled these science applications to achieve up to 4.6 times speed-up and prepare for future exascale architectures. # Goal/Relevance of Session The roofline model [1,2] is a powerful tool for analyzing the performance of applications with respect to the theoretical peak achievable on a given computer architecture. It allows one to graphically represent the performance of an application in terms of operational intensity, i.e. the ratio of flops performed and bytes moved from memory in order to guide optimization efforts. Given the scale and complexity of modern science applications, it can often be a tedious task for the user to perform the analysis on the level of functions or loops to identify where performance gains can be made. With new Intel tools, it is now possible to automate this task, as well as base the estimates of peak performance on measurements rather than vendor specifications. The goal of this session is to demonstrate how the roofline feature of Intel Advisor can be used to balance memory vs. computation related optimization efforts and effectively identify performance bottlenecks. A series of typical optimization techniques: cache blocking, structure refactoring, data alignment, and vectorization illustrated by the kernel cases will be addressed. # Description of the codes ## XGC1 The XGC1 code [3] is a magnetic fusion Particle-In-Cell code that uses an unstructured mesh for its Poisson solver that allows it to accurately resolve the edge plasma of a magnetic fusion device. After recent optimizations to its collision kernel [4], most of the computing time is spent in the electron push (pushe) kernel, where these optimization efforts have been focused. The kernel code scaled well with MPI+OpenMP but had almost no automatic compiler vectorization, in part due to indirect memory addresses and in part due to low trip counts of low-level loops that would be candidates for vectorization. Particle blocking and sorting have been implemented to increase trip counts of low-level loops and improve memory locality, and OpenMP directives have been added to vectorize compute-intensive loops that were identified by Advisor. The optimizations have improved the performance of the pushe kernel 2x on Haswell processors and 1.7x on KNL. The KNL node-for-node performance has been brought to within 30% of a NERSC Cori phase I Haswell node and we expect to bridge this gap by reducing the memory footprint of compute intensive routines to improve cache reuse. ## PICSAR is a Fortran/Python high-performance Particle-In-Cell library targeting at MIC architectures first designed to be coupled with the PIC code WARP for the simulation of laser-matter interaction and particle accelerators. PICSAR also contains a FORTRAN stand-alone kernel for performance studies and benchmarks. A MPI domain decomposition is used between NUMA domains and a tile decomposition (cache-blocking) handled by OpenMP has been added for shared-memory parallelism and better cache management. The so-called current deposition and field gathering steps that compose the PIC time loop constitute major hotspots that have been rewritten to enable more efficient vectorization. Particle communications between tiles and MPI domain has been merged and parallelized. All considered, these improvements provide speedups of 3.1 for order 1 and 4.6 for order 3 interpolation shape factors on KNL configured in SNC4 quadrant flat mode. Performance is similar between a node of cori phase 1 and KNL at order 1 and better on KNL by a factor 1.6 at order 3 with the considered test case (homogeneous thermal plasma).

97 MATHEMATICS AND COMPUTING↗

A Review of Finite Element Method Models for Nuclear Graphite Applications

Graphite components serve as structural support and as a neutron moderator in multiple nuclear power stations and is also a prospective material in Generation IV reactor concepts. The environment (e.g. high temperature and intense radiation) within a graphite core changes the geometry, microstructure and material properties of the graphite components, that could limit the lifetime (i.e. safe operation) of a reactor. Some of the assessments made for graphite reactor core components include Finite Element Method (FEM) simulations. These simulations are used to estimate the stress state and strain state and probability of failure of irradiated graphite components. FEM also forms part of the validation of measurements or as a complementary technique to experimental procedures in graphite research. A comprehensive overview of FEM methodologies is provided in this review, including: FEM software, basics steps to set up a FEM model, and analytical solutions used to validate simulations for this field of study. In conclusion, the objective of this paper is to inform researchers in the different available FEM techniques and tools used in the research of nuclear graphite.

42 ENGINEERING↗

Effect of Whole Pelvic Radiotherapy for Patients With Locally Advanced Prostate Cancer Treated With Radiotherapy and Long-Term Androgen Deprivation Therapy

Purpose: To evaluate the effect of whole pelvic radiotherapy (WPRT) in prostate cancer patients treated with RT and long-term (>1 year) androgen deprivation therapy (ADT). Methods and materials: Prostate cancer patients with high-risk features (Stage T3-T4 and/or Gleason score {>=}7 and/or prostate-specific antigen level {>=}20 ng/mL) who had undergone RT and long-term ADT were included in the present analysis. Patients with bowel inflammatory disease, colon diverticula, and colon diverticulitis were excluded from WPRT and treated with prostate-only radiotherapy (PORT). Patients were grouped according to nodal risk involvement as assessed by the Roach formula using different cutoff levels (15%, 20%, 25%, and 30%). Biochemical disease-free survival (bDFS) was analyzed in each group according to the RT type (WPRT or PORT). Results: A total of 358 patients treated between 1994 and 2007 were included in the analysis (46.9% with WPRT and 53.1% with PORT). The median duration of ADT was 24 months (range, 12-38). With a median follow-up of 52 months (range, 20-150), the overall 4-year bDFS rate was 90.5%. The 4-year bDFS rate was similar between the patients who had undergone WPRT or PORT (90.4% vs. 90.5%; p = NS). However, in the group of patients with the greatest nodal risk (>30%), a significant bDFS improvement was recorded for the patients who had undergone WPRT (p = .03). No differences were seen in acute toxicity among the patients treated with WPRT or PORT. The late gastrointestinal toxicity was similar in patients treated with PORT or WPRT (p = NS). Conclusions: Our analysis has supported the use of WPRT in association with long-term ADT for patients with high-risk nodal involvement (>30%), although a definitive recommendation should be confirmed by a randomized trial.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Early Proctoscopy is a Surrogate Endpoint of Late Rectal Toxicity in Prostate Cancer Treated With Radiotherapy

Purpose: To predict the grade and incidence of late clinical rectal toxicity through short-term (1 year) mucosal alterations. Methods and Materials: Patients with prostate adenocarcinoma treated with curative or adjuvant radiotherapy underwent proctoscopy a year after the course of radiotherapy. Mucosal changes were classified by the Vienna Rectoscopy Score (VRS). Late toxicity data were analyzed according to the Kaplan-Meier method. Comparison between prognosis groups was performed by log-rank analysis. Results: After a median follow-up time of 45 months (range, 18-99), the 3-year incidence of grade {>=}2 rectal late toxicity according to the criteria of the European Organization for Research and Treatment of Cancer and the Radiation Therapy Oncology Group was 24%, with all patients (24/24; 100%) experiencing rectal bleeding. The occurrence of grade {>=}2 clinical rectal late toxicity was higher in patients with grade {>=}2 (32% vs. 15 %, p = 0.02) or grade {>=}3 VRS telangiectasia (47% vs. 17%, p {<=} 0.01) and an overall VRS score of {>=}2 (31% vs. 16 %, p = 0.04) or {>=}3 (48% vs. 17%, p = 0.01) at the 1-year proctoscopy. Conclusions: Early proctoscopy (1 year) predicts late rectal bleeding and therefore can be used as a surrogate endpoint for late rectal toxicity in studies aimed at reducing this frequent complication.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Electron tomography of unirradiated and irradiated nuclear graphite

Graphite is the moderator material of several Generation IV nuclear reactor concepts, as well as the British Advanced Gas-cooled Reactors (AGR). Porosity can heavily influence the material properties, mechanical irradiation response, and neutron induced shrinkage or swelling of nuclear-grade graphite. Due to the sub-micron size of several types of pores found in graphite, only a high-resolution imaging technique such as electron tomography are capable of visualizing these features in three dimensions. In this research, we used electron tomography to characterize as-received and neutron irradiated samples of IG-110 nuclear-grade graphite to show for the first time the 3D structure of both native and irradiation-induced nano-cracks. This technique also reveals unique characteristics of graphite such as the structure that surrounds pores and could be used to inform molecular dynamic simulations of irradiated graphite and experimental techniques such as gas-absorption. This work also shows the utility of this technique for the study of other nuclear porous carbon-based materials.

36 MATERIALS SCIENCE↗

Postoperative Intensity Modulated Radiation Therapy in High Risk Prostate Cancer: A Dosimetric Comparison

The aim of this study was to compare intensity-modulated radiation therapy (IMRT) with 3D conformal technique (3D-CRT), with respect to target coverage and irradiation of organs at risk for high dose postoperative radiotherapy (PORT) of the prostate fossa. 3D-CRT and IMRT treatment plans were compared with respect to dose to the rectum and bladder. The dosimetric comparison was carried out in 15 patients considering 2 different scenarios: (1) exclusive prostate fossa irradiation, and (2) pelvic node irradiation followed by a boost on the prostate fossa. In scenario (1), a 3D-CRT plan (box technique) and an IMRT plan were calculated and compared for each patient. In scenario (2), 3 treatment plans were calculated and compared for each patient: (a) 3D-CRT box technique for both pelvic (prophylactic nodal irradiation) and prostate fossa irradiation (3D-CRT only); (b) 3D-CRT box technique for pelvic irradiation followed by an IMRT boost to the prostatic fossa (hybrid 3D-CRT and IMRT); and (c) IMRT for both pelvic and prostate fossa irradiation (IMRT only). For exclusive prostate fossa irradiation, IMRT significantly reduced the dose to the rectum (lower Dmean, V50%, V75%, V90%, V100%, EUD, and NTCP) and the bladder (lower Dmean, V50%, V90%, EUD and NTCP). When prophylactic irradiation of the pelvis was also considered, plan C (IMRT only) performed better than plan B (hybrid 3D-CRT and IMRT) as respect to both rectum and bladder irradiation (reduction of Dmean, V50%, V75%, V90%, equivalent uniform dose [EUD], and normal tissue complication probability [NTCP]). Plan (b) (hybrid 3D-CRT and IMRT) performed better than plan (a) (3D-CRT only) with respect to dose to the rectum (lower Dmean, V75%, V90%, V100%, EUD, and NTCP) and the bladder (Dmean, EUD, and NTCP). Postoperative IMRT in prostate cancer significantly reduces rectum and bladder irradiation compared with 3D-CRT.

61 RADIATION PROTECTION AND DOSIMETRY↗

Identification of Family Non-universal Gauge Bosons in High-energy Electron-positron Collisions

We examine effects on measurable observables in e{sup +}e{sup -} collisions resulting from the existence of additional neutral gauge bosons originating in extensions of the standard model. In particular, we consider family non-universal neutral gauge bosons occurring in the minimal supersymmetric standard model and in the Sp(6){sub L} x U(1){sub Y} model, as well as other theoretically motivated popular neutral gauge bosons. We show how the proper employment of the generation-dependent couplings of the extra gauge boson, and the appropriate adjustment of the beam polarization, not only improved the identification of the models but also enhanced the discovery potential of the family non-universal extra gauge bosons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Acceleration of the IMplicit–EXplicit nonhydrostatic unified model of the atmosphere on manycore processors

We present the acceleration of an IMplicit-EXplicit (IMEX) nonhydrostatic atmospheric model on manycore processors such as graphic processing units (GPUs) and Intel's Many Integrated Core (MIC) architecture. IMEX time integration methods sidestep the constraint imposed by the Courant-Friedrichs-Lewy condition on explicit methods through corrective implicit solves within each time step. In this work, we implement and evaluate the performance of IMEX on manycore processors relative to explicit methods. Using 3D-IMEX at Courant number C = 15, we obtained a speedup of about 4x relative to an explicit time stepping method run with the maximum allowable C = 1. Moreover, the unconditional stability of IMEX with respect to the fast waves means the speedup can increase significantly with the Courant number as long as the accuracy of the resulting solution is acceptable. We show a speedup of 100x at C = 150 using 1D-IMEX to demonstrate this point. Several improvements on the IMEX procedure were necessary in order to outperform our results with explicit methods: (a) reducing the number of degrees of freedom of the IMEX formulation by forming the Schur complement, (b) formulating a horizontally explicit vertically implicit 1D-IMEX scheme that has a lower workload and better scalability than 3D-IMEX, (c) using high-order polynomial preconditioners to reduce the condition number of the resulting system, and (d) using a direct solver for the 1D-IMEX method by performing and storing LU factorizations once to obtain a constant cost for any Courant number. Without all of these improvements, explicit time integration methods turned out to be difficult to beat. We discuss in detail the IMEX infrastructure required for formulating and implementing efficient methods on manycore processors. Several parametric studies are conducted to demonstrate the gain from each of the abovementioned improvements. Lastly, we validate our results with standard benchmark problems in numerical weather prediction and evaluate the performance and scalability of the IMEX method using up to 4192 GPUs and 16 Knights Landing processors.

97 MATHEMATICS AND COMPUTING↗

Determination of the specific area of liquid gas and the velocity of weak pressure waves in aqueous foams

The specific area of an aqueous foam contained in a tube is obtained by determining the fraction of the quantity of light emitted by a source and diffused by the specific area of the column of the foam. The velocity of the waves of weak pressure (which propagate in the air with the velocity of sound) is measured by noticing the moment when the wave penetrates the column and the moment when, having reached the opposite side, it determines a variation of the light diffused by the area of the latter.

Saint-Cloud, Jean↗

Characterization of the Irradiation Effects in Nuclear Graphite

Graphite serves as the neutron moderator of British Gas-Cooled Reactors and is a candidate material for the Generation IV of nuclear power stations. The neutron irradiation affects the crystalline structure of graphite; this damage results in volumetric changes of the bulk material and generation of porosity. Dimensional changes induced by irradiation in graphite bricks can make difficult the refueling operations and reduce the stability of the graphite core. Some of the neutron irradiation damage mechanisms in the crystal lattice of graphite are not well understood. Results of TEM, SEM, FIB-SEM tomography, gas adsorption and optical microscopy will be combined to help understand how the irradiation affects graphite at multiple length scales. Furthermore, some new evidence of porosity and microstructural changes generated by irradiation will be presented.

Arregui Mena, Jose'↗

Microstructural characterization of the oxidation of nuclear graphite under chronic and accident conditions via XCT and synchrotron XCT

Graphite is a candidate material to moderate fast neutrons and for structural components in the US next-generation graphite-moderated reactors. A graphite core is conceived as a large formation of interconnected bricks that primarily serves as the moderator of fast neutrons, holds instrumentation, fuel elements, control rods, and is a receptacle for molten salts for Molten Salt Reactors (MSRs) designs. During the operation of a graphite-moderated power plant, graphite components might be subjected to chronic oxidation during normal operating conditions or aggressive oxidation as a result of accidental ingress of air that reacts vigorously with the graphite core. Using synchrotron, x-ray computed tomography (XCT), this research systematically characterized microstructural changes that accompanied these two oxidation scenarios. Chronic oxidation was studied by characterizing IG-110, PCEA, and NBG-18 specimens that were gradually oxidized in air at a low temperature (520°C). The accidental ingress of air into the graphite reactor was simulated by subjecting the grades mentioned above to high-temperature oxidation at approximately 750°C. This research is the first, in situ, systematic characterization of nuclear graphite microstructural evolution that can be associated with the two possible oxidation scenarios and provide insight into related repercussions. The results show that the microstructure and pore connectivity influence the rate of oxidation and evolution of the microstructure under the two oxidation regimes. These results are crucial to understanding which grades are more resilient to each type of oxidation and documenting the damage created in the graphite.

Arregui Mena, Jose'↗