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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

18 records

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Acoustic fluidization and the scale dependence of impact crater morphology

A phenomenological Bingham plastic model has previously been shown to provide an adequate description of the collapse of impact craters. This paper demonstrates that the Bingham parameters may be derived from a model in which acoustic energy generated during excavation fluidizes the rock debris surrounding the crater. Experimental support for the theoretical flow law is presented. Although the Bingham yield stress cannot be computed without detailed knowledge of the initial acoustic field, the Bingham viscosity is derived from a simple argument which shows that it increases as the 3/2 power of crater diameter, consistent with observation. Crater collapse may occur in material with internal dissipation Q as low as 100, comparable to laboratory observations of dissipation in granular materials. Crater collapse thus does not require that the acoustic field be regenerated during flow.

Melosh, H. J.↗

DOE OSTI2013

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↗

NASA NTRS1978

Yield strengths of flows on the earth, Mars, and moon

Dimensions of flows on the earth, Mars, and moon and their topographic gradients obtained from remote measurements are used to calculate yield strengths with a view to explore the validity of the Bingham plastic model and determine whether there is a relation between yield strengths and silica contents. Other factors are considered such as the vagaries of natural phenomena that might contribute to erroneous interpretations and measurements. Comparison of yield strengths of Martian and lunar flows with terrestrial flows suggests that the Martian and lunar flows are more akin to terrestrial basalts than they are to terrestrial andesites, trachytes, and rhyolites.

Moore, H. J.↗

NASA NTRS1991

Fluid volcanism on Miranda and Ariel - Flow morphology and composition

Several types of volcanic units have been recognized on the icy Uranian satellites Miranda and Ariel. On Miranda, ridges characterized by crest grooves are up to 10 km wide and 500 m high. A continuous flat-topped flow band also 10 km wide and 500 m high forms the outer southern margin of Elsinore Corona, which appears to comprise coalesced flow bands and ridges. On Ariel, in addition to at least one ridge unit similar to those on Miranda, flood plains material has covered the floors of deep chasmata (grabens) and an irregular depression. Flows on both satellites are characterized by linear vent geometries and distinct topographic margins, which indicate extrusion of a relatively viscous material. The topography of the flows can be used to estimate flow viscosity or yield strength using a Bingham plastic model. Extrusion viscosity estimates, incorporating plausible volcanologically based emplacement time scales and a rigid crust correction, range from 10 MP to 1 GP (10 TP in the unlikely absence of a chilled crust). Viscosity estimates are dependent on the assumed emplacement time scale, however, and could be as high as 10 PP, if a solid-state-based time scale is assumed.

Schenk, Paul M.↗

DOE OSTI2012

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↗

DOE OSTI2022

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↗

DOE OSTI2013

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↗

DOE OSTI2018

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↗

DOE OSTI2019

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↗

NASA NTRS1992

Fluid outflows from Venus impact craters - Analysis from Magellan data

Many impact craters on Venus have unusual outflow features originating in or under the continuous ejecta blankets and continuing downhill into the surrounding terrain. These features clearly resulted from flow of low-viscosity fluids, but the identity of those fluids is not clear. In particular, it should not be assumed a priori that the fluid is an impact melt. A number of candidate processes by which impact events might generate the observed features are considered, and predictions are made concerning the rheological character of flows produce by each mechanism. A sample of outflows was analyzed using Magellan images and a model of unconstrained Bingham plastic flow on inclined planes, leading to estimates of viscosity and yield strength for the flow materials. It is argued that at least two different mechanisms have produced outflows on Venus: an erosive, channel-forming process and a depositional process. The erosive fluid is probably an impact melt, but the depositional fluid may consist of fluidized solid debris, vaporized material, and/or melt.

Asimow, Paul D.↗

NASA NTRS1984

Dynamics of lava flow - Thickness growth characteristics of steady two-dimensional flow

The thickness growth characteristics of flowing lava are investigated using a heat balance model and a two-dimensional model for flow of a Bingham plastic fluid down an inclined plane. It is found that yield strength plays a crucial role in the thickening of a lava flow of given flow rate. To illustrate this point, downstream thickness profiles and yield strength distributions were calculated for flows with mass flow rates of 10,000 and 100,000 kg/m-sec. Higher flow rates led to slow cooling rates which resulted in slow rate of increase of yield strength and thus greater flow lengths.

Park, S.↗

DOE OSTI2015

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↗

DOE OSTI2021

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↗

NASA NTRS1982

A schematic model of crater modification by gravity

The morphology of craters found on planets and moons of the solar system is examined and a development model which can account for the observed crater characteristics is discussed. The prompt collapse of craters to form flat floors, terraced walls, and central peak structures is considered to be the result of an approximate Bingham plastic rheology of the material surrounding the crater. This rheology is induced dynamically by the strong incoherent acoustic 'noise' accompanying excavation of the crater. Central pits, peak rings, and other multiple symmetric-profile rings originate by oscillation of this fluid. Large craters with transient depths comparable to the lithosphere thickness are subject to collapse by fragmentation of the lithosphere as well as fluidization. The considered concepts are developed mathematically. A model emerges which appears capable of explaining most of the qualitative features of large impact structures.

Melosh, H. J.↗

DOE OSTI2018

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↗

NASA NTRS1993

The origin of Venusian channels: Modelling of thermal erosion by lava

Magellan imagery has revealed that channels, apparently volcanic in origin, are abundant on the surface of Venus. There has been much debate about the origin of these channels. Are they the result of erosional (either thermal or mechanical) or constructional processes? A common characteristic of the simple sinuous channels is that they show evidence of erosion near their source and then become purely constructional, forming levees and in some cases roofing over completely. One method of showing that thermal erosion is capable of producing the type of channels seen is to use computer modeling incorporating the physical conditions on Venus and the physical characteristics of the different types of lava that may have been erupted. It is possible to calculate, relatively easily, two channel parameters. The first is the erosion rate, which combined with eruption duration, gives depth. The second is for how long after leaving the source the erupted lava will continue to be capable of thermal erosion before constructional processes dominate. Making assumptions about the rheology of the lava (e.g., assume it behaves as a Bingham plastic) along with the slope angle yields a flow velocity and therefore a distance over which thermal erosion will take place. Due to the resolution (both vertical and horizontal) of the Magellan altimetric data, the distance from the source that the channel is erosional can be much more accurately measured than the depth of the channel. This will remain the case until stereo imagery becomes available for large areas of the planet.

Bussey, D. B. J.↗

DOE OSTI2018

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↗

DOE OSTI2023

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)↗