Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Large deviations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,477 records · Page 82

Single-site Lennard-Jones models via polynomial chaos surrogates of Monte Carlo molecular simulation

In this work, two Polynomial Chaos (PC) surrogates were generated to reproduce Monte Carlo (MC) molecular simulation results of the canonical (single-phase) and the NVT-Gibbs (two-phase) ensembles for a system of normalized structureless Lennard-Jones (LJ) particles. The main advantage of such surrogates, once generated, is the capability of accurately computing the needed thermodynamic quantities in a few seconds, thus efficiently replacing the computationally expensive MC molecular simulations. Benefiting from the tremendous computational time reduction, the PC surrogates were used to conduct large-scale optimization in order to propose single-site LJ models for several simple molecules. Experimental data, a set of supercritical isotherms, and part of the two-phase envelope, of several pure components were used for tuning the LJ parameters (ε, σ). Based on the conducted optimization, excellent fit was obtained for different noble gases (Ar, Kr, and Xe) and other small molecules (CH{sub 4}, N{sub 2}, and CO). On the other hand, due to the simplicity of the LJ model used, dramatic deviations between simulation and experimental data were observed, especially in the two-phase region, for more complex molecules such as CO{sub 2} and C{sub 2} H{sub 6}.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Clustering of local group distances: publication bias or correlated measurements? I. The large Magellanic cloud

The distance to the Large Magellanic Cloud (LMC) represents a key local rung of the extragalactic distance ladder yet the galaxy's distance modulus has long been an issue of contention, in particular in view of claims that most newly determined distance moduli cluster tightly—and with a small spread—around the 'canonical' distance modulus, (m – M){sub 0} = 18.50 mag. We compiled 233 separate LMC distance determinations published between 1990 and 2013. Our analysis of the individual distance moduli, as well as of their two-year means and standard deviations resulting from this largest data set of LMC distance moduli available to date, focuses specifically on Cepheid and RR Lyrae variable-star tracer populations, as well as on distance estimates based on features in the observational Hertzsprung-Russell diagram. We conclude that strong publication bias is unlikely to have been the main driver of the majority of published LMC distance moduli. However, for a given distance tracer, the body of publications leading to the tightly clustered distances is based on highly non-independent tracer samples and analysis methods, hence leading to significant correlations among the LMC distances reported in subsequent articles. Based on a careful, weighted combination, in a statistical sense, of the main stellar population tracers, we recommend that a slightly adjusted canonical distance modulus of (m – M){sub 0} = 18.49 ± 0.09 mag be used for all practical purposes that require a general distance scale without the need for accuracies of better than a few percent.

79 ASTRONOMY AND ASTROPHYSICS↗

TU-AB-202-06: Quantitative Evaluation of Deformable Image Registration in MRI-Guided Adaptive Radiation Therapy

Purpose: To assess the performance of the deformable image registration algorithm used for MRI-guided adaptive radiation therapy using image feature analysis. Methods: MR images were collected from five patients treated on the MRIdian (ViewRay, Inc., Oakwood Village, OH), a three head Cobalt-60 therapy machine with an 0.35 T MR system. The images were acquired immediately prior to treatment with a uniform 1.5 mm resolution. Treatment sites were as follows: head/neck, lung, breast, stomach, and bladder. Deformable image registration was performed using the ViewRay software between the first fraction MRI and the final fraction MRI, and the DICE similarity coefficient (DSC) for the skin contours was reported. The SIFT and Harris feature detection and matching algorithms identified point features in each image separately, then found matching features in the other image. The target registration error (TRE) was defined as the vector distance between matched features on the two image sets. Each deformation was evaluated based on comparison of average TRE and DSC. Results: Image feature analysis produced between 2000–9500 points for evaluation on the patient images. The average (± standard deviation) TRE for all patients was 3.3 mm (±3.1 mm), and the passing rate of TRE<3 mm was 60% on the images. The head/neck patient had the best average TRE (1.9 mm±2.3 mm) and the best passing rate (80%). The lung patient had the worst average TRE (4.8 mm±3.3 mm) and the worst passing rate (37.2%). DSC was not significantly correlated with either TRE (p=0.63) or passing rate (p=0.55). Conclusions: Feature matching provides a quantitative assessment of deformable image registration, with a large number of data points for analysis. The TRE of matched features can be used to evaluate the registration of many objects throughout the volume, whereas DSC mainly provides a measure of gross overlap. We have a research agreement with ViewRay Inc.

60 APPLIED LIFE SCIENCES↗

Search for heavy particles decaying into top-quark pairs using lepton-plus-jets events in proton–proton collisions at $\sqrt{s} = 13$ $\text {TeV}$ with the ATLAS detector

Here, a search for new heavy particles that decay into top-quark pairs is performed using data collected from proton–proton collisions at a centre-of-mass energy of 13 TeV by the ATLAS detector at the Large Hadron Collider. The integrated luminosity of the data sample is 36.1 fb –1 . Events consistent with top-quark pair production are selected by requiring a single isolated charged lepton, missing transverse momentum and jet activity compatible with a hadronic top-quark decay. Jets identified as likely to contain b-hadrons are required to reduce the background from other Standard Model processes. The invariant mass spectrum of the candidate top-quark pairs is examined for local excesses above the background expectation. No significant deviations from the Standard Model predictions are found. Exclusion limits are set on the production cross-section times branching ratio for hypothetical Z' bosons, Kaluza–Kein gluons and Kaluza–Klein gravitons that decay into top-quark pairs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for heavy charged long-lived particles in the ATLAS detector in 31.6 fb &#x2212; 1 of proton-proton collision data at s = 13 TeV

CERN-LHC. Search for heavy charged long-lived particles in the ATLAS detector in 36.1/fb of proton-proton collision data at sqrt(s) = 13 TeV. A search for heavy charged long-lived particles is performed using a data sample of 36.1/fb of proton-proton collisions at sqrt(s) = 13 TeV collected by the ATLAS experiment at the Large Hadron Collider. The search is based on observables related to ionization energy loss and time of flight, which are sensitive to the velocity of heavy charged particles traveling significantly slower than the speed of light. Multiple search strategies for a wide range of lifetimes, corresponding to path lengths of a few meters, are defined as model-independently as possible, by referencing several representative physics cases that yield long-lived particles within supersymmetric models, such as gluinos/squarks (R-hadrons), charginos and staus. No significant deviations from the expected Standard Model background are observed. Upper limits at 95% confidence level are provided on the production cross sections of long-lived R-hadrons as well as directly pair-produced staus and charginos. These results translate into lower limits on the masses of long-lived gluino, sbottom and stop R-hadrons, as well as staus and charginos of 2000 GeV, 1250 GeV, 1340 GeV, 430 GeV and 1090 GeV, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Measurement and interpretation of same-sign W boson pair production in association with two jets in $pp$ collisions at $\sqrt{s}$ = 13 TeV with the ATLAS detector

This paper presents the measurement of fiducial and differential cross sections for both the inclusive and electroweak production of a same-sign W-boson pair in association with two jets (W ± W ± jj) using 139 fb –1 of proton-proton collision data recorded at a centre-of-mass energy of $\sqrt{s}$ = 13 TeV by the ATLAS detector at the Large Hadron Collider. The analysis is performed by selecting two same-charge leptons, electron or muon, and at least two jets with large invariant mass and a large rapidity difference. The measured fiducial cross sections for electroweak and inclusive W ± W ± jj production are 2.92 ± 0.22 (stat.) ± 0.19 (syst.) fb and 3.38 ± 0.22 (stat.) ± 0.19 (syst.) fb, respectively, in agreement with Standard Model predictions. The measurements are used to constrain anomalous quartic gauge couplings by extracting 95% confidence level intervals on dimension-8 operators. A search for doubly charged Higgs bosons H ±± that are produced in vector-boson fusion processes and decay into a same-sign W boson pair is performed. The largest deviation from the Standard Model occurs for an H ±± mass near 450 GeV, with a global significance of 2.5 standard deviations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fast-ion transport in quasisymmetric equilibria in the presence of a resonant Alfvénic perturbation

Significant progress has been made in designing magnetic fields that provide excellent confinement of the guiding-centre trajectories of alpha particles using quasisymmetry (QS). Given the reduction in this transport channel, we assess the impact of resonant Alfvén eigenmodes (AEs) on the guiding-centre motion. The AE amplitudes are chosen to be consistent with experimental measurements and large-scale simulations. We evaluate the drift resonance condition, phase-space island width and island overlap criterion for quasisymmetric configurations. Kinetic Poincaré plots elucidate features of the transport, including stiff transport above a critical perturbation amplitude. Our analysis highlights key departures from the AE-driven transport in tokamaks, such as the avoidance of phase-space island overlap in quasihelical configurations and the enhanced transport due to wide phase-space islands in low magnetic shear configurations. In configurations that are closer to QS, with QS deviations $\delta B/B_0 \lesssim 10^{-3}$ , the transport is primarily driven by the AE, while configurations that are further from QS, $\delta B/B_0 \sim 10^{-2}$ , experience significant transport due to the QS-breaking fields in addition to the AE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Anti-site disorder and improved functionality of Mn₂NiX (X = Al, Ga, In, Sn) inverse Heusler alloys

Recent first-principles calculations have predicted Mn₂NiX (X = Al, Ga, In, Sn) alloys to be magnetic shape memory alloys. Moreover, experiments on Mn₂NiGa and Mn₂NiSn suggest that the alloys deviate from the perfect inverse Heusler arrangement and that there is chemical disorder at the sublattices with tetrahedral symmetry. In this work, we investigate the effects of such chemical disorder on phase stabilities and magnetic properties using first-principles electronic structure methods. We find that except Mn₂NiAl, all other alloys show signatures of martensitic transformations in presence of anti-site disorder at the sublattices with tetrahedral symmetry. This improves the possibilities of realizing martensitic transformations at relatively low fields and the possibilities of obtaining significantly large inverse magneto-caloric effects, in comparison to perfect inverse Heusler arrangement of atoms. We analyze the origin of such improvements in functional properties by investigating electronic structures and magnetic exchange interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurement of simplified template cross sections of the Higgs boson produced in association with W or Z bosons in the H → b b ¯ decay channel in proton-proton collisions at s = 13 TeV

Differential cross sections are measured for the standard model Higgs boson produced in association with vector bosons (W, Z) and decaying to a pair of b quarks. Measurements are performed within the framework of the simplified template cross sections. The analysis relies on the leptonic decays of the W and Z bosons, resulting in final states with 0, 1, or 2 electrons or muons. The Higgs boson candidates are either reconstructed from pairs of resolved b-tagged jets, or from single large-radius jets containing the particles arising from two b quarks. Proton-proton collision data at $\sqrt{s}$ =13 TeV, collected by the CMS experiment in 2016–2018 and corresponding to a total integrated luminosity of 138 fb -1 , are analyzed. The inclusive signal strength, defined as the product of the observed production cross section and branching fraction relative to the standard model expectation, combining all analysis categories, is found to be μ = $1.15^{+0.22}_{-0.20}$. This corresponds to an observed (expected) significance of 6.3 (5.6) standard deviations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Evaluation of the dosimetric impact of interfractional anatomical variations on prostate proton therapy using daily in-room CT images

Purpose: To quantify interfractional anatomical variations and their dosimetric impact during the course of fractionated proton therapy (PT) of prostate cancer and to assess the robustness of the current treatment planning techniques. Methods: Simulation and daily in-room CT scans from ten prostate carcinoma patients were analyzed. PT treatment plans (78 Gy in 39 fractions of 2 Gy) were created on the simulation CT, delivering 25 fractions to PTV1 (expanded from prostate and seminal vesicles), followed by 14 boost fractions to PTV2 (expanded from prostate). Plans were subsequently applied to daily CT, with beams aligned to the prostate center in the sagittal plane. For five patients having a sufficiently large daily imaging volume, structure contours were manually drawn, and plans were evaluated for all CT sets. For the other five patients, the plans were evaluated for six selected fractions. The daily CT was matched to the simulation CT through deformable registration. The registration accuracy was validated for each fraction, and the three patients with a large number of accurately registered fractions were used for dose accumulation. Results: In individual fractions, the coverage of the prostate, seminal vesicles, and PTV1 was generally maintained at the corresponding prescription dose. For PTV2, the volume covered by the fractional prescription dose of 2 Gy (i.e., V2) was, on average, reduced by less than 3% compared to the simulation plan. Among the 225 (39 x 5 + 6 x 5) fractions examined, 15 showed a V2 reduction larger than 5%, of which ten were caused by a large variation in rectal gas, and five were due to a prostate shift in the craniocaudal direction. The fractional dose to the anterior rectal wall was found to increase for one patient who had large rectal gas volume in 25 of the 39 fractions, and another who experienced significant prostate volume reduction during the treatment. The fractional bladder dose generally increased with decreasing fullness. In the total accumulated dose for the three patients after excluding a few fractions with inaccurate registration due to a large amount of rectal gas (a condition inconsistent with RTOG protocol), 98.5%, 96.6%, and 98.2% of the PTV2 received the prescription dose of 78 Gy. The V75 and V70 of the anterior rectal wall and bladder both remained within tolerance. Conclusions: The results confirm that the PT planning techniques and dose constraints used at our institution ensure that target coverage to the prescription dose is maintained in the presence of interfractional anatomical variations. Dose coverage in individual fractions can be compromised, and normal tissue dose increased, due to deviations in the bladder and rectal volume compared to the simulation plans or progressive changes in the prostate volume during the treatment. Deviations from the plan can be reduced with efforts aimed at maintaining consistent daily patient anatomy.

61 RADIATION PROTECTION AND DOSIMETRY↗

WE-D-BRA-01: FEATURED PRESENTATION and BEST IN PHYSICS (THERAPY): Predicting Potentially Problematic VMAT Treatment Plans Before Patient Specific QA Measurements

Purpose: Several promising IMRT QA tools have been developed in recent years to combat problems found in the lack of sensitivity in planar dose measurements analyzed using consensus gamma analysis criteria. The increased complexity and added information with such devices adds not only increased time, but new challenges in determining endpoints for pass/fail criteria. Using a large cohort of previously measured planar IMRT QA data, it may be possible to correlate potentially problematic plans with calculated plan metrics that can be done a priori, such that these tools can be used only in clinically relevant situations. Methods: 90 previously measured, clinically delivered VMAT plans were exported in DICOM RT format. Using a Matlab program, plan metrics were computed based on a previously developed set of equations (Du et al. 2014). These metrics included MU-weighted beam irregularity, which quantifies an MLC shape’s deviation from that of a circle. Machine delivery parameters such as MU delivered per degree and leaf movement in millimeters per MU were also calculated. Based on a previous analysis of 394 IMRT QA measurements, a “failing” plan was defined as one with less than 85% gamma pass rate when computed at 1% dose difference and 2 mm distance to agreement; 16 of the 90 plans were identified as failing Results: A ROC curve was generated with an AUC of 0.9409. All 16 outliers were detected with a specificity of 85.1% when using a threshold value based on a linear combination of the MU-weighted plan irregularity and the leaf speed in mm/degree (average MU delivered per degree multiplied by average leaf movement in mm per MU). Conclusion: This data supports that potentially problematic VMAT plans can possibly be predicted before measurement by assessing the MU-weighted irregularity of the MLC shapes combined with the averaged leaf speed of the MLC.

60 APPLIED LIFE SCIENCES↗

SU-G-206-07: Dual-Energy CT Inter- and Intra-Scanner Variability Within One Make and Model

Purpose: It can be logistically quite difficult to scan patients on the same exact device for their repeat visits in multi-scanner facilities. The reliability between dual-energy CT scanners’ quantitative results is not known, nor is their individual repeatability. Therefore, we evaluated inter- and intra-scanner variability with respect to several key clinical quantitative metrics specific to dual-energy CT. Methods: Eleven identical GE HD-750 CT scanners in a busy clinical environment were used to perform dual-energy (DE) CT scans of a large elliptical quality control (QC) phantom (Gammex, Inc.; Middleton, WI) which contains many standard insert materials. The DE-QC phantom was scanned bi-weekly during 2016; 3 to 4 scans were obtained from each scanner (a total of 35 data sets were used for analysis). Iodine accuracy for the 2mg/ml, 5mg/ml and 15mg/ml rods (from the Iodine(Water) image set) and soft tissue HU (40 HU based on NIST constants) from the 50keV data set were used to assess inter- and intra-scanner variability (standard deviation). Results: Intra-scanner variability average for 2mg/ml Iodine was 0.10 mg/ml (range 0.05–0.15 mg/ml), for 5mg/ml Iodine was 0.12 mg/ml (range 0.07–0.16 mg/ml), for 15 mg/ml Iodine was 0.25 mg/ml (range 0.16–0.37 mg/ml), and for the soft tissue inserts was 2.1 HU (range 1.8–2.6 HU). Inter-scanner variability average for 2mg/ml Iodine was 0.16 mg/ml (range 0.11–0.19 mg/ml), for 5mg/ml Iodine was 0.18 mg/ml (range 0.11–0.22 mg/ml), for 15 mg/ml Iodine was 0.35 mg/ml (range 0.23–0.44 mg/ml), and for the soft tissue inserts was 3.8 HU (range 3.1–4.5 HU). Conclusion: Intra-scanner variability for the iodine and soft tissue inserts averaged 3.1% and 5.2% respectively, and inter-scanner variability for these regions analyzed averaged 5.0% and 9.5%, respectively. Future work will include determination of smallest measurable change and acceptable limits for DE-CT scanner variability over longer time intervals. This research has been supported by funds from Dr. William Murphy, Jr., the John S. Dunn, Sr. Distinguished Chair in Diagnostic Imaging at MD Anderson Cancer Center.

60 APPLIED LIFE SCIENCES↗

Search for doubly and singly charged Higgs bosons decaying into vector bosons in multi-lepton final states with the ATLAS detector using proton-proton collisions at $\sqrt{s}$ = 13 TeV

A search for charged Higgs bosons decaying into W ± W ± or W ± Z bosons is performed, involving experimental signatures with two leptons of the same charge, or three or four leptons with a variety of charge combinations, missing transverse momentum and jets. A data sample of proton-proton collisions at a centre-of-mass energy of 13 TeV recorded with the ATLAS detector at the Large Hadron Collider between 2015 and 2018 is used. The data correspond to a total integrated luminosity of 139 fb -1 . The search is guided by a type-II seesaw model that extends the scalar sector of the Standard Model with a scalar triplet, leading to a phenomenology that includes doubly and singly charged Higgs bosons. Two scenarios are explored, corresponding to the pair production of doubly charged H ±± bosons, or the associated production of a doubly charged H ±± boson and a singly charged H ± boson. No significant deviations from the Standard Model predictions are observed. H ±± bosons are excluded at 95% confidence level up to 350 GeV and 230 GeV for the pair and associated production modes, respectively.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A versatile machine learning workflow for high-throughput analysis of supported metal catalyst particles

Accurate and efficient characterization of nanoparticles (NPs), particularly regarding particle size distribution, is essential for advancing our understanding of their structure-property relationship and facilitating their design for various applications. In this study, we introduce a novel two-stage artificial intelligence (AI)-driven workflow for NP analysis that leverages prompt engineering techniques from state-of-the-art single-stage object detection and large-scale vision transformer (ViT) architectures. This methodology is applied to transmission electron microscopy (TEM) and scanning TEM (STEM) images of heterogeneous catalysts, enabling high-resolution, high-throughput analysis of particle size distributions for supported metal catalyst NPs. The model's performance in detecting and segmenting NPs is validated across diverse heterogeneous catalyst systems, including various metals (Ru, Cu, PtCo, and Pt), supports (silica (SiO 2 ), γ-alumina (γ-Al 2 O 3 ), and carbon black), and particle diameter size distributions with mean and standard deviations ranging from 1.6 ± 0.2 nm to 9.7 ± 4.6 nm. The proposed machine learning (ML) methodology achieved an average F1 overlap score of 0.91 ± 0.01 and demonstrated the ability to disentangle overlapping NPs anchored on catalytic support materials. The segmentation accuracy is further validated using the Hausdorff distance and robust Hausdorff distance metrics, with the 90th percent of the robust Hausdorff distance showing errors within 0.4 ± 0.1 nm to 1.4 ± 0.6 nm. In conclusion, our AI-assisted NP analysis workflow demonstrates robust generalization across diverse datasets and can be readily applied to similar NP segmentation tasks without requiring costly model retraining.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fine-Scale Variability of Observed and Simulated Surface Albedo Over the Southern Great Plains

Surface albedo can be highly variable in both space and time. The Department of Energy’s Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field study provides a unique opportunity to characterize the variability over the Southern Great Plains of the United States using integrated tower, aircraft, and satellite observations of the surface albedo and Normalized Difference Vegetation Index (NDVI). The primary advantage of the aircraft and satellite observations is the ability to examine the spatial distribution of surface albedo and NDVI over a large area, while the tower measurements have the ability to examine both diurnal and day-to-day variability at a single location. In general, consistency was found between both the broadband (BB) albedo and NDVI measured from the surface, air and space. There was a small increase from 0.186 to 0.194 in the aircraft BB surface albedo between May and September (about 4% change), while the MODIS black-sky BB surface albedo increased from 0.151 to 0.166 over the same period (about 10% change) while the standard deviations in the aircraft and MODIS BB albedo were similar. The largest seasonal differences in both the aircraft BB albedo and NDVI were found for areas with winter wheat or forest, while areas with pasture or grasses showed a smaller seasonal diversity. The Weather Research and Forecasting (WRF) model was used to simulate the BB surface albedo. In comparison with the aircraft and satellite observations, the WRF-simulated BB surface albedo had no seasonal change and a much narrower distribution.

Berg, Larry K.↗

Tip-based source of femtosecond electron pulses at 30 keV

We present a nano-scale photoelectron source, optimized for ultrashort pulse durations and well-suited for time-resolved diffraction and advanced laser acceleration experiments. A tungsten tip of several-ten-nanometers diameter mounted in a suppressor-extractor electrode configuration allows the generation of 30 keV electron pulses with an estimated pulse duration of 9 fs (standard deviation; 21 fs full width at half maximum) at the gun exit. We infer the pulse duration from particle tracking simulations, which are in excellent agreement with experimental measurements of the electron-optical properties of the source in the spatial domain. We also demonstrate femtosecond-laser triggered operation of the apparatus. The temporal broadening of the pulse upon propagation to a diffraction sample can be greatly reduced by collimating the beam. Besides the short electron pulse duration, a tip-based source is expected to feature a large transverse coherence and a nanometric emittance.

07 ISOTOPE AND RADIATION SOURCES↗

Probing asymmetric structures in the outskirts of galaxies

Upcoming large imaging surveys will allow detailed studies of the structure and morphology of galaxies aimed at addressing how galaxies form and evolve. Computational approaches are needed to characterize their morphologies over large samples. We introduce an automatic method to quantify the outer structure of galaxies. The key to our approach is the division of a galaxy image into two sections delineated by the isophote, which encloses half the total brightness of the galaxy. We call the central section the inner half-flux region (IHR) and the outer section the outer half-flux region (OHR). From this division, we derive two parameters: A {sub o}, which measures the asymmetry of the OHR, and D {sub o}, which measures the deviation of the intensity weighted centroid of the OHR from that of the IHR relative to the effective radius. We derive the two parameters from HST/ACS z {sub 850}-band images for a sample of 764 galaxies with z {sub 850} < 22 mag and 0.35 < z < 0.9 selected from the GEMS and GOODS-South surveys. We show that the sample galaxies having strong asymmetric structures, particularly tidal tails, are well-separated from those with regular morphologies in the A {sub o}-D {sub o} space. Meanwhile, the widely used CAS and Gini-M {sub 20} methods turn out to be insensitive to such morphological features. We stress that the A {sub o}-D {sub o} method is an efficient way to select galaxies with significant asymmetric features like tidal tails and study galaxy mergers in the dynamical phase traced by these delicate features.

79 ASTRONOMY AND ASTROPHYSICS↗

IAEA Coordinated Research Project on HTGR Physics, Thermal Hydraulics, and Depletion Uncertainty Analysis: PHISICS/RELAP5-3D Results for the Phase III Coupled Core Exercises

This report details the Parallel and Highly Innovative Simulation for INL Code System (PHISICS)/Reactor Excursions and Leak Analysis Program (RELAP5)-3D results obtained for the neutronics stand-alone core exercises defined for Phase II of the International Atomic Energy Agency (IAEA) Coordinated Research Project (CRP) on high-temperature gas cooled reactor (HTGR) uncertainty analysis in modeling (UAM). The Phase II models and results are linked to the earlier Standardized Computer Analyses for Licensing Evaluation (SCALE)/Sampler/New ESC-based Weighting Transport (NEWT) data generated for the lattice physics (lattice) stage Phase I of the CRP. The focus of this report is the propagation of cross-section uncertainties from the lattice to the core calculation (core) phase of the prismatic 350 megawatt (MW) General Atomics (GA) Modular HTGR (MHTGR). The Phase II neutronics exercises consist of two sub-cases (II-2a and II-2b) specified for the fresh and mixed core loading of the MHTGR-350 prismatic high-temperature reactor design. Exercises II-2a and II-2b are performed without temperature feedback using isothermal temperature profiles. For the nominal lattice physics (i.e., best-estimate) model, a series of fresh and depleted single block and super cell SCALE/KENO and SCALE/NEWT models were developed at Idaho National Laboratory (INL) to assess the impact of spectral effects on the few-group cross sections. A combination of these lattice cells were used in a super cell calculation to construct 8-group cross-section libraries for use in the INL-developed PHISICS/RELAP5-3D code suite. The block-level homogenized libraries were assigned to the 220 fuel zones utilized in the MHTGR-350 PHISICS core model, according to the fresh and mixed fresh/depleted core loading pattern. In order to evaluate the effect of uncertainty propagation, several sets consisting of 1,000 perturbed NEWT lattice and PHISICS/RELAP5-3D core calculations were performed using the SCALE 6.2 ENDF/B-VII.1 cross-section co-variance data and the stochastic sampling module Sampler in SCALE 6.2.0. The Figures of Merit (FOM) investigated include the mean and standard deviation uncertainty values obtained for the eigenvalues, control rod worth values, block-averaged power profiles, power axial offsets, and isothermal temperature feedback coefficients. A sensitivity study on the main nuclide reaction contributors to the FOM uncertainties has not yet been performed since the SCALE/PHISICS/RAVEN coupling with ENDF/B-VII co-variance data is still being developed. It was found that uncertainties in the cross sections lead to uncertainties in the local power values of up to 1.5%. The variances between the fresh and mixed/depleted core models are not as large as the axial variations between the top and bottom within each core model. The core eigenvalue uncertainty is approximately 0.5%, which is significant in comparison to typical eigenvalue difference targets (less than ~0.1% for LWRs). The power peaking and axial offset uncertainties are smaller (0.10% – 0.15%) and probably below the level of significance in typical core design metrics. The use of an 8-group energy structure slightly under-estimates some of the uncertainties, but in absolute terms, the uncertainty results are reasonably close to the 26-group results. The 26- and 8-group mean and nominal results differ by a larger degree, and if it can be afforded in terms of calculation resources, the 26-group structure should be used for HTGR calculations. The 8-group structure can, however, still be useful for reasonably accurate uncertainty assessments of transient cases where computational time is important. It was furthermore found that the use of a super cell during the cross-section generation process for better representation of the spectral environment at the core periphery did not lead to significant changes in the core eigenvalues, axial offset, and power peaking uncertainties. The impact of the super cell cross sections on the nominal and mean eigenvalues is, however, important enough (0.6% – 1.1% lower) to justify the use of these super cells for HTGR core simulation to represent the softer spectral environment in the core peripheral regions for nominal best-estimate calculations. The use of super cells for HTGR uncertainty studies is, therefore, not strictly necessary, but the improvements in nominal results could justify their use.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗