SEARCH · Engineering Papers

Results for “Plating”

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 667 records · Page 37

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Three-dimensional full-field velocity measurements in shock compression experiments using stereo digital image correlation

Shock compression plate impact experiments conventionally rely on point-wise velocimetry measurements based on laser-based interferometric techniques. This study presents an experimental methodology to measure the free surface full-field particle velocity in shock compression experiments using high-speed imaging and three-dimensional (3D) digital image correlation (DIC). The experimental setup has a temporal resolution of 100 ns with a spatial resolution varying from 90 to 200 μm/pixel. Experiments were conducted under three different plate impact configurations to measure spatially resolved free surface velocity and validate the experimental technique. First, a normal impact experiment was conducted on polycarbonate to measure the macroscopic full-field normal free surface velocity. Second, an isentropic compression experiment on Y-cut quartz–tungsten carbide assembly is performed to measure the particle velocity for experiments involving ramp compression waves. To explore the capability of the technique in multiaxial loading conditions, a pressure shear plate impact experiment was conducted to measure both the normal and transverse free surface velocities under combined normal and shear loading. The velocities measured in the experiments using digital image correlation are validated against previous data obtained from laser interferometry. Numerical simulations were also performed using established material models to compare and validate the experimental velocity profiles for these different impact configurations. The novel ability of the employed experimental setup to measure full-field free surface velocities with high spatial resolutions in shock compression experiments is demonstrated for the first time in this work.

47 OTHER INSTRUMENTATION↗

DOE OSTI2023

Validation of a detonation product equation of state for an insensitive high explosive via slab geometry expansion tests

Slab expansion (SLABEX) tests are conducted to validate a process for calibrating the detonation products equation of state (EOS) of a high explosive (HE). The SLABEX tests use rectangular slabs of PBX 9502, a polymer-bonded HE formulation consisting of 95 wt.% 1,3,5-triamino-2,4,6-trinitrobenzene bound with Kel F-800, a co-polymer of chlorotrifluoroethylene and vinylidene-fluoride. Three PBX 9502 slab thicknesses are examined, each confined symmetrically by two rectangular copper (Cu) plates approximately one-tenth the thickness of the HE slab. For the duration of each experiment, the detonation flow along the central axis of the PBX 9502 slab remains two-dimensional. The lateral flow velocity component of the outer surfaces of the expanding Cu plates is measured, along with the steady axial detonation speed along the central axis of the SLABEX. Hydrodynamic simulations of the Cu plate expansion in the SLABEX geometry, driven by the energy stored in the detonation products by the detonation combustion event, are conducted using a Jones–Wilkins–Lee EOS for the detonation products. This EOS form was recently parameterized for PBX 9502 in the cylinder expansion test geometry using a newly developed calibration technique [Voelkel et al., Combust. Flame 246, 112373 (2022)]. Good agreement between the experiment and prediction is found in each SLABEX test, demonstrating that the detonation product EOS calibration technique produces EOSs that are predictive when applied to other geometries.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2025

Multifrequency eddy-current detection of fast transient thermal signatures for in situ monitoring applications

Eddy-current (EC) nondestructive evaluation has a long history of use in a variety of ex situ defect monitoring applications because of its exquisite sensitivity to local material variations. Due to the relationship between a material's conductivity and its temperature, EC methods have also been used to investigate quasistatic, long-range temperature variations in casting applications. However, these techniques remain underutilized for the measurement of rapidly varying, spatially nonuniform temperature distributions. In this work, we construct a model system capable of generating repeatable temperature transients in steel plates and measure real-time eddy-current signals with millisecond time resolution and spatial resolution of the order of 1 mm. Using a combination of Multiphysics simulations, fast thermal imaging, and time-resolved holographic interferometry, we tease apart contributions to the eddy-current signals arising from temperature variations and transient plate deformation. Finally, we perform a systematic study in which we vary the plate thickness and the eddy-current excitation frequency to demonstrate that eddy-current techniques can provide information about a three-dimensional, time-varying, subsurface thermal distribution, which is inaccessible to the traditional thermal imaging techniques.

Rosenberg, Ethan R. [Lawrence Livermore National L↗

DOE OSTI2026

Calculating the space-charge-limited current density for nonplanar geometries by simulating the charge-free electric field

Calculating the space-charge-limited-current density (SCLCD) for a complicated diode geometry often requires computationally expensive particle-in-cell (PIC) codes. Here, this paper addresses this issue by using the charge-free electric field $E_0$ calculated using COMSOL Multiphysics to determine local and global SCLCD. The SCLCD obtained by using the surface average of $|E_0|^2$ on the cathode recovers theoretical results for one-dimensional (1D) planar, cylindrical, and tip-to-tip geometries in appropriate limits. We next compared tip-to-tip calculations with the SCLCD obtained using the PIC code Empire. The SCLCD calculated using COMSOL agreed well with Empire for flatter 1D tip-to-tip geometries and diverged with increasing sharpness. Physically, Empire predicts lower SCLCD than COMSOL because the electrons spread due to concentrated space-charge at the tip, whereas theory assumes that the electrons follow the charge-free electric field lines. We further assess the behavior of the SCLCD for tips protruding from the centers of flat, circular plates of various areas. Larger plate areas with constant tip size recover the 1D planar SCLCD globally and 1D tip-to-tip SCLCD locally, while reducing the difference between Empire and COMSOL calculations since larger plates capture more of the emitted electrons, reducing SCLCD suppression due to beam spreading. These results show that charge-free electric field simulations can be used to determine the SCLCD without needing to simulate particle dynamics in PIC.

Wright, Jack K. [Purdue Univ., West Lafayette, IN ↗

DOE OSTI2024

Simulation of plasma and neutral transport in PISCES-RF using SOLPS-ITER

In this research, the fluid plasma transport code SOLPS-ITER is applied and validated against experimental data from the plasma interaction surface component experimental station (PISCES)-RF linear plasma device to establish a physics basis for plasma and neutral transport in its two magnetic field (B-field) geometry setups-(1)the cusp and (2) non-cusp or linear B-field. The main focus of this study is to understand (1) radial plasma transport (2) heat and particle loads on the upstream dump and downstream target plate, and (3) the physics of plasma-neutral interactions in PISCES-RF. The simulation setup adheres to typical PISCES-RF experimental conditions, with a 2D helicon power deposition profile as an input heating source. SOLPS-ITER simulations reproduce experimental conditions with the Bohm diffusion model for both B-field configurations of the PISCES-RF experiment. Major energy loss channels include neutral radiation and power deposited on the wall and dump plate, with only 1% of the input power reaching the target. The ionization front is well confined near the dump plate due to the heating and puffing regions. Additionally, SOLPS-ITER simulation results are also found to be in very good agreement with the particle-in-cell calculations using the code-PICOS++ which supports the validity of SOLPS in low collisionality regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Simulation of the SPARC plasma boundary with the UEDGE code

Here in this work, the UEDGE edge transport code is used to examine conditions in the SPARC divertor and edge plasma for various levels of carbon impurity and power from the core (P SOL ). A double-null magnetic configuration is simulated assuming up-down symmetry in geometry and physics. The anomalous heat and particle transport coefficients are tuned to match empirical predictions for SPARC's midplane density profiles, target plate heat flux profiles, and inner/outer divertor power sharing. Convective transport is included on the low-field side, while on the high-field side the transport is modeled as purely diffusive. Hydrogen neutrals are modeled as a fluid with inertial effects, and a carbon impurity is included using the fixed-fraction model. We find that detachment induced by impurity seeding could significantly reduce the heat flux to the divertor surfaces in the SPARC tokamak. At P SOL = 28 MW (the value predicted for SPARC's full-power H-mode scenario) cases with both divertor legs detached were obtained with a carbon impurity fraction between 0.3%–1.4%, far below Z eff limits for SPARC. When the plasma in the outer leg is detached, the peak heat flux density perpendicular to the target plate is below 1 MW m -2 , electron and ion temperatures are less than 1.5 eV, and momentum detachment is observed. However, the detachment state is found to be sensitive to the side-wall boundary conditions, the level of neutral pumping, and the target plate tilt. Finally, a broadly similar SOLPS simulation of SPARC is used to assess the appropriateness of the simpler impurity and neutral models used in UEDGE.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Plasma–surface interaction in the stellarator W7-X: conclusions drawn from operation with graphite plasma-facing components

W7-X completed its plasma operation in hydrogen with island divertor and inertially cooled test divertor unit (TDU) made of graphite. A substantial set of plasma-facing components (PFCs), including in particular marker target elements, were extracted from the W7-X vessel and analysed post-mortem. The analysis provided key information about underlying plasma–surface interactions (PSI) processes, namely erosion, transport, and deposition as well as fuel retention in the graphite components. The net carbon (C) erosion and deposition distribution on the horizontal target (HT) and vertical target (VT) plates were quantified and related to the plasma time in standard divertor configuration with edge transform ι = 5/5, the dominant magnetic configuration of the two operational phases (OP) with TDU. The operation resulted in integrated high net C erosion rate of 2.8 mg s -1 in OP1.2B over 4809 plasma seconds. Boronisations reduced the net erosion on the HT by about a factor 5.4 with respect to OP1.2A owing to the suppression of oxygen (O). In the case of the VT, high peak net C erosion of 11μm at the strike line was measured during OP1.2B which converts to 2.5 nm s -1 or 1.4 mg s -1 when related to the exposed area of the target plate and the operational time in standard divertor configuration. PSI modelling with ERO2.0 and WallDYN-3D is applied in an interpretative manner and reproduces the net C erosion and deposition pattern at the target plates determined by different post-mortem analysis techniques. This includes also the 13 C tracer deposition from the last experiment of OP1.2B with local 13 CH 4 injection through a magnetic island in one half module. The experimental findings are used to predict the C erosion, transport, and deposition in the next campaigns aiming in long-pulse operation up to 1800 s and utilising the actively cooled carbon-fibre composite (CFC) divertor currently being installed. The CFC divertor has the same geometrical design as the TDU and extrapolation depends mainly on the applied plasma boundary. Extrapolation from campaign averaged information obtained in OP1.2B reveals a net erosion of 7.6 g per 1800 s for a typical W7-X attached divertor plasma in hydrogen.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Controlling the size of non-axisymmetric magnetic footprints using resonant magnetic perturbations

The structure of the non-axisymmetric heat load distribution at the divertor plates is determined not only by the toroidal but also from the poloidal spectrum of non-axisymmetric field perturbations. Whether they are intrinsic, like error fields, or they are applied through 3D coils, the non-axisymmetric fields produce complex 3D edge magnetic topologies (footprints) that alter the properties of the heat and particle flux distributions on the divertor target plates. In this manuscript, a study of the impact of applied 3D field poloidal spectrum on the footprint size and structure is done for the DIII-D tokamak using the resistive MHD code M3D-C1 coupled with the field line tracing code TRIP3D. To resolve the impact of the poloidal spectrum of the magnetic perturbation, the relative phase of the two rows of in-vessel 3D coils used to produce both a n = 2 and a n = 3 perturbation is varied, where n is the toroidal harmonic of the magnetic perturbation. This shows that the largest footprint is predicted when the relative phase of the two rows is close to zero, which is also where the resonant coupling with the plasma is maximized. These results suggest that it will be challenging to decouple the footprint size from the requisite resonant coupling for RMP–ELM control. The correlation between the measured heat load and particle flux distributions at the outer divertor plates in DIII-D and the magnetic measurements is in good agreement with the predicted dependence of the magnetic footprint size on the amplitude of the resonant component of the plasma response.

3D fields↗

DOE OSTI2024

Bayesian parameter estimation and evaluation of the K -ω shear stress transport model for plane impinging jets

Numerical simulations with semi-empirical turbulence models are commonly used to model impinging jets, often used for cooling solid surfaces. In this work, the constants in the k-ω shear stress transport model in ANSYS FLUENT are calibrated to experimental velocity and heat transfer data for a plane turbulent impinging air jet to determine if Kennedy-O'Hagan calibration (Kennedy and O'Hagan 2001 J. R. Stat. Soc. B 63 425–64) can improve predictions of near-surface velocities and surface Nusselt numbers for similar flows. Impinging jets have been proposed to cool the target plates of the divertor in future magnetic fusion energy reactors, where simulations are used to estimate divertor performance. The flat-plate divertor (Wang et al 2009 Fusion Sci. Technol .56 1023–7) uses a plane jet of helium issuing from a B = 0.5 mm slot to cool a surface with radius of curvature of 44 B at a distance 4 B from the slot. Predictions from the calibrated numerical model are compared with independent experimental data at different flow conditions, as well as surface temperature data for a flat plate divertor test section. The contribution of this work is evaluation of the accuracy of a calibrated turbulence model for modest extrapolations in flow geometry and flow conditions for a plane impinging jet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Seismic imaging of the Ecuadorian forearc and arc from joint ambient noise, local, and teleseismic tomography: catching the Nazca slab in the act of flattening

SUMMARY The Ecuadorian Andes are a complex region characterized by accreted oceanic terranes driven by the ongoing subduction of the oceanic Nazca plate beneath South America. Present-day tectonics in Ecuador are linked to the downgoing plate geometry featuring the subduction of the aseismic, oceanic Carnegie Ridge, which is currently entering the trench. Using seismic tomography, we jointly invert arrival times of P and S waves from local and teleseismic earthquakes with surface wave dispersion curves to image the structure of the forearc and magmatic arc of the Ecuadorian Andes. Our data set includes > 100 000 traveltimes recorded at 294 stations across Ecuador. Our images show the basement of the central forearc is composed of accreted oceanic terranes with high elastic wave speeds. Inboard of the Carnegie Ridge, the westernmost forearc and coastal cordilleras display relatively low Vp and Vs and high Vp/Vs values, which we attribute to the increased hydration and fracturing of the overriding plate due to the subduction of the thick oceanic crust of the Carnegie Ridge. We additionally image across-arc differences in magmatic architecture. The frontal volcanic arc overlies accreted terranes and is characterized by low velocities and high Vp/Vs indicative of partial melt reservoirs which are limited to the upper crust. In contrast, the main arc displays regions of partial melt across a wider range of depths. The Subandean zone of Ecuador has two active volcanoes built on continental crust suggesting the arc is expanding eastwards. The mid to lower crust does not show indications of being modified from the magmatic process. We infer that the slab is in the process of flattening as a consequence of early-stage subduction of the buoyant Carnegie Ridge.

58 GEOSCIENCES↗

DOE OSTI2021

Blast from the past: constraining progenitor models of SN 1972E

ABSTRACT We present a novel technique to study Type Ia supernovae (SNe Ia) by constraining surviving companions of historical extragalactic SN by combining archival photographic plates and Hubble Space Telescope(HST) imaging. We demonstrate this technique for Supernova 1972E, the nearest known SN Ia in 125 yr. Some models of SNe Ia describe a white dwarf with a non-degenerate companion that donates enough mass to trigger thermonuclear detonation. Hydrodynamic simulations and stellar evolution models show that these donor stars will survive the explosion, and show increased luminosity for at least a 1000 yr. Thus, late-time observations of the exact location of a supernova can constrain the presence of a surviving donor star and progenitor models. We find the explosion site of SN 1972E by analysing 17 digitized photographic plates taken with the European Southern Observatory 1-m Schmidt and 1 plate taken with the Cerro Tololo Inter-American Observatory 1.5-m telescope. Using the Gaia eDR3 catalogue to determine Supernova 1972E’s location yields: α = 13h39m52${_{.}^{\rm s}}$708 ± 0${_{.}^{\rm s}}$004 and δ = −31°40’9${_{.}^{\prime\prime}}$00 ± 0${_{.}^{\prime\prime}}$04 (ICRS). In 2005, HST/ACS imaged the host galaxy of SN 1972E with the F435W, F555W, and F814W filters covering the explosion site. The nearest detected source is offset by 3.0 times our positional precision, and is inconsistent with the colours expected of a surviving donor star. Thus, the limiting magnitude of the HST observation (F555W > 28 mag) rules out all He star companion models and the most luminous main-sequence companion model currently in the literature. The remaining main-sequence companion models could be tested with a 10 orbit HST exposure in the F606W filter.

Do, Aaron↗

DOE OSTI2020

CRAGE-mediated insertion of fluorescent chromosomal markers for accurate and scalable measurement of co-culture dynamics in Escherichia coli

Monitoring population dynamics in co-culture is necessary in engineering microbial consortia involved in distributed metabolic processes or biosensing applications. However, it remains difficult to measure strain-specific growth dynamics in high-throughput formats. This is especially vexing in plate-based functional screens leveraging whole-cell biosensors to detect specific metabolic signals. Here, we develop an experimental high-throughput co-culture system to measure and model the relationship between fluorescence and cell abundance, combining chassis-independent recombinase-assisted genome engineering (CRAGE) and whole-cell biosensing with a PemrR-green fluorescent protein (GFP) monoaromatic reporter used in plate-based functional screening. CRAGE was used to construct Escherichia coli EPI300 strains constitutively expressing red fluorescent protein (RFP) and the relationship between RFP expression and optical density (OD600) was determined throughout the EPI300 growth cycle. A linear equation describing the increase of normalized RFP fluorescence during deceleration phase was derived and used to predict biosensor strain dynamics in co-culture. Measured and predicted values were compared using flow cytometric detection methods. Induction of the biosensor lead to increased GFP fluorescence normalized to biosensor cell abundance, as expected, but a significant decrease in relative abundance of the biosensor strain in co-culture and a decrease in bulk GFP fluorescence. Taken together, these results highlight sensitivity of population dynamics to variations in metabolic activity in co-culture and the potential effect of these dynamics on the performance of functional screens in plate-based formats. The engineered strains and model used to evaluate these dynamics provide a framework for optimizing growth of synthetic co-cultures used in screening, testing and pathway engineering applications.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Computational and Experimental Study of Film-Cooling Effectiveness With and Without Downstream Vortex Generators

Downstream vortex generators that involve a pair of rectangular plates arranged in an open V-shape placed just downstream of each film-cooling hole were shown to create flow and vortical structures that entrain lifted film-cooling flow back to the surface and increase its lateral spreading on the surface (GT2020–14317). In this study, computations and measurements were performed to examine the flow mechanisms this vortex generator induces to improve film-cooling effectiveness of a flat plate with the cooling flow emanating from one row of inclined holes. Parameters studied include blowing ratio (BR = 0.75 and 1.0) and temperature ratio (TR = 1.07 and 1.9). Here, the computational study is based on steady Reynolds-averaged Navier–Stokes (RANS) closed by the shear-stress transport (SST) turbulence model with and without conjugate analysis. The experimental study was conducted by using a conjugate heat transfer test rig with a plenum, where cooling flow is introduced. Measurements made include velocity and temperature profiles upstream and downstream of the film-cooling holes as well as the temperature at several locations on the hot and cold sides of the film-cooled flat plate. The computational study was validated by comparing computed results with those from measurements at BR = 0.75 and 1.0 and TR = 1.9. Computational and experimental results are presented to show the effects of BR and TR on the flow structures and how those structures improve the effectiveness of film cooling with and without the downstream vortex generators and with and without conjugate heat transfer.

42 ENGINEERING↗

DOE OSTI2024

Effect of Substrate Temperature on the Electrochemical and Supercapacitance Properties of Pulsed Laser-Deposited Titanium Oxynitride Thin Films

Electrocatalytically active titanium oxynitride (TiNO) thin films were fabricated on commercially available titanium metal plates using a pulsed laser deposition method for energy storage applications. The elemental composition and nature of bonding were analyzed using X-ray photoelectron spectroscopy (XPS) to reveal the reacting species and active sites responsible for the enhanced electrochemical performance of the TiNO electrodes. Symmetric supercapacitor devices were fabricated using two TiNO working electrodes separated by an ion-transporting layer to analyze their real-time performance. The galvanostatic charge–discharge studies on the symmetric cell have indicated that TiNO films deposited on the polycrystalline titanium plates at lower temperatures are superior to TiNO films deposited at higher temperatures in terms of storage characteristics. For example, TiNO films deposited at 300 °C exhibited the highest specific capacity of 69 mF/cm 2 at 0.125 mA/cm 2 with an energy density of 7.5 Wh/cm 2 . The performance of this supercapacitor (300 °C TiNO) device is also found to be ~22% better compared to that of a 500 °C TiNO supercapacitor with a capacitance retention ability of 90% after 1000 cycles. Furthermore, the difference in the electrochemical storage and capacitance properties is attributed to the reduced leaching away of oxygen from the TiNO films by the Ti plate at lower deposition temperatures, leading to higher oxygen content in the TiNO films and, consequently, a high redox activity at the electrode/electrolyte interface.

30 DIRECT ENERGY CONVERSION↗

DOE OSTI2024

Iterative Stress Reconstruction Algorithm to Estimate Three-Dimensional Residual Stress Fields in Manufactured Components

Residual stress (RS) significantly impacts the mechanical performance of components. Measurement of RS often provides incomplete data in terms of components of stress and spatial density. Employing such fields in finite element simulations results in significant modification of the field to achieve equilibrium and compatibility among strains. To overcome this, an iterative stress reconstruction algorithm (ISRA) is developed to estimate 3D RS fields that satisfy equilibrium, are stress component-wise complete, and represent the characterized data sampled. An Al 7075-T651 plate and an additively manufactured (AM) A36 steel wall are considered for RS reconstruction using measurement data from the literature. A maximum variation of ~2.5 MPa in the Al plate, and ~10 MPa in the steel wall are observed between the reconstructed and measured stresses. Furthermore, unknown stress components emerge and reach significant magnitudes (upto ~2.3 MPa in the Al plate and ~45 MPa in the AM wall) during ISRA. Indeed, it is found that minor errors in measurement or data processing are eliminated through the physical requirements during ISRA. Employing a reconstructed RS field is hence not just more accurate given its compatibility, but it additionally corrects for minor errors in measurement. Furthermore, it is found that spatially dense measurement data result in convergence with fewer iterations. Finally, although ISRA yields a nonunique solution dependent on boundary conditions, measurement errors, fitting errors, and mesh density, it accommodates for uncertainties and inaccuracies in measurement, as opposed to failing to reach a physically realistic converged solution.

42 ENGINEERING↗

DOE OSTI2022

ORGANIC RANKINE CYCLE TURBINE AND HEAT EXCHANGER SIZING FOR LIQUID AIR COMBINED CYCLE

Cryogenic energy storage offers several opportunities to design turbomachinery and other equipment for novel cycles. This paper presents the design and analysis of turbomachinery and heat exchangers for an Organic Rankine Cycle (ORC) subsystem for a hybrid energy storage concept. The Liquid Air Combined Cycle is an energy storage system that stores air at cryogenic conditions at times with high variable renewable energy to be dispatched along with a gas turbine to recover the exhaust heat. In order to re-vaporize the air, the liquid air is coupled with an ORC as an additional bottoming cycle. The ORC turbine is expected to expand the fluid with a pressure ratio of nearly 30 and a flow rate of approximately 45 kg/s. Sizing calculations for both a radial and axial turbine solution were performed over a range of speeds and stages to determine the optimal design point. The results show that either an axial (8- or 9-stage) or radial (four stages at two shaft speeds) turbine are capable of handling the pressure ratios. Further trades of the two configurations would be required to determine the best option. The ORC system also incorporates five heat exchangers to distribute heat, vaporize the liquid air, or recover exhaust heat from the gas turbine. Three heat exchangers were analyzed to understand the size of heat exchangers and pressure drop for the overall system. Different types of heat exchangers were explored for the different purposes, including plate-fin heat exchangers, gasketed plate heat exchangers and shell-in-tube heat exchangers. It was determined that the ORC recuperator, liquid-air vaporizer, and vaporized air pre-heater would be counter-flow heat exchangers using a gasketed plate design. Keywords: Energy Storage, Liquid Air Energy Storage, Organic Rankine Cycle

Pryor, Owen↗

DOE OSTI2022

Computational Modeling of Multi-Pass Rolling Parameters Effect on Resulting Fuel Foil Shape

A focus of the U.S. Department of Energy is to improve production yield and reduce the cost of Low Enriched Uranium (LEU)-molybdenum alloy (U-10Mo) monolithic fuel plates that will be replacing High Enriched Uranium (HEU) oxide dispersion fuels used currently in the United States High Performance Research Reactors (USHPRR). One area of improvement is lowering the high transverse waviness and longitudinal waviness currently present within rolled foils prior to cladding to produce fuel plates. Traditional rolling manufacturing techniques for other metal foils use winders to pull and straighten the foil as it is rolled back and forth to the final thickness. This approach cannot be used to roll thin U-10Mo foils (0.008-0.025” thick) because only small castings can be rolled due to nuclear criticality safety concerns. As a result, the fuel foils are too short (1 m in length) to use traditional winders. Therefore, it is crucial to identify other rolling parameters (i.e., roller friction, axial tension load, roller diameter, and roll pass reduction percent) that might reduce transverse waviness and longitudinal waviness in the rolled fuel foil and develop a high-yield, low-cost multi-pass rolling manufacturing process. This report documents a systematic finite element modeling study to investigate the effects of numerous rolling parameters to reduce resulting transverse waviness and longitudinal waviness in the fuel foil during multi-pass rolling of U-10Mo foils. The rolling of a U-10Mo plate with initial dimensions of 1”x1”x 0.048” is modeled using Abaqus CAE. This rolling is modeled to undergo eight 20% reduction roll passes to a final fuel foil thickness of 0.01”. The elastic-plastic constitutive model of the U-10Mo alloy was input to the fuel foil rolling model. The rollers were modeled as rigid bodies. A comparison of rolling friction coefficients of 0.3 and 0.7 over a wide range of applied axial tension loads were investigated in order to evaluate the effect of using a lubricant during rolling. The effect of roller diameter on the resulting transverse waviness and longitudinal waviness of the fuel foil over a wide range of axial tension loads were also investigated by modeling rollers 7/8” and 3.75” in diameter. The results of this systematic finite element method study will aid manufacturers in producing low transverse waviness and reduced longitudinal waviness in U-10Mo fuel foils.

U-10Mo, FEA, Rolling, Residual Stress, Fuel Foil↗

DOE OSTI2020

Mechano-Electrochemical Interaction in Solid-State Lithium Batteries

A solid-state lithium (Li) battery primarily consists of Li metal anode, solid electrolyte separator, and cathode. The asymmetric volume changes, originating from ion transport and interfacial Li growth during plating, lead to significant stresses in the layered architecture. In this study, we develop a coupled mechanics-electrochemistry formalism for polymer electrolyte based solid-state batteries, in particular, focusing on the stress effect on electrochemical performance. By means of a coupling coefficient, it is found that stress-assisted ion transport in the electrolyte results in a delayed Sand’s time and increased critical current density of unstable electrodeposition, and consequently, alleviates the propensity of dendrite formation. Stress at the Li metal-electrolyte interface affects the electrochemical reaction kinetics, and the influences from the deviatoric stress and hydrostatic pressure vary with Li plating time. In addition, a low restraint stiffness to the layered structure could elastically buffer the volumetric changes and thus reduce the stress during Li plating. This fundamental study provides guidance for the design of solid-state batteries, aimed at stable electrodeposition and mechanical integrity.

25 ENERGY STORAGE↗