Engineering topics

Khan, Monsurul

Publications and source records attributed to Khan, Monsurul.

7 records

DOE OSTI2024

Rheology of bi-disperse dense fiber suspensions

While significant progress has been made in the modeling and simulation of uniform fiber suspensions, no existing model has been validated for industrially-relevant concentrated suspensions containing fibers of multiple aspect ratios. In the present work, we investigate bi-disperse suspensions with two fiber populations in varying aspect ratios in a steady shear flow using direct numerical simulations. Moreover, we measure the suspension viscosity by creating a controlled length bidispersity for nylon fibers suspended in a Newtonian fluid. Here the results showed good agreement between the experimentally measured and numerically predicted viscosity for bi-disperse suspensions. The ratio between the aspect ratio of large to small fibers (size ratio) and the volume fraction of large fibers (composition) in bi-disperse systems strongly affected the rheological behavior of the suspension. The increment of relative viscosity associated with size ratio and composition can be explained by the decrease in the maximum flowable limit or jamming volume fraction. Moreover, the relative viscosity of bi-disperse suspensions collapses, when plotted against the reduced volume fraction, demonstrating the controlling influence of the jamming fraction in bi-disperse fiber suspensions.

42 ENGINEERING↗

DOE OSTI2023

Rheology of dense fiber suspensions: Origin of yield stress, shear thinning, and normal stress differences

We use high-fidelity computational model of fiber suspension to investigate suspension rheology. We investigate the interplay between the hydrodynamic, noncontact attractive and repulsive, and interfiber contact interactions. The shear-thinning viscosity and finite yield stress obtained from the Immerse Boundary Method simulations align quantitatively with experimental findings from the literature. The study demonstrates that attractive interactions lead to both yield stress and shear thinning behavior in rigid fiber suspensions. This discovery holds significance as it contributes to the ongoing debate on the source of yield stress in fiber suspensions. The proposed model is used to quantify normal stresses in addition to shear thinning and yield stress.

Ardekani, Arezoo↗

DOE OSTI2023

Rheology of dense fiber suspensions: Origin of yield stress, shear thinning, and normal stress differences

We explain the origins of yield stress, shear thinning, and normal stress differences in rigid fiber suspensions. We investigate the interplay between the hydrodynamic, noncontact attractive and repulsive, and interfiber contact interactions. The shear-thinning viscosity and finite yield stress obtained from the Immerse Boundary Method simulations are in quantitative agreement with experiential results from the literature. In this study, we show that attractive interactions result in yield stress and shear thinning rheology in the suspensions of rigid fibers. This is an important finding, given the ongoing discussion regarding the origin of the yield stress for suspensions of fibers. Finally, the ability of the proposed model to quantitatively predict the rheology is not limited to only shear thinning and yield stress but also extends to normal stresses.

42 ENGINEERING↗

DOE OSTI2023

Rheology of concentrated fiber suspensions with a load-dependent friction coefficient

Here we investigate the effects of fiber aspect ratio, roughness, flexibility, and volume fraction on the rheology of concentrated suspensions in a steady shear flow using direct numerical simulations. We model the fibers as inextensible continuous flexible slender bodies with the Euler-Bernoulli beam equation governing their dynamics suspended in an incompressible Newtonian fluid. The fiber dynamics and fluid flow coupling is achieved using the immersed boundary method. In addition, the fiber surface roughness might lead to interfiber contacts, resulting in normal and tangential forces between the fibers, which follow Coulomb’s law of friction. The surface roughness is modeled as hemispherical protrusions on the fiber surfaces. Their deformation results in a normal load-dependent friction coefficient. Our simulations accurately predict the experimentally observed shear thinning in fiber suspensions. Furthermore, we find that the suspension viscosity η increases with increasing the volume fraction, roughness, fiber rigidity, and aspect ratio. The increase in η is the macroscopic manifestation of a similar increase in the microscopic contact contribution to the total stress with these parameters. In addition, we observe positive and negative first N 1 and second N 2 normal stress differences, respectively, with |N 2 | < |N 1 |, in agreement with previous experiments. Last, we propose a modified Maron-Pierce law to quantify the reduction in the jamming volume fraction by increasing the fiber aspect ratio and roughness. Our results and analysis establish the use of fiber surface tribology to tune the suspension flow behavior.

42 ENGINEERING↗

DOE OSTI2023

A constitutive model for sheared dense fiber suspensions

Here we propose a constitutive model to predict the viscosity of fiber suspensions, which undergoes shear thinning, at various volume fractions, aspect ratios, and shear stresses/rates. We calibrate the model using the data from direct numerical simulation and prove the accuracy by predicting experimental measurements from the literature. We use a friction coefficient decreasing with the normal load between the fibers to quantitatively reproduce the experimentally observed shear thinning in fiber suspensions. In this model, the effective normal contact force, which is directly proportional to the bulk shear stress, determines the effective friction coefficient. A rise in the shear stress reduces the effective friction coefficient in the suspension. As a result, the jamming volume fraction increases with the shear stress, resulting in a shear thinning in the suspension viscosity. Moreover, we extend the model to quantify the effects of fiber volume fraction and aspect ratio in the suspension. We calibrate this model using the data from numerical simulations for the rate-controlled shear flow. Once calibrated, we show that the model can be used to predict the relative viscosity for different volume fractions, shear stresses, and aspect ratios. The model predictions are in excellent agreement with the available experimental measurements from the literature. The findings of this study can potentially be used to tune the fiber size and volume fraction for designing the suspension rheology in various applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Multi-fidelity modeling to predict the rheological properties of a suspension of fibers using neural networks and Gaussian processes

Unveiling the rheological properties of fiber suspensions is of paramount interest to many industrial applications. There are multiple factors, such as fiber aspect ratio and volume fraction, that play a significant role in altering the rheological behavior of suspensions. Three-dimensional (3D) numerical simulations of coupled differential equations of the suspension of fibers are computationally expensive and time-consuming. Machine learning algorithms can be trained on the available data and make predictions for the cases where no numerical data are available. However, some widely used machine learning surrogates, such as neural networks, require a relatively large training dataset to produce accurate predictions. Multi-fidelity models, which combine high-fidelity data from numerical simulations and less expensive lower fidelity data from resources such as simplified constitutive equations, can pave the way for more accurate predictions. Here, we focus on neural networks and the Gaussian processes with two levels of fidelity, i.e., high and low fidelity networks, to predict the steady-state rheological properties, and compare them to the single-fidelity network. High-fidelity data are obtained from direct numerical simulations based on an immersed boundary method to couple the fluid and solid motion. The low-fidelity data are produced by using constitutive equations. Multiple neural networks and the Gaussian process structures are used for the hyperparameter tuning purpose. Results indicate that with the best choice of hyperparameters, both the multi-fidelity Gaussian processes and neural networks are capable of making predictions with a high level of accuracy with neural networks demonstrating marginally better performance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Shear thinning rheology in a concentrated suspension of fibers: the role of attractive interactions

Fiber-reinforced composites are ubiquitously encountered in the engineering, automobile, and aerospace industries. Fabrication of these composites requires the mixing of fibers dispersed in a liquid, and the final structure is affected by fiber properties, interactions, and flow fields. A good understanding of the rheology of fiber suspensions can aid in the design and optimization of the fabrication processes. To this end, we simulate the simple shear flow of fiber suspensions by accounting for short range lubrication, van der Waals attractive, electrostatic repulsion interactions, as well as friction between fibers. Direct numerical simulations are performed using the Immersed Boundary Method where the fibers are modeled as flexible slender bodies governed by the Euler-Bernoulli beam theory and the fluid flow is resolved using the Navier-Stokes equations. The simulation results for the suspension viscosity and yield stress are consistent with the experimental data from the literature for polyamide fiber suspensions. The shear-thinning behavior becomes stronger as we increase the magnitude of attractive interactions. Furthermore, we perform a parametric study varying fiber flexibility, volume fraction, aspect ratio, inertia and examine suspension viscosity and yield stress.

Khan, Monsurul↗