Past exam of the advanced fluid dynamics course of the University of Oxford / 2017 / 2

Take $\mathbf{R}={\mathbf{r}}_2-{\mathbf{r}}_1$ and label each dumbbell once. Locally assume its probability density and number density vary negligibly over a connector length. The oriented connector crossing count follows geometrically: a connector crosses a small plane element precisely when its midpoint lies within a slab of normal thickness $|\mathbf{R}\cdot \mathbf{n}|$. For configurations in $d^{3}R$, the unsigned crossing number per unit area is thereforeThe expression without the absolute value is an oriented crossing count, or an unsigned count restricted to $\mathbf{R}\cdot \mathbf{n}\ge 0$. It cannot be an unrestricted nonnegative number: a connector $\mathbf{R}=-R\mathbf{n}$ gives a negative printed expression despite physically intersecting the plane.

To obtain traction, keep the force orientation as well. For $\mathbf{R}\cdot \mathbf{n}>0$, the spring force on the bead below the plane exerted by the bead above it is $H\mathbf{R}$. For the opposite orientation it is $-H\mathbf{R}$. Multiplying that force by the unsigned crossing count in either case givesIntegrating over all labelled configurations yieldsThere is no extra factor of two: each dumbbell was labelled and counted once, and the orientation-dependent force has already supplied the sign. This is the local continuum stress approximation; resolving appreciable spatial variation across a finite connector would require a nonlocal stress construction.

Let $p_T=n{k}_{B}T$, $C=⟨\mathbf{R}{\mathbf{R}}^T⟩$, and $M=I+\alpha \sigma$, so the single-bead mobility is $M/\zeta$. Assume an incompressible solvent, constant $n,H,T,\zeta$, negligible bead inertia, and a symmetric positive-definite $M$. The last requirement makes the resistance invertible and Brownian diffusion physically admissible. At a given position and time, $M$ depends on the moment $C$ but not on the individual configuration $\mathbf{R}$.

Use the source's gradient convention $K_{ij}={\partial }_i{u}_j$, and write $G=K^T$ for the component-first gradient. Then $\mathbf{u}({\mathbf{r}}_2)-\mathbf{u}({\mathbf{r}}_1)\simeq G\mathbf{R}$. The spring forces on the two beads are ${\mathbf{F}}_1=H\mathbf{R}$ and ${\mathbf{F}}_2=-H\mathbf{R}$. Their overdamped probability currents areHere the Brownian term written as $-k_{B}T\nabla \psi$ is a force weighted by probability density. Its corresponding conditional osmotic force is $-k_{B}T\nabla \log \psi$; treating the given term as an unweighted deterministic force per bead would not produce the stated Fokker–Planck equation.

With midpoint $\mathbf{x}=({\mathbf{r}}_1+{\mathbf{r}}_2)/2$, the derivatives are ${\nabla }_{r_1}=\frac{1}{2}{\nabla }_x-{\nabla }_R$ and ${\nabla }_{r_2}=\frac{1}{2}{\nabla }_x+{\nabla }_R$. Subtracting the bead currents gives the configuration currentThe midpoint current is $\mathbf{u}\psi$ after neglecting the specified spatial Brownian diffusion. Configuration probability conservation therefore givesSince $G\mathbf{R}$ has components $R_j{\partial }_j{u}_i$, this agrees with the source's $\mathbf{R}\cdot \nabla \mathbf{u}$. The factor two in both spring relaxation and configuration diffusivity comes from relative motion of both beads.

The second moment $C$ is the unnormalized conformation tensor. Now multiply by $R_i{R}_j$ and integrate over configuration space. Assume sufficiently decaying tails so every boundary term at infinite separation vanishes. The affine current gives $GC+C{G}^T$; spring drift gives $-(2H/\zeta )(MC+CM)$; two integrations by parts of diffusion give $(4k_{B}T/\zeta )M$. HenceBecause $\sigma =nHC-p_{T}I$, $M$ is a polynomial in $C$ and commutes with it. ThusThis is the Hookean dumbbell conformation equation with anisotropic mobility, an exact second-moment closure within the stated model; no Gaussian approximation to $\psi$ was needed.

The upper-convected derivative of the identity is $I^{\nabla }=-(G+G^T)$. Applying it to the polymer stress therefore givesConsequently the coefficients and the closed stress equation arewhere ${\sigma }^{\nabla }=D_t\sigma -\sigma K-K^T\sigma$ in the printed convention. Here $\lambda$ is a polymer relaxation time, unrelated to the inverse-length parameter of question c; $\mu$ is the polymer viscosity scale. The coefficient $\alpha$ has inverse-stress units. In the dimensional mobility parameter in the Giesekus model, the dimensionless factor is $b=\alpha p_T$, not the dimensional $\alpha$ printed in this problem.

The upper-convected rate removes the affine stretching and rotation of material connectors. If $\dot{F}=GF$, then a tensor transported purely as $C=F{C}_0{F}^T$ has $C^{\nabla }=0$. Under a time-dependent rigid rotation $Q$, the rate transforms as $Q{C}^{\nabla }{Q}^T$: this is material frame indifference, not an ordinary componentwise time derivative. The resulting equation is the Giesekus model in these dimensional conventions.

The isotropic subtraction is fixed by zero polymer extra stress at thermal equilibrium. With no flow and positive mobility, the stationary distribution isso $\boxed{{\sigma }_{\mathrm{eq}}=0}$. The constitutive equation correspondingly has a zero-stress quiescent state and a Newtonian low-rate limit $\sigma \simeq \mu (K+K^T)$. Without subtracting exactly $p_{T}I$, the spring contribution would retain an equilibrium isotropic stress and the equation would acquire corresponding reference terms. Because mobility itself depends on $\sigma$, that reference cannot simply be shifted while retaining the same model. The separate fluid pressure enforcing incompressibility is not included in this polymer extra stress.

For steady uniaxial extension of the Giesekus model, the component-first velocity gradient is diagonal, $G=\mathrm{diag}(\gamma ,-\gamma /2,-\gamma /2)$. For a steady spatially uniform diagonal stress, the material derivative vanishes and the convected terms are $-2\gamma {\sigma }_{11}$ and $+\gamma {\sigma }_{22}$. The off-diagonal equations vanish identically, and symmetry permits equal transverse components. The remaining equations areSet $W=\lambda \gamma \ge 0$, $b=\alpha p_T$ and $s_i={\sigma }_{ii}/p_T$, with $\mu =\lambda p_T$. For $b>0$, the branches continuously connected to zero stress at zero rate areA sufficient nondegenerate range for physical solutions at every nonnegative rate is $0<b<1$. Indeed $M=(1-b)I+bA$ with $A=HC/(k_{B}T)$, so positive conformation implies positive mobility. The longitudinal quadratic equation is nonpositive at $s_1=0$ and has a nonnegative root of a polynomial. The transverse quadratic equation has values $b-1<0$ at $s_2=-1$ and $W\ge 0$ at $s_2=0$, so its larger root of a polynomial lies in $(-1,0]$. Thus the displayed branches satisfy $s_i>-1$, giving $A_i=1+s_i>0$. The other algebraic roots do not provide a physical branch connected to the equilibrium distribution. For constant positive diagonal $A$, the associated stationary configuration distribution is a Gaussian with covariance $C=(k_{B}T/H)A$.

For positive $\alpha$ on the physical branch, the stretching stress remains finite through $W=1/2$. At that rate,At large extension rate,Thus the stress continues to grow, but the uniaxial extensional viscosity ${\eta }_E=({\sigma }_{11}-{\sigma }_{22})/\gamma$ approaches the finite plateau $2\lambda /\alpha$ rather than developing a finite-rate pole. The transverse stress approaches $-p_T$ as its conformation tensor variance becomes small. At small rate the two stresses are $2\mu \gamma$ and $-\mu \gamma$, so ${\eta }_E\to 3\mu$.

When $\alpha =0$, the physical steady branch instead isThe stretching stress diverges as $W\to 1/2$ from below. Above that threshold its algebraic continuation has $1+s_1=1/(1-2W)<0$ and cannot be a conformation covariance. There is no physical finite steady stretching state. The transient equationgives linear growth at $W=1/2$ and exponential growth at larger rates, the Hookean coil-stretch failure of the upper-convected Maxwell limit.

Parameter qualification. Positivity of $\alpha$ alone, without a restriction on $b=\alpha n{k}_{B}T$, does not ensure real physical steady solutions for the entire tensor. For example $b=2$, $W=1$ makes the transverse equation $2s_2^2+2s_2+1=0$, which has no real root. Its diffusion/resistance interpretation is already inadmissible for some conformations. The physically positive-mobility regime stated above avoids that counterexample; the source's positive-$\alpha$ case must not be taken as a guarantee for arbitrary parameters.