Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Selective capture”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A direct numerical simulation study for confined non-isothermal jet impingement at moderate nozzle-to-plate distances: Capturing jet-to-ambient density effects

A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Re τ = 505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15 K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Furthermore, radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Understanding Gas Absorption in Multivalent Ionic Liquids via Solute-Solvent Interaction Analyses

Solute-solvent interactions are often evaluated to understand and explain gas solubilities. This information is regularly used to interpret gas solubility and separation within ionic liquids (ILs). Here, we use several different quantum chemical and previous molecular dynamics simulation results to show that strong solute-solvent interactions neither lead to high solubility nor to high selectivity, with respect to other solutes. Instead, the solute-solvent interaction strength with different isolated IL ions is affected by the electrostatic extremes. Also, the solute-solvent interaction with different IL molecules is strongly correlated to the molecular volumes or polarizabilities (within a series of same-anion based ILs).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Factors controlling injection-induced rupture of intersecting faults during geological sequestration of CO 2

This study addresses coupled multiphase fluid flow and geomechanics effects on potential fault activation associated with subsurface CO 2 injection around intersecting faults. An enhanced fault-representation model is used to capture geomechanical responses of two intersecting faults with finite length during CO 2 injection. The faults are embedded in a strike-slip stress regime of a caprock-reservoir-basement system with the faults represented by zero-thickness interfaces with adjacent finite-thickness damage zones. A sensitivity analysis is conducted to study the effect of fault permeability, slip-weakening behavior, well location relative to the orientation of faults, and well placement (the number and location of injection wells). Five metrics (pressure, CO 2 plume, shear state on the fault, as well as shear displacement and stress path at selected fault monitoring points) are selected to assess CO 2 migration and reactivation of intersecting faults. The results show that induced ruptures are favored by low permeability faults due to high pressure buildup and by slip-weakening behavior resulting from fault strength reduction. The location of one injection well relative to fault orientation determines the magnitude of changes in effective normal stress and shear stress, affecting the location of induced ruptures. Well placement (two injection wells used in the paper) dominates pressure diffusion around the intersection and tips of faults. This redistributes changes in effective normal stress caused by each injection well, influencing the spatial distribution of ruptures along faults. A larger injection volume induces far-field ruptures that are controlled by stress transfer within the injection layer. The findings presented here can provide valuable insights into engineering operations for a long-term, safe, and reliable geologic CO 2 storage.

Fault permeability↗
DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Systematic Validation of Ensemble Cloud-Process Simulations using Polarimetric Radar Observations and Simulator over the NASA Wallops Flight Facility

The BiLateral Operational Storm-Scale Observation and Modeling (BLOSSOM) project was initiated to establish a long-term supersite to improve understanding of cloud physical states and processes as well as to support satellite and climate model programs over the WFF site via a bilateral approach of storm-scale observations and process modeling. This study highlights a noble systematic validation framework of the BLOSSOM ensemble cloud-process simulations through mixed-phase, light-rain, and deep-convective precipitation cases. The framework consists of creating a domain-shifted ensemble of large-scale forcing datasets, and configuring and performing cloud-process simulations with three different bulk microphysics schemes. Validation uses NASA S-band dual-POLarimetric radar (NPOL) observations in the form of statistical composites and skill scores via a polarimetric radar simulator and newly developed CfRad Data tool (CfRAD). While the simulations capture the overall structures of the reflectivity composites, polarimetric signals are still poorly simulated, mainly due to a lack of representation of ice microphysics diversity in shapes, orientation distributions, and their complex mixtures. Despite the limitation, this new ensemble-based validation framework demonstrates that 1) no particular forcing or microphysics scheme outperforms the rest and 2) the skill scores of coarse- and fine-resolution ensemble simulations with different domain-shifted forcing and microphysics schemes are highly correlated with each other with no clear improvement. On the other hand, this suggests that coarse-resolution ensemble simulations are relevant for selecting the best meteorological forcing and microphysics scheme before conducting computationally demanding large eddy simulations (LESs) in support of aircraft and satellite instrument development as well as cloud-precipitation-convection parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Infrared-Fused Vision-Based Thermoregulation Performance Estimation for Personal Thermal Comfort-Driven HVAC System Controls

Thermal comfort is one of the primary factors influencing occupant health, well-being, and productivity in buildings. Existing thermal comfort systems require occupants to frequently communicate their comfort vote via a survey which is impractical as a long-term solution. Here, we present a novel thermal infrared-fused computer vision sensing method to capture thermoregulation performance in a non-intrusive and non-invasive manner. In this method, we align thermal and visible images, detect facial segments (i.e., nose, eyes, face boundary), and accordingly read the temperatures from the appropriate coordinates in the thermal image. We focus on the human face since it is often clearly visible to cameras and is not merged into a hot background (unlike hands). We use a regularized Gaussian Mixture model to track the thermoregulation changes over time and apply a heuristic algorithm to extract hot and cold indices. We present a personalized and a generalized comfort modeling method, selected based on the availability of the occupant historical indices measurements in a neutral environment, and use the time-series of the hot and cold indices to define corrections to HVAC system operations in the form of setpoint constraints. To evaluate the efficacy of our proposed approach in responding to thermal stimuli, we designed a series of controlled experiments to simulate exposure to cold and hot environments. While applying personalized modeling showed an acceptable average accuracy of 91.3%, the generalized model’s average accuracy was only 65.2%. This shows the importance of having access to physiological records in modeling and assessing comfort. We also found that individual differences should be considered in selecting the cooling and heating rates when some knowledge of the occupant’s overall thermal preference is available.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

A continuum model for heat and mass transfer in moving-bed reactors for thermochemical energy storage

In this work, a continuum heat and mass transfer model coupling transport phenomena and high-temperature thermochemical reactions is developed for stationary packed-bed and counter-flow moving-bed reactors. After presenting the general modeling framework, we focus on the 2D axisymmetric version of the model for which validation is conducted with experimental results for a packed-bed reactor in the literature for manganese-iron oxide reduction/oxidation and an in-house counter-flow moving-bed reactor for magnesium-manganese oxide reduction up to 1450 °C. Transient simulation results including the local distributions of gas/solid temperatures, oxygen concentration and the extent of reaction, as well as the various energy flow components and energy conversion efficiencies are reported. The results based on the 2D axisymmetric model are also compared with those obtained from a previous 1D model. The comparison shows that capturing the radial variation is critical in reactor modeling and the 2D results demonstrate improved agreement with experiments. Specifically, large temperature variations along the radial direction are observed especially in the reaction zone; this non-uniform radial temperature distribution has a significant effect on the chemical reaction extent due to its strong dependence on temperature; and the overall oxygen concentration at the reactor exit and the predicted system efficiency are slightly lower in the 2D model compared to the 1D model. Finally, the present heat and mass transfer model can provide valuable insights into reactor design, scale-up, and operating conditions selection to maximize system energy storage efficiency.

25 ENERGY STORAGE↗
DOE OSTI2022

An antibody-free platform for multiplexed, sensitive quantification of protein biomarkers in complex biomatrices

Sensitive, multiplexed protein quantification remains challenging despite recent advancements in LC-MS assays for targeted protein biomarker quantification. High-sensitivity protein biomarker measurements usually require immuno-affinity enrichment of target protein; a process which is highly dependent on capture reagent and limited in capability to measure multiple analytes. Herein, we report a novel antibody-free platform, which measures multiple biomarkers from complex matrices employing a strategically optimized solid-phase extraction cleanup and orthogonal multidimensional LC-MS. Eight human protein biomarkers with different specifications were spiked into canine plasma as a model investigation system. The developed strategy achieved the desired sensitivity, robustness, and throughput via the following steps: (1) post digestion mixed-mode cation exchange-reverse phase SPE enrichment cleaned up the sample initially; (2) rapid, high-pH peptide fractionation further eliminated background components efficiently while selectively enriched signature peptides (SP) to provide sufficient sensitivity for multiple targets; and (3) trapping-micro-LC-MS analysis delivered high sensitivity comparable to a nano-LC-MS method but with much better robustness and throughput for the final analysis. Compared with a conventional LC-MS assay with direct protein digestion and limited clean-up, analysis with this antibody-free platform improved the LLOQ by 1–2 orders of magnitude for the eight protein biomarkers, reaching as low as 5 ng/mL in plasma, with feasible robustness and throughput. In conclusion, this platform was applied for the quantification of biomarkers of respiratory conditions in patients with various lung diseases, demonstrating real-world applicability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Secondary Crash Identification using Crowdsourced Waze User Reports

Secondary crashes are crashes that occur as a result of the nonrecurrent congestion originating from primary crashes, and always have a greater impact on safety and traffic than a single crash. A better understanding of secondary crashes would benefit traffic incident management, and this requires accurate identification of secondary crashes. This study explores using crowdsourced Waze user reports to identify secondary crashes. Here, a network-based clustering algorithm is proposed to extract the primary crash cluster, including all user reports originating from the primary crash, and any crash that occurred within the cluster would be a secondary crash. This method works as a filter to select accurate primary–secondary relationships, thus precisely identifying secondary crashes. A case study is performed with crashes occurring from June to December 2019 on a 30-mi stretch of I-40 in Knoxville, TN. A static threshold method (crash duration and 10 mi) was used to preselect the potential primary–secondary crash pairs, and 75 out of 708 crashes were identified as potential secondary crashes. Based on the preselected primary–secondary crash pairs, 17 secondary crashes were obtained with the proposed method and the results were compared with one of the commonly used methods, the speed contour plot method. Though the proposed method captured fewer secondary crashes, it did identify several secondary crashes that could not be observed with the speed contour plot method. The results showed the applicability of the method and the potential of crowdsourced Waze user reports in secondary crash identification.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Synthesis of Nitro Compounds from Nitrogen Dioxide Captured in a Metal-Organic Framework

Increasing levels of air pollution are driving the need for the development of new processes that take “waste-to-chemicals”. Herein, we report the capture and conversion under ambient conditions of a major air pollutant, NO 2 , using a robust metal-organic framework (MOF) material, Zr-bptc (H 4 bptc = 3,3',5,5'-biphenyltetracarboxylic acid), comprising {Zr 6 (μ 3 -O) 4 (μ 3 -OH) 4 (COO) 12 } clusters linked by 4-connected bptc 4– ligands in an ftw topology. At 298 K, Zr-bptc shows exceptional stability and adsorption of NO 2 at both low (4.9 mmol g –1 at 10 mbar) and high pressures (13.8 mmol g –1 at 1.0 bar), as measured by isotherm experiments. Dynamic breakthrough experiments have confirmed the selective retention of NO 2 by Zr-bptc at low concentrations under both dry and wet conditions. The immobilized NO 2 can be readily transformed into valuable nitro compounds relevant to construction, agrochemical, and pharmaceutical industries. In situ crystallographic and spectroscopic studies reveal strong binding interactions of NO 2 to the {Zr 6 (μ 3 -O) 4 (μ 3 -OH) 4 (COO) 12 } cluster node. This study paves a circular pathway to enable the integration of nitrogen-based air pollutants into the production of fine chemicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Cooperative and inhibitory ion transport in functionalized angstrom-scale two-dimensional channels

Significant success has been achieved in fabricating angstrom-scale artificial solid ionic channels aiming to replicate the biological ion channels (BICs). Besides high selectivity, BICs also exhibit sophisticated ion gating and interplay. However, such behavior and functionality are seldomly recreated in the artificial counterparts due to the insufficient understanding of the molecular origin. Here we report cooperative and inhibitory ion transport in angstrom-scale acetate functionalized MoS 2 two-dimensional channels. For cooperative ion transport, the permeability of K + is doubled in the presence of only 1% Pb 2+ (versus K + by molarity), while the permeability of Pb 2+ is independent of K + . Molecular dynamics simulations reveal complex interplay among K + , Pb 2+ , and the anions in governing the cooperativity, such that Pb 2+ ions capture and slow down the anions via long-range interaction, which leads to the synchronization of anions with K + to transport as ion pairs with reduced interaction with the channel surface. For inhibitory ion transport, divalent Co 2+ (or Ba 2+ ) and Pb 2+ can replace each other in the confined channel and compete for the limited transport cross section. Our work reveals ion transport phenomena in extreme confinement and highlights the potential of manipulating ion interplay in confinement for achieving advanced functionalities.

Functionalized Angstrom↗
DOE OSTI2023

Rational design of stable functional metal–organic frameworks

Functional porous metal–organic frameworks (MOFs) have been explored for a number of potential applications in catalysis, chemical sensing, water capture, gas storage, and separation. MOFs are among the most promising candidates to address challenges facing our society related to energy and environment, but the successful implementation of functional porous MOF materials are contingent on their stability; therefore, the rational design of stable MOFs plays an important role towards the development of functional porous MOFs. In this Focus article, we summarize progress in the rational design and synthesis of stable MOFs with controllable pores and functionalities. The implementation of reticular chemistry allows for the rational top-down design of stable porous MOFs with targeted topological networks and pore structures from the pre-selected building blocks. We highlight the reticular synthesis and applications of stable MOFs: (1) MOFs based on high valent metal ions (e.g., Al 3+ , Cr 3+ , Fe 3+ , Ti 4+ and Zr 4+ ) and carboxylate ligands; (2) MOFs based on low valent metal ions (e.g., Ni 2+ , Cu 2+ , and Zn 2+ ) and azolate linkers. We envision that the synthetic strategies, including modulated synthesis and post-synthetic modification, can potentially be extended to other more complex systems like metal-phosphonate framework materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Reactive capture and electrochemical conversion of CO 2 with ionic liquids and deep eutectic solvents

Ionic liquids (ILs) and deep eutectic solvents (DESs) have tremendous potential for reactive capture and conversion (RCC) of CO 2 due to their wide electrochemical stability window, low volatility, and high CO 2 solubility. There is environmental and economic interest in the direct utilization of the captured CO 2 using electrified and modular processes that forgo the thermal- or pressure-swing regeneration steps to concentrate CO 2 , eliminating the need to compress, transport, or store the gas. The conventional electrochemical conversion of CO 2 with aqueous electrolytes presents limited CO 2 solubility and high energy requirement to achieve industrially relevant products. Additionally, aqueous systems have competitive hydrogen evolution. In the past decade, there has been significant progress toward the design of ILs and DESs, and their composites to separate CO 2 from dilute streams. In parallel, but not necessarily in synergy, there have been studies focused on a few select ILs and DESs for electrochemical reduction of CO 2 , often diluting them with aqueous or non-aqueous solvents. The resulting electrode–electrolyte interfaces present a complex speciation for RCC. In this review, we describe how the ILs and DESs are tuned for RCC and specifically address the CO 2 chemisorption and electroreduction mechanisms. Critical bulk and interfacial properties of ILs and DESs are discussed in the context of RCC, and the potential of these electrolytes are presented through a techno-economic evaluation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Aqueous Carbon Capture Using Guanidinium-Functionalized Hollow Fiber Sorbent Contactors

As part of the growing suite of technologies aimed at combatting rising temperatures, negative emissions technologies have become a powerful tool in the global effort to minimize the consequences of human-induced climate change. Among these, carbon removal from aqueous sources, which contain much higher carbon concentrations than the atmosphere, remains largely unexplored. Indeed, developing robust and efficient carbon capture materials for usage in complex aqueous environments remains a significant challenge. Here, we explore the potential of functionalizing polyvinylidene fluoride (PVDF) hollow fiber contactors grafted with a guanidinium-derived polymer sorbent for carbon removal from aqueous sources, including saline waters. Computational screening against amine-based analogs is utilized to identify guanidinium as a promising motif for bicarbonate (HCO 3 – ) ions binding. To leverage this finding, synthesis of a guanidinium polymer and subsequent covalent grafting onto PVDF hollow fibers is employed to structured polymer–sorbent–grafted hollow fiber contactors. Our prototype achieves an initial HCO 3 – removal of 34% with an increase to 98% after four cycles. The functionalized fibers demonstrate aqueous stability over 13 adsorption/desorption cycles in model NaHCO 3 solutions where regeneration is facilitated by a mild pH swing. Importantly, the system maintains selective performance in the presence of competitive chloride ions over multiple cycles; carbon removal remained above 10% even at high (10:1) NaCl/NaHCO 3 ratios. These findings demonstrate the feasibility of sorbent-based aqueous carbon removal and highlight its potential as a promising approach for negative emissions.

carbon capture↗
DOE OSTI2020

Quantifying Progress Across Different CMIP Phases With the ESMValTool

Abstract More than 40 model groups worldwide are participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6), providing a new and rich source of information to better understand past, present, and future climate change. Here, we use the Earth System Model Evaluation Tool (ESMValTool) to assess the performance of the CMIP6 ensemble compared to the previous generations CMIP3 and CMIP5. While CMIP5 models did not capture the observed pause in the increase in global mean surface temperature between 1998 and 2013, the historical CMIP6 simulations agree well with the observed recent temperature increase, but some models have difficulties in reproducing the observed global mean surface temperature record of the second half of the twentieth century. While systematic biases in annual mean surface temperature and precipitation remain in the CMIP6 multimodel mean, individual models and high‐resolution versions of the models show significant reductions in many long‐standing biases. Some improvements are also found in the vertical temperature, water vapor, and zonal wind speed distributions, and root‐mean‐square errors for selected fields are generally smaller with reduced intermodel spread and higher average skill in the correlation patterns relative to observations. An emerging property of the CMIP6 ensemble is a higher effective climate sensitivity with an increased range between 2.3 and 5.6 K. A possible reason for this increase in some models is improvements in cloud representation resulting in stronger shortwave cloud feedbacks than in their predecessor versions.

Bock, L.↗
DOE OSTI2020

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↗
DOE OSTI2020

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

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

Crank angle-resolved PIV↗
DOE OSTI2022

Benchmark Exercise for the Control Rod Swelling Evaluation

The VTR core has six reactivity control assemblies and three safety assemblies. The control assemblies or primary control rods are adjusted during the normal operation to balance the core reactivity and to control the reactor power. A typical control assembly radial layout is presented in Figure 1. The figure shows the swelled absorber (B 4 C) rod. Initially, helium gas fills the gap between the pin and the cladding before irradiation swelling takes place. For VTR, HT9 steel was selected as the cladding and duct material. The main neutron absorbing material used in the VTR is B 4 C. When residing in the core, the neutronics, thermophysical, and mechanical properties of the materials used in a control assembly will degrade due to accumulated neutron damage. Material degradation limits how long a control assembly can reside in the core. Many phenomena affect the control assembly lifetime, such as the loss of reactivity worth due to B 4 C depletion, the mechanical interaction of the absorber rod and the cladding due to B 4 C swelling, the helium gas buildup in the pin due to B-10 capture, etc. B 4 C swelling, which causes closure of the gap between the absorber rod and the cladding, is usually considered as the main limiting factor from past experience. An initial study was conducted at PNNL to evaluate the irradiation behavior of a VTR control assembly. The evaluation was performed using the CNRD2 code that was initially developed for the FFTF. The study also included an assessment of the VTR control assembly and focused on a 61-pin control assembly design, which is different from that used (37-pin design) in the core design study. The study conducted by PNNL was reviewed independently by ANL. A Python script referred to as the Control Assembly Evaluation Script (CAES) was developed for the independent review and additional assessment of 37-pin control assembly design. The script has focused on the assessment of the absorber rod swelling for its importance in determining the control assembly lifetime. CAES uses geometry, neutronics, materials data as input to predict the swelling of the absorber rod during its residence in the reactor core. The results from CAES showed some non-negligible differences against the PNNL results. Some of the differences can be attributed to the different interpretation of the control rod assembly dimensions. To resolve this issue, a benchmark exercise was proposed. The benchmark specification was developed by PNNL. The benchmark exercise was performed independently at PNNL and ANL using different codes/scripts (CRND2 and CAES). This memo documents the results calculated using the different codes. However, this report is limited to presenting the results obtained. Further investigation of the cause of the observed difference will be performed as part of future activities, pending continuation of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗