Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solution multiplicity”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Unsteady panel method for flows with multiple bodies moving along various paths

A potential flow based three-dimensional panel method was modified to treat time-dependent conditions in which several submerged bodies can move within the fluid along different trajectories. This modification was accomplished by formulating the momentary solution in an inertial frame of reference, attached to the undisturbed stationary fluid. Consequently, the numerical interpretation of the multiple-body, solid-surface boundary condition and the viscous wake rollup was considerably simplified. The usteady capability of this code was calibrated and validated by comparing computed results with closed-form analytical results available for an airfoil, which was impulsively set into a constant speed forward motion. To demonstrate the multicomponent capability, computations were made for two wings following closely intersecting paths (i.e., simulations aimed at avoiding mid-air collisions) and for a flowfield with relative rotation (i.e., the case of a helicopter rotor rotating relative to the fuselage). Computed results for the cases were compared to experimental data, when such data was available.

Richardson, Thomas F.↗
NASA NTRS1978

Aerodynamic preliminary analysis system. Part 1: Theory

A comprehensive aerodynamic analysis program based on linearized potential theory is described. The solution treats thickness and attitude problems at subsonic and supersonic speeds. Three dimensional configurations with or without jet flaps having multiple non-planar surfaces of arbitrary planform and open or closed slender bodies of non-circular contour may be analyzed. Longitudinal and lateral-directional static and rotary derivative solutions may be generated. The analysis was implemented on a time sharing system in conjunction with an input tablet digitizer and an interactive graphics input/output display and editing terminal to maximize its responsiveness to the preliminary analysis problem. Nominal case computation time of 45 CPU seconds on the CDC 175 for a 200 panel simulation indicates the program provides an efficient analysis for systematically performing various aerodynamic configuration tradeoff and evaluation studies.

Bonner, E.↗
DOE OSTI2023

Role of reinforcement learning for risk-based robust control of cyber-physical energy systems

Critical infrastructures such as cyber-physical energy systems (CPS-E) integrate information flow and physical operations that are vulnerable to natural and targeted failures. Safe, secure, and reliable operation and control of CPS-E is critical to ensure societal well-being and economic prosperity. Automated control is key for real-time operations and may be mathematically cast as a sequential decision-making problem under uncertainty. Emergence of data-driven techniques for decision making under uncertainty, such as reinforcement learning (RL), have led to promising advances for addressing sequential decision-making problems for risk-based robust CPS-E control. However, existing research challenges include understanding the applicability of RL methods across diverse CPS-E applications, addressing the effect of risk preferences across multiple RL methods, and development of open-source domain-aware simulation environments for RL experimentation within a CPS-E context. This article systematically analyzes the applicability of four types of RL methods (model-free, model-based, hybrid model-free and model-based, and hierarchical) for risk-based robust CPS-E control. Problem features and solution stability for the RL methods are also discussed. We demonstrate and compare the performance of multiple RL methods under different risk specifications (risk-averse, risk-neutral, and risk-seeking) through the development and application of an open-source simulation environment. Motivating numerical simulation examples include representative single-zone and multizone building control use cases. Finally, six key insights for future research and broader adoption of RL methods are identified, with specific emphasis on problem features, algorithmic explainability, and solution stability.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1986

Transonic flow solutions using a composite velocity procedure for potential, Euler and RNS equations

Solutions for transonic viscous and inviscid flows using a composite velocity procedure are presented. The velocity components of the compressible flow equations are written in terms of a multiplicative composite consisting of a viscous or rotational velocity and an inviscid, irrotational, potential-like function. This provides for an efficient solution procedure that is locally representative of both asymptotic inviscid and boundary layer theories. A modified conservative form of the axial momentum equation that is required to obtain rotational solutions in the inviscid region is presented and a combined conservation/nonconservation form is applied for evaluation of the reduced Navier-Stokes (RNS), Euler and potential equations. A variety of results is presented and the effects of the approximations on entropy production, shock capturing, and viscous interaction are discussed.

Gordnier, R. E.↗
NASA NTRS1989

Transonic flow solutions using a composite velocity procedure for potential, Euler and RNS equations

Solutions for transonic viscous and inviscid flows using a composite velocity procedure are presented. The velocity components of the compressible flow equations are written in terms of a multiplicative composite consisting of a viscous or rotational velocity and an inviscid, irrotational, potential-like function. This provides for an efficient solution procedure that is locally representative of both asymptotic inviscid and boundary layer theories. A modified conservative form of the axial momentum equation that is required to obtain rotational solutions in the inviscid region is presented and a combined conservation/nonconservation form is applied for evaluation of the reduced Navier-Stokes (RNS), Euler and potential equations. A variety of results is presented and the effects of the approximations on entropy production, shock capturing, and viscous interaction are discussed.

Gordnier, R. E.↗
NASA NTRS1995

An analytic radiative transfer model for a coupled atmosphere and leaf canopy

A new analytical radiative transfer model of a leaf canopy is developed that approximates multiple-scattering radiance by a four-stream formulation. The canopy model is coupled to a homogeneous atmospheric model as well as a non-Lambertian lower boundary soil surface. The same four-stream formulation is also used for the calculation of multiple scattering in the atmosphere. Comparisons of radiance derived from the four-stream model with those calculated by an iterative numerical solution of the radiative transfer equation show that the analytic model has a very high accuracy, even with a turbid atmosphere and a very dense canopy in which multiple scattering dominates. Because the coupling of canopy and atmospheric models fully accommodates anisotropic surface reflectance and atmospheric scattering and its effect on directional radiance, the model is especially useful for application to directional radiance and measurements obtained by remote sensing. Retrieval of biophysical parameters using this model is under investigation.

Liang, Shunlin↗
NASA NTRS1993

Unsteady panel method for flows with multiple bodies moving along various paths

A potential flow based three-dimensional panel method was modified to treat time dependent conditions in which several submerged bodies can move within the fluid along different trajectories. This modification was accomplished by formulating the momentary solution in an inertial frame-of-reference, attached to the undisturbed stationary fluid. Consequently, the numerical interpretation of the multiple-body, solid-surface boundary condition and the viscous wake rollup was considerably simplified. The unsteady capability of this code was validated by comparing computed and experimental results for a finite wing undergoing pitch oscillations. In order to demonstrate the multicomponent capability, computations were made for two wings following closely intersecting paths (e.g., to avoid mid air collisions) and for a flow field with relative rotation (e.g., helicopter-rotor/fuselage interaction). Results were compared to experimental data when such data was available.

Richason, Thomas F.↗
NASA NTRSDate not supplied

Processes and Methods used for NASA Communications Services Project Conceptualization and Formulation

The U.S. National Aeronautics and Space Administration (NASA) is partnering with the commercial SATCOM industry to demonstrate the industry’s capabilities to provide services to NASA’s near-Earth missions. Leveraging commercial capabilities will become increasingly necessary for NASA, as its incumbent Government-owned-and-operated near-Earth satellite network, the Tracking and Data Relay Satellite System (TDRSS), will approach retirement towards the end of the current decade. NASA is actively taking the first steps to making the vision of using commercial SATCOM a reality. For more than two years, NASA’s Communication Services Project (CSP) has been evaluating the feasibility of employing commercial satellite communication (SATCOM) networks for near-Earth operations. This past April, CSP awarded US$278.5 million in Funded Space Act Agreements to select commercial SATCOM providers, so they may begin developing and demonstrating their specific near-Earth space communication services that may support future agency missions. Each company will match or exceed agency contributions during the five-year development and demonstration period, totaling more than $1.5 billion of cost-share investment. The companies will complete technology development and in-space demonstrations by 2025, after which NASA intends to seek multiple long-term contracts to acquire these kinds of services. The six companies (Inmarsat Government Inc., Kuiper Government Solutions (KGS), SES Government Solutions, Space Exploration Technologies (SpaceX), Telesat U.S. Services and Viasat Inc.) will shortly begin developing and demonstrating their capabilities to provide services, ranging from launch to standard on-orbit operations. The demonstrations will be across multiple spectrum regimes from L-Band to optical, and will utilize Geostationary Orbit (GEO), Medium-Earth Orbit (MEO) and Low-Earth Orbit (LEO) constellations. This paper discusses methodologies and processes used over the past two years to reach the point of entering into these partnerships. It also includes CSP’s goals and objectives going forward, following the award of six Funded Space Act Agreements to commercial SATCOM providers to develop and demonstrate these capabilities.

Communications Services Project↗
NASA NTRSDate not supplied

Processes and Methods used for Nasa Communications Services Project Conceptualization and Formulation

The U.S. National Aeronautics and Space Administration (NASA) is partnering with the commercial SATCOM industry to demonstrate the industry’s capabilities to provide services to NASA’s near-Earth missions. Leveraging commercial capabilities will become increasingly necessary for NASA, as its incumbent Government-owned-and-operated near-Earth satellite network, the Tracking and Data Relay Satellite System (TDRSS), will approach retirement towards the end of the current decade. NASA is actively taking the first steps to making the vision of using commercial SATCOM a reality. For more than two years, NASA’s Communication Services Project (CSP) has been evaluating the feasibility of employing commercial satellite communication (SATCOM) networks for near-Earth operations. This past April, CSP awarded US$278.5 million in Funded Space Act Agreements to select commercial SATCOM providers, so they may begin developing and demonstrating their specific near-Earth space communication services that may support future agency missions. Each company will match or exceed agency contributions during the five-year development and demonstration period, totaling more than $1.5 billion of cost-share investment. The companies will complete technology development and in-space demonstrations by 2025, after which NASA intends to seek multiple long-term contracts to acquire these kinds of services. The six companies (Inmarsat Government Inc., Kuiper Government Solutions (KGS), SES Government Solutions, Space Exploration Technologies (SpaceX), Telesat U.S. Services and Viasat Inc.) will shortly begin developing and demonstrating their capabilities to provide services, ranging from launch to standard on-orbit operations. The demonstrations will be across multiple spectrum regimes from L-Band to optical, and will utilize Geostationary Orbit (GEO), Medium-Earth Orbit (MEO) and Low-Earth Orbit (LEO) constellations. This paper discusses methodologies and processes used over the past two years to reach the point of entering into these partnerships. It also includes CSP’s goals and objectives going forward, following the award of six Funded Space Act Agreements to commercial SATCOM providers to develop and demonstrate these capabilities.

Communications Services Project↗
DOE OSTI2025

Integration of multiple coinflip devices for high-quality random sampling

Artificial intelligence, scientific computing, and probabilistic computing use random sampling to approximate solutions to various problems, with larger models requiring a substantial quantity of random numbers. To generate the required vast quantity of random numbers at high rates, we explore so-called “coinflip” devices, which are stochastic microelectronic devices ideally capable of independently generating random bits with a tunable weight at a high rate. However, coinflip devices are inherently analog and demonstrate nonidealities, like temperature dependence and drift, that can introduce determinism into the outputs. We present important considerations for building systems of multiple coinflip devices to produce high-quality bitstreams with low error and little dependency on previous bits. Using tunnel diodes as coinflip devices, we implement a control loop to adapt to temperature dependence and generate fair bitstreams with each device. While this can lead to dependencies between bits in a single bitstream, we demonstrate that combining results generated in parallel with individual tunnel diodes can produce fair and unpredictable bitstreams. The suitability of these bitstreams for use in probabilistic computing is then demonstrated through a Monte Carlo approximation of π.

Taylor, Brady Garland [Sandia National Laboratorie↗
NASA NTRS2019

Cloud Computing Methods for Near Rectilinear Halo Orbit Trajectory Design

Complicated mission design problems require innovative computational solutions. As spacecraft depart from a proposed Gateway in a Near Rectilinear Halo Orbit (NRHO), recontact analysis is required to avoid risk of collision and ensure safe operations. Escape dynamics from NRHOs are governed by multiple gravitational bodies, yielding a trajectory design space that is exhaustively large. This paper summarizes the recontact analysis for departure from the NRHO and describes how the Deep Space Trajectory Explorer (DSTE) trajectory design software incorporates high performance cloud computing to compute and visualize the orbit design space. Recent focus on exploration missions to cislunar space has kindled accelerated interest in multibody orbit solutions. Trajectory analysis in the presence of multiple gravity fields is complex, and innovative computational tools are needed to simplify complicated design spaces, to generate large quantities of data quickly, and to visualize the output for user accessibility. The Gateway mission is a prime example. The Gateway1 is proposed as a human outpost in deep space. The current baseline orbit for the Gateway is a Near Rectilinear Halo Orbit (NRHO) near the Moon.2 The NRHO exists in a regime that experiences the gravitational effects of the Earth and the Moon simultaneously, complicating orbit analysis. The mission design process benefits greatly from updated computational tools for multibody missions like the Gateway. As an example, consider the problem of assessing the risk of collision in an NRHO. As a staging location to missions to the lunar surface and beyond the Earth-Moon system, the Gateway will experience spacecraft and other objects regularly arriving and departing. Departing objects potentially include spent logistics modules, visiting crew vehicles, debris objects, wastewater particles, and cubesats. Each departure is governed by the dynamics of the Gateway orbit and the surrounding dynamical environment. Over time, any unmaintained object in such an orbit eventually departs due to the small instabilities associated with the NRHOs. A separation maneuver speeds the departure from the NRHO, but the effects of the maneuver on the spacecraft behavior depend on the location, magnitude, and direction of the burn. Escape dynamics from the NRHO with regard to these maneuver options open up an enormous potential trajectory design space where subtle changes in input can produce dramatically large changes in the results. Any departing object must avoid recontacting the Gateway as it leaves the lunar vicinity, and a recontact analysis thus involves a significant number of computations and extensive output data. To explore the dynamics of this extensive design space, the Deep Space Trajectory Explorer3 (DSTE) trajectory design software incorporates new High Performance Computing (HPC) services and novel interactive visualizations. This paper details the HPC and cloud infrastructure techniques that are implemented in the DSTE, applying the new capabilities to analysis of recontact risk with the Gateway in NRHO. NEAR RECTILINEAR HALO ORBITS The Gateway is planned to fly in a lunar NRHO as its baseline orbit. The NRHO families of orbits are subsets of the larger halo families, which originate from planar orbits near the L1 and L2 libration points; the Earth-Moon L2 halo family appears in Figure 1. Each halo orbit is perfectly periodic in the Circular Restricted 3-Body Problem (CR3BP) and becomes a quasi-periodic orbit in a higher fidelity ephemeris force model. The NRHOs are defined as those members of the halo family with bounded stability properties;2 they pass near the Moon at perilune and are nearly polar. Families exist with apolunes located both above the lunar north pole and above the lunar south pole; the Gateway is planned to reside in a southern L2 NRHO in a 9:2 resonance with the lunar synodic period. The 9:2 NRHO is characterized by a period of about 6.5 days, a perilune radius of about 3,500 km, and an apolune radius of about 71,000 km; it is strongly affected by the gravity of both the Earth and the Moon simultaneously. This NRHO offers extended communications with assets on the south pole of the Moon,4 as well as low-cost orbit maintenance and attitude control,5 favorable eclipse avoidance properties,6 and inexpensive transfers from Earth and to other destinations.5,7 The NRHO portion of the southern L2 halo family is highlighted in black in Figure 1, and the 9:2 NRHO appears in blue.

Phillips, Sean M.↗
NASA NTRS1979

A two-dimensional unsteady Euler-equation solver for flow regions with arbitrary boundaries

A new technique is described for solving supersonic fluid dynamic problems containing multiple regions of continuous flow, each bounded by a permeable or impermeable surface. Region boundaries are, in general, arbitrarily shaped and time dependent. Discretization of such a region for solution by conventional finite difference procedures is accomplished using an elliptic solver which alleviates the dependence on a particular base coordinate system. Multiple regions are coupled together through the boundary conditions. The technique has been applied to a variety of problems including the shock diffraction and pointed wedge with detached bow shock.

Hindman, R. G.↗
NASA NTRS1981

Two-Dimensional Unsteady Euler-Equation Solver for Arbitrarily Shaped Flow Regions

A new technique is described for solving supersonic fluid dynamics problems containing multiple regions of continuous flow, each bounded by a permeable or impermeable surface. Region boundaries are, in general, arbitrarily shaped and time dependent. Discretization of such a region for solution by conventional finite difference procedures is accomplished using an elliptic solver which alleviates the dependence on a particular base coordinate system. Multiple regions are coupled together through the boundary conditions. The technique has been applied to a variety of problems including a shock diffraction problem and supersonic flow over a pointed ogive.

Hindman, R. G.↗
DOE OSTI2020

Low-Temperature Plasma Treatment for Enhanced Recovery of Highly Valued Critical REEs from Coal (Final Report)

The principle objective of this project was to develop a novel process flowsheet using low-temperature plasma treatment (LTP) integrated with conventional hydrometallurgical processes to recover rare earth elements (REEs), especially highly valued REEs (i.e. scandium and critical REEs), from coal and coal byproducts. LTP treatment was used to pretreat the solid feed material prior to leaching in an effort to improve the leaching characteristics of the REEs and potentially reduce acid consumption. Initially, a detailed experimental program was performed on a laboratory LTP furnace to identify a set of process parameter values to be used for assessing the enhanced leaching performance achieved on four specific gravity fractions of West Kentucky No. 13 and Fire Clay coarse coal refuse materials. High-temperature oxidation (HTO)) was also evaluated as a means for improving REE leachability on the same fractions of both coal sources. The lixiviants used to assess leachability improvements included: i) untreated water, ii) 0.1 M ammonium sulfate solution and iii) 1.2 M sulfuric acid solution at a temperature of 75°C. The laboratory data was used to develop multiple process flow sheets using HTO only to produce a rare earth oxide mix concentrate or utilizing both LTP and HTO treatment processes to generate a high purity heavy REO mix concentrate and a light REO mix concentrate. A techno-economic analysis was performed on each circuit. Low-temperature plasma oxidation was found to enhance REE leaching characteristics for all specific gravity (SG) fractions of the West Kentucky No. 13 coarse refuse material. The leaching recovery improvement was preferentially achieved for the heavy REEs with higher values obtained when treating the lower density fraction material. For example, the leaching recovery for the heavy REEs in the untreated 1.6 SG Float fraction material was 8% in a 0.1 M ammonium sulfate solution. After LTP treatment for 32 hours, heavy REE recovery increased to 33%. The corresponding recovery value for the light REEs was 10%. These findings indicate a strong affinity of the heavy REEs with the organic material, which is the component being reduced in content during the LTP treatment. Scanning electron micrographs showed the development of a honeycomb structure with increased LTP treatment with higher surface area as confirmed from BET measurements. Higher surface area promotes improved leaching efficiency due to enhance chemical interactions. When using 1.2 M sulfuric acid, the LTP treatment was found to provide a heavy REE recovery increase from 23% to 53% for the 1.6 SG Float fraction which corresponded to an ash content increase from around 20% to 35%. However, LTP treatment provided no improvement in leaching characteristic for the 2.2 SG sink material in the acid solution, which was likely due to the low organic content. The LTP treatment of the Fire Clay SG fractions provided no improvement in REE leaching characteristic, which may be due to the association with REE minerals such as monazite and REE-associated mineral matter such as clays. High temperature oxidation using temperatures in the range of 600-750OC was found to significantly enhance the leaching characteristics of both the light and heavy REEs for all density fractions in both coarse coal refuse materials. Generally, the increase in light REE recovery for the West Kentucky material was from 15% for the untreated material to 90% for the HTO treated material at 600°C. Recovery values dropped significantly when using temperatures below 600°C and above 750°C. The impact of the excessive temperatures was most likely a result of sintering. For the heavy REEs, HTO significantly improved recovery in the lighter SG fractions, which may be due to an affinity with the organic content. Only limited heavy REE recovery improvement was obtained for the denser SG fractions. Another benefit of HTO is the conversion of calcium-based minerals to a very soluble calcium oxide, which is removable using a mild acid solution, and the conversion of Fe minerals to relatively insoluble minerals such as hematite and magnetite. The pregnant leaching solution produced from the leaching of all SG fractions obtained from both coal sources was treated using two stages of solvent extraction procedurally modeling a rougher-cleaner circuit followed by oxalic acid precipitation. Rare earth oxide products having a purity greater than 97% was produced from nearly all of the SG fractions. A novel process flowsheet was developed based on the selective leaching of heavy REEs that was made possible through LTP treatment of West Kentucky No. 13 coarse refuse followed by salt leaching. In the flowsheet, coal refuse material is ground to a top size of 1mm and then treated by LTP. The treated material is mixed into a 0.1 M ammonium sulfate solution and then filtered to produce a heavy REE enriched pregnant leach solution. The heavy REEs are concentrated by two stages of selective precipitation to produce a heavy REE mix oxalate and then roasted to obtain a heavy REO mix concentrate. The solid residues remaining after salt leaching is dried and then subjected to HTO treatment followed by concentration using selective precipitation and roasting to produce a light REO mix concentrate. An alternative circuit was developed using HTO treatment only to treat 100% of the coarse refuse material followed by selective precipitation and roasting to generate a high purity REO mix product. A techno-economic analysis performed for the treatment of West Kentucky No. 13 and Fire Clay coarse refuse based on a hypothetical 500-tph plant revealed that the net present value for all scenarios was significantly negative. Primary reasons include: i) a substantial decrease in the basket value of the final product relative to the feed and ii) high chemical costs relative to the amount of product generated. The basket value reduction was primarily due to little to no scandium recovery due to poor leaching characteristics when using low acid concentrations as necessitated by the low feed grades. Operating costs were primarily influenced by the chemical costs relative to the kilograms of concentrate produced and feed grade. For the Fire Clay coarse refuse material, the total operating cost was estimated to be $354 per kg of REO produced which was higher than the basket value of the REO product ($276/kg). Approximately 37% of the total cost was the sulfuric acid used during the leaching process. A comparative analysis performed to assess the economic viability of using low-temperature plasma furnaces found that market value enhancement for the generation of both heavy REO and light REO concentrates will be needed in addition to significant technology advancements to reduce the cost of plasma furnaces.

01 COAL, LIGNITE, AND PEAT↗
NASA NTRS2013

A Floating Node Method for the Modelling of Discontinuities Within a Finite Element

This paper focuses on the accurate numerical representation of complex networks of evolving discontinuities in solids, with particular emphasis on cracks. The limitation of the standard finite element method (FEM) in approximating discontinuous solutions has motivated the development of re-meshing, smeared crack models, the eXtended Finite Element Method (XFEM) and the Phantom Node Method (PNM). We propose a new method which has some similarities to the PNM, but crucially: (i) does not introduce an error on the crack geometry when mapping to natural coordinates; (ii) does not require numerical integration over only part of a domain; (iii) can incorporate weak discontinuities and cohesive cracks more readily; (iv) is ideally suited for the representation of multiple and complex networks of (weak, strong and cohesive) discontinuities; (v) leads to the same solution as a finite element mesh where the discontinuity is represented explicitly; and (vi) is conceptually simpler than the PNM.

Pinho, Silvestre T.↗
NASA NTRS1987

Thermosolutal convection during directional solidification. II - Flow transitions

The influence of thermosolutal convection on solute segregation in crystals grown by vertical directional solidification of binary metallic alloys or semiconductors is studied. Finite differences are used in a two-dimensional time-dependent model which assumes a planar crystal-melt interface to obtain numerical results. It is assumed that the configuration is periodic in the horizontal direction. Consideration is given to the possibility of multiple flow states sharing the same period. The results are represented in bifurcation diagrams of the nonlinear states associated with the critical points of linear theory. Variations of the solutal Rayleigh number can lead to the occurrence of multiple steady states, time-periodic states, and quasi-periodic states. This case is compared to that of thermosolutal convection with linear vertical gradients and stress-free boundaries.

Mcfadden, G. B.↗
DOE OSTI2020

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↗
DOE OSTI2022

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↗