Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiphase”

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 1,567 records · Page 87
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Fault-related CO 2 degassing, geothermics, and fluid flow in southern California basins---Physiochemical evidence and modeling

Our studies have had an important impact on societal issues. Experimental and field observations show that CO 2 degassing, such as might occur from stored CO 2 reservoir gas, can result in significant stable isotopic disequilibrium. In the offshore South Ellwood field of the Santa Barbara channel, we show how oil production has reduced natural seep rates in the area, thereby reducing greenhouse gases. Permeability is calculated to be ~20-30 millidarcys for km-scale fault-focused fluid flow, using changes in natural gas seepage rates from well production, and poroelastic changes in formation pore-water pressure. In the Los Angeles (LA) basin, our characterization of formation water chemistry, including stable isotopic studies, allows the distinction between deep and shallow formations waters. Our multiphase computational-based modeling of petroleum migration demonstrates the important role of major faults on geological-scale fluid migration in the LA basin, and show how petroleum was dammed up against the Newport-Inglewood fault zone in a “geologically fast” interval of time (less than 0.5 million years). Furthermore, these fluid studies also will allow evaluation of potential cross-formational mixing of formation fluids. Lastly, our new study of helium isotopes in the LA basin shows a significant leakage of mantle helium along the Newport Inglewood fault zone (NIFZ), at flow rates up to 2 cm/yr. Crustal-scale fault permeability (~60 microdarcys) and advective versus conductive heat transport rates have been estimated using the observed helium isotopic data. The NIFZ is an important deep-seated fault that may crosscut a proposed basin decollement fault in this heavily populated area, and appears to allow seepage of helium from the mantle sources about 30 km beneath Los Angeles. The helium study has been widely cited in recent weeks by the news media, both in radio and on numerous web sites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Advanced Envelope Research for Factory Built Housing, Phase 3. Whole-House Prototyping

The Advanced Envelope Research effort will provide factory homebuilders with high performance, cost-effective envelope designs that can be effectively integrated into the plant production process while meeting the thermal requirements of the 2012 IECC standards. This work is part of a multiphase effort. Phase 1 identified seven envelope technologies and provided a preliminary assessment of three methods for building high performance walls. Phase 2 focused on developing viable product designs, manufacturing strategies, addressing code and structural issues, and cost analysis of the three selected options. An industry advisory committee helped narrow the research focus to perfecting a stud wall design with exterior continuous insulation (CI). This report describes Phase 3, which was completed in two stages and continued the design development effort, exploring and evaluating a range or methods for applying CI to factory built homes. The scope also included material selection, manufacturing and cost analysis, and prototyping and testing. During this phase, a home was built with CI, evaluated, and placed in service. The experience of building a mock up wall section with CI and then constructing on line a prototype home resolved important concerns about how to integrate the material into the production process. First steps were taken toward finding least expensive approaches for incorporating CI in standard factory building practices and a preliminary assessment suggested that even at this early stage the technology is attractive when viewed from a life cycle cost perspective.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2015

Uniaxial Creep as a Control on Mercury Intrusion Capillary Pressure in Consolidating Rock Salt

The nature of geologic disposal of nuclear waste in salt formations requires validated and verified two-phase flow models of transport of brine and gas through intact, damaged, and consolidating crushed salt. Such models exist in other realms of subsurface engineering for other lithologic classes (oil and gas, carbon sequestration etc. for clastics and carbonates) but have never been experimentally validated and parameterized for salt repository scenarios or performance assessment. Models for waste release scenarios in salt back-fill require phenomenological expressions for capillary pressure and relative permeability that are expected to change with degree of consolidation, and require experimental measurement to parameterize and validate. This report describes a preliminary assessment of the influence of consolidation (i.e. volume strain or porosity) on capillary entry pressure in two phase systems using mercury injection capillary pressure (MICP). This is to both determine the potential usefulness of the mercury intrusion porosimetry method, but also to enable a better experimental design for these tests. Salt consolidation experiments are performed using novel titanium oedometers, or uniaxial compression cells often used in soil mechanics, using sieved run-of-mine salt from the Waste Isolation Pilot Plant (WIPP) as starting material. Twelve tests are performed with various starting amounts of brine pore saturation, with axial stresses up to 6.2 MPa (~900 psi) and temperatures to 90°C. This corresponds to UFD Work Package 15SN08180211 milestone “FY:15 Transport Properties of Run-of-Mine Salt Backfill – Unconsolidated to Consolidated”. Samples exposed to uniaxial compression undergo time-dependent consolidation, or creep, to various degrees. Creep volume strain-time relations obey simple log-time behavior through the range of porosities (~50 to 2% as measured); creep strain rate increases with temperature and applied stress as expected. Mercury porosimetry is used to determine characteristic capillary pressure curves from a series of consolidation tests and show characteristic saturation-capillary pressure curves that follow the common van Genuchten (1978, 1980) formulation at low stresses. Higher capillary pressure data are suspect due to the large potential for sample damage, including fluid inclusion decrepitation and pore collapse. Data are supportive of use of the Leverett “J” function (Leverett, 1941) to use for scaling characteristic curves at different degrees of consolidation, but better permeability determinations are needed to support this hypothesis. Recommendations for further and refined testing are made with the goal of developing a self- consistent set of constitutive laws for granular salt consolidation and multiphase (brine-air) flow.

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

EOS7CA Version 1.0: TOUGH2 Module for Gas Migration in Shallow Subsurface Porous Media Systems

EOS7CA is a TOUGH2 module for mixtures of a non-condensible gas (NCG) and air with or without a gas tracer, an aqueous phase, and water vapor. The user can select the NCG as being CO 2 , N 2 , or CH 4 . EOS7CA uses a cubic equation of state with a multiphase version of Darcy’s Law to model flow and transport of gas and aqueous phase mixtures over a range of pressures and temperatures appropriate to shallow subsurface porous media systems. The limitation to shallow systems arises from the use of Henry’s Law for gas solubility which is appropriate for low pressures but begins to over-predict solubility starting at pressures greater than approximately 1 MPa (10 bar). The components modeled in EOS7CA are water, brine, NCG, gas tracer, air, and optional heat. The real gas properties module (ZEVCA) has options for Peng-Robinson, Redlich-Kwong, or Soave-Redlich-Kwong equations of state to calculate gas mixture density, enthalpy departure, and viscosity. Transport of the gaseous and dissolved components is by advection and Fickian molecular diffusion. This user guide provides instructions for use and two sample problems as verification and demonstration of EOS7CA.

03 NATURAL GAS↗
DOE OSTI2015

V&V Of CFD Modeling Of The Argonne Bubble Experiment: FY15 Summary Report

In support of the development of accelerator-driven production of the fission product Mo 99, computational fluid dynamics (CFD) simulations of an electron-beam irradiated, experimental-scale bubble chamber have been conducted in order to aid in interpretation of existing experimental results, provide additional insights into the physical phenomena, and develop predictive thermal hydraulic capabilities that can be applied to full-scale target solution vessels. Toward that end, a custom hybrid Eulerian-Eulerian-Lagrangian multiphase solver was developed, and simulations have been performed on high-resolution meshes. Good agreement between experiments and simulations has been achieved, especially with respect to the prediction of the maximum temperature of the uranyl sulfate solution in the experimental vessel. These positive results suggest that the simulation methodology that has been developed will prove to be suitable to assist in the development of full-scale production hardware.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2015

Implementation and Re nement of a Comprehensive Model for Dense Granular Flows

Dense granular ows are ubiquitous in both natural and industrial processes. They manifest three di erent ow regimes, each exhibiting its own dependence on solids volume fraction, shear rate, and particle-level properties. This research project sought to develop continuum rheological models for dense granular ows that bridges multiple regimes of ow, implement them in open-source platforms for gas-particle ows and perform test simulations. The rst phase of the research covered in this project involved implementation of a steady- shear rheological model that bridges quasi-static, intermediate and inertial regimes of ow into MFIX (Multiphase Flow with Interphase eXchanges - a general purpose computer code developed at the National Energy Technology Laboratory). MFIX simulations of dense granular ows in hourglass-shaped hopper were then performed as test examples. The second phase focused on formulation of a modi ed kinetic theory for frictional particles that can be used over a wider range of particle volume fractions and also apply for dynamic, multi- dimensional ow conditions. To guide this work, simulations of simple shear ows of identical mono-disperse spheres were also performed using the discrete element method. The third phase of this project sought to develop and implement a more rigorous treatment of boundary e ects. Towards this end, simulations of simple shear ows of identical mono-disperse spheres con ned between parallel plates were performed and analyzed to formulate compact wall boundary conditions that can be used for dense frictional ows at at frictional boundaries. The fourth phase explored the role of modest levels of cohesive interactions between particles on the dense phase rheology. The nal phase of this project focused on implementation and testing of the modi ed kinetic theory in MFIX and running bin-discharge simulations as test examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Enhancements to the TOUGH2 Simulator as Implemented in iTOUGH2

iTOUGH2 is a program for parameter estimation, sensitivity analysis, and uncertainty propagation analysis. It is based on the TOUGH2 simulator for non-isothermal multiphase, multicomponent flow and transport in fractured and porous media [Pruess, 1987, 1991, 2005, 2011; Falta et al., 1995; Pruess et al., 1999, 2002, 2012; Doughty, 2013]. The core of iTOUGH2 contains slightly modified versions of TOUGH2 modules. Most code modifications are editorial and do not affect the simulation results. As a result, standard TOUGH2 input files can be used in iTOUGH2, and identical results are obtained if iTOUGH2 is run in forward mode. However, a number of modifications have been made as described in this report. They enhance the functionality, flexibilitu, and eas-of-use of the forward simulator. This report complements the reports iTOUGH2 User's Guide, iTOUGH2 Command Referecne, and the collection of tutorial examples in iTOUGH2 Sample Problems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Final report of “A Detailed Study of the Physical Mechanisms Controlling CO2-Brine Capillary Trapping in the Subsurface” (University of Arizona, DE-SC0006696)

Carbon capture and storage (CCS) of carbon dioxide emissions generated by production or combustion of fossil fuels is a technologically viable means to reduce the build-up of CO2 in the atmosphere and oceans. Using advantages of scale and location, CCS is particularly suitable for large point sources near ubiquitous deep saline aquifers, depleted gas reservoirs, or at production reservoirs for enhanced oil recovery (EOR). In the BES-funded research project, Oregon State University (OSU) carried out capillary trapping experiments with proxy fluids that mimic the properties of the scCO2/brine system under ambient temperatures and pressures, and successfully developed a unique and novel x-ray compatible, high-pressure, elevated temperature setup to study the scCO2/brine system under challenging reservoir conditions. Both methodologies were applied to a variety of porous media, including synthetic (glass bead) and geologic (Bentheimer sandstone) materials. The University of Arizona (UA) developed pore-scale lattice Boltzmann (LB) models which are able to handle the experimental conditions for proxy fluids, as well as the scCO2/brine system, that are capable of simulating permeability in volumes of tens of millions of fluid elements. We reached the following summary findings (main institute indicated): 1. (OSU/UA) To understand capillary trapping in a multiphase fluid-porous medium system, the system must be analyzed from a pore-scale force balance perspective; trapping can be enhanced by manipulating wetting and nonwetting phase fluid properties. 2. (OSU) Pore-scale fluid connectivity and topology has a clear and direct effect on nonwetting phase capillary trapping efficiency. 3. (OSU) Rock type and flow regime also have a pronounced effects on capillary trapping. 4. (OSU/UA) There is a predictable relationship between NWP connectivity and NWP saturation, which allows for development of injection strategies that optimize trapping. The commonly used Land model (Land, 1968) does not predict amount of trapped NWP accurately. 5. (UA) There are ambiguities regarding the segmentation of large-volume gray-scale CT data into pore-volumes suitable for pore-scale modeling. Simulated permeabilities vary by three orders of magnitude and do not resemble observed values very well. Small-volume synchrotron-based CT data (such as produced by OSU) does not suffer significantly from segmentation ambiguities. 6. (UA) A standard properly parameterized Shan-Chen model LB model is useful for simulating porous media with proxy fluids as well as the scCO2/brine system and produces results that are consistent with tomographic observations. 7. (UA) A LB model with fluid-interactions defined by a (modified) Peng-Robinson Equation of State is able to handle the scCO2/brine system with variable solid phase wettability. This model is numerically stable at temperatures between 0 and 250 °C and pressures between 3 and 50 MPa, and produces appropriate densities above the critical point of CO2 and exhibits three-phase separation below. Based on above findings OSU and UA have proposed continued experimentation and pore-scale modeling of the scCO2/brine system. The reported research has extensively covered capillary trapping using proxy fluids, but due to limited beam-time availability we were unable to apply our high-pressure CO2 setup to sufficient variation in fluid properties, and initial scCO2 connectivity. New data will also allow us to test, calibrate and apply our LB models to reservoir conditions beyond those that are currently feasible experimentally. Such experiments and simulations will also allow us to provide information how suitable proxy fluids are for the scCO2/brine system. We believe it would be worthwhile to pursue the following new research questions: 1. What are the fundamental differences in the physics underlying capillary trapping at ambient vs. supercritical conditions? 2. Do newly developed pore-scale trapping interactions and relationships translate to continuum scales? A motivation for these questions was elaborated in “Capillary Trapping of Super-Critical CO2: Linking Pore and Continuum Scales to Verify new Relationships” that was submitted to DOE-BES in 2015.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2016

Crenulation-­1 Flash Report

We briefly report on the first of PHELIX driven Crenulation experiments diagnosed with proton radiography. PHELIX is a 300 kJ capacitor bank located at the LANL LANSCE pRad facility. It is capable of delivering a 4 MA, 10 us current pulse to a low inductance cylindrical load. Using a magnetically driven cylindrical liner, we have shocked a concentric tin cylindrical shell to melt-on-release and compared to theory and calculations. The inner surface of the tin cylinder has three sectors of single-mode perturbations, which are subject to Richtmyer-Meshkov Instability (RMI). Recent theoretical work on the EOS of tin has modified both the Hugoniot and the isentropes for release into various states in tabular data. The new multiphase EOS for tin, SESAME 2161, includes the beta and gamma solid phases as well as a liquid phase. It predicts a lower pressure boundary for release to pure solid (~20 GPa) and a higher pressure boundary for release to pure liquid (~35 GPa) than the existing SESAME 2160 table.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Investigation of Coal-biomass Catalytic Gasification using Experiments, Reaction Kinetics and Computational Fluid Dynamics

A collaborative effort involving experiments, kinetic modeling, and computational fluid dynamics (CFD) was used to understand co-gasification of coal-biomass mixtures. The overall goal of the work was to determine the key reactive properties for coal-biomass mixed fuels. Sub-bituminous coal was mixed with biomass feedstocks to determine the fluidization and gasification characteristics of hybrid poplar wood, switchgrass and corn stover. It was found that corn stover and poplar wood were the best feedstocks to use with coal. The novel approach of this project was the use of a red mud catalyst to improve gasification and lower gasification temperatures. An important results was the reduction of agglomeration of the biomass using the catalyst. An outcome of this work was the characterization of the chemical kinetics and reaction mechanisms of the co-gasification fuels, and the development of a set of models that can be integrated into other modeling environments. The multiphase flow code, MFIX, was used to simulate and predict the hydrodynamics and co-gasification, and results were validated with the experiments. The reaction kinetics modeling was used to develop a smaller set of reactions for tractable CFD calculations that represented the experiments. Finally, an efficient tool was developed, MCHARS, and coupled with MFIX to efficiently simulate the complex reaction kinetics.

09 BIOMASS FUELS↗
DOE OSTI2016

Modeling Manufacturing Impacts on Aging and Reliability of Polyurethane Foams

Polyurethane is a complex multiphase material that evolves from a viscous liquid to a system of percolating bubbles, which are created via a CO2 generating reaction. The continuous phase polymerizes to a solid during the foaming process generating heat. Foams introduced into a mold increase their volume up to tenfold, and the dynamics of the expansion process may lead to voids and will produce gradients in density and degree of polymerization. These inhomogeneities can lead to structural stability issues upon aging. For instance, structural components in weapon systems have been shown to change shape as they age depending on their molding history, which can threaten critical tolerances. The purpose of this project is to develop a Cradle-to-Grave multiphysics model, which allows us to predict the material properties of foam from its birth through aging in the stockpile, where its dimensional stability is important.

36 MATERIALS SCIENCE↗
DOE OSTI2016

TOGA: A TOUGH code for modeling three-phase, multi-component, and non-isothermal processes involved in CO 2 -based Enhanced Oil Recovery

TOGA is a numerical reservoir simulator for modeling non-isothermal flow and transport of water, CO 2 , multicomponent oil, and related gas components for applications including CO 2 -enhanced oil recovery (CO 2 -EOR) and geologic carbon sequestration in depleted oil and gas reservoirs. TOGA uses an approach based on the Peng-Robinson equation of state (PR-EOS) to calculate the thermophysical properties of the gas and oil phases including the gas/oil components dissolved in the aqueous phase, and uses a mixing model to estimate the thermophysical properties of the aqueous phase. The phase behavior (e.g., occurrence and disappearance of the three phases, gas + oil + aqueous) and the partitioning of non-aqueous components (e.g., CO 2 , CH 4 , and n-oil components) between coexisting phases are modeled using K-values derived from assumptions of equal-fugacity that have been demonstrated to be very accurate as shown by comparison to measured data. Models for saturated (water) vapor pressure and water solubility (in the oil phase) are used to calculate the partitioning of the water (H 2 O) component between the gas and oil phases. All components (e.g., CO 2 , H 2 O, and n hydrocarbon components) are allowed to be present in all phases (aqueous, gaseous, and oil). TOGA uses a multiphase version of Darcy’s Law to model flow and transport through porous media of mixtures with up to three phases over a range of pressures and temperatures appropriate to hydrocarbon recovery and geologic carbon sequestration systems. Transport of the gaseous and dissolved components is by advection and Fickian molecular diffusion. New methods for phase partitioning and thermophysical property modeling in TOGA have been validated against experimental data published in the literature for describing phase partitioning and phase behavior. Flow and transport has been verified by testing against related TOUGH2 EOS modules and CMG. The code has also been validated against a CO 2 -EOR experimental core flood involving flow of three phases and 12 components. Results of simulations of a hypothetical 3D CO 2 -EOR problem involving three phases and multiple components are presented to demonstrate the field-scale capabilities of the new code. This user guide provides instructions for use and sample problems for verification and demonstration.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

STAR-CCM+ Verification and Validation Plan

The commercial Computational Fluid Dynamics (CFD) code STAR-CCM+ provides general purpose finite volume method solutions for fluid dynamics and energy transport. This document defines plans for verification and validation (V&V) of the base code and models implemented within the code by the Consortium for Advanced Simulation of Light water reactors (CASL). The software quality assurance activities described herein are port of the overall software life cycle defined in the CASL Software Quality Assurance (SQA) Plan [Sieger, 2015]. STAR-CCM+ serves as the principal foundation for development of an advanced predictive multi-phase boiling simulation capability within CASL. The CASL Thermal Hydraulics Methods (THM) team develops advanced closure models required to describe the subgrid-resolution behavior of secondary fluids or fluid phases in multiphase boiling flows within the Eulerian-Eulerian framework of the code. These include wall heat partitioning models that describe the formation of vapor on the surface and the forces the define bubble/droplet dynamic motion. The CASL models are implemented as user coding or field functions within the general framework of the code. This report defines procedures and requirements for V&V of the multi-phase CFD capability developed by CASL THM. Results of V&V evaluations will be documented in a separate STAR-CCM+ V&V assessment report. This report is expected to be a living document and will be updated as additional validation cases are identified and adopted as part of the CASL THM V&V suite.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

HPC Institutional Computing Project: W15_lesreactiveflow KIVA-hpFE Development: A Robust and Accurate Engine Modeling Software

KIVA-hpFE is a high performance computer software for solving the physics of multi-species and multiphase turbulent reactive flow in complex geometries having immersed moving parts. The code is written in Fortran 90/95 and can be used on any computer platform with any popular complier. The code is in two versions, a serial version and a parallel version utilizing MPICH2 type Message Passing Interface (MPI or Intel MPI) for solving distributed domains. The parallel version is at least 30x faster than the serial version and much faster than our previous generation of parallel engine modeling software, by many factors. The 5th generation algorithm construction is a Galerkin type Finite Element Method (FEM) solving conservative momentum, species, and energy transport equations along with two-equation turbulent model k-ω Reynolds Averaged Navier-Stokes (RANS) model and a Vreman type dynamic Large Eddy Simulation (LES) method. The LES method is capable modeling transitional flow from laminar to fully turbulent; therefore, this LES method does not require special hybrid or blending to walls. The FEM projection method also uses a Petrov-Galerkin (P-G) stabilization along with pressure stabilization. We employ hierarchical basis sets, constructed on the fly with enrichment in areas associated with relatively larger error as determined by error estimation methods. In addition, when not using the hp-adaptive module, the code employs Lagrangian basis or shape functions. The shape functions are constructed for hexahedral, prismatic and tetrahedral elements. The software is designed to solve many types of reactive flow problems, from burners to internal combustion engines and turbines. In addition, the formulation allows for direct integration of solid bodies (conjugate heat transfer), as in heat transfer through housings, parts, cylinders. It can also easily be extended to stress modeling of solids, used in fluid structure interactions problems, solidification, porous media modeling and magneto hydrodynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Ni and Cr addition to alloy waste forms to reduce radionuclide environmental releases

Reference alloy waste forms (RAW) were fabricated and underwent hybrid corrosion/immersion testing to parameterize the ANL analytical oxidative-dissolution model to enable the calculation of fractional release rates and to determine the effectiveness of Ni and Cr trim additions in reducing release rates of radionuclide surrogates. Figure 1 shows the prototypical multiphase microstructure of the alloys with each phase type contributing about equally to the exposed surface area. The waste forms tested at SRNL were variations of the RAW-6 formulation that uses HT9 as the main alloy component, and are meant to enable evaluation of the impact of Ni and Cr trim additions on the release rates of actinides and Tc-99. The test solutions were deaerated alkaline and acidic brines, ranging in pH 3 to pH 10, representing potential repositories with those conditions. The testing approach consisted of 4 major steps; 1) bare surface corrosion measurements at pH values of 3, 5, 8, and 10, 2) hybrid potentiostatic hold/exposure measurements at pH 3, 3) measurement of radionuclide concentrations and relations to anodic current from potentiostatic holds, and 4) identification of corroding phases using SEM/EDS of electrodes.

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

High Fidelity BWR Fuel Simulations

This report describes the Consortium for Advanced Simulation of Light Water Reactors (CASL) work conducted for completion of the Thermal Hydraulics Methods (THM) Level 3 milestone THM.CFD.P13.03: High Fidelity BWR Fuel Simulation. High fidelity computational fluid dynamics (CFD) simulation for Boiling Water Reactor (BWR) was conducted to investigate the applicability and robustness performance of BWR closures. As a preliminary study, a CFD model with simplified Ferrule spacer grid geometry of NUPEC BWR Full-size Fine-mesh Bundle Test (BFBT) benchmark has been implemented. Performance of multiphase segregated solver with baseline boiling closures has been evaluated. Although the mean values of void fraction and exit quality of CFD result for BFBT case 4101-61 agreed with experimental data, the local void distribution was not predicted accurately. The mesh quality was one of the critical factors to obtain converged result. The stability and robustness of the simulation was mainly affected by the mesh quality, combination of BWR closure models. In addition, the CFD modeling of fully-detailed spacer grid geometry with mixing vane is necessary for improving the accuracy of CFD simulation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2017

CFD-DEM Onset of Motion Analysis for Application to Bed Scour Risk Assessment

This CFD study with DEM was done as a part of the Federal Highway Administration’s (FHWA’s) effort to improve scour design procedures. The Computational Fluid Dynamics-Discrete Element Method (CFD-DEM) model, available in CD-Adapco’s StarCCM+ software, was used to simulate multiphase systems, mainly those which combine fluids and solids. In this method the motion of discrete solids is accounted for by DEM, which applies Newton's laws of motion to every particle. The flow of the fluid is determined by the local averaged Navier–Stokes equations that can be solved using the traditional CFD approach. The interactions between the fluid phase and solids phase are modeled by use of Newton's third law. The inter-particle contact forces are included in the equations of motion. Soft-particle formulation is used, which allows particles to overlap. In this study DEM was used to model separate sediment grains and spherical particles laying on the bed with the aim to analyze their movement due to flow conditions. Critical shear stress causing the incipient movement of the sediment was established and compared to the available experimental data. An example of scour around a cylindrical pier is considered. Various depths of the scoured bed and flow conditions were taken into account to gain a better understanding of the erosion forces existing around bridge foundations. The decay of these forces with increasing scour depth was quantified with a ‘decay function’, which shows that particles become increasingly less likely to be set in motion by flow forces as a scour hole increases in depth. Computational and experimental examples of the scoured bed around a cylindrical pier are presented.

97 MATHEMATICS AND COMPUTING↗