Search NASA⌕ Search

SEARCH · Search NASA

Results for “Bingham Plastic Model”

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.

A study of grout flow pattern analysis

A new disposal unit, designated as Salt Disposal Unit no. 6 (SDU6), is being designed for support of site accelerated closure goals and salt nuclear waste projections identified in the new Liquid Waste System plan. The unit is cylindrical disposal vault of 380 ft diameter and 43 ft in height, and it has about 30 million gallons of capacity. Primary objective was to develop the computational model and to perform the evaluations for the flow patterns of grout material in SDU6 as function of elevation of grout discharge port, and slurry rheology. A Bingham plastic model was basically used to represent the grout flow behavior. A two-phase modeling approach was taken to achieve the objective. This approach assumes that the air-grout interface determines the shape of the accumulation mound. The results of this study were used to develop the design guidelines for the discharge ports of the Saltstone feed materials in the SDU6 facility. The focusing areas of the modeling study are to estimate the domain size of the grout materials radially spread on the facility floor under the baseline modeling conditions, to perform the sensitivity analysis with respect to the baseline design and operating conditions such as elevation of discharge port, discharge pipe diameter, and grout properties, and to determine the changes in grout density as it is related to grout drop height. An axi-symmetric two-phase modeling method was used for computational efficiency. Based on the nominal design and operating conditions, a transient computational approach was taken to compute flow fields mainly driven by pumping inertia and natural gravity. Detailed solution methodology and analysis results are discussed here.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Saltstone SDU6 Modeling Study

A new disposal unit, designated as Saltstone Disposal Unit 6 (SDU6), is being designed for support of site accelerated closure goals and salt waste projections identified in the new Liquid Waste System Plan. The unit is a cylindrical disposal cell of 375 ft in diameter and 43 ft in height, and it has a minimum 30 million gallons of capacity. SRNL was requested to evaluate the impact of an increased grout placement height on the flow patterns radially spread on the floor and to determine whether grout quality is impacted by the height. The primary goals of the work are to develop the baseline Computational Fluid Dynamics (CFD) model and to perform the evaluations for the flow patterns of grout material in SDU6 as a function of elevation of grout discharge port and grout rheology. Two transient grout models have been developed by taking a three-dimensional multiphase CFD approach to estimate the domain size of the grout materials radially spread on the facility floor and to perform the sensitivity analysis with respect to the baseline design and operating conditions such as elevation height of the discharge port and fresh grout properties. For the CFD modeling calculations, air-grout Volume of Fluid (VOF) method combined with Bingham plastic and time-dependent grout models were used for examining the impact of fluid spread performance for the initial baseline configurations and to evaluate the impact of grout pouring height on grout quality. The grout quality was estimated in terms of the air volume fraction for the grout layer formed on the SDU6 floor, resulting in the change of grout density. The study results should be considered as preliminary scoping analyses since benchmarking analysis is not included in this task scope. Transient analyses with the Bingham plastic model were performed with the FLUENT TM code on the high performance parallel computing platform in SRNL. The analysis coupled with a transient grout aging model was performed by using ANSYS-CFX code in the parallel computing platform in Mercer University. Recommended operational guidance was developed assuming that local shear rates and flow patterns related to radial spread along the SDU floor can be used as a measure of grout performance and spatial dispersion affected by the grout height and viscosity. The main conclusions drawn from the grout modeling and calculations are as follows: • The 5 ft height baseline results show that when the 150 gpm grout flow with a 5 Pa yield stress and a 60 cp viscosity is poured down through a 3 inch discharge port, the grout is spread radially up to about 64 ft distance from the pouring center after 2 hours’ pouring time. The air volume fraction of the grout layer is about 29% at 5 minutes’ transient time, and it is reduced by about 9% in 2 hours’ pouring time, resulting in the grout density consisting of about 80% grout and 20% air volume fractions. • The sensitivity results show that when the discharge port is located at a higher position, a larger amount of air is trapped inside the layer formed below the discharge port at the early transient time of less than 30 minutes because of the higher impinging momentum of the grout flow on the floor, resulting in the formation of less smooth layer. • The results clearly indicate that the radial spread for the 43 ft discharge port is about 10% faster than that of the 5 ft discharge port for the early transient period of 5 minutes. However, for the pouring time longer than half an hour, the discharge port height does not affect the radial distance spread on the disposal floor. • When grout quality is related to grout volume fraction, the grout volume fraction for the 43 ft discharge port has lower volume fraction than the 5 ft discharge port for the transient period of the first 5 minutes. However, for the pouring time longer than half an hour, the discharge port height does not affect the grout volume fraction for the layer accumulated on the disposal floor. • A modified Bingham plastic model coupled with time-dependent viscosity behavior was developed for conducting the initial scoping calculations to assess the impact of fluid residence time on radial spreading and basic flow patterns. • The results for the transient viscosity model show that when grout material becomes more viscous, the thickness of the grout layer accumulated on the floor becomes higher, but the radial distance spread on the horizontal floor becomes smaller. The early transient results for the grout density with about 32% air volume fractions are in reasonable agreement with those of the idealized Bingham plastic model. • It is recommended that the current models developed here be benchmarked against the experimental results for critical applications of the modeling results.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Modeling gas release from a Bingham plastic slurry and deconvoluting measured data

Physics-based models for predicting the off-gassing characteristics of agitated slurries are presented. This approach decomposes the system into two separate but coupled slurry and headspace models. The physics driving bubble transport through the slurry and gas mixing within the headspace are discussed. An analytical expression for predicting the time evolution of the headspace concentration using first principles theory is also presented. Predictions from the numerical models, as well as expectations from the analytical solution, both agree with measured off gassing data for two different experimental operating conditions. After benchmarking the numerical model predictions against experimental data, the model was used to make predictions for gas release from a full-scale vessel and to perform sensitivity analyses to examine the sensitivity of gas release to parameters such as yield stress, consistency index, slurry density, bubble size, bubble concentration, and impeller speed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Review Of Rheology Models For Hanford Waste Blending

The area of rheological property prediction was identified as a technology need in the Hanford Tank Waste – waste feed acceptance initiative area during a series of technical meetings among the national laboratories, Department of Energy-Office of River Protection, and Hanford site contractors. Meacham et al. delivered a technical report in June 2012, RPP-RPT-51652 “One System Evaluation of Waste Transferred to the Waste Treatment Plant” that included estimating of single shell tank waste Bingham plastic rheological model constants along with a discussion of the issues inherent in predicting the rheological properties of blended wastes. This report was selected as the basis for moving forward during the technical meetings. The report does not provide an equation for predicting rheological properties of blended waste slurries. The attached technical report gives an independent review of the provided Hanford rheological data, Hanford rheological models for single tank wastes, and Hanford rheology after blending provided in the Meacham report. The attached report also compares Hanford to SRS waste rheology and discusses some SRS rheological model equations for single tank wastes, as well as discussing SRS experience with the blending of waste sludges with aqueous material, other waste sludges, and frit slurries.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Short communication: Estimation of yield stress/viscosity of molten octol

Explosive HMX particles are similar in morphology and chemistry to RDX particles, the main constituent of Composition B-3 (Comp B-3). This suggests molten HMX-TNT formulations may show Bingham plasticity, much like recent studies have shown for Comp B-3. Here a Bingham plastic viscosity model, including yield stress and shear thinning, is presented for octol (70/30wt% HMX/TNT) as a function of HMX particle volume fraction. The effect of HMX dissolution into molten TNT is included in this analysis.

36 MATERIALS SCIENCE↗

Slurry rheology of Hanford sludge

Here we evaluate the slurry viscosity of radioactive waste. Hanford waste is a motely mixture of particle sizes and densities composed of metal and oxide particles with constituents spanning much of the periodic table stored at high ionic strength and high pH. Hanford wastes have traditionally been represented as simply an ideal Bingham plastic, but this rheological approximation often provides a poor fit from 0 to 200 1/s, precisely the range of interest for engineering design of waste transport and mixing systems. Here we use a modified Cross formulation inclusive of a yield stress to fit stress-strain rate data for reconstituted Hanford REDOX sludge waste, measured by increasing shear-rate (up-ramp) and decreasing shear-rate (down-ramp) sweeps. We find that this continuous single equation formulation quantitatively matches experimental data better than other rheological models including the Bingham plastic and Herschel-Bulkley approximations. The increase of viscosity and yield stress from up ramp to down ramp indicates the formation of particle microstructure during shear due to attractive forces.

radioactive waste, high ionic strength slurries, y↗

Modeling a Newtonian Fluid with a Rate-Dependent, Von Mises Plasticity Model for Solid Mechanics Applications

This memo derives a particular use of a rate-dependent, perfectly plastic Von Mises Plasticity constitutive model that represents a Newtonian (or Bingham) fluid. Only two parameters are needed: the viscosity and bulk modulus to complete the model parameterization. Under certain deformations, the model exactly represents a Newtonian fluid. Under others, the model shear thins. A Newtonian fluid representation may be used in confined flow applications in which a fluid-solid interaction is desired but code coupling is not.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Computational Fluid Dynamics Simulation of Compressible Non-Newtonian Biomass in a Compression-Screw Feeder

Compression-screw feeders play a critical role in biorefineries to transport lignocellulosic biomass feedstocks from the feeding hoppers to biomass-conversion reactors in order to pretreat and convert biomass to hydrocarbon liquid biofuels and other power and energy resources. One of the main challenges in the operation of screw feeder is plugging and jamming of compressed biomass with high concentration of insoluble solids. The focus of this paper is to numerically investigate the screw feeder at these challenging operating conditions and help with the optimization of the screw feeder design to avoid operation failure. In this work a customized CFD model based on open-source OpenFOAM package [1] was developed to simulate the concentrated biomass as a highly viscous non-Newtonian fluid in the screw feeder. The biomass is modeled as a single-phase compressible Bingham fluid with a plastic viscosity as well as a density-dependent yield stress. The compressibility formulation (pressure-dependent density) and the density-dependent yield stress formulation in the governing equations follow the suggestions from a recent study by Duncan et al. [2]. A pilot-scale hopper/screw feeding system at NREL [3] is used to compare the experimental observations with our simulation results. The auger is 280 mm long and tapered with outer diameter changing from 80 mm to 35 mm. The auger rotates from 10 to 60 rpm in a conical throat which contains anti-rotational bars. The simulations predicted the required torque for the screw feeder and the pressure increase at the exit for biomass feedstocks with various fluid viscosity properties and auger rotating speeds. The analysis of the stress forces helped to identify the critical conditions were the screw feeder excessive wear or jamming could occur.

biofuels↗

Integrated rheology model: Explosive Composition B-3

Composition B-3 (Comp B-3) is a high explosive formulation composed of 60/40wt% RDX (1,3,5-trinitroperhydro-1,3,5-triazine) /TNT (2,4,6 trinitrotoluene). Above approximately 78°C this formulation partially melts to form a multiphase system with solid RDX particles in a molten TNT matrix. This multiphase system presents a number of phenomena that influence its apparent viscosity. In an earlier study explosive Composition B-3 (Comp B-3, 60/40wt% RDX/TNT) was examined for evidence of yield stress using a non-isothermal falling ball viscometer and a yield stress model was proposed in this paper. An integrated viscosity model suitable for use in computational fluid dynamics (CFD) simulations is developed to capture the transition from a heterogeneous solid to a Bingham viscoplastic fluid. This viscosity model is used to simulate the motion of imbedded spheres falling through molten Comp B-3. Finally, comparison of the simulations to physical tests show agreement between the positions predicted by the model and the measured locations of the spheres as a function of temperature between 90C and 165C.

36 MATERIALS SCIENCE↗

Highly parallelisable simulations of time-dependent viscoplastic fluid flow with structured adaptive mesh refinement

We present the extension of an efficient and highly parallelisable framework for incompressible fluid flow simulations to viscoplastic fluids. The system is governed by incompressible conservation of mass, the Cauchy momentum equation, and a generalised Newtonian constitutive law. In order to simulate a wide range of viscoplastic fluids, we employ the Herschel-Bulkley model for yield-stress fluids with nonlinear stress-strain dependency above the yield limit. We utilise Papanastasiou regularisation in our algorithm to deal with the singularity in apparent viscosity. The resulting system of partial differential equations is solved using the IAMR (Incompressible Adaptive Mesh Refinement) code, which uses second-order Godunov methodology for the advective terms and semi-implicit diffusion in the context of an approximate projection method to solve adaptively refined meshes. By augmenting the IAMR code with the ability to simulate regularised Herschel-Bulkley fluids, we obtain efficient numerical software for time-dependent viscoplastic flow in three dimensions, which can be used to investigate systems not considered previously due to computational expense. We validate results from simulations using this new capability against previously published data for Bingham plastics and power-law fluids in the two-dimensional lid-driven cavity. In doing so, we expand the range of Bingham and Reynolds numbers which have been considered in the benchmark tests. Moreover, extensions to time-dependent flow of Herschel-Bulkley fluids and three spatial dimensions offer new insights into the flow of viscoplastic fluids in this test case, and we provide missing benchmark results for these extensions.

42 ENGINEERING↗

Understanding the transient large amplitude oscillatory shear behavior of yield stress fluids

A full understanding of the sequence of processes exhibited by yield stress fluids under large amplitude oscillatory shearing is developed using multiple experimental and analytical approaches. A novel component rate Lissajous curve, where the rates at which strain is acquired unrecoverably and recoverably are plotted against each other, is introduced and its utility is demonstrated by application to the analytical responses of four simple viscoelastic models. Using the component rate space, yielding and unyielding are identified by changes in the way strain is acquired, from recoverably to unrecoverably and back again. The behaviors are investigated by comparing the experimental results with predictions from the elastic Bingham model that is constructed using the Oldroyd–Prager formalism and the recently proposed continuous model by Kamani, Donley, and Rogers in which yielding is enhanced by rapid acquisition of elastic strain. The physical interpretation gained from the transient large amplitude oscillatory shear (LAOS) data is compared to the results from the analytical sequence of physical processes framework and a novel time-resolved Pipkin space. The component rate figures, therefore, provide an independent test of the interpretations of the sequence of physical processes analysis that can also be applied to other LAOS analysis frameworks. Each of these methods, the component rates, the sequence of physical processes analysis, and the time-resolved Pipkin diagrams, unambigiously identifies the same material physics, showing that yield stress fluids go through a sequence of physical processes that includes elastic deformation, gradual yielding, plastic flow, and gradual unyielding.

Kamani, Krutarth M. (ORCID:0000000338975420)↗