Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple solutions”

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

Verification of the DIF3D Software to Support Fast Reactor Analysis (Rev. 3)

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.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

HIPPO – A Software Platform for Electricity Market Research and Development

The goal of this project is to provide Regional transmission organizations (RTOs) and independent system operators (ISOs) a market design and prototyping software, High-Performance Power-Grid Optimization (HIPPO), that they can evaluate electricity market design options, calculate market planning strategies and operational performance. With the high standards and strict reliability requirements for operating power systems, impacts of new technologies need to be fully investigated prior to any consideration for adoption. A market design and prototyping software tool which can be used to prototype electricity market design options, to calculate market planning strategies and operational performance with high precision, and to investigate the impacts for integrating future power grid technologies will be valuable to RTOs/ISOs who operate power systems, to vendors like GE and ABB who provide the market solvers, and to market participants and researchers who are actively doing market research. HIPPO is a such tool that can be used to improve the current market operations and provide capabilities for rigorous forward-looking design and prototyping of next-generation energy markets. HIPPO has a high-resolution model for the day-ahead SCUC, which was validated with MISO and GE-Grid Solutions. HIPPO is built with parallel and distributed computing capabilities and can be executed in both multi-thread and high-performance computing (HPC) settings. This capability provides fast solution speed necessary to handle the larger and more complex SCUC problems of real-world cases and the potentially growing size and complexity of future scenarios. In addition, HIPPO has a concurrent optimizer (CO) which manages multiple algorithm executions simultaneously and leverages the advantages from different algorithms. This structure provides flexibility to better benchmark competing approaches. Highly accurate market model, fast solution technologies and flexible model and algorithm control are the features which will make HIPPO an extensible platform for developing and testing multiple approaches to meet a wide range of future market needs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Path-Integrated Stereo X-Ray Digital Image Correlation: Resolving the Violation of Conservation of Intensity

In this study, x-ray imaging addresses many challenges with visible light imaging in extreme environments, where optical diagnostics such as digital image correlation (DIC) and particle image velocimetry (PIV) suffer biases from index of refraction changes and/or cannot penetrate visibly occluded objects. However, conservation of intensity—the fundamental principle of optical image correlation algorithms—may be violated if ancillary components in the X-ray path besides the surface or fluid of interest move during the test. The test series treated in this work sought to characterize the safe use of fiber-epoxy composites in aerospace and aviation industries during fire accident scenarios. Stereo X-ray DIC was employed to measure test unit deformation in an extreme thermal environment involving a visibly occluded test unit, incident surface heating to temperatures above 600°C, and flames and soot from combusting epoxy decomposition gasses. The objective of the current work is to evaluate two solutions to resolve the violation of conservation of intensity that resulted from both the test unit and the thermal chamber deforming during the test. The first solution recovered conservation of intensity by pre-processing the path-integrated X-ray images to isolate the DIC pattern of the test unit from the thermal chamber components. These images were then correlated with standard, optical DIC software. The second solution, called Path-Integrated Digital Image Correlation (PI-DIC), modified the fundamental matching criterion of DIC to account for multiple, independently-moving components contributing to the final image intensity. The PI-DIC algorithm was extended from a 2D framework to a stereo framework and implemented in a custom DIC software. Both solutions accurately measured the cylindrical shape of the undeformed test unit, recovering radii values of R = 76.20±0.12 mm compared to the theoretical radius of R theor = 76.20 mm. During the test, the test unit bulged asymmetrically as decomposition gasses pressurized the interior and eventually burned in a localized jet. Both solutions measured the heterogeneous radius of this bulge, which reached a maximum of approximately R = 91 mm, with a discrepancy of 2–3% between the two solutions. Two solutions that resolve the violation of conservation of intensity for path-integrated X-ray images were developed. Both were successfully employed in a stereo X-ray DIC configuration to measure deformation of an aluminum-skinned, fiber-epoxy composite test unit in a fire accident scenario.

47 OTHER INSTRUMENTATION↗
NASA NTRS2011

Cartesian Off-Body Grid Adaption for Viscous Time- Accurate Flow Simulation

An improved solution adaption capability has been implemented in the OVERFLOW overset grid CFD code. Building on the Cartesian off-body approach inherent in OVERFLOW and the original adaptive refinement method developed by Meakin, the new scheme provides for automated creation of multiple levels of finer Cartesian grids. Refinement can be based on the undivided second-difference of the flow solution variables, or on a specific flow quantity such as vorticity. Coupled with load-balancing and an inmemory solution interpolation procedure, the adaption process provides very good performance for time-accurate simulations on parallel compute platforms. A method of using refined, thin body-fitted grids combined with adaption in the off-body grids is presented, which maximizes the part of the domain subject to adaption. Two- and three-dimensional examples are used to illustrate the effectiveness and performance of the adaption scheme.

Buning, Pieter G.↗
DOE OSTI2022

Flow and transport in three-dimensional discrete fracture matrix models using mimetic finite difference on a conforming multi-dimensional mesh

Here, we present a comprehensive workflow to simulate single-phase flow and transport in fractured porous media using the discrete fracture matrix approach. The workflow has three primary parts: (1) a method for conforming mesh generation of and around a three-dimensional fracture network, (2) the discretization of the governing equations using a second-order mimetic finite difference method, and (3) implementation of numerical methods for high-performance computing environments. A method to create a conforming Delaunay tetrahedralization of the volume surrounding the fracture network, where the triangular cells of the fracture mesh are faces in the volume mesh, that addresses pathological cases which commonly arise and degrade mesh quality is also provided. Our open-source subsurface simulator uses a hierarchy of process kernels (one kernel per physical process) that allows for both strong and weak coupling of the fracture and matrix domains. We provide verification tests based on analytic solutions for flow and transport, as well as numerical convergence. We also provide multiple expositions of the method in complex fracture networks. In the first example, we demonstrate that the method is robust by considering two scenarios where the fracture network acts as a barrier to flow, as the primary pathway, or offers the same resistance as the surrounding matrix. In the second test, flow and transport through a three-dimensional stochastically generated network containing 257 fractures is presented.

97 MATHEMATICS AND COMPUTING↗