skip to main content

Attention:

The NSF Public Access Repository (PAR) system and access will be unavailable from 11:00 PM ET on Friday, December 13 until 2:00 AM ET on Saturday, December 14 due to maintenance. We apologize for the inconvenience.


Title: Geometric decomposition of the conformation tensor in viscoelastic turbulence
This work introduces a mathematical approach to analysing the polymer dynamics in turbulent viscoelastic flows that uses a new geometric decomposition of the conformation tensor, along with associated scalar measures of the polymer fluctuations. The approach circumvents an inherent difficulty in traditional Reynolds decompositions of the conformation tensor: the fluctuating tensor fields are not positive definite and so do not retain the physical meaning of the tensor. The geometric decomposition of the conformation tensor yields both mean and fluctuating tensor fields that are positive definite. The fluctuating tensor in the present decomposition has a clear physical interpretation as a polymer deformation relative to the mean configuration. Scalar measures of this fluctuating conformation tensor are developed based on the non-Euclidean geometry of the set of positive definite tensors. Drag-reduced viscoelastic turbulent channel flow is then used an example case study. The conformation tensor field, obtained using direct numerical simulations, is analysed using the proposed framework.  more » « less
Award ID(s):
1652244
PAR ID:
10084482
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
842
ISSN:
0022-1120
Page Range / eLocation ID:
395 to 427
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We consider the problem of formulating perturbative expansions of the conformation tensor, which is a positive definite tensor representing polymer deformation in viscoelastic flows. The classical approach does not explicitly take into account that the perturbed tensor must remain positive definite – a fact that has important physical implications, e.g. extensions and compressions are represented similarly to within a negative sign, when physically the former are unbounded and the latter are bounded from below. Mathematically, the classical approach assumes that the underlying geometry is Euclidean, and this assumption is not satisfied by the manifold of positive definite tensors. We provide an alternative formulation that retains the conveniences of classical perturbation methods used for generating linear and weakly nonlinear expansions, but also provides a clear physical interpretation and a mathematical basis for analysis. The approach is based on treating a perturbation as a sequence of successively smaller deformations of the polymer. Each deformation is modelled explicitly using geodesics on the manifold of positive definite tensors. Using geodesics, and associated geodesic distances, is the natural way to model perturbations to positive definite tensors because it is consistent with the manifold geometry. Approximations of the geodesics can then be used to reduce the total deformation to a series expansion in the small perturbation limit. We illustrate our approach using direct numerical simulations of the nonlinear evolution of Tollmien–Schlichting waves. 
    more » « less
  2. Summary

    We propose and investigate an additive regression model for symmetric positive-definite matrix-valued responses and multiple scalar predictors. The model exploits the Abelian group structure inherited from either of the log-Cholesky and log-Euclidean frameworks for symmetric positive-definite matrices and naturally extends to general Abelian Lie groups. The proposed additive model is shown to connect to an additive model on a tangent space. This connection not only entails an efficient algorithm to estimate the component functions, but also allows one to generalize the proposed additive model to general Riemannian manifolds. Optimal asymptotic convergence rates and normality of the estimated component functions are established, and numerical studies show that the proposed model enjoys good numerical performance, and is not subject to the curse of dimensionality when there are multiple predictors. The practical merits of the proposed model are demonstrated through an analysis of brain diffusion tensor imaging data.

     
    more » « less
  3. ABSTRACT

    Transport coefficients in turbulence are comprised of correlation functions between turbulent fluctuations and efficient methods to calculate them are desirable. For example, in mean-field dynamo theories used to model the growth of large-scale magnetic fields of stars and galaxies, the turbulent electromotive force is commonly approximated by a series of tensor products of turbulent transport coefficients with successively higher order spatial derivatives of the mean magnetic field. One ingredient of standard models is the kinematic coefficient of the zeroth-order term, namely the averaged kinetic pseudo-tensor $\boldsymbol \alpha$, that converts toroidal to poloidal fields. Here we demonstrate an efficient way to calculate this quantity for rotating stratified turbulence, whereby the pre-averaged quantity is calculated for the motion of a single plume, and the average is then taken over an ensemble of plumes of different orientations. We calculate the plume dynamics in the most convenient frame, before transforming back to the lab frame and averaging. Our concise configuration space calculation gives essentially identical results to previous lengthier approaches. The present application exemplifies what is a broadly applicable method.

     
    more » « less
  4. Abstract

    Recent developments of burst-mode lasers and imaging systems have opened new realms of simultaneous diagnostics for velocity and density fields at a rate of 1 kHz–1 MHz. These enable the exploration of previously unimaginable shock-driven turbulent flow fields that are of significant importance to problems in high-energy density physics. The current work presents novel measurements using simultaneous measurements of velocity and scalar fields at 60 kHz to investigate Richtmyer-Meshkov instability (RMI) in a spatio-temporal approach. The evolution of scalar fields and the vorticity dynamics responsible for the same are shown, including the interaction of shock with the interface. This temporal information is used to validate two vorticity-deposition models commonly used for initiation of large scale simulations, and have been previously validated only via simulations or integral measures of circulation. Additionally, these measurements also enable tracking the evolution and mode merging of individual flow structures that were previously not possible owing to inherently random variations in the interface at the smallest scales. A temporal evolution of symmetric vortex merging and the induced mixing prevalent in these problems is presented, with implications for the vortex paradigms in accelerated inhomogenous flows.

     
    more » « less
  5. From flocks of birds to biomolecular assemblies, systems in which many individual components independently consume energy to perform mechanical work exhibit a wide array of striking behaviors. Methods to quantify the dynamics of these so-called active systems generally aim to extract important length or time scales from experimental fields. Because such methods focus on extracting scalar values, they do not wring maximal information from experimental data. We introduce a method to overcome these limitations. We extend the framework of correlation functions by taking into account the internal headings of displacement fields. The functions we construct represent the material response to specific types of active perturbation within the system. Utilizing these response functions we query the material response of disparate active systems composed of actin filaments and myosin motors, from model fluids to living cells. We show we can extract critical length scales from the turbulent flows of an active nematic, anticipate contractility in an active gel, distinguish viscous from viscoelastic dissipation, and even differentiate modes of contractility in living cells. These examples underscore the vast utility of this method which measures response functions from experimental observations of complex active systems.

     
    more » « less