Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Plumes”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Western sector in-situ chemical oxidation project: Supplemental results after injection activities

The purpose of the Western Sector ISCO project is to inject chemical oxidants (i.e., potassium permanganate and sodium persulfate) into the eight injections wells (i.e., WSI001B, WSI001C, WSI002B, WSI002C, WSI003B, WSI003C, WSI004B, and WSI004C) and evaluate the effectiveness of the chemical oxidants at degrading indiscriminate DNAPL that is present in the subsurface. Fifteen monitoring wells (i.e., WSM001BB, WSM001B, WSM001CC, WSM001C, WSM002BB, WSM002B, WSM002CC, WSM002C, WSM003BB, WSM003B, WSM003CC, WSM003C, MSB107B, MSB107CC, and MSB107C) are currently being used by SRS to monitor the effect of the oxidant injections on the VOC plume of the Western Sector of the M-Area HWMF. Figure 1 shows the project location in reference to the M-Area Settling Basin. This is the second of three planned interim reports to track the performance of the project. The purpose of this document is to report groundwater monitoring data collected during the baseline monitoring before oxidant injections occurred and the fifteen months after oxidant injections were completed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Selection of Vadose Zone Flow and Transport Properties with Gravel Fraction Corrections for the Hanford Site Composite Analysis and Cumulative Impact Evaluation

This environmental calculation file (ECF) is a compilation of hydrologic and radiological/chemical properties to be used for the updated Hanford Site Composite Analysis (CA) and Cumulative Impact Evaluation (CIE) vadose zone (VZ) modeling. The CA will provide an all-pathways dose projection to a hypothetical future member of the public from all planned low-level radioactive waste disposal facilities and potential contributions from all other projected end-state sources of radioactive material left at the Hanford Site following site closure. Its primary purpose is to support the decision-making process of the U.S. Department of Energy (DOE) under DOE O 435.1 Chg 1, Radioactive Waste Management, related to managing low-level waste disposal facilities at the Hanford Site. The CIE evaluates the effects of cleanup decisions regarding groundwater quality in the Hanford Site Central Plateau (DOE/RL-2018-69, Cumulative Impact Evaluation Technical Approach Document). Due to the complexity and large number of waste sites in source operable units (OUs), the computational tools used for the CIE must be capable of representing a range of site conditions and source terms in the VZ while also efficiently computing the impact that cleanup decisions have on the underlying aquifer. Waste-site proximity between and within source OUs has resulted in contaminants commingling in the vadose and saturated zones in complex ways. Plume commingling requires cleanup decisions to be evaluated considering the surrounding waste sites and existing groundwater contamination, therefore demonstrating the need to evaluate cumulative impacts from the VZ to groundwater.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Probabilistic predictions and uncertainty estimation for radiological and nuclear effects modeling

Better characterization of meteorological uncertainty is needed to predict the effects of radiological and nuclear releases in the atmosphere. Multi-physics ensembles of atmospheric simulations can estimate uncertainty of physical processes, but they are often limited for computational reasons. Because weather uncertainty is important for emergency response applications, but the conventional way of running multi-physics ensembles is infeasible, we have developed new machine learning-based methods that quickly estimate and identify key sources of meteorological uncertainty. Machine learning accelerates the analysis of weather uncertainty by iteratively learning about model physics options that affect plume predictions. We have demonstrated the power of machine learning by creating a large multi-physics weather and dispersion ensemble dataset containing 1,200 simulations for two hypothetical radiological release scenarios. By training on less than 10% of the ensemble, we can quickly and accurately predict different outcomes, including the mass of radiological material deposited and spatial distribution and extent of contaminated areas.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Hotspot 2.07.1 doses (FGR-13) from releases of 1 g Pu-239 (SBK Calcs)

This calculation is to be used for accident analysis at B239, B331, B332, and B334. The objective of this analysis is to determine: The dose per gram of Pu-239 Solubility Class S (FGR-13) for various release heights, receptor heights and distances from 100 to 1000 meters. This calculation can be used to evaluate particulate plume dispersion releases for the co-located worker and public receptors.

61 RADIATION PROTECTION AND DOSIMETRY↗
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

E-Area Low-Level Waste Facility GoldSim System Model

This report documents the development of the E-Area Low-Level Waste Facility (ELLWF) trench system model. The GoldSim® Monte Carlo simulation software is utilized to model the release and transport of radiological inventory disposed (both currently and in the future) within Engineered and Slit Trenches. This model is in support of the sensitivity and uncertainty analysis for the ELLWF Performance Assessment. The ELLWF system model utilizes a hybrid-approach to accurately describe the disposal system. The Hydrologic Evaluation of Landfill Performance model provides the infiltration data to both PORFLOW and GoldSim. PORFLOW is used to calibrate the GoldSim model to ensure confidence in the stochastic results. Finally, the concentrations from GoldSim transport simulations are fed into the SRNL Dose Toolkit to calculate dose impacts and assess plume interaction

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Sensor Recommendations for Long Term Monitoring of the F-Area Seepage Basins

In mid-2018, a new paradigm for long-term monitoring was developed after of decade of applied research projects funded by the Department of Energy’s office of Environmental Management Technology Development program. The program at SRNL was focused on transitioning complex environmental waste sites from active to passive remediations strategies. A key result of these studies was that the use of enhanced attenuation approaches at radiologically contaminated sites will result in the creation of secondary source areas in the subsurface that will require monitoring for decades. Alternative monitoring approaches are being developed and tested at the Savannah River Site’s F-Area Hazardous Waste Management Facility, the new paradigm provides innovative solutions that will significantly lower costs of monitoring through the coupling of data collection, machine learning and deterministic groundwater modeling. The foundation of this approach is a well-optimized network of sensors for measuring hydrogeochemical master variables that control, and therefore act as indicators of groundwater contaminant transport. By monitoring changes in the controlling master variables over time arising from geological and environmental shifts, predictive modelling can assist with identifying new strategies for ensuring regulatory requirements are met if trends toward conditions for potential remobilization of attenuated contaminants are detected. In this report, we evaluated commercially available single parameter sensor platforms (e.g., temperature/depth) and configurable multi-parameter sensor platforms (e.g., pH, oxidation-reduction potential, temperature, depth, dissolved oxygen, and conductivity). Each was scored using an optimization function based on how well the system supports the proposed long-term monitoring paradigm, in general, and the site-specific conditions at F-Area, in particular. Several viable sensor systems were identified. Of these, a combined platform including the In-Situ Aqua TROLL 500 multi-parameter sensor platform and the In-Situ temperature/depth sensor had the highest rating and was identified as the most suitable candidate for installation and monitoring of the master variables and potentiometric surface that control groundwater contaminant plumes emanating from the F-Area Seepage Basins. The discussion of recommended potential deployment locations builds upon recommendations made by Denham et al (2019).

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Laser Plasma Interaction in the MG Magnetic Field (Final Report)

Strong magnetic fields change properties of plasma. Employment of the magnetic fields can improve conditions for inertial confinement fusion. Compression and heating of the magnetized plasma is a base of the MagLIF pulsed power approach to the controlled fusion. Study of expansion and heating of plasma in the external magnetic fields, development of plasma instabilities, parametric and other effects are important for basic physics and applications. The main goal of the research for the grant was experimental investigation of laser produced plasma in the well characterized and controlled MG magnetic fields. We developed a platform for investigation of high intensity laser-plasma interaction (LPI) in the MG fields using a university-scale pulsed power machine. A Zebra machine at UNR generated transverse magnetic fields of 3 MG on the surface of the rod load and longitudinal fields of 1.4 MG in the coil loads. The magnetic fields were characterized with two-color Faraday rotation diagnostics. The pulsed power machine provided almost static magnetic fields for LPI compared to the much faster laser driven magnetic fields. We found that eddy currents in the rising magnetic field generate plasma on the metal targets. Eddy currents did not impact Si and CH dielectric targets. A Leopard laser with intensity in the focus up to 10 19 W/cm 2 produced plasma in the MG magnetic fields. The modified reflected Faraday diagnostics was used to study laser driven B-field at the MTW laser at Laboratory for Laser Energetics, UR. We demonstrated a slow dependence of the magnetic field in the coil on the applied laser intensity. The magnetic field increased by a factor of 2 if laser intensity increased by a factor of 30. This scaling is important for comparison of experiments at different laser facilities. Expansion of the plasma in the 2-3 MG external B-field was studied. Laser-produced plasma in the azimuthal magnetic field took the unique form of a thin disc expanding radially with a velocity of 250 km/s and confined in the vertical direction. A HYDRA MHD program at the University of Rochester was applied for simulations of LPI in the magnetic field B >3 MG. Simulations for the plasma disc parameters and the expansion dynamics were in a good agreement with experiments. Generation of narrow dense plasma jets in the longitudinal magnetic field of 0.6-0.8 MG was studied. Narrow plasma jets reached a length of 3-4 mm with the electron plasma density of (0.2-1.2)x10 20 cm -3 . A jet tip propagated with the velocity of 160-200 km/s. MHD simulations showed a good agreement of the dynamics of the formation of plasma jets with experiments. These jets are relevant to astrophysical jets. Simulations of K-shell spectra of Si plasma was performed for LPI in the B-field. The laser prepulse generated plasma near the laser target. A PrismSPECT model with a MeV electron beam produced by the laser was used. The spectral modelling showed the increased plasma density by a factor by 2-3 in the magnetic field. This was in agreement with the observed dynamics of the plasma plume confined by the axial magnetic field. The two-plasmon decay was studied in the 2.5-3 MG transversal magnetic fields. Strong 2-3.5 nm spectral widening and a 2-4 nm shifts of “red” and “blue” 3/2ω 0 spectral components were observed. The large shift and widening exceed the expected temperature and magnetic shift. PIC simulations are performing to clarify physics of this effect. The experimental research program for the grant is completed. New plasma effects in the MG fields were studied at the university-scale pulsed power generator. Two graduate and two undergraduate UNR and UR students carried experiments at the Zebra generator, supported and developed plasma diagnostics, provided data processing, and performed MHD simulations of LPI in the strong magnetic field. Research results are published in 9 referred papers and presented in 12 conferences. Program manager: Kramer U. Akli; Collaborators: R. C. Mancini and H. Sawada, University of Nevada, Reno (UNR); R. Betti and A. V. Maximov, University of Rochester (UR).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Source Term of Accidents Involving Vehicle Crashes and Fuel Fires

This calculation is to be used when developing accident analysis scenarios involving fuel fires at the WSF. The objectives of the calculation are to determine (1) the Source Term from releases of TRU waste involved in impacts, impacts with fuel fire, and fuel fire with no impact; and (2) the plume sensible heat from various sizes of gasoline pool fires.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Source Term of Accidents Involving Aircraft Crashes and Fuel Fires

This calculation is to be used when developing accident analysis scenarios involving aircraft fuel fires at the Waste Storage Facility (WSF). The objectives of the calculation are to determine (1) the Source Term from releases of TRU waste involved in impacts, impacts with fuel fire, and fuel fire with no impact; and (2) the plume sensible heat from various sizes of gasoline pool fires.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

AniMACCS User Guide

AniMACCS is a utility code in the MELCOR Accident Consequence Code System (MACCS) software suite that allows for certain MACCS output information to be visually displayed and overlaid onto a geospatial map background. AniMACCS was developed by Sandia National Laboratories for the U.S. Nuclear Regulatory Commission. MACCS is designed to calculate health and economic consequences following a release of radioactive material in the atmosphere. MACCS accomplishes this by modeling the atmospheric dispersion, deposition, and consequences of the release, which depend on several factors including the source term, weather, population, economic, and land-use characteristics of the impacted geographical area. From these inputs, MACCS determines the characteristics of the plume, as well as ground and air concentrations as a function of time and radionuclide.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2020

Determining the Appropriate Emergency Planning Attributes for Microreactors

Micro-reactor designs incorporate design features that provide very low reactor decay heat power at 24 hours after shutdown that is manageable in comparison with larger reactors. This attribute translates into a low probability of core damage negligible offsite dose. This paper provides the conceptual framework for appropriately structuring emergency planning requirements for micro-reactors while ensuring the U.S. Nuclear Regulatory Commission’s (NRC) commitment to safety is met. Due to the diversity in reactor designs comprising the micro-reactor community, this paper identifies the necessary concepts that should be considered to ensure emergency plans for micro-reactors are appropriate for the risk. This paper does not mandate specific design features for micro-reactors but provides reactor designers a conceptual framework for a scalable, graded approach to Emergency Planning Standards that should be considered in their individual designs in order to support simplified emergency plans. Given that accident source terms associated with micro-reactors are essentially negligible when compared with those for large light water reactors, revisions to emergency planning standards are justified. A graded approach to implementing emergency planning guidance can be used to appropriately structure micro-reactor emergency plans and reduce the size of the emergency planning zone and plume exposure pathway. Appropriately structuring emergency planning requirements will better optimize licensee and offsite agencies emergency planning resources and reduce the resources associated with emergency planning.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

Distribution of Infiltration in the 216-U-10 and 216-B-3 Pond Systems 1944-1997

A modeling system to simulate fate and transport of contaminant plumes in groundwater beneath the Central Plateau to support the Cumulative Impact Evaluation (CIE) and Composite Analysis (CA) is being developed. The CIE will evaluate the cumulative effect of multiple Central Plateau radiological and chemical sources on groundwater quality (DOE/RL-2018-69, Cumulative Impact Evaluation Technical Approach Document ). This will enable cleanup decisions to be evaluated in the context of all Central Plateau waste sites and existing groundwater conditions. The CA evaluates the cumulative effect of all sources of radioactive materials on radiological dose to future members of the public (DOE/RL-2018-60, Annual Status Report (FY 2018): Composite Analysis of Low Level Waste Disposal in the Central Plateau of the Hanford Site ). Fate and transport modeling will be used to forecast future groundwater concentrations for both the CIE and CA. The CIE/CA subsurface modeling system is divided into source (inventory), vadose zone, and saturated zone facets. Source inventory is a necessary input to the vadose zone facet, and releases from the vadose zone to groundwater are inputs to the saturated zone facet. The vadose zone facet simulates historical and future releases from the vadose zone to groundwater for both liquid waste and solid waste disposal sites. Effluent release volumes and radionuclide inventory estimates for liquid waste sites are contained in ECFHANFORD- 17-0079, Hanford Soil Inventory Model (SIM-v2) Calculated Radionuclide Inventory of Direct Liquid Discharges to Soil in the Hanford Site’s 200 Areas , herein referred to as SIM-v2. The liquid waste sites in SIM-v2 include the 216-B-3 Pond near 200 East Area and the 216-U-10 Pond in 200 West Area. The 216-B-3 Pond system consisted of a main pond and 3 expansion lobes that were operational during overlapping time periods. It also included influent ditches, but these are not listed as liquid waste sources in SIM-v2 and are not included in the CIE/CA modeling. The total effluent volume released to the pond system was measured, but the partitioning of this volume to the various expansion lobes was not determined and is not accounted for in SIM-v2. Regarding the 216-U-10 Pond, it received water from several influent ditches. The pond and the influent ditches are listed as liquid waste sites in SIM-v2 and are included as sources in the CIE/CA modeling. Because water in the ditches flowed into the pond, only a portion of the effluent volumes and contaminant inventories assigned to the ditches in SIM-v2 infiltrated from the ditches. Thus, partitioning of infiltration and contaminant inventories between the ditches and the pond is not fully accounted for in SIM-v2. Further, some of the water in the 216-U-10 Pond entered overflow ditches and the amount of this overflow was not recorded and is not accounted for in SIM-v2. Thus, the purpose of this Environmental Calculation File (ECF) is to partition the infiltration of wastewater and contaminant inventory between the main pond and each expansion lobe of the 216-B-3 Pond system, and to partition the infiltration of wastewater and contaminant inventory between the ditches and pond of the 216-U-10 Pond system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Charged Wellbore Casing Controlled Source Electromagnetics (CWC-CSEM) for Reservoir Imaging and Monitoring (Final report)

This project addresses the needs of the U.S. Department of Energy (DOE) to develop advanced monitoring technologies and protocols to track the fate of subsurface carbon dioxide (CO2) plumes for carbon storage. Specifically, the project seeks to develop and test a unique and novel system of technologies consisting of electromagnetic data acquisition, coupled multiphysics imaging, and reservoir model enhancement to understand the migration and long-term distribution of CO2 in the subsurface. The overarching objective is to develop an integrated approach for long term monitoring of carbon storage. The two main components of the project include the methodology development and the test of the method at a field site. The methodology component consists of 1) developing the field procedure and protocol for collecting time-lapse controlled-source electromagnetic (CSEM) data with source electric current injected into the subsurface through wellbore casings; 2) building of background 3D electrical conductivity utilizing multiple sources of data such as supplemental surface transient EM (TEM) surveys, well-logs, and seismic structural information, for enhancing CSEM signal from reservoir depths; and 3) coupled multiphysics simulations and inversion of CSEM data constrained by production data and by structural information from seismic imaging of the reservoir and overlying formations. The testing component used the field site of Bell Creek Oil Field, which served both as a field laboratory for the method development as well as a test site to evaluate the CSEM signal strengths and the methodology developed in this research project. We have accomplished all the proposed tasks and developed the methodology as planned. These include the procedure for time-lapse CSEM data acquisition, data processing techniques, integration with 3D conductivity model building, fast reservoir simulation for history matching using machine learning, and interpreting CSEM data with coupling to the reservoir modeling. Collectively, the outcome of these tasks form a coherent workflow that can be applied to monitor dedicated carbon storage in saline reservoirs. The testing component evaluated the applicability and limitations of the method, and concluded that the method would be ideal for monitoring dedicated carbon storage sites utilizing saline reservoirs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluation of Dose- and Risk-Based Groundwater Cleanup Levels for Low Energy Beta Radioisotopes

Low-energy beta-emitting radionuclides that were released historically during reactor operations and plutonium separations activities at some U.S. Department of Energy sites have migrated into the groundwater, forming contaminant plumes that are subject to federally regulated remediation actions. At the Hanford Site, the low-energy beta-emitting radionuclides include iodine-129, technetium-99, chlorine-36, carbon-14, and tritium (H-3). All are highly mobile in the subsurface, and except for tritium, and have very long half-lives—thousands to millions of years. The geochemistry and transport behavior of these contaminants in the subsurface present significant challenges for remediation of groundwater to federal drinking water standards (DWS)—the appropriate or relevant and applicable requirements (ARARs) for cleanup. For some of the low-energy beta-emitter contaminants, particularly iodine-129, cleanup and restoration of groundwater to DWS may not be attainable within a reasonable timeframe using currently available treatment technologies.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗
DOE OSTI2020

Dynamic Behavior of Natural Seep Vents: Analysis of Field and Laboratory Observations and Modeling (Final Scientific/Technical Report)

In this project, we have analyzed data collected by the U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) in a high pressure water tunnel (HPWT) and data from two research cruises to natural seeps in the Gulf of Mexico to adapt and validate a numerical model to predict the dynamics of natural seeps in the deep oceans. The HPWT data include video observations of the shrinkage rate of individual methane and natural gas bubbles under simulated deep-water conditions. Field data were collected during two cruises by the Gulf Integrated Spill Research (GISR) Consortium led by Texas A&M University and funded by the Gulf of Mexico Research Initiative (GoMRI). These data included in situ observations from a remotely operated vehicle (ROV) of gas bubbles at two natural seep sites in the Gulf and acoustic observations of the natural seep bubble flares in the ocean water column. The acoustic data were from multibeam echosounders, one mounted in a forward-looking orientation on the ROV and another mounted down-looking in the haul of the ship. All of these laboratory and field data were focused on the dynamics of natural gas bubbles at temperatures and pressures favorable for clathrate hydrate formation between the gas and water. Our analyses of this data focused on understanding the mechanisms responsible for gas bubble dissolution within the hydrate stability zone (HSZ) of the oceans. Ice-like hydrate shells may form on the bubble-water interface under these conditions, and it was unknown how this might affect the mass transfer of gas into the ocean. We were able to extract bubble shrinkage rates from the HPWT datasets. Using this data we determined that mass transfer coefficients with and without a hydrate shell match empirical values for bubbles in contaminated systems (so-called dirty bubbles contaminated by naturally occurring surfactants). We also showed that free gas, and not gas hydrate, is the dominant dissolving phase when the hydrate sub-cooling is below 11 degree Celsius (temperature difference between hydrate the hydrate formation temperature and ambient temperator) or the pressure is reducing as bubbles rise through the ocean water column. Using this mass transfer model, our numerical model of bubble dissolution matched the over 200 HPWT experiments with an average error of 10% for predicting the bubble size at the end of an experiment. From field data in the literature, we also observed that gas bubbles dissolve faster when they are initially released, following mass transfer coefficients for so-called clean-bubbles (those not yet contaminated by surfactants). Shortly after release within the HSZ, a hydrate shell forms on the bubble-water interface, and the mass transfer reduces to rates matching those of dirty bubbles. We correlated this transition time from clean to dirty bubble behavior with the initial bubble surface area and the hydrate sub-cooling. With this model for hydrate formation time and using the mass transfer coefficients deduced from the HPWT data, we validated our numerical model for predicting the rise heights of natural seep flares in the oceans. Flare heights are commonly observed in haul-mounted acoustic multibeam data. The numerical model predicts bubbles to rise high in the ocean water column owing to the slower mass transfer rates for dirty bubbles that accompany the majority of their rise time. We found that the numerical model predictions matched the observed flare heights within 5% to 10% accuracy when we compared the rise heights of the largest bubbles released from the seafloor with the bubbles acoustically visible in the multibeam data. Bubbles become acoustically transparent as they shrink to sizes of order 1 mm in diameter for the multibeam frequencies used in the field. The forward-looking multibeam on the ROV also provided data on the lateral spreading of bubbles in natural seep flares. Our analysis of this data showed that spreading follows a diffusion process, with the effective diffusivity correlating with the wobbling length scale of these ellipsoidal bubbles. When we apply this diffusivity in a random displacement model of bubble spreading, our numerical simulations match closely the lateral spread observed by the M3 in the ocean water column. Finally, we compared the seep model predictions for the acoustic properties of these natural seep plumes with that observed by the acoustic instruments in the field. The M3 and EM 302 observations were converted to relative values of target strength using a calibration we obtained in the laboratory for the M3 and using an algorithm from the manufacturer for the EM 302. Comparing the numerical seep model to these data, we obtain good agreement over the whole height of rise of these bubble flares. This further validates the numerical model. Overall, our validated seep model captures the key dynamics of gas bubbles released from natural seeps in the oceans and helps to predict the fate of methane in the water column.

03 NATURAL GAS↗
DOE OSTI2020

Emergency Planning Considerations for Advanced Nuclear Power Reactors

The purpose of this scoping study is to develop an approach for establishing emergency planning requirements for advanced nuclear power reactors and other new reactor technologies. The approach considers existing emergency planning requirements and guidance. More specifically the study focuses on establishing criteria and process to determine the size of the plume and ingestion exposure pathway emergency planning zone. The review of emergency planning in place for existing licensed nuclear facilities provides insight and informs the suggested process for establishing this Emergency Planning Zone (EPZ) process.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗