Integrated Radiation Transport and Nuclear Fuel Performance for Assembly-Level Simulations

The Advanced Multi-Physics (AMP) Nuclear Fuel Performance code (AMPFuel) is focused on predicting the temperature and strain within a nuclear fuel assembly to evaluate the performance and safety of existing and advanced nuclear fuel bundles within existing and advanced nuclear reactors. AMPFuel was extended to include an integrated nuclear fuel assembly capability for (one-way) coupled radiation transport and nuclear fuel assembly thermo-mechanics. This capability is the initial step toward incorporating an improved predictive nuclear fuel assembly modeling capability to accurately account for source-terms and boundary conditions of traditional (single-pin) nuclear fuel performance simulation, such as the neutron flux distribution, coolant conditions, and assembly mechanical stresses. A novel scheme is introduced for transferring the power distribution from the Scale/Denovo (Denovo) radiation transport code (structured, Cartesian mesh with smeared materials within each cell) to AMPFuel (unstructured, hexagonal mesh with a single material within each cell), allowing the use of a relatively coarse spatial mesh (10 million elements) for the radiation transport and a fine spatial mesh (3.3 billion elements) for thermo-mechanics with very little loss of accuracy. In addition, a new nuclear fuel-specific preconditioner was developed to account for the high aspect ratio of each fuel pin (12 feet axially, but 1 4 inches in diameter) with many individual fuel regions (pellets). With this novel capability, AMPFuel was used to model an entire 17 17 pressurized water reactor fuel assembly with many of the features resolved in three dimensions (for thermo-mechanics and/or neutronics), including the fuel, gap, and cladding of each of the 264 fuel pins; the 25 guide tubes; the top and bottom structural regions; and the upper and lower (neutron) reflector regions. The final, full assembly calculation was executed on Jaguar using 40,000 cores in under 10 hours to model over 162 billion degrees of freedom for 10 loading steps. The single radiation transport calculation required about 50% of the time required to solve the thermo-mechanics with a single loading step, which demonstrates that it is feasible to incorporate, in a single code, a high-fidelity radiation transport capability with a high-fidelity nuclear fuel thermo-mechanics capability and anticipate acceptable computational requirements. The results of the full assembly simulation clearly show the axial, radial, and azimuthal variation of the neutron flux, power, temperature, and deformation of the assembly, highlighting behavior that is neglected in traditional axisymmetric fuel performance codes that do not account for assembly features, such as guide tubes and control rods

Multicenter evaluation of the clinical utility of laparoscopy-assisted ERCP in patients with Roux-en-Y gastric bypass

Background and Aims The obesity epidemic has led to increased use of Roux-en-Y gastric bypass (RYGB). These patients have an increased incidence of pancreaticobiliary diseases yet standard ERCP is not possible due to surgically altered gastroduodenal anatomy. Laparoscopic-ERCP (LA-ERCP) has been proposed as an option but supporting data are derived from single center small case-series. Therefore, we conducted a large multicenter study to evaluate the feasibility, safety, and outcomes of LA-ERCP. Methods This is retrospective cohort study of adult patients with RYGB who underwent LA-ERCP in 34 centers. Data on demographics, indications, procedure success, and adverse events were collected. Procedure success was defined when all of the following were achieved: reaching the papilla, cannulating the desired duct and providing endoscopic therapy as clinically indicated. Results A total of 579 patients (median age 51, 84% women) were included. Indication for LA-ERCP was biliary in 89%, pancreatic in 8%, and both in 3%. Procedure success was achieved in 98%. Median total procedure time was 152 minutes (IQR 109-210) with median ERCP time 40 minutes (IQR 28-56). Median hospital stay was 2 days (IQR 1-3). Adverse events were 18% (laparoscopy-related 10%, ERCP-related 7%, both 1%) with the clear majority (92%) classified as mild/moderate whereas 8% were severe and 1 death occurred. Conclusion Our large multicenter study indicates that LA-ERCP in patients with RYGB is feasible with a high procedure success rate comparable with that of standard ERCP in patients with normal anatomy. ERCP-related adverse events rate is comparable with conventional ERCP, but the overall adverse event rate was higher due to the added laparoscopy-related events

BOTTOM-UP CONSTRUCTION AND 2:1 BALANCE REFINEMENT OF LINEAR OCTREES IN PARALLEL ∗

Abstract. In this article, we propose new parallel algorithms for the construction and 2:1 balance refinement of large linear octrees on distributed memory machines. Such octrees are used in many problems in computational science and engineering, e.g., object representation, image analysis, unstructured meshing, finite elements, adaptive mesh refinement and N-body simulations. Fixed-size scalability and isogranular analysis of the algorithms, using an MPI-based parallel implementation, was performed on a variety of input data and demonstrated good scalability for different processor counts (1 to 1024 processors) at the Pittsburgh Supercomputing Center’s TCS-1 AlphaServer. The results are consistent for different data distributions. Octrees with over a billion octants were constructed and balanced in less than a minute on 1024 processors. Like other existing algorithms for constructing and balancing octrees, our algorithms have O(n log n) work and O(n) storage complexity. Under reasonable assumptions on the distribution of octants and the work per octant, the parallel time complexity is O(n/np log(n/np) + np log np), were n is the final number of leaves and np is the number of processors. Key words. Linear octrees, Balance refinement, Morton encoding, large scale parallel computing, space filling curves AMS subject classifications. 65N50, 65Y05, 68W10, 68W15 1. Introduction. Spatia

Low-constant parallel algorithms for finite element simulations using linear octrees

In this article we propose parallel algorithms for the construction of conforming finite-element discretization on linear octrees. Existing octree-based discretizations scale to billions of elements, but the complexity constants can be high. In our approach we use several techniques to minimize overhead: a novel bottom-up tree-construction and 2:1 balance constraint enforcement; a Golomb-Rice encoding for compression by representing the octree and element connectivity as an Uniquely Decodable Code (UDC); overlapping communication and computation; and byte alignment for cache efficiency. The cost of applying the Laplacian is comparable to that of applying it using a direct indexing regular grid discretization with the same number of elements. Our algorithm has scaled up to four billion octants on 4096 processors on a Cray XT3 at the Pittsburgh Supercomputing Center. The overall tree construction time is under a minute in contrast to previous implementations that required several minutes; the evaluation of the discretization of a variable-coefficient Laplacian takes only a few seconds. 1