Search NASA⌕ Search

SEARCH · Search NASA

Results for “Homogenization”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Characterizing Turbulence at a Forest Edge: Comparing Sub-Filter Scale Turbulence Models in Simulations of Flow over a Canopy

In wildfires, atmospheric turbulence plays a major role in the transfer of turbulent kinetic energy. Understanding how turbulence feeds back into a dynamical system is important, down to the varying small scales of fuel structures (i.e. pine needles, grass). Large eddy simulations (LES) are a common way of numerically representing turbulence. The Smagorinsky model (1963) serves as one of the most studied sub-grid scale representations in LES. In this investigation, the Smagorinsky model was implemented in HIGRAD/FIRETEC, LANL’s coupled fire-atmosphere model. This study was motivated by the need to quantitatively investigate the vorticity budget equation in HIGRAD/FIRETEC. The Smagorinsky turbulent kinetic energy (TKE) was compared to FIRETEC’s 1.5-order TKE eddy-viscosity subgrid-scale model, known as the Linn turbulence model. This was done in simulations of flow over flat terrain with a homogeneous, cuboidal canopy in the center of the domain. Examinations of the modeled vertical TKE profile and turbulent statistics at the leading edge, and throughout the canopy, show that the Smagorinsky model provides comparable results to that of the original closure model posed in FIRETEC.

58 GEOSCIENCES↗
DOE OSTI2024

Evidence of Superconductivity in Electrical Resistance Measurements of Hydrides Under High Pressure

In the standard van der Pauw four-probe configuration commonly used for electrical resistance measurements, including those of hydrogen-rich samples at high pressures, the application of electrical current through one pair of leads while measuring voltage difference across another pair is a fundamental practice (see the inset in Fig. 1c). In the scenario described in Ref. [5], where electrical resistance vanishes due to disconnection between current and voltage probes or due to an onset of giant magnetoresistance in parts of the sample, it is crucial to consider the implications across different probe orientations: while one orientation may indeed result in vanishing electrical resistance, a divergence of the electrical resistance towards infinity will be observed in an alternative orientation where the applied current pattern is rotated a quarter turn (see the scheme on the inset in Fig. 1c). It is imperative to acknowledge that experimental data encompassing all possible probe orientations are essential for eliminating artifacts and assessing sample homogeneity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Fractionation of UF6 and daughter progeny in storage cylinders from external heating

Abstract Uranium hexafluoride (UF 6 ) is a significant concern for material accountancy and verification in the international safeguards community. Verification of the contents of UF 6 cylinders is generally attempted with gamma spectroscopy but the current methods assume a uniform, homogeneous UF 6 mass distribution within the cylinder. In this work, it was found experimentally and confirmed via modeling, that under an external heat load (the sun), the UF 6 and its daughter products undergo fractionation in the cylinder. This fractionation of the UF 6 and daughter products can cause an errant measurement of the enrichment of the cylinder when using the current verification methods.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2025

Polymer-assisted deposition of nanoparticle feedstocks for target fabrication

Traditional polymer-assisted deposition has been shown to produce highly uniform thin films of metal oxides, including actinide oxides. Furthermore, while producing thicker films for nuclear targets is possible through repeated coating application, we exchanged the dissolved metal species with nanoparticles to maximize the thickness that can be achieved with an individual layer. Using CeO 2 nanoparticles in a polyethyleneimine matrix, we produced targets with single-layer areal densities of 0.22 ± 0.01 mg·cm −2 (1σ) and thicknesses of 690 ± 80 nm (1σ). A custom 3D-printed spin coating chuck attachment with an inlay improved target homogeneity and will streamline future work with radioactive materials.

and nuclear chemistry↗
DOE OSTI2025

Uncertainty quantification for nuclear forensics with population analyses

Although neural networks offer cutting-edge predictive power, their deployment in high-consequence nuclear forensic applications is limited, partly because of their black-box nature. Incorporating robust uncertainty quantification methods into the predictive frameworks of neural networks is progress towards their future deployment in such scenarios. This work integrates uncertainty quantification into neural networks for nuclear reactor core-average burnup estimation from simulated environmental samples. We test two regimes (homogeneous and heterogeneous events) on DeepSets and Set Transformer architectures, we find both quantify predictive uncertainty effectively, but Set Transformer excels in partitioning latent events, offering superior predictive power and more informative uncertainty estimates.

Hatton, Conner [ORNL] (ORCID:0009000804970959)↗
DOE OSTI2025

Initial Conditions of Planet Formation: Time Constraints from Small Bodies and the Lifetime of Reservoirs in the Solar Protoplanetary Disk

This review explores the timescales of the initial phase of planet formation, from nebular dust (CAIs and chondrules) to planetesimal accretion and differentiation, using evidence from meteorite research. Aluminium-Mg systematics of CAIs indicate either an extended period of CAI formation (∼0.3 Ma) or an initial 26 Al heterogeneity, with evidence supporting a homogeneous 26 Al abundance in the protoplanetary disk. Thermal and aqueous alteration on the parent body can disturb the U-Pb and Al-Mg chronometers in chondrules. Focusing on relatively robust isochron data from plagioclase of pristine (types ≤3.05) chondrites indicates a shift in chondrule formation locations, moving from the inner to the outer disk over time. Ages of basaltic achondrites show that silicate differentiation on small bodies was well underway within the first few million years (Ma) of our solar system. Their age record, however, reveals inconsistencies between different chronometers, partly caused by (i) secondary disturbances due to thermal metamorphism, aqueous alteration, or impacts, (ii) the presence of xenolithic minerals, and (iii) potentially variable initial 26 Al abundances due to disturbances at the mineral scale. Nucleosynthetic isotope data indicate that parent bodies of iron and stony meteorites formed in two distinct regions within the protoplanetary disk: the inner, non-carbonaceous (NC) and the outer, carbonaceous (CC) region. Based on Hf-W chronometry it has been demonstrated that NC and CC parent bodies of magmatic iron meteorites segregated their cores within ∼1–3 Ma after CAI formation, implying that parent body accretion occurred within <1 Ma in both reservoirs. Combining accretion ages with nucleosynthetic data further reveals that, at first order, NC and CC reservoirs in the solar protoplanetary disk were established within 1 Ma and existed over several Ma with limited exchange between them. In the CR chondrite accretion region of the disk, planetary bodies formed over at least 3 Ma, while in most other regions, formation spanned at least 1 Ma, with minimal changes in nucleosynthetic isotope compositions. Aerodynamical size sorting of dust likely introduced or amplified some of these variations.

26Al distribution↗
DOE OSTI2025

Quantifying precursors to void nucleation and coalescence in aluminum

Ductile rupture is a common failure mode for engineering alloys. It is generally comprised of three mechanisms: void nucleation at secondary particles, void growth, and coalescence of voids. Conventional models for ductile rupture have limited precision for predicting macroscopic failure. This is partially because they have been corroborated using ex-situ observations of these three mechanisms and calibrated by macroscopic strain or stress metrics. In this study, in-situ high-energy X-ray characterization was conducted to correlate sites of void nucleation via particle cracking and sites of void growth with grain-scale metrics. Here it is shown that particle cracking is not predicated by elevated stress metrics, which disagrees with classical models. Instead, particle cracking tended to occur for the largest, least spherical particles. This trend persisted on the aggregate-scale and within individual neighborhoods around particles. Furthermore, an Eshelby analysis illuminated a statistically significant increase in maximum principal stress within a cracked particle, compared to an uncracked particle. In-situ observations also revealed two separate examples of intragranular void growth comprised of flat, crack-like features that grew showing alignment with slip systems of the highest resolved shear stress. Post-mortem fractography revealed a secondary population of dimples on these crack-like features, implying that void sheeting may be occurring during this crack-like void growth. Furthermore, these in-situ observations imply that classical models for void growth, that assume a homogenous medium, should be extended to account for the anisotropic grain-level behavior to correctly capture distinct features of void nucleation and growth.

Al-2219↗
DOE OSTI2026

Quantitative characterization of gradient microstructures: A study on friction stir spot processing of pure cobalt

Heterogeneous microstructures in polycrystalline metals can enhance the strength and ductility, outperforming homogeneous structures of similar composition. This study investigates deformed cobalt via friction stir spot processing with varying dwell times to uncover the effects of plastic deformation and heat generation on the formation of morphological, phase, and grain boundary character gradients. A new approach to quantify the morphological gradients in materials, which describes grain morphology in terms of density followed by parametric regression, enables direct quantification of processing depth and gradient sharpness. Results show that longer processing times increase the steepness of morphological gradients and reduce the deformation depth for friction stir spot processing with low plunge depths and high tool rotational speeds. The amount of retained FCC is increased in the shorter processing conditions, primarily due to refined grain size, increased defect content, and reduced heat generation. Crystallographic texture analysis of the HCP phase indicated a dominant B-fiber described by (0001) ∥ shear plane normal in the extreme processing conditions and the formation of a P-fiber, shear direction ∥ ⟨11$\bar2$0⟩ for intermediate dwell times. The texture of the FCC phase for low processing times was a C texture {100}⟨011⟩ where longer processing times were dominated by a [001] fiber texture with a main {110}⟨100⟩ orientation and emergence of a slight [111] fiber in the longest processing condition. The approaches outlined in this work give insight into quantifying gradients and improve the understanding of highly deformed cobalt.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Measuring thermomechanical response of large-format printed polymer composite structures via digital image correlation

Large-format additive manufacturing (LFAM) is a branch of additive manufacturing (AM) research with the ability to create large structures typically measuring several meters in scale. LFAM is advantageous for tooling applications, not only because it offers the ability to create complex geometries not easily made using subtractive manufacturing processes, but the cost savings of pelletized feedstock used by these systems result in larger parts printed at faster speeds than traditional AM systems. Fiber reinforced polymer (FRP) is a commonly used feedstock material in LFAM structures because it reduces the distortion experienced during printing. However, FRP introduces highly anisotropic thermomechanical properties and contributes to a nonhomogeneous microstructure that can result in critical distortion of dimensions during tooling. Measuring the global thermomechanical response of LFAM structures requires a more representative method that accounts for not only anisotropic properties but also the nonhomogeneous nature of the final part. This is where traditional techniques to measure thermomechanical response, such as thermomechanical analysis (TMA), fall short as they assume homogeneity. This study evaluated the coefficient of thermal expansion (CTE) of LFAM structures as measured by TMA as compared to a novel digital image correlation oven (DIC Oven) system. The LFAM structures were made from 20 % by weight carbon fiber reinforced acrylonitrile butadiene styrene (CF-ABS). TMA measurements showed significant variations in CTE across a single LFAM bead, confirming the need for a global technique that captures overall thermomechanical response. The CTE values measured using the DIC Oven compared well to average TMA values obtained from localized measurements across the sample. The DIC Oven was also used to quantify the effects of different layer orientations on thermomechanical properties, which cannot be easily captured using TMA. A predictive model was also developed by using localized TMA values across an LFAM bead to predict the overall thermomechanical response of an LFAM structure.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dissolution zone model of the oxide structure in additively manufactured dispersion-strengthened alloys

The structural evolution of oxides in dispersion-strengthened superalloys during laser-powder bed fusion is considered in detail. Alloy chemistry and process parameter effects on oxide structure are assessed through a parameter study on the model alloy Ni-20Cr, doped with varying concentrations of Y 2 O 3 and Al. Small angle neutron scattering measurements of the dispersoid size distribution show the dispersoid size increases with higher laser power, slower scan speed, and increasing Y 2 O 3 and Al content. Complementary electron microscopy measurements reveal reactions between Y 2 O 3 and Al, even in nanoscale dispersoids, and the presence of micron-scale oxide slag inclusions in select specimens. A scaling analysis of mass and momentum transport within the melt pool, presented here, establishes that diffusional structural evolution mechanisms dominate for nanoscale dispersoids, while fluid forces and advection become significant for larger slag inclusions. These findings are developed into a theory of dispersoid structural evolution, integrating quantitative models of diffusional processes – dispersoid dissolution, nucleation, growth, coarsening – with a reduced order model of time-temperature trajectories of fluid parcels within the melt pool. Calculations of the dispersoid size in single-pass melting reveal a zone in the center of the melt track in which the oxide feedstock fully dissolves. Within this zone the final Y 2 O 3 size is independent of feedstock size and determined by nucleation and growth kinetics. If the dissolution zones of adjacent melt tracks overlap sufficiently with each other to dissolve large oxides, formed during printing or present in the powder feedstock, then the dispersoid structure throughout the build volume is homogeneous and matches that from a single pass within the dissolution zone. Gaps between adjacent dissolution zones result in oxide accumulation into larger slag inclusions. Predictions of final dispersoid size and slag formation using this dissolution zone model match the present experimental data and explain process-structure linkages speculated in the open literature.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Infrared triggered dwell and active cooling thermal control effects on microstructural uniformity in DED

Directed Energy Deposition (DED) offers rapid large scale fabrication, but difficulty in delivering consistent microstructures and properties hinders the use of DED fabricated components in safety or performance critical applications. Variability stems from the complex thermal cycles generated by the toolpath used to print the required geometry. Several practical methods have become established in DED to regulate overheating, such as active cooling of the baseplate structure or the use of an infrared camera to inject interlayer pauses to ensure the top layer of the component cools to a set temperature, which have been shown to affect microstructure. However, no critical assessment has been performed as to how effective these controls are in promoting microstructural uniformity in the context of complex layer timing commonly generated by non-prismatic geometries. Here we show how controls influence the thermal field, phase transformations, and dynamic annealing of a low-temperature transformation steel using infrared imaging and operando neutron diffraction. Counterintuitively, common thermal homogenization process controls can reduce microstructural uniformity because these approaches stabilize peak temperature while overlooking temperatures near the solid-state phase transformation fronts. Instead, the cyclic reheating induces spatially-variant dynamically annealed regions which can be modulated via control parameters. We show that these controls have spatially linked effects centimeters away from the active weld, which implies that microstructure control must co-optimize thermal input across many subsequent layers. In conclusion, our results demonstrate the pressing need for higher order controls that integrate predictive elements of simulation data to stabilize printed properties for future qualification of DED components.

Directed energy deposition↗
DOE OSTI2026

A meshing framework for digital twins for extrusion based additive manufacturing

Additive manufacturing (AM) allows for manufacturing of complex three-dimensional geometries not typically realizable with standard manufacturing practices. The internal microstructure of AM components has a significant impact on mechanical, vibrational, and shock properties and permits richer design space when this is controllable. Due to complex interactions of internal geometry of an extrusion-based AM component, it is common practice to assume homogeneous behavior or to perform characterization testing on specific toolpath configurations. To avoid testing or material waste, it is necessary to develop a consistently accurate numerical simulation framework with relevant boundary value problems that can handle the complicated geometry of internal material microstructure present in AM components. Herein, a framework is proposed to directly create computational meshes suitable for finite element analysis (FEA) of the fine-scale features generated from extrusion-based AM tool paths to maintain a strong process–structure–property-performance linkage. This mesh can be manually or automatically analyzed using standard FEA simulations such as quasi-static preloading or modal analysis. The framework allows an in-silico assessment of a target AM geometry where fine-scale features greatly impact quantities of design interest such as in soft elastomeric lattices where toolpath infill can greatly influence the self-contact of a structure in compression, which we use as a motivating exemplar. This approach greatly reduces both time and resource waste present in traditional build and test design cycles for non-intuitive design spaces, and acts as a tool for use in the production of a key component of a digital twin, a mesh suitable for finite element analysis. In conclusion, it also further allows for the exploration of toolpath infill to optimize component properties beyond simple linear properties such as density and stiffness.

Additive manufacturing↗
DOE OSTI2024

A Modeling framework for flocculated cohesive sediment transport in the current bottom boundary layer

Cohesive sediment transport, where its settling velocity is controlled by the flocculation process, is a crucial component in determining biochemical cycles, fate of pollutants, and morphodynamics in many aquatic ecosystems. In this study, a modeling framework is presented to investigate how flocculation influences cohesive sediment transport in the current bottom boundary layer in dilute conditions, consistent with the calibration range of the flocculation model. From a local analysis of floc dynamics in homogenous turbulence, we identify that the floc size distribution is mainly controlled by floc cohesion and yield strength. The uncertainty in fractal dimension plays a minor role for the floc size but it influences the resulting floc density and settling velocity. The transport analysis in the current boundary layer shows that the flocculation process alters the vertical distribution of the settling velocity and hence the sediment concentration with a strong dependence on cohesion, floc yield strength, and floc structure. When the flocs are more susceptible to breaking, a well-mixed concentration profile is obtained. In contrast, for flocs with higher cohesion or yield strength, higher concentration with a sharp gradient is observed close to the bed. Overall, the settling velocity exhibits a low vertical variability within 20 % of the depth-averaged value except near the bed. Further, this suggests that using a depth-averaged settling velocity yields acceptable predictions of the sediment concentration profiles, especially for flocs with lower cohesion.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Noise-aware optimization in nominally identical manufacturing and measuring systems for high-throughput parallel workflows

Device-to-device variability in experimental noise critically impacts reproducibility, especially in automated, high-throughput systems like additive manufacturing farms. While manageable in small labs, such variability can escalate into serious risks at larger scales, such as architectural 3D printing, where noise may cause structural or economic failures. This contribution presents a noise-aware decision-making algorithm that quantifies and models device-specific noise profiles to manage variability adaptively. It uses distributional analysis and pairwise divergence metrics with clustering to choose between single-device and robust multi-device Bayesian optimization strategies. Unlike conventional methods that assume homogeneous devices or enforce generic robustness, the proposed framework explicitly determines whether shared optimization across devices is appropriate based on the degree of inter-device noise heterogeneity. This enables improved performance, reproducibility, and efficiency. An experimental case study involving three nominally identical 3D printers (same brand, model, and close serial numbers) demonstrates reduced redundancy, lower resource usage, and improved reliability, along with improved convergence stability and solution quality through the selection of the appropriate optimization strategy based on the degree of inter-device noise heterogeneity. Overall, this framework establishes a general approach for precision- and resource-aware optimization in scalable, automated experimental platforms, demonstrated here on a representative multi-device 3D printing case study.

Schenk, Christina↗
DOE OSTI2025

Spectral distortions to momentum and scalar exchanges by non-turbulent motion and patchy landscape variability

Modifications to the spectra of turbulent velocity and scalars and co-spectra of vertical fluxes of momentum and scalars due to patchy landscape heterogeneity and non-stationarity are explored for a Mediterranean ecosystem. About 9 months of high frequency measurements of the three velocity components, water vapor concentration, carbon dioxide concentration, and air temperature were analyzed for different seasons (spring/summer) and prevalent wind directions (southeast/northwest). The two wind directions sampled a contrast of clumped and patchy landscape comprised of olive trees (southeast) and wall bounded flow disturbed by the presence of few upwind trees (northwest). The measured spectra and co-spectra were also compared to theoretical scaling forms from stationary, planar homogeneous flow, in the absence of subsidence as derived from the Kansas experiment. To assess the role of low frequency non-turbulent motion on the spectral and co-spectral content, a 5-min Fourier cutoff was introduced and the analysis was limited to near-neutral conditions where the boundary layer depth is shallow compared to its unstable counterpart. It was shown that the velocity statistics were not appreciably impacted by the low-frequency motion causing non-stationarity. Moreover, the turbulent scalar fluxes were also shown not to be significantly impacted by such low frequency motion. The scalar variances were impacted, especially the water vapor variance and its concomitant spectral shape. When the non-turbulent motion was filtered, the scalar spectra at low wavenumbers followed expectations from the so-called attached eddy hypothesis (i.e. exhibited a $k^{-1}_x$ scaling with $k_x$ defining the longitudinal wavenumber) applicable for near-neutral conditions. For momentum co-spectra, the canonical shapes from the Kansas experiment appear to describe well the measurements here and in both dominant directions and seasons with some adjustment to the integral time scales based on wind direction. For the scalar co-spectra, deviations from the Kansas experiment were prevalent. The most noticeable and surprising deviations were their slow decay with increased sampling frequency at inertial subrange scales. This slow decay was shown not to contribute appreciably to the overall scalar fluxes. At those fine scales, predictions from local isotropy were expected to hold. The scalar co-spectral deviations from local isotropy were then discussed using a simplified co-spectral budget model where scalar–scalar co-spectra naturally emerged and the interplay between landscape heterogeneity and a scale-dependent pressure-scalar de-correlation time was postulated. It is also envisaged that the findings here offer a preliminary template for analyzing eddy-covariance data in situations that deviate from ideal conditions, especially regarding low-frequency modulations of scalar spectra and vertical scalar flux co-spectra.

Canopy turbulence↗
DOE OSTI2026

Impact of vertical and seasonal variation in leaf traits on simulating soybean canopy photosynthesis via 1D and 3D modeling

Accurate modeling of photosynthesis is crucial for predicting crop productivity and quantifying the carbon cycle in agroecosystems. Leaf traits are essential inputs for modeling canopy photosynthesis. Yet, many existing models still use fixed plant functional type (PTF)-based values to parameterize leaf traits under a big-leaf or two-big-leaf assumption, neglecting their vertical profiles and seasonal changes. This simplification may introduce significant uncertainties in estimating gross primary productivity (GPP). In this study, we simulated soybean GPP and tested the effects of vertical and seasonal variation in three key leaf photosynthetic traits: the maximum carboxylation rate at 25 °C (Vcmax 25 ), leaf chlorophyll content (LCC), and leaf mass per area (LMA) in the 1D-SCOPE and 3D-Helios models. Weekly field measurements were conducted during the growing season of 2024 to support the simulation. We designed ten leaf trait parameterization schemes by incorporating different combinations of vertical profiles and seasonal changes, while assuming homogeneous canopy architecture in both models. Our results revealed that Vcmax 25 vertical and seasonal variation had the strongest influence on simulated GPP in both 1D and 3D models, while LCC and LMA effects were minimal. Particularly, the scheme with an empirically parameterized Vcmax 25 profile achieved comparable performance to the scheme with the measured Vcmax 25 profile. Both 1D-SCOPE and 3D-Helios accurately modeled GPP (SCOPE: R 2 = 0.87, Bias = 0.55 µmol m⁻² s⁻¹; Helios: R 2 = 0.9, Bias = 0.22 µmol m⁻² s⁻¹) under the most complex scheme, and their responses to vertical and seasonal variation in leaf traits were consistent, demonstrating the robustness of our findings. Based on our findings, we propose a scalable framework for parameterizing leaf traits to improve GPP simulations. This study contributes to improving the representation of leaf trait dynamics in canopy-level photosynthesis models, potentially enhancing our ability to predict crop productivity and understand agroecosystem carbon dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Griffin: A MOOSE-based reactor physics application for multiphysics simulation of advanced nuclear reactors

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor physics application for multiphysics simulations of advanced reactor designs jointly developed by Idaho National Laboratory and Argonne National Laboratory. This paper summarizes the motivation, significance, architecture, design, and features of Griffin. Griffin offers flexible and extensible features to address the challenges associated with advanced reactor designs. These features range from fundamental particle transport to specific reactor physics tasks. The features cover a wide range including on-the-fly and traditional two-step cross-section generation methods, steady-state and transient transport solvers suitable for both heterogeneous and homogeneous models, high-fidelity depletion where thousands of isotopes can be tracked and low-fidelity depletion characterized by burnup, etc. The most fundamental aspect that sets Griffin apart from other reactor analysis codes is that it is developed based on the MOOSE framework. A modular development approach is strongly enforced, with multiphysics being an essential element considered since the beginning of Griffin’s development. Griffin links various MOOSE physics modules and couples to other MOOSE-based applications and non-MOOSE-based applications for multiphyiscs simulations. Griffin includes three modules: ISOXML for preparing and managing multigroup cross sections, radiation transport for solving the neutron transport equation, and reactor analysis for user-oriented reactor physics analysis functionalities. Griffin uses various finite element methods for spatial discretization, multigroup approximation for energy discretization and discrete ordinates method, spherical harmonics expansion method, and diffusion approximation for streaming direction discretization to solve the neutron transport equation. Griffin’s flexibility is evidenced through Griffin’s various applications to fast reactor, high-temperature reactor, pebble bed reactor, molten salt reactor, and microreactor designs. Griffin development follows the software quality assurance procedure for MOOSE-based applications and with software requirements consistent with the ASME NQA-1 standard. Griffin has been adopted into the reactor analysis system for the U.S. NRC and is in use at U.S. companies, universities and national laboratories.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Development of a Griffin model of the advanced test reactor

In the pursuit of a higher fidelity deterministic simulation capability of the Advanced Test Reactor, it is important to have a fast yet accurate deterministic neutronics model. Here, to achieve this, we employed an advanced two-step method. The first step involves generating homogenized cross sections using OpenMC, a cutting-edge Monte Carlo neutron transport code. OpenMC offers excellent modular capabilities, allowing for easy component integration and flexibility in incorporating new designs into the model. The second step involves deterministic transport calculations, which are performed using Griffin, a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). To ensure the accurate spatial resolution and assignment of material cross sections, a Cubit-generated mesh for the Advanced Test Reactor is utilized as an intermediate step between the OpenMC and Griffin models; Griffin utilizes the mesh for its finite element solution, while OpenMC material identifications are written to the mesh file to be used in Griffin material assignments. Additionally, a Python-based script converts the cross sections generated by OpenMC into the ISOXML format required by Griffin. Initial comparisons using the Griffin diffusion solver indicated good agreement between the neutron multiplication factors obtained from the standalone OpenMC model and the Griffin model, with differences of less than 10 pcm in the 2D geometry configuration; it was later determined that this agreement was likely due to compensating effect and was more likely on the order of –700 pcm relative to the OpenMC solution. However, in three-dimensional calculations, an unacceptably large error (almost 8,000 pcm) was found in the Griffin solution with the diffusion solver. Subsequent calculations using Griffin’s discrete ordinates solver demonstrated substantially improved agreement, within 116 pcm of the OpenMC solution used to generate the cross sections for Griffin. Building on this capability, future work will seek to perform more detailed validation calculations. The ultimate goal is to evaluate both transient and multiphysics simulations of the reactor.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗