SEARCH · Engineering Papers

Results for “Radiative capture”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS2011

Assessment of Radiative Heating Uncertainty for Hyperbolic Earth Entry

This paper investigates the shock-layer radiative heating uncertainty for hyperbolic Earth entry, with the main focus being a Mars return. In Part I of this work, a baseline simulation approach involving the LAURA Navier-Stokes code with coupled ablation and radiation is presented, with the HARA radiation code being used for the radiation predictions. Flight cases representative of peak-heating Mars or asteroid return are de ned and the strong influence of coupled ablation and radiation on their aerothermodynamic environments are shown. Structural uncertainties inherent in the baseline simulations are identified, with turbulence modeling, precursor absorption, grid convergence, and radiation transport uncertainties combining for a +34% and ..24% structural uncertainty on the radiative heating. A parametric uncertainty analysis, which assumes interval uncertainties, is presented. This analysis accounts for uncertainties in the radiation models as well as heat of formation uncertainties in the flow field model. Discussions and references are provided to support the uncertainty range chosen for each parameter. A parametric uncertainty of +47.3% and -28.3% is computed for the stagnation-point radiative heating for the 15 km/s Mars-return case. A breakdown of the largest individual uncertainty contributors is presented, which includes C3 Swings cross-section, photoionization edge shift, and Opacity Project atomic lines. Combining the structural and parametric uncertainty components results in a total uncertainty of +81.3% and ..52.3% for the Mars-return case. In Part II, the computational technique and uncertainty analysis presented in Part I are applied to 1960s era shock-tube and constricted-arc experimental cases. It is shown that experiments contain shock layer temperatures and radiative ux values relevant to the Mars-return cases of present interest. Comparisons between the predictions and measurements, accounting for the uncertainty in both, are made for a range of experiments. A measure of comparison quality is de ned, which consists of the percent overlap of the predicted uncertainty bar with the corresponding measurement uncertainty bar. For nearly all cases, this percent overlap is greater than zero, and for most of the higher temperature cases (T >13,000 K) it is greater than 50%. These favorable comparisons provide evidence that the baseline computational technique and uncertainty analysis presented in Part I are adequate for Mars-return simulations. In Part III, the computational technique and uncertainty analysis presented in Part I are applied to EAST shock-tube cases. These experimental cases contain wavelength dependent intensity measurements in a wavelength range that covers 60% of the radiative intensity for the 11 km/s, 5 m radius flight case studied in Part I. Comparisons between the predictions and EAST measurements are made for a range of experiments. The uncertainty analysis presented in Part I is applied to each prediction, and comparisons are made using the metrics defined in Part II. The agreement between predictions and measurements is excellent for velocities greater than 10.5 km/s. Both the wavelength dependent and wavelength integrated intensities agree within 30% for nearly all cases considered. This agreement provides confidence in the computational technique and uncertainty analysis presented in Part I, and provides further evidence that this approach is adequate for Mars-return simulations. Part IV of this paper reviews existing experimental data that include the influence of massive ablation on radiative heating. It is concluded that this existing data is not sufficient for the present uncertainty analysis. Experiments to capture the influence of massive ablation on radiation are suggested as future work, along with further studies of the radiative precursor and improvements in the radiation properties of ablation products.

Johnston, Christopher O.↗

DOE OSTI2021

A temperature profile diagnostic for radiation waves on OMEGA-60

Predicting and matching radiation wave propagation with computational models has proven difficult. Information provided by experiments studying radiation flow has been limited when only radiation breakout is measured. We have developed the COAX (co-axial) diagnostic platform to provide spatial temperature profiles of a radiation wave through low density foams as a more detailed constraint for simulations. COAX uses a standard, laser-driven OMEGA-60 halfraum to drive radiation down a titanium-laden silicon oxide foam. Point-projection X-ray absorption spectroscopy perpendicular to the radiation flow measures the spatial profile of titanium ionization. The spectroscopic measurement utilizes a broadband capsule backlighter. Imaging and streak spectroscopy are used to characterize the size and spectrum of this source. Radiography provides an additional constraint by capturing the developing shock as the radiation flow becomes subsonic. The DANTE diagnostic is used to measure the halfraum temperature. Here, we provide a spectroscopic analysis of COAX data to determine temperature, and we describe experimental sources of uncertainty. The temperature is obtained by comparison to multi-temperature synthetic spectra post-processed from radiation-hydrodynamics simulations. Quantitative comparison between data and synthetic spectra generated from temperature profiles at relevant simulation times enable determination of a peak temperature of 114 ± 8 eV at 265 ± 22.4 μm from the halfraum. This represents an improvement over the temperature uncertainties of previous radiation flow experiments. Further refinements to the spectroscopic analysis could achieve ±4 eV. The combination between space-resolved spectroscopy and radiography enables us to determine the distance from the halfraum of both the radiation front and the shock front at the time of measurement. For the example shown in this paper the radiation front position is 600–630 μm at 3.43 ± 0.16 ns and the shock front position is 633 μm at 3.3 ± 0.24 ns.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Theoretical estimates of flammability bounds for thin condensed fuel diffusion flames in microgravity using detailed models of chemistry and radiation

Recently, U-shape flammability maps have been constructed showing minimal oxygen vs. flame strain for opposed flame spread in micro-gravity by Olson and Ferkul. The U-shape defines the limiting flammability bounds from radiative extinction and flame blow-off. Here, the minimum of the U corresponds to the minimum possible oxidizer concentration where burning can occur, and is an important quantity of interest for fire safety. While high strain extinction bounds have been well analyzed, low strain radiative extinction has not. To estimate low strain extinction, in this study an analytical theory is developed based on thin flame theory coupled with a heat and mass transfer model for solid fuels. A reaction progress variable based on the Damköhler number is adapted in the theory to account for incomplete combustion at high strain rates and enable the capturing of the full flammability map. The analytical model is compared to a one dimensional numerical model w/ detailed chemical kinetics and coupled radiation heat transfer in planar and spherical geometries. The flammability maps are then qualitatively compared to experimental extinguishment data compiled by Olson and Ferkul for cylindrical rods of PMMA showing similar trends. The results show the newly developed analytics capture the radiative extinction bound compared to the numerical model and qualitatively agrees with microgravity data.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

NASA NTRS2017

Device and Method of Scintillating Quantum Dots for Radiation Imaging

A radiation imaging device includes a radiation source and a micro structured detector comprising a material defining a surface that faces the radiation source. The material includes a plurality of discreet cavities having openings in the surface. The detector also includes a plurality of quantum dots disclosed in the cavities. The quantum dots are configured to interact with radiation from the radiation source, and to emit visible photons that indicate the presence of radiation. A digital camera and optics may be used to capture images formed by the detector in response to exposure to radiation.

Burke, Eric R.↗

NASA NTRS2011

The Properties of Large Amplitude Whistler Mode Waves in the Magnetosphere: Propagation and Relationship with Geomagnetic Activity

Wepresent resultsof a studyof the characteristicsof very large amplitude whistler mode waves inside the terrestrial magnetosphere at radial distances of less than 15 RE using waveform capture data from the Wind spacecraft. We observed 247 whistler mode waves with at least one electric field component (105/247 had !80 mV/m peak!to!peak amplitudes) and 66 whistler mode waves with at least one search coil magnetic field component (38/66 had !0.8 nT peak!to!peak amplitudes). Wave vectors determined from events with three magnetic field components indicate that 30/46 propagate within 20 of the ambient magnetic field, though some are more oblique (up to "50 ). No relationship was observed between wave normal angle and GSM latitude. 162/247 of the large amplitude whistler mode waves were observed during magnetically active periods (AE > 200 nT). 217 out of 247 total whistler mode waves examined were observed inside the radiation belts. We present a waveform capture with the largest whistler wave magnetic field amplitude (^8 nT peak!to!peak) ever reported in the radiation belts. The estimated Poynting flux magnitude associated with this wave is ^300 mW/m2, roughly four orders of magnitude above estimates from previous satellite measurements. Such large Poynting flux values are consistent with rapid energization of electrons.

Wilson, L. B., III↗

NASA NTRS2005

A New Perspective on Trapped Radiation Belts in Planetary Atmospheres

The charged particle fluxes trapped in the magnetic dipole fields of certain planets in our Solar System are interesting signatures of planetary properties in space physics. They also represent a source of potentially hazardous radiation to spacecraft during planetary and interplanetary exploration. The Earth s trapped radiation belts have been studied for years and the physical mechanisms by which primary radiation from the Sun and Galaxy is captured is well understood. The higher-energy particles collide with molecules in the planetary atmosphere and initiate large cascades of secondary radiation which itself becomes trapped by the magnetic dipole field of the planet. Some of it is even backscattered as albedo neutrons.

Diaz, A.↗

DOE OSTI2021

Reduction of the neutron imaginary potential off the stability line and its possible impact on neutron capture rates

The effect of the isovector imaginary optical potential has been studied experimentally by using the particle-evaporation technique for the 11 B + 48 Ca reaction with a 21.8 MeV 11 B beam. Spectra of neutron, proton, and α particles emitted from the neutron-rich compound nucleus 59 Mn have been measured and analyzed with traditional optical model potentials with their original parametrizations as well as with adjusted isovector imaginary components. It is shown that the isovector component of the imaginary potential is indispensable in the reproduction of proton and α-particle yields from this reaction and even needs to be enhanced compared with the suggestions of the original model parametrizations. Furthermore, this can lead to important consequences for astrophysical reaction-rate calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Nuclear Reaction Data Evaluations for Nuclear Applications Physics and modeling [Slides]

This presentation provides an evaluation workflow scheme for the project. It also discusses nuclear data libraries, including an overview of worldwide nuclear data libraries and current funding agencies. The presentation also discusses the theoretical framework commonly used in the evaluated work, including the R-matrix theory (≈ up to a few keVs for heavy nuclei) as well as pre-equilibrium (PE) and Hauser-Feshbach (HF) theory (⪆1 MeV for medium/heavy nuclei). Inclusion of experimental effects are also discussed. Examples of evaluated work, including theory, experiments, and applications are discussed, specifically n+ 28,29,30 Si evaluation (direct capture component), n+ 53 Cr evaluation (strong multiple scattering effects), n+ 239 Pu(n,γf) reaction (fluctuations in neutron multiplicities), and α+ 17,18 O evaluations (emerging neutron spectra).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRSDate not supplied

Comparison between Hayabusa 2 Spectral Measurements and Simulations

Recent Hayabusa 2 spectral measurements are compared with state-of-the-art shock-layer radiation simulations resulting from the LAURA/HARA code suite. These simulations include coupled ablation, which accounts for the injection of ablation products into the flowfield, and coupled radiation, which accounts for radiative energy loss in the flowfield. To enable the cou- pled ablation simulations, a best-estimate model is developed for Hayabusa’s carbon-phenolic ablator based on available published information. The comparison between the simulations and measurements focuses on three atomic nitrogen lines, two atomic oxygen lines, one atomic hydrogen line, and the CN Violet band system. Parametric uncertainties are evaluated for these features using a recently developed capability in LAURA/HARA. These uncertainties capture the impact of flowfield and radiation modeling uncertainties on the observed radiation simulations, and include the impact of coupled radiation and ablation product emission. The resulting parametric uncertainties range from 50 to over 100% for the atomic lines, and 80% for CN Violet. Uncertainties in flowfield kinetics and non-Boltzmann modeling provide the dominant contribution to these uncertainties. For the atomic nitrogen and oxygen lines, com- parisons between the nominal simulations and measurements show agreement well within these simulation uncertainty bounds, with the measurements within 20% of the nominal simulation over most of the trajectory, and within 10% at peak emission. This excellent agreement is un- precedented for observed radiation measurements of atomic lines, where previous Stardust and Hayabusa 1 comparisons were significantly worse. The improved agreement for these Hayabusa 2 comparisons is both the result of improved measurement quality and enhanced flowfield and radiation modeling. For the ablation product emission from the 656 nm atomic hydrogen line, the agreement between measurements and simulations is within ±20% for all trajectory points except one outlier. This good comparison is surprising considering this species depends on the ablation rate and is shown to have strong wake contribution, which both contribute to the large ≈±70% parametric uncertainty. For the ablation product emission from CN Violet, the simulations over-predict the measurements by over a factor of two early in the trajectory, with the comparison improving later in the trajectory. This disagreement for CN Violet is consistent with the ≈±80% parametric uncertainty bounds evaluated for the simulations. To utilize the excellent agreement between the atomic line simulations and measurements in the assessment of the radiative heating margin for a future flight vehicle, the flowfield property binning (FPB) approach is applied. Using the Mars Sample Return (MSR) Earth Entry System (EES) as an example, the combination of the excellent Hayabusa comparisons and the FPB analysis justify nearly a 10% reduction in the radiative heating margin from an assumed baseline margin of 30%. This direct quantitative link between observed radiation measurements and the radia- tive heating margin of future flight vehicles demonstrates the value of these observed radiation measurements.

Shock-Layer Radiative Heating↗

DOE OSTI2020

Discriminating Uranium Isotopes Based on Fission Signatures Induced by Delayed Neutrons

The use of active interrogation to induce delayed neutron emission is a well-established technique for the characterization of special nuclear materials. Delayed neutrons have isotope-characteristic spectral and temporal signatures, which provide the basis for isotope identification. However, in bulk materials that contain an appreciable fissile (e.g., 235 U or 233 U) fraction, such as highly enriched uranium (HEU), delayed neutrons have a high probability of inducing additional fissions. As a result, the overall delayed neutron signature consists of two distinct components: the “primary” delayed neutrons (emitted directly by fission fragments) and the “secondary prompt” fission neutrons produced in fission induced by primary delayed neutrons. These prompt products differ from “primary” delayed neutrons both in their energy spectra and in the presence of coincident radiation released by the parent fission event. The presence and relative quantity of prompt products from delayed fission depend on the cross section of the material in the energy range of delayed neutrons, which may differ significantly between isotopes, thus providing an exploitable means for isotope differentiation. Here, we demonstrate two experimental approaches for discriminating between 235 U and 238 U isotopes based on the measurement of delayed neutron-induced fission products. First, HEU and depleted uranium objects are differentiated through the detection of high-energy prompt neutrons from delayed fission using both recoil-based organic liquid scintillators and thermalization spectra from a custom-built capture-gated composite detector. Secondly, coincident radiation measurements are used as the basis for discrimination by comparing the overall rates and time evolution of fission events when delayed neutrons are present.

07 ISOTOPE AND RADIATION SOURCES↗

DOE OSTI2024

Electron tail suppression and effective collisionality due to synchrotron emission and absorption in mildly relativistic plasmas

Synchrotron radiation losses are a significant cause of concern for high-temperature aneutronic fusion reactions such as proton–Boron 11. The fact that radiation losses occur primarily in the high-energy tail, where the radiation itself has a substantial impact on the electron distribution, necessitates a self-consistent approach to modeling the diffusion and drag induced by synchrotron absorption and emission. Furthermore, an accurate model must account for the fact that the radiation emission spectrum is momentum-dependent, and the plasma opacity is frequency-dependent. Here, we present a simple Fokker–Planck operator, built on a newly solved-for blackbody synchrotron diffusion operator, which captures all relevant features of the synchrotron radiation. Focusing on magnetic mirror fusion plasmas, we show that significant suppression of the electron distribution occurs for relativistic values of the perpendicular electron momentum, which therefore emit much less radiation than predicted under the assumption of a Maxwell–Jüttner distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Fusion burn-propagation simulations using the collisional and radiative particle-in-cell code TRIFORCE

The ability to accurately model burn propagation in inertial confinement fusion plasmas is crucial for advancing fusion energy research. This work presents enhancements to the triforce hybrid fluid-kinetic multiphysics code, focusing on its kinetic half, which employs the particle-in-cell (PIC) method with Monte Carlo collisions (MCC). We use a moment-preserving collision model that mitigates numerical noise, particularly in spherical geometries where particle weights vary significantly. Additionally, we refine the treatment of inverse bremsstrahlung to account for electron–ion collision frequency reductions in degenerate plasmas and incorporate a blackbody radiation source to enable realistic photon injection. These improvements enable the simulation of 1-dimensional (1D) spherical fusion burn propagation in deuterium–tritium plasmas. Benchmark comparisons with the hydra radiation-hydrodynamics code confirm that triforce accurately captures the dynamics of hot-spot expansion and burn propagation, demonstrating sensitivity to ignition thresholds consistent with theoretical models. Findings show the ignition cliff to be less steep in our work compared to radiation-hydrodynamic modeling. These results highlight the role of kinetic effects in fusion ignition physics and underscore the necessity of hybrid fluid-kinetic models for advancing predictive capabilities in high-energy-density plasma systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

The cost savings potential of controlling solar thermal collectors with storage for time-of-use electricity rates

Homes and commercial buildings account for roughly 40% of U.S. energy consumption. A significant portion of this energy consumption, and therefore utility cost, comes from domestic hot-water use. Solar collectors such as flat plate collectors, evacuated tubes, and combined photovoltaic/thermal panels, are renewable options that collect thermal energy and could help reduce utility costs. Hot water is generated by pumping a working fluid through the solar collector, where the fluid captures thermal energy via absorbed solar radiation incident on the collector surface. The collection of heat coincides with the availability of solar radiation, but the best time to use this thermal energy depends on the utility rate at a given time, which often changes throughout the day. In this work, we implement and simulate a control strategy using Transient System Simulation Tool (TRNSYS) that reduces utility costs by controlling when the system stores and uses thermal energy based on solar availability and time-of-use electricity rates. Our modeled results show that the control strategy can reduce utility costs associated with using hot water between 3% and 30%. Here, we also show how the on-peak utility rates, the size and type of the solar collectors, and the location's climate impact the results.

14 SOLAR ENERGY↗

DOE OSTI2022

Natural and anthropogenic contributions to the hurricane drought of the 1970s–1980s

Atlantic hurricane activity experienced a pronounced lull during the 1970s and 1980s. The current explanation that anthropogenic aerosol radiative forcing cooled the sea surface locally fails to capture the magnitude of this large decrease in activity. To explain this hurricane drought, we propose that the radiative effects of sulfate aerosols from Europe and North-America decreased precipitation in the Sahara-Sahel region, leading to an enhancement of dust regional emissions and transport over the Atlantic. This dust in turn enhanced the local decrease of sea-surface temperature and of hurricane activity. Here, we show that dust emissions from the Sahara peaked in phase with regional sulfate aerosol optical thickness and Sahel drought conditions, and that dust optical depth variations alone can explain nearly half of the sea-surface temperature depression in the 1970s and 1980s.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Physics insights from a large-scale 2D UEDGE simulation database for detachment control in KSTAR

A large-scale database of two-dimensional UEDGE simulations has been developed to study detachment physics in KSTAR and to support surrogate models for control applications. Nearly 70 000 steady-state solutions were generated, systematically scanning upstream density, input power, plasma current, impurity fraction, and anomalous transport coefficients, with magnetic and electric drifts across the magnetic field included. The database identifies robust detachment indicators, with strike-point electron temperature at detachment onset consistently T e,target ~ 3-4 eV, largely insensitive to upstream conditions. Scaling relations reveal weaker impurity sensitivity than one-dimensional models and show that heat flux widths follow Eich’s scaling only for uniform, low D and χ. Distinctive in–out divertor asymmetries are observed in KSTAR, differing qualitatively from DIII-D. Complementary time-dependent simulations quantify plasma response to gas puffing, with delays of 5 - 15 ms at the outer strike point and ∼ 40 ms for the low-magnetic-field-side radiation front. These dynamics are well captured by first-order-plus-dead-time models and are consistent with experimentally observed detachment-control behavior in KSTAR (Gupta et al 2025 Plasma Phys. Control. Fusion (submitted)).

Physics - Plasma physics↗

DOE OSTI2022

An iron evaluation story: from TALYS model parameters to validation on the ASPIS benchmark with the Monte-Carlo code TRIPOLI-4

Iron is a widespread structural material in nuclear applications. However some benchmark studies show that for many libraries, including TENDL-2019, producing good iron evaluations is still an open challenge. This work aims at improving TENDL-2019 {sup 56}Fe and {sup 54}Fe evaluations through feedbacks from experimental data and Monte Carlo results of the iron benchmark ASPIS. Part of the T6 code infrastructure was used to produce the evaluations. In the new evaluations, particular attention to resonance parameters and formalism was paid to address the importance of elastic scattering. The total inelastic scattering and inelastic continuum cross sections were shaped by modifying the input parameters of TALYS, one of the T6 codes. Additionally, it was necessary to normalise some reaction channels to those of JENDL-4.0. With the new files ASPIS results, obtained with the Monte Carlo code TRIPOLI-4, improved noticeably compared to TENDL 2019. The evaluations and parameter setting modifications are to be included into TENDL-2021. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

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↗