Nitrogen Evaluation [Slides]
This presentation includes a R-Matrix Analysis of 15N system with SAMMY, Discussion of Boundary Conditions in R-Matrix Analyses and a discussion of Analyses with Multiple Incident Channels [Slides]
SEARCH · Search NASA
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.
This presentation includes a R-Matrix Analysis of 15N system with SAMMY, Discussion of Boundary Conditions in R-Matrix Analyses and a discussion of Analyses with Multiple Incident Channels [Slides]
Software Health Management (SWHM) must make sure that the software always remains in safe and healthy regions of the state space. Boundaries between healthy and unhealthy regions are important for the detection of violations and health management.In this position paper, we present a framework, which employs techniques from Bayesian statistical modeling and active learning to efficiently characterize health boundaries in high-dimensional spaces. We will discuss, how this framework supports SWHM during design time and during operation of learning/adapting software systems
During the late 1950s, the Atomic Energy Commission (AEC) (now the U.S. Department of Energy [DOE]) conducted a series of safety experiments to determine if a nuclear device subjected to a large conventional explosives detonation would result in a nuclear yield. The AEC obtained temporary use of a large portion of western Emigrant Valley from the U.S. Department of Defense (DOD) for one of these experiments: Project 57. Following the Project 57 safety experiment, the AEC fenced the contaminated area and returned control of the surrounding land to the DOD. Emigrant Valley is part of the Nevada Test and Training Range (NTTR). For safety and security reasons, access to the NTTR is controlled using both physical (i.e., fences) and administrative (e.g., signs and postings) controls. Therefore, the public cannot access the Project 57 site and there are no known human receptors that routinely access the site. Project 57 was detonated on April 24, 1957, in Emigrant Valley approximately 13 mi (21 km) northeast of the north end of Yucca Flat. This test was undertaken to develop: (1) a means of estimating immediate distribution and long-term redistribution of plutonium dispersed during a nonnuclear detonation; (2) biomedical evaluation techniques for use in likely plutonium-laden environments; (3) methods of decontaminating ground areas, pavements, and building materials; and (4) alpha survey instruments and field monitoring procedures to promptly estimate contaminant deposition. Although the test did not result in the fission of nuclear materials, it did disseminate plutonium across the ground surface. Various radiological surveys have been performed in the area since Project 57 was conducted. The AEC delineated the area containing radioactive material based on radioactivity surveys performed shortly after the test was conducted. A designated Contamination Area (CA), marked by a fence, was established later. The distribution of americium-241 (Am-241) in the area was determined again in a 1997 flyover, which showed Am-241 ranging from as much as 70,000 counts per second (cps) at ground zero to background (<70 cps). This survey also documented Am-241 radioactivity on the ground surface beyond the east side CA fence at levels of up to 150 cps. In 2007, the DOE expanded the CA by posting “Contamination Area” signs 200 ft to 400 ft (60 m to 120 m) beyond the original fence, which formed a new CA boundary.
The overall objective of the project was to develop a computationally feasible and user-friendly mechanized process to integrate traditional probabilistic risk assessment (PRA) and dynamic PRA (DPRA) results. Starting with the systematic identification of items in an existing PRA that need dynamic augmentation, the project used a generic 4-loop pressurized reactor (PWR) and 3-loop PWR as example plants. Station blackout (SBO) and large break loss of coolant accident SBLOCA) were selected as the example initiating events. Using the traditional event-tree (ET)/fault-tree (FT) methodology augmented by dynamic evet tree approach, the potential consequences of the initiating events were simulated with RELAP-3D and MELCOR/RASCAL codes to cover Level 1 through Level 3 of PRA. RAVEN and ADAPT software were used to generate Level 1 simulations with RELAP-3D and Level 2/3 simulations with MELCOR (Level 2)/RASCAL (Level 3), respectively. Example branching conditions (BCs) for SBO included AC power recovery time, valve repair failure time, reactor coolant pump leak time/break size and emergency power supply duration to a total of 9. Example BCs for LOCA included off-site power recovery time, diesel generator power recovery time, auxiliary feed water system operation time, safety relief valve failure to open upon demand, reactor coolant pump seal break time and size to a total of 21. Each RELAP-3D simulation (9,587 scenarios) was labelled OK or Core Damage based on the maximum allowed peak clad temperature (2,100oF). Each MELCOR simulation (4610 scenarios) was labeled as Bin over 10rem or Bin 0-10rem based on the dose at the site boundary. The scenarios were clustered based on the criteria above using the mean shift methodology. Classical PRA (CPRA) and DPRA results were compared to identify the ET sequences that need DPRA augmentation. Several approaches were proposed for the incorporation of these sequences into CPRA using clustering with the mean shift methodology, restructuring the CPRA ETs by adding new BCs/sequences, and using the concept of a limit surface. Procedures for decision making regarding the possible consequences of an initiating event (e.g. core damage or not, site evacuation or not) were developed using a convolutional neural network (CNN), a recurrent neural network (RNN) and a transformer neural network (TNN). The project has led to two PhD degrees, three archival journal papers and five refereed conference proceedings.
The Seventh Summer Program of the Center for Turbulence Research took place in the four-week period, July 5 to July 31, 1998. This was the largest CTR Summer Program to date, involving thirty-six participants from the U. S. and nine other countries. Thirty-one Stanford and NASA-Ames staff members facilitated and contributed to most of the Summer projects. A new feature, and perhaps a preview of the future programs, was that many of the projects were executed on non-NASA computers. These included supercomputers located in Europe as well as those operated by the Departments of Defense and Energy in the United States. In addition, several simulation programs developed by the visiting participants at their home institutions were used. Another new feature was the prevalence of lap-top personal computers which were used by several participants to carry out some of the work that in the past were performed on desk-top workstations. We expect these trends to continue as computing power is enhanced and as more researchers (many of whom CTR alumni) use numerical simulations to study turbulent flows. CTR's main role continues to be in providing a forum for the study of turbulence for engineering analysis and in facilitating intellectual exchange among the leading researchers in the field. Once again the combustion group was the largest. Turbulent combustion has enjoyed remarkable progress in using simulations to address increasingly complex and practically more relevant questions. The combustion group's studies included such challenging topics as fuel evaporation, soot chemistry, and thermonuclear reactions. The latter study was one of three projects related to the Department of Energy's ASCI Program (www.llnl.gov/asci); the other two (rocket propulsion and fire safety) were carried out in the turbulence modeling group. The flow control and acoustics group demonstrated a successful application of the so-called evolution algorithms which actually led to a previously unknown forcing strategy for jets yielding increased spreading rate. A very efficient algorithm for flow in complex geometries with moving boundaries based on the immersed boundary forcing technique was tested with very encouraging results. Also a new strategy for the destruction of aircraft trailing vortices was introduced and tested. The Reynolds Averaged Modeling (RANS) group demonstrated that the elliptic relaxation concept for RANS calculations is also applicable to transonic flows with shocks; however, prediction of laminar/turbulent transition remains an important pacing item. A large fraction of the LES effort was devoted to the development and testing of a new algorithmic procedure (as opposed to phenomenological model) for subgrid scale modeling based on regularized de-filtering of the flow variables. This appears to be a very promising approach, and a significant effort is currently underway to assess its robustness in high Reynolds number flows and in conjunction with numerical methods for complex flows. As part of the Summer Program two review tutorials were given on Turbulent structures in hydrocarbon pool fires (Sheldon Tieszen), and Turbulent combustion modeling: from RANS to LES via DNS (Luc Vervisch); and two seminars entitled Assessment of turbulence models for engineering applications (Paul Durbin) and Subgrid-scale modeling for non-premixed, turbulent reacting flows (James Riley) were presented. A number of colleagues from universities, government agencies, and industry attended the final presentations of the participants on July 31 and participated in the discussions. There are twenty-six papers in this volume grouped in five areas. Each group is preceded with an overview by its coordinator.
In this study, the capabilities of a physics-informed convolutional neural network (CNN) for reconstructing the temperature field from a limited set of measurements taken at the boundaries of internal flows are demonstrated. Such an approach enables the development of less invasive monitoring methods for real-time plant diagnostics. As a test case, a Molten Salt Fast Reactor (MSFR) design was selected. This circulating fuel reactor has received interest from both scientific and industrial communities due to its intrinsic safety and sustainability. Molten salt flows in such reactors, however, can present highly localized temperature peaks that can induce significant thermal stresses onto the vessel walls. At these local maxima, the salt temperature may exceed a thousand kelvins, which makes a direct measurement challenging or even unfeasible. The proposed CNN algorithm allows one to detect indirectly such discontinuities through an accurate, albeit indirect, temperature measurement method during reactor operation. The datasets employed to train and test the machine learning models in the present work were generated with Nek5000, a computational fluid dynamics (CFD) code developed at Argonne National Laboratory. The CNN algorithm is trained with CFD results that span a set of MSFR operational power and flow ranges. Here, to demonstrate the efficacy of the algorithm, predictions are made for test cases contained within the training range but for which the CFD data were not used when training. Results demonstrate that the proposed technique properly characterizes temperature peaks and distributions within the domain for a broad range of scenarios.
Uncertainties in the manufacturing, fabrication process, material behavior, loads, and boundary conditions results in the variation of the stresses and strains induced in the flexures and its fatigue life. Past experience and the test data at material coupon levels revealed a significant amount of scatter of the fatigue life. Owing to these facts, the design of the flexure, using conventional approaches based on safety factor or traditional reliability based on similar equipment considerations does not provide a direct measure of reliability. Additionally, it may not be feasible to run actual long term fatigue tests due to cost and time constraints. Therefore it is difficult to ascertain material fatigue strength limit. The objective of the paper is to present a methodology and quantified results of numerical simulation for the reliability of flexures used in the Stirling convertor for their structural performance. The proposed approach is based on application of finite element analysis method in combination with the random fatigue limit model, which includes uncertainties in material fatigue life. Additionally, sensitivity of fatigue life reliability to the design variables is quantified and its use to develop guidelines to improve design, manufacturing, quality control and inspection design process is described.
This report describes a confirmatory calculation for a proposed consolidated interim storage facility (CISF) to support the US Nuclear Regulatory Commission’s review of the license application submitted by the Interim Storage Partners. The scope of this report is limited to an assessment of the annual dose from direct radiation associated with Phase 1 of the Waste Control Specialists (WCS) CISF and the determination of the minimum controlled area boundary based on the 10 CFR 72.104 limit of 25 mrem annual dose to the whole body. The confirmatory dose rate calculations used the source term and shielding calculation capabilities of the SCALE 6.2.3 computer code system. The calculation method employed in the confirmatory calculations uses a detailed Monte Carlo radiation transport simulation from source to dose rate. This method differs from the two-step method used in the safety analysis report (SAR), which requires determination of the photon and neutron energy and angular distributions on the cask external surface and then the use of this surface source in a new radiation transport calculation to determine dose rate as a function of distance from the storage facility. A simulation was made using the complete site geometry (all casks present) but with only one cask containing source, and a total of 467 independent calculations were performed to obtain the dose rate maps produced by each storage cask. Basic input data and assumptions used in the SAR (e.g., cask design parameters and design basis assembly characteristics) were used in this confirmatory calculation for consistency with the SAR. Local atmospheric conditions for skyshine calculations were simulated using US standard atmosphere data, which describe average air temperature, pressure, and density as a function of altitude. The determined minimum required distances from the WCS Phase 1 CISF to a contour defined by the 25 mrem annual dose limit were approximately 623 m (0.387 miles) in the N-NE direction, 594 m (0.369 miles) in the SSW direction, and 533 m (0.331 miles) in the E-SE direction and the W-NW direction. The annual dose results of the confirmatory calculation are in relatively good agreement with the annual dose values presented in the WCS SAR. The annual dose produced by the Phase 1 CISF at the nearest site boundary, 0.75 miles away from the CISF, was estimated by this study to be approximately 6.02E-02 ± 2.0E-04 mrem, which is significantly smaller than the 25 mrem annual dose limit provided in 10 CFR 72.104. The SAR annual dose rate value for the nearest site boundary was 7.52E-2 mrem, which included 7.73E-03 mrem due to postulated leakage of the FO-, FC-, and FF-canisters. A comparison with the dose rate values reported in the SAR for 16 detector locations, D1 through D16, on the protected area boundary is also included in this report. This comparison shows that the SAR determined higher annual dose values for detectors D1 through D14 by 12% to 45% and lower annual dose values for detectors D15 and D16 by 18% and 8%, respectively, as compared with this ORNL confirmatory calculation. Sensitivity studies have been performed to support a comprehensive review of the application. In these studies, the effects on dose rate of the modeled air volume, density, and humidity were analyzed. The sensitivity calculations showed that the air volume used in the Monte Carlo simulations was adequate and that a larger air volume would not further increase the dose rate due to skyshine contributions. The sensitivity analysis showed that simulation of local atmospheric conditions is important for the accuracy of dose rate estimates at large distances from the facility. Air humidity was demonstrated to have negligible effects on dose rate values for the simulated local atmospheric conditions.
Abstract. For the underground disposal of high-level nuclear waste in rock salt formations, the safety concept includes the backfilling of open cavities with crushed salt. For the prognosis of the sealing function of the backfill for the safe containment of the radioactive waste, it is crucial to have a comprehensive process understanding of the crushed-salt compaction behavior. The crushed-salt compaction process is influenced by internal properties (e.g., grain size, mineralogy, and moisture content) and boundary conditions (e.g., temperature, stress state, and compaction rate) and, therefore, involves several coupled thermal–hydro–mechanical (THM) processes (Hansen et al., 2014; Kröhn et al., 2017). With the paradigm shift from the limited release of radionuclides to safe containment due to the German Repository Site Selection Act passed in 2017, the importance of crushed salt as geotechnical barrier has increased, with a focus on the evolution of its hydraulic properties. Based on the knowledge gaps in the current process understanding, the “Compaction of crushed salt for safe containment” (KOMPASS) projects were initiated to improve the scientific basis behind using crushed salt for the long-term isolation of high-level nuclear waste within rock salt repositories. The efforts to improve the prediction of crushed-salt compaction begun during the first phase of the KOMPASS projects (Czaikowski et al., 2020) and were followed up in a second phase ending in June 2023. The primary achievements of the projects are as follows (Czaikowski et al., 2020; Friedenberg et al., 2022): specification of the KOMPASS reference material, an easily available and reproducible synthetic crushed-salt material, for generic investigations; development of pre-compaction methods and successful production of samples in the short term and under in situ loading conditions; formulation of an extended laboratory program addressing the isolated investigation of known relevant factors influencing the compaction behavior of crushed salt (Düsterloh et al., 2022); execution of long-term compaction tests addressing isotropic and deviatoric load changes, temperature, and compaction state; construction of a backfill body using the KOMPASS reference material in the Sondershausen mine through collaboration with the SAVER (Entwicklung eines salzgrusbasierten Versatzkonzepts unter der Option Rückholbarkeit) project (Schaarschmidt and Friedenberg, 2022); advancement of the tools for microstructure investigation methods (Svensson and Laurich, 2022); generation (first stages) of a microphysical process list combining literature research with our own findings; benchmarking of long-term compaction test for model development and optimization of various existing models as well as the development of new models; application of a virtual demonstrator (2D model representing a backfilled drift in rock salt) for the visualization of developments and the quantification of the models (Rabbel, 2022). In summary, the KOMPASS projects contributed to the reduction of uncertainties and the strengthening of the safety case for using crushed salt within rock salt repositories.
This paper extends the analysis first presented in Jardin et al. [Phys. Rev. Lett. 128, 245001 (2022)] to more thoroughly examine the stability of spherical torus equilibrium to ideal magnetohydrodynamic (MHD) infernal modes and their nonlinear consequences. Here we demonstrate that in a 3D resistive magnetohydrodynamic (MHD) simulation of a NSTX discharge, anomalous transport can occur due to these instabilities. Here we generate a family of equilibrium of differing β and use this to show that these instabilities could explain the experimentally observed flattening of the electron temperature profile at modest β. The modes studied in this paper are found to occur with poloidal mode number m and toroidal mode number n when the ratio m/ n is in the range of 1.2–1.5, when the central safety factor is in this range or slightly lower, and when the central region has very low magnetic shear. Our analysis gives some insight as to why the unstable linear growth rates are oscillatory functions of the toroidal mode number n. We present a simulation of an initially stable configuration that passes through a stability boundary at a critical β as it is heated. We also show that a particular NSTX discharge is unstable to these modes over a timescale of several hundred ms. We conclude that these modes must be taken into account when performing predictive modeling. An appendix shows that similar modes can be found in R / a = 4 tokamaks for certain q-profiles and β values.
The primary objective of the NASA Advanced Air Transportation Technologies (AATT) program is to fully explore the possibilities of the "Free Flight" concept. AVIATION products will enable substantial increases in the effectiveness of national and global air transportation systems. These increases will be achieved by developing and testing automation aids that can assist in the decision-making process among pilots, air traffic controllers, and dispatchers. The AATT program is responsible for defining, exploring, and developing advanced air traffic system concepts to a level suitable for preproduction prototype assessment by the FAA which, if successful, will result in full-scale deployment. These decision support tools will allow all airspace users to choose the best flight path for their own purpose within the constraints of safety and the needs of other users. To do this, several goals must be met: allow users to minimize operating costs by making trade-offs between time and routing; improve the effectiveness of high-density operations on the ground and in the air, enable safe operation in a smooth and efficient manner across boundaries of free-flight and capacity-constrained flight regions; provide system improvements that are easily deployable anywhere in the world; and improve the ability to simulate advanced capabilities in the airspace system.
The Fourth Convection and Moisture Experiment (CAMEX 4) was a scientific field experiment based in Florida in summer 2001 focused on the study of hurricanes off the east coast of the United States. Sponsored by the National Aeronautics and Space Administration's Office of Earth Science, and conducted in collaboration with the National Oceanic and Atmospheric Administration's annual hurricane research program, CAMEX 4 utilized aircraft, ground-based and satellite instrumentation to obtain unprecedented, three dimensional characterizations of these important storms. The Aerosonde UAV was selected by NASA to participate in CAMEX 4 because it provided a unique capability to obtain measurements in the atmospheric boundary layer in and around the storms, unattainable by other platforms or measurement capabilities. This talk focuses on the NASA review process that was followed to coordinate the UAV activity with the conventional aircraft operations, as well as with the other participating agencies and the FAA. We will discuss how Aerosonde addressed the issues of safety, coordination and communication and summarize the lessons learned.
Future human exploration of Mars will require detailed knowledge of the surface and upper several meters of the subsurface in potential landing sites. Likewise, many of the Planetary Science Decadal Survey science goals, such as understanding the history of Mars climate change, determining how the surface was altered through processes like volcanism and fluvial activity, and locating regions that may have been hospitable to life in the past, would be significantly advanced through mapping of the upper meters of the surface. Synthetic aperture radar (SAR) is the only remote sensing technique capable of penetrating through meters of material and imaging buried surfaces at high (meters to tens-of-meters) spatial resolution. SAR is capable of mapping the boundaries of buried units and radar polarimetry can provide quantitative information about the roughness of surface and subsurface units, depth of burial of stratigraphic units, and density of materials. Orbital SAR systems can obtain broad coverage at a spatial scale relevant to human and robotic surface operations. A polarimetric SAR system would greatly increase the safety and utility of future landed systems including sample caching.
The James Webb Space Telescope (JWST), set to launch in 2018, is NASAs next-generation flagship telescope. The Optical Telescope Element (OTE) and Integrated Science Instrument Module (ISIM) contain all of the optical surfaces and instruments to capture and analyze the telescopes infrared targets. The integrated OTE and ISIM are denoted as OTIS, and will be tested as a single unit in a critical thermal-vacuum test in mid-2017 at NASA Johnson Space Centers Chamber A facility. The payload will be evaluated for workmanship and functionality in a 20K simulated flight environment during this thermal-vacuum test. However, the sheer thermal mass of the OTIS payload as well as the restrictive gradient, rate, and contamination-related constraints placed on test components precludes rapid cooldown or warmup to its steady-state cryo-balance condition. Hardware safety considerations precludes injection of helium gas for free molecular heat transfer. Initial thermal analysis predicted that transient radiative cooldown from ambient temperatures, while meeting all limits and constraints, would take 33.3 days; warmup similarly would take 28.4 days. This paper discusses methods used to reduce transition times from the original predictions through modulation of boundary temperatures and environmental conditions. By optimizing helium shroud transition rates and heater usage, as well as rigorously re-examining previously imposed constraints, savings of up to three days on cooldown and up to a week on warmup can be achieved. The efficiencies gained through these methods allow the JWST thermal test team to create faster cooldown and warmup profiles, thus reducing the overall test duration and cost, while keeping all of the required test operations.
The James Webb Space Telescope (JWST), set to launch in 2018, is NASAs next-generation flagship telescope. The Optical Telescope Element (OTE) and Integrated Science Instrument Module (ISIM) contain all of the optical surfaces and instruments to capture and analyze the telescopes infrared targets. The integrated OTE and ISIM are denoted as OTIS, and will be tested as a single unit in a critical thermal-vacuum test in mid-2017 at NASA Johnson Space Centers Chamber A facility. The payload will be evaluated for workmanship and functionality in a 20K simulated flight environment during this thermal-vacuum test. However, the sheer thermal mass of the OTIS payload as well as the restrictive gradient, rate, and contamination-related constraints placed on test components precludes rapid cooldown or warmup to its steady-state cryo-balance condition. Hardware safety considerations precludes injection of helium gas for free molecular heat transfer. Initial thermal analysis predicted that transient radiative cooldown from ambient temperatures, while meeting all limits and constraints, would take 33.3 days; warmup similarly would take 28.4 days. This paper discusses methods used to reduce transition times from the original predictions through modulation of boundary temperatures and environmental conditions. By optimizing helium shroud transition rates and heater usage, as well as rigorously re-examining previously imposed constraints, savings of up to three days on cooldown and up to a week on warmup can be achieved. The efficiencies gained through these methods allow the JWST thermal test team to create faster cooldown and warmup profiles, thus reducing the overall test duration and cost, while keeping all of the required test operations.
When developing this course several years back I asked myself, "What do I think the TFAWS community needs to know about lithium-ion batteries and why should they care?" I think the answers to these questions can be boiled down to the following statements: 1. Lithium-ion (Li-ion) battery electrical performance and efficiency are heavily driven by thermal conditions. 2. Li-ion battery assemblies can experience single cell thermal runaway events which can lead to cell-to-cell propagation; these are thermally driven failure events that can be controlled with effective thermal management systems. 3. Knowledge of BOTH Li-ion battery fundamentals AND traditional thermal design principles are required to develop safe battery assemblies that are both gravimetrically and volumetrically optimized.This year the Short Course on Lithium-ion Batteries will be presented in a 2-part series. The first portion, "Fundamentals, Battery Safety, and Thermal Runaway," will focus on educating participants on the fundamental aspects of Li-ion batteries and on battery safety related topics (i.e. thermal runaway, cell-to-cell propagation, safe handling practices, et...). The second portion, "Practical Thermal Simulation Techniques," will provide participants with a real time demonstration of how to use Thermal Desktop to model of Li-ion battery assemblies. This part of the lesson will cover geometry simplification, mesh development, boundary conditions, and heat loads for both nominal operations (charge and discharge) and for thermal runaway events.
In a previous publication (K. P. Travis, D. Burley and F. G. F. Gibb, WM2017 Conference, Phoenix Arizona paper no. 17480) we introduced a numerical model to estimate the time it would take for a pressure perturbation to subside following its creation from the construction of a deep geological borehole in water-saturated rock. The model was based on the notion of a point source of momentum - a solution of the time-dependent pressure-diffusion equation. Using the model, we estimated that it would take on the order of 10 k years for physical equilibrium to become re-established following the sinking of a 5 km borehole in granite. The significance of such a model is that it places an upper bound on the lifetime required of a borehole sealing system - a necessary input to a borehole post-closure safety case assessment. No engineered sealing system has ever been devised which is capable of retaining its sealing properties for 10 half-lives of long-lived radioisotopes in spent fuel or high-level waste. One of the key advantages of disposing of nuclear waste via Deep Borehole Disposal (DBD) is the natural sealing provided by density stratified groundwater. Upon sinking of a borehole, and subsequent filling with fresh water, brine or drilling mud, this natural barrier may be temporarily damaged. Over time, fresh brine from the far-field will flow towards or away from the hole (driven by a pressure gradient) and will eventually re-establish the original salinity gradient. It follows that the engineered seals need last only as long as the time required for this salinity gradient to reset itself. One of the limitations of our previous model was the use of a static boundary condition on the borehole wall. The model used a boundary pressure which varied quadratically with depth (arising from differences between the pressure of a column of fresh water in the borehole and that of a column of brine in the host rock). However as brine replaces fresh water in the borehole (driven by a pressure gradient), the boundary function must change with time. We now introduce an improved model which takes this time dependent boundary condition into account. The model also takes into account the time taken for a mixture of brine and fresh water to re-establish chemical equilibrium through the process of diffusion. The paper contains the mathematical details of our new iterative model as well as results showing the time taken to reach steady state, and concentration profiles for the components of the brine-filled borehole as a function of time together with a discussion on the implications of the results for developing a post-closure safety assessment for DBD. (authors)
For economic reasons, the US nuclear industry is renewing efforts to build a technical basis to extend rod average burnup limits above the current regulatory burnup limit of 62 GWd/MTU. The primary driver is to increase pressurized water reactor cycle lengths to 24 months, reducing the number of fresh fuel assemblies and core design constraints, thereby making core energy utilization more efficient. However, fuel pellet fragmentation and pulverization, termed high burnup fuel fragmentation (HBFF), has been observed in the high burnup (>90 GWd/MTU) Halden loss-of-coolant-accident (LOCA) integral test series. The issue gained attention when fuel fragmentation and pulverization were also observed closer to the current US regulatory limit during the US Nuclear Regulatory Commission (NRC) sponsored out-of-core integral test at Studsvik Nuclear in early 2011. This led to NRC concerns with potential changes to fuel and core designs relative to fuel pellet pulverization. In a letter to the NRC Commissioners, the staff specifically identified a need to “…define the boundary of safe operation for key fuel design and operating parameters,” stating that “the staff is challenged to evaluate the acceptability of future fuel design advancements and fuel utilization changes.” As such, it can be concluded that HBFF and potential dispersal into the reactor coolant system introduces additional complications in light-water reactor (LWR) fuel safety evaluations. However, it is not clear how much fuel will be susceptible to HBFF; nor has there been a methodology developed to evaluate fuel susceptibility to HBFF. To that end, this paper proposes an analysis methodology to assess fuel susceptibility to HBFF during LOCA scenarios. The work presented here uses the BISON fuel performance code to evaluate a representative pressurized water reactor fuel rod exposed to a rod average burnup of 75 GWd/MTU. Sensitivity studies investigated the impact of the peak cladding temperature, transient fission gas released, and pre-transient fission gas release on cladding ballooning and burst timing. Subsequently, a methodology to assess fuel susceptibility to HBFF will be developed based on experimental data published in the open literature. The methodology will then be demonstrated by calculating the mass of fuel susceptibility to HBFF. The BISON results conclude that increasing peak cladding temperature drastically decreased time to failure, and decreased balloon size both of which have been confirmed experimentally. Additionally, the effect of pre-transient and transient fission gas release affected cladding balloon size and burst timing. Finally, fuel susceptibility to HBFF significantly decreased as a function of peak cladding temperature.