Search NASA⌕ Search

SEARCH · Search NASA

Results for “test particle simulations”

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 73 records · Page 4

The structure of dissipative dark matter halos

Dissipative dark matter self-interactions can affect halo evolution and change its structure. We perform a series of controlled N-body simulations to study impacts of the dissipative interactions on halo properties. The interplay between gravitational contraction and collisional dissipation can significantly speed up the onset of gravothermal collapse, resulting in a steep inner density profile. For reasonable choices of model parameters controlling the dissipation, the collapse timescale can be a factor of 10–100 shorter than that predicted in purely elastic self-interacting dark matter. The effect is maximized when energy loss per collision is comparable to characteristic kinetic energy of dark matter particles in the halo. Our simulations provide guidance for testing the dissipative nature of dark matter with astrophysical observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Knowing our Sun: data fusion for optimizing space weather forecasts

Outline: Space weather: model for predicting solar wind & polarity; Particle filtering: optimization with Monte Carlo; Simulation & observation: twin tests and real data; Back to Earth: independent work in inertial confinement fusion (ICF)

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling Turbulence in Galactic Centers

Turbulence is a prevalent phenomenon in the interstellar medium, and in particular, the environment at the centers of galaxies. For example, detailed observations of the Milky Way’s Central Molecular Zone (CMZ) revealed that it has a complex and turbulent structure. Turbulence on galactic scales is often modeled using star formation and feedback. However, these effects do not appear to be sufficient for explaining the high-velocity dispersion observed in the CMZ, indicating that additional gas-stirring processes are likely to be operating. Here we introduce a proof-of-concept method to drive turbulence in gas that orbits under the influence of a galactic potential. Instead of relying on a particular physical mechanism, we have adopted a Fourier forcing module and have applied it using a smoothed particle hydrodynamics code. To test our method, we performed simulations of a simplistic model of the CMZ. Our turbulence injection method is capable of balancing the self-gravity of the gas, which allows us to run the simulations for long timescales and thereby follow the evolution of the CMZ. Our results show that turbulence induces a flocculent spiral pattern in our model, analogous to that found in galactic-scale simulations. Furthermore, we find that our turbulence injection method induces inward migration of gas, a result consistent with previous numerical simulations. We submit that this injection method is a promising new tool to simulate turbulence in galactic centers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Possibilities and Limitations of Kinematically Identifying Stars from Accreted Ultra-faint Dwarf Galaxies

Abstract The Milky Way has accreted many ultra-faint dwarf galaxies (UFDs), and stars from these galaxies can be found throughout our Galaxy today. Studying these stars provides insight into galaxy formation and early chemical enrichment, but identifying them is difficult. Clustering stellar dynamics in 4D phase space ( E , L z , J r , J z ) is one method of identifying accreted structure that is currently being utilized in the search for accreted UFDs. We produce 32 simulated stellar halos using particle tagging with the Caterpillar simulation suite and thoroughly test the abilities of different clustering algorithms to recover tidally disrupted UFD remnants. We perform over 10,000 clustering runs, testing seven clustering algorithms, roughly twenty hyperparameter choices per algorithm, and six different types of data sets each with up to 32 simulated samples. Of the seven algorithms, HDBSCAN most consistently balances UFD recovery rates and cluster realness rates. We find that, even in highly idealized cases, the vast majority of clusters found by clustering algorithms do not correspond to real accreted UFD remnants and we can generally only recover 6% of UFDs remnants at best. These results focus exclusively on groups of stars from UFDs, which have weak dynamic signatures compared to the background of other stars. The recoverable UFD remnants are those that accreted recently, z accretion ≲ 0.5. Based on these results, we make recommendations to help guide the search for dynamically linked clusters of UFD stars in observational data. We find that real clusters generally have higher median energy and J r , providing a way to help identify real versus fake clusters. We also recommend incorporating chemical tagging as a way to improve clustering results.

79 ASTRONOMY AND ASTROPHYSICS↗

Super‐Droplet Method to Simulate Lagrangian Microphysics of Nuclear Fallout in a Homogeneous Cloud

Abstract Nuclear detonations produce hazardous local and global particles or fallout. Predicting fallout size, chemical components, and location is necessary to inform officials and determine immediate guidance for the public. However, existing nuclear detonation fallout models prescribe the particle size distributions based on limited observations. In this work, we apply the super‐droplet method, which is a numerical modeling technique developed for cloud microphysics, to simulate size distributions of particles in a mushroom cloud formed post‐detonation of a nuclear device. We model fallout formation and evolution with homogeneous nucleation and condensation of a single species and a Monte Carlo coagulation algorithm. We verify the numerical methods representing coagulation and condensation processes against analytical test problems. Additionally, we explore several scenarios for the integral system mass and yield in equivalent kilotons (kt) of TNT (trinitrotoluene). The fallout size distribution median diameter d pg follows a scaling law based on the integral system mass m v 0 kg and yield Y kt: nm. We test the effect of cloud turbulence, enhanced nucleation and growth, and vapor volatility with a sensitivity study. The range in median diameter predictions for simulations of historical tests performed over the Pacific encompass the measurements of particles sampled from the cloud caps. Predicted median particle size ranges up to 217, 123, 86, and 35 nm for historical tests with yields of 0.2, 0.7, 2, and 10 Mt, respectively. This work can be expanded in many different directions to build a more predictive model for fallout formation.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Abrasion wear at high temperature in particle receiver type Concentrating Solar Power Systems

The use of solid particles as a heat transfer fluid in receivers for concentrating solar powerplants has garnered much attention in the past decade. Particles are not only stable at temperatures above 700 °C but are also able to absorb solar energy directly, which gives them a distinct advantage over conventional molten salts. The use of particles can, however, cause material degradation from the flow of hot or cold particles through discharge hoppers or along the inner receiver surfaces and other system components (e.g. tubes, valves etc.), depending on operating mode of the receiver. Stainless steel alloy 316, nickel-based alloy Inconel® 740H®, Haynes 230, and a 60% alumina, low-cement refractory castable from Allied Mineral Products were tested for wear from first-of-a-kind high-temperature low-velocity abrasion-erosion experiments. The tests are designed to simulate potential conditions in Generation 3 falling particle Concentrating Solar Power (CSP) systems. Simultaneous experiments were conducted at room temperature as well to draw comparison on the role of temperature on abrasion wear rates. It was noted that the abrasion wear rate in candidate materials was driven by the oxidation resistance of the material. The findings from these experiments will be used to determine life-time of critical components in particle receiver concentrated solar powerplants.

14 SOLAR ENERGY↗

Prediction of the energetic particle redistribution by an improved critical gradient model and analysis of the transport threshold

Based on the theory of critical gradient model (CGM) and following the simulation method proposed by Waltz et al., a combination of TGLFEP and EPtran code is employed to predict the energetic particle (EP) transport induced by Alfvén eigenmodes (AEs). To be consistent with the experiment, recent improvements to the simulation method include consideration of threshold evolution and orbit loss due to finite orbit width. The revised CGM is applied to simulate two DIII-D experimental discharges (#142111 and #153071). It well reproduces the experimental profiles with multiple unstable AEs and large-scale EP transport. Discharge #142111 had previously been simulated using a nonlinear MHD-kinetic code MEGA with a transport mechanism based on stochasticity induced by overlapping AE. By comparing the simulated EP profiles, we find that the AE transport threshold is approximated by both the MEGA nonlinear stability threshold and the proposed CGM threshold (error <5% for single n and <17% for multiple n simulation). Both of them are larger than the linear stability threshold of the most unstable AE mode by a quantity of the order of the flux needed to sustain EP transport by the background turbulence. Additionally, we have also applied the improved CGM to simulate the α particle redistribution for a China Fusion Engineering Test Reactor steady state scenario. Because of the clear separation between the AE unstable region and the loss cone, only a moderate α particle loss of ~9.6% is predicted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Ab initio investigation of a hypersonic double cone experiment

This article presents a direct molecular simulation (DMS) of a reactive Mach 8.2 oxygen flow over a double cone geometry. The free stream conditions and article configuration generate a flow with thermal and chemical nonequilibrium, which are common attributes of hypersonic flight. This scenario was first studied experimentally at Calspan University of Buffalo Research Center’s test facility. DMS is a particle method that uses quantum mechanically derived interaction potentials to simulate molecular collisions within a flow field. Since these interaction potentials are the only modeling inputs used in the simulation, all flow features can solely be attributed to the ab initio potential energy surfaces. Hence, providing a comparison of a hypersonic ground test and numerical data anchored to quantum mechanics. The fundamental nature of DMS is leveraged to investigate molecular level mechanisms prevalent in the flow, and comparisons with lower fidelity simulations are presented to highlight the role of these first principles calculations as benchmark solutions.

Science & Technology - Other Topics↗

Small to Full Height Scale Comparisons of Cesium Ion Exchange Performance with Crystalline Silicotitanate

The U.S. Department of Energy’s (DOE) Hanford Site houses 56 million gallons of high-level radioactive waste generated from plutonium production from 1944 to 1988. The supernatant waste, currently stored in underground tanks, is intended to be vitrified following filtration and 137 Cs removal at the Hanford Waste Treatment and Immobilization Plant (WTP) Low-Activity Waste (LAW) Vitrification Facility. The WTP Pretreatment Facility will not be operational for several years. The Tank Side Cesium Removal (TSCR) system is a technology demonstration that will remove cesium from tank waste supernate to support directly feeding LAW to the vitrification facility. The 137 Cs removal is important to meet the WTP LAW contract specification and ultimately for creating a contact-handled waste form. The waste acceptance criteria (WAC) limit for the WTP LAW Facility is <3.18E-5 Ci 137 Cs per mole Na. The TSCR system is skid mounted and employs two key technologies: dead end filtration for solids removal and ion exchange (IX) for cesium removal. Filtration is necessary to protect the functionality of the IX columns. The IX process utilizes sodium form crystalline silicotitanate (CST) IX media produced by UOP Honeywell, LLC (Des Plaines, IL) under the product name IONSIV R9140-B, 18 x 50 mesh, in a lead-lag-polish column configuration. Each column contains a CST bed height of 92 inches and a volume of approximately 157 gallons of CST IX media. The full-scale TSCR operation will run at a residence time of 1.9 bed volumes (BV)/h, which results in a superficial velocity of 7.3 cm/min and a flowrate of ~5 gal/min. Column testing at small (2.5% of the full bed height), medium (12% of the full bed height), and full-height scales has been previously conducted to evaluate process variables and scale up performance of Cs exchange onto the CST. Cesium load performances from various sieve cuts at the small scale indicated that a 30-mesh sieve cut be tested to determine if it better reflects the 12% and full-height column performances at the small scale. Two process flowrates were tested in the small-column configuration with <30-mesh CST and simulant solutions. These tests were compared to full-height column tests at the same residence times to assess CST particle size effect on column scaling. Table S.1 summarizes the observed column performance determined for the two flowrates juxtaposed to the previous work with 5.6 M Na simulant at higher scales (used as benchmarks). The WAC breakthroughs between the small and full-height tests at 1.3 BV/h differed by ~59 BVs. The WAC breakthroughs for small, medium, and full-height tests at 1.8 BV/h were consistent at 240 BVs. The 50% Cs breakthroughs were nominally equivalent for all column tests. The common transition zones and onsets of Cs breakthrough at the 1.8 BV/h tests indicated that the Cs mass transfers were equivalent and thus the <30-mesh CST Cs load performance at the small scale successfully modeled that of the full-height system. It is recommended that the <30-mesh CST be used in subsequent 10-mL CST bed tests.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

BISON Simulated and Experimental Fission Product Release Comparisons from Reradiated AGR-3/4 Compacts During High Temperature Heating Tests

The fuel performance modeling code BISON was used to predict the release of fission products iodine-131 (131I), xenon-133 (133Xe), and krypton-85 (85Kr) from four re-irradiated AGR-3/4 fuel compacts containing tristructural isotropic (TRISO) coated particles during high-temperature isothermal heating tests. The AGR-3/4 fuel compacts were irradiated in the Advanced Test Reactor (ATR) as part of the third and fourth series of planned experiments to support the Advanced Gas Reactor (AGR) Program. They were subsequently stored and re-irradiated in the Neutron Radiography (NRAD) reactor for approximately five days and then stored for another five to eight days before being subjected to isothermal heating tests in the Fuel Accident Condition Simulation (FACS furnace) for 200 to 300 hours at temperatures between 1000°C and 1600°C to evaluate fission product release at elevated temperatures. New nuclide-specific fission product source term models for the three nuclides of interest were developed using the reactor multiphysics code Griffin and implemented into BISON to support this work. The new source term models were incorporated into coupled compact- and particle-scale BISON simulations, which predict spatially- and temporally-resolved radionuclide generation, radioactive decay, transport, and release throughout the entire irradiation history, including the initial ATR irradiation, NRAD re-irradiations, FACS heating tests, and intermediate periods spent in storage. The experimentally measured fission product release from the heating tests were compared to modeling release predictions calculated by BISON to evaluate how well the code compares to experimental results. Overall, the experimental measured and BISON predicted comparative release results varied but generally agreed to within 5 particle equivalents. Comparative release results identified general observations to take into consideration to help refine future models and reduce uncertainties associated with both the measurement results and predictive results. This includes developing new uranium oxycarbide (UCO) specific kernel diffusivities for the three isotopes examined to more accurately reflect the material properties of the fuel form. Deriving new diffusivities will aid in producing a more informed BISON model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Testing dark matter halo properties using self-similarity

ABSTRACT We use self-similarity in N-body simulations of scale-free models to test for resolution dependence in the mass function and two-point correlation functions of dark matter haloes. We use 10243 particle simulations performed with abacus , and compare results obtained with two halo finders: friends-of-friends (fof ), and rockstar . The fof mass functions show a systematic deviation from self-similarity which is explained by resolution dependence of the fof mass assignment previously reported in the literature. Weak evidence for convergence is observed only starting from haloes of several thousand particles, and mass functions are overestimated by at least as much as $20-25{{\ \rm per\ cent}}$ for haloes of 50 particles. The mass function of the default rockstar halo catalogue (with bound virial spherical overdensity mass), on the other hand, shows good convergence of the order of 50 to 100 particles per halo, with no detectable evidence at the few percent level of any systematic dependence for larger particle number. Tests show that the mass unbinding procedure in rockstar is the key factor in obtaining this much improved resolution. Applying the same analysis to the halo–halo two point correlation function, we find again strong evidence for convergence only for rockstar haloes, at separations sufficiently large so that haloes do not overlap. At these separations, we can exclude dependence on resolution at the $5-10{{\ \rm per\ cent}}$ level once haloes have of the order of 50 to 100 particles. At smaller separations results are not converged even at significantly larger particle number, and bigger simulations would be required to establish the resolution required for convergence.

79 ASTRONOMY AND ASTROPHYSICS↗

Hybrid Solar System (Final Scientific/Technical Report)

GTI Energy (GTI) teamed with the University of California at Merced (UCM) to scaleup the hybrid solar system (HSS) technology for demonstrating its performance at the US Gypsum (USG) plant in Plaster City, California. The technology integrates two-stage concentrating solar collector with matching particle thermal transport and storage (TSS) system to deliver cost-effective, and on-demand distributed high temperature industrial process heat up to 600°C with solar thermal, in this case to a gypsum kettle, to reduce its fuel use and carbon footprint. Current solar technologies, which reach these temperatures, are not distributable (towers) or cost-effective (dish). The research team developed a conceptual system design for host site retrofit, including preliminary heat balance, process flow diagram, particle to process heat exchanger and equipment placements at the site. Subsequently, parallel efforts were carried out at UCM to design, build and test a 12 m long commercial scale prototype concentrating thermal-only collector system and at GTI to design, build and test a matching 650°C capable particle TTS system. The nominal 50 kWth collector consists of a parabolic trough and three 4 m long two-stage receivers in series. Prior to on-sun testing, a 4 m long receiver was fabricated and successfully tested at 650 °C in a laboratory setting for 100 hrs of continuous operation showing less than 15% radiation loss. A 7 m wide x 17 m long parabolic trough was then installed at UCM for on-sun testing of the 12 m long receiver, and concurrently several 4 m long receivers were built. The optics of the parabolic trough were calibrated, and on-sun test were carried out on 12 m long receivers. During tests, the intense solar radiation (53x) caused the absorber tubes in the receivers to bend, reducing the overall optical efficiency. To address the bending issue, a self-consistent algorithm that includes ray tracing, thermal and deformation models was developed to perform thermal stress analysis on absorbers for parabolic solar collectors. Results obtained with this algorithm showed a dramatic rise in deformation as absorber tube length increases. A combined efficiency parameter that includes the occluded area for the mounts was developed to obtain an optimized tube length obtained. Based on the results, a length of 2.7 m for the absorber + 0.2 m for the coupler was chosen to minimize any bending and optimize optical efficiency while maintaining ease of mounting. The associated particle TTS system was designed, built and successfully tested at GTI. It includes storage, receiving and lock hoppers and piping that simulates the transfer of captured solar energy to an actual industrial furnace. Tests over 77 charge-discharge cycles demonstrated <2% particle degradation, with no problematic particle accumulations and no flow interruptions. The piping pressure drop was about 5 psi. The team also worked with Stanley Consultants (Stanley) to prepare conceptual and preliminary engineering packages to facilitate follow-on development and commercialization efforts. These include process and instrumentation diagram’s (P&ID’s), general arrangements, electrical one-line, project definitions document, equipment data sheets, schedule, and construction cost estimate for 2 MWth system. Updated HSS technology commercialization and customer engagement plans and detailed costs and evaluated market trade-offs and manufacturing.

03 NATURAL GAS↗

Developing Low Noise, Rad-hard Detectors for Parity Violating Experiments

As parity-violating electron scattering (PVeS) experiments continue to push the frontiers of precision electro-weak asymmetry measurements, the demands on the experimental techniques and apparatus are also pushed. In particular, the need for higher resolution and radiation-hardness of the main integrating detectors has increased over the successive generations of PVeS experiments. The demand for precise measurement of GHz event rates has pushed the field toward the use of thin quartz Cerenkov light based detectors. We have been designing and testing such thin quartz detectors specifically for the upcoming experiments at Jefferson Lab, PREX-II and CREX, but this work also greatly influences the detector development for several future experiments including MOLLER at Jefferson Lab. These are all PVeS experiments that will use the new thin quartz Cerenkov detector design concept for their main asymmetry measurements (as well as for beam monitoring). The new design concept gives not only significant performance improvements compared to its predecessor from PREX-I, but also we have now thoroughly characterized its operational design and performance using a combination of test-beam data and detailed particle and optical Monte Carlo Geant4 simulations. These activities have culminated in the development and implementation of a "bench-marked" Monte Carlo package, QSIM, which constitutes a powerful design tool for present and future PVeS quartz detectors. The benchmarked simulation can replicate real photoelectron distributions (RMS and Mean) from testbeam data with ~5% precision – limited mainly by the systematic uncertainty of the PMT gain measurements.

Villarreal, Carlos Bula Villarreal↗

Electrochemical Corrosion of SIMFUEL: Effects of dissolved H2 and noble metal particles

This report describes the results from electrochemical corrosion tests that were conducted using simulated spent fuel (SIMFUEL) materials comprised of UO2 and surrogate fission products. Two SIMFUEL compositions were tested to quantify the effect of noble metal inclusions and dissolved H2 on the UO2 dissolution rate. One material consisted of UO2 with added lanthanide oxides (UO2 N) and the other consisted of UO2 with added lanthanide oxides and noble metals (UO2-H) at concentrations to simulate high burnup fuel. The electrochemical corrosion tests on the SIMFUEL materials were conducted in aqueous electrolyte solutions (pH 10) that were either saturated with air or purged with a H2/Ar gas mixture to maintain a constant dissolved H2 concentration. The results from electrochemical corrosion tests on SIMFUEL can be applied to qualitatively understand and also quantify the separate effects of water chemistry and fuel composition on the degradation behavior of the UO2 matrix. Open circuit potential (OCP) measurements on a SIMFUEL material of known composition (i.e., known fraction of NMPs at the fuel surface) immersed in a known solution chemistry (i.e., pH, Eh, [O2], [H2]) at known temperature provide qualitative insight into the degradation behavior of the UO2 matrix. If the OCP is above the threshold potential at which the oxidative dissolution of U(IV) to U(VI) occurs, the SNF is expected to degrade by oxidative processes under the experimental conditions. The net currents that are measured during potentiostatic tests, during which the surface potential of the SIMFUEL material is fixed by a potentiostat, can be used to quantify the UO2 degradation rate and optimize the rate constant values used in the Fuel Matrix Degradation Model (FMDM) for half reactions that occur on the SNF surface. Specifically, electrochemical measurements enable the estimation of the total anodic current at ECORR–which is the surface potential at which the total anodic and total cathodic currents are equal–so that the rate constant values for key reactions can be calculated. The open circuit potential measurements, potentiostatic tests, surface property measurements (scanning electron microscopy images and electrochemical impedence spectroscopy plots), and solution elemental composition analyses are being performed to update and optimize the FMDM. A case study is presented herein to show how electrochemical corrosion test results and accompanying characterization results for the UO2 N material in air-saturated solution can be used to validate the SNF surface reaction module of the FMDM. Future electrochemical tests will be conducted to provide quantitative information on the effects of NM content (burnup), H2 concentration, water chemistry, temperature, and galvanic couples with cladding and EBS alloys on UO2 degradation that can be used to improve the accuracy and functionality of the FMDM.

Thomas, Sara↗

Ordered Particle Packing in Dense TRISO/SiC Fuel Elements and Preliminary Assessment of Neutronic and Thermomechanical Characteristics

Detailed analysis of the particle distribution in Transformational Challenge Reactor fuel elements indicates that particle packing is not random; instead, it follows a relatively ordered structure near fuel element surfaces. Discrete particle neutronic simulations indicate that the core reactivity is not impacted when assuming homogenization of particles with the silicon carbide matrix. However, the neutronic power distribution resulting from the ordered packing structure indicates that the highest-power particles reside at the top and bottom of the fuel elements and nearest the YH 1.85 moderator rods. The power distribution results were applied to thermo-mechanical simulations using mesh-based power distributions. Previous results indicated high stress at the bottom of the fuel element, where packing is most ordered. Additionally, to reduce this stress concentration, additively manufactured protrusions were added to the bottom of a test fuel element to disrupt dense particle packing. These protrusions reduced the overall power peaking, but the thermomechanical simulations did not indicate a significant change in the fuel element’s maximum stress or failure probability.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Computational model development and validation of fuel dispersal phenomena

As the US nuclear industry is renewing efforts to extend the rod average burnup limits in existing pressurized water reactors, it becomes imperative to investigate the safety implications associated with high burnup fuel fragmentation and the potential dispersal of fragmented fuel into the reactor coolant system. Additionally, as there is a growing interest in transitioning to new nuclear fuel designs to accommodate higher fuel burnup levels, it becomes essential to conduct studies on fuel behavior during design basis accident scenarios, particularly during a postulated loss-of-coolant accident, in order to ensure the safe operation of light-water reactors. The present study develops a computational model to simulate the complex three-phase flow of the fuel dispersal phenomena that occur following a breach in the cladding, when the stream of solid fuel particles and carrier fission gas stream may interact with the surrounding liquid medium. Furthermore, the proposed framework considers the solid phase as discrete Lagrangian particles while the interpenetrating gas–liquid continuum is modeled using a Eulerian framework. The developed simulation framework is validated with experimental results performed in a separate-effect test facility by comparing predicted particle settlement locations with the experimental observations reported in the literature. Key highlights of this study include the interpretation of the high-pressure boundary conditions and the three-phase flow coupling strategy. Also discussed are the gas–liquid hydrodynamics and the behavior of particles as they are transported through the gas and liquid phases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Assessment of Buffer-IPyC Thermomechanical Debonding Behavior Using New Experiment Strength Data in BISON

TRIstructural ISOtropic (TRISO) single sub-millimeter-diameter nuclear fuel particle consists of a spherical fuel kernel surrounded by four coating layers. During reactor operation, the TRISO coating layers are subjected to irradiation-induced dimensional changes and the associated thermomechanical behavior of each layer. One of the observed behaviors is gap formation between the buffer and IPyC layer due to the porous buffer’s irradiation-induced shrinkage exceeding that of the IPyC layer. To predict TRISO thermomechanical performance, multi-physics models have been built that are continually updated and refined. The BISON code is a finite element-based nuclear fuel performance code that may be used for TRISO particle simulations. One of the recent additions to the BISON code is the ability to model the process of layer debonding. This paper will focus on the simulation results of the improved BISON debonding model that will utilize updated strengths measured from irradiated TRISO fuel particles. The new experimental strength data from micromechanical tests of irradiated TRISO fuel samples were exercised in the BISON simulations and compared to baseline strength data to assess their applicability in the models. This also includes updated buffer-IPyC bond strengths to simulate layer delamination. Based on current experimental observations it is noted that the buffer-IPyC separation occurs not exactly at the junction of these two layers, but more on the side of the buffer layer. This improved modeling approach using experimental strength data to characterize buffer-IPyC debonding and its potential subsequent cracking will be presented in the paper along with comparisons to available experimental observations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Validation of Computational Fluid Dynamics Simulations for Biological Performance Assessment in Hydropower units (Final Report)

The biological performance assessment (BioPA) toolset developed by Pacific Northwest National Laboratory (PNNL) estimates the relative biological performance of fish passage at a hydroelectric power turbine unit. The tool is based on the use of computational fluid dynamics (CFD) and fish biological response relationships. The recent release, BioPA-v3, is based on directly computed trajectory and collision of material Lagrangian particles using CFD simulation codes rather than the prior version that relies on Tecplot to compute streamtrace trajectories. Before modifying the toolset, a series of validation tests were performed at the various steps of modification in the toolset. Validation is a critical step of any numerical investigation that reflects the accuracy and reliability of the predicted results. It raises the confidence level of the user to use the modified version of the BioPA toolset. Several test cases were simulated and compared, where available, to observed data. The trajectory and collision of the small spherical and cylindrical particles in a water flume were compared to in-house experiments. The CFD predicted collision rate and flow field compared well with experimental observation for vane array and large cylinder as target bodies. Next, the CFD-predicted flow field and hydraulic performance of a laboratory-scale model of a Francis turbine was also successfully validated. Note that the trajectory of the particles is significantly affected by the flow field in such extreme conditions. In addition to the particle trajectories and flow field, the collision detection method employed in the CFD simulations was also successfully validated. The CFD predicted impact velocity, collision time, velocity, and trajectory of a sphere excellently matched with analytical value for a bouncing ball in the elastic collision. A similar approach was also tested and successfully validated for a collision of sphere with a 45° inclined plane. After successfully validating different cases, the BioPA toolset was modified to use direct output of the CFD prediction and the new version can be used in evaluating biological performance at hydroelectric turbines.

13 HYDRO ENERGY↗