Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model checking”

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 811 records · Page 45

Deciding Between an X-Ray and 137Cs Irradiator – It’s Not Just About Energy Spectra

Irradiators utilizing radioactive cesium-137 or cobalt-60 have been used for medical applications for many decades. This includes the irradiation of much of the nation’s blood supply, and radiation biology research. In 2005, the U.S. Nuclear Regulatory Commission was assigned to evaluate and prevent the misuse of radioactive materials with malicious intentions, with gamma-ray sources being high on the list. Licensees who possessed higher quantities of radioactive material of concern were required to set up increased security controls such as constant surveillance, controlled access, personnel background checks, and other security controls. As a result of introducing such regulations, and problems with disposing of high activity sources, many users of gamma-ray sources started to look into alternative technologies to replace their gamma-ray irradiators – the main alternative being an X-ray irradiator. To make an educated decision on whether a particular type of X-ray irradiator is equivalent to a particular type of cesium-137 irradiator for specific applications, one has to rely on relevant published comparison studies from other researchers, or perform the comparison studies on their own. This publication focuses on the authors’ technical performance comparison of two cesium-137 irradiator models and two X-ray irradiator models, for the purpose of determining whether the X-ray irradiator models could validly replace the cesium-137 irradiator models for certain applications. Utilizing 47 gram and 120 gram tissue-equivalent rodent dosimetry phantoms, our results indicate that, depending upon the user’s dose uncertainty budget and maximum size of specimens to be irradiated, the Model RS2000 160 kVp X-ray irradiator, Model XRAD-160 X-ray irradiator, or Model XRAD-320 X-ray irradiator, could successfully replace a cesium-137 irradiator. It is expected such equivalency could be extrapolated to other irradiated specimens of known sizes and density. Results also reveal that differences in inherent source geometry, field geometry and irradiation geometry can counter some of the influence due to differences in energy spectrum. It is hoped that this publication can be used as a guide for other similar studies, providing investigators guidance on important details that can make the difference between strong and weak comparison conclusions.

Cs-137 irradiator, X-ray Irradiator, self-containe↗

THE INFRARED COLORS OF THE SUN

Solar infrared colors provide powerful constraints on the stellar effective temperature scale, but they must be measured with both accuracy and precision in order to do so. We fulfill this requirement by using line-depth ratios to derive in a model-independent way the infrared colors of the Sun, and we use the latter to test the zero point of the Casagrande et al. effective temperature scale, confirming its accuracy. Solar colors in the widely used Two Micron All Sky Survey (2MASS) JHK{sub s} and WISE W1-4 systems are provided: (V - J){sub Sun} = 1.198, (V - H){sub Sun} = 1.484, (V - K{sub s} ){sub Sun} = 1.560, (J - H){sub Sun} = 0.286, (J - K{sub s} ){sub Sun} = 0.362, (H - K{sub s} ){sub Sun} = 0.076, (V - W1){sub Sun} = 1.608, (V - W2){sub Sun} = 1.563, (V - W3){sub Sun} = 1.552, and (V - W4){sub Sun} = 1.604. A cross-check of the effective temperatures derived implementing 2MASS or WISE magnitudes in the infrared flux method confirms that the absolute calibration of the two systems agrees within the errors, possibly suggesting a 1% offset between the two, thus validating extant near- and mid-infrared absolute calibrations. While 2MASS magnitudes are usually well suited to derive T{sub eff}, we find that a number of bright, solar-like stars exhibit anomalous WISE colors. In most cases, this effect is spurious and can be attributed to lower-quality measurements, although for a couple of objects (3% {+-} 2% of the total sample) it might be real, and may hint at the presence of warm/hot debris disks.

79 ASTRONOMY AND ASTROPHYSICS↗

A TYPE II RADIO BURST WITHOUT A CORONAL MASS EJECTION

Type II radio bursts are thought to be a signature of coronal shocks. In this paper, we analyze a short-lived type II burst that started at 07:40 UT on 2011 February 28. By carefully checking white-light images, we find that the type II radio burst is not accompanied by a coronal mass ejection, only by a C2.4 class flare and narrow jet. However, in the EUV images provided by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory, we find a wave-like structure that propagated at a speed of ∼600 km s{sup −1} during the burst. The relationship between the type II radio burst and the wave-like structure is, in particular, explored. For this purpose, we first derive the density distribution under the wave by the differential emission measure method, which is used to restrict the empirical density model. We then use the restricted density model to invert the speed of the shock that produces the observed frequency drift rate in the dynamic spectrum. The inverted shock speed is similar to the speed of the wave-like structure. This implies that the wave-like structure is most likely a coronal shock that produces the type II radio burst. We also examine the evolution of the magnetic field in the flare-associated active region and find continuous flux emergence and cancellation taking place near the flare site. Based on these facts, we propose a new mechanism for the formation of the type II radio burst, i.e., the expansion of the strongly inclined magnetic loops after reconnecting with a nearby emerging flux acts as a piston to generate the shock wave.

79 ASTRONOMY AND ASTROPHYSICS↗

The impact of tillage row orientation on physical and chemical sediment enrichment

Abstract This study aimed to better understand how tillage row orientation with respect to dominant flow‐pathway along hillslope impacts runoff and the transport of different sediment size fractions. Experimental plots were constructed in contour ridge till (CRT) and parallel ridge till (PRT) sites to monitor runoff and sediment fluxes. Particle size fractions of the sediment, along with organic C and N contents were measured to quantify physical and chemical enrichment ratios for the two tillage orientations. In the CRT plot, tillage produced large oriented roughness elements along the contours, which acted as little check dams, while in the PRT site, the roughness elements helped confine and concentrate the runoff. In the CRT, the “check dams” resulted in a runoff coefficient of .03, while in the PRT, the flow confinement between rows produced a runoff coefficient of .82. Moreover, the erosion rates at the CRT site were 97% less than those in the PRT plot. Large contours produced finer sediment fractions due to the selective sorting in ponded furrows. More aggregate sediment fractions were present in the PRT site, which was dominated by rill erosion. Physical enrichment revealed selective entrainment of finer sediment particles during erosion. Finer‐grain particles with higher specific surface areas, attached more organic C resulting in chemical enrichment. Physical and chemical enrichment methods were found to be in good agreement. These findings suggest transport models that can simulate size fraction updates to the soil active layer can be used to estimate SOC redistribution and hence more accurate C budgets.

Wacha, Kenneth M.↗

PDBx/mmCIF Ecosystem: Foundational Semantic Tools for Structural Biology

PDBx/mmCIF, Protein Data Bank Exchange (PDBx) macromolecular Crystallographic Information Framework (mmCIF), has become the data standard for structural biology. With its early roots in the domain of small-molecule crystallography, PDBx/mmCIF provides an extensible data representation that is used for deposition, archiving, remediation, and public dissemination of experimentally determined three-dimensional (3D) structures of biological macromolecules by the Worldwide Protein Data Bank (wwPDB, wwpdb.org). Extensions of PDBx/mmCIF are similarly used for computed structure models by ModelArchive (modelarchive.org), integrative/hybrid structures by PDB-Dev (pdb-dev.wwpdb.org), small angle scattering data by Small Angle Scattering Biological Data Bank SASBDB (sasbdb.org), and for models computed generated with the AlphaFold 2.0 deep learning software suite (alphafold.ebi.ac.uk). Community-driven development of PDBx/mmCIF spans three decades, involving contributions from researchers, software and methods developers in structural sciences, data repository providers, scientific publishers, and professional societies. Having a semantically rich and extensible data framework for representing a wide range of structural biology experimental and computational results, combined with expertly curated 3D biostructure data sets in public repositories, accelerates the pace of scientific discovery. Herein, we describe the architecture of the PDBx/mmCIF data standard, tools used to maintain representations of the data standard, governance, and processes by which data content standards are extended, plus community tools/software libraries available for processing and checking the integrity of PDBx/mmCIF data. Use cases exemplify how the members of the Worldwide Protein Data Bank have used PDBx/mmCIF as the foundation for its pipeline for delivering Findable, Accessible, Interoperable, and Reusable (FAIR) data to many millions of users worldwide.

59 BASIC BIOLOGICAL SCIENCES↗

ON THE SECULAR BEHAVIOR OF DUST PARTICLES IN AN ECCENTRIC PROTOPLANETARY DISK WITH AN EMBEDDED MASSIVE GAS GIANT PLANET

We investigate the dust velocity and spatial distribution in an eccentric protoplanetary disk under the secular gravitational perturbation of an embedded planet of about 5 Jupiter masses. We first employ the FARGO code to obtain the two-dimensional density and velocity profiles of the eccentric gas disk exterior to the gap opened up by the embedded planet in the quasi-steady state. We then apply the secular perturbation theory and incorporate the gas drag to estimate the dust velocity and density on a secular timescale. The dust-to-gas ratio of the unperturbed disk is simply assumed to be 0.01. In our fiducial disk model with the planet at 5 AU, we find that 0.01 cm to 1 m sized dust particles are well coupled to the gas. Consequently, the particles behave similarly to the gas and exhibit asymmetric dynamics as a result of eccentric orbits. The dust surface density is enhanced around the apocenter of the disk. However, for the case of a low-density gaseous disk (called the 'transition disk' henceforth in this work) harboring the planet at 100 AU, the azimuthal distributions of dust of various sizes can deviate significantly. Overall, the asymmetric structure exhibits a phase correlation between the gas velocity fields and dust density distribution. Therefore, our study potentially provides a reality check as to whether an asymmetric disk gap detected at submillimeter and centimeter wavelengths is a signpost of a massive gas giant planet.

79 ASTRONOMY AND ASTROPHYSICS↗

Inputs to GCAM-USA: IM3 Phase 2 Experiments

Overview This dataset contains XML input files for the IM3 Phase 2 version of GCAM-USA. The files are organized into two categories: Scenario-specific inputs represent hydroclimate and socioeconomic effects on water availability, heating and cooling degree-hours, and agricultural productivity. They support eight IM3 canonical scenarios: rcp45cooler_ssp3 rcp45cooler_ssp5 rcp45hotter_ssp3 rcp45hotter_ssp5 rcp85cooler_ssp3 rcp85cooler_ssp5 rcp85hotter_ssp3 rcp85hotter_ssp5 Model-improvement inputs extend GCAM-USA v5.3 with updated representations of coal and nuclear power plant retirements, electricity trade among U.S. interconnections, offshore carbon storage costs, and groundwater depletion constraints. Data structure Scenario-specific inputs rcp45_runoff/ and rcp85_runoff/XML files describing water availability by HUC2 basin under the RCP 4.5 and RCP 8.5 scenarios. rcp45_hdcd/ and rcp85_hdcd/XML files containing monthly-day and monthly-night heating and cooling degree-hours at the U.S. state level for different RCP-SSP combinations. rcp45_agyields/ and rcp85_agyields/XML files describing changes in agricultural productivity at the intersection of GCAM regions and HUC2 water basins for different RCP-SSP combinations. rcp45_emissions_pathway/The emissions-constraint XML file used to represent the RCP 4.5 pathway. Model-improvement inputs core_retire/Updates coal-fired power plant retirement schedules based on New England ISO. GCAMUSA_IM3_elec_trade_interconnect.xmlRestricts electricity trade to occur within the ERCOT, WECC, and IE interconnections. nuclear_USA.xmlUpdates the retirement schedules of the Diablo Canyon and Palisades nuclear power plants. high_cost_offshore_carbon.xmlUpdates the assumed cost of offshore carbon storage. water_supply_constrained_gleeson_5pct.xmlReplaces WaterGAP historical groundwater-depletion estimates with data from the Gleeson dataset and limits groundwater extraction to 5% of the available groundwater in each Superwell grid cell. How to use the data This dataset is designed for use with the IM3 version of GCAM-USA. Download or clone GCAM-USA from the IM3 GCAM GitHub repository at https://github.com/IMMM-SFA/gcam-core and check out the gcam-usa-im3 branch. Place the downloaded folder im3scenarios in the gcam-core/input directory while preserving the provided folder structure.

Energy↗

Determination of 2D poloidal maps of the intrinsic W density for transport studies in JET-ILW

The experimental method developed at ASDEX Upgrade for the determination of the intrinsic tungsten (W) density profile coupling data from the soft X-ray (SXR) diagnostic and vacuum-ultra-violet (VUV) spectroscopy has been upgraded for application to JET plasmas. The strong poloidal asymmetries in the SXR emission are modeled assuming a ln ln ( ϵ ( ρ , R ) / ϵ ( ρ , R 0 ) ) = λ ( ρ ) ( R 2 − R 0 2 ) distribution, where ρ is the flux coordinate, R is the major radius, and λ is the fit parameter. The W density is calculated from the resulting 2D SXR emissivity maps accounting for contributions from a low-Z impurity (typically beryllium) and main ion with the assumption that their contributions are poloidally symmetric. Comparing the result with the independent W concentration measurement of VUV spectroscopy, a recalibration factor for the SXR emissivity is calculated making the method robust against the decrease in the sensitivity of the SXR diodes which has been observed across multiple campaigns. The final 2D W density map is checked for consistency versus the time-evolution of the W concentration measurement from VUV spectroscopy, toroidal rotation measurements from charge exchange recombination spectroscopy, and tomographic reconstructions of bolometry data. The method has been found to be robust for W concentrations above a few 10 -5 and in cases where the contributions from other medium-Z impurities such as Ni are negligible.

47 OTHER INSTRUMENTATION↗

A Runtime Verification Framework for Control System Simulation

n a standard workflow for the validation of a control system, the control system is implemented as an extension to a simulator. Such simulators are complex software systems, and engineers may unknowingly violate constraints a simulator places on extensions. As such, errors may be introduced in the implementation of either the control system or the simulator leading to invalid simulation results. This paper presents a novel runtime verification approach for verifying control system implementations within simulators. The major contribution of the approach is the two-tier specification process. In the first tier, engineers model constraints using a domain-specific language tailored to modeling a controller’s response to changes in its input. The language is high-level and effectively hides the implementation details of the simulator, allowing engineers to specify design-level constraints independent of low-level simulator interfaces. In the second tier, simulator developers provide mapping rules for mapping design-level constraints to the implementation of the simulator. Using the rules, an automated tool transforms the design-level specifications into simulator-specific runtime verification specifications and generates monitoring code which is injected into the implementation of the simulator. During simulation, these monitors observe the input and output variables of the control system and report changes to the verifier. The verifier checks whether these changes follow the constraints of the control system. We describe application of this approach to the verification of the constraints of an HVAC control system implemented with the power grid simulator GridLAB-D.

control system↗

Cross section measurements for neutron inelastic scattering and the ( n , 2 n γ ) reaction on Pb 206

We measured excitation functions for γ production associated with the neutron inelastic scattering and the (n, 2n) reactions on 206 Pb from threshold up to 18 MeV for about 40 transitions. Two independent measurements were performed using different samples and acquisition systems to check consistency of the results. Moreover, the neutron flux was determined with a 235 U fission chamber and a procedure that were validated against a fluence standard. For incident energy higher than the threshold for the first excited level and up to 3.5 MeV, estimates are provided for the total inelastic and level cross sections by combining the present γ production cross sections with the level and decay data of 206 Pb reported in the literature. The uncertainty common to all incident energies is 3.0% allowing overall uncertainties from 3.3% to 30% depending on transition and neutron energy. Finally, the present data agree well with earlier work, but significantly expand the experimental database while comparisons with model calculations using the talys reaction code show good agreement over the full energy range.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

EXPLAINING THE MOST ENERGETIC SUPERNOVAE WITH AN INEFFICIENT JET-FEEDBACK MECHANISM

We suggest that the energetic radiation from core-collapse super-energetic supernovae (SESNe) is due to a long-lasting accretion process onto the newly born neutron star (NS), resulting from an inefficient operation of the jet-feedback mechanism (JFM). The jets that are launched by the accreting NS or black hole maintain their axis due to a rapidly rotating pre-collapse core and do not manage to eject core material from near the equatorial plane. The jets are able to eject material from the core along the polar directions and reduce the gravity near the equatorial plane. The equatorial gas expands, and part of it falls back over a timescale of minutes to days to prolong the jet-launching episode. According to the model for SESNe proposed in the present paper, the principal parameter that distinguishes between the different cases of core-collapse supernova (CCSN) explosions, such as between normal CCSNe and SESNe, is the efficiency of the JFM. This efficiency, in turn, depends on the pre-collapse core mass, envelope mass, core convection, and, most of all, the angular momentum profile in the core. One prediction of the inefficient JFM for SESNe is the formation of a slow equatorial outflow in the explosion. The typical velocity and mass of this outflow are estimated to be v {sub eq} ≈ 1000 km s{sup −1} and M {sub eq} ≳ 1 M {sub ⊙}, respectively, though quantitative values will have to be checked in future hydrodynamic simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Chandrasekhar's relation and stellar rotation in the Kepler field

According to the statistical law of large numbers, the expected mean of identically distributed random variables of a sample tends toward the actual mean as the sample increases. Under this law, it is possible to test the Chandrasekhar's relation (CR), (V) = (π/4){sup –1}(Vsin i), using a large amount of Vsin i and V data from different samples of similar stars. In this context, we conducted a statistical test to check the consistency of the CR in the Kepler field. In order to achieve this, we use three large samples of V obtained from Kepler rotation periods and a homogeneous control sample of Vsin i to overcome the scarcity of Vsin i data for stars in the Kepler field. We used the bootstrap-resampling method to estimate the mean rotations ((V) and (Vsin i)) and their corresponding confidence intervals for the stars segregated by effective temperature. Then, we compared the estimated means to check the consistency of CR, and analyzed the influence of the uncertainties in radii measurements, and possible selection effects. We found that the CR with (sin i) = π/4 is consistent with the behavior of the (V) as a function of (Vsin i) for stars from the Kepler field as there is a very good agreement between such a relation and the data.

79 ASTRONOMY AND ASTROPHYSICS↗

Two-dimensional forms of robust CO 2 reduction photocatalysts

Photoelectrocatalysts that use sunlight to power the CO 2 reduction reaction will be crucial for carbon-neutral power and energy-efficient industrial processes. Scalable photoelectrocatalysts must satisfy a stringent set of criteria, such as stability under operating conditions, product selectivity, and efficient light absorption. Two-dimensional materials can offer high specific surface area, tunability, and potential for heterostructuring, providing a fresh landscape of candidate catalysts. From a set of promising bulk CO 2 reduction photoelectrocatalysts, we screen for candidate monolayers of these materials, then study their catalytic feasibility and suitability. For stable monolayer candidates, we verify the presence of visible-light band gaps, check that band edges can support CO 2 reduction, determine exciton binding energies, and compute surface reactivity. We find visible light absorption for SiAs, ZnTe, and ZnSe monolayers, and that due to a lack of binding, CO selectivity is possible. We thus identify SiAs, ZnTe, and ZnSe monolayers as targets for further investigation, expanding the chemical space for CO 2 photoreduction candidates.

36 MATERIALS SCIENCE↗

STOCHASTIC DESCRIPTION OF THE HIGH-FREQUENCY CONTENT OF DAILY SUNSPOTS AND EVIDENCE FOR REGIME CHANGES

The irregularity index λ is applied to the high-frequency content of daily sunspot numbers ISSN. This λ is a modification of the standard maximal Lyapunov exponent. It is computed here as a function of embedding dimension m, within four-year time windows centered at the maxima of Schwabe cycles. The λ(m) curves form separate clusters (pre-1923 and post-1933). This supports a regime transition and narrows its occurrence to cycle 16, preceding the growth of activity leading to the Modern Maximum. The two regimes are reproduced by a simple autoregressive process AR(1), with the mean of Poisson noise undergoing 11 yr modulation. The autocorrelation a of the process (linked to sunspot lifetime) is a ≈ 0.8 for 1850-1923 and ≈0.95 for 1933-2013. The AR(1) model suggests that groups of spots appear with a Poisson rate and disappear at a constant rate. We further applied the irregularity index to the daily sunspot group number series for the northern and southern hemispheres, provided by the Greenwich Royal Observatory (RGO), in order to study a possible desynchronization. Correlations between the north and south λ(m) curves vary quite strongly with time and indeed show desynchronization. This may reflect a slow change in the dimension of an underlying dynamical system. The ISSN and RGO series of group numbers do not imply an identical mechanism, but both uncover a regime change at a similar time. Computation of the irregularity index near the maximum of cycle 24 will help in checking whether yet another regime change is under way.

79 ASTRONOMY AND ASTROPHYSICS↗

Interpretation of Lyman opacity measurements in JET with the ITER-like wall using a particle balance approach

Recent measurements of the Ly β /D α ratio in the JET ITER-like wall (ILW) divertor show comparatively more pronounced Ly β reabsorption relative to previous results in JET with the carbon wall (JET-C). At the outer horizontal target ion current rollover point a Ly β reabsorption rate of 60% is measured, increasing to 80% in more pronounced detachment, implying a Ly α reabsorption rate of 90%–98%. The radially resolved Ly β /D α measurements are used to constrain Lyman opacity corrections to atomic rate coefficients using the population escape factor technique in order to capture the local changes to the excited state population structure in the high Lyman opacity regions at and outboard of the outer strike point. To check the self-consistency of the Lyman opacity measurements, a detailed spectroscopic interpretation of the outer divertor particle balance is presented, in which the impact of opacity corrections to the Ly α inverse photon efficiency coefficients is assessed. A five-fold deficit in the estimated D + source rate obtained with optically thin plasma assumptions is reconciled once the opacity corrections are factored into the Ly α photon rate to ionization rate conversion. The experiment results are reproduced in EDGE2D-EIRENE density scan simulations in which an imposed ad hoc Lyman reabsorption rate is recovered spectroscopically using synthetic measurements, and a similar shortfall in the D + source rate estimates is reconciled using Lyman opacity corrected atomic data. The model limitations prevent a more detailed self-consistent analysis of the Lyman opacity impact on divertor parameters and detachment evolution, and hence motivate renewed efforts to re-establish routine exploitation of the photon transport modelling capabilities in the EIRENE code package. Lyman opacity corrections to atomic data coefficients are necessary for spectroscopic interpretation of the JET-ILW divertor plasma, with significant influence on the divertor plasma also a likely consequence of the presence of strong Ly α reabsorption. Finally, a scan of auxiliary heating powers spanning L-mode and H-mode conditions reveals a strong correlation of the measured Ly β opacity with the outer target temperature, suggesting new possibilities for using Lyman opacity measurements in establishing detachment scalings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Spring Boot based REST API to Improve Data Quality Report Generation for Big Scientific Data: ARM Data Center Example

Web application technologies are growing rapidly with continuous innovation and improvements. This paper focuses on the popular Spring Boot [1] java-based framework for building web and enterprise applications and how it provides the flexibility for service-oriented architecture (SOA). One challenge with any Spring-based applications is its level of complexity with configurations. Spring Boot makes it easy to create and deploy stand-alone, production-grade Spring applications with very little Spring configuration. Example, if we consider Spring Model-View-Controller (MVC) framework [2], we need to configure dispatcher servlet, web jars, a view resolver, and component scan among other things. To solve this, Spring Boot provides several Auto Configuration options to setup the application with any needed dependencies. Another challenge is to identify the framework dependencies and associated library versions required to develop a web application. Spring Boot offers simpler dependency management by using a comprehensive, but flexible, framework and the associated libraries in one single dependency, which provides all the Spring related technology that you need for starter projects as compared to CRUD web applications. This framework provides a range of additional features that are common across many projects such as embedded server, security, metrics, health checks, and externalized configuration. Web applications are generally packaged as war and deployed to a web server, but Spring Boot application can be packaged either as war or jar file, which allows to run the application without the need to install and/or configure on the application server. In this paper, we discuss how Atmospheric Radiation Measurement (ARM) Data Center (ADC) at Oak Ridge National Laboratory, is using Spring Boot to create a SOA based REST [4] service API, that bridges the gap between frontend user interfaces and backend database. Using this REST service API, ARM scientists are now able to submit reports via a user form or a command line interface, which captures the same data quality or other important information about ARM data.

Guntupally, Kavya↗

Verification of the DIF3D Software to Support Fast Reactor Analysis

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the polynomial order used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Verification of the DIF3D Software to Support Fast Reactor Analysis

Ongoing design activities at Argonne National Laboratory are requiring a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system. The driver for this effort requires the 3D Cartesian, triangular-Z, and hexagonal-Z core geometry options of DIF3D be verified. Previous work identified the DIF3D features required to be verified to support current design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying DIF3D’s ability to correctly translate the user’s model in to DIF3D’s preferred format, verifying that options planned for use have the desired effect, and verifying the correctness of the eigenvalue, fixed-source, forward, and adjoint solvers in DIF3D-FD and DIF3D-VARIANT. This manuscript provides the verification tasks and their results with respect to the features needed for current design activities. Since analytic solutions of the neutron diffusion and transport equations are either limited in scope or not possible, multiple tiers of problems unique to each solver and geometry type were implemented. Each of these tiers tests features independent and complementary arguments for why the separate testing of functionalities is acceptable. Finally, this separate testing was also supplemented with a high-level integral check of each the diffusion and transport capabilities and applicable geometries. To accommodate cases which an analytic solution is not feasible, MCNP6.2 was relied upon to provide a higher-order reference solution. This therefore required that the capabilities within MCNP6.2 which were relied upon for this work are also verified in this work. No MCNP discrepancies were noted in this effort. Note that the MCNP6.2 verification included in this work does not stand as a full verification of MCNP6.2, but merely verifies the features used in verifying DIF3D. The verification effort identified no issues that are debilitating or otherwise impactful to design usage of DIF3D, and thus DIF3D version 11.0, release 3012 is considered verified. As some additional changes have been made to the ARC software since this point all versions between release 3012 and 3266 can be considered verified as version 3253 was used for all updates in this revision. The types of issues that were identified were predominantly in the areas of: unclear documentation, software bugs which were inconsequential to final results, editing options which were ignored in favor of printing more information than requested, bugs in the outputs of intermediate results, or secondary output binary file information which was not present. While not a bug, this verification report also identified that the algorithm used to evaluate the peak fast flux in a nodal transport solution can be quite unreliable due to the methodology used and the location of the peak within the mesh. The authors of the report therefore recommend the usage of the EvaluateFlux software (distributed with ARC) as a more robust alternative noting that DIF3D will properly notify the user when the peaking values it is providing are potentially incorrect.

97 MATHEMATICS AND COMPUTING↗