Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “boundary element method”

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 199 records · Page 11

An Efficient B-Spline Lagrangian/Eulerian Method for Compressible Flow, Shock Waves, and Fracturing Solids

This study presents a new method for modeling the interaction between compressible flow, shock waves, and deformable structures, emphasizing destructive dynamics. Extending advances in time-splitting compressible flow and the Material Point Methods (MPM), we develop a hybrid Eulerian and Lagrangian/Eulerian scheme for monolithic flow-structure interactions. We adopt the second-order WENO scheme to advance the continuity equation. To stably resolve deforming boundaries with sub-cell particles, we propose a blending treatment of reflective and passable boundary conditions inspired by the theory of porous media. The strongly coupled velocity-pressure system is discretized with a new mixed-order finite element formulation employing B-spline shape functions. Shock wave propagation, temperature/density-induced buoyancy effects, and topology changes in solids are unitedly captured.

97 MATHEMATICS AND COMPUTING↗

Correlative STEM-APT characterization of radiation-induced segregation and precipitation of in-service BWR 304 stainless steel

In this study, radiation induced segregation and precipitation phenomena in an in-service boiling water reactor 304 stainless steel component were investigated using directly correlated 3D-atom probe tomography and scanning transmission electron microscopy. Significant quantitative differences in measured segregation at grain boundaries were found between the atom probe and energy dispersive spectroscopy measurements of the exact same locations. In particular, a much stronger Si segregation (~10 atomic% via atom probe versus ~4 atomic% via electron microscopy) and different Cr profile shapes were detected that are critical to models of radiation induced segregation and stress corrosion cracking behavior. These quantitative differences highlight the need for comparative microscopy and critical evaluation of limitations in each analytical method. Elemental segregation to dislocations and conjoined-clusters were also highlighted by atom probe; confirming and expanding upon what has been observed in test reactor neutron and accelerator-based ion irradiations.

36 MATERIALS SCIENCE↗

Seismic response of vertical dry storage casks under three-dimensional earthquake motions

Ensuring the long-term seismic safety of dry storage casks (DSCs) is becoming increasingly critical as these systems evolve from temporary to de facto permanent repositories for spent nuclear fuels. Traditional seismic soil–structure interaction (SSI) assessment methods use one-dimensional deconvolution or simplified boundary conditions to model incident waves. Although computationally appealing, simplifying assumptions may alter the seismic risk by neglecting the full complexity of three-dimensional (3D) wave propagation effects. To address this challenge, this paper introduces a novel high-fidelity computational framework that leverages the Domain Reduction Method (DRM) with perfectly matched layers (PML) to accurately transfer complex, 3D seismic wavefields from regional-scale fault-rupture simulations into local-scale finite element models of DSCs. Using broadband, physics-based ground motions from a generic M w 7.0 strike-slip event, both single-cask and multi-cask configurations were investigated under near- and far-field conditions. Emphasis is placed on capturing complex SSI, spatial variability in the ground motion, and nonlinear phenomena such as cask rocking and sliding. Numerical results demonstrate that near-field conditions, where forward directivity and fling-step effects dominate, lead to significantly higher DSC rocking and sliding. Far-field cases, by contrast, generally exhibit modest responses. Incorporating SSI tends to amplify or alter DSC response spectra and introduce response variability, which underscores the need for site-specific evaluations and robust modeling approaches to ensure the seismic integrity of DSCs in interim spent fuel storage installations.

Das, Tonmoy↗

Sample shape and boundary dependence of measured transverse thermal properties

In spite of increased interest in thermal Hall measurements for the analysis of insulating quantum materials, there remains large uncertainty in such measurements due to contact misalignment. In this paper, we propose that sample geometry and uncertain boundary conditions provide an additional source of uncertainty in the measurement of $D_{xy}$ or $κ_{xy}$. By running simple simulations in an open source finite-element solver, we demonstrate that measured $D_{xy}$ can be changed by a non-negligible fraction in samples with similar width and length. This geometric corrective factor depends on the distinction between a uniform heat flow and constant temperature boundary couplings to a bath. Sample geometry and boundary conditions can be accounted for through simulation or by using rectangular samples to make thermal Hall measurements more reliable and reproducible. Finally, we detail a contactless optical method for measuring $D_{xy}$ based on an existing photothermal microscope technique. This method is insensitive to the longitudinal diffusivity pollution caused by contact misalignment.

36 MATERIALS SCIENCE↗

Interpreting strain tensor data to characterize and monitor reservoirs for CO2 storage and other applications

Recent advances in instrumentation have made it feasible to measure the transient strain tensor caused by small changes in fluid volume or pressure in the subsurface and this has opened the door to new opportunities for characterization and monitoring. We have deployed strainmeters and then conducted injection well tests in an underlying reservoir at 530m depth. The resulting data indicated that the horizontal strain at shallow strainmeters (30 to 40m depth) was tensile and the vertical strain was compressive. The radial strain was less than the horizontal strain, and the strain rates decreased from 100 nanostrain/day to roughly 10 n/d over a few days. The signal at two strainmeters at shallow depth were consistent, although the magnitude of the horizontal strains were different reflecting the different radial directions from the well. The signal at a deep strainmeter deployed at reservoir depth was much different, with tensile vertical strains and compressive horizontal strains. These data can be interpreted by inverting poroelastic forward models developed using numerical and analytical methods. The average horizontal strain in the caprock resembles the transient pressure in the underlying reservoir and classic type-curve methods from transient well testing can be used for preliminary interpretations of strain data. We have developed closed-form analytical solutions to a pressurized poroelastic inclusion and inhomogeneity in a half-space. This model is fast and can be inverted to estimate reservoir stiffness and geometry. Numerical models developed using finite element methods allow more details of the subsurface to be included in the inversion, but they require much longer run times and this makes inversion cumbersome using standard methods. We have developed an inversion approach that uses a proxy model created using machine learning to do most of the forward calculations. The proxy model is periodically updated and refined using the finite element model to ensure accuracy. This approach significantly reduces the computational requirements and makes it feasible to use Bayesian inversion with large numerical models. We have shown that the strain tensor in the caprock is sensitive to pressure in the reservoir, boundaries in the reservoir, and pressure in the caprock caused by leaks. These results indicate that coupling strain tensor data with inversion has the potential to help evaluate reservoirs during initial characterization, and to monitor them during the CO2 injection and storage process.

Murdoch, Larry↗

Local and long-range atomic/magnetic structure of non-stoichiometric spinel iron oxide nanocrystallites

Spinel iron oxide nanoparticles of different mean sizes in the range 10–25 nm have been prepared by surfactant-free up-scalable near- and super-critical hydrothermal synthesis pathways and characterized using a wide range of advanced structural characterization methods to provide a highly detailed structural description. The atomic structure is examined by combined Rietveld analysis of synchrotron powder X-ray diffraction (PXRD) data and time-of-flight neutron powder-diffraction (NPD) data. The local atomic ordering is further analysed by pair distribution function (PDF) analysis of both X-ray and neutron total-scattering data. It is observed that a non-stoichiometric structural model based on a tetragonal γ-Fe 2 O 3 phase with vacancy ordering in the structure (space group P 4 3 2 1 2) yields the best fit to the PXRD and total-scattering data. Detailed peak-profile analysis reveals a shorter coherence length for the superstructure, which may be attributed to the vacancy-ordered domains being smaller than the size of the crystallites and/or the presence of anti-phase boundaries, faulting or other disorder effects. The intermediate stoichiometry between that of γ-Fe 2 O 3 and Fe 3 O 4 is confirmed by refinement of the Fe/O stoichiometry in the scattering data and quantitative analysis of Mössbauer spectra. The structural characterization is complemented by nano/micro-structural analysis using transmission electron microscopy (TEM), elemental mapping using scanning TEM, energy-dispersive X-ray spectroscopy and the measurement of macroscopic magnetic properties using vibrating sample magnetometry. Notably, no evidence is found of a Fe 3 O 4 /γ-Fe 2 O 3 core-shell nanostructure being present, which had previously been suggested for non-stoichiometric spinel iron oxide nanoparticles. Finally, the study is concluded using the magnetic PDF (mPDF) method to model the neutron total-scattering data and determine the local magnetic ordering and magnetic domain sizes in the iron oxide nanoparticles. The mPDF data analysis reveals ferrimagnetic collinear ordering of the spins in the structure and the magnetic domain sizes to be ∼60–70% of the total nanoparticle sizes. The present study is the first in which mPDF analysis has been applied to magnetic nanoparticles, establishing a successful precedent for future studies of magnetic nanoparticles using this technique.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Measurement of grain boundary strength of Inconel X-750 superalloy using in-situ micro-tensile testing techniques in FIB/SEM system

Grain boundaries (GBs), known as two-dimensional defects, are omnipresent in polycrystalline metallic alloys and thus influence a wide range of mechanical properties under different environmental conditions like irradiation and corrosion. Therefore, quantifying the strength of individual GBs is critical for understanding the degradation of mechanical properties of materials under different conditions. In this study we developed an efficient approach for the fabrication of micro-tensile specimens with a GB almost perpendicular to the tensile direction, which is expected to advance the development of individual GB tensile testing at micro or nanoscale in a wide scope of materials. Further, an in-situ cantilever micro-tensile testing method was developed and used to quantify the strength of a Σ3 GB in Inconel X-750 with the combination of finite element modeling. The average ultimate tensile strength (UTS) of a non-irradiated Σ3 GB is estimated at around 1.4 GPa, comparable to that of a neutron-irradiated Σ3 GB with a dose of ~1.5 dpa (1.3 GPa). Moreover, the in-situ push-to-pull micro-tensile testing technique developed in this work provides valuable insights into the high-angle GB deformation and fracture behavior. This method generates qualitatively similar ductility behavior before and after neutron irradiation as the bulk material testing. However, the ductility and UTS values obtained from this method are different from bulk measurements due to vastly different specimen dimensions.

36 MATERIALS SCIENCE↗

Stable, Ductile and Strong Ultrafine HT-9 Steels via Large Strain Machining

Beyond the current commercial materials, refining the grain size is among the proposed strategies to manufacture resilient materials for industrial applications demanding high resistance to severe environments. Here, large strain machining (LSM) was used to manufacture nanostructured HT-9 steel with enhanced thermal stability, mechanical properties, and ductility. Nanocrystalline HT-9 steels with different aspect rations are achieved. In-situ transmission electron microscopy annealing experiments demonstrated that the nanocrystalline grains have excellent thermal stability up to 700 °C with no additional elemental segregation on the grain boundaries other than the initial carbides, attributing the thermal stability of the LSM materials to the low dislocation densities and strains in the final microstructure. Nano-indentation and micro-tensile testing performed on the LSM material pre- and post-annealing demonstrated the possibility of tuning the material’s strength and ductility. The results expound on the possibility of manufacturing controlled nanocrystalline materials via a scalable and cost-effective method, albeit with additional fundamental understanding of the resultant morphology dependence on the LSM conditions.

36 MATERIALS SCIENCE↗

Multilevel Monte Carlo methods for the Grad-Shafranov free boundary problem

The equilibrium configuration of a plasma in an axially symmetric reactor is described mathematically by a free boundary problem associated with the celebrated Grad-Shafranov equation. The presence of uncertainty in the model parameters introduces the need to quantify the variability in the predictions. This is often done by computing a large number of model solutions on a computational grid for an ensemble of parameter values and then obtaining estimates for the statistical properties of solutions. In this study, we explore the savings that can be obtained using multilevel Monte Carlo methods, which reduce costs by performing the bulk of the computations on a sequence of spatial grids that are coarser than the one that would typically be used for a simple Monte Carlo simulation. We examine this approach using both a set of uniformly refined grids and a set of adaptively refined grids guided by a discrete error estimator. Numerical experiments show that multilevel methods dramatically reduce the cost of simulation, with cost reductions typically on the order of 60 or more and possibly as large as 200. Furthermore, adaptive griding results in more accurate computation of geometric quantities such as x-points associated with the model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mechanical Characterization of the NIF Ignition Target TMPSA Bonding Flexure

The Thermo-Mechanical Package Sub-Assembly (TMPSA) provides a critical mechanical, thermal, and electrical interface between a silicon arm and a TMP aluminum can in a NIF ignition target. During assembly, sixteen silicon pads are bonded to the aluminum using a fixture that positions the components and applies a repeatable prescribed displacement through a compliant flexure. The flexure converts fixture interference into displacement and reaction force. Because bondline thickness variation must be maintained within +/-4 µm, consistent flexure behavior is important to the assembly process. With a recent string of TMPSAs exhibiting low bond strength, the flexures were inspected. Despite being manufactured to the same specifications, flexures were found to exhibit variation in measured stiffness. Additionally, the measured stiffness values did not always follow the presumed beam mechanics model. This work addresses two questions: whether the bonding fixture is working as intended, and whether the previously made stiffness measurements are accurate representations of the use case. The flexures are analyzed using Euler-Bernoulli beam theory, measured beam dimensions, finite element analysis, and tolerance stack-up calculations. The analysis shows that the fixed-guided beam mechanics model appropriately represents the flexure during TMPSA bonding, but the chisel-tip stiffness measurement method introduces a deformation to the inner ring of the flexure that is not represented during use. The measured stiffness values should therefore be interpreted as test-condition stiffness values rather than direct measurements of operational flexure stiffness. The discrepancy is therefore attributed primarily to the measurement boundary condition rather than to failure of the fixed-guided beam model. Recommendations are provided for GD&T, dimensional inspection, and a representative stiffness testing method to better control bondline variation.

42 ENGINEERING↗

A numerical-homogenization based phase-field fracture modeling of linear elastic heterogeneous porous media

Most porous media, such as geomaterials and biomaterials are highly heterogeneous in nature, and they contain large variations of microscopic pore structures, such as pore sizes, pore distribution, and pore shapes. The oscillation of microscopic structures is a substantial challenge in theoretical characterization and is usually ignored in continuous modeling. However, mechanical behavior of porous media such as deformation and failure, are essentially impacted by the microscopic heterogeneity which needs to be considered in modeling a porous media. Here, this research proposes a numerical modeling framework with a capability to investigate the effect of microscopic heterogeneity on the macroscopic fracture behavior in porous media by using a numerical homogenization technique, combined with the phase-field fracture modeling method. This numerical modeling strategy computes a homogenized elasticity tensor based on microscopic heterogeneous pore structures heterogenous porous domain by solving boundary value problems at microscopic domain. The strain energy and subsequent propagation of macroscopic fractures will be updated using homogenized stiffness information. Using this numerical scheme, the microscopic pore structure’s impact on the fracture behavior through the homogenized elastic tensor will be taken into account. This multiscale technique is benchmarked against classical problems. Finally, the results highlight the importance of the underlying pore structure and reveal that both fracture strength and propagation path can be influenced by the microscopic heterogeneity.

36 MATERIALS SCIENCE↗

Modeling low cycle fatigue (LCF) of additively manufactured Hastelloy X using An accelerated crystal plasticity fatigue damage model

This paper presents a microstructure-based model for low cycle fatigue (LCF) behavior and life of Nickel-based alloy Hastelloy X manufactured using laser-powder bed fusion (L-PBF) additive manufacturing (AM). AM Hastelloy X, a solution-strengthened alloy, is tested at elevated temperature under fully reversed LCF conditions at different strain levels. A generalized plane strain finite element model is generated from electron backscatter diffraction (EBSD) characterization. The constitutive behavior of the material under fatigue is modeled using crystal plasticity and calibrated with both monotonic tensile and cyclic stress–strain data. The fatigue micro-crack initiation and propagation in the microstructure is modeled using a modified Chaboche fatigue damage model. An embedded boundary condition with a homogenous medium is used to apply the cyclic deformation and prevent numerically introduced over-constraints during fatigue simulation. A ‘cycle-jump’ method is used to accelerate the fatigue simulation and reduce the computational cost. The simulation results are compared to LCF experiments, showing satisfactory matches in cyclic stress behavior and number of cycles to macro-crack initiation for all applied strain ranges. In addition, the model illustrates the potential for quantifying microscale fatigue life impacting factors such as microstructure and surface roughness, which is needed to accurately quantify the reliability of AM components in service.

36 MATERIALS SCIENCE↗

Tunable phononic bandgap materials designed via topology optimization

Topology optimization is used to design phononic bandgap materials that are tunable by mechanical deformation. A periodic media is considered, which due to the assumption of length scale separation, allows the dispersion relations to be obtained by analyzing a single unit cell subjected to Floquet–Bloch boundary conditions. A finite macroscopic deformation is applied to the unit cell to affect its geometry and hence dispersion. We tune the dispersion–deformation relation to our liking by solving a topology optimization problem using nonlinear programming. The adjoint method is employed to compute the sensitivities, and the non-differentiability of degenerate eigenvalues is avoided using symmetric polynomials. Several tunable phononic crystal designs are presented. Also, a verification analysis is performed, wherein the optimized design is interpreted and analyzed using a conforming finite element mesh.

42 ENGINEERING↗

Chemical order transitions within extended interfacial segregation zones in NbMoTaW

Interfacial segregation and chemical short-range ordering influence the behavior of grain boundaries in complex concentrated alloys. In this study, we use atomistic modeling of a NbMoTaW refractory complex concentrated alloy to provide insight into the interplay between these two phenomena. Hybrid Monte Carlo and molecular dynamics simulations are performed on columnar grain models to identify equilibrium grain boundary structures. Here our results reveal extended near-boundary segregation zones that are much larger than traditional segregation regions, which also exhibit chemical patterning that bridges the interfacial and grain interior regions. Furthermore, structural transitions pertaining to an A2-to-B2 transformation are observed within these extended segregation zones. Both grain size and temperature are found to significantly alter the widths of these regions. An analysis of chemical short-range order indicates that not all pairwise elemental interactions are affected by the presence of a grain boundary equally, as only a subset of elemental clustering types are more likely to reside near certain boundaries. The results emphasize the increased chemical complexity that is associated with near-boundary segregation zones and demonstrate the unique nature of interfacial segregation in complex concentrated alloys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Identification of High-Frequency Transverse Acoustic Modes in Multi-Nozzle Can Combustors

Abstract High frequency thermoacoustic instabilities are problematic for lean-premixed gas turbines. Identifying which acoustic mode is being excited is important, in that it provides insight into potential mitigation measures and mechanical stress/life. However, the frequency spacing between modes becomes quite close for high frequency instabilities in a can combustor. This makes it difficult to distinguish between the modes (e.g., the first transverse mode vs. a higher order axial/mixed mode) based upon frequency calculations alone, which inevitably have uncertainties in boundary conditions, temperature profiles, and combustion response. This paper presents a methodology to simultaneously identify the acoustic mode shapes in the axial and azimuthal directions from acoustic pressure measurements. Multiple high temperature pressure transducers, located at distinct axial and azimuthal positions, are flush mounted in the combustor wall. The measured pressure oscillations from each sensor are then used to reconstruct the pressure distributions by using a least squares method in conjunction with a solution of a three dimensional wave equation. In order to validate the methodology, finite element method (FEM) calculations with estimated post-flame temperature is used to provide the candidate frequencies and corresponding mode shapes. The results demonstrate the reconstructed mode shapes and standing/spinning character of transverse waves, as well as the associated frequencies, both of which are consistent with the FEM predictions. Nodal line location was also extracted from the experimental data during the instabilities in the pressure data.

Kim, J.↗

Magneto-Structural Combined Dimensional and Topology Optimization of Interior Permanent Magnet Synchronous Machine Rotors

A magneto-structural combined dimensional and topology optimization approach for interior permanent magnet synchronous machine (IPMSM) rotor design is proposed using the Solid Isotropic Material Penalization (SIMP) density-based topology optimization method. This method optimizes the location and dimensions of the permanent magnets (PMs) embedded in the rotor, while simultaneously optimizing the lamination flux barriers using topology optimization.To allow for the location and dimensions of the PMs to change while retaining the same number of mesh elements and connectivity, the mesh must be deformed. Mesh deformation is accomplished using a Laplace smoothing technique. To increase the range of PM displacement and size changes, the non-PM areas of the rotor design domain are divided into sub-regions along whose boundaries the mesh nodes are redistributed also using a Laplace smoothing algorithm. Three examples of IPMSM rotor optimizations are presented to demonstrate the proposed approach.

33 ADVANCED PROPULSION SYSTEMS↗

Analysis of the SBP-SAT Stabilization for Finite Element Methods Part I: Linear Problems

In the hyperbolic community, discontinuous Galerkin (DG) approaches are mainly applied when finite element methods are considered. As the name suggested, the DG framework allows a discontinuity at the element interfaces, which seems for many researchers a favorable property in case of hyperbolic balance laws. On the contrary, continuous Galerkin methods appear to be unsuitable for hyperbolic problems and there exists still the perception that continuous Galerkin methods are notoriously unstable. To remedy this issue, stabilization terms are usually added and various formulations can be found in the literature. However, this perception is not true and the stabilization terms are unnecessary, in general. In this paper, we deal with this problem, but present a different approach. We use the boundary conditions to stabilize the scheme following a procedure that are frequently used in the finite difference community. Here, the main idea is to impose the boundary conditions weakly and specific boundary operators are constructed such that they guarantee stability. This approach has already been used in the discontinuous Galerkin framework, but here we apply it with a continuous Galerkin scheme. No internal dissipation is needed even if unstructured grids are used. Further, we point out that we do not need exact integration, it suffices if the quadrature rule and the norm in the differential operator are the same, such that the summation-by-parts property is fulfilled meaning that a discrete Gauss Theorem is valid. This contradicts the perception in the hyperbolic community that stability issues for pure Galerkin scheme exist. In numerical simulations, we verify our theoretical analysis.

97 MATHEMATICS AND COMPUTING↗

Identifying Challenges in Safeguards for Metallic Fuel Fabrication Facilities

As new advanced reactors gain popularity, there is an increasing interest in metallic fuel fabrication for fast reactors. While metallic fuels themselves are not a new idea, as many of the first reactors employed metallic fuels, new designs, compositions, and fabrication methods are appearing throughout the nuclear community. As the interest grows and facilities are constructed, both domestic and international safeguards will need to be heavily involved to support safeguards-by-design (SBD) measures from the start. This work compiles a review of historical and modern fuel types and fabrication methods, fabrication processes, safeguards gaps, and potential safeguards solutions. Metallic nuclear fuel types have been around for many decades and were included in some of the first reactors including the Experimental Breeder Reactor (EBR)-I and -II, the Fermi 1 reactor, the Integral Fast Reactor (IFR), and the Dounreay Fast Reactor (DFR). These reactors used various compositions including pure uranium (U) metal, U-zirconium (Zr) alloys, plutonium (Pu)-aluminum (Al) alloys, U-fissium (Fs) alloys, U-Pu-Zr alloys, and U-molybdenum (Mo) alloys [1, 2, 3, 4, 5]. These small alloying additions are included to improve the material properties of the pure U metal. The alpha-phase U (stable below 661C) suffers elongation in one direction causing grain boundary cracking and increasing creep rate due to irradiation growth, thermal cycling, and preferential crystal orientation. It is ideal to utilize the gamma-phase U (typically stable above 769C) by adding small amounts of alloying elements such as Zr or Mo to stabilize this phase down to room temperature [3]. Additionally, some research has been focused on U with transuranic (TRU) elements present, typically coming from the used fuel recycling process. Including these elements in fast reactor fuel can aid in the reduction of nuclear waste by burning minor long-lived actinides. However, the additions of TRU elements can cause concerns to arise when trying to fabrication or safeguard metallic fuels. A typical metallic fuel element is shown in Figure 1. Sodium is added into the cladding to create a thermal bond between the fuel slug and cladding wall. The fuel slug is then inserted and the end plug is welded on to the top of the fuel element. A gas plenum is left to create a headspace for gaseous fission products to escape rather than continue to build in the fuel itself [1, 5]. Other fuel element geometries exist as well, such as the Lightbridge twisted cruciform geometry shown in Figure 2 [6]. This design allows for better cooling performance and provides room for fuel rod swelling without impacting the fuel rod diameter. There are many different fabrication methods for metallic fuels, which is one of the many benefits of these fuel types. Many of these fabrication methods are relatively easy and cost-efficient. The most popular fabrication method is injection casting, sometimes called vacuum induction melting (VIM), shown in Figure 3 [4, 8, 9, 7, 10]. This method was largely used for EBR-II fuel fabrication. The injection casting system is contained inside of a vessel consisting of a Y2O3-coated graphite crucible surrounded by an induction coil with ZrO2-coated quartz molds suspended above the crucible. The fuel feedstock is placed inside of the graphite crucible and melted using the induction furnace. The induction furnace utilizes a dual frequency with the high frequency melting the feedstock and the low frequency causing stirring of the melted feedstock to form a homogeneous mixture. The mixture is heated to approximately 1600C in an argon environment. The vessel is evacuated and then the quartz molds are lowered into the graphite crucible containing the molten metal and the vessel is repressurized to inject the metal fuel upwards into the molds. The molds are removed and then shattered to release the fuel slugs. This fabrication method was used to fabricate 39,000 metallic fuel pins for EBR-II. While injection casting has been the most common metallic fuel fabrication method throughout the decades, many other methods have been explored including low-pressure gravity casting, microwave casting, continuous casting, centrifugal casting, coextrusion, and many others [11, 12, 8, 13, 14, 15]. Some of these methods aim to mitigate challenges that arise with americium (Am) volatilization during the casting process for TRU-containing fuel feedstocks, an issue with injection casting. Coextrusion is one of the methods explored at the Idaho National Laboratory (INL) and has been utilized for the initial fabrication tests of Lightbridge's unique fuels, as well as other metallic fuels with cladding coextruded. In this process, large billets are formed and machined and then inserted into a molten salt bath for approximately 30 minutes. The billets are then loaded into the extrusion press and extruded. This process can be seen in Figure 4 [15].

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗