Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sensitivity study”

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 829 records · Page 46

Load Mitigation on Floating Offshore Wind Turbines with Advanced Controls and Tuned Mass Dampers

General Electric, the National Renewable Energy Laboratory (NREL), the University of Massachusetts Amherst (UMass), and Glosten have recently completed a US Department of Energy (DOE)-funded research program to study technologies for mitigating loads on floating offshore wind turbines through the use of advanced turbine controls and tuned mass dampers (TMDs). The analysis was based upon the Glosten PelaStar tension leg platform (TLP) with GE Haliade 150 turbine, a system developed in a previous front end engineering design (FEED) study funded by the Energy Technology Institute (ETI) in the UK. The platform was designed for the WaveHub wave energy research site, with a mean water depth of 59-m. Loads were analyzed by running time-domain simulations in four 50-year return period (50-YRP) ultimate load state (ULS) conditions and 77 fatigue load state (FLS) environmental conditions. In 50-YRP conditions advanced controls are not active. The influence of TMDs on ULS loads have been reported previously (Park et al. [2]). In FLS conditions advanced controls and TMDs afford dramatic reductions in fatigue damage, offering the potential of significant savings in tower structural requirements. Simulations in turbine idling conditions were run in OrcaFlex, and simulations in operating conditions were run in FAST. Simulations were run with a baseline turbine controller, representative of the current state of the art, and an advanced controller developed by NREL to use collective and individual blade pitch control to maintain rotor speed and reduce tower loads. UMass developed a number of TMD types, with varying system configurations, including passive nonlinear dampers and semi-actively controlled dampers with an inverse velocity groundhook control algorithm. Loads and accelerations in FLS conditions were evaluated on the basis of damage equivalent loads (DELs), and fatigue damage was computed by Miner's summations of stress cycles at the tower base. To study sensitivity to water depth, loads were analyzed at both the 59-m WaveHub depth and a more commercially realistic depth of 100 m. TMDs reduce fatigue damage at the tower-column interface flange by up to 52% in 59-m water depth and up to 28% in 100 m water depth. Advanced controls reduce fatigue damage at the tower-column flange by up to 22% in 59-m water depth and up to 40% in 100 m water depth. The most effective load-mitigation strategy is combining advanced controls with TMDs. This strategy affords a 71% reduction in fatigue damage in both 59-m and 100-m water depths.

17 WIND ENERGY↗

Effects of the shape of the foil corners on the irradiation performance of U10Mo alloy based monolithic mini-plates

Monolithic plate-type fuel is a fuel form being developed for high performance research and test reactors to minimize the use of enriched material. These fuel elements are comprised of a high density, low enrichment, U-Mo alloy based fuel foil, sandwiched between Zirconium liners and encapsulated in Aluminum cladding. The use of a high density fuel in a foil form presents a number of fabrication and operational concerns, such as: foil centering, flatness of the foil, fuel thickness variation, geometrical tilting, foil corner shape etc. To benchmark this new design, effects of various geometrical and operational variables on irradiation performance have been evaluated. As a part of these series of sensitivity studies, the shape of the foil corners were studied. To understand the effects of the corner shapes of the foil on thermo-mechanical performance of the plates, a behavioral model was developed for a selected plate from RERTR-12 experiments (Plate L1P785). Both fabrication and irradiation processes were simulated. Once the thermo-mechanical behavior the plate is understood for the nominal case, the simulations were repeated for two additional corner shapes to observe the changes in temperature, displacement and stress-strain fields. The results from the fabrication simulations indicated that the foil corners do not alter the post-fabrication stress-strain magnitudes. Furthermore, the irradiation simulations revealed that post-fabrication stresses of the foil would be relieved very quickly in operation. While, foils with chamfered and filleted corners yielded stresses with comparable magnitudes, they are slightly lower in magnitudes, and provided a more favorable mechanical response compared with the foil with sharp corners.

Finite Element Analysis↗

WE-E-17A-01: Characterization of An Imaging-Based Model of Tumor Angiogenesis

Purpose: Understanding the transient dynamics of tumor oxygenation is important when evaluating tumor-vasculature response to anti-angiogenic therapies. An imaging-based tumor-vasculature model was used to elucidate factors that affect these dynamics. Methods: Tumor growth depends on its doubling time (Td). Hypoxia increases pro-angiogenic factor (VEGF) concentration which is modeled to reduce vessel perfusion, attributing to its effect of increasing vascular permeability. Perfused vessel recruitment depends on the existing perfused vasculature, VEGF concentration and maximum VEGF concentration (VEGFmax) for vessel dysfunction. A convolution-based algorithm couples the tumor to the normal tissue vessel density (VD-nt). The parameters are benchmarked to published pre-clinical data and a sensitivity study evaluating the changes in the peak and time to peak tumor oxygenation characterizes them. The model is used to simulate changes in hypoxia and proliferation PET imaging data obtained using [Cu- 61]Cu-ATSM and [F-18]FLT respectively. Results: Td and VD-nt were found to be the most influential on peak tumor pO2 while VEGFmax was marginally influential. A +20 % change in Td, VD-nt and VEGFmax resulted in +50%, +25% and +5% increase in peak pO2. In contrast, Td was the most influential on the time to peak oxygenation with VD-nt and VEGFmax playing marginal roles. A +20% change in Td, VD-nt and VEGFmax increased the time to peak pO2 by +50%, +5% and +0%. A −20% change in the above parameters resulted in comparable decreases in the peak and time to peak pO2. Model application to the PET data was able to demonstrate the voxel-specific changes in hypoxia of the imaged tumor. Conclusion: Tumor-specific doubling time and vessel density are important parameters to be considered when evaluating hypoxia transients. While the current model simulates the oxygen dynamics of an untreated tumor, incorporation of therapeutic effects can make the model a potent tool for analyzing anti-angiogenic therapies.

60 APPLIED LIFE SCIENCES↗

The s process in massive stars, a benchmark for neutron capture reaction rates

A clear definition of the contribution from the slow neutron-capture process (s process) to the solar abundances between Fe and the Sr-Zr region is a crucial challenge for nuclear astrophysics. Robust s-process predictions are necessary to disentangle the contribution from other stellar processes producing elements in the same mass region. Nuclear uncertainties are affecting s-process calculations, but most of the needed nuclear input are accessible to present nuclear experiments or they will be in the near future. Neutron-capture rates have a great impact on the s process in massive stars, which is a fundamental source for the solar abundances of the lighter s-process elements heavier than Fe (weak s-process component). In this work we present a new nuclear sensitivity study to explore the impact on the s process in massive stars of 86 neutron-capture rates, including all the reactions between C and Si and between Fe and Zr. We derive the impact of the rates at the end of the He-burning core and at the end of the C-burning shell, where the 22 Ne(α,n) 25 Mg reaction is is the main neutron source. We confirm the relevance of the light isotopes capturing neutrons in competition with the Fe seeds as a crucial feature of the s process in massive stars. For heavy isotopes we study the propagation of the neutron-capture uncertainties, finding a clear difference of the impact of Fe and Co isotope rates with respect to the rates of heavier stable isotopes. The local uncertainty propagation due to the neutron-capture rates at the s-process branching points is also considered, discussing the example of 85 Kr. The complete results of our study for all the 86 neutron-capture rates are available online. Finally, we present the impact on the weak s process of the neutron-capture rates included in the new ASTRAL library (v0.2).

79 ASTRONOMY AND ASTROPHYSICS↗

The CUPID-Mo experiment for neutrinoless double-beta decay: performance and prospects

CUPID-Mo is a bolometric experiment to search for neutrinoless double-beta decay ($0\nu\beta\beta$) of $^{100}$Mo. In this article, we detail the CUPID-Mo detector concept, assembly, installation in the underground laboratory in Modane in 2018, and provide results from the first datasets. The demonstrator consists of an array of 20 scintillating bolometers comprised of $^{100}$Mo-enriched 0.2 kg Li$_2$MoO$_4$ crystals. The detectors are complemented by 20 thin cryogenic Ge bolometers acting as light detectors to distinguish $\alpha$ from $\gamma$/$\beta$ events by the detection of both heat and scintillation light signals. We observe good detector uniformity, facilitating the operation of a large detector array as well as excellent energy resolution of 5.3 keV (6.5 keV) FWHM at 2615 keV, in calibration (physics) data. Based on the observed energy resolutions and light yields a separation of $\alpha$ particles at much better than 99.9\% with equally high acceptance for $\gamma$/$\beta$ events is expected for events in the region of interest for $^{100}$Mo $0\nu\beta\beta$. We present limits on the crystals' radiopurity ($\leq$3 $\mu$Bq/kg of $^{226}$Ra and $\leq$2 $\mu$Bq/kg of $^{232}$Th). Based on these initial results we also discuss a sensitivity study for the science reach of the CUPID-Mo experiment, in particular, the ability to set the most stringent half-life limit on the $^{100}$Mo $0\nu\beta\beta$ decay after half a year of livetime. In conclusion, the achieved results show that CUPID-Mo is a successful demonstrator of the technology - developed in the framework of the LUMINEU project - selected for the CUPID experiment, a proposed follow-up of CUORE, the currently running first tonne-scale cryogenic $0\nu\beta\beta$ experiment.

100Mo↗

Simulating the Impact of Dynamic Rerouting on Metropolitan-scale Traffic Systems

The rapid introduction of mobile navigation aides that use real-time road network information to suggest alternate routes to drivers is making it more difficult for researchers and government transportation agencies to understand and predict the dynamics of congested transportation systems. Computer simulation is a key capability for these organizations to analyze hypothetical scenarios; however, the complexity of transportation systems makes it challenging for them to simulate very large geographical regions, such as multi-city metropolitan areas. In this article, we describe enhancements to the Mobiliti parallel traffic simulator to model dynamic rerouting behavior with the addition of vehicle controller actors and vehicle-to-controller reroute requests. The simulator is designed to support distributed-memory parallel execution using discrete event simulation and be scalable on high-performance computing platforms. We demonstrate the potential of the simulator by analyzing the impact of varying the population penetration rate of dynamic rerouting on the San Francisco Bay Area road network. Using high-performance parallel computing, we can simulate a day in the San Francisco Bay Area with 19 million vehicle trips with 50 percent dynamic rerouting penetration over a road network with 0.5 million nodes and 1 million links in less than three minutes. We present a sensitivity study on the dynamic rerouting parameters, discuss the simulator’s parallel scalability, and analyze system-level impacts of changing the dynamic rerouting penetration. Furthermore, we examine the varying effects on different functional classes and geographical regions and present a validation of the simulation results compared to real-world data.

97 MATHEMATICS AND COMPUTING↗

Sensitivity of surface temperature to oceanic forcing via q-Flux Green’s function experiments. Part III: asymmetric response to warming and cooling

Climate response is often assumed to be linear in climate sensitivity studies. However, by examining the surface temperature (TS) response to pairs of oceanic forcings of equal amplitude but opposite sign in a large set of local q-flux perturbation experiments with CAM5 coupled to a slab, we find strong asymmetry in TS responses to the heating and cooling forcings, indicating a strong nonlinearity intrinsic to the climate system examined. Regardless of where the symmetric forcing is placed, the cooling response to the negative forcing always exceeds the warming to the positive forcing, implying an intrinsic inclination towards cooling of our current climate. Thus, the ongoing global warming induced by increasing greenhouse gases may have already been alleviated by the asymmetric component of the response. The common asymmetry in TS response peaks in high latitudes, especially along sea ice edges, with notable seasonal dependence. Decomposition into different radiative feedbacks through a radiative kernel indicates that the asymmetry in the TS response is realized largely through lapse rate and albedo feedbacks. Further process interference experiments disabling the seasonal cycle and/or sea ice reveal that the asymmetry originates ultimately from the presence of the sea ice component and is further amplified by the seasonal cycle. The fact that a pair of opposite tropical q-flux forcings can excite very similar asymmetric response as a pair placed at 55°S strongly suggests the asymmetric response is a manifestation of an internal mode of the climate model system.

Liu, Fukai↗

A Bimodal Diagnostic Cloud Fraction Parameterization. Part I: Motivating Analysis and Scheme Description

Cloud fraction parameterizations are beneficial to regional, convection-permitting numerical weather prediction. For its operational regional midlatitude forecasts, the Met Office uses a diagnostic cloud fraction scheme that relies on a unimodal, symmetric subgrid saturation-departure distribution. This scheme has been shown before to underestimate cloud cover and hence an empirically based bias correction is used operationally to improve performance. This first of a series of two papers proposes a new diagnostic cloud scheme as a more physically based alternative to the operational bias correction. The new cloud scheme identifies entrainment zones associated with strong temperature inversions. Additionally, for model grid boxes located in this entrainment zone, collocated moist and dry Gaussian modes are used to represent the subgrid conditions. The mean and width of the Gaussian modes, inferred from the turbulent characteristics, are then used to diagnose cloud water content and cloud fraction. It is shown that the new scheme diagnoses enhanced cloud cover for a given gridbox mean humidity, similar to the current operational approach. It does so, however, in a physically meaningful way. Using observed aircraft data and ground-based retrievals over the southern Great Plains in the United States, it is shown that the new scheme improves the relation between cloud fraction, relative humidity, and liquid water content. An emergent property of the scheme is its ability to infer skewed and bimodal distributions from the large-scale state that qualitatively compare well against observations. A detailed evaluation and resolution sensitivity study will follow in Part II.

54 ENVIRONMENTAL SCIENCES↗

Learning curves for drug response prediction in cancer cell lines

Motivated by the size and availability of cell line drug sensitivity data, researchers have been developing machine learning (ML) models for predicting drug response to advance cancer treatment. As drug sensitivity studies continue generating drug response data, a common question is whether the generalization performance of existing prediction models can be further improved with more training data. We utilize empirical learning curves for evaluating and comparing the data scaling properties of two neural networks (NNs) and two gradient boosting decision tree (GBDT) models trained on four cell line drug screening datasets. The learning curves are accurately fitted to a power law model, providing a framework for assessing the data scaling behavior of these models. The curves demonstrate that no single model dominates in terms of prediction performance across all datasets and training sizes, thus suggesting that the actual shape of these curves depends on the unique pair of an ML model and a dataset. The multi-input NN (mNN), in which gene expressions of cancer cells and molecular drug descriptors are input into separate subnetworks, outperforms a single-input NN (sNN), where the cell and drug features are concatenated for the input layer. In contrast, a GBDT with hyperparameter tuning exhibits superior performance as compared with both NNs at the lower range of training set sizes for two of the tested datasets, whereas the mNN consistently performs better at the higher range of training sizes. Moreover, the trajectory of the curves suggests that increasing the sample size is expected to further improve prediction scores of both NNs. These observations demonstrate the benefit of using learning curves to evaluate prediction models, providing a broader perspective on the overall data scaling characteristics. A fitted power law learning curve provides a forward-looking metric for analyzing prediction performance and can serve as a co-design tool to guide experimental biologists and computational scientists in the design of future experiments in prospective research studies.

60 APPLIED LIFE SCIENCES↗

A New Equivalence Theory Method for Treating Doubly Heterogeneous Fuel - I. Theory

A new methodology has been developed to treat resonance self-shielding in doubly heterogeneous very high temperature gas-cooled reactor systems in which the fuel compact region of a reactor lattice consists of small fuel grains dispersed in a graphite matrix. This new method first homogenizes the fuel grain and matrix materials using an analytically derived disadvantage factor from a two-region problem with equivalence theory and intermediate resonance method. This disadvantage factor accounts for spatial self-shielding effects inside each grain within the framework of an infinite array of grains. Then the homogenized fuel compact is self-shielded using a Bondarenko method to account for interactions between the fuel compact regions in the fuel lattice. In the final form of the equations for actual implementations, the double-heterogeneity effects are accounted for by simply using a modified definition of a background cross section, which includes geometry parameters and cross sections for both the grain and fuel compact regions. With the new method, the doubly heterogeneous resonance self-shielding effect can be treated easily even with legacy codes programmed only for a singly heterogeneous system by simple modifications in the background cross section for resonance integral interpolations. This paper presents a detailed derivation of the new method and a sensitivity study of double-heterogeneity parameters introduced during the derivation. The implementation of the method and verification results for various test cases are presented in the companion paper.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

An economics-by-design approach to a radiant integrated thermophotovoltaic microreactor system

This paper presents an application of the economics-by-design approach to the Radiant Integrated Thermo-photovoltaic Microreactor System (RITMS). The RITMS design is unique in that it directly couples a critical fission reactor with thermo-photovoltaic (TPV) panels for high efficiency energy conversion. This significant shift from electric conversion using traditional dynamic heat pumps leads to a simpler and more reliable system without turbomachinery and high pressurization. In working towards wrapping up the early design work, the economics-by-design approach, which centers economic competitiveness as the optimization parameter, is well suited in making a final determination on viability. This paper describes the computational sequence that was developed to couple the radiative and conductive heat transfer and feed operational performance to cost estimation. The framework was applied to maximize the power of the system, while minimizing the fuel enrichment. This method is applied to a reference RITMS design as part of a parametric sensitivity study, which revealed that the system is under moderated and that single unit plants could produce power as low as 300 dollar/MWh. This price point supports the notion that early in the RITMS design implementation, adoption into niche markets as a first-of-a-kind technology is possible. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Kimberlina - Regional Static 3D Geologic Model of the Southern San Joaquin Basin

Regional static 3D geologic model of the southern San Joaquin Basin for the Westcarb Kimberlina demonstration project created in Earthvision (DGI, Inc.). Also contains a smaller-scale model from the larger basin domain for the geometry used in Birkholzer et al. (2011). Data sets derived from: Wagoner, J., 2009, 3D geologic modeling of the Southern San Joaquin Basin for the Westcarb Kimberlina demonstration project—A status report: Lawrence Livermore National Laboratory LLNL-TR-410813, 27 p. Birkholzer, J.T., Zhou, Q., Cortis, A. and Finsterle, S., 2011, A sensitivity study on regional pressure buildup from large-scale CO2 storage projects. Energy Procedia, 4, pp.4371-4378. LLNL contact and file inventory creator: Kayyum Mansoor (mansoor1@llnl.gov)

3D Geological Model↗

Hydraulic-Fracture-Height Growth Under the Combined Influence of Stress Barriers and Natural Fractures

A fully coupled finite-element/finite-volume code is used here to model 3D hydraulically driven fractures under the influence of strong vertical variations in closure stress interacting with natural fractures. Previously unknown 3D interaction mechanisms on fracture-height growth are revealed. Slipping of a natural fracture, triggered by elevated fluid pressure from an intersecting hydraulic fracture, can induce both increases and decreases of normal stress in the minimum-horizontal-stress direction, toward the center and tip of the natural fracture, respectively. Consequently, natural fractures are expected to be able to both encourage and inhibit the progress of hydraulic fractures propagating through stress barriers, depending on the relative locations between the intersecting fractures. Once the hydraulic fracture propagates above the stress barrier through the weakened segment near a favorably located natural fracture, a configuration consisting of two opposing fractures cuts the stress barrier from above and below. The fluid pressure required to break the stress barrier under such opposing-fracture configurations is substantially lower than that required by a fracture penetrating the same barrier from one side. Sensitivity studies of geologic conditions and operational parameters have also been performed to explore the feasibility of controlled fracture height. The interactions between hydraulic fractures, natural fractures, and geologic factors such as stress barriers in three dimensions are shown to be much more complex than in two dimensions. Although it is impossible to exhaust all the possible configurations, the ability of a 3D, fully coupled numerical model to naturally capture these processes is well-demonstrated.

42 ENGINEERING↗

Low-Energy Solvents For Carbon Dioxide Capture Enabled By A Combination Of Enzymes And Vacuum Regeneration

An integrated bench-scale system combining the attributes of the bio-renewable enzyme carbonic anhydrase (CA) with low-enthalpy CO 2 absorption solvents and vacuum regeneration was designed, built and operated for 500 hours using simulated flue gas. The objective was to develop a CO 2 capture process with improved efficiency and sustainability when compared to NETL Case 10 monoethanolamine (MEA) scrubbing technology. The use of CA accelerates inter-conversion between dissolved CO 2 and bicarbonate ion to enhance CO 2 absorption, and the use of low enthalpy CO 2 absorption solvents makes it possible to regenerate the solvent at lower temperatures relative to the reference MEA-based solvent. The vacuum regeneration-based integrated bench-scale system operated successfully for an accumulated 500 hours using aqueous 23.5 wt% K 2 CO 3 -based solvent containing 2.5 g/L enzyme to deliver an average 84% CO 2 capture when operated with a 20% enzyme replenishment rate per ~7 hour steady-state run period. The total inlet gas flow was 30 standard liters per minute with 15% CO 2 and 85% N 2 . The absorber temperature was 40°C and the stripper operated under 35 kPa pressure with an approximate 77°C stripper bottom temperature. Tests with a 30°C absorber temperature delivered >90% capture. On- and off-line operational measurements provided a full process data set, with recirculating enzyme, that allowed for enzyme replenishment and absorption/desorption kinetic parameter calculations. Dissolved enzyme replenishment and conventional process controls were demonstrated as straightforward approaches to maintain system performance. Preliminary evaluation of a novel flow-through ultrasonically enhanced regeneration system was also conducted, yet resulted in CO2 release within the range of temperature-dependent release, and further work would be needed to validate the benefits of ultrasonic enhanced stripping. A full technology assessment was completed in which four techno-economic cases for enzyme-enhanced aqueous K 2 CO 3 solvent with vacuum stripping were considered and a corresponding set of sensitivity studies were developed. The cases were evaluated using bench-scale and laboratory-based observations, AspenPlus® process simulation and modeling, AspenTech’s CCE® Parametric Software, current vendor quotations, and project partners’ know-how of unit operations. Overall, the DOE target of 90% CO 2 capture could be met using the benign enzyme-enhanced aqueous K 2 CO 3 -based alternative to NETL Case 10. The model-predicted plant COE performance, scaled to 550 MWe net output, was 9% higher than NETL Case 10 for an enzyme-activated case with minimized technical risk and highest confidence in physical system performance utilizing commercially available equipment. A COE improvement of 2.8% versus NETL Case 10 was predicted when favorable features of improved enzyme longevity and additional power output from a very low pressure (VLP) turbine were combined, wherein corresponding high capital and operational costs limited the level of COE benefit. The environmental, health and safety (EH&S) profile of the system was found to be favorable and was compliant with the Federal EH&S legislation reviewed. Further work on a larger scale test unit is recommended to reduce the level of uncertainty inherent in extrapolating findings from a bench-scale unit to a full scale PCC plant, and to further investigate several identified opportunities for improvement. Production feasibility and suitability of carbonic anhydrases for scale-up testing was confirmed both through the current project and through parallel efforts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nuclear Energy Advanced Modeling and Simulation (NEAMS) Accident Tolerant Fuels High Impact Problem: FeCrAl Modeling Capabilities

Iron-chromium-aluminum alloys are a potential candidate replacement for zirconium-based claddings in light water reactor fuel rods due to their increased strength (lower creep rates and higher yield stress) and slower oxidation kinetics with steam. Their larger neutron absorption cross-section necessitates thinner claddings and larger fuel pellets to maintain the same as-fabricated fuel-to-clad gap. The aggressive development schedule of the United States Department of Energy's (DOE) Advanced Fuels Campaign (AFC) to have potential accident tolerant fuels (ATF) concepts in lead test rods or assemblies by 2022 has necessitated the need for advanced modeling and simulation techniques to investigate the behavior of these materials under normal operating and accident conditions. This report presents the modeling capabilities for FeCrAl claddings at the engineering (continuum) scale developed as part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) high impact problem (HIP). Sensitivity studies of select material parameters have been completed to determine their significance in predictions of important fuel performance parameters such as fuel centerline and surface temperatures, cladding stresses and strains, fission gas released, and rod internal pressure. Comparative analyses of FeCrAl and Zircaloy-4 claddings for separate effects conditions are presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A status review on DNB prediction using CASL baseline boiling model and possible suggestions regarding wall boiling closures

The main goal of FY18 THM’s focus task is to evaluate where the current CASL boiling model stands for DNB predictive performance, to identify which components should be improved or calibrated to enhance DNB performance. Over the past three year’s validation tests, Gen1A boiling model has presented an under-predictive behavior for DNB at the high pressure of interest for PWR application. With this in mind, this year’s task is to focus on 1) best practice guidelines on DNB detection criteria, 2) interfacial momentum transfer closure tests, and 3) wall boiling closure sensitivity on DNB. Best practice guidelines on DNB detection criteria are demonstrated by the following three monitoring variables: 1) max wall temperature, 2) K-value and 3) energy balance. Interfacial momentum transfer closure govern the bubble distribution, especially in lateral direction, as well as bubble mean diameter. A clear correlation between interfacial momentum closure and DNB prediction is not observed in the current tests. However, void distributions are changed with different interfacial momentum closure tests which requires further validation with corresponding experimental measurement. From a variety of wall boiling closure tests, bubble departure diameter seems the most sensitive for DNB performance while fulfilling numerical robustness in DNB simulation. Tolubinski’s bubble departure diameter model was calibrated to see if any improvement is observed in DNB performance. Notable DNB improvement was observed as the bubble departure diameter is reduced. However, reduced bubble departure diameter by calibrated Tolubinski’s model seemed to lead to relatively high wall superheat temperature prediction which requires further validation with associated experimental measurement It should be noted that lack of experimental measurement close to DNB condition hampers this boiling closure sensitivity study significantly. This situation directly reflects some of limitations of the CASL boiling model development. With this in mind, the implementation of improved boiling closure on 5x5 NMVG and MVG tests will be pursued in FY19, the associated DNB best practice guidelines and possible limitations will be addressed for industry DNB application.

42 ENGINEERING↗

Development of Preliminary VERA-CS Crud Induced Localized Corrosion Modeling Capability

The objective of this milestone was to develop a crud-induced localized corrosion (CILC) screening capability in VERA-CS using existing codes and capabilities that had been developed in the Consortium for Advanced Simulation of Light Water Reactors (CASL) as part of the efforts to resolve the crud challenge problem. The milestone work focused on improving existing CASL software and implementing clad corrosion modeling into VERA-CS. Clad corrosion modeling was implemented by creating an application interface (API) in the Cicada software so that its clad heat transfer and corrosion models could be used in VERA-CS. Cicada was coupled to the subchannel code CTF for solving crud and corrosion growth. Additionally, a clad corrosion model was added directly to CTF to serve as an alternate approach. The rod thermal hydraulic reconstruction (ROTHCON) process, which is used for capturing computational fluid dynamics (CFD)-predicted rod-surface heat and turbulence behavior that can impact crud-induced power shift (CIPS) and CILC prediction, was also implemented into VERA-CS. Usability improvements were made to this feature, which was developed over several previous milestones, so that the VERA-CS user can link spacer grids in a VERA-CS model to an applicable pregenerated spacer grid data file if available. Several assessments were performed to demonstrate the proper functioning of the new capability. Sensitivity studies were performed to demonstrate the effect of crud lithium content, oxide thermal conductivity, Cicada mesh refinement, and ROTHCON coupling mesh refinement on predicted amounts of crud and oxide growth. A full-scale demonstration was performed by running Seabrook Cycle 5 with the CILC capability. Some general observations were made from the assessments. It was observed that predicted corrosion thickness is less than anticipated for CILC cases, so further work is needed to better understand and improve the corrosion modeling. Refining the coupling mesh in the ROTHCON process leads to increased corrosion and crud thickness, which can affect where crud and corrosion grow, how thick it gets, and how many rods experience crud and corrosion growth. Results indicate that applying the ROTHCON process and refining the coupling mesh leads to only a slight increase in total core maximum crud and corrosion thickness, but it leads to a significant increase in the number of rods with a maximum thickness that surpasses a set threshold value (i.e., more rods grow thicker corrosion and crud). Computational performance of the simulation is degraded by a factor of 5{6 with high coupling mesh refinement. The impact of the ROTHCON coupling mesh refinement on total core crud mass is less significant; however, using the ROTHCON process with CFD-generated data over the Yao{Hochreiter{Leech (YHL) model in CTF leads to significantly different amounts of total core crud mass. Future work on the CILC screening capability will include model improvement and validation, development and testing of a CILC-risk screening criteria, and performance improvements.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Preliminary Evaluation of FeCrAl Cladding for Accident Tolerant Fuel Concepts

Since the accident at the Fukushima Daiichi Nuclear Power Station, enhancing the accident tolerance of light water reactors (LWRs) has become an important research topic. In particular, the community is actively developing enhanced fuels and cladding for LWRs to improve safety in the event of accidents in the reactor or spent fuel pools. Fuels with enhanced accident tolerance are those that, in comparison with the standard UO2-zirconium alloy system, can tolerate loss of active cooling in the reactor core for a considerably longer time period during design-basis and beyond design-basis events while maintaining or improving the fuel performance during normal operations and operational transients. This report presents early work in developing thermal and mechanical models for FeCrAl, a proposed cladding material. BISON, the finite element-based fuel performance code in development at Idaho National Laboratory, was used to compare results from normal operation conditions with UO2/Zr-4 behavior. In addition, sensitivity studies are presented for evaluating the relative importance of material parameters such as elastic constants and thermal conductivity in order to provide guidance on future experiments for this material.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗