Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model benchmarking”

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 343 records · Page 19

Evaluated benchmark experiments at critical assemblies simulating features of an HTHR at the ASTRA facility

The design of the ASTRA facility and critical assemblies that simulate physics features of modular high-temperature reactors (HTHR-Ms) with a graphite moderator and reflectors loaded with fuel particles having multilayer ceramic coatings is described in detail. Geometrical dimensions of the main elements and regions of the critical assemblies, composition of the materials used, and experimental results for various configurations of the critical assemblies are presented. A detailed computational benchmark model allowing for the structural and compositional features of the critical assembly configurations in question is developed on the basis of all the above data. The results are to be used for verification of the neutronics codes used for calculations of high-temperature helium-cooled reactors.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Validation and bias quantification of criticality safety codes for SRNS operations

The validation and bias quantification methods employed by the Criticality Safety Engineering group at Savannah River Nuclear Solutions (SRNS) are discussed. Both SCALE and MCNP are validated against benchmark experiments, where these benchmark models have been divided into subgroups representative of the types of fissile material operations at SRNS. A quantification of the bias is made for each subgroup for each code and the three methods employed for quantification are discussed. The bias information is employed directly in the establishment of the Upper Safety Limit (USL) value of the neutron multiplication factor when using standardized documented SRNS methods. A simple example of validation and bias quantification is given and the procedure for setting the USL using the derived bias information is also presented. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Development of an aerosol microphysical module: Aerosol Two-dimensional bin module for foRmation and Aging Simulation (ATRAS)

Number concentrations, size distributions, and mixing states of aerosols are essential parameters for accurate estimation of aerosol direct and indirect effects. In this study, we developed an aerosol module, designated Aerosol Two-dimensional bin module for foRmation and Aging Simulation (ATRAS), that can represent these parameters explicitly by considering new particle formation (NPF), black carbon (BC) aging, and secondary organic aerosol (SOA) processes. A two-dimensional bin representation is used for particles with dry diameters from 40 nm to 10 µm to resolve both aerosol size (12 bins) and BC mixing state (10 bins) for a total of 120 bins. The particles with diameters from 1 to 40 nm are resolved using an additional 8 size bins to calculate NPF. The ATRAS module was implemented in the WRF-chem model and applied to examine the sensitivity of simulated mass, number, size distributions, and optical and radiative parameters of aerosols to NPF, BC aging and SOA processes over East Asia during the spring of 2009. BC absorption enhancement by coating materials was about 50% over East Asia during the spring, and the contribution of SOA processes to the absorption enhancement was estimated to be 10 – 20% over northern East Asia and 20 – 35% over southern East Asia. A clear north-south contrast was also found between the impacts of NPF and SOA processes on cloud condensation nuclei (CCN) concentrations: NPF increased CCN concentrations at higher supersaturations (smaller particles) over northern East Asia, whereas SOA increased CCN concentrations at lower supersaturations (larger particles) over southern East Asia. Application of ATRAS to East Asia also showed that the impact of each process on each optical and radiative parameter depended strongly on the process and the parameter in question. The module can be used in the future as a benchmark model to evaluate the accuracy of simpler aerosol models and examine interactions between NPF, BC aging, and SOA processes under different meteorological conditions and emissions.

Matsui, H.↗

Comparative analysis of energy deposition modes available in Serpent 2 within the framework of the supercritical water reactor - Fuel qualification test reactor physics benchmark

A joint European Canadian Chinese development of a supercritical water-cooled small modular reactor (SCW-SMR) technology is in progress since September 2020 in the framework of a Horizon 2020 project called ECC-SMART. As a main purpose of the project, proper estimates of energy deposition and its spatial distribution are prerequisites for the accurate analysis of safety related parameters of the SCW-SMR concept under development. A supercritical water reactor fuel computational benchmark model, provided by Canadian Nuclear Laboratories, was applied for detailed comparison of different energy deposition calculation options available in the Serpent 2 Monte Carlo code. The effect of energy deposition options on the normalization of the results as well as on the spatial distribution of the energy deposition are discussed. Consistent energy deposition calculation methods are presented between three Monte Carlo codes, viz., Serpent 2, MCNP6 and OpenMC. Although resource-intensive, the use of the coupled neutron-photon transport mode of Serpent 2 is recommended for accurate spatial and quantitative characterization of energy deposition in the SCW-SMR fuel assemblies, accounting for both neutron and photon heating of all the materials. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Major improvements on the longwave radiative interactions between surface and clouds in the polar regions in atmospheric global circulation model

Motivations: One of the most important tasks in the study of physical and socioeconomic aspects of climate change is the simulation of polar climate and how it changes in response to global warming. One import component of the polar energy budget (and thus of polar climate) is the radiative interaction between surface and clouds. As of today, the vast majority of state-of-the-art global climate models (GCMs), including all three US flagship GCMs (NCAR CESM, GFDL CM3, and GISS E2-H/E2-R), still assume constant surface emissivity overall spectral bands in their longwave radiation treatment. Similarly, a majority of state-of-the-art GCMs assume non-scattering clouds in the longwave spectrum, including both the NCAR CESM and GFDL CM3. The issues of ignoring spectral variation in surface emissivity and scattering by clouds manifest themselves in the far-infrared (IR) over the polar continents because (1) there is so little water vapor in such areas that surface far-IR emission can reach clouds; and (2) the imaginary part of the refractive index of ice has a local minimum over 350-550 cm-1, which implies possible strong scattering effects over this spectral region both within and between the snow surface and ice clouds. Our sensitivity studies using CloudSat-retrieved hydrometeor profiles over the high-elevation Antarctic continent show that, in winter, an appropriate inclusion of snow surface spectral emissivity and ice cloud scattering in radiative transfer calculations can noticeably reduce the monthly-mean surface net downward far-IR flux and net atmospheric far-IR emission over the entire region. For the far-IR alone, the magnitude of these effects is ~1Wm-2 or even larger. Projection Objectives: We propose to carry out the following studies: (1) Develop or update schemes that can be used in the CESM to compute the optical properties of the snow surface and ice clouds consistently across the longwave and the shortwave spectrum; (2) Modify the longwave radiation scheme in the CESM such that it can incorporate surface spectral emissivity and scattering of clouds; (3) Validate schemes developed for (1) and (2) against measurements or benchmark models; (4) Evaluate the effects of these new treatments on the simulated polar energy budget and climate by the CESM. Technical Approaches:1.Use a detailed snowpack radiative transfer model to simulate spectrally resolved emissivity of snow surfaces as a function of snow effective radius. The model will be validated against available high spectral-resolution measurements of snow spectral emissivity (ASTER spectra library v2.0). Currently the CESM prognoses vertically-resolved snow effective radius, and snow shortwave albedo is parameterized based on snow effective radius. Our treatment of snow spectral emissivity here will be consistent with the treatment of snow albedo in the shortwave and seamlessly connected with the snow physical module in the CESM. 2. Update the parameterization of ice cloud optical properties in the CESM, for both longwave and shortwave, based on the latest ice-cloud scattering database developed by co-I Prof. Yang. 3. Modify RRTMG_LW, the current longwave radiation scheme in the CESM, to include the surface spectral emissivity and the scattering of clouds in the longwave calculation. Validate modified scheme against a line-by-line radiative transfer model with multiple scattering capability (LBLRTM/CHARTS). 4. Implement products from Steps 1-3 into the CESM and conduct simulations to evaluate their effects on the simulated polar climate and energy budget. These runs will be compared to control runs and evaluated using routine model-observation diagnostics, including data from the ARM NSA (North Slope of Alaska) site. Moreover, a novel aspect in this step is that we will evaluate the simulated longwave band-by-band fluxes and cloud radiative effect against observational counterparts over both polar continents, using a product recently derived by the lead PI Prof. Huang from collocated AIRS and CERES satellite observations from 2002 to 2013. Benefits and outcomes: This project will significantly improve the longwave radiation scheme in the CESM and its successors and provide more realistic treatment of radiative interactions between polar surfaces and clouds. In the future, it will also be straightforward to extend this framework to include spectral emissivity of all surface types. This is responsive to identified DoE Earth System Modeling Science Drivers to improve the performance of the fully coupled climate system modeling and to help address questions on how rapid changes in cryospheric systems interact with the climate system.

54 ENVIRONMENTAL SCIENCES↗

Di-Higgs production in the 4b channel and gravitational wave complementarity

We present a complementarity study of gravitational waves and double Higgs production in the 4 b channel, exploring the gauge singlet scalar extension of the SM. This new physics extension serves as a simplified benchmark model that realizes a strongly first-order electroweak phase transition necessary to generate the observed baryon asymmetry in the universe. In calculating the signal-to-noise ratio of the gravitational waves, we incorporate the effect of the recently discovered significant suppression of the gravitational wave signals from sound waves for strong phase transitions, make sure that supercooled phase transitions do complete and adopt a bubble wall velocity that is consistent with a successful electroweak baryogenesis by solving the velocity profiles of the plasma. The high-luminosity LHC sensitivity to the singlet scalar extension of the SM is estimated using a shape-based analysis of the invariant 4 b mass distribution. We find that while the region of parameter space giving detectable gravitational waves is shrunk due to the new gravitational wave simulations, the qualitative complementary role of gravitational waves and collider searches remain unchanged.

79 ASTRONOMY AND ASTROPHYSICS↗

Search for a heavy resonance decaying into a top quark and a W boson in the lepton+jets final state at $\sqrt{s}$ = 13 TeV

A search for a heavy resonance decaying into a top quark and a W boson in proton-proton collisions at $\sqrt{s}$ = 13 TeV is presented. The data analyzed were recorded with the CMS detector at the LHC and correspond to an integrated luminosity of 138 fb -1 . The top quark is reconstructed as a single jet and the W boson, from its decay into an electron or muon and the corresponding neutrino. A top quark tagging technique based on jet clustering with a variable distance parameter and simultaneous jet grooming is used to identify jets from the collimated top quark decay. The results are interpreted in the context of two benchmark models, where the heavy resonance is either an excited bottom quark b* or a vector-like quark B. A statistical combination with an earlier search by the CMS Collaboration in the all-hadronic final state is performed to place upper cross section limits on these two models. The new analysis extends the lower range of resonance mass probed from 1.4 down to 0.7 TeV. For left-handed, right-handed, and vector-like couplings, b* masses up to 3.0, 3.0, and 3.2 TeV are excluded at 95% confidence level, respectively. The observed upper limits represent the most stringent constraints on the b* model to date.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

How resonance-continuum interference changes 750 GeV diphoton excess: Signal enhancement and peak shift

A hypothetical new scalar resonance, a candidate explanation for the recently observed 750 GeV diphoton excess at the LHC 13 TeV, necessarily interferes with the continuum background gg → γγ. The interference has two considerable effects: (1) enhancing or suppressing diphoton signal rate due to the imaginary-part interference and (2) distorting resonance shape due to the real-part interference. We study them based on the best-fit analysis of two benchmark models: two Higgs doublets with ~50 GeV width (exhibiting the imaginary-part interference effect) and a singlet scalar with 5 GeV width (exhibiting the real-part one), both extended with vector-like fermions. Furthermore, we find that the resonance contribution can be enhanced by a factor of 2 (1.6) for 3 (6) fb signal rate, or the 68% CL allowed mass region is shifted by O (1) GeV. If the best-fit excess rate decreases in the future data, the interference effects will become more significant.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for leptoquarks decaying into the bτ final state in pp collisions at $\sqrt{\textrm{s}}$ = 13 TeV with the ATLAS detector

A search for leptoquarks decaying into the $bτ$ final state is performed using Run 2 proton-proton collision data from the Large Hadron Collider, corresponding to an integrated luminosity of 139 fb —1 at $\sqrt{\textrm{s}}$ = 13 TeV recorded by the ATLAS detector. The benchmark models considered in this search are vector leptoquarks with electric charge of 2/3e and scalar leptoquarks with an electric charge of 4/3e. No significant excess above the Standard Model prediction is observed, and 95% confidence level upper limits are set on the cross-section times branching fraction of leptoquarks decaying into bτ. For the vector leptoquark production two models are considered: the Yang-Mills and Minimal coupling models. In the Yang-Mills (Minimal coupling) scenario, vector leptoquarks with a mass below 1.58 (1.35) TeV are excluded for a gauge coupling of 1.0 and below 2.05 (1.99) TeV for a gauge coupling of 2.5. In the case of scalar leptoquarks, masses below 1.28 (1.53) TeV are excluded for a Yukawa coupling of 1.0 (2.5). Finally, an interpretation of the results with minimal model dependence is performed for each of the signal region categories, and limits on the visible cross-section for beyond the Standard Model processes are provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for magnetic monopoles and stable particles with high electric charges in $\sqrt{s}$ = 13 TeV pp collisions with the ATLAS detector

We present a search for magnetic monopoles and high-electric-charge objects using LHC Run 2 $\sqrt{s}$ = 13 TeV proton-proton collisions recorded by the ATLAS detector. A total integrated luminosity of 138 fb –1 was collected by a specialized trigger. No highly ionizing particle candidate was observed. Considering the Drell-Yan and photon-fusion pair production mechanisms as benchmark models, cross-section upper limits are presented for spin-0 and spin-1/2 magnetic monopoles of magnetic charge 1g D and 2g D and for high-electric-charge objects of electric charge 20 ≤ |z| ≤ 100, for masses between 200 GeV and 4000 GeV. The search improves by approximately a factor of three the previous cross-section limits on the Drell-Yan production of magnetic monopoles and high-electric charge objects. Also, the first ATLAS limits on the photon-fusion pair production mechanism of magnetic monopoles and high-electric-charge objects are obtained.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Implementation of the direct S ( α , β ) method in the KENO Monte Carlo code

The Monte Carlo code KENO contains thermal scattering data for a wide variety of thermal moderators. These data are processed from Evaluated Nuclear Data Files (ENDF) by AMPX and stored as double differential probability distribution functions. The method examined in this study uses S(α,β) probability distribution functions derived from the ENDF data files directly instead of being converted to double differential cross sections. This allows the size of the cross section data on the disk to be reduced substantially amount. KENO has also been updated to allow interpolation in temperature on these data so that problems can be run at any temperature. Results are shown for several simplified problems for a variety of moderators. In addition, benchmark models based on the KRITZ reactor in Sweden were run, and the results are compared with the previous versions of KENO without the direct S(α,β) method. Results from the direct S(α,β) method compare favorably with the original results obtained using the double differential cross sections. Finally, sampling the data increases the run-time of the Monte Carlo calculation, but memory usage is decreased substantially.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scalable DPG multigrid solver for Helmholtz problems: A study on convergence

This paper presents a scalable multigrid preconditioner targeting large-scale systems arising from discontinuous Petrov–Galerkin (DPG) discretizations of high-frequency wave operators. This work is built on previously developed multigrid preconditioning techniques of Petrides and Demkowicz (Comput. Math. Appl. 87 (2021) pp. 12–26) and extends the convergence results from $\mathscr{O}$(10 7 ) degrees of freedom (DOFs) to $\mathscr{O}$(10 9 ) DOFs using a new scalable parallel MPI/OpenMP implementation. Novel contributions of this paper include an alternative definition of coarse-grid systems based on restriction of fine-grid operators, yielding superior convergence results. In the uniform refinement setting, a detailed convergence study is provided, demonstrating h and p robust convergence and linear scaling with respect to the wave frequency. Finally, the paper concludes with numerical results on hp -adaptive simulations including a large-scale seismic modeling benchmark problem with high material contrast.

97 MATHEMATICS AND COMPUTING↗

A two-stage optical fusion framework for wildfire severity mapping across the conterminous United States

Accurate wildfire severity mapping (WSM) is essential for post-fire recovery planning, erosion risk assessment, ecosystem monitoring, and disaster risk reduction. Although Landsat and Sentinel optical imagery have been widely used for burn severity assessment, the added value of fusing multiple optical sensors has not been sufficiently quantified across diverse fire events, particularly since the launch of Landsat-9. This study evaluates whether multisensor optical fusion improves wildfire severity mapping relative to single-sensor baselines using Sentinel-2, Landsat-8, and Landsat-9 imagery across 40 wildfire events in the conterminous United States. We tested a two-stage fusion framework that combines feature-level fusion with pixel-level dimensionality reduction. First, feature-level fused datasets were created through early fusion by combining standardized post-fire bands from each sensor into a single predictor stack. Both raw reflectance bands and pairwise spectral transforms were retained to capture within- and cross-sensor spectral interactions. Second, Linear Discriminant Analysis was applied to both single-sensor and fused datasets to produce comparable low-dimensional feature spaces. Six machine-learning classifiers were then used to benchmark model performance with repeated spatially buffered train–test splits. Results show that Landsat-9 was the strongest single-sensor baseline. Among the fusion strategies, Sentinel-2 + Landsat-9 produced the most consistent improvement and reduced performance variability. Landscape-condition analysis further showed that this fusion was most beneficial in shrubland-dominated and high-terrain fires, where it achieved the highest overall mean accuracy and the fewest failures. In contrast, its benefits were less reliable in evergreen forests, mixed vegetation, and low- to moderate-elevation terrain. In operational settings, the Sentinel-2 + Landsat-9 configuration offers a practical solution for post-fire recovery planning, erosion-risk assessment, watershed management, and ecological monitoring when field observations are available and timely satellite-based information is needed.

Landsat↗

Modernizing the active neutron collar poison rod correction using simulations to enhance PWR fresh fuel verification

The declared inventory of 235 U in fresh, low enriched uranium, nuclear fuel assemblies is routinely verified by nuclear safeguards inspectorates using the UNCL (Uranium Neutron Collar – Light Water Reactor Fuel) nondestructive assay instrument. The trend in modern fuels is towards higher initial enrichment, which in turn requires a larger number of burnable poison pins with higher Gadolinia concentration to hold down the initial reactivity. UNCL assay error for modern fuel assemblies typically exceeds 10%. The traditional algorithm used to correct the response for burnable poison content needs revision to achieve accuracy comparable to the $\sigma_\text{R}$ of 4.5% expected for poison-free assemblies. Here we review the problem and use a large set of 287 Monte Carlo simulations based on an experimentally benchmarked model corresponding to the 16×16 PWR array used in Brazilian Angra type II and III fuel. Simulated relative responses are used to evaluate the functional form of the poison correction and determine parameters optimal for 16×16 PWR assemblies. This update is conducted using between 4–24 Gadolinia pins ranging 2–11 wt% Gd 2 O 3 with assemblies having mean enrichments ranging 2.5%–5%. Experimental validation of the updated coefficients using seven measurements shows bias is reduced by an order of magnitude and $\sigma_\text{R}$R reduced from 10% to less than 2%. Two updated sets of coefficients are provided, the most accurate of which is directly useable in the existing analysis code (INCC) in use by inspectorate at fuel fabrication facilities.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Search for resonances decaying into photon pairs in 139 fb -1 of pp collisions at s=13 TeV with the ATLAS detector

Searches for new resonances in the diphoton final state, with spin 0 as predicted by theories with an extended Higgs sector and with spin 2 using a warped extra-dimension benchmark model, are presented using 139 fb$^{-1}$ of $\sqrt{s} = $ 13 TeV $pp$ collision data collected by the ATLAS experiment at the LHC. No significant deviation from the Standard Model is observed and upper limits are placed on the production cross-section times branching ratio to two photons as a function of the resonance mass.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Search for a new heavy resonance decaying to a top quark and a neutral scalar boson in proton-proton collisions at $\sqrt{s}$ = 13 TeV

A first search at the LHC for a new heavy resonance decaying to a top quark and a neutral scalar boson ϕ in the fully hadronic final state is presented, where the ϕ boson is identified by its decay into a bottom quark-antiquark pair. The search is focused on final states in which the decay products of the highly Lorentz boosted top quark and ϕ boson are each reconstructed as a single, large-radius jet with distinct substructure. The analysis is performed using proton-proton collision data at a center-of-mass energy of 13 TeV, corresponding to an integrated luminosity of 138 fb −1 , recorded with the CMS detector at the LHC in 2016–2018. The single production of a vector-like top quark, T′, is used as a benchmark model for the signal process. The results of this search are combined with those of a previous CMS search in which semileptonic decays of the top quark were used. No significant excess of data is observed with respect to the background prediction. For the case where the neutral scalar is a standard model Higgs boson and the T′ quark width is 5% of its mass, T′ quark masses between 0.85 and 1.3 TeV are excluded at 95% confidence level and the most stringent limits to date are set for masses above 2 TeV . For other ϕ boson masses, upper limits as low as 0.1 fb are set on the product of the T′ quark production cross section and branching fraction for its decay to a top quark and a ϕ boson.

Boosted↗

Depletion benchmark for a high-assay low-enriched uranium fuel experiment in the advanced test reactor

Reactor physics depletion benchmarks for high-assay low-enriched uranium (HALEU) fuel are limited in number. In particular, there is limited data for HALEU benchmarks for U-10Mo (uranium-10% molybdenum) plate fuel that is being developed for use in the United States’ high performance research reactors including the Advanced Test Reactor (ATR), Advanced Test Reactor Critical Facility (ATR-C), High Flux Isotope Reactor (HFIR), Massachusetts Institute of Technology Reactor (MITR), University of Missouri Research Reactor (MURR), National Bureau of Standards Reactor (NBSR). These six reactors currently operate with highly enriched uranium dispersed fuel in an aluminum matrix. In support of conversion to a HALEU fuel, qualification of U-10Mo formed into a monolithic foil is being performed. Fuel qualification involves irradiating fuel specimens in the ATR. The irradiation tests provide an opportunity to benchmark depletion capabilities of reactor physics codes in support of the ATR operation, as well as develop benchmarks that can be used by other institutions to benchmark other reactor physics codes. This paper documents the development of a benchmark model of the irradiation of the ATR Full-size plate In center flux trap Position 7 (AFIP-7) experiment using the depletion codes MC21 and Advanced Dimensional Depletion for Engineering of Reactors (ADDER).

Nielsen, Joseph W. [Idaho National Laboratory (INL↗