Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TiP”

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 775 records · Page 43

Development of Solid Synchronous Reluctance Rotors With Multi-Material Additive Manufacturing

Synchronous reluctance (SynR) machines are promising rare-earth material-free alternatives to permanent magnet machines. However, structural challenges limit their operating speed and power density. This paper proposes and investigates multi-material additive manufacturing (MMAM) as a key-enabler to realize power-dense and high-speed SynR machines. It does so by proposing designs that guide magnetic flux through solid rotors realized by selective placement of magnetic and non-magnetic materials. To explore this concept, first, material samples are additively manufactured and experimentally characterized to assess the structural and magnetic properties that can be expected for the proposed rotors. Second, the design space of each rotor type is explored using these measured properties within finite element analysis. The results reveal that MMAM can enable fabrication of SynR motors with power density levels that are at the leading edge of all conventional electric machine topologies. It is shown that tip speeds in excess of 300 m/s can be achieved, resulting in 3-4x improvement in power density over conventional SynR motors. A solid SynR rotor is printed in an experimental MMAM laser powder bed fusion system. The rotor is paired with an existing stator to create a functional SynR motor with a saliency ratio of 2.59 and torque rating of 4.15 Nm. This is the first publication of a SynR rotor prototype constructed via MMAM.

36 MATERIALS SCIENCE↗

Modeling of Magnetic Properties of Magnetorheological Elastomers Using JA Hysteresis Model

Magnetorheological elastomers (MREs) are composite materials that consist of magnetically permeable particles in a nonmagnetic polymeric matrix. Under the influence of an external magnetic field, a reversible deformation change occurs in the mechanical properties of these materials. Due to their coupled magnetomechanical response, these materials have been found suitable for a range of applications including tunable vibration absorbers, sensors, and actuators. Notably, improvement of such devices are prerequisites to efficient energy conversion systems, hence the need to understand further the MRE technology. The Jiles-Atherton (JA) theory takes into consideration the magneto-coupling experienced by effective domains in a magnetic material. Algorithm based on the theory yields five model parameters; saturation magnetization (M s ), domain density (a), domain coupling (α), loss coefficient (k), and reversibility (c). Using JA theory, model parameters were calculated and linked to the physical attributes of Fe powder and isotropic MRE. The results show that the calculated Ms for the MRE is reasonably related to that of the Fe powder by a factor of the particle's volume fraction used in the MRE. Further, the calculated k, a, and α provided support for the reduced pinning factor, domain density, and increased domain coupling in the MRE due to the changes in the domain structure between the two materials. From the calculated JA parameters, finite-element modeling (FEM) of the MRE hysteresis loop was performed. The analysis showed that the modeled magnetic properties including coercivity, remanence, and coordinates of the hysteresis loop tip vary with geometric position.

42 ENGINEERING↗

Design and Testing of a Prototype Eddy Current Actuated Valve for the ITER Shattered Pellet Injection System

Reliably mitigating disruptions is essential for ITER to meet its long-term operational research plan without damage to the in-vessel components. Currently, the shattered pellet injection (SPI) technique is the most effective radiator of thermal energy and has been chosen for the baseline disruption mitigation system (DMS) for ITER. The SPI process uses cryogenic temperatures to desublimate material into the barrel of a pipe gun forming a solid cylindrical pellet. Pellets for ITER will initially be hydrogen and hydrogen-neon mixtures. Once formed, pellets are dislodged and accelerated using high-pressure gas (40-60 bar) delivered by a fast-opening valve. The solenoid valves currently used for SPI experiments will not operate in an ITER environment due to the large background magnetic field. An ITER prototype fast-opening valve, called a flyer plate valve (FPV), has been designed and has undergone a wide range of testing. The FPV operates by pulsing current through a pancake coil that is closely coupled with a ``flyer plate.'' The flyer plate is an aluminum plate in which eddy currents are generated creating a repulsive force from the pancake coil. The force generated in the flyer plate rapidly lifts the valve tip off the seat and delivers a pulse of gas to the rear of the pellet, breaking it free from the barrel and accelerating the pellet downstream to its intended target. The design of the valve has been iterated on over the lifetime of this project, as the DMS for ITER shifted from massive gas injection (MGI) to SPI. The most recent design has been tested, and operational ranges have been mapped. The valve must survive 3000+ cycles in an ITER-like magnetic field. The principal functional requirement of this valve is to reliably dislodge and accelerate hydrogen (or H-Ne mixture) pellets into ITER. The valve was mated with an ITER SPI test stand and has been shown to be capable of launching pellets reliably. The valve and power supply design will be discussed in this article, along with the various testing setups used to determine the feasibility of this valve for use on ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental and computational studies of load ratio effect on fatigue crack growth rates in Alloy 709 at elevated temperatures

Abstract Experimental testing and finite element simulations of fatigue crack growth were performed in austenitic 20Cr‐25Ni (Alloy 709) steel at different load ratios and elevated temperatures. The experimental tests were performed using compact tension specimens, and crack growth rates were measured at stress intensity factors ranging between 5 and 35 MPa√m. Fractographic analysis using scanning electron microscopy indicated crack surface roughness and secondary cracking depending on testing temperature and load ratio. Finite element simulations of fatigue crack growth were performed to compute plasticity‐induced crack opening loads and predict crack growth rates. Predictions of fatigue crack growth rates using finite element simulations were performed using the computed crack‐tip opening loads, and they are shown to match well the experimental measurements.

Engineering↗

Climate regime shift and forest loss amplify fire in Amazonian forests

Frequent Amazonian fires over the last decade have raised the alarm about the fate of the Earth's most biodiverse forest. The increased fire frequency has been attributed to altered hydrological cycles. However, observations over the past few decades have demonstrated hydrological changes that may have opposing impacts on fire, including higher basin-wide precipitation and increased drought frequency and severity. Here, we use multiple satellite observations and climate reanalysis datasets to demonstrate compelling evidence of increased fire susceptibility in response to climate regime shifts across Amazonia. In this work, we show that accumulated forest loss since 2000 warmed and dried the lower atmosphere, which reduced moisture recycling and resulted in increased drought extent and severity, and subsequent fire. Extremely dry and wet events accompanied with hot days have been more frequent in Amazonia due to climate shift and forest loss. Simultaneously, intensified water vapor transport from the tropical Pacific and Atlantic increased high-altitude atmospheric humidity and heavy rainfall events, but those events did not alleviate severe and long-lasting droughts. Amazonia fire risk is most significant in the southeastern region where tropical savannas undergo long seasonally dry periods. We also find that fires have been expanding through the wet-dry transition season and northward to savanna-forest transition and tropical seasonal forest regions in response to increased forest loss at the "Arc of Deforestation." Tropical forests, which have adapted to historically moist conditions, are less resilient and easily tip into an alternative state. Our results imply forest conservation and fire protection options to reduce the stress from positive feedback between forest loss, climate change, and fire.

54 ENVIRONMENTAL SCIENCES↗

From Depletion to Restoration: Lessons From Long‐Term Monitoring of Carbon Gains and Losses in Cropping Systems

As global atmospheric CO 2 rapidly approaches a key tipping point, there is an urgent need to implement strategies to reverse this pattern. A generally accepted understanding of carbon (C) in agricultural fields includes: (H1) substantial C loss occurs when natural vegetation is converted to crops, (H2) soils typically reach a steady-state C concentration under contemporary practices, and (H3) improved management or crop selection can enhance soil C stocks over time. Significant variability exists, but studies consistently show large C losses from agricultural ecosystems, supporting H1. Although steady-state C levels (H2) are commonly assumed, measuring C gains or losses in mature agroecosystems is challenging. Efforts to increase soil C storage (H3) have limited data due to the diversity of potential practices, compounded by substantial variability in soil C measurements. Here, long-term (7–17 year) ecosystem C flux data from diverse cropping systems revealed that conventionally tilled annual row crops (maize and soybean) act as significant long-term atmospheric C sources, challenging H2. Furthermore, conservation tillage practices reduced C losses compared with conventional tillage but showed minimal evidence for long-term ecosystem C storage, even after 20+ years. This indicates that no-till practices reduce C losses but imply that no soil C is added, challenging H3. By contrast, perennial Miscanthus × giganteus, Panicum virgatum, and restored tallgrass prairie systems store C at the ecosystem scale more effectively than minimally tilled annual row crops. Analysis over multiple years demonstrates significant ecosystem C storage with perennial crops, varying by species, starting in the first year of transition. These findings, although focused on one region, suggest that the assumptions of steady-state C levels and increased storage from conservation practices do not universally apply and that significant changes to agroecosystems are required to increase C storage.

59 BASIC BIOLOGICAL SCIENCES↗

Structure-mechanical properties correlation in bulk LiPON glass produced by nitridation of metaphosphate melts

The glassy solid electrolyte Lithium phosphorous oxynitride (LiPON) has been widely researched in thin film solid state battery format due to its outstanding stability when cycled against lithium. In addition, recent reports show thin film LiPON having interesting mechanical behaviors, especially its ability to resist micro-scale cracking via densification and shear flow. Here, in the present study, we have produced bulk LiPON glasses with varying nitrogen contents by ammonolysis of LiPO 3 melts. The resulting compositions were determined to be LiPO 3-3z/2 N z , where 0 ≤ z ≤ 0.75, and the z value of 0.75 is among the highest ever reported for this series of LiPON glasses. The short-range order structures of the different resulting compositions were characterized by infrared, Raman, 31 P magic angle spinning nuclear magnetic resonance, and X-ray photoelectron spectroscopies. Instrumented nano-indentation was used to measure mechanical properties. It was observed that similar to previous studies, both trigonally coordinated (N t ) and doubly bonded (N d ) N co-exist in the glasses in about the same amounts for z ≤ 0.36, the limit of N content in most previous studies. For glasses with z > 0.36, it was found that the fraction of the N t increased significantly while the fraction of N d correspondingly decreased. The incorporation of nitrogen increased both the elastic modulus and hardness of the glass by approximately a factor of 1.5 when N/P ratio reaches 0.75. At the same time, an apparent embrittlement of the glass was observed due to nitridation, which was revealed by nanoindentation with an extra sharp nanoindenter tip.

25 ENERGY STORAGE↗

BUZZ : an essential gene for postinitiation root hair growth and a mediator of root architecture in Brachypodium distachyon

Here, we discover a player in root development. Recovered from a forward-genetic screen in Brachypodium distachyon, the buzz mutant initiates root hairs but they fail to elongate. In addition, buzz roots grow twice as fast as wild-type roots. Also, lateral roots show increased sensitivity to nitrate, whereas primary roots are less sensitive to nitrate. Using whole-genome resequencing, we identified the causal single nucleotide polymorphism as occurring in a conserved but previously uncharacterized cyclin-dependent kinase (CDK)-like gene. The buzz mutant phenotypes are rescued by the wild-type B. distachyon BUZZ coding sequence and by an apparent homolog in Arabidopsis thaliana. Moreover, T-DNA mutants in A. thaliana BUZZ have shorter root hairs. BUZZ mRNA localizes to epidermal cells and develops root hairs and, in the latter, partially colocalizes with the NRT1.1A nitrate transporter. Based on qPCR and RNA-Seq, buzz overexpresses ROOT HAIRLESS LIKE SIX-1 and -2 and misregulates genes related to hormone signaling, RNA processing, cytoskeletal, and cell wall organization, and to the assimilation of nitrate. Overall, these data demonstrate that BUZZ is required for tip growth after root hair initiation and root architectural responses to nitrate.

59 BASIC BIOLOGICAL SCIENCES↗

How deep should we go to understand roots at the top of the world?

Harsh environmental conditions and the short summers of northern, high-latitude biomes impose unique constraints on the plants that live in the arctic tundra and the boreal forest. To escape the harsh aboveground environment, plants in these habitats often allocate a large portion of their biomass belowground to facilitate nutrient acquisition. In turn, the proximity of living plant roots to vast stores of sequestered soil carbon in these biomes means that shifts in rooting depth distribution and the size of the root–soil interface could significantly contribute to ongoing climate change. Indeed, plant ‘priming’ of rhizosphere decomposition via root exudation, particularly from shallowly distributed roots, can lead to losses of carbon from tundra soils. While we have a hard-won understanding of the distribution of plant communities across the arctic tundra and the boreal forest from direct field observations scaled to the landscape level using climate-informed mapping techniques (i.e. the Circumpolar Arctic Vegetation Map (CAVM); Walker et al. , 2005), these vegetation maps are only the tip of the iceberg (Iversen et al. , 2015). Root form and function remain hidden beneath the land surface. In an article recently published in New Phytologist , Blume-Werry et al. (2023, 10.1111/nph.18998) asked whether rooting depth distribution, and ensuing carbon emissions, could be inferred from commonly used vegetation mapping classifications across the pan-Arctic. An important question to guide our understanding, mapping, and prediction of belowground characteristics and ecosystem feedbacks at the top of the world. Unfortunately, they found that the answer was ‘not quite’. While rooting depth distribution varied demonstrably, in turn causing substantial changes in modeled carbon emissions via rhizosphere ‘priming’, variation across rooting depth profiles did not correspond with vegetation mapping classes. If we are unable to predict belowground rooting depth distributions across large spatial scales by leveraging aboveground vegetation community distributions, how then should belowground researchers proceed?

59 BASIC BIOLOGICAL SCIENCES↗

The UPR regulator IRE1 promotes balanced organ development by restricting TOR-dependent control of cellular differentiation in Arabidopsis

Proteostasis of the endoplasmic reticulum (ER) is controlled by sophisticated signaling pathways that are collectively called the unfolded protein response (UPR) and are initiated by specialized ER membrane-associated sensors. The evidence that complete loss-of-function mutations of the most conserved of the UPR sensors, inositol-requiring enzyme 1 (IRE1), dysregulates tissue growth and development in metazoans and plants raises the fundamental question as to how IRE1 is connected to organismal growth. To address this question, we interrogated the Arabidopsis primary root, an established model for organ development, using the tractable Arabidopsis IRE1 mutant ire1a ire1b, which has marked root development defects in the absence of exogenous stress. We demonstrate that IRE1 is required to reach maximum rates of cell elongation and root growth. We also established that in the actively growing ire1a ire1b mutant root tips the Target of Rapamycin (TOR) kinase, a widely conserved pro-growth regulator, is hyperactive, and that, unlike cell proliferation, the rate of cell differentiation is enhanced in ire1a ire1b in a TOR-dependent manner. By functionally connecting two essential growth regulators, these results underpin a novel and critical role of IRE1 in organ development and indicate that, as cells exit an undifferentiated state, IRE1 is required to monitor TOR activity to balance cell expansion and maturation during organ biogenesis.

59 BASIC BIOLOGICAL SCIENCES↗

Energy Release Rate Approximation for Small Surface Cracks in Three-Dimensional Domains Using the Topological Derivative

The topological derivative describes the variation of a response functional with respect to infinitesimal changes in topology, such as the introduction of an infinitesimal crack or hole. In this three-dimensional fracture mechanics work, we propose an approximation of the energy release rate field associated with a small surface crack of any boundary location, direction, and orientation combination using the topological derivative. Furthermore, this work builds on the work of Silva et al. (“Energy Release Rate Approximation for Small Surface-Breaking Cracks Using the Topological Derivative,” J. Mech. Phys. Solids 59(5), pp. 925–939), in which the authors proposed an approximation of the energy release rate field which was limited to two-dimensional domains. The proposed method is computationally advantageous because it only requires a single analysis. By contrast, current boundary element and finite element-based methods require an analysis for each crack length-location-direction-orientation combination. Furthermore, the proposed method is evaluated on the non-cracked domain, obviating the need for refined meshes in the crack tip region.

42 ENGINEERING↗

Additively Manufactured Compliant Hybrid Gas Thrust Bearing for Supercritical Carbon Dioxide Turbomachinery: Design and Proof of Concept Testing

The following paper presents a new type of gas lubricated thrust bearing that utilizes additive manufacturing or also known as direct metal laser melting (DMLM) to fabricate the bearing. The motivation for the new bearing concept is derived from the need for highly efficient supercritical carbon dioxide (sCO2) turbomachinery in the mega-watt power range. The paper provides a review of existing gas thrust-bearing technology, outlines the need for the new DMLM concept, and discusses proof-of-concept testing results. The new concept combines hydrostatic pressurization with individual tilting pads that are flexibly mounted using hermetic squeeze film dampers (HSFD) in the bearing-pad support. This paper describes the thrust-bearing concept and discusses the final design approach. Proof-of-concept testing in air for a 6.8 in. (173 mm) outer diameter thrust gas bearing was performed; with thrust loading, up to 1500 lbs (6.67 kN) and a thrust runner speed of 10krpm (91 m/s tip speed). The experiments were performed with a bent shaft resulting in thrust runner axial vibration magnitudes of 2.9 mils (74 μm) p-p and dynamic thrust loads of 270 lbs (1.2 kN) p-p. In addition, force deflection characteristics and stiffness coefficients of the bearing system are presented for an inlet hydrostatic pressure of 380 psi (2.62 MPa). Results at 10 krpm show that the pad support architecture was able to sustain high levels of dynamic misalignment equaling 6 times the nominal film clearance while demonstrating a unit load-carrying capacity of 55 psi (0.34 Mpa). Gas-film force deflection tests portrayed nonlinear behavior like a hardening spring, while the bearing pad support stiffness was measured to be linear and independent of gas film thickness.

Engineering↗

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗

A Phenomenological Thermal Spray Wall Interaction Modeling Framework Applied to a High Temperature Ignition Assistant Device

Airborne compression ignition engines must operate with reliable ignition systems to achieve proper ignition at every cycle, particularly at high altitudes. Glow-plug-based ignition-assistant (IA) devices can provide the necessary energy to preheat the fuel and ensure ignitability of the fuel-air mixture. Ignitability of liquid sprays can be facilitated via direct impingement onto the hot IA surface, however this comes with adverse effects on the IA durability. Therefore, optimizing an IA’s design requires detailed understanding of the physics of fuel spray impingement of superheated surfaces. While spray impingement on relatively low wall temperatures has been extensively studied and appropriate numerical models have been proposed through the years, fundamental understanding of high-speed liquid spray impingement on superheated walls is still elusive. This work aims to formulate a phenomenological thermal spray-wall interaction framework for modeling the film-boiling-induced heat transfer, atomization, and dispersion of fuel spray droplets impinging onto a superheated IA device. A qualitative comparison of the new phenomenological model is performed against optical experiments from the literature of an F-24 fuel spray injected onto an IA device located 12 mm away from the injector tip. The temperature of the IA was set at 1400 K. The fuel injection pressure was 400 bar, while the ambient gas pressure and temperature were 30 bar and 800 K, respectively. The performance of the phenomenological model is evaluated in comparison with two other state-of-art models from the literature. A qualitative analysis of the different spray and fuel-air mixture characteristics is performed to outline the differences in the predictions offered by the new phenomenological model and the two state-of-art spray-wall interaction models.

droplet dispersion↗

Numerical Investigation of Ammonia/n-Heptane Dual-Fuel Spray Flames Using Large Eddy Simulations

Diesel engines are extensively used in heavy-duty transportation, power generation, and marine vehicles due to their superior thermal efficiency and extended high-load operability compared to spark ignition (SI) engines. However, combustion in diesel engines is generally characterized by locally rich fuel–air mixtures and high combustion temperatures, causing significant amounts of soot and NO x emissions from these engines. Utilizing carbon-free alternative fuels and enhancing fuel efficiency represent promising strategies to mitigate greenhouse gas (GHG) and other emissions in the heavy-duty transportation sector. In this context, ammonia (NH 3 ), as a hydrogen carrier, has received significant attention as a viable substitute for hydrocarbon fuels due to its carbon-free composition, relatively high energy density, and well-established infrastructure. Many previous studies have considered combustion and emission characteristics of ammonia-hydrocarbon fuel blends in engines and simplified flames. But, detailed investigations on the effects of ammonia on the performance of hydrocarbon fuels under engine conditions are lacking. In the present study, we perform large eddy simulations (LES) of the ignition and flame processes in a constant-volume combustion reactor, where n-heptane is injected in an ammonia/air ambient mixture in a diesel-like environment. A detailed and validated reaction mechanism containing 302 species and 1981 reactions is employed. The Engine Combustion Network Spray H experimental data is used to validate the spray model under both non-reacting and reacting conditions. Dual-fuel combustion is simulated using the well-stirred reactor (WSR) approach. Results are presented for two spray cases: (1) single fuel (SF) with n-heptane injected into a mixture of air and combustion products and (2) dual-fuel (DF) with the injection of n-heptane in a mixture of air, ammonia, and combustion products. It is observed that the presence of ammonia has a significant effect on the ignition and flame development processes. With ammonia addition, both the first- and second-stage ignition delay times increase, but the effect of ammonia on the second-stage ignition is significantly more prominent. In addition, the ignition kernel size and growth rate decrease noticeably. For SF spray, the main ignition is characterized by multiple ignition kernels near the spray tip, whereas for DF spray, a single relatively small ignition kernel forms and grows slowly in the downstream direction. The flame development and the final quasi-steady flame structure are also modified due to ammonia. Here, the outcome of this research would enable a better understanding of ammonia–diesel dual-fuel spray flame behavior and guide the development of associated engine combustion strategies.

ammonia↗

SYSTEMATIC MULTI-INDEX VALIDATIONS OF SIDI ENGINE FLOW FIELD LES COMPUTATIONS USING CRANK ANGLE-RESOLVED PIV MEASUREMENTS

In-cylinder flow fields make significant impacts on the fuel atomization, fuel mixture formation, and combustion process in spark ignition direct injection (SIDI) engines. In recent years, model-based simulation approaches are preferred in regard to investigating the transient in-cylinder flow field characteristics. Most commonly, the simulation models are validated using single representative flow field at a typical crank angle measured by particle image velocimetry (PIV). However, it provides only limited knowledge about the flow field which is highly three-dimensional and of transient nature. In this study, crank angle-resolved Ply measurements are conducted on three distinct planes inside the cylinder to capture the transient process of flow field characteristics which vary with the crank angle. These three planes consist of one tumble plane through the spark plug tip, one tumble plane along two intake ports, and one swirl plane at 30 mm below the cylinder head. Large eddy simulation (LES) is employed for the numerical computations using the CONVERGE codes. On the basis of large datasets for both temporal and spatial domains, a multi-index systematic validation approach is conducted. Crank angle-resolved calculations of global indices and local indices are implemented using the flow fields velocity data obtained from both PIV and LES on select planes. Global indices reveal the trends in similarities of different crank angle degrees and locations, while local indices give the detail comparison results. In summary, with the systematic multi-index validation approach, the crucial crank angle degrees and locations for model verification will be detected. Furthermore, the corresponding critical flow features are analyzed. Practical guideline of flow field validation is proposed.

Crank angle-resolved PIV↗

AN ESTIMATE OF SPENT NUCLEAR FUEL MECHANICAL LOADS IN THE GENERAL 30 CM PACKAGE DROP SCENARIO

The US Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is performing research to determine the mechanical loading conditions applied to spent nuclear fuel (SNF) during normal conditions of transport to inform mechanical tests of SNF and close an important knowledge gap related to the practical disposition of SNF in the US. Researchers at Pacific Northwest National Laboratory (PNNL) have completed an extensive finite element study to characterize and estimate the potential mechanical loads on SNF during a hypothetical 30 cm drop of an SNF transportation package. This modeling study is validated with test data collected by the SFWST program during a physical test campaign that included one-third scale package drop tests and full-scale single fuel assembly drop tests. The test campaign was led by Sandia National Laboratories (SNL) and included international collaboration with Equipos Nucleares S.A, S.M.E (ENSA) and Bundesanstalt für Materialforschung und -prüfung (BAM). The key contribution of the modeling study is to go beyond the limitations of the limited number of physical tests to estimate the impact response to variations in impact angle, initial gap conditions, fuel assembly design, burnup and other parameters that affect the mechanical loads. The methodology of this study included a classic parametric study to calculate the impact response of highly detailed fuel assemblies over many combinations of parameters. Models of a 17x17 pressurized water reactor fuel assembly and a generic 10x10 boiling water reactor fuel assembly were both used in this study to cover the major fuel assembly types in the US inventory. Over 2,000 impact responses were calculated. The results of the parametric study were evaluated using traditional methods and basic statistics. The results were also used to construct a damage model using multiple nonlinear regression techniques to predict the mechanical loads over the full range of all input parameters. The damage model was found to work very well for all impact angle cases where the cask came to rest on its side. It was concluded that end drop cases where the cask remained vertical (instead of tipping over onto its side) were not sufficiently characterized by the current set of parametric study cases to include in the damage model, but it was not a priority to fully investigate that range because the highest mechanical loads were observed in the broader range of side impact cases. This modeling work provides sufficient insight into the mechanical loads on SNF during a hypothetical 30 cm package drop that, when considered along with the physical test data collected by the SNL-led team, the SFSWT program can consider the knowledge gap closed.

Klymyshyn, Nicholas A.↗

MODELING A SPENT NUCLEAR FUEL CASK SEISMIC TEST

The US Department of Energy Spent Fuel and Waste Science and Technology (SFWST) program is planning to conduct a series of full-scale shake table tests to simulate hypothetical earthquake conditions and record the response of surrogate spent nuclear fuel (SNF) assemblies in a canister dry storage system mockup. The shake table motions will represent a range of hypothetical earthquake conditions at hypothetical locations in the continental US to generally define the range of mechanical loads that SNF can be expected to experience during extended dry storage periods. The earthquake conditions will represent seismic hazards in the 2,000-to-20,000-year return period range. The test will use instrumented pressurized water reactor fuel assemblies with surrogate mass inside the fuel rods instead of radioactive fuel pellets. Dummy assemblies with similar mass and dimensions to actual fuel assemblies will occupy the rest of the fuel assembly locations in a SNF canister that holds 32 fuel assemblies. The canister will be located inside a mockup vertical concrete cask. Instrumentation will record the motion of the major components of the complex dynamic system, and strain gauges will be used to record the cladding strain at select locations. Preparations for the test require modeling predictions to identify the range of response of the system and to help select specific earthquake cases to be simulated on the shake table from a large set of potential cases. This paper describes the pretest nonlinear finite element modeling efforts that have been completed to date, including cask system level modeling and fuel assembly modeling in LS-DYNA. The cask system level models are critical for anticipating sliding or tipping of an unanchored cask during the test. The fuel assembly model is needed to estimate the range of cladding strain response and fuel assembly structural response to be expected during the test and cladding strain measurements will be key metrics for model validation and the development of modeling best practices after the test is concluded.

Klymyshyn, Nicholas A.↗