Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fuel elements”

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 451 records · Page 25

Exploring the Nature of f-Element Soft Donor Interactions Using Electronically Tunable Azolate Ionic Liquids

This project was undertaken to advance the understanding of how f-elements interact with moderately soft donors, a heavily investigated yet open question which is of prime importance in spent nuclear fuel processing and fundamental inorganic chemistry. During the course of the project, based on exciting results, a stretch goal was developed to study the hydrolysis products of transuranic actinide metals, a somewhat understudied field even with its significance in nuclear fuel cycle and impacts in environmental chemistry. The stretch goal was to take our serendipitous discovery of an easy route to isolation of crystalline multinuclear ƒ-element hydrolysis products, and apply it to gaining a mechanistic understanding of Pu(III/IV) colloid formation. The simplicity of our techniques should lend themselves to the remote handling required for study of many transuranic elements. We developed several methodologies using azolium azolate chemistry to overcome ƒ-element hydrolysis problems that prohibit the isolation of ƒ-element soft donor complexes and to build a crystallographic library of ƒ-element N-donor complexes as a means to understand the fundamental differences between actinide and lanthanide interactions with moderately soft donor ligands. Our next major endeavor will be to transfer this chemistry from 4ƒ elements to transuranic elements, particularly in the study of hydrolysis of Pu(III/IV). While our work is fundamental in nature, applications of the knowledge we are generating should be felt in such diverse fields as catalysis, separations in general, nuclear waste remediation specifically, and many other applications in f-element magnetic and luminescent properties. The potential ramifications of the consistent and robust formation of hydrolysis controlled hexanuclear lanthanide structures are enormous, with future uses being catalyst formation, higher-nuclearity structure synthesis using our hexanuclear motif as a fundamental building block, Pu waste remediation, separations, and many other potential applications resulting from characteristic magnetic and luminescent properties of lanthanide polynuclear structures. Three synthetic methodologies (direct mixing with variable stoichiometries, use of volatile solvent, metathesis) were developed starting with 7 acidic and 6 basic azoles to obtain ionic liquids suitable for f-element coordination. Proton transfer by acidic/basic azole combination led to suitable low melting salts and two cocrystals. Acid/base reaction of azoles with soft-donor permanent cations of ([X 4444 ][OH] (where [X 4444 ] + = tetrabutylammonium [N 4444 ] + or tetrabutylphosphonium [P 4444 ] + ) with weakly acidic azoles including imidazole, 1,2,3-triazole, 1,2,4-triazole, 5-aminotetrazole, 4,5-dicyanoimidazole, and 2-amino-4,5-dicyanoimidazole) revealed several suitable low-melting salts. Metathesis reactions of Na(azolate) were conducted by first using weakly acidic azoles including 4,5-dicyanoimidazole, 2-amino-4,5-dicyanoimidazole, 5-aminotetrazole, and 1,2,4-triazole to form sodium or lithium salts using group(I) hydroxides in methanolic solutions. The best results were obtained by reacting the basic and acidic azoles directly in 1:1 or 3:1 ratios at elevated temperatures. Twenty-two azole mixtures were identified which are either low melting solids or room temperature liquids. Each of the low melting solids was confirmed to be a new solid phase, each of which is being further characterized. The liquids and solids are anticipated to be ILs, eutectics, or partially ionized systems, all of which will be suitable for the dissolution of f-element salts. Five new synthetic methodologies were developed to finding suitable crystallization conditions for f-element complexation with the goal of finding simple one pot reaction syntheses and crystallization strategies that could be used under the demanding conditions of transuranic chemistries. These synthetic methods yield many new crystalline phases which were studied by single crystal X-ray diffraction.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Large-scale parametric modeling of spent nuclear fuel dynamics in the 30 cm package drop scenario

Packages used to transport spent nuclear fuel (SNF) are required by the U.S. Code of Federal Regulations 10 CFR 71.71 to demonstrate satisfactory performance during a drop scenario. While the CFR is meant to ensure safe package function, it does not evaluate survival of the SNF within. The U.S. Department of Energy Spent Fuel and Waste Science and Technology program is working on closing the knowledge gap related to the response of SNF to external mechanical loads, including the hypothetical 30 cm package drop scenario in the CFR. In support of this effort, LS-DYNA finite element simulations were developed by Pacific Northwest National Laboratory (PNNL) to model generic drop scenarios at both the package and fuel assembly level. The models were validated against one-third scale package and full scale fuel assembly drop test data and were exercised to predict fuel cladding strains in a narrow range of model configurations. This work describes a large-scale parametric study conducted by PNNL using the previously developed and validated PWR finite element model, with the addition of a new generic BWR assembly model. The motivation for the parametric study was to characterize the broad range of SNF responses in the 30 cm package drop scenario. This was accomplished by varying the drop orientation, fuel assembly type (17x17 PWR and 10x10 BWR), burnup, cladding temperature, spacer grid buckling load, package mass, impact limiter stiffness, and mechanical gap conditions within the basket. A MATLAB framework was developed to automate LS-DYNA model generation and execution on PNNL institutional computing resources. In total, over 2000 simulations were performed. For each simulation, the SNF response was quantified in terms of permanent grid deformation, fuel rod contact pressure, and strains within the fuel rods, guide tubes, and water rods. The results provide valuable insight into the range of responses that could be reasonably expected from SNF in the hypothetical drop scenario, as well as the sensitivity to each input parameter. The results of this parametric study are a key component of the testing and modeling strategy the Spent Fuel and Waste Science and Technology program is using to close the external loads knowledge gap.

Kadooka, Kevin↗

Predicting the evolution of biomass bulk density through feedstock preprocessing: Discrete element modeling, regression analysis, and pilot-scale validation

Bulk density is an important material property of biomass feedstocks, influencing handling, storage, transport costs, and conversion efficiency. In this study, predictive regression models for loose and tapped bulk densities of Alamo and Cave-in-Rock switchgrass are developed using a comprehensive dataset generated via calibrated bonded-sphere discrete element method (DEM) simulations. Here, a key contribution of this study is the use of a DEM-based approach, which correlates density with moisture content and particle size distribution parameters and enables analysis across a continuous particle size range, overcoming limitations of purely experimental data. For comparison, regression models are also developed using only experimental data from pilot-scale runs at the Biomass Feedstock National User Facility at Idaho National Laboratory. Validation against pilot-scale data showed reasonable prediction accuracy for both model types, particularly for smaller particle sizes (post-secondary grinding). While the experimental model showed slightly better performance matching the validation data in some cases, the DEM-based model benefits from a much larger dataset, reduced predictor multicollinearity, and continuous parameter coverage, highlighting the utility of validated simulation models for developing robust predictive tools for biomass preprocessing applications.

09 - BIOMASS FUELS↗

Thermal Model for Horizontal In-Water and In-Air Mini-Plate Experiments in the Advanced Test Reactor

Mini-Plate (MP) irradiation tests are a fueled experiment designed for irradiation in multiple test locations in the Advanced Test Reactor (ATR) at the Idaho National Laboratory (INL). MP experiments are considered a non-instrumented drop-in test where small aluminum-clad fuel plate samples are cooled directly by the ATR Primary Coolant System (PCS) water. These fueled experiments contain aluminum-clad fuel mini plates consisting of monolithic U-10Mo. Previous thermal analyses have been documented concerning oxide growth[1] and thermal safety margin[2] for the MP-2 experiment. A thermal analysis has been performed on the MP-2 experiment using a new technique for calculating the necessary time after the experiment has been irradiated and reactor shutdown necessary to not exceed the 95% lower bound blister threshold temperature for the fuel in a horizontal in air configuration. A horizontal in water model has also been analyzed and will be discussed. After the experiment is irradiated in the ATR, it is necessary to calculate the fuel meat temperatures at 12 hours after shutdown. This is done to assure that at the 12 hour mark, if this experiment is removed from the reactor and accidentally dropped onto the working tray and is in the most severe geometrical configuration of horizontal (minimum heat transfer) that fission products will not be released. The blister threshold temperature is considered this safety limit. The horizontal in air analysis is also performed considering that the experiment is horizontal in a drained canal and filled with air and natural convection heat transfer with air on the outside. The initial mini-plate testing campaign known as MP-1 [3] has been irradiated and its purpose was to form the base fuel specimens fabricated by various processes so that a single process can be down-selected for qualification. MP-2 will then serve to populate a comprehensive fuel performance data set for the down-selected fabrication process. The MP-2 experiment will be irradiated in several locations of the ATR as shown in Figure 1. The fuel plate thickness and geometry is identical to the MP-1 experiment [3], with the only variable being the fuel meat thickness. This model of the B-12 position is treated as encompassing the I-21, I-24, and B-9 positions since it has more fuel and higher total heat than the other positions. Since the blister threshold temperature is a function of fission density, a series of fission densities is calculated for the B-12 position to cover the other three positions. The B 12 position specifically uses fuel meat that is 0.016 in thick by 3.25 in long and 0.75 in tall. The total plate thickness is 0.050 in. Aluminum cladding forms the barrier between the PCS and the fuel meat. The purpose of this paper is to demonstrate this new technique developed to calculate the necessary decay time after shutdown for a series of fission densities using the ABAQUS [4] finite element and heat transfer code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Catalysts and methods for the conversion of carbonaceous materials to liquid fuels

The present disclosure relates to a composition that includes copper (Cu), aluminum (Al), oxygen, and an element (M) that includes at least one of magnesium, cerium, and/or a transition metal, where the copper and the element are present at a first molar ratio relative to the aluminum between about 0.1:1 and about 30:1 ((Cu+M):Al), and the copper and the element are present at a second molar ratio between about 0.1:4 and about 20:1 (Cu:M).

Vardon, Derek Richard↗

Reproducible benchmark for the SNAP 8 experimental reactor at dry conditions

Here, this work provides fully reproducible benchmark models of the Systems for Nuclear Auxiliary Power (SNAP) 8 Experimental Reactor. Validation benchmarks of the criticality configuration experiments under dry conditions with no coolant are presented. In addition, relevant reactivity worth experiments, and measurement of power distributions are validated and presented here. Discrepancies between modeled and experimental results are at most 110 pcm for fuel and poison worths, and critical configurations that do not manipulate control elements are within 50 pcm. Larger discrepancies found in reflector element worth experiments are due to experimental methodology which was noted as being simplified due to limited calculational capabilities at the time, while discrepancy in the control element manipulated critical configuration is due the absence of structural components for modeling simplification. Power distributions closely follow what was observed in experiment with expected peaking in reflecting elements. All models are well documented with cited references describing and justifying assumptions, simplifications, and adjustments to reproduce the models with any nuclear safety codes; model inputs and outputs are stored in the snapReactors GitHub repository.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Efficient computational search for lanthanide-binding additive dopants for advanced U-Zr based fuels

A rational design strategy combining theory and experiment is highly desirable for the development of U-Zr based fuels that are resistant to lanthanide-induced fuel-cladding chemical interaction (FCCI). Here, we performed rapid computational screening of FCCI-mitigating (via strong lanthanide-binding) elements across the periodic table by utilizing three simple criteria and density functional theory (DFT) calculated total energies of relevant elements and compounds retrieved from high-throughput DFT databases. In addition to successfully identifying previously demonstrated additives, our search leads to the discovery of Bi as a highly effective element among all 68 elements screened, and its effectiveness was experimentally confirmed in this study.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Biomanufacturing and bioprocessing of lunar regolith

Microbial biomanufacturing is important to accelerate lunar construction because it can leverage lunar material and waste streams as feedstocks to create a circular production system. In-space bio-mining and biomanufacturing using moon and asteroidal source material will enable the creation of infrastructure, produce industrial fuels and lubricants, and enable recovery of actinides and rare-earth elements (REEs) present in trace concentrations. Moreover, biomanufacturing in closed-loop systems (recycling and reuse of resources toward the establishment of a circular economy) will enable long-term lunar activities by recycling waste (CO 2 , gray water) and producing oxygen and biomaterials. Our response focuses on the use of lunar regolith and waste streams as feedstocks for protein and microbial-enabled biomining and bioprocessing to extract actinides and REEs, and to create biocomposites for lunar infrastructure. We envision an enclosed process that initiates with (1a) extracting actinides and REEs from lunar regolith using immobilized proteins, followed by (1b) creating biocomposites from the post-extracted lunar regolith for infrastructure, and (1c) cultivating diatoms and other microalgae on waste streams to harvest silica shells for incorporating into biocomposites and to generate O 2 for human respiration and/or producing refinable feedstocks. LLNL has significant expertise in all three processes and provides facilities, personnel, and expertise at the intersection of metal (lanthanide, actinide, transition) separations, purifications, biohydrometallurgy, radiobiochemistry, synthetic and systems biology, and materials science and engineering. Importantly, all three processes are relatively well-studied for Earth-based workflows and can be derisked for demonstration on the lunar surface by 2029.

59 BASIC BIOLOGICAL SCIENCES↗

LDRD poster - Identifying the speciation of salt-based actinides in the presence of contaminants

Molten salt reactors (MSRs) are among the most promising Generation IV nuclear reactors. The salts used in these reactors exhibit a high affinity for moisture and oxygen, and heating does not eliminate the oxides that are formed. Additionally, due to the high operating temperatures, reactor parts experience pronounced corrosion, and the resulting corrosion products integrate into the fuel salt system. It is hypothesized that the speciation of f-elements in molten salts would change in the presence of oxides and corrosion products. These changes in local structure can alter the rheological properties, which are critical to the operational and safety parameters of the reactors. Therefore, an in-depth understanding of how corrosion products influence lanthanide/actinide speciation is crucial. To test this hypothesis, lanthanides (as surrogates for actinides) and actinide chlorides were mixed with oxides (such as nickel oxide and lithium oxide) and/or corrosion products (e.g., nickel chloride) and reacted in molten salts. The resulting phases were analyzed using techniques such as powder X-ray diffraction, Raman spectroscopy, UV-Visible spectroscopy, and nuclear magnetic resonance spectroscopy. In the presence of oxides, due to the high oxophilicity of f-elements, these elements capture oxygen from the oxides, forming lanthanide/actinide oxide or oxychloride phases. Specifically, with cerium chloride, the introduction of lithium oxide to the formation of various final products depending on the oxide concentration. At lower oxide concentrations, cerium oxychloride (CeOCl) in the +3 oxidation state was observed, while an increase in the oxide concentration led to the formation of cerium dioxide (CeO2), where cerium is in the +4 oxidation state. In contrast, with nickel oxide, CeOCl formation was consistently observed, regardless of the oxide concentration. This project successfully determined that the speciation of f-elements change in the presence of oxides and different oxides in molten salts result in different final products. When the influence of nickel chloride was examined, changes in f-element speciation were not observed. However, it reacted with lithium chloride to generate a lithium-nickel-chloride (Li6NiCl8) phase, which was also observed when nickel oxide was reacted. Computational models were developed to simulate the phases anticipated under experimental conditions.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗

Development of Innovative Non-Concentrated Alloy (NCA) Fuel Cladding for Advanced Nuclear Reactors

As the nuclear energy sector advances toward next-generation reactors, the need for high-performance fuel cladding materials has become increasingly urgent. Traditional alloys like zirconium and stainless steel are reaching their performance limits under higher temperatures, more corrosive coolants, and extended irradiation. This report presents the development of a new class of fuel cladding materials based on Non-Concentrated Alloys (NCAs) ? multi-element systems designed to deliver enhanced mechanical strength, corrosion resistance, and radiation tolerance. Through a combination of computational modeling (Computer Coupling of Phase Diagrams and Thermochemistry: CALPHAD), simulation-guided alloy selection, and experimental fabrication via arc melting and spark plasma sintering (SPS), three strategic alloy design paths were explored: (1) FeCrAl-based NCAs, (2) refractory-lean neutron-efficient alloys, and (3) equimolar high-entropy compositions. Microstructural analysis confirmed the formation of stable body-centered cubic BCC_A2 phases, while mechanical testing demonstrated hardness values significantly exceeding those of conventional cladding materials. The results highlight the tunability of NCA systems and their potential for balancing strength and ductility ? a critical consideration for in-reactor performance. Looking forward, future work will focus on thermomechanical optimization, CALPHAD refinement, and benchmarking against industry standards to enable scalable deployment. This work not only advances the science of nuclear materials but also supports broader goals in nuclear safety, performance, and nuclear energy innovation.

36 - MATERIALS SCIENCE↗

GPHS Module Pre-Impact test Temperature Prediction with ANSYS [PowerPoint]

The purpose of the project was to: Develop a finite element thermal model of the General Purpose Heat Source (GPHS) Module loaded with four fueled clads; Perform steady state analysis to validate the thermal model via energy balance; Perform transient thermal analysis to predict clad temperatures at various times; and, Transient model results will be used to inform module soak time needed to reach a desired clad temperature prior to impact.

42 ENGINEERING↗

Comparing Sensor Fusion and Multimodal Chemometric Models for Monitoring U(VI) in Complex Environments Representative of Irradiated Nuclear Fuel

Optical sensors and chemometric models were leveraged for the quantification of uranium(VI) (0–100 μg mL –1 ), europium (0–150 μg mL –1 ), samarium (0–250 μg mL –1 ), praseodymium (0–350 μg mL –1 ), neodymium (0–1000 μg mL –1 ), and HNO 3 (2–4 M) with varying corrosion product (iron, nickel, and chromium) levels using laser fluorescence, Raman scattering, and ultraviolet–visible–near-infrared absorption spectra. In this paper, an efficient approach to developing and evaluating tens of thousands of partial least-squares regression (PLSR) models, built from fused optical spectra or multimodal acquisitions, is discussed. Each PLSR model was optimized with unique preprocessing combinations, and features were selected using genetic algorithm filters. The 7-factor D-optimal design training set contained just 55 samples to minimize the number of samples. The performance of PLSR models was evaluated by using an automated latent variable selection script. PLS1 regression models tailored to each species outperformed a global PLS2 model. PLS1 models built using fused spectra data and a multimodal (i.e., analyzed separately) approach yielded similar information, resulting in percent root-mean-square error of prediction values of 0.9–5.7% for the seven factors. Further, the optical techniques and data processing strategies established in this study allow for the direct analysis of numerous species without measuring luminescence lifetimes or relying on a standard addition approach, making it optimal for near-real-time, in situ measurements. Nuclear reactor modeling helped bound training set conditions and identified elemental ratios of lanthanide fission products to characterize the burnup of irradiated nuclear fuel. Leveraging fluorescence, spectrophotometry, experimental design, and chemometrics can enable the remote quantification and characterization of complex systems with numerous species, monitor system performance, help identify the source of materials, and enable rapid high-throughput experiments in a variety of industrial processes and fundamental studies.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Recovering high-purity uranyl nitrate from simulated used nuclear fuel dissolver solutions by crystallization: rejecting technetium

The separation of U from Tc and other problematic fission product elements like Mo and Ru, along with Sr, Zr, Cs, and Nd, has been achieved via the crystallization of uranyl nitrate hexahydrate (UNH). Rejection of technetium as pertechnetate anion ( 99 TcO 4 – ) is an especially important feature of this system, as it otherwise tends to follow U (VI) in extractive separations. It also raises the salient question regarding why this oxoanion cannot replace nitrate within the crystalline lattice of UNH. Results showed high-yield (>90 %), high-purity (>99 %) recovery of U as UNH from solutions containing 99 TcO 4 – by simple reduction of temperature from 60°C to 20°C. There was no observable interaction of 99 TcO 4 – with UO 2 2+ . The addition of other cations like, Sr 2+ , Zr 4+ , Cs + , and Nd 3+ , also did not form secondary, contaminant solid phases, leaving the > 99 % of the fission product elements in the mother liquor, while the U was recovered at > 90 %. Similarly, Mo and Ru, when added to the mixture, were shown to behave as the other fission-product elements, remaining in the mother liquor during crystallization. As a result, DFT calculations showed that, despite the higher binding strength of TcO 4 – , HMoO 4 – , and BiO 3 – with the UO 2 2+ cation compared to NO 3 – , the hydrogen-bonding network of the two coordinated ions and four waters of hydration in the UNH crystal structure is the driving force for the high specificity of this separation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Assessment Approach to Advanced Fuel Models

This document constitutes completion of the NEAMS milestone M3MS-21IN0201016, which is titled: Assessment of advanced fuel improvements. In this report we present: (1) simulations of separate effects creep tests for uranium alloy and HT-9; (2) examples of standardizing and updating metallic fuel assessments to include the latest models and key output figures of merit, which is coordinated with LANL to better evaluate advanced models; (3) examples of how the FIPD database at ANL is integrated into metallic fuel assessment cases; and finally (4) documentation of fuel-specific figures of merit and the BISON review process for assessments with best practices for preparing and reviewing finite element simulations in the appendix.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Marmot V2

MARMOT is a robust numerical tool for mesoscale modeling of fuel performance developed under the NEAMS Fuels technical area to predict the coevolution of microstructure and properties in fuel and cladding materials. MARMOT accomplishes this using the phase field method coupled with finite strain mechanics and heat conduction. MARMOT is based on the open source Multiphysics Object-Oriented Simulation Environment (MOOSE) and solves the coupled partial differential equations defining the physics using the finite element method. MARMOT is being developed in order to facilitate the development of improved materials models for fuel performance, but it is also being developed as a powerful tool in and of itself for the simulation of mesoscale fuel performance.

Aagesen, LarryK.↗

Field Test Report Neutron Scintillator Array Dry Storage Cask Scanner FY2024

During two weeks of Field Testing at the Idaho National Laboratory INTEC Cask Farm in July and August 2024, the LLNL Dry Storage Cask Scanner Array was lifted on top of an MC-10 dry storage fuel cask and operated to acquire neutron and gamma-ray data from the 24 fuel bundle positions. Neutron and gamma-ray data acquisition scans across the top of the cask of varying dwell times were performed July 15-18, 2024 and August 19-22, 2024 to evaluate the ability of the scanner data to reveal asymmetries in the fuel positions that reflect asymmetries in the MC-10 cask fuel bundle loading. The MC-10 cask 24 position fuel bundle loading at the INTEC Cask Farm is well documented, including the locations of six empty fuel bundle positions. This loading presents an opportunity to test the ability of the scanner system to detect diversion of spent fuel bundles as well as to validate the MC-10 cask MCNP modeling. The cask scanner array consists of six Stilbene crystal scintillator detectors and a linear actuator frame that moves the six detectors across the MC-10 dry storage cask to obtain data above each of the 24 fuel bundle positions. The detectors are connected to a pulse-shape discrimination data acquisition system capable of generating separate neutron and gamma-ray spectra for each detector and for each scan position. From the prior single detector Field Test in 2021 and iteration with MCNP modeling, the neutron and gamma-ray data were analyzed in multiple energy regions to identify an analysis method that would provide the strongest and most consistent signature of the asymmetric MC-10 cask fuel loading1 . From both the 2021 Field Test and the current Field Test results, the neutron capture gamma-ray count rate around 2.2 MeV provides the strongest signature of the asymmetric MC-10 cask fuel loading and has qualitative agreement with MCNP calculations. Counting all gamma-rays produces a similar signature. Neutrons emerging from the cask top are moderated and captured by the hydrogen in the polyethylene moderator and scintillator detector, producing a 2.2 MeV gamma ray which is seen in the scintillator gamma-ray spectrum. The count rate in the 2.2 MeV gamma-ray region is ~50 c/s, which is ~1000x higher than the ~0.05 n/s rate in the > 4MeV neutron region, and ~50x greater than the ~1 n/s rate in the neutrons > 500 keV region. Analysis of the 2.2 MeV neutron-capture Compton-scattered gamma-rays produces a statistically significant signature of the INTEC Cask Farm MC-10 asymmetric fuel loading. MCNP simulations indicate that the average neutron energy spectrum offers the potential to detect a large asymmetry from several missing bundles as well as individual missing fuel bundles. Testing this feature will require measurements on a cask with single missing elements.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Implement and Test 3D Mortar Contact in BISON

We leverage the extension of the generation of mortar segment meshes to three dimensions in MOOSE’s framework to extend thermomechanical modeling capabilities to problems with three dimensions. A modular approach to gap heat transfer physics using the mortar finite element method was created and documented, mechanical contact was extended to three dimensions—including frictional behavior, performance and ease of use were improved, and steps towards scalability of solid mechanics problems involving contact were taken. Many of these new developments are demonstrated in the simulation of 3D light-water reactor (LWR) problems, where the thermomechanical interface problem is solved using the mortar finite element method. Usage of the mortar framework has improved convergence in 2D problems and has enabled employing friction in 3D problems, of which we show results of a short, local stack of 3D pellets. Consequently, the benefits of mortar in terms of solution convergence and quality are extended to three dimensions. Section 2 discusses fundamental developments that enabled the simulation of practical mortar problems in three dimensions and other general improvements, including the reduction of the derivative container size, the modification of dual basis computations when edge dropping (lack of secondary element projection) takes place, the improvement of conditioning when employing the VCP in-edge dropping conditions, and code usability and quality improvements. These latter code enhancements include the migration of tests using “old” mortar contact constraints to using dual mortar with a semi-smooth Newton solution strategy and the reuse of lower dimensional domains for straightforwardly setting up a mortar thermomechanical LWR problem, i.e. the MOOSE action is employed for mechanical contact and the thermal LWR action is employed to capture the gas conductance, contact, and radiation components of gap heat transfer physics. Independently of the mortar LWR thermal action, we developed a modular approach to gap heat transfer that resides in MOOSE and can be leveraged, e.g., in metallic fuel problems. This approach, whose code design based on MOOSE’s user objects to model specific physics was proposed by the maintenance activity, is detailed in Section 3. Based on the dual mortar finite element method, the frictional contact constraints were extended to three dimensions. A block sheared in two directions in and out of contact with a rigid plane is employed in Section 4 to show the way the approach handles changes in frictional states (e.g. stick to slip) within a competitive number of Newton iterations. Equations and numerical results on the use of the VCP with Cartesian Lagrange multipliers, whose combination enables their direct condensation, are described in Section 5.3. Two-dimensional and three-dimensional BISON LWR simulations are discussed in Section 6. Particularly, a stack of five eccentric pellets with a surface defect is simulated and the effect of pellet-cladding friction is assessed. Finally, conclusions are outlined in Section 7.

42 ENGINEERING↗