Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Viscous suspensions”

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.

Machine learning methods for particle stress development in suspension Poiseuille flows

Numerical simulations are used to study the dynamics of a developing suspension Poiseuille flow with monodispersed and bidispersed neutrally buoyant particles in a planar channel, and machine learning is applied to learn the evolving stresses of the developing suspension. The particle stresses and pressure develop on a slower time scale than the volume fraction, indicating that once the particles reach a steady volume fraction profile, they rearrange to minimize the contact pressure on each particle. Here we consider how the stress development leads to particle migration, time scales for stress development, and present a new physics-informed Galerkin neural network that allows for learning the particle stresses when direct measurements are not possible. The particle fluxes are compared with the Suspension Balance Model with good agreement. We show that when stress measurements are possible, the MOR-physics operator learning method can also capture the particle stresses.

97 MATHEMATICS AND COMPUTING↗

Bidisperse supension balance model

The suspension balance model (SBM) for viscous Stokes flow has been well studied for the case of identical, monodisperse spherical particles in a channel, however, more work remains to be done to explore the bidisperse and polydisperse SBM. Here, we present a simple extension of the SBM that allows for modeling of suspensions of particles of bidisperse size. The comparison with avaiable experiments and direct simulations is found to be good. Additionally, we present a range of simulations to justify the assumptions made in the polydisperse suspension balance model for non-Brownian, neutrally buoyant particles that vary in size with a bidisperse distribution. The aim is to study the effects of moderate size variation on the rheology of the suspension and the distribution of the particles across the flow, especially in the near-wall regions. It is shown that, at the size ratios considered, bidispersity in the particle size does not significantly affect rheological parameters including the particle volume fraction and relative viscosity in the bulk of a Couette flow. The particle phase stresses are distributed between small and large particles in proportion to the relative volume fractions of each species. In bidisperse suspensions, the large particles in the wall layer develop a spatial structure and form chains in the streamwise direction.

97 MATHEMATICS AND COMPUTING↗

Caustic Neutralization and Precipitation of Acidic Dissolved Simulated Stainless Steel–Clad Plutonium and Plutonium/Uranium Nuclear Fuel

Simulated dissolved stainless steel (SS) clad Pu and Pu/U nuclear fuel in HNO 3 was neutralized to a free hydroxide (OH - ) concentration of 0.6 M. A thermal neutron poison, Gd, was added to the simulants at concentrations of either ~3 - 6 g/L or ~37 - 38 g/L. The supernate Pu concentration the day of neutralization ranged from 0.48 to 8.75 mg/L. The supernate Pu concentration of a simplified simulant neutralized to 0.6 M OH - above precipitated solids containing Pu was demonstrated to decrease over 18 days. A significant portion of precipitated Pu was found to be insoluble in 8 M HNO 3 at ambient temperature, but essentially quantitative Pu dissolution was achieved in 11.5 M HNO 3 /0.1 M KF at 100 °C. The difficulty in dissolving the Pu precipitate is believed to be due to the formation of refractory PuO 2 •xH 2 O during the neutralization process. Initial Gd concentrations of ~37 – 38 g/L were found to result in a greater Al precipitation when neutralized to 0.6 M OH - than initial Gd concentrations of ~3 – 6 g/L. Physical properties of the resultant slurries were measured and used to calculate limiting flowrates and slurry velocities by gravity only in transfer piping between the Savannah River Site’s H-Canyon Facility and the Concentration, Storage, and Transfer Facility (CSTF). These results were compared to calculated deposition velocities to predict if solids would settle during the transfer. The Newtonian model was found to be reasonable for each diluted slurry evaluated. Deposition velocities of Pu containing slurries are lower than nuclear fuel slurries primarily composed of U due to the high density of Pu solids. In conclusion, dilution of slurries reduces the margin between the slurry and deposition velocities due to the reduction in viscosity because higher viscous forces on the particles promote maintained suspension.

Actinide Neutralization↗

Simulation of Particulate Transport for Delivery of Solid Amendments into the Subsurface: FY24 Status Report

For particulate-based amendments to be viable for field-scale remediation at the Hanford Site (e.g., 200 DV-1 Operable Unit), particles need to be delivered a sufficient radial distance from an injection well and retained at concentrations high enough for effective treatment. An accurate description of the particle radius of influence (ROI) is critical for developing an overall remediation strategy. However, field-scale particle simulations are currently limited due to insufficient simulation capabilities and a lack of experimental data to validate and parameterize particle transport models. To help build toward field-scale deployment, this fiscal year (FY) we have (1) developed a pre screening tool to estimate particle transport, (2) implemented particle transport models within PFLOTRAN, and (3) conducted preliminary estimations of particle ROI. While field-scale numerical simulations will ultimately be necessary before remedy design and field implementation, we have developed a pre-screening tool that offers valuable estimations of expected particle injectability and ROI in a 1-D system. The advantage of the tool is that it does not require extensive laboratory experiments and instead makes predictions based solely on routine laboratory measurements. This tool can assist in down-selection and decision-making by identifying which particle amendment systems are worth pursuing in future laboratory experiments, such as 1-D column tests and beyond. With any system, scaling up from the lab to the field presents challenges. Currently, there is no field data available for model calibration or validation. However, the theoretical particle models being developed herein are the best tools available to guide progress toward field deployment. To help bridge this gap and verify model predictions, larger-scale lab experiments are being proposed. To advance simulation capabilities, six particle transport models are being integrated into the reactive transport simulator PFLOTRAN. These include colloid filtration theory (CFT) and five additional particle transport models (M1-M5). Each model, from M1 to M5, progressively incorporates additional particle transport and retention processes. Ultimately, the simplest model capable of accurately describing 1-D column data will be selected and parameterized. During FY24, the CFT and M1 model have been fully implemented within PFLOTRAN. Using an existing 1 D column experiment, the two currently implemented particle transport models (CFT and M1), and associated parameters, were fit to this experiment. While simpler model formulations are helpful for estimations, these formulations could not fully describe particle transport and retention behavior in the previous 1-D column experiment. Thus, additional complexities will need to be considered, which will be accounted for in the M2-M5 model formulations. Additionally, because a viscous, shear thinning fluid was required to keep particles in suspension, considerations for flow will also need to also be accounted for. Therefore, a new immiscible two-phase flow mode is currently being implemented in PFLOTRAN. With some modifications, this new flow module could also support simulation of non-Newtonian liquid amendments, foams, and emulsions. We also estimated the expected ROI of solid amendments using 1-D simulations. The average predicted ROI was approximately 15 ft for micron-sized zero valent iron (mZVI) suspended in xanthan gum (XG). Using the pre screening tool and ROI estimates, additional amendment-delivery laboratory characterization and experiments are proposed. The results from additional experiments can be used to validate and parametrize particulate transport model formulations, which will ultimately provide predictive capabilities for field amendment-delivery systems.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Thixotropy, antithixotropy, and viscoelasticity in hysteresis

Thixotropy, antithixotropy, and viscoelasticity are three types of time-dependent dynamics that involve fundamentally different underlying physical processes. Here, we show that the three dynamics exhibit different signatures in hysteresis by examining the fingerprints of the simplest thixotropic kinetic model, a new antithixotropic model that we introduce here, and the Giesekus model. We start by showing that a consistent protocol to generate hysteresis loops is a discrete shear-rate controlled ramp that begins and ends at high shear rates, rather than at low shear rates. Using this protocol, we identify two distinguishing features in the resulting stress versus shear rate loops. The first is the direction of the hysteresis loops: clockwise for thixotropy, but counterclockwise for viscoelasticity and antithixotropy. A second feature is achieved at high ramping rates where all responses lose hysteresis: the viscoelastic response shows a stress plateau at low shear rates due to lack of stress relaxation, whereas the thixotropic and antithixotropic responses are purely viscous with minimal shear thinning or thickening. We establish further evidence for these signatures by experimentally measuring the hysteresis of Laponite suspensions, carbon black suspensions, and poly(ethylene oxide) solutions, each representing a historically accepted example of each class of material behavior. The signatures measured in experiments are consistent with those predicted by the three models. This study reveals different fingerprints in hysteresis loops associated with thixotropy, antithixotropy, and viscoelasticity, which may be helpful in distinguishing the three time-dependent responses.

Mechanics↗

Filament extension atomizers: A study of integration potential into spray dryers for high-viscosity and non-Newtonian fluids

Filament extension atomizers (FEAs) are an emerging class of spray nozzles designed to atomize high-viscosity and non-Newtonian fluids that are challenging for conventional pressure or two-fluid nozzles. In this study, we investigate the influence of roller geometry, surface velocity, roller material, and fluid rheology on the atomization performance of FEA systems using concentrated whey protein suspensions (50–70 wt. %). Extensional and shear rheology experiments, along with high-speed imaging and particle image velocimetry, reveal that filament breakup dynamics are governed by competition between inertial, capillary, and viscoelastic stresses. High roller rotational velocity leads to narrower spray cones, contradicting rheology experiments and suggests a significant inertial contribution to filament breakup. Smaller rollers operating at the same rotational velocity led to broader spray cones consistent with expectations. An FEA nozzle was integrated into a conventional dryer producing particles in the 100 μm range. Results suggest that FEA technology enables atomization of highly viscous fluids at industrially relevant spray cone angles similar to those generated by pressure nozzles, offering a pathway to improve energy efficiency in spray drying by enabling higher solids loading feedstocks. Furthermore, these insights provide critical guidance for optimizing FEA nozzle designs and process parameters across a range of applications.

Energy efficiency↗