skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Macro–Micro-Coupled Simulations of Dilute Viscoelastic Fluids
Modeling the flow of polymer solutions requires knowledge at various length and time scales. The macroscopic behavior is described by the overall velocity, pressure, and stress. The polymeric contribution to the stress requires knowledge of the evolution of polymer chains. In this work, we use a microscopic model, the finitely extensible nonlinear elastic (FENE) model, to capture the polymer’s behavior. The benefit of using microscopic models is that they remain faithful to the polymer dynamics without information loss via averaging. Their downside is the computational cost incurred in solving the thousands to millions of differential equations describing the microstructure. Here, we describe a multiscale flow solver that utilizes GPUs for massively parallel, efficient simulations. We compare and contrast the microscopic model with its macroscopic counterpart under various flow conditions. In particular, significant differences are observed under nonlinear flow conditions, where the polymers become highly stretched and oriented.  more » « less
Award ID(s):
1751339
PAR ID:
10520139
Author(s) / Creator(s):
;
Publisher / Repository:
MPDI
Date Published:
Journal Name:
Applied Sciences
Volume:
13
Issue:
22
ISSN:
2076-3417
Page Range / eLocation ID:
12265
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Poole, Robert J (Ed.)
    The FENE-P (Finitely-Extensible Nonlinear Elastic) dumbbell constitutive equation is widely used in simulations and stability analyses of free and wall-bounded viscoelastic shear flows due to its relative simplicity and accuracy in predicting macroscopic properties of dilute polymer solutions. The model contains three independent material parameters, which expressed in dimensionless form correspond to a Weissenberg number (Wi), i.e., the ratio of the dumbbell relaxation time scale to a characteristic flow time scale, a finite extensibility parameter (L), corresponding to the ratio of the fully extended dumbbell length to the root mean square end-to-end separation of the polymer chain under equilibrium conditions, and a solvent viscosity ratio. An exact solution for the rheological predictions of the FENE-P model in steady simple shear flow is available [Sureshkumar et al., Phys Fluids (1997)], but the resulting nonlinear and nested set of equations do not readily reveal the key shear-thinning physics that dominates at high Wi as a result of the finite extensibility of the polymer chain. In this note we review a simple way of evaluating the steady material functions characterizing the nonlinear evolution of the polymeric contributions to the shear stress and first normal stress difference as the shear rate increases, provide asymptotic expansions as a function of Wi , and show that it is in fact possible to construct universal master curves for these two material functions as well as the corresponding stress ratio. Steady shear flow experiments on three highly elastic dilute polymer solutions of different finite extensibilities also follow the identified master curves. The governing dimensionless parameter for these master curves is Wi/L and it is only in strong shear flows exceeding Wi/L > 1 that the effects of finite extensibility of the polymer chains dominate the evolution of polymeric stresses in the flow field. We suggest that reporting the magnitude of Wi/L when performing stability analyses or simulating shear-dominated flows with the FENE-P model will help clarify finite extensibility effects. 
    more » « less
  2. Arrested soft materials such as gels and glasses exhibit a slow stress relaxation with a broad distribution of relaxation times in response to linear mechanical perturbations. Although this macroscopic stress relaxation is an essential feature in the application of arrested systems as structural materials, consumer products, foods, and biological materials, the microscopic origins of this relaxation remain poorly understood. Here, we elucidate the microscopic dynamics underlying the stress relaxation of such arrested soft materials under both quiescent and mechanically perturbed conditions through X-ray photon correlation spectroscopy. By studying the dynamics of a model associative gel system that undergoes dynamical arrest in the absence of aging effects, we show that the mean stress relaxation time measured from linear rheometry is directly correlated to the quiescent superdiffusive dynamics of the microscopic clusters, which are governed by a buildup of internal stresses during arrest. We also show that perturbing the system via small mechanical deformations can result in large intermittent fluctuations in the form of avalanches, which give rise to a broad non-Gaussian spectrum of relaxation modes at short times that is observed in stress relaxation measurements. These findings suggest that the linear viscoelastic stress relaxation in arrested soft materials may be governed by nonlinear phenomena involving an interplay of internal stress relaxations and perturbation-induced intermittent avalanches. 
    more » « less
  3. In this article, we introduce a new method for discretizing micro-macro models of dilute polymeric fluids by integrating a finite element method (FEM) discretization for the macroscopic fluid dynamics equation with a deterministic variational particle scheme for the microscopic Fokker-Planck equation. To address challenges arising from micro-macro coupling, we employ a discrete energetic variational approach to derive a coarse-grained micro-macro model with a particle approximation first and then develop a particle-FEM discretization for the coarse-grained model. The accuracy of the proposed method is evaluated for a Hookean dumbbell model in a Couette flow by comparing the computed velocity field with existing analytical solutions. We also use our method to study nonlinear FENE dumbbell models in different scenarios, such as extensional flow, pure shear flow, and lid-driven cavity flow. Numerical examples demonstrate that the proposed deterministic particle approach can accurately capture the various key rheological phenomena in the original FENE model, including hysteresis and δ-function-like spike behavior in extensional flows, velocity overshoot phenomenon in pure shear flows, symmetries breaking, vortex center shifting, and vortices weakening in lid-driven cavity flows, with a small number of particles. 
    more » « less
  4. The adaptive mechanical properties of soft and fibrous biological materials are relevant to their functionality. The emergence of the macroscopic response of these materials to external stress and intrinsic cell traction from local deformations of their structural components is not well understood. Here, we investigate the nonlinear elastic behavior of blood clots by combining microscopy, rheology, and an elastic network model that incorporates the stretching, bending, and buckling of constituent fibrin fibers. By inhibiting fibrin cross-linking in blood clots, we observe an anomalous softening regime in the macroscopic shear response as well as a reduction in platelet-induced clot contractility. Our model explains these observations from two independent macroscopic measurements in a unified manner, through a single mechanical parameter, the bending stiffness of individual fibers. Supported by experimental evidence, our mechanics-based model provides a framework for predicting and comprehending the nonlinear elastic behavior of blood clots and other active biopolymer networks in general. 
    more » « less
  5. The phenomenon of the ‘coil-stretch’ (C-S) transition, wherein a long-chain polymer initially in a coiled state undergoes a sudden configuration change to become fully stretched under steady elongational flows, has been widely recognized. This transition can display intricate hysteresis behaviours under specific critical conditions, giving rise to unique rheological characteristics in dilute polymer solutions. Historically, microscopic stochastic models and Brownian dynamics simulations have shed light on the underlying mechanisms of the transition by uncovering bistable configurations of polymer chains. Following the initial work by Cerf (J. Chem. Phys., vol. 20, 1952, pp. 395–402), we introduce a continuum model in this study to investigate the C-S transition in a constant uniaxial elongational flow. Our approach involves approximating the unfolding process of the polymer chain as an axisymmetric deformation of an elastic particle. We make the assumption that the particle possesses uniform material properties, which can be represented by a nonlinear, strain-hardening constitutive equation to replicate the finite extensibility of the polymer chain. Subsequently, we analytically solve for the steady-state deformation using a polarization method. By employing this reduced model, we effectively capture the C-S transition and establish its specific correlations with material and geometric properties. The hysteresis phenomena can be comprehended through a force-balance analysis, which involves comparing the externally applied viscous forces with the intrinsic elastic responsive forces. We demonstrate that our model, while simple, unveils rich elastohydrodynamics of the C-S transition. 
    more » « less