Search NASA⌕ Search

SEARCH · Search NASA

Results for “NUCLEATION PROCESS”

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 217 records · Page 12

A new method for operating a continuous-flow diffusion chamber to investigate immersion freezing: assessment and performance study

Abstract. Glaciation in mixed-phase clouds predominantly occurs through the immersion-freezing mode where ice-nucleating particles (INPs) immersed within supercooled droplets induce the nucleation of ice. Model representations of this process currently are a large source of uncertainty in simulating cloud radiative properties, so to constrain these estimates, continuous-flow diffusion chamber (CFDC)-style INP devices are commonly used to assess the immersion-freezing efficiencies of INPs. This study explored a new approach to operating such an ice chamber that provides maximum activation of particles without droplet breakthrough and correction factor ambiguity to obtain high-quality INP measurements in a manner that previously had not been demonstrated to be possible. The conditioning section of the chamber was maintained at −20 ∘C and water relative humidity (RHw) conditions of 113 % to maximize the droplet activation, and the droplets were supercooled with an independently temperature-controlled nucleation section at a steady cooling rate (0.5 ∘C min−1) to induce the freezing of droplets and evaporation of unfrozen droplets. The performance of the modified compact ice chamber (MCIC) was evaluated using four INP species: K-feldspar, illite-NX, Argentinian soil dust, and airborne soil dusts from an arable region that had shown ice nucleation over a wide span of supercooled temperatures. Dry-dispersed and size-selected K-feldspar particles were generated in the laboratory. Illite-NX and soil dust particles were sampled during the second phase of the Fifth International Ice Nucleation Workshop (FIN-02) campaign, and airborne soil dust particles were sampled from an ambient aerosol inlet. The measured ice nucleation efficiencies of model aerosols that had a surface active site density (ns) metric were higher but mostly agreed within 1 order of magnitude compared to results reported in the literature.

54 ENVIRONMENTAL SCIENCES↗

The hierarchical energy landscape of screw dislocation motion in refractory high-entropy alloys

High-entropy alloys (HEAs) with concentrated solid solutions are conceived to possess a rugged atomic and energy landscape in which dislocation motion necessarily proceeds to accommodate mechanical deformation. Fundamental questions remain as to how rough the energy landscape is and to what extent it can be influenced by the local ordering of their constituent elements. Here, we construct and report the potential energy landscape (PEL) governing screw dislocation motion in refractory HEAs that reveals a hierarchical and multilevel structure with a collection of small basins nested in large metabasin. This striking feature pertaining to HEAs exerts a trapping force and back stress on saddle point activations, retarding dislocation movement. Further, by introducing chemical short-range order, the energy landscape is smoothed but skewed to different degrees that shifts the rate-liming process from kink-glide to kink-pair nucleation. The chemical disorder-roughened PEL in HEAs, analogous to structural disorder induced in metallic glasses, signifies the role of various barrier-hopping processes underlying the extraordinary mechanical behaviors of the two distinct groups of materials.

36 MATERIALS SCIENCE↗

Atmospheric ice nuclei concentration measurements over a high altitude-station in the Western Ghats, India

Mixed-phase and ice clouds play a significant role in modulating the Indian summer monsoon rainfall, and modeling studies suggest that inaccurate representation of ice nucleating particles (INP) is one of the factors, which have resulted in a deficiency in the prediction of cloud properties and hydrological cycle. We do not understand the ice nucleation mechanisms to represent the ice formation in the model, and therefore more INP measurements at various supercooled temperatures are needed. This study reports for the first time real-time INP concentration at a high altitude station in the Western Ghats region of India, and these results will be used to constrain the ice nucleation parameterizations in the model. The variation of INP for any particular cloud condition can have strong microphysical responses of mixed-phase precipitation formation processes. Also, to improve the representation of heterogeneous ice nucleation, INP measurements of various types of aerosols commonly observed over the Indian region are needed. This study presents, the INP concentrations and daily concentration of non-refractory chemical composition (organics, sulfate, nitrate, ammonium, and chloride), cloud condensation nuclei (CCN), and aerosol size distribution. Spectrometer for Ice Nuclei (SPIN) was operated to measure INP concentration at the mountain site from August to December 2018. INP measurements were performed at three different temperatures (-25, -30, and - 34 °C) in the immersion freezing mode relevant for mixed-phase cloud conditions. The average INP concentrations approximately varied from 0.18 to 12.4 L-1, 0.39 to 24 L-1, 1.1 to 40.2 L-1 at -25, -30, and - 34 °C, respectively. These concentrations are within the concentrations reported in the literature from different regions across the globe. The air mass back trajectory analysis showed that continental air mass contains more INP compared to maritime air mass. The information on INP concentration and chemical characteristics of aerosol will help to improve heterogeneous ice nucleation parameterization in numerical models.

Kumar, Anil V.↗

Molecular-Dynamics Analysis of the Mechanical Behavior of Plasma-Facing Tungsten

We report a systematic computational analysis of the mechanical behavior of plasma-facing component (PFC) tungsten focusing on the impact of void and helium (He) bubble defects on the mechanical response beyond the elastic regime. Specifically, we explore the effects of porosity and He atomic fraction on the mechanical properties and structural response of PFC tungsten, at varying temperature and bubble size. We find that the Young modulus of defective tungsten undergoes substantial softening that follows an exponential scaling relation as a function of matrix porosity and He atomic content. Beyond the elastic regime, our high strain rate simulations reveal that the presence of nanoscale spherical defects (empty voids and He bubbles) reduces the yield strength of tungsten in a monotonically decreasing fashion, obeying an exponential scaling relation as a function of tungsten matrix porosity and He concentration. Our detailed analysis of the structural response of PFC tungsten near the yield point reveals that yielding is initiated by emission of dislocation loops from bubble/matrix interfaces, mainly 1/2<111> shear loops, followed by gliding and growth of these loops and reactions to form <100> dislocations. Furthermore, dislocation gliding on the <111>{211} twin systems nucleates 1/6<111> twin regions in the tungsten matrix. These dynamical processes reduce the stress in the matrix substantially. Subsequent dislocation interactions and depletion of the twin phases via nucleation and propagation of detwinning partials lead the tungsten matrix to a next deformation stage characterized by stress increase during applied straining. Furthermore, our structural analysis reveals that the depletion of twin boundaries (areal defects) is strongly impacted by the density of He bubbles at higher porosities. After the initial stress relief upon yielding, increase in the dislocation density in conjunction with decrease in the areal defect density facilitates the initiation of dislocation-driven deformation mechanisms in the PFC crystal.

36 MATERIALS SCIENCE↗

A 1D Model for Nucleation of Ice From Aerosol Particles: An Application to a Mixed‐Phase Arctic Stratus Cloud Layer

Abstract Mixed‐phase clouds (MPCs) have been identified as significant contributors to uncertainties in climate projections, attributable to model representation of processes controlling the formation and loss of supercooled water droplets and ice particles from the atmosphere. Arctic MPCs are commonly widespread and long‐lived, with sustained ice crystal formation processes that challenge current understanding. This study examines the ice‐nucleating particle (INP) reservoir dynamics governing immersion‐mode heterogeneous freezing in an observed case of Arctic MPCs using a simplified 1D aerosol‐cloud model. The model setup includes prescribed dynamical forcings and thermodynamic profiles, and represents INPs as multicomponent and polydisperse particle size distributions. Diagnostic and prognostic approaches to immersion freezing parameterization are compared, including time‐independent (singular) number‐ and surface area‐based descriptions and a time‐dependent description following classical nucleation theory (CNT). The choice of freezing parameterization defines the size of the INP reservoir. The CNT‐based description yields an orders of magnitude larger INP reservoir than the singular parameterizations, which is the dominant factor for sustained ice crystal formation. The efficiency of the freezing process and cloud cooling are of secondary importance. A diagnostic treatment neglecting INP loss is only accurate when the INP reservoir size is large and INP depletion weak. Since a larger INP reservoir sustains ice crystal formation substantially longer, and ice water path scales with ice crystal concentrations for the conditions considered, resolving the source of differences in INP reservoir dynamics due to model implementation is a high priority for advancing climate model physics.

54 ENVIRONMENTAL SCIENCES↗

Nucleation Rate Theory for Coordination Number: Elucidating Water-Mediated Formation of a Zigzag Na 2 SO 4 Morphology

Predicting and controlling nanostructure formation during nucleation can pave the way to synthesizing novel energy materials via crystallization. However, such control over nucleation and crystallization remains challenging due to an inadequate understanding of critical factors that govern evolving atomistic structures and dynamics. Here, in this paper, we utilize coordination number as a reaction coordinate and rate theory to investigate how sodium sulfate, commonly known as a phase-change energy material, nucleates in a supersaturated aqueous solution. In conjunction with ab initio and force field-based molecular dynamics simulation, the rate theoretical analysis reveals that sodium sulfate from an initially dissolved metastable state transits to a heterogeneous mixture of prenucleated clusters and finally to a large cylindrical zigzag morphology. Measurements of Raman spectra and their ab initio modeling confirm that this nucleated morphology contains a few waters for every sulfate. Rate processes such as solvent exchange and desolvation exhibit high sensitivity to the evolving prenucleation/nucleation structures, providing a means to distinguish between critical nucleation precursors. Desolvation and forming the first-shell interionic coordination structure via monomer-by-monomer addition around sulfates are found to explain the formation of large nuclei. Thus, a detailed understanding of the step-by-step structure formation across scales has been achieved. This can be leveraged to predict nucleation-related structures and dynamics and potentially control the synthesis of novel phase-change materials for energy applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Novel estimation of aerosol processes with particle size distribution measurements: a case study with the TOMAS algorithm v1.0.0

Atmospheric aerosol microphysical processes are a significant source of uncertainty in predicting climate change. Specifically, aerosol nucleation, emissions, and growth rates, which are simulated in chemical transport models to predict the particle size distribution, are not understood well. However, long-term size distribution measurements made at several ground-based sites across Europe implicitly contain information about the processes that created those size distributions. This work aims to extract that information by developing and applying an inverse technique to constrain aerosol emissions as well as nucleation and growth rates based on hourly size distribution measurements. We developed an inverse method based upon process control theory into an online estimation technique to scale aerosol nucleation, emissions, and growth so that the model–measurement bias in three measured aerosol properties exponentially decays. The properties, which are calculated from the measured and predicted size distributions, used to constrain aerosol nucleation, emission, and growth rates are the number of particles with a diameter between 3 and 6 nm, the number with a diameter greater than 10 nm, and the total dry volume of aerosol ( N 3–6 , N 10 , V dry ), respectively. In this paper, we focus on developing and applying the estimation methodology in a zero-dimensional “box” model as a proof of concept before applying it to a three-dimensional simulation in subsequent work. The methodology is first tested on a dataset of synthetic and perfect measurements that span diverse environments in which the true particle emissions, growth, and nucleation rates are known. The inverse technique accurately estimates the aerosol microphysical process rates with an average and maximum error of 2 % and 13 %, respectively. Next, we investigate the effect that measurement noise has on the estimated rates. The method is robust to typical instrument noise in the aerosol properties as there is a negligible increase in the bias of the estimated process rates. Finally, the methodology is applied to long-term datasets of in situ size distribution measurements in western Europe from May 2006 through June 2007. At Melpitz, Germany, and Hyytiälä, Finland, the average diurnal profiles of estimated 3 nm particle formation rates are reasonable, having peaks near noon local time with average peak values of 1 and 0.15 cm –3 s –1 , respectively. The normalized absolute error in estimated N 3–6 , N 10 , and V dry at three European measurement sites is less than 15 %, showing that the estimation framework developed here has potential to decrease model–measurement bias while constraining uncertain aerosol microphysical processes.

54 ENVIRONMENTAL SCIENCES↗

Effect of PVD-coated chromium on the subcooled flow boiling performance of nuclear reactor cladding materials

Here we elucidate the separate effect of a thin Cr coating deposited by physical vapor deposition (PVD) on the subcooled flow boiling performance of zircaloy-4. First, we run flow boiling experiments on prototypical zircaloy-4 surfaces mimicking the scratch pattern and surface roughness of nuclear reactor claddings. Then, we PVD-coat a 0.3 µm thick chromium layer on the same exact surface and repeat the same flow boiling investigations. All experiments are run using deionized water at atmospheric pressure, flowing on a 1 × 3 cm 2 rectangular cross section channel at a rate of 1000 kg/m 2 /s and a subcooling of 10 K. We measure the average temperature of the boiling surface at increasing surface heat fluxes, covering a wide range of heat transfer regimes, from single-phase forced convection to the boiling crisis. We also record high-speed videos of the boiling process, which we postprocess to measure bubble nucleation site density, growth time, departure diameter and frequency. The surface analysis reveals that, while the chromium coating does not seem change the surface roughness and morphology, it improves surface wettability. However, it decreases the critical heat flux. The chromium coating causes an increase of nucleation temperature, bubble departure diameter and growth time, and a reduction of the nucleation site density. The concurrence of these observations indicates that a size reduction of the nucleation sites, conformally covered by the chromium coating, may be the cause of the boiling performance deterioration. We confirm this hypothesis repeating the same analysis on a FeCrAl sample prepared and tested using the same protocol as the zircaloy-4 sample, but with a different initial surface texture.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

An Efficient Reactive Force Field without Explicit Coordination Dependence for Studying Caustic Aluminum Chemistry

Reactive force fields (RFFs) are an expedient approach to sample chemical reaction paths in complex systems, relative to density functional theory. However, there is continued need to improve efficiencies, specifically in systems that have slow transverse degrees of freedom, as in highly viscous and superconcentrated solutions. Here, we present an RFF that is differentiated from current models (e.g., ReaxFF) by omitting explicit dependence on the atom coordination and employing a small parameter set based on Lennard-Jones, Gaussian, and Stillinger–Weber potentials. In this study, the model was parametrized from AIMD simulation data and is used to model aluminate reactivity in sodium hydroxide solutions with extensive validation against experimental radial distribution functions, computed free energy profiles for oligomerization, and formation energies. In conclusion, the model enables simulation of early stage Al(OH) 3 nucleation which has significant relevance to industrial processing of aluminum and has a computational cost that is reduced by 1 order of magnitude relative to ReaxFF.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A unifying criterion of the boiling crisis

Abstract We reveal and justify, both theoretically and experimentally, the existence of a unifying criterion of the boiling crisis. This criterion emerges from an instability in the near-wall interactions of bubbles, which can be described as a percolation process driven by three fundamental boiling parameters: nucleation site density, average bubble footprint radius and product of average bubble growth time and detachment frequency. Our analysis suggests that the boiling crisis occurs on a well-defined critical surface in the multidimensional space of these parameters. Our experiments confirm the existence of this unifying criterion for a wide variety of boiling surface geometries and textures, two boiling regimes (pool and flow boiling) and two fluids (water and liquid nitrogen). This criterion constitutes a simple mechanistic rule to predict the boiling crisis, also providing a guiding principle for designing boiling surfaces that would maximize the nucleate boiling performance.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Equal abundance of summertime natural and wintertime anthropogenic Arctic organic aerosols

Aerosols play an important yet uncertain role in modulating the radiation balance of the sensitive Arctic atmosphere. Organic aerosol is one of the most abundant, yet least understood, fractions of the Arctic aerosol mass. Here we use data from eight observatories that represent the entire Arctic to reveal the annual cycles in anthropogenic and biogenic sources of organic aerosol. We show that during winter, the organic aerosol in the Arctic is dominated by anthropogenic emissions, mainly from Eurasia, which consist of both direct combustion emissions and long-range transported, aged pollution. In summer, the decreasing anthropogenic pollution is replaced by natural emissions. These include marine secondary, biogenic secondary and primary biological emissions, which have the potential to be important to Arctic climate by modifying the cloud condensation nuclei properties and acting as ice-nucleating particles. Their source strength or atmospheric processing is sensitive to nutrient availability, solar radiation, temperature and snow cover. Our results provide a comprehensive understanding of the current pan-Arctic organic aerosol, which can be used to support modelling efforts that aim to quantify the climate impacts of emissions in this sensitive region.

58 GEOSCIENCES↗

A Deep Learning Pipeline for Optimizing Large-scale Phase Field Simulations

Phase field (PF) simulations are computationally expensive but remain a key analysis tool to understand the complex mechanisms of additive manufacturing (AM) processes. Each PF simulation-aided analysis requires thousands of node hours on leadership-class supercomputers. One of the main goals of these analyses is the study of microstructure evolution during the build process which begins with the onset of nucleation. Nucleation occurs under certain thermomechanical conditions which are not known a priori and many PF simulations are required to identify ranges of input thermo-mechanical parameters that can result in the onset of nucleation. Since many of the simulations do not result in nucleation, an analysis campaign often ends up wasting tremendous amounts of precious computing resources executing nucleation-absent simulations. The goal of this work is to design and train deep learning models to inform a PF simulation about the likelihood of the occurrence of nucleation in a future simulation time-step based on the state summary over a finite number of past time-steps of a running simulation. If the prediction determines that the running simulation is unlikely to reach nucleation in the allotted time, then its execution is stopped immediately ultimately resulting in vast reduction in wasted computations when accrued over all the PF simulations typically performed in a single or multiple analysis campaign(s). The paper presents the performance of a machine learning pipeline that uses a convolutional neural network (CNN) model to learn an embedding which is then used with a self-attention network to build a multi-task deep learning model to predict the likelihood of nucleation. The model also predicts the input parameters used in a simulation. Performance is compared with a baseline pipeline that uses an off-the-shelf LeNet-5 model to learn the initial embedding. Despite their smaller size, performance results indicate significant improvement in accuracy of the proposed models compared to the larger baseline models.

Kannan, Ramakrishnan {ramki}↗

A dynamic parameterization of sulfuric acid–dimethylamine nucleation and its application in three-dimensional modeling

Sulfuric acid (SA) is a governing gaseous precursor for atmospheric new particle formation (NPF), a major source of global ultrafine particles, in environments studied around the world. In polluted urban atmospheres with high condensation sinks (CSs), the formation of stable SA–amine clusters, such as SA–dimethylamine (DMA) clusters, usually initializes intense NPF events. Coagulation scavenging and cluster evaporation are dominant sink processes of SA–amine clusters in urban atmospheres, yet these loss processes are not quantitatively included in the present parameterizations of SA–amine nucleation. We herein report a parameterization of SA–DMA nucleation, based on cluster dynamic simulations and quantum chemistry calculations, with certain simplifications to greatly reduce the computational costs. Compared with previous SA–DMA nucleation parameterizations, this new parameterization was able to reproduce the dependences of particle formation rates on temperature and CSs. We then incorporated it in a three-dimensional (3-D) chemical transport model to simulate the evolution of the particle number size distributions. Simulation results showed good consistency with the observations in the occurrence of NPF events and particle number size distributions in wintertime Beijing and represented a significant improvement compared to that using a parameterization without coagulation scavenging. Quantitative analysis shows that SA–DMA nucleation contributes significantly to nucleation rates and aerosol population during the 3-D simulations in Beijing (>99 % and >60 %, respectively). These results broaden the understanding of NPF in urban atmospheres and stress the necessity of including the effects of coagulation scavenging and cluster stability in simulating SA–DMA nucleation in 3-D simulations. Representing these processes is thus likely to improve model performance in particle source apportionment and quantification of aerosol effects on air quality, human health, and climate.

Li, Yuyang↗

Unveiling real‐time crystallization with nucleators and thickeners for zinc nitrate hexahydrate as a phase change material

Abstract The primary challenge of salt hydrates as phase change materials (PCMs) is their high degree of supercooling (∆ T Sup ). Past studies have attempted to alleviate ∆ T Sup by incorporating nucleators possessing similar lattice structure, without a thorough analysis of the nucleator/PCM interactions. This work presents a novel in‐operando study to visually observe the crystallization process of zinc nitrate hexahydrate (ZNH) in real‐time using suitable nucleators and a thickener. We have introduced a new dataset where zinc acetate dihydrate (ZnAc.2H) has been employed as an additive to ZNH for thermal energy storage (TES). According to our in‐operando study, ZnAc.2H underwent hydrolysis in molten ZNH to precipitate needle‐like ZnO particles, which acted as nucleators for the PCM. In the absence of a thickener, the crystal propagation of ZNH exhibited a preferential directional tendency, whereas incorporation of thickener made the phase change process uniform along all directions. Self‐hydrolysis of ZnAc.2H generated excess water which caused the undesirable effect of broadening the endothermic peak. The addition of carboxymethyl cellulose (CMC) thickener could restrict the activity of excess water and reduce the broadening of the endothermic peak. Differential scanning calorimetry (DSC) revealed that 5 wt% ZnAc.2H additive and 2 wt% CMC thickener added to ZNH can considerably improve the TES properties of ZNH PCM (∆ H Fusion 125.9 J·g −1 , ∆ T Sup 3.0°C, melting point [M.P.] 29.0°C). Our in‐operando studies can unveil real‐time phase change behaviors to better design PCM systems with desired characteristics.

Chakraborty, Anirban↗

A phase-field model for quasi-dynamic nucleation, growth, and propagation of rate-and-state faults

Despite its critical role in the study of earthquake processes, numerical simulation of the entire stages of fault rupture remains a formidable task. The main challenges in simulating a fault rupture process include the complex evolution of fault geometry, frictional contact, and off-fault damage over a wide range of spatial and temporal scales. Here, in this paper, we develop a phase-field model for quasi-dynamic fault nucleation, growth, and propagation, which features two standout advantages: (i) it does not require any sophisticated algorithms to represent fault geometry and its evolution; and (ii) it allows for modeling fault nucleation, propagation, and off-fault damage processes with a single formulation. Built on a recently developed phase-field framework for shear fractures with frictional contact, the proposed formulation incorporates rate- and state-dependent friction, radiation damping, and their impacts on fault mechanics and off-fault damage. We show that the numerical results of the phase-field model are consistent with those obtained from well-verified approaches that model the fault as a surface of discontinuity, without suffering from the mesh convergence issue in the existing continuous approaches to fault rupture (e.g., the stress glut method). Further, through numerical examples of fault propagation in various settings, we demonstrate that the phase-field approach may open new opportunities for investigating complex earthquake processes that have remained overly challenging for the existing numerical methods.

58 GEOSCIENCES↗

Water Interplays during Dysprosium Electrodeposition in Pyrrolidinium Ionic Liquid: Deconvoluting the Pros and Cons for Rare Earth Metallization

The electrochemical production of rare earth metals (REMs) in ionic liquids (ILs) has received much attention as a promising, sustainable and replacement to the molten salt electrolysis. Water additives have been recently suggested as a strategy to overcome the limitations of the ionic liquid process, but without complete assessment of its overall viability for the REMs production. In this regard, a full investigation of water impact on dysprosium (Dy) electrodeposition in pyrrolidinium triflate (BMPyOTf) ionic liquid was carried out. Water introduction was revealed to involve an interplay of implications on the electrodeposition process, including coordination, speciation, reduction pathways, interfacial dynamics, nucleation and metal stability and purity. In highly dry conditions, the reduction occurs at very negative potentials (-3.3V) in a consecutive pathway, resulting in negligible metal electrodeposition (low rate and efficiency) at the electrode surface. Initial water introduction (<30% molar fraction) leads to partitioning of the Dy complex between water and IL coordinated speciation, giving rise to an additional wave at more positive potentials (-2.4V). Probing and quantification of the heterogeneous Dy speciation by spectroscopic analyses enabled the uncovering of the reduction mechanism and accurate evaluation of the mass transport properties. While lowering the reduction thermodynamics, water addition also improved the nucleation, deposition rate and faradaic efficiency. In spite of these benefits, stripping voltammetric analysis also indicated substantial chemical reactivity of the deposited Dy metal with water additives and/or electrolyte components, under extended timescales. Surface characterization of the obtained product confirmed limited Dy metal stability (mostly oxidized/fluorinated) and purity (~60%) in the studied system. Moreover, high water introduction (>30% molar fraction) was shown to trigger a fast hydrogen evolution reaction (HER), downgrading the robustness of system efficiency. We report the overall impact of water additives seems to engender both promoting and mitigating effects on the Dy electrodeposition process, making the case for more innovative strategies to be sought.

36 MATERIALS SCIENCE↗

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗

A new method for operating a continuous flow diffusion chamber to investigate immersion freezing: assessment and performance study

Glaciation in mixed-phase clouds predominately occurs through the immersion freezing mode where ice nucleating particles (INPs) immersed within supercooled droplets induce nucleation of ice. Currently, model representations of this process are a large source of uncertainty in simulating cloud radiative properties, and to constrain these estimates, continuous flow diffusion chamber (CFDC)-style INP devices are commonly used to assess the immersion freezing efficiencies of INPs. In this study, a new approach was explored to operating such an ice chamber that provides maximum activation of particles without droplet breakthrough and correction factor ambiguity to obtain high-quality INP measurements in a manner that has not been demonstrated as possible previously. The conditioning section of the chamber was maintained at ~ –20 °C and water relative humidity (RH w ) ~ 113 % conditions to maximize the droplet activation, and the droplets were supercooled with an independently temperature-controlled nucleation section at a steady cooling rate (0.5 °C min –1 ) to induce the freezing of droplets and evaporation of unfrozen droplets. The performance of the modified ice chamber was evaluated using four INP species: K-feldspar, illite-NX, Argentinian soil dust, and airborne arable dust that had shown ice nucleation over a wide span of supercooled temperatures. Dry dispersed and size-selected K-feldspar particles were generated in the laboratory. Illite-NX and soil dust particles were sampled during the second phase of the Fifth International Ice Nucleation Workshop (FIN-02) campaign, and airborne arable dust particles were sampled from the aerosol inlet located on the rooftop of the laboratory. The measured ice nucleation efficiencies of model aerosols with a surface active site density (n s ) metric were higher, but mostly agreed within one order of magnitude compared to literature results.

54 ENVIRONMENTAL SCIENCES↗