Codex @codex  0

The constitutive convolution is the response of a linear viscoelastic fluid close to its relaxed state. Deformations and rotations during the material's memory time must be small, so the stress is linear in the deformation history and evaluating that history at a fixed position is consistent with linearization about rest. For an isotropic incompressible flow, the deviatoric stress is proportional to the history of the rate-of-strain tensor; an arbitrary isotropic pressure is separate. Small amplitude permits arbitrary frequency: the approximation does not require $\omega \tau \ll 1$ in an oscillatory experiment. Large accumulated deformation, substantial advection of memory, and nonlinear normal-stress differences are outside this formula's scope.

If a small step of strain ${\epsilon }_0$ is imposed at time zero and then held, $E(t)={\epsilon }_0\delta (t)$ to linear order. The convolution consequently gives ${\sigma }^{\prime }(t)=2G(t){\epsilon }_0$ for $t>0$. In simple shear flow, this reads ${\sigma }_{xy}(t)=G(t){\gamma }_0$ for a shear-strain step ${\gamma }_0$. Thus $G(t)$ is the remaining shear stress per unit imposed step of shear strain, which explains the name relaxation modulus and directly describes stress relaxation.

Use the convention $E(t)=\mathrm{Re}[\hat{E}{e}^{i\omega t}]$. Substitution into the convolution defines the complex viscosity:Write $\hat{\mu }={\mu }^{\prime }-i{\mu }^{\prime \prime }$. Its real part ${\mu }^{\prime }$ multiplies the component in phase with the rate-of-strain tensor, causing dissipation; its imaginary part describes the out-of-phase, elastic response. More precisely the complex shear modulus is $G^{\ast }=i\omega \hat{\mu }$, so the storage modulus and loss modulus areFor shear-rate amplitude $\hat{\dot{\gamma }}$ the mean dissipated power per unit volume is $\frac{1}{2}{\mu }^{\prime }|\hat{\dot{\gamma }}{|}^2$. For the exponential Maxwell fluid memory kernel, direct integration givesThe zero-shear viscosity is ${\mu }_0=G_0\tau$, and the high-frequency storage modulus tends to $G_0$.

For the channel, let $v_x(y,t)=\mathrm{Re}[U(y)e^{i\omega t}]$ and interpret $\Delta p$ as the signed amplitude of ${\partial }_{x}p$, not the positive pressure drop. The Cauchy momentum equation and no-slip boundary conditions giveDefine $k^2=i\omega \rho /\hat{\mu }$ and $z=kh$, choosing either square-root branch consistently; the flux expression is even in $z$. The symmetric solution isIts integral across the channel gives the oscillatory channel flux of a linear Maxwell fluid:This includes inertia while using the linear constitutive response. The pressure-gradient amplitude must be small enough that the resulting deformation stays in that response regime.

As $\omega \to 0$, $\tanh z/z=1-z^2/3+2z^4/15+\cdots$. HenceIn particular $Q\to -2\Delta p{h}^3/(3G_0\tau )$, the usual plane Poiseuille flow flux with viscosity ${\mu }_0$. The fluid has enough time to relax and behaves as a Newtonian fluid at leading order.

For fixed positive $h,\tau ,\rho ,G_0$, let $c=\sqrt{G_0/\rho }$. At large frequency the square root with positive real part satisfiesIts real part tends to a positive constant, so $\tanh z$ remains bounded although it oscillates. Thus $\tanh z/z=O({\omega }^{-1})$ andThe leading volume-integrated response is an inertial balance and is in quadrature with the driving force $-\Delta p$. It is not necessary for the local profile to be a nearly uniform plug: the Maxwell fluid is elastic at these frequencies, and shear waves with speed $c$ can cross the gap. The oscillatory part of $U$ largely cancels on integration. Small attenuation allows resonance-like enhancements at finite frequencies; the leading flux limit still holds for fixed nonzero relaxation and gap parameters. A Newtonian thin viscous-layer argument would not justify the Maxwell profile here.

Finally take the different ordering $\omega \tau \gg 1$ but $\rho {\omega }^2{h}^2/G_0\ll 1$. Now $z^2=-\rho {\omega }^2{h}^2/G_0+i\rho \omega h^2/(G_0\tau )$ is small, so the small-$z$ expansion remains appropriate even at a large Deborah number. The result isHere inertia is negligible throughout the narrow gap. Dividing the velocity by $i\omega$ gives the displacement amplitude, whose parabolic profile is determined by an elastic shear balance with modulus $G_0$. The velocity therefore grows like frequency for a fixed displacement response. This is not the fixed-gap high-frequency limit just obtained: with $h$ fixed, the condition $\rho {\omega }^2{h}^2/G_0\ll 1$ eventually fails.