1 Introduction
In the present study, we address the following questions that arise in the context of viscoelastic turbulent flows. (i) Given a turbulent flow whose dynamics is partially governed by state variables that are positive definite tensors representing material deformation, what is an appropriate method to decompose the flow into a mean, or nominal, component and a deviation about that mean that preserves the physical character of the state variables? (ii) Are there corresponding scalar measures of the turbulence associated with these positive definite state variables? The conformation tensor is the relevant positive definite state variable in viscoelastic turbulence.
Dilute polymer solutions, viscoelastic flows obtained by adding small amounts of polymers to an incompressible Newtonian solvent, are the focus of the present work. The added polymers impart elasticity to the solvent which then causes the fluid to react not only to the deformation rate but also to the deformation history. As a result, a complete physical description of viscoelastic turbulence requires characterization of both the velocity,
$\boldsymbol{u}$
, and the conformation tensor,
$\unicode[STIX]{x1D63E}$
, which together form the state variables. The conformation tensor is a second-order positive definite tensor that encapsulates the polymer deformation history and is obtained by averaging, over molecular realizations, the dyad formed by the polymer end-to-end vector (Bird, Armstrong & Hassager Reference Bird, Armstrong and Hassager1987).
The conformation tensor affects the velocity field through the polymer stress,
$\unicode[STIX]{x1D64F}=\unicode[STIX]{x1D64F}(\unicode[STIX]{x1D63E})$
, while gradients in the velocity field are responsible for polymer stretching. Characterizing the mean polymer stress, or stress deficit, was the focus of early work because it was found to be necessary for closing the mean momentum balance (Willmarth, Wei & Lee Reference Willmarth, Wei and Lee1987). A non-vanishing stress deficit suggests the possibility of maintaining a turbulent velocity profile in the absence of Reynolds stresses, in which case the polymer dynamics would sustain turbulence. The experiments of Warholic, Massah & Hanratty (Reference Warholic, Massah and Hanratty1999) showed that a turbulent mean profile can, indeed, be maintained in the near absence of Reynolds stresses in channel flow. However, the polymer deformation itself is not readily accessible experimentally. Therefore, much of the work in understanding the mechanisms that lead to behaviour such as that found by Warholic et al. (Reference Warholic, Massah and Hanratty1999) has resorted to analytical treatments or direct numerical simulations (DNS), which we will briefly review below.
1.1 Previous approaches to quantifying polymer deformation and its effects
The main approach to analysing the polymer dynamics has been to utilize the statistics of the polymer forces and torques, or the normal stresses. The polymer force is the divergence of the polymer stress and the polymer torque is the curl of the polymer force. For example, de Angelis, Casciola & Piva (Reference de Angelis, Casciola and Piva2002) and Dubief et al. (Reference Dubief, Terrapon, White, Shaqfeh, Moin and Lele2005) showed that the cross-stream polymer force in turbulent channel flow counteracts spanwise variations in the velocity while enhancing streamwise advection, consistent with drag-reducing behaviour. The polymer torque acts in lockstep with the polymer force and counteracts streamwise vortices. It also inhibits generation of the heads of hairpin vortices (Kim et al. Reference Kim, Li, Sureshkumar, Balachandar and Adrian2007; Kim & Sureshkumar Reference Kim and Sureshkumar2013). Recent theoretical work proposed a vorticity–polymer torque formulation of the linearized governing equations and used it to reveal a reverse Orr mechanism for turbulence production in viscoelastic parallel shear flows (Page & Zaki Reference Page and Zaki2014, Reference Page and Zaki2015). Min, Choi & Yoo (Reference Min, Choi and Yoo2003a ), Min et al. (Reference Min, Yoo, Choi and Joseph2003b ), on the other hand, studied viscoelastic turbulent channel flow but used the elastic energy, defined there as proportional to the sum of the normal polymer stresses, to posit a theory of drag reduction that relied on an active exchange of elastic and kinetic energies in the flow.
A more appropriate quantity to probe the polymer deformation itself, and one that is also a state variable, is the conformation tensor
$\unicode[STIX]{x1D63E}$
. The trace of
$\unicode[STIX]{x1D63E}$
, denoted here as
$\text{tr}\,\unicode[STIX]{x1D63E}$
, is commonly used in the literature to analyse
$\unicode[STIX]{x1D63E}$
since it is equal to the sum of its principal stretches and is therefore a measure of the polymer deformation. For example, Sureshkumar, Beris & Handler (Reference Sureshkumar, Beris and Handler1997) considered first-order statistics of
$\text{tr}\,\unicode[STIX]{x1D63E}$
in their pioneering paper on the DNS of viscoelastic turbulent channel flows. The quantity
$\text{tr}\,\unicode[STIX]{x1D63E}$
is frequently used because it is proportional to the elastic energy in purely Hookean constitutive models of the polymers (Beris & Edwards Reference Beris and Edwards1994; Min et al.
Reference Min, Yoo, Choi and Joseph2003b
). However, it is often not a sufficiently complete descriptor of the polymer deformation; even if
$\text{tr}\,\unicode[STIX]{x1D63E}$
is held constant, the polymer may undergo a volumetric deformation.
Housiadas & Beris (Reference Housiadas and Beris2003) evaluated
$\text{tr}\,\unicode[STIX]{x1D63E}$
for a wide range of flow parameters in a viscoelastic turbulent channel flow and found a surprising result for certain parameter ranges: the mean of
$\text{tr}\,\unicode[STIX]{x1D63E}$
can increase with increasing elasticity without a commensurate effect on the mean velocity profile. A similar trend was also reported by Xi & Graham (Reference Xi and Graham2010) in minimal flow unit simulations. This trend is not inconsistent with the mean momentum balance as the mean of
$\text{tr}\,\unicode[STIX]{x1D63E}$
in turbulent channel flows can increase without effecting the mean velocity profile. This behaviour arises because the (mean) stress deficit is not a function of any of the normal components of
$\unicode[STIX]{x1D63E}$
. The normal components of
$\unicode[STIX]{x1D63E}$
do affect the mean momentum balance through the dynamical coupling between the different components of
$\unicode[STIX]{x1D63E}$
, but this relationship cannot be captured by
$\text{tr}\,\unicode[STIX]{x1D63E}$
. The situation described above highlights the importance of simultaneously considering all of the components of
$\unicode[STIX]{x1D63E}$
in order to arrive at a complete picture of the polymer deformation and its effect on the velocity field. Additionally, mean quantities such as those used by Housiadas & Beris (Reference Housiadas and Beris2003) and others are in themselves insufficient descriptors of the fluctuating polymer deformation; higher-order statistical quantities associated with
$\unicode[STIX]{x1D63E}$
are required to describe the fluctuations and their deviation away from the mean.
The fluctuating conformation tensor,
$\unicode[STIX]{x1D63E}^{\prime }$
, and its various moments provide one method to obtain pertinent higher-order statistical descriptions of the full conformation tensor,
$\unicode[STIX]{x1D63E}$
(see also Lee & Zaki Reference Lee and Zaki2017). The tensor
$\unicode[STIX]{x1D63E}^{\prime }$
is obtained by subtracting the mean conformation tensor from
$\unicode[STIX]{x1D63E}$
, in analogy with the Reynolds decomposition of
$\boldsymbol{u}$
. However, this fluctuating tensor is not guaranteed to be physically realizable; at least one realization of
$\unicode[STIX]{x1D63E}^{\prime }$
implies negative material deformation since it is guaranteed to lose positive-definiteness (
$\text{tr}\,\unicode[STIX]{x1D63E}^{\prime }$
must be
${\leqslant}0$
for at least one sample pulled from a statistical ensemble). Furthermore, it is not clear which scalar functions of
$\unicode[STIX]{x1D63E}^{\prime }$
provide mathematically consistent measures of the turbulence intensity associated with
$\unicode[STIX]{x1D63E}$
.
Although
$\unicode[STIX]{x1D63E}^{\prime }$
is not a physical conformation tensor, it can still be used for modelling Reynolds-averaged quantities in turbulence that do not necessarily require physically realizable fluctuating quantities. Indeed, it has been used with varying degrees of success in recent work to develop turbulence models (Li et al.
Reference Li, Gupta, Sureshkumar and Khomami2006; Iaccarino, Shaqfeh & Dubief Reference Iaccarino, Shaqfeh and Dubief2010; Resende et al.
Reference Resende, Kim, Younis, Sureshkumar and Pinho2011; Masoudian et al.
Reference Masoudian, Kim, Pinho and Sureshkumar2013) and to quantify subgrid stress contributions in large-eddy simulations (LES) (Masoudian, da Silva & Pinho Reference Masoudian, da Silva and Pinho2016). However, characterizing the polymer fluctuations using physically meaningful quantities is advantageous in that physical interpretations aid in modelling and provide a greater understanding of mechanism.
A physically motivated description of the velocity and conformation tensor field can be obtained using Karhunen–Loève or proper orthogonal decomposition (POD), as recently shown by Wang et al. (Reference Wang, Graham, Hahn and Xi2014). A POD is a global decomposition of a field quantity that yields an orthonormal basis that is optimally ordered in the sense of the best representation of the Euclidean norm (see Lumley Reference Lumley1970, for more details). For the velocity field this norm is the square-root of the kinetic energy but in a straightforward POD of the conformation tensor field the norm is not directly related to the elastic energy. Wang et al. (Reference Wang, Graham, Hahn and Xi2014) showed that a POD of the square-root of
$\unicode[STIX]{x1D63E}$
instead ensures best representation in terms of the elastic energy. Their approach crucially assumed that
$\text{tr}\,\unicode[STIX]{x1D63E}$
is proportional to the elastic energy. When this assumption is satisfied, their approach provides a valuable tool to extract the spatial structure of the dominant energetic components of viscoelastic turbulence. One can also use the individual modes to construct positive definite tensors that represent modal polymer deformation, since the POD basis is orthogonal with respect to a Frobenius inner product integrated over the spatial domain. However, these tensors cannot be used to construct a local decomposition of the conformation tensor into mean and fluctuating components because the sum of squares is not equal to the square of the sum: the cross-contributions of the tensors only vanish when we take the trace and integrate over the spatial domain. Reynolds decomposing the square-root of the conformation tensor is a local approach in which the cross-contribution similarly does not vanish instantaneously.
In the following subsection, we outline a local decomposition the conformation tensor into mean and fluctuating components, which overcomes the previously described difficulties with earlier approaches and which is one of the contributions of the present work. Another contribution is the development of scalar measures of the fluctuating component that depend on the Riemannian geometric structure of the set of positive definite tensors. Although the results we invoke for the latter contribution are well established and broadly applicable, they have not yet been widely used in fluid mechanics.
1.2 The present approach: motivation and summary of the framework
In the present work, we develop a formalism that allows us to evaluate the instantaneous deviation of the polymer deformation away from mean which respects the mathematical structure, and physical interpretation, of
$\unicode[STIX]{x1D63E}$
. Such a deviation yields an associated conformation tensor that can be used in analysis and modelling, e.g. the approach of Wang et al. (Reference Wang, Graham, Hahn and Xi2014) can be adapted to the analysis of this new conformation tensor. We also develop scalar measures of the turbulence intensity in the polymer deformation that reflect the distance, in the mathematically precise sense of a distance metric on a manifold, of the instantaneous deviation away from the mean.
In order to motivate our proposed formalism, we consider the following analogue encountered in the study of the stretching of material lines in turbulence (Batchelor Reference Batchelor1952). In this case, the normalized squared length of a material line,
$\ell ^{2}(t)>0$
, serves as the scalar analogue of
$\unicode[STIX]{x1D63E}$
. Since material lines cannot vanish,
$\ell ^{2}\neq 0$
. Let us assume that a statistically stationary state is possible and that
$\langle \ell ^{2}\rangle$
is then the expected value of the squared length of a material line. A Reynolds decomposition of
$\ell ^{2}=\langle \ell ^{2}\rangle +(\ell ^{2})^{\prime }$
yields a fluctuation
$(\ell ^{2})^{\prime }(t)$
that is not always positive, which implies a negative normalized squared length. One may, for the sake of argument, side-step the physical ambiguity implied by the negative squared length by considering only
$|(\ell ^{2})^{\prime }|$
but this does not solve the problem of asymmetry of
$|(\ell ^{2})^{\prime }|$
with respect to the direction of the stretching: when
$\ell ^{2}/\langle \ell ^{2}\rangle \in (1,\infty )$
, the material line is expanded with respect to
$\langle \ell ^{2}\rangle$
and when
$\ell ^{2}/\langle \ell ^{2}\rangle \in (0,1)$
it is compressed, which means that similarly probable states (expansion and contraction) would be described by fluctuations with very different magnitudes. A meaningful way to study the fluctuations in
$\ell ^{2}$
is by instead considering
$\log (\ell ^{2}/\langle \ell ^{2}\rangle )$
. Our goal is to generalize this latter type of construction to the conformation tensor, where one must take into account the tensorial nature of
$\unicode[STIX]{x1D63E}$
which encodes directional information not included in a scalar such as
$\ell ^{2}$
.
Following the scalar case described above, one approach to evaluating fluctuations in
$\unicode[STIX]{x1D63E}$
is to use
$\log \unicode[STIX]{x1D63E}$
, in lieu of
$\unicode[STIX]{x1D63E}$
, where
$\log$
here refers to the matrix logarithm. This approach is appealing because the logarithm of a positive definite matrix is a symmetric matrix and the set of symmetric matrices form a vector space, which therefore allows for a Reynolds decomposition analogous to that of
$\boldsymbol{u}$
, i.e.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn1.gif?pub-status=live)
where
$\langle \log \unicode[STIX]{x1D63E}\rangle$
is the expected value of
$\log \unicode[STIX]{x1D63E}$
. To the authors’ best knowledge, such a decomposition has not been previously used to characterize fluctuations in
$\unicode[STIX]{x1D63E}$
. However,
$\log \unicode[STIX]{x1D63E}$
itself has been an object of some interest in the viscoelastic literature. Fattal & Kupferman (Reference Fattal and Kupferman2004) introduced an approach for simulating viscoelastic flows that relied on evolving
$\log \unicode[STIX]{x1D63E}$
instead of
$\unicode[STIX]{x1D63E}$
. Fattal & Kupferman (Reference Fattal and Kupferman2004) and Hulsen, Fattal & Kupferman (Reference Hulsen, Fattal and Kupferman2005) then provided closed-form evolution equations for
$\log \unicode[STIX]{x1D63E}$
which explicitly depended on both
$\log \unicode[STIX]{x1D63E}$
as well as its spectral decomposition. Recently, Knechtges, Behr & Elgeti (Reference Knechtges, Behr and Elgeti2014) and Knechtges (Reference Knechtges2015) eliminated the explicit dependence on the spectral decomposition but at the expense of either imposing restrictions on the spectral radius of
$\unicode[STIX]{x1D63E}$
or introducing Dunford–Taylor-type integrals into the equations.
At least two additional difficulties arise in using
$\log \unicode[STIX]{x1D63E}$
. The first is that the expected value of
$\unicode[STIX]{x1D63E}$
is not equal to
$\text{e}^{\langle \log \unicode[STIX]{x1D63E}\rangle }$
. This fact implies that evaluating the effect of the polymer stress on the mean momentum balance requires all statistical moments of
$\log \unicode[STIX]{x1D63E}$
, even when the polymer stress is a linear function of
$\unicode[STIX]{x1D63E}$
. The second difficulty is that, in general,
$\text{e}^{\langle \log \unicode[STIX]{x1D63E}\rangle +(\log \unicode[STIX]{x1D63E})^{\prime }}\neq \text{e}^{\langle \log \unicode[STIX]{x1D63E}\rangle }\boldsymbol{\cdot }\text{e}^{(\log \unicode[STIX]{x1D63E})^{\prime }}$
which means that there is no way to associate
$(\log \unicode[STIX]{x1D63E})^{\prime }$
with a conformation tensor or a physical polymer deformation. It also means that there is no clear way to separate the effect of
$\langle \log \unicode[STIX]{x1D63E}\rangle$
in the fluctuating momentum balance.
In this paper, we derive a new conformation tensor,
$\unicode[STIX]{x1D642}$
, from a physical decomposition of the polymer deformation. Instead of the traditional additive decomposition
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn2.gif?pub-status=live)
the proposed geometric decomposition of
$\unicode[STIX]{x1D63E}$
is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn3.gif?pub-status=live)
where
$\overline{\unicode[STIX]{x1D641}}$
is a deformation gradient tensor that can be calculated directly from a base-flow conformation tensor,
$\overline{\unicode[STIX]{x1D63E}}$
, such that
$\overline{\unicode[STIX]{x1D641}}\boldsymbol{\cdot }\overline{\unicode[STIX]{x1D641}}^{\mathsf{T}}=\overline{\unicode[STIX]{x1D63E}}$
(this choice will be justified in § 3 below). This conformation tensor
$\unicode[STIX]{x1D642}$
is analogous to the scalar fluctuating quantity,
$\ell ^{2}/\langle \ell ^{2}\rangle$
.
The tensor
$\unicode[STIX]{x1D642}$
represents turbulent deviations from the mean conformation tensor and can be analysed by resorting to the curved, Riemannian geometry of the manifold of positive definite tensors. Interestingly, the first two moment invariants of the tensorial equivalent of
$\log (\ell ^{2}/\langle \ell ^{2}\rangle )$
, i.e.
$\log \unicode[STIX]{x1D642}$
, then appear as the relevant scalar measures for the fluctuations in
$\unicode[STIX]{x1D642}$
. The first moment invariant is the logarithm of the ratio of the volume of
$\unicode[STIX]{x1D63E}$
to the volume of
$\overline{\unicode[STIX]{x1D63E}}$
, where the volume of the conformation tensor refers to its determinant. The latter is proportional to the squared volume of the ellipsoid representing the coarse-grained polymer (Truesdell & Noll Reference Truesdell and Noll2004). The determinant also corresponds to the sphericity or conformational probability of the molecular structure (Beris & Edwards Reference Beris and Edwards1994). The second moment invariant is the metric distance of
$\unicode[STIX]{x1D642}$
away from
$\unicode[STIX]{x1D644}$
on the manifold of second-order positive definite tensors. Finally, we also propose a measure of the anisotropy of
$\unicode[STIX]{x1D63E}$
relative to the mean, based on the work of Moakher & Batchelor (Reference Moakher and Batchelor2006). This measure is equal to the metric distance of
$\unicode[STIX]{x1D642}$
to the closest isotropic tensor.
Finally, we use the proposed framework and direct numerical simulations to gain insight into the dynamics of viscoelastic (FENE-P) turbulent channel flow. Such flows are known to exhibit greatly reduced drag relative to an equivalent Newtonian flow, up to 60 % or more reduction in some cases (Toms Reference Toms1948; White & Mungal Reference White and Mungal2008). For such turbulent flows, separating the mean and fluctuating components of the conformation tensor in a physically consistent manner is an important step towards developing a quantitative understanding of the dynamics and isolating the relevant mechanisms at play.
The rest of the paper is organized as follows. The governing equations used in viscoelastic flows are reviewed in § 2. Section 3 presents the geometric decomposition of the conformation tensor, the associated evolution equations and the relation between the decomposition and the elastic energy. A review of a geometry constructed specifically for the set of positive definite tensors along with scalar measures of the polymer deformation based on this geometry and associated evolution equations are presented in § 4. In § 5, we present an example case study of viscoelastic turbulent channel flow to illustrate the concepts developed in the paper.
2 Governing equations
The non-dimensional governing equations for the velocity,
$\boldsymbol{u}$
, and conformation tensor,
$\unicode[STIX]{x1D63E}$
, in a viscoelastic flow are
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn4.gif?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn5.gif?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn6.gif?pub-status=live)
where
$\text{D}(\cdot )/\text{D}t=\unicode[STIX]{x2202}_{t}(\cdot )+u_{k}\unicode[STIX]{x2202}_{k}(\cdot )$
is the convective derivative,
$\text{sym}(\unicode[STIX]{x1D63C})=(1/2)(\unicode[STIX]{x1D63C}+\unicode[STIX]{x1D63C}^{\mathsf{T}})$
is the symmetric part of a second-order tensor
$\unicode[STIX]{x1D63C}$
,
$p$
is the pressure,
$Re$
is the Reynolds number,
$\unicode[STIX]{x1D6FD}\in [0,1]$
is the viscosity ratio and the polymer stress,
$\unicode[STIX]{x1D64F}$
, is a function of the conformation tensor,
$\unicode[STIX]{x1D63E}$
. The left-hand side of (2.3) is equal to the upper-convected Maxwell derivative, or the Lie derivative with respect to
$\boldsymbol{u}$
, of
$\unicode[STIX]{x1D63E}$
. Although we restrict our focus to the upper-convected Maxwell derivative, it can be replaced in (2.3) with any other co-rotational derivative, or objective rate.
The functional form of
$\unicode[STIX]{x1D64F}(\unicode[STIX]{x1D63E})$
depends on the particular constitutive model and strain measure used. Although we will not use a particular model in the theoretical development that will follow, we note for completeness that in the absence of inherent directionality in the polymers,
$\unicode[STIX]{x1D64F}$
is an isotropic function of
$\unicode[STIX]{x1D63E}$
defined locally at each
$(\boldsymbol{x},t)$
. Therefore, by the representation theorem (Truesdell & Noll Reference Truesdell and Noll2004) we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn7.gif?pub-status=live)
where the three characteristic tensor invariants of
$\unicode[STIX]{x1D63E}$
are defined as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn8.gif?pub-status=live)
and the Weissenberg number
$Wi$
is the polymer relaxation time normalized by the convective time scale.
Table 1 lists the coefficient functions,
$\unicode[STIX]{x1D707}_{i}$
for
$i=1,2,3$
, for two polymer models that are popular in the viscoelastic turbulence literature. The parameter
$L_{\max }$
is the maximum polymer extensibility.
In the following (primarily in §§ 3.1, 3.2 and 4.3), angled brackets,
$\langle \cdot \rangle$
, denote Reynolds spatio-temporal filtering (see appendix A in Sagaut Reference Sagaut2006), i.e. for a variable
$\unicode[STIX]{x1D719}(\boldsymbol{x},t)$
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn9.gif?pub-status=live)
where
$\mathscr{G}$
is a filtering kernel that is normalized so that
$\langle 1\rangle =1$
, and is defined such that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn10.gif?pub-status=live)
for any two integrable functions
$f_{1}=f_{1}(\boldsymbol{x},t)$
,
$f_{2}=f_{2}(\boldsymbol{x},t)$
. The mean of a quantity
$\unicode[STIX]{x1D719}$
is then
$\langle \unicode[STIX]{x1D719}\rangle$
and the
$n$
th moment of
$\unicode[STIX]{x1D719}$
is
$\langle \unicode[STIX]{x1D719}^{n}\rangle$
. The properties, (2.7), further imply that
$\langle \mathscr{F}(\langle \unicode[STIX]{x1D719}\rangle )\rangle =\mathscr{F}(\langle \unicode[STIX]{x1D719}\rangle )$
for any analytic function,
$\mathscr{F}$
. While the example case study presented in § 5 uses traditional Reynolds time-averaging, we present definitions using the filtering formulation since the approach is also expected to be valid more generally.
We use an overline symbol within the present text to denote the nominal or base-flow quantity associated with the symbol, which may be distinct from the averaged or filtered quantity.
Table 1. Coefficients
$\unicode[STIX]{x1D707}_{i}$
for two common models of polymers. Note that only
$\unicode[STIX]{x1D707}_{1}$
in the FENE-P model depends on an invariant of
$\unicode[STIX]{x1D63E}$
.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_tab1.gif?pub-status=live)
3 Decomposition of the conformation tensor
In the following, we will denote the general linear group of degree
$n$
, i.e. the set of
$n\times n$
matrices with non-zero determinant, as
$\boldsymbol{G}\boldsymbol{L}_{n}$
. We define the structure-preserving group action of
$\boldsymbol{G}\boldsymbol{L}_{n}$
on a set
$\boldsymbol{W}_{n}\subseteq \mathbb{R}^{n\times n}$
as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn11.gif?pub-status=live)
where
$\unicode[STIX]{x1D63C}\in \boldsymbol{G}\boldsymbol{L}_{n}$
and
$\unicode[STIX]{x1D63D}\in \boldsymbol{W}_{n}$
and by definition, we require
$\boldsymbol{W}_{n}$
to be invariant under the action.
From the perspective of continuum mechanics,
$\unicode[STIX]{x1D63E}>0$
is the left Cauchy–Green tensor associated with the deformation of the polymers (Beris & Edwards Reference Beris and Edwards1994; Rajagopal & Srinivasa Reference Rajagopal and Srinivasa2000; Cioranescu, Girault & Rajagopal Reference Cioranescu, Girault and Rajagopal2016), i.e.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn12.gif?pub-status=live)
where
$\unicode[STIX]{x1D641}$
is the deformation gradient with respect to an equilibrium configuration, also known as the distortion tensor. If the spatial coordinates in the micro-structure are given by
$\boldsymbol{a}=\boldsymbol{a}(\boldsymbol{a}_{0},t)$
where
$\boldsymbol{a}_{0}$
are the material coordinates, then
$\unicode[STIX]{x1D641}=\unicode[STIX]{x1D735}_{\boldsymbol{a}_{0}}\boldsymbol{a}=\unicode[STIX]{x2202}\boldsymbol{a}/\unicode[STIX]{x2202}\boldsymbol{a}_{0}$
so that a material line
$\text{d}\boldsymbol{a}_{0}$
deforms to
$\text{d}\boldsymbol{a}=\unicode[STIX]{x1D641}\boldsymbol{\cdot }\text{d}\boldsymbol{a}_{0}$
under the deformation represented by
$\unicode[STIX]{x1D63E}$
. When
$\unicode[STIX]{x1D641}$
is restricted to be symmetric, (3.2) reduces to the factorization proposed by Balci et al. (Reference Balci, Thomases, Renardy and Doering2011) to improve numerical schemes for evolving the conformation tensor equations.
Let
$\overline{\unicode[STIX]{x1D63E}}$
be a nominal conformation tensor such as the mean or laminar base-flow conformation tensor. The only requirement we impose on
$\overline{\unicode[STIX]{x1D63E}}$
is that it must be defined according to a rule that ensures that
$\unicode[STIX]{x1D63E}$
and
$\overline{\unicode[STIX]{x1D63E}}$
cannot be arbitrarily rotated with respect to each other. In other words, if
$\unicode[STIX]{x1D63E}$
transforms to
$[\unicode[STIX]{x1D63E}]_{\unicode[STIX]{x1D64D}}$
then
$\overline{\unicode[STIX]{x1D63E}}$
must transform to
$[\overline{\unicode[STIX]{x1D63E}}]_{\unicode[STIX]{x1D64D}}$
for any
$\unicode[STIX]{x1D64D}\in \boldsymbol{S}\boldsymbol{O}_{3}$
, where
$\boldsymbol{S}\boldsymbol{O}_{n}$
denotes the
$n\times n$
special orthogonal group (or rotation matrices). Define
$\overline{\unicode[STIX]{x1D641}}\in \boldsymbol{G}\boldsymbol{L}_{3}$
with
$\text{det}\,\overline{\unicode[STIX]{x1D641}}>0$
as the tensor that satisfies
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn13.gif?pub-status=live)
Such an
$\overline{\unicode[STIX]{x1D641}}$
is non-unique as it can be parametrized as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn14.gif?pub-status=live)
for any
$\unicode[STIX]{x1D64D}\in \boldsymbol{S}\boldsymbol{O}_{3}$
and where
$\overline{\unicode[STIX]{x1D63E}}^{1/2}$
is the unique matrix square-root of
$\overline{\unicode[STIX]{x1D63E}}$
. Since the polar decomposition of
$\overline{\unicode[STIX]{x1D641}}$
and the square-root of
$\overline{\unicode[STIX]{x1D63E}}$
(up to a
$\pm$
sign change) are both unique, (3.4) is a parametrization of all possible
$\overline{\unicode[STIX]{x1D641}}$
. The tensor
$\overline{\unicode[STIX]{x1D641}}$
serves as a deformation gradient associated with the mean configuration.
The
$n$
th power of a positive definite tensor
$\unicode[STIX]{x1D63C}$
is a tensor with the same eigenvectors as
$\unicode[STIX]{x1D63C}$
and associated eigenvalues equal to the corresponding eigenvalues of
$\unicode[STIX]{x1D63C}$
raised to the
$n$
th power. In practice, since these
$n$
th powers are isotropic functions of
$\unicode[STIX]{x1D63C}$
, one need not explicitly perform a spectral decomposition to calculate them. For example, an application of the representation theorem can be used to express
$\unicode[STIX]{x1D63C}^{1/2}$
and
$\unicode[STIX]{x1D63C}^{-1/2}$
solely in terms of
$\unicode[STIX]{x1D63C}$
and its invariants (Hoger & Carlson Reference Hoger and Carlson1984; Ting Reference Ting1985).
Given a specific
$\overline{\unicode[STIX]{x1D641}}$
one satisfying (3.4), we then decompose the full distortion tensor
$\unicode[STIX]{x1D641}$
about
$\overline{\unicode[STIX]{x1D641}}$
by considering successive transformations on the material line
$\text{d}\boldsymbol{a}_{0}$
, i.e.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn15.gif?pub-status=live)
where
$\unicode[STIX]{x1D647}=\overline{\unicode[STIX]{x1D641}}^{\,-1}\boldsymbol{\cdot }\unicode[STIX]{x1D641}$
is the fluctuating distortion tensor. This decomposition is illustrated in figure 1.
Substituting
$\unicode[STIX]{x1D641}=\overline{\unicode[STIX]{x1D641}}\boldsymbol{\cdot }\unicode[STIX]{x1D647}$
in (3.2), we then arrive at a geometric decomposition of the conformation tensor
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn16.gif?pub-status=live)
where
$\unicode[STIX]{x1D642}=\unicode[STIX]{x1D647}\boldsymbol{\cdot }\unicode[STIX]{x1D647}^{\mathsf{T}}$
is a left Cauchy–Green tensor that is analogous to
$\unicode[STIX]{x1D63E}$
. Comparing (3.6) and (1.2), we can relate
$\unicode[STIX]{x1D63E}^{\prime }$
to
$\unicode[STIX]{x1D642}$
as follows:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn17.gif?pub-status=live)
From this point of view, the geometric decomposition provides a framework for interpreting the fluctuating tensor,
$\unicode[STIX]{x1D63E}^{\prime }$
, obtained from the Reynolds decomposition.
Although a specific
$\unicode[STIX]{x1D642}$
, and in particular only its set of principal axes, depends on
$\unicode[STIX]{x1D64D}\in \boldsymbol{S}\boldsymbol{O}_{3}$
chosen in (3.4), any function of only the invariants of
$\unicode[STIX]{x1D642}$
is independent of the choice of
$\unicode[STIX]{x1D64D}$
. This class of functions includes all objective scalar functions of
$\unicode[STIX]{x1D642}$
; indeed, the scalar characterizations of the fluctuations that we develop later are also independent of
$\unicode[STIX]{x1D64D}$
. With respect to the full tensor,
$\unicode[STIX]{x1D642}$
, we will later find that
$\unicode[STIX]{x1D64D}=\unicode[STIX]{x1D644}$
is a natural choice.
The decomposition of
$\unicode[STIX]{x1D641}$
into successive deformations, as in (3.5), is reminiscent of the multiplicative decomposition in large deformation theory that has found numerous applications over the last few decades (Casey Reference Casey2015; Sadik & Yavari Reference Sadik and Yavari2017). For example, in elasto-plasticity theory, the deformation gradient is decomposed into successive plastic and elastic deformations with the objective of formulating constitutive laws for each of the deformations somewhat independently. Similar constructions are used in thermo-elasticity and biomechanics (Lubarda Reference Lubarda2004). A full review of that literature is beyond the scope of the present work, but it suffices to note that the present case is greatly simplified because the constitutive laws are already specified and the focus is on the analysis of the polymer deformation due to turbulence.
We next present the equations for mean and fluctuating quantities in the geometric decomposition when the nominal conformation tensor is obtained by averaging or Reynolds filtering.
3.1 Evolution equations in the Reynolds-filtered case
In this section, we will consider the case when the nominal tensor is obtained using Reynolds filtering. We choose to restrict our attention to Reynolds filters but the development can be generalized to other filters, e.g. for applications in large-eddy simulations (LES). We thus have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn18.gif?pub-status=live)
By the properties (2.7), the associated
$\overline{\unicode[STIX]{x1D641}}$
satisfies
$\langle \overline{\unicode[STIX]{x1D641}}\rangle =\overline{\unicode[STIX]{x1D641}}$
. Applying the averaging operation to (3.6) yields
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn19.gif?pub-status=live)
where
$\unicode[STIX]{x1D64D}(\boldsymbol{x},t)$
is the rotation tensor field given in (3.4). Henceforth, we will restrict the rotation tensor field so that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn20.gif?pub-status=live)
By (2.7), we then have
$\langle \unicode[STIX]{x1D642}\rangle =\unicode[STIX]{x1D644}$
.
The Reynolds decomposition is applied to
$p$
and
$\boldsymbol{u}$
while
$\unicode[STIX]{x1D63E}$
is decomposed using (3.6) with
$\overline{\unicode[STIX]{x1D63E}}$
defined according to (3.8). We thus have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn21.gif?pub-status=live)
where
$\overline{\boldsymbol{u}}=\langle \boldsymbol{u}\rangle$
and
$\overline{p}=\langle p\rangle$
and the primes denote fluctuating quantities obtained via the Reynolds decomposition. In general,
$\overline{p}=\overline{p}(\boldsymbol{x},t)$
,
$\overline{\boldsymbol{u}}=\overline{\boldsymbol{u}}(\boldsymbol{x},t)$
,
$\overline{\unicode[STIX]{x1D63E}}=\overline{\unicode[STIX]{x1D63E}}(\boldsymbol{x},t)$
. Note that
$\overline{\unicode[STIX]{x1D641}}\neq \langle \unicode[STIX]{x1D641}\rangle$
, in general.
Following the standard procedure, we can then decompose the momentum equation as follows:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn22.gif?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn23.gif?pub-status=live)
where
$\overline{\unicode[STIX]{x1D64F}}=\langle \unicode[STIX]{x1D64F}\rangle$
,
$\unicode[STIX]{x1D64F}^{\prime }=\unicode[STIX]{x1D64F}-\overline{\unicode[STIX]{x1D64F}}$
,
$\overline{\boldsymbol{u}^{\prime }\boldsymbol{u}^{\prime }}=\langle \boldsymbol{u}^{\prime }\boldsymbol{u}^{\prime }\rangle$
and
$(\boldsymbol{u}^{\prime }\boldsymbol{u}^{\prime })^{\prime }=\boldsymbol{u}^{\prime }\boldsymbol{u}^{\prime }-\overline{\boldsymbol{u}^{\prime }\boldsymbol{u}^{\prime }}$
.
The precise form of
$\overline{\unicode[STIX]{x1D64F}}$
, which appears in the mean momentum equation in (3.13), depends on the constitutive model used. In the Oldroyd-B model,
$\overline{\unicode[STIX]{x1D64F}}$
only depends on
$\overline{\unicode[STIX]{x1D641}}$
:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn24.gif?pub-status=live)
In models that are nonlinear in
$\unicode[STIX]{x1D63E}$
, the fluctuating tensor
$\unicode[STIX]{x1D642}$
cannot be eliminated or factored out of
$\overline{\unicode[STIX]{x1D64F}}$
. For example, in the FENE-P model,
$\overline{\unicode[STIX]{x1D64F}}$
can be expressed as a series in which the dominant term is equal to (3.14) while the remaining terms depend on higher-order moments of
$\unicode[STIX]{x1D642}$
. In general, we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn25.gif?pub-status=live)
Substituting (3.11) into (2.3) and applying the filtering operation
$\langle \cdot \rangle$
defined in (2.6) yields the following equations for
$\overline{\unicode[STIX]{x1D63E}}$
:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn26.gif?pub-status=live)
Term (a) is the averaged turbulent transport and terms (b), (c) describe the mean stretching and rotation of the polymer arising due to the gradients in the fluctuating velocity field. The right-hand side of (3.16) is the cumulative effect of the turbulent fluctuations on the mean balance. The mean balance (3.16) can also be written as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn27.gif?pub-status=live)
where
$\overline{\unicode[STIX]{x1D640}}=\langle \unicode[STIX]{x1D640}\rangle$
,
$\unicode[STIX]{x1D640}^{\prime }=\unicode[STIX]{x1D640}-\overline{\unicode[STIX]{x1D640}}$
,
$\overline{\unicode[STIX]{x1D648}}=\langle \unicode[STIX]{x1D648}\rangle$
, and
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn28.gif?pub-status=live)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn29.gif?pub-status=live)
Here, the tensors
$\unicode[STIX]{x1D640}$
and
$\unicode[STIX]{x1D648}$
serve as modified velocity gradient and polymer stress tensors. Note that the invariants of
$\unicode[STIX]{x1D735}\boldsymbol{u}$
and
$\unicode[STIX]{x1D640}$
coincide. The equation (3.17) shows that the mean modified stress,
$\overline{\unicode[STIX]{x1D648}}$
, is a function of the mean velocity gradient,
$\overline{\unicode[STIX]{x1D640}}$
, the mean stretching due to turbulent velocity gradients
$\langle \unicode[STIX]{x1D642}\boldsymbol{\cdot }\unicode[STIX]{x1D640}^{\prime }\rangle$
, and the turbulent advection of the fluctuating polymer deformation,
$\langle \boldsymbol{u}^{\prime }\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D642}\rangle$
. Additionally, new terms appear due to the time-rate of change of
$\overline{\unicode[STIX]{x1D641}}$
, and modified advection of
$\overline{\unicode[STIX]{x1D641}}$
.
We can find the evolution equation for
$\unicode[STIX]{x1D642}$
by substituting (3.11) into (2.3) and simplifying, which yields
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn30.gif?pub-status=live)
where
$\unicode[STIX]{x1D646}$
is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn31.gif?pub-status=live)
and represents the modified velocity gradient augmented with an additional stretching that arises due to the decomposition. The expression (3.20) is more general than considered here; an equivalent expression can be derived when the nominal tensor,
$\overline{\unicode[STIX]{x1D63E}}$
, is not equal to the mean conformation tensor.
The higher-dimensional nature of
$\unicode[STIX]{x1D642}$
makes the quantification of the fluctuating turbulent polymer deformation a more difficult task. We will examine the elastic potential energy as a method to evaluate this deformation in the next subsection and then introduce more general scalar characterizations of
$\unicode[STIX]{x1D642}$
in the next section.
3.2 Elastic energy and its relation to
$\unicode[STIX]{x1D642}$
The turbulent mean polymer configuration is not the thermodynamic equilibrium state, and thus
$\unicode[STIX]{x1D642}$
alone is not sufficient to fully determine thermodynamic quantities such as the elastic potential energy,
$\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D63E})$
. For example, Beris & Edwards (Reference Beris and Edwards1994) define
$\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D63E})$
for an Oldroyd-B model as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn32.gif?pub-status=live)
where we have rewritten the expression in terms of
$\unicode[STIX]{x1D642}$
and
$\overline{\unicode[STIX]{x1D63E}}$
by setting
$\overline{\unicode[STIX]{x1D641}}=\overline{\unicode[STIX]{x1D63E}}^{1/2}$
and using the cyclic property of the trace to obtain
$\text{tr}\,\unicode[STIX]{x1D63E}=\text{tr}(\overline{\unicode[STIX]{x1D63E}}\boldsymbol{\cdot }\unicode[STIX]{x1D642})$
. Here,
$\unicode[STIX]{x1D6FA}$
is the spatial domain, and the scalar function
$\unicode[STIX]{x1D713}(\boldsymbol{x})$
is proportional to the polymer elastic constant times the elasticity density.
The mean elastic potential,
$\langle \unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D63E})\rangle$
, for the Oldroyd-B model has the convenient property that it can be written solely in terms of the mean conformation tensor:
$\langle \unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D63E})\rangle =\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\overline{\unicode[STIX]{x1D63E}})$
. However, the contribution of
$\unicode[STIX]{x1D642}$
in
$\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D63E})$
cannot be fully separated from that of
$\overline{\unicode[STIX]{x1D63E}}$
because
$\text{tr}\,\unicode[STIX]{x1D63C}\boldsymbol{\cdot }\unicode[STIX]{x1D63D}\neq \text{tr}\,\unicode[STIX]{x1D63C}\,\text{tr}\,\unicode[STIX]{x1D63D}$
. Nonetheless, insight into the role of the different contributions can be obtained by using a trace inequality proved by Mori (Reference Mori1988), which yields
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn33.gif?pub-status=live)
where
$\unicode[STIX]{x1D70E}_{i}(\unicode[STIX]{x1D63C})$
denotes the
$i$
th largest eigenvalue of a tensor
$\unicode[STIX]{x1D63C}$
. In terms of the bounds in (3.23), the contribution of
$\overline{\unicode[STIX]{x1D63E}}$
to
$\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D63E})$
is equivalent to a modification of the local elasticity density or the elastic constant.
In other constitutive models, the contribution to the elastic potential energy from the mean polymer deformation is more difficult to separate. For example, the elastic potential energy for the FENE-P model (Beris & Edwards Reference Beris and Edwards1994) is
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn34.gif?pub-status=live)
where
$L_{\max }$
,
$\unicode[STIX]{x1D6FA}$
and
$\unicode[STIX]{x1D713}$
are as defined before. Here, the mean elastic potential energy,
$\langle \unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D63E};L_{\max })\rangle$
, cannot be separated from
$\unicode[STIX]{x1D642}$
as in the Oldroyd-B model because, according to (3.24),
$\langle \unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D63E})\rangle \neq \unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\overline{\unicode[STIX]{x1D63E}})$
. However, we can again bound
$\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D63E};L_{\max })$
as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn35.gif?pub-status=live)
In terms of the bounds in (3.25), the contribution of
$\overline{\unicode[STIX]{x1D63E}}$
in
$\unicode[STIX]{x1D700}_{\unicode[STIX]{x1D713}}(\unicode[STIX]{x1D63E};L_{\max })$
is equivalent to a modification of the local polymer maximum extensibility.
Elastic energy may itself be insufficient to fully characterize the polymer deformation. For example, in both Oldroyd-B and FENE-P models, the elastic energy is equal for all conformation tensors that are given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn36.gif?pub-status=live)
where
$\text{diag}(\unicode[STIX]{x1D719}_{1},\unicode[STIX]{x1D719}_{2},\unicode[STIX]{x1D719}_{3})$
denotes a diagonal tensor with the
$i$
th diagonal component given by
$\unicode[STIX]{x1D719}_{i}$
, and
$\unicode[STIX]{x1D6FC}_{1}\geqslant \unicode[STIX]{x1D6FC}_{2}\geqslant \unicode[STIX]{x1D6FC}_{3}>0$
and
$\unicode[STIX]{x1D64C}\in \boldsymbol{S}\boldsymbol{O}_{3}$
are fixed. Even though the trace is fixed, the volume of the deformation ellipsoid changes with
$\unicode[STIX]{x1D6FF}$
, and is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn37.gif?pub-status=live)
In addition, since the governing equations are not Hamiltonian (Beris & Edwards Reference Beris and Edwards1994), the elastic potential energy only provides a partial characterization of the dynamics underlying the polymer deformation. Due to the above limitations of the elastic energy, and its dependence on the choice of the particular constitutive model, we instead develop an approach to characterizing the polymer deformation using the inherent geometric structure underlying
$\unicode[STIX]{x1D642}$
. This approach, introduced in the next section, is mathematically rigorous and can be applied to any positive definite tensor.
4 A Riemannian approach to the fluctuating polymer deformation
Any scalar characterization of
$\unicode[STIX]{x1D642}$
obeying the principle of objectivity can be a function only of its invariants,
$\text{I}_{\unicode[STIX]{x1D642}}$
,
$\text{II}_{\unicode[STIX]{x1D642}}$
and
$\text{III}_{\unicode[STIX]{x1D642}}$
. The invariants can be interpreted in terms of the fluctuating deformation ellipsoid, i.e. the ellipsoid associated with
$\unicode[STIX]{x1D642}$
. The first invariant,
$\text{I}_{\unicode[STIX]{x1D642}}$
, is proportional to the average radius of the ellipsoid, the second invariant,
$\text{II}_{\unicode[STIX]{x1D642}}$
, is proportional to a lower bound for the surface area (Klamkin Reference Klamkin1971) and the third invariant
$\text{III}_{\unicode[STIX]{x1D642}}$
is the volume of the deformation ellipsoid. Note that the eigenvalues (or principal stretches) of a conformation tensor are equal to the squared polymer stretches.
In practice, multiple difficulties arise in naively using the invariants of
$\unicode[STIX]{x1D642}$
to characterize the conformation tensor. For example, consider the isotropic case with
$\unicode[STIX]{x1D63E}=a\unicode[STIX]{x1D644}$
and
$\overline{\unicode[STIX]{x1D63E}}=b\unicode[STIX]{x1D644}$
. We then have
$\unicode[STIX]{x1D642}=(a/b)\unicode[STIX]{x1D644}$
and the three invariants reduce to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn38.gif?pub-status=live)
which implies that the invariants are bounded between 0 and 1 for compressions with respect to
$\overline{\unicode[STIX]{x1D63E}}$
and between 1 and
$+\infty$
for expansions with respect to
$\overline{\unicode[STIX]{x1D63E}}$
. This inherent asymmetry in the characterization is undesirable. The statistical moments of the invariants also vary over several orders of magnitude, rendering these moments uninformative predictors of the level of turbulent stretching in the polymers.
The problems discussed above arise because the set of
$n\times n$
positive definite matrices, denoted
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
, for
$n>0$
, does not form a vector space and thus the Euclidean notions of translation and shortest distances between points are not valid. For example, let
$\unicode[STIX]{x1D63C},\unicode[STIX]{x1D63D}\in \boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
, and define
$\unicode[STIX]{x1D653}$
as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn39.gif?pub-status=live)
One may wish to use the parameter
$r$
to denote ‘distance of
$\unicode[STIX]{x1D653}$
to
$\unicode[STIX]{x1D63C}$
’ along the ‘direction between
$\unicode[STIX]{x1D63C}$
and
$\unicode[STIX]{x1D63D}$
’. However,
$\unicode[STIX]{x1D653}$
is then guaranteed to be positive definite only if
$r\in [0,1]$
. While
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
is not a vector space, it has a Riemannian geometric structure that can exploited to formulate alternative scalar measures of
$\unicode[STIX]{x1D642}$
that do not suffer from the problems mentioned above. We introduce this geometry in § 4.1, including definitions of shortest paths and distances between tensors. Subsequently, in § 4.2, we introduce scalar measures based on the development in § 4.1 that can be used to quantify the turbulent fluctuations in the polymers. In § 4.3, we derive the Reynolds-filtered evolution equations for the scalar measures.
4.1 Geodesic curves and distances between positive definite tensors
The set
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
is a Cartan–Hadamard manifold: it is a simply connected, geodesically complete Riemannian manifold with seminegative curvature (Lang Reference Lang2001). We summarize this characterization in the present section in order to develop a notion of distances between positive definite tensors that will be used to formulate appropriate scalar measures of the fluctuating conformation tensor
$\unicode[STIX]{x1D642}$
. Details on the Riemannian structure of
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
and theorems leading to the results used in this section are presented in the appendix.
Consider two matrices
$\unicode[STIX]{x1D653},\unicode[STIX]{x1D654}\in \boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
. In the set of all curves along the manifold
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
connecting
$\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D654}$
, there exists a unique curve that minimizes the distance between
$\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D654}$
with respect to the Riemannian metric on
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
, i.e. there exists a
$\unicode[STIX]{x1D64B}(r)$
with
$\unicode[STIX]{x1D64B}(0)=\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D64B}(1)=\unicode[STIX]{x1D654}$
that uniquely minimizes the distance traversed along the manifold between
$\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D654}$
. We call this curve the geodesic curve along the manifold and it is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn40.gif?pub-status=live)
The geodesic distance associated with the geodesic curve between
$\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D654}$
is the minimum separation between them along the manifold and is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn41.gif?pub-status=live)
The distance,
$d(\unicode[STIX]{x1D653},\unicode[STIX]{x1D654})$
, is affine invariant, i.e.
$d(\unicode[STIX]{x1D653},\unicode[STIX]{x1D654})=d([\unicode[STIX]{x1D653}]_{\unicode[STIX]{x1D63C}},[\unicode[STIX]{x1D654}]_{\unicode[STIX]{x1D63C}})$
for all
$\unicode[STIX]{x1D63C}\in \boldsymbol{G}\boldsymbol{L}_{3}$
.
Geodesic curves and distances along a Riemannian manifold are analogous to straight lines and distances in Euclidean space. In the case of
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
, the analogy can be taken quite far because
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
is geodesically complete; a geodesic connecting any two points on the manifold parametrized by
$r$
can be arbitrarily extended letting
$r\in \mathbb{R}$
. For example,
$\unicode[STIX]{x1D653}\#_{r}\unicode[STIX]{x1D654}$
with
$r\in [0,a]$
, is a geodesic between
$\unicode[STIX]{x1D653}\#_{0}\unicode[STIX]{x1D654}$
and
$\unicode[STIX]{x1D653}\#_{a}\unicode[STIX]{x1D654}$
for all
$[0,a]\subseteq \mathbb{R}$
. Furthermore, for each
$r>0$
, we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn42.gif?pub-status=live)
We now illustrate the geodesic distance derived above using two specific examples.
-
(i) Isotropic tensors: Let
$\unicode[STIX]{x1D653}=a\unicode[STIX]{x1D644}$ and
$\unicode[STIX]{x1D654}=b\unicode[STIX]{x1D644}$ be elements of the one-dimensional sub-manifold of
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$ consisting of the isotropic tensors. The geodesic path joining
$\unicode[STIX]{x1D653}$ and
$\unicode[STIX]{x1D654}$ is given by
(4.6)and the geodesic distance is given by$$\begin{eqnarray}\displaystyle \unicode[STIX]{x1D653}\#_{r}\unicode[STIX]{x1D654}=(a^{1-r}b^{r})\unicode[STIX]{x1D644} & & \displaystyle\end{eqnarray}$$
$d(\unicode[STIX]{x1D653},\unicode[STIX]{x1D654})=\sqrt{3}\log (b/a)$ .
Notice that
$a^{1-r}b^{r}$ , which appears in (4.6), is a generalized geometric mean of
$a$ and
$b$ with the classical definition realized at
$r=1/2$ . It can be shown that a similar interpretation is admissible when
$\unicode[STIX]{x1D653}$ and
$\unicode[STIX]{x1D654}$ are not isotropic (Bhatia Reference Bhatia2015). This fact has formed the basis of efforts to formulate alternative definitions of statistical quantities such as means and covariances so that they conform to the geometric structure of
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$ (Pennec, Fillard & Ayache Reference Pennec, Fillard and Ayache2006; Fletcher & Joshi Reference Fletcher and Joshi2007).
-
(ii) Tensors differing by a rotation: Consider
$\unicode[STIX]{x1D653}$ and
$\unicode[STIX]{x1D654}=[\unicode[STIX]{x1D653}]_{\unicode[STIX]{x1D64D}}$ for
$\unicode[STIX]{x1D64D}\in \boldsymbol{S}\boldsymbol{O}_{3}$ . The geodesic joining
$\unicode[STIX]{x1D653}$ and
$\unicode[STIX]{x1D654}$ is given by
(4.7)The distance between$$\begin{eqnarray}\displaystyle \unicode[STIX]{x1D653}\#_{r}\unicode[STIX]{x1D654}=[([\unicode[STIX]{x1D653}]_{\unicode[STIX]{x1D653}^{-1/2}\boldsymbol{\cdot }\unicode[STIX]{x1D64D}})^{r}]_{\unicode[STIX]{x1D653}^{1/2}}. & & \displaystyle\end{eqnarray}$$
$\unicode[STIX]{x1D653}$ and
$\unicode[STIX]{x1D654}$ is then bounded as
(4.8)where the lower bound is achieved for$$\begin{eqnarray}\displaystyle 0\leqslant d(\unicode[STIX]{x1D653},\unicode[STIX]{x1D654})\leqslant \sqrt{3}\min \,\left\{\max _{i}\left\{\left|\log \left(\frac{\unicode[STIX]{x1D70E}_{i}(\unicode[STIX]{x1D654})}{\unicode[STIX]{x1D70E}_{3}(\unicode[STIX]{x1D653})}\right)\right|\right\},\max _{i}\left\{\left|\log \left(\frac{\unicode[STIX]{x1D70E}_{1}(\unicode[STIX]{x1D654})}{\unicode[STIX]{x1D70E}_{i}(\unicode[STIX]{x1D653})}\right)\right|\right\}\right\}, & & \displaystyle \nonumber\\ \displaystyle & & \displaystyle\end{eqnarray}$$
$\unicode[STIX]{x1D64D}=\unicode[STIX]{x1D644}$ . The upper bound in (4.8) suggests that a differential rotation of a second-order tensor,
$\unicode[STIX]{x1D653}$ , leads to an excursion along
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$ with a path length that depends on the anisotropy of
$\unicode[STIX]{x1D653}$ . For isotropic
$\unicode[STIX]{x1D653}$ , the path length is zero and it otherwise increases with increasing anisotropy.
The geometry of
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
and the properties discussed above are next used to define scalar measures that characterize the turbulent fluctuations in
$\unicode[STIX]{x1D642}$
.
4.2 Scalar measures of the fluctuating conformation tensor
In this subsection we introduce scalar measures that can be used to quantify the fluctuating polymer deformation represented by
$\unicode[STIX]{x1D642}$
. In what follows, we will denote the matrix logarithm of
$\unicode[STIX]{x1D642}$
as
$\boldsymbol{{\mathcal{G}}}$
, i.e.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn46.gif?pub-status=live)
The matrix logarithm is guaranteed to exist and is unique since
$\unicode[STIX]{x1D642}$
is positive definite. A key point to note is that the eigenvalues of
$\boldsymbol{{\mathcal{G}}}$
are the logarithms of the eigenvalues of
$\unicode[STIX]{x1D642}$
.
4.2.1 Logarithmic volume ratio,
$\unicode[STIX]{x1D701}$
Let
$\unicode[STIX]{x1D6E4}_{i}=\unicode[STIX]{x1D70E}_{i}(\unicode[STIX]{x1D642})$
, for
$i=1,2,3$
, be the eigenvalues of
$\unicode[STIX]{x1D642}$
. Then
$\log (\text{det}\,\unicode[STIX]{x1D642})=\log (\prod _{i=1}^{3}\unicode[STIX]{x1D6E4}_{i})=\sum _{i=1}^{3}\log \unicode[STIX]{x1D6E4}_{i}$
. We thus define the logarithmic volume ratio of the fluctuation,
$\unicode[STIX]{x1D701}$
, as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn47.gif?pub-status=live)
When
$\unicode[STIX]{x1D701}=0$
, the mean and the instantaneous conformation tensors have the same volume; when
$\unicode[STIX]{x1D701}$
is negative (positive), the instantaneous conformation tensor has a smaller (larger) volume than the volume of the mean. The logarithm ensures that there is no asymmetry between compressions and expansions with respect to the mean.
4.2.2 Squared distance from the mean,
$\unicode[STIX]{x1D705}$
When
$\unicode[STIX]{x1D63E}=\overline{\unicode[STIX]{x1D63E}}$
, we have
$\unicode[STIX]{x1D642}=\unicode[STIX]{x1D644}$
. When
$\unicode[STIX]{x1D63E}\neq \overline{\unicode[STIX]{x1D63E}}$
, we wish to consider the (appropriately defined) shortest distance between
$\unicode[STIX]{x1D644}$
and
$\unicode[STIX]{x1D642}$
as a measure of the magnitude of the fluctuation. The shortest path between
$\unicode[STIX]{x1D644}$
and
$\unicode[STIX]{x1D642}$
along the manifold
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
is given by the geodesic along
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
connecting
$\unicode[STIX]{x1D644}$
and
$\unicode[STIX]{x1D642}$
,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn48.gif?pub-status=live)
The squared geodesic distance associated with this path is then
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn49.gif?pub-status=live)
where (4.12) follows from (4.4). Using (4.4), one can verify that
$d^{2}(\unicode[STIX]{x1D644},\unicode[STIX]{x1D642})=d^{2}(\unicode[STIX]{x1D644},\unicode[STIX]{x1D642}^{-1})$
and thus the squared distance measure treats both expansions and compressions with respect to the mean similarly. The affine-invariance property, furthermore, ensures that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn50.gif?pub-status=live)
for all
$\unicode[STIX]{x1D63C}\in \boldsymbol{G}\boldsymbol{L}_{3}$
. In particular, with
$\unicode[STIX]{x1D63C}=\overline{\unicode[STIX]{x1D641}}$
, we obtain
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn51.gif?pub-status=live)
which exhibits the highly desirable property that the squared distance between
$\unicode[STIX]{x1D63E}$
and
$\overline{\unicode[STIX]{x1D63E}}$
is equal to the squared distance between
$\unicode[STIX]{x1D644}$
and
$\unicode[STIX]{x1D642}$
. A further consequence of the affine-invariance property is that
$d^{2}(\unicode[STIX]{x1D644},\unicode[STIX]{x1D642})$
is independent of the choice of the rotation
$\unicode[STIX]{x1D64D}\in \boldsymbol{S}\boldsymbol{O}_{3}$
in (3.4).
The path between
$\overline{\unicode[STIX]{x1D63E}}$
and
$\unicode[STIX]{x1D63E}$
along
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn52.gif?pub-status=live)
which reduces to
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn53.gif?pub-status=live)
when
$\unicode[STIX]{x1D64D}=\unicode[STIX]{x1D644}$
in (3.4). The choice
$\unicode[STIX]{x1D64D}=\unicode[STIX]{x1D644}$
is then natural in the sense that it allows the path along the manifold between
$\overline{\unicode[STIX]{x1D63E}}$
and
$\unicode[STIX]{x1D63E}$
, whose distance is a measure of the fluctuation, to be described using only
$\overline{\unicode[STIX]{x1D641}}$
and
$\unicode[STIX]{x1D642}$
.
We next consider realizability in the
$(\unicode[STIX]{x1D701},\unicode[STIX]{x1D705})$
plane. Since
$\boldsymbol{{\mathcal{G}}}$
is symmetric, its eigenvalues must be real, i.e. the eigenvalues must together belong in
$\mathbb{R}^{3}$
. In the
$\mathbb{R}^{3}$
space of eigenvalues of
$\boldsymbol{{\mathcal{G}}}$
, surfaces of constant
$\unicode[STIX]{x1D701}$
are planes, and surfaces of constant
$\unicode[STIX]{x1D705}$
are spheres. A particular choice of
$\unicode[STIX]{x1D701}$
and
$\unicode[STIX]{x1D705}$
is realizable only if the plane and sphere intersect. The coordinates along the intersecting circle of the sphere satisfy
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn54.gif?pub-status=live)
where
$\unicode[STIX]{x1D719}$
is the inclination angle and
$\unicode[STIX]{x1D703}$
is the azimuthal angle, in a spherical coordinate representation of
$\mathbb{R}^{3}$
. The physically realizable region in the
$(\unicode[STIX]{x1D705},\unicode[STIX]{x1D701})$
plane is thus given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn55.gif?pub-status=live)
When
$\unicode[STIX]{x1D705}=(1/3)\unicode[STIX]{x1D701}^{2}$
, the circle of intersection reduces to a point and
$\boldsymbol{{\mathcal{G}}}$
consequently has only two independent tensor invariants,
$\unicode[STIX]{x1D701}$
and
$\unicode[STIX]{x1D705}$
. The angles that maximize the left-hand side of (4.17) are given by
$\unicode[STIX]{x1D703}_{\max }=\unicode[STIX]{x03C0}/4$
,
$\unicode[STIX]{x1D719}_{\max }=\arctan \sqrt{2}$
. At
$(\unicode[STIX]{x1D703},\unicode[STIX]{x1D719})=(\unicode[STIX]{x1D703}_{\max },\unicode[STIX]{x1D719}_{\max })$
, it is readily verified that
$\unicode[STIX]{x1D705}=(1/3)\unicode[STIX]{x1D701}^{2}$
and also that the eigenvalues of
$\boldsymbol{{\mathcal{G}}}$
are all equal. Thus
$\boldsymbol{{\mathcal{G}}}$
, and hence
$\unicode[STIX]{x1D642}$
, is isotropic at the realizability bounds.
4.2.3 Anisotropy index,
$\unicode[STIX]{x1D709}$
Following the approach taken by Moakher & Batchelor (Reference Moakher and Batchelor2006), we define the anisotropy index,
$\unicode[STIX]{x1D709}$
, of
$\unicode[STIX]{x1D642}$
as the squared geodesic distance between
$\unicode[STIX]{x1D642}$
and the closest isotropic tensor,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn56.gif?pub-status=live)
By differentiation, we find that
$a^{3}=\prod _{i=1}^{3}\unicode[STIX]{x1D70E}_{i}(\unicode[STIX]{x1D642})=\text{det}\,\unicode[STIX]{x1D642}$
is a minimizing stationary point of (4.19) and hence the closest isotropic tensor to
$\unicode[STIX]{x1D642}$
along
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
is
$(\sqrt[3]{\text{det}\,\unicode[STIX]{x1D642}})\unicode[STIX]{x1D644}$
. We then have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn57.gif?pub-status=live)
Notice that
$\unicode[STIX]{x1D712}=0$
if and only if
$\unicode[STIX]{x1D701}^{2}=3\unicode[STIX]{x1D705}$
. But since we already showed that
$\boldsymbol{{\mathcal{G}}}$
, and hence
$\unicode[STIX]{x1D642}$
, are isotropic at the bound
$\unicode[STIX]{x1D701}^{2}=3\unicode[STIX]{x1D705}$
, it follows that
$\unicode[STIX]{x1D709}=0$
only for isotropic tensors.
Batchelor et al. (Reference Batchelor, Moakher, Atkinson, Calamante and Connelly2005) first introduced the index
$\sqrt{\unicode[STIX]{x1D709}}$
for characterizing positive definite diffusion tensors measured in magnetic resonance imaging. The index is analogous to the ‘fractional anisotropy index’ that is commonly used in turbulence and which provides the Euclidean distance to the closest isotropic tensor,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn58.gif?pub-status=live)
where
$\Vert \unicode[STIX]{x1D63C}\Vert _{F}=\text{tr}\,(\unicode[STIX]{x1D63C}^{\mathsf{T}}\boldsymbol{\cdot }\unicode[STIX]{x1D63C})$
indicates the Frobenius norm of matrix
$\unicode[STIX]{x1D63C}$
. A review of anisotropy measures is available in Moakher & Batchelor (Reference Moakher and Batchelor2006).
The three scalar measures presented above can be used together to obtain a better understanding of the fluctuations in the conformation tensor. The logarithmic volume ratio,
$\unicode[STIX]{x1D701}$
, is positive (negative) for volumetric expansions (contractions) with respect to the mean. However,
$\unicode[STIX]{x1D701}=0$
does not necessarily imply no deformation, since
$\text{det}(\text{det}(\unicode[STIX]{x1D642})\,\unicode[STIX]{x1D63C})=\text{det}(\unicode[STIX]{x1D642})$
for all
$\unicode[STIX]{x1D63C}$
with unit determinant. The squared geodesic distance to the identity,
$\unicode[STIX]{x1D705}$
, helps distinguish such cases since
$\unicode[STIX]{x1D705}=0$
only when
$\unicode[STIX]{x1D642}=\unicode[STIX]{x1D644}$
(
$\unicode[STIX]{x1D63E}=\overline{\unicode[STIX]{x1D63E}}$
). Finally, the anisotropy index,
$\unicode[STIX]{x1D709}$
, provides a quantification of the deviation of the shape of the polymer from the shape of the mean conformation tensor because it is a measure of the distance from
$\unicode[STIX]{x1D642}$
to the closest isotropic tensor, or equivalently, the minimizing distance between
$\unicode[STIX]{x1D63E}$
and
$a\overline{\unicode[STIX]{x1D63E}}$
over all
$a>0$
.
We next derive evolution equations for the scalar measures presented above and for the particular case when
$\overline{\unicode[STIX]{x1D63E}}=\langle \unicode[STIX]{x1D63E}\rangle$
.
4.3 Evolution equations for
$\unicode[STIX]{x1D701}$
,
$\unicode[STIX]{x1D705}$
and
$\unicode[STIX]{x1D709}$
Since
$\unicode[STIX]{x1D701}$
,
$\unicode[STIX]{x1D705}$
and
$\unicode[STIX]{x1D709}$
are scalar characterizations of
$\unicode[STIX]{x1D642}$
, one need only evolve
$\unicode[STIX]{x1D642}$
(or equivalently,
$\unicode[STIX]{x1D63E}$
) to obtain the field-valued
$\unicode[STIX]{x1D701}$
,
$\unicode[STIX]{x1D705}$
and
$\unicode[STIX]{x1D709}$
. Nevertheless, it is of interest to mathematically evaluate the evolution equations of these scalar measures separately in order to find the quantities that contribute to their dynamics.
Using (3.20) and the relationship
$\text{tr}\,\boldsymbol{{\mathcal{G}}}^{n}=\sum _{i=1}^{3}(\log \unicode[STIX]{x1D6E4}_{i})^{n}$
, we can derive the following equations for the fluctuating scalar measures:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn59.gif?pub-status=live)
where
$\text{dev}\,\boldsymbol{{\mathcal{G}}}=\boldsymbol{{\mathcal{G}}}-(\text{tr}\,\boldsymbol{{\mathcal{G}}}/3)\unicode[STIX]{x1D644}$
is the deviatoric part of
$\boldsymbol{{\mathcal{G}}}$
, and
$\boldsymbol{{\mathcal{D}}}$
is defined as
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn60.gif?pub-status=live)
The derivation of (4.22), omitted here for brevity, closely follows the procedure used by Vaithianathan & Collins (Reference Vaithianathan and Collins2003) to obtain evolution equations for the continuous eigendecomposition of
$\unicode[STIX]{x1D63E}$
.
The Cauchy–Schwarz inequality can be used to show that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn61.gif?pub-status=live)
The bounds in (4.24) illustrate the role of the stretching and relaxation balance,
$\boldsymbol{{\mathcal{D}}}$
, in bounding the growth of
$\unicode[STIX]{x1D705}$
and
$\unicode[STIX]{x1D709}$
. In the Reynolds-filtered case,
$\boldsymbol{{\mathcal{D}}}$
can be simplified using (3.17) so that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn62.gif?pub-status=live)
which shows that the turbulence intensity of the fluctuating conformation tensor, as measured by
$\unicode[STIX]{x1D705}$
, is not directly affected by the mean velocity gradient tensor
$\unicode[STIX]{x1D735}\overline{\boldsymbol{u}}$
. The contribution of
$\unicode[STIX]{x1D735}\overline{\boldsymbol{u}}$
to
$\unicode[STIX]{x1D705}$
is captured indirectly through
$\overline{\unicode[STIX]{x1D641}}$
, which is determined based on the mean balance.
According to (4.25), the tensor
$\boldsymbol{{\mathcal{D}}}$
consists of a stretching component:
$2\,\text{sym}(\unicode[STIX]{x1D640}^{\prime }-\langle \unicode[STIX]{x1D642}\boldsymbol{\cdot }\unicode[STIX]{x1D640}^{\prime }\rangle )$
, a component that arises due to gradients in
$\overline{\unicode[STIX]{x1D641}}$
and represents modified advection of
$\overline{\unicode[STIX]{x1D641}}$
:
$-2\,\text{sym}[(\boldsymbol{u}^{\prime }-\langle \unicode[STIX]{x1D642}\boldsymbol{\cdot }\boldsymbol{u}^{\prime }\rangle )\boldsymbol{\cdot }\unicode[STIX]{x1D735}\overline{\unicode[STIX]{x1D641}}]$
, a component that comprises mean advection of
$\unicode[STIX]{x1D642}$
by the fluctuating velocity field:
$\langle \boldsymbol{u}^{\prime }\boldsymbol{\cdot }\unicode[STIX]{x1D735}\unicode[STIX]{x1D642}\rangle$
and finally a component that resembles a fluctuating relaxation contribution:
$-(\unicode[STIX]{x1D648}\boldsymbol{\cdot }\unicode[STIX]{x1D642}^{-1}-\overline{\unicode[STIX]{x1D648}})$
.
4.3.1 Reynolds filtering the evolution equations
As a first-order statistical characterization of the fluctuating quantities,
$\unicode[STIX]{x1D701}$
,
$\unicode[STIX]{x1D705}$
and
$\unicode[STIX]{x1D709}$
, we will consider their filtered or averaged values,
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn63.gif?pub-status=live)
Reynolds filtering (4.22) and using the expression (4.25), we obtain the filtered evolution equations for
$\unicode[STIX]{x1D701}$
and
$\unicode[STIX]{x1D705}$
. The filtered equation for
$\unicode[STIX]{x1D701}$
is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn64.gif?pub-status=live)
where each term can be compared to those in (4.25) that were described in the previous subsection. Similarly, the filtered equation for
$\unicode[STIX]{x1D705}$
is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn65.gif?pub-status=live)
The equation for
$\unicode[STIX]{x1D709}$
, which we omit here for brevity, can be similarly derived.
5 Case study: viscoelastic turbulent channel flow
The general framework we have developed can be applied to a variety of flows. We focus on the classical problem of viscoelastic turbulent channel flow as a case study and use direct numerical simulations (DNS) to investigate the turbulent dynamics. The algorithmic details of the simulation are identical to that of Lee & Zaki (Reference Lee and Zaki2017) with the exception of the treatment of the conformation tensor which is documented in the appendix for the interested reader. The code employed was validated against linear growth rates of Tollmien–Schlichting waves (see Lee & Zaki Reference Lee and Zaki2017, for a study of natural transition in viscoelastic flows) and also against the results of Agarwal, Brandt & Zaki (Reference Agarwal, Brandt and Zaki2014) for the evolution of a localized disturbance.
We define
$x$
,
$y$
and
$z$
as the streamwise, wall-normal and spanwise coordinates, respectively. The flow is homogeneous in
$x$
and
$z$
and all the coordinates are normalized with respect to the channel half-height, with the channel walls located at
$y=\pm 1$
. The parameters for the calculation are listed in table 2. The Reynolds number,
$Re$
, is defined based on the channel half-height and bulk velocity, while
$Re_{\unicode[STIX]{x1D70F}}$
is the friction Reynolds number defined based on the friction velocity, calculated using the slope of the mean velocity at the wall, and channel half-height. The flow is driven by a pressure gradient which is adjusted in time to maintain a constant mass flow rate. The symbol
$\langle \cdot \rangle$
denotes averaging over
$x$
,
$z$
and
$t$
. Therefore, all of the averaged quantities are functions of only
$y$
.
The computational grid is uniform in the
$(x,z)$
directions and employs hyperbolic tangent stretching in the
$y$
direction with a Planck taper (McKechan, Robinson & Sathyaprakash Reference McKechan, Robinson and Sathyaprakash2010) applied such that grid spacing very close to the wall is constant. The maximum change of the grid spacing in the
$y$
direction is less than 3 % throughout in the domain. The resolution in friction units is listed in table 2. The initial turbulent state was generated from a separate simulation that followed the evolution of a Tollmien–Schlichting wave to the fully turbulent state (Lee & Zaki Reference Lee and Zaki2017). A snapshot from the fully turbulent state of Lee & Zaki (Reference Lee and Zaki2017) was used as an initial condition and first run for at least 150 convective time units before any statistics were collected. The evolution of the friction Reynolds number,
$Re_{\unicode[STIX]{x1D70F}}$
, was used to check whether the simulation had reached a statistically stationary state.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_fig2g.gif?pub-status=live)
Figure 2. Mean velocity profile from a FENE-P drag-reduced channel flow simulation. The solid line (——) is the mean streamwise velocity, the red dotted (lower) line (- - - - - -) is the von Kármán log law,
$\overline{u}_{von\,K\acute{a} rm\acute{a} n}^{+}=2.5y^{+}+5.5$
, and the red dashed (upper) line (– – –) is Virk’s asymptote,
$\overline{u}_{Virk}^{+}=11.7y^{+}-17.0$
.
Table 2. Parameters of the simulation of viscoelastic turbulent channel flow. The length scale is the channel half-height and velocity scale is the bulk flow speed. In the
$p$
th direction, the size of the domain is
$L_{p}$
, the number of grid points is
$N_{p}$
, and the spatial resolution, in friction units, is
$\unicode[STIX]{x1D6E5}_{p}^{+}$
.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_tab2.gif?pub-status=live)
5.1 Mean profiles and comparisons with the laminar profiles
The statistics presented in this section were obtained by averaging in space and over 750 time units, and by exploiting the symmetry of the flow about the centreline. Halving the number of samples maintained the trends and caused only minor deviations in the statistics, with no impact on the conclusions.
The mean streamwise velocity profile is shown in figure 2 as a function of
$y^{+}=Re_{\unicode[STIX]{x1D70F}}(y+1)$
, where
$Re_{\unicode[STIX]{x1D70F}}$
is always taken to be the turbulent frictional Reynolds number given in table 2. Also shown are the von Kármán log law and Virk’s maximum drag reduction asymptote. The mean velocity lies in between these two lines, indicating a drag-reduced state. Using Dean’s correlation for the skin friction (Dean Reference Dean1978), we obtain a friction Reynolds number of approximately 284 for a Newtonian flow with
$Re=4667$
and thus the drag reduction percentage is
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn66.gif?pub-status=live)
The non-zero components of
$\langle \unicode[STIX]{x1D63E}\rangle$
calculated for the same parameter values, are shown in figure 3. All the components of
$\langle \unicode[STIX]{x1D63E}\rangle$
are even functions of
$y$
except
$\langle \unicode[STIX]{x1D60A}_{xy}\rangle$
, which is an odd function of
$y$
. The streamwise stretch
$\langle \unicode[STIX]{x1D60A}_{xx}\rangle$
is four times larger than the laminar case (not shown) near the wall. It is also an order of magnitude larger than
$\langle \unicode[STIX]{x1D60A}_{yy}\rangle$
and
$\langle \unicode[STIX]{x1D60A}_{zz}\rangle$
. The remaining normal components of the conformation tensor are also larger in the turbulent case than the laminar: figure 3 shows that
$\max _{y}\langle \unicode[STIX]{x1D60A}_{yy}\rangle \approx 45$
and
$\max _{y}\langle \unicode[STIX]{x1D60A}_{zz}\rangle \approx 120$
, while
$\unicode[STIX]{x1D60A}_{yy}=\unicode[STIX]{x1D60A}_{zz}\approx 1$
throughout the channel when the flow is laminar. The peak values of each of the components
$\langle \unicode[STIX]{x1D60A}_{xx}\rangle$
,
$\langle \unicode[STIX]{x1D60A}_{yy}\rangle$
and
$\langle \unicode[STIX]{x1D60A}_{zz}\rangle$
occur at different locations in the channel. Figure 3 also shows that
$\langle \unicode[STIX]{x1D60A}_{zz}\rangle \geqslant \langle \unicode[STIX]{x1D60A}_{yy}\rangle$
throughout the channel. The trends above are consistent with those reported in the literature (Dallas, Vassilicos & Hewitt Reference Dallas, Vassilicos and Hewitt2010).
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_fig3g.gif?pub-status=live)
Figure 3. Mean conformation tensor profiles from a FENE-P drag-reduced channel flow simulation. (a)
$\overline{\unicode[STIX]{x1D60A}}_{xx}$
. (b) The solid line with star symbols (— * —) is
$\overline{\unicode[STIX]{x1D60A}}_{yy}$
, the dashed line with square symbols (– – ▫ – –) is
$\overline{\unicode[STIX]{x1D60A}}_{zz}$
and the dashed-dot line with circle symbols (– - ○ – -) is
$\overline{\unicode[STIX]{x1D60A}}_{xy}$
. The remaining components of the mean conformation tensor are 0. Note that the symbols in (b) are identifiers and are thus only a small subset of all the data points used in the line plots.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_fig4g.gif?pub-status=live)
Figure 4. Scalar measures applied to nominal conformation tensors, in equivalent dimensions, plotted as functions of
$y^{+}$
. (a) The solid line with star symbols (— * —) is
$\log \det \,\overline{\unicode[STIX]{x1D63E}}$
, the logarithmic volume; the dashed line with square symbols (– – ▫ – –) is
$d(\unicode[STIX]{x1D644},\overline{\unicode[STIX]{x1D63E}})$
, the geodesic distance between
$\overline{\unicode[STIX]{x1D63E}}$
and
$\unicode[STIX]{x1D644}$
on
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
. (b) Anisotropy index,
$\sqrt{(d(\unicode[STIX]{x1D644},\overline{\unicode[STIX]{x1D63E}}))^{2}-(1/3)(\log \det \,\overline{\unicode[STIX]{x1D63E}})^{2}}$
, which is the geodesic distance from the closest isotropic tensor. For both (a) and (b), black lines are for
$\overline{\unicode[STIX]{x1D63E}}=\langle \unicode[STIX]{x1D63E}\rangle$
and grey lines are when
$\overline{\unicode[STIX]{x1D63E}}$
is equal to the FENE-P laminar conformation tensor. The red dotted line (- - - - - -) in (b) is
$-1.375\log y^{+}+7.925$
. Note that the symbols in (a) are identifiers and are thus only a small subset of all the data points used in the line plots.
Figure 4(a) shows the logarithmic volume of the mean and laminar conformation tensors along with the distance from the origin on the manifold of positive definite tensors. The figure shows that both the logarithmic volume,
$\log \det \,\overline{\unicode[STIX]{x1D63E}}$
, and the distance from the origin,
$d(\unicode[STIX]{x1D644},\overline{\unicode[STIX]{x1D63E}})$
are larger in the turbulent case compared to the laminar. Furthermore, these two quantities are monotonically decreasing in the laminar case but have peaks around
$y^{+}\approx 60$
in the turbulent case. Both quantities in the two cases asymptote to a constant at locations very close to the wall,
$y^{+}\leqslant 2$
. The weak growth in
$d(\unicode[STIX]{x1D644},\overline{\unicode[STIX]{x1D63E}})$
in the turbulent case despite a rapid increase in
$\log \det \,\overline{\unicode[STIX]{x1D63E}}$
is due to the increase in isotropy (sphericity) as we move away from the wall, since for a given volume the tensor closest to
$\unicode[STIX]{x1D644}$
is an isotropic tensor. In figure 3 we see that the mean normal stretches in the
$y$
and
$z$
directions in the turbulent case peak between
$y^{+}\approx 40$
and
$y^{+}\approx 70$
and are at least an order of magnitude larger than the stretches in the laminar case, where
$\overline{\unicode[STIX]{x1D60A}}_{yy}=\overline{\unicode[STIX]{x1D60A}}_{zz}\approx 1$
. The term
$\overline{\unicode[STIX]{x1D60A}}_{xx}$
is decreasing towards the channel centre in both the turbulent and laminar case but is accompanied by an increase in
$\overline{\unicode[STIX]{x1D60A}}_{yy}$
,
$\overline{\unicode[STIX]{x1D60A}}_{zz}$
in the turbulent flow which leads to increased isotropy for locations sufficiently removed from the wall.
Figure 4(b) shows the anisotropy, the geodesic distance to the closest isotropic tensor, of the mean and laminar conformation tensors. The anisotropy index is approximately constant in the vicinity of the wall for both the laminar and turbulent cases, and decays away from the wall. In the turbulent flow, the decay starts very close to the wall – approximately three friction units away from the wall – and then shows a remarkable logarithmic decay that proceeds all the way to very close to the centreline where it sharply turns and forms a stationary point. The increased isotropy in the turbulent case may be explained by the fact that, although more stretching occurs in this case, the stretching in the cross-stream directions is much larger than in the laminar case. Overall, this leads to a more isotropic mean conformation tensor.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20180323123059-63880-mediumThumb-S0022112018001180_fig5g.jpg?pub-status=live)
Figure 5. Wall-parallel
$(x,z)$
planes of isocontours of instantaneous
$\text{I}_{\unicode[STIX]{x1D642}}=\text{tr}\,\unicode[STIX]{x1D642}$
: (a)
$y^{+}=15$
and (b)
$y^{+}=180$
(centreline).
5.2 Invariants of the fluctuating conformation tensor
In order to motivate the scalar measures proposed in the present work, we consider the invariants of
$\unicode[STIX]{x1D642}$
as alternatives in this subsection. Figure 5 shows isocontours of instantaneous
$\text{I}_{\unicode[STIX]{x1D642}}$
at a given time at two wall-parallel planes,
$y^{+}=15$
and
$y^{+}=180$
(centreline). The isocontours of instantaneous
$\text{II}_{\unicode[STIX]{x1D642}}$
and
$\text{III}_{\unicode[STIX]{x1D642}}$
are qualitatively similar to those of
$\text{I}_{\unicode[STIX]{x1D642}}$
and are thus not shown here. The instantaneous
$\text{I}_{\unicode[STIX]{x1D642}}$
can vary over several orders of magnitude. As a result, obtaining reliable statistics for the invariants is challenging. We found that the peak root-mean-square (r.m.s.) of the invariants (not shown) are at least an order of magnitude larger than the corresponding mean values. This large spread in the instantaneous invariants of
$\unicode[STIX]{x1D642}$
suggests that
$\log \unicode[STIX]{x1D642}$
is a more appropriate quantity to consider, and reinforces the need for the geometrically consistent scalar measures introduced in § 4.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary-alt:20180323123059-11875-mediumThumb-S0022112018001180_fig6g.jpg?pub-status=live)
Figure 6. Wall-parallel
$(x,z)$
planes of isocontours of instantaneous (a,b) logarithmic volume ratio,
$\unicode[STIX]{x1D701}$
, (c,d) geodesic distance from the identity,
$\unicode[STIX]{x1D705}$
and (e,f) anisotropy index,
$\unicode[STIX]{x1D709}$
. (a,c,e)
$y^{+}=15$
, (b,d,f)
$y^{+}=180$
(centreline).
5.3 The scalar measures:
$\unicode[STIX]{x1D701}$
,
$\unicode[STIX]{x1D705}$
and
$\unicode[STIX]{x1D709}$
Figure 6 shows isocontours of instantaneous values of
$\unicode[STIX]{x1D701}$
,
$\unicode[STIX]{x1D705}$
and
$\unicode[STIX]{x1D709}$
for two wall-parallel planes,
$y^{+}=15$
and
$y^{+}=180$
. As a comparison, the fluctuating tensor
$\unicode[STIX]{x1D60A}_{xx}^{\prime }$
obtained by the Reynolds decomposition (1.2) and normalized by the local
$\overline{\unicode[STIX]{x1D60A}}_{xx}$
is shown in figure 7. We normalized
$\unicode[STIX]{x1D60A}_{xx}^{\prime }$
so that the fluctuations near the wall could be compared to those at the centreline, since
$\overline{\unicode[STIX]{x1D60A}}_{xx}$
differs by an order of magnitude between the two locations. The isocontours of
$\unicode[STIX]{x1D60A}_{xx}^{\prime }/\overline{\unicode[STIX]{x1D60A}}_{xx}$
, and in particular the negative values, are difficult to interpret since the correspondence to a physical deformation or a mathematical metric is unclear.
Figure 6(a,b) shows the logarithmic volume ratio,
$\unicode[STIX]{x1D701}$
. This quantity is the logarithm of
$\text{III}_{\unicode[STIX]{x1D642}}$
, which itself is qualitatively similar to
$\text{I}_{\unicode[STIX]{x1D642}}$
and hence we observe a strong visual resemblance between figures 5 and 6(a,b). The colour scale in the former is logarithmic and is thus consistent with the linear scale in figure 6. Both figures 6(a) and 6(b) have predominantly negative values, indicating that the instantaneous volume is smaller than the volume of the mean conformation. We also find regions of very high
$\unicode[STIX]{x1D701}$
adjacent to regions of very low values, especially at the centreline. This is partially a result of the lack of diffusion in the polymers since there is no direct mechanism for smoothing out shocks in the conformation tensor field.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_fig7g.gif?pub-status=live)
Figure 7. Wall-parallel
$(x,z)$
planes of isocontours of instantaneous
$\unicode[STIX]{x1D60A}_{xx}^{\prime }/\overline{\unicode[STIX]{x1D60A}}_{xx}$
at (a)
$y^{+}=15$
, (b)
$y^{+}=180$
(centreline), where
$\overline{\unicode[STIX]{x1D60A}}_{xx}(y^{+}=15)=2.22\times 10^{3}$
and
$\overline{\unicode[STIX]{x1D60A}}_{xx}(y^{+}=180)=1.02\times 10^{2}$
. The limits of the divergent colour map are set at the planar maxima and minima of
$\unicode[STIX]{x1D60A}_{xx}^{\prime }/\overline{\unicode[STIX]{x1D60A}}_{xx}$
.
Another important effect is that of memory: polymers are stretched near the wall where the shear is significant, and are then transported out to the centreline. If the half-channel transit time, the ratio of the channel half-height to the r.m.s. wall-normal velocity, is smaller than the polymer relaxation time, we expect to observe a footprint of the near-wall stretching all the way out to the centreline. In the present case, the relaxation time is two to three times the half-channel transit time. Since material points that are initially close are exponentially diverging in a turbulent flow (see Johnson et al. Reference Johnson, Hamilton, Burns and Meneveau2017, for a study of Lagrangian stretching in Newtonian turbulent channel flow), it is unsurprising to find adjacent regions of strongly and weakly stretched polymers.
The reductionist explanation given above is useful for a basic understanding but is insufficient to account for other observed features of the flow. For example, the present considerations would suggest that the streamwise elongated shape of the isocontours of
$\unicode[STIX]{x1D701}$
near the wall would lead to a similar shape at the centreline. However, this is manifestly not the case. Instead, the
$\unicode[STIX]{x1D701}$
field appears to generate, on the whole, highly curved isocontours at the centre of the channel.
The measure
$\unicode[STIX]{x1D701}$
does not distinguish between volume-preserving deformations. For example,
$\unicode[STIX]{x1D701}$
does not distinguish between
$\unicode[STIX]{x1D63E}$
and
$(\det \,\unicode[STIX]{x1D63E})\unicode[STIX]{x1D63C}$
for any
$\unicode[STIX]{x1D63C}$
with determinant
$=1$
. In particular,
$\unicode[STIX]{x1D701}=0$
, does not imply
$\unicode[STIX]{x1D63E}=\overline{\unicode[STIX]{x1D63E}}$
. In order to identify regions where
$\unicode[STIX]{x1D63E}=\overline{\unicode[STIX]{x1D63E}}$
is true, and quantify the deviation when it is not, we use the squared distance away from the origin (
$\unicode[STIX]{x1D644}$
) along the manifold,
$\unicode[STIX]{x1D705}$
. Figure 6(c,d) shows isocontours of instantaneous
$\unicode[STIX]{x1D705}$
. Most of the conformation tensor field is significantly far away, in the sense of distance along the manifold, from
$\overline{\unicode[STIX]{x1D63E}}$
. However, regions where
$\overline{\unicode[STIX]{x1D63E}}$
is a good representation of
$\unicode[STIX]{x1D63E}$
are interspersed between regions where
$\unicode[STIX]{x1D705}$
is large. This behaviour is true both in the near-wall region as well as the channel centre but more so in the latter. In contrast, it is well known that kinetic energy fluctuations are weakest at the centreline in a Newtonian channel flow. A different behaviour for the polymers is unsurprising since, due to the strong memory effect,
$\unicode[STIX]{x1D63E}$
at each point is strongly dependent on the Lagrangian path that is obtained by a pull-back of the particular Eulerian point of interest.
Figure 6(e,f) shows isocontours of instantaneous
$\unicode[STIX]{x1D709}$
, the anisotropy index. This index shows how close the shape of instantaneous conformation tensor is to the shape of the mean conformation tensor, irrespective of volumetric changes. The visual resemblance of
$\unicode[STIX]{x1D705}$
and
$\unicode[STIX]{x1D709}$
suggests that deformations to the mean conformation are largely anisotropic, or in other words, lead to shape change.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_fig8g.gif?pub-status=live)
Figure 8. Mean scalar measures, plotted in equivalent dimensions, as functions of
$y^{+}$
. The solid line with star symbols (— * —) is minus the mean volume ratio,
$-\overline{\unicode[STIX]{x1D701}}=-\langle \unicode[STIX]{x1D701}\rangle$
; the dashed line with square symbols (– – ▫ – –) is the square-root of the mean geodesic distance from the identity,
$\overline{\unicode[STIX]{x1D705}}^{1/2}=\sqrt{\langle \unicode[STIX]{x1D705}\rangle }$
; the dashed-dot line with circle symbols (– - ○ – -) is the square-root of the mean anisotropy index,
$\overline{\unicode[STIX]{x1D709}}^{1/2}=\sqrt{\langle \unicode[STIX]{x1D709}\rangle }$
. Red dotted lines (- - - - - -) are logarithmic fits to the profiles: the fit to the
$\overline{\unicode[STIX]{x1D705}}^{1/2}$
profile (– – ▫ – –) is given by
$0.725\log y^{+}+2.15$
, the fits to the
$\overline{\unicode[STIX]{x1D709}}^{1/2}$
profile (– - ○ – -) are given by
$0.65\log y^{+}+1.20$
and
$-0.9\log y^{+}+7.6$
. Note that the symbols are identifiers and are thus only a small subset of all the data points used in the line plots.
Figure 8 shows the mean values of
$\unicode[STIX]{x1D701}$
,
$\unicode[STIX]{x1D705}$
and
$\unicode[STIX]{x1D709}$
in dimensions of distance along the manifold. These statistics were generated using 225 convective time units. We checked for convergence by halving the number of samples. This process led to only minor deviations in the results, with no material significance to the discussion that follows. As was inferred earlier, the average logarithmic volume ratio is negative throughout the channel and is monotonically decreasing towards the centreline where it becomes roughly constant, similar to the behaviour very near the wall
$y^{+}\lesssim 2$
. This behaviour of the mean logarithmic volume being smaller than the volume of the mean is consistent with the ‘swelling’ problem associated with the arithmetic mean of positive definite tensors that has been previously reported in the literature (Arsigny et al.
Reference Arsigny, Fillard, Pennec and Ayache2007). It also suggests that, although an arithmetic mean of the conformation tensor may be unavoidable for modelling in the averaged equations, it may not be the most representative conformation tensor for deducing the most likely physical deformation of the polymers. A more extensive study, beyond the scope of the present work, is required to determine better alternatives to the arithmetic mean.
The square-root mean squared distance from the origin along the manifold,
$\overline{\unicode[STIX]{x1D705}}^{1/2}$
, is logarithmically increasing up to close to the centreline but peaks at
$y^{+}\approx 100$
. The anisotropy,
$\overline{\unicode[STIX]{x1D709}}^{1/2}$
, shows logarithmic increase over a small range
$3\lesssim y^{+}\lesssim 20$
, peaking at
$y^{+}\approx 60$
, but then shows a logarithmic decrease towards the centreline. The logarithmic behaviour in these quantities, especially in
$\overline{\unicode[STIX]{x1D705}}^{1/2}$
where the behaviour extends over a significant range, resembles the behaviour that appears in the mean velocity as well as in the statistical moments and two-point correlations of the velocity fluctuations in wall-bounded shear flows (Meneveau & Marusic Reference Meneveau and Marusic2013; Yang, Marusic & Meneveau Reference Yang, Marusic and Meneveau2016).
If we momentarily accept the simplified picture of polymers being deformed closer to the wall in an ‘active’ region and then passively transported out to the centreline of the channel, then these results indicate that the active region of the channel exists all the way up to
$y^{+}\approx 100$
. Beyond this region, the stretching of polymers weakens and thus
$\overline{\unicode[STIX]{x1D705}}^{1/2}$
decreases monotonically. The active region involves a region of logarithmic increase, and so we can write
$\overline{\unicode[STIX]{x1D705}}^{1/2}=a_{1}\int y^{-1}\,\text{d}y+a_{0}$
for constants
$a_{1}$
and
$a_{0}$
. If the stretching at each wall-normal station is actually additive, the polymers undergo deformation that on average leads to a diminishing increment in the deformation. This behaviour is consistent with the velocity gradients weakening with wall-normal distance. At
$y^{+}\approx 100$
, the velocity gradients then either weaken or act on such long time scales that polymers relax quickly enough not to retain any additional deformation.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_fig9g.gif?pub-status=live)
Figure 9. Joint probability density functions (JPDF) of the logarithmic volume ratio,
$\unicode[STIX]{x1D701}$
, and the geodesic distance from the identity,
$\unicode[STIX]{x1D705}$
, at four different wall-normal locations: (a)
$y^{+}=2$
, (b)
$y^{+}=15$
, (c)
$y^{+}=100$
and (d)
$y^{+}=180$
(centreline). Dotted lines (- - - - -) are isocontours of the anisotropy index,
$\unicode[STIX]{x1D709}$
. The thick red dashed line (– – –) denotes the realizability bound,
$\unicode[STIX]{x1D705}=(1/3)\unicode[STIX]{x1D701}^{2}$
, derived in (4.18), which coincides with the zero anisotropy index isocontour (
$\unicode[STIX]{x1D709}=0$
).
In order to quantify the fluctuations observed in figure 6 in more detail, we calculated the joint probability density function (JPDF) for
$\unicode[STIX]{x1D701}$
and
$\unicode[STIX]{x1D705}$
using 12 snapshots evenly spaced over 120 convective time units. The JPDFs, at four different wall-normal locations, are shown in figure 9 along with isocontours of
$\unicode[STIX]{x1D709}$
, which is purely a function of
$\unicode[STIX]{x1D701}$
and
$\unicode[STIX]{x1D705}$
.
The JPDFs are non-zero primarily on the lower half of the realizability region, which is consistent with the isocontours in figure 6(a,b) and
$\overline{\unicode[STIX]{x1D701}}<0$
throughout the channel in figure 8. In addition, the isocontours tend to concentrate along the isotropy line (thick red dashed line) that was derived in (4.18) as a realizability bound. However, with the exception of the centreline, the most probable
$(\unicode[STIX]{x1D701},\unicode[STIX]{x1D705})$
are located away from the isotropy line. This implies that the most likely conformation tensor away from the centreline does not have the same shape as the local mean conformation tensor. Although the most likely conformation tensor at the centreline assumes the shape of the mean conformation tensor, the JPDF at the centreline occupies a greater area in the
$(\unicode[STIX]{x1D701},\unicode[STIX]{x1D705})$
plane than at
$y^{+}=2$
or
$y^{+}=15$
, which implies a greater degree of uncertainty. In addition, the most likely conformation tensor, as determined by the peak of the JPDF, is further away from the mean than at any other wall-normal location.
The JPDF indicate that the most intermittent region of the flow, determined by the most extreme excursions away from the identity on
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
, do not occur near the wall or at the centreline. This can be seen in figure 9, where the JPDF at
$y^{+}=100$
shows events with up to
$\unicode[STIX]{x1D705}=100$
. This behaviour is consistent with the peak
$\overline{\unicode[STIX]{x1D705}}$
occurring away from the centreline in figure 8.
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_fig10g.gif?pub-status=live)
Figure 10. Root-mean-square profiles of
$\unicode[STIX]{x1D63E}^{\prime }$
,
$[\unicode[STIX]{x1D60A}_{ij}]_{rms}=\sqrt{\langle \unicode[STIX]{x1D60A}_{ij}^{2}\rangle -\langle \unicode[STIX]{x1D60A}_{ij}\rangle ^{2}}$
based on the Reynolds decomposition (1.2). (a)
$[\unicode[STIX]{x1D60A}_{xx}]_{rms}$
. (b) The solid line with star symbols (— * —) is
$[\unicode[STIX]{x1D60A}_{yy}]_{rms}$
, the dashed line with square symbols (– – ▫ – –) is
$[\unicode[STIX]{x1D60A}_{zz}]_{rms}$
and the dashed-dot line with circle symbols (– - ○ – -) is
$[\unicode[STIX]{x1D60A}_{xy}]_{rms}$
. Note that the symbols in (b) are identifiers and are thus only a small subset of all the data points used in the line plots.
Finally the r.m.s. of
$\unicode[STIX]{x1D63E}^{\prime }$
, defined according to the Reynolds decomposition (1.2), is shown in figure 10 for comparison to the present approach. The protocol used to obtain these quantities was the same as that used to obtain the mean conformation tensor in figure 3. The r.m.s. values of
$\unicode[STIX]{x1D63E}^{\prime }$
are of similar magnitude to the components of
$\langle \unicode[STIX]{x1D63E}\rangle$
. In fact, the peak r.m.s. of
$\unicode[STIX]{x1D60A}_{yy}$
,
$\unicode[STIX]{x1D60A}_{zz}$
and
$\unicode[STIX]{x1D60A}_{xy}$
are larger than their respective mean values. Interestingly, the peak fluctuating deformation found at
$y^{+}\approx 100$
using our present framework is not discernible from the r.m.s. fluctuations. Different components of the r.m.s. tensor peak at different locations in the channel with
$[\unicode[STIX]{x1D60A}_{yy}]_{rms}$
showing a peak that is closest to
$y^{+}=100$
. The r.m.s. quantities only show the component-wise behaviour of the conformation tensor, and hence are not indicative of the total polymer deformation. A more appropriate quantity to evaluate in the context of
$\unicode[STIX]{x1D63E}^{\prime }$
would then be the JPDF of all six independent components of
$\unicode[STIX]{x1D63E}^{\prime }$
. Owing to high-dimensionality, this characterization is more difficult to both calculate and analyse. The scalar measures suggested in the present work provide a good alternative to such a characterization.
6 Conclusion
We have developed a geometric decomposition, given in (3.6), that overcomes the difficulties associated with the traditional Reynolds decomposition of
$\unicode[STIX]{x1D63E}$
. The geometric decomposition yields a conformation tensor,
$\unicode[STIX]{x1D642}$
, that describes the deformation of the polymer with respect to the mean deformation. We characterized the fluctuations in
$\unicode[STIX]{x1D642}$
by using a geometry specifically constructed for
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$
and obtained three scalar measures: the logarithmic volume ratio,
$\unicode[STIX]{x1D701}$
, given in (4.10), the squared geodesic distance of the perturbation conformation tensor away from the origin,
$\unicode[STIX]{x1D705}$
, given in (4.12) and the anisotropy index,
$\unicode[STIX]{x1D709}$
, given in (4.20), defined as the squared geodesic distance to the closest isotropic tensor. The average values and JPDFs of these scalar measures provided interesting insights about the fluctuating polymer deformation that are not readily available from a Reynolds decomposition of
$\unicode[STIX]{x1D63E}$
. These insights include the following.
-
(i) The anisotropy in
$\overline{\unicode[STIX]{x1D63E}}$ , measured as geodesic distance away from
$\unicode[STIX]{x1D644}$ on
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{3}$ , decreases logarithmically from
$y^{+}=5$ to close to the centreline.
-
(ii) The mean conformation tensor tends to be significantly different to the most likely conformation tensor observed in the flow.
-
(iii) The mean polymer deformation, measured in terms of
$\overline{\unicode[STIX]{x1D705}}$ , increases logarithmically from
$y^{+}=10$ and peaks at
$y^{+}\approx 100$ .
-
(iv) As evidenced by the JPDF of
$\unicode[STIX]{x1D705}$ , the peak turbulence intensity in the polymers occurs in between the wall and centreline, at approximately
$y^{+}\approx 100$ .
The universality of the trends mentioned above, and others documented in the present work, and their connection to larger issues in viscoelastic turbulence are open questions. The framework we have developed can be used to probe the dynamics in viscoelastic turbulence beyond channel flow and can also be exploited for developing and benchmarking reduced-order models for viscoelastic turbulence. The approach can also be adapted to other similar problems, for example in the analysis of deforming droplets in turbulence using a model based on the droplet conformation tensor (Maffettone & Minale Reference Maffettone and Minale1998; Biferale, Meneveau & Verzicco Reference Biferale, Meneveau and Verzicco2014).
An important, open question that needs to be resolved in future work is the relationship between the fluctuating conformation tensor
$\unicode[STIX]{x1D642}$
and elastic energy of the polymers. In contrast to the clear meaning of the kinetic energy associated with the fluctuating velocity field, a deeper understanding of the elastic energy and its relation to
$\unicode[STIX]{x1D642}$
and the scalar measures introduced in the present work is unavailable. The attainment of such an understanding is partially hindered by the myriad of constitutive models prevalent in the literature (Beris & Edwards Reference Beris and Edwards1994). Instead of using the details of a particular constitutive model, the aim of the present work was to maintain as much generality as possible by exploiting the mathematical structure of
$\unicode[STIX]{x1D642}$
to characterize the fluctuating polymer deformation.
Acknowledgements
The authors would like to acknowledge the Maryland Advanced Research Computing Centre (MARCC) for computing time, S. J. Lee and J. Lee for help in the validation of the numerical code, and the referees for their valuable comments. This work is sponsored in part by the National Science Foundation under grants CBET–1652244, CBET–1511937 and PHY–1125915.
Appendix A. Riemannian structure of the set of positive definite matrices
The theoretical results presented in this section on the geometric structure of
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
are standard with detailed accounts available in pp. 322–339 of Lang (Reference Lang2001) and also pp. 201–235 of Bhatia (Reference Bhatia2015).
A.1 Riemannian metric
We can define an inner product
$(\cdot ,\cdot ):\mathbb{R}^{n\times n}\times \mathbb{R}^{n\times n}\rightarrow \mathbb{R}$
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn67.gif?pub-status=live)
for
$\unicode[STIX]{x1D63C},\unicode[STIX]{x1D63D}\in \mathbb{R}^{n\times n}$
where
$\unicode[STIX]{x1D653}\in \boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
is fixed. When
$\unicode[STIX]{x1D653}=\unicode[STIX]{x1D644}$
, (A 1) reduces to the definition of the standard Frobenius inner product. The guaranteed factorization
$\unicode[STIX]{x1D653}=\unicode[STIX]{x1D653}^{1/2}\boldsymbol{\cdot }\unicode[STIX]{x1D653}^{1/2}$
and the cyclical property of the trace ensures the positivity of
$\Vert \unicode[STIX]{x1D63C}\Vert _{\unicode[STIX]{x1D653}}$
, while the remaining properties of the inner product and the norm follow on from the standard Frobenius theory. The space
$\mathbb{R}^{n\times n}$
is a Hilbert space when equipped with such an inner product and the norm induced by it:
$\Vert \unicode[STIX]{x1D63C}\Vert _{\unicode[STIX]{x1D653}}=\sqrt{(\unicode[STIX]{x1D63C},\unicode[STIX]{x1D63C})_{\unicode[STIX]{x1D653}}}$
. The subset of
$\mathbb{R}^{n\times n}$
consisting of symmetric matrices forms a vector space,
$\boldsymbol{S}\boldsymbol{y}\boldsymbol{m}_{n}$
, and can also be Hilbertized under the inner product (A 1).
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
is an open subset of
$\mathbb{R}^{n\times n}$
(and also
$\boldsymbol{S}\boldsymbol{y}\boldsymbol{m}_{n}$
) in the
$\Vert \cdot \Vert _{\unicode[STIX]{x1D653}}$
metric and is thus a (smooth) manifold.
The tangent bundle of
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
, which consists of the manifold
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
equipped with a tangent space
$T_{\unicode[STIX]{x1D653}}\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
at each point
$\unicode[STIX]{x1D653}$
of
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
, provides a natural projection that can be used to study the geometry of
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
. A simple argument shows that the tangent space at each point of
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
coincides with
$\boldsymbol{S}\boldsymbol{y}\boldsymbol{m}_{n}$
. (Let
$\unicode[STIX]{x1D654}$
be defined by
$\unicode[STIX]{x1D654}=\unicode[STIX]{x1D653}+\unicode[STIX]{x1D700}\unicode[STIX]{x1D64E}$
, for some
$\unicode[STIX]{x1D653}\in \boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
,
$\unicode[STIX]{x1D700}\in \mathbb{R}$
,
$\unicode[STIX]{x1D64E}\in \boldsymbol{S}\boldsymbol{y}\boldsymbol{m}_{n}$
. By Weyl’s inequality, there exists some
$\unicode[STIX]{x1D700}>0$
sufficiently small such that
$\unicode[STIX]{x1D654}\in \boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
. This implies that
$\boldsymbol{S}\boldsymbol{y}\boldsymbol{m}_{n}\subseteq T_{\unicode[STIX]{x1D653}}\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
. Since
$\unicode[STIX]{x1D654}\notin \boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
for any
$\unicode[STIX]{x1D700}\neq 0$
and
$\unicode[STIX]{x1D64E}\notin \boldsymbol{S}\boldsymbol{y}\boldsymbol{m}_{n}$
, we have
$\boldsymbol{S}\boldsymbol{y}\boldsymbol{m}_{n}=T_{\unicode[STIX]{x1D653}}\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
.) The latter result is the geometric underpinning of numerical algorithms that time-march the conformation tensor by translations of
$\unicode[STIX]{x1D63E}$
by symmetric matrices (the right-hand side of the evolution equation for
$\unicode[STIX]{x1D63E}$
).
A manifold
$\boldsymbol{M}$
equipped with a scalar product over
$T_{\unicode[STIX]{x1D653}}\boldsymbol{M}$
for each
$\unicode[STIX]{x1D653}\in \boldsymbol{M}$
is a Riemannian manifold. The set of such scalar products is called the Riemannian metric of the manifold.
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
is a Riemannian manifold with Riemannian metric given by (Lang Reference Lang2001; Bhatia Reference Bhatia2015)
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn68.gif?pub-status=live)
with the understanding that the scalar product on
$T_{\unicode[STIX]{x1D653}}\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
is
$(\cdot ,\cdot )_{\unicode[STIX]{x1D653}}\in g$
and the domain of
$(\cdot ,\cdot )_{\unicode[STIX]{x1D653}}$
is restricted to
$\boldsymbol{S}\boldsymbol{y}\boldsymbol{m}_{n}=T_{\unicode[STIX]{x1D653}}\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
. An infinitesimal distance around the point
$\unicode[STIX]{x1D653}$
on the manifold is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn69.gif?pub-status=live)
The metric given by (A 2) ensures that distances between points
$\unicode[STIX]{x1D653},\unicode[STIX]{x1D654}\in \boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
along the manifold calculated using (A 3) are invariant under the action
$[\cdot ]_{\unicode[STIX]{x1D63C}}$
of any
$\unicode[STIX]{x1D63C}\in \boldsymbol{G}\boldsymbol{L}_{n}$
, i.e. invariant to transformations such as (3.6).
A.2 Geodesic curves and distances
Consider a parametrized curve on
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
connecting points
$\unicode[STIX]{x1D653},\unicode[STIX]{x1D654}\in \boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
, i.e.
$\unicode[STIX]{x1D64B}:[0,1]\rightarrow \boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
with
$\unicode[STIX]{x1D64B}(0)=\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D64B}(1)=\unicode[STIX]{x1D654}$
. The distance, in the sense of the metric
$g$
, traversed on the manifold along the curve
$\unicode[STIX]{x1D64B}=\unicode[STIX]{x1D64B}(r)$
is given by
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn70.gif?pub-status=live)
$\ell _{\unicode[STIX]{x1D64B}}$
has an attractive property in that it is invariant under affine transformations.
Lemma 1 (Affine invariance).
For each positive definite
$\unicode[STIX]{x1D63C}$
and differentiable path
$\unicode[STIX]{x1D64B}$
on the Riemannian manifold of positive definite matrices, we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn71.gif?pub-status=live)
Proof. § 6.1.1 in Bhatia (Reference Bhatia2015).
We call a curve
$\unicode[STIX]{x1D64B}(r)$
on
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
that minimizes
$\ell _{\unicode[STIX]{x1D64B}}(1)$
a geodesic curve connecting
$\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D654}$
. In general, the existence and/or uniqueness of a geodesic curve is not guaranteed. We also define
$d(\unicode[STIX]{x1D653},\unicode[STIX]{x1D654})$
, the geodesic distance between
$\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D654}$
, as the infimum of
$\ell _{\unicode[STIX]{x1D64B}}(1)$
over all possible curves
$\unicode[STIX]{x1D64B}$
connecting
$\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D654}$
:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn72.gif?pub-status=live)
A corollary of affine invariance is that
$d(\unicode[STIX]{x1D653},\unicode[STIX]{x1D654})=d([\unicode[STIX]{x1D653}]_{\unicode[STIX]{x1D63C}},[\unicode[STIX]{x1D654}]_{\unicode[STIX]{x1D63C}})$
.
It turns out that the existence and uniqueness of geodesics is guaranteed on
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
. Furthermore, we can obtain analytical expressions for these geodesics. Following Bhatia (Reference Bhatia2015), we present three key theorems that allow this construction.
Theorem 1 (Exponential metric increasing property).
For any two real symmetric
$\boldsymbol{{\mathcal{X}}}$
and
$\boldsymbol{{\mathcal{Y}}}$
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn73.gif?pub-status=live)
where we note that
$\text{e}^{\boldsymbol{{\mathcal{X}}}}$
,
$\text{e}^{\boldsymbol{{\mathcal{Y}}}}$
are positive definite matrices.
Proof. § XII.2 in Lang (Reference Lang2001) and § 6.1.4 in Bhatia (Reference Bhatia2015).
Equality is achieved in (A 7) when
$\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D654}$
commute and we can also parametrize the geodesic in this case, as expressed in the proposition below.
Proposition 1. Let
$\unicode[STIX]{x1D653}=\text{e}^{\boldsymbol{{\mathcal{X}}}}$
and
$\unicode[STIX]{x1D654}=\text{e}^{\boldsymbol{{\mathcal{Y}}}}$
be positive definite matrices such that
$\unicode[STIX]{x1D653}\boldsymbol{\cdot }\unicode[STIX]{x1D654}=\unicode[STIX]{x1D654}\boldsymbol{\cdot }\unicode[STIX]{x1D653}$
. Then, the exponential function maps the line segment
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn74.gif?pub-status=live)
in the Euclidean space of symmetric matrices to the geodesic between
$\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D654}$
on the Riemannian manifold of positive definite matrices and
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn75.gif?pub-status=live)
Proof. Chapter 6 in Bhatia (Reference Bhatia2015).
Finally, using the affine-invariance property of the Riemannian metric and noting that
$\unicode[STIX]{x1D644}$
commutes with every element of
$\boldsymbol{P}\boldsymbol{o}\boldsymbol{s}_{n}$
, one can prove the following theorem.
Theorem 2. Let
$\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D654}$
be positive definite matrices. There exists a unique geodesic
$\unicode[STIX]{x1D653}\#_{r}\unicode[STIX]{x1D654}$
on the Riemannian manifold of positive definite matrices that joins
$\unicode[STIX]{x1D653}$
and
$\unicode[STIX]{x1D654}$
with the following parametrization:
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn76.gif?pub-status=live)
which is natural in the sense that
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn77.gif?pub-status=live)
for each
$r\in \mathbb{R}$
. Furthermore, we have
![](https://static.cambridge.org/binary/version/id/urn:cambridge.org:id:binary:20180323121557109-0154:S0022112018001180:S0022112018001180_eqn78.gif?pub-status=live)
Proof. Chapter 6 in Bhatia (Reference Bhatia2015).
Appendix B. Details of the numerical evolution of the conformation tensor equations
The conformation tensor field is solved on a temporal grid that is staggered by half a time step with respect to that of the velocity field. We use linear interpolation to transfer between the two temporal grids. The conformation tensor is time-marched using an equal sub-step second-order Runge–Kutta scheme. The advection and stretching terms employ an Adams–Bashforth discretization at each sub-step. Following Dubief et al. (Reference Dubief, Terrapon, White, Shaqfeh, Moin and Lele2005), we ensure that the polymer does not exceed its maximum extensibility by using a semi-implicit approach for the relaxation term in the conformation tensor.
The advection term in the conformation tensor is second-order accurate in space with the exception of a small number points where it is only first-order accurate:
${\sim}0.1$
–0.3 % in the
$x$
and
$y$
directions and
${\sim}3$
–5 % in the
$z$
direction. This change in order is due to the special treatment of the advection term in the conformation tensor that is needed to avoid numerical issues which lead to the conformation tensor losing positive-definiteness. The special treatment is an adaptation of the slope-limiting approach of Vaithianathan et al. (Reference Vaithianathan, Robert, Brasseur and Collins2006) (see also Dallas et al. (Reference Dallas, Vassilicos and Hewitt2010) for an implementation). The approach of Vaithianathan et al. (Reference Vaithianathan, Robert, Brasseur and Collins2006) requires the evaluation of three schemes, utilizing forward, backward and centred stencils, to approximate the flux at the boundaries of each computational cell, at each time step. The scheme that maximizes the eigenvalues of the conformation tensor at the boundaries is then chosen. We modify this approach by evaluating the following schemes in order, and choosing the first one that yields a positive definite tensor at the boundaries: (i) centred stencil, (ii) upwind biased stencil, (iii) downwind biased stencil and (iv) first-order approximation that equates the conformation tensor at the cell centre and the boundary. Our approach is mathematically consistent with Vaithianathan et al. (Reference Vaithianathan, Robert, Brasseur and Collins2006) but is computationally more efficient and defers to the unbiased approximation when possible. In practice, we find that case (i) is sufficient for the vast majority of the points in the domain (
${>}90$
%).
The evolution equation for the conformation tensor has no associated boundary conditions since it is hyperbolic. The conformation tensor at the walls, where no derivatives of the conformation tensor are needed, is then explicitly marched in time.