Glimpses of Flow Development and Degradation During Type B Drag Reduction by Aqueous Solutions of Polyacrylamide B1120

Flow development and degradation during Type B turbulent drag reduction by 0.10 to 10 wppm solutions of a partially-hydrolysed polyacrylamide B1120 of MW \(=\) 18x10⁶ was studied in a smooth pipe of ID \(=\) 4.60 mm and L/D \(=\) 210 at Reynolds numbers from 10000 to 80000 and wall shear stresses Tw from 8 to 600 Pa. B1120 solutions exhibited facets of a Type B ladder, including segments roughly parallel to, but displaced upward from, the P-K line; those that attained asymptotic maximum drag reduction at low Re√ f but departed downwards into the polymeric regime at a higher retro-onset Re√ f; and segments at MDR for all Re√ f. Axial flow enhancement profiles of S\(^{\prime }\) vs L/D reflected a superposition of flow development and polymer degradation effects, the former increasing and the latter diminishing S\(^{\prime }\) with increasing distance downstream. Solutions that induced normalized flow enhancements S\(^{\prime }\)/S\(^{\prime }_{\mathrm {m}} <\) 0.4 developed akin to solvent, with L_e,p/D \(=\) L_e,n/D \(<\) 42.3, while those at maximum drag reduction showed entrance lengths L_e,m/D \(\sim \) 117, roughly 3 times the solvent L_e,n/D. Degradation kinetics were inferred by first detecting a falloff point (Re√f^∧, S\(^{{\prime }\wedge }\)), of maximum observed flow enhancement, for each polymer solution. A plot of S\(^{{\prime }\wedge }\)vs C revealed S\(^{{\prime }\wedge }\)linear in C at low C, with lower bound [S\(^{\prime }\)] \(=\) 5.0 wppm^??1, and S\(^{{\prime }\wedge }\) independent of C at high C, with upper bound S\(^{\prime }_{\mathrm {m}} =\) 15.9. The ratio S\(^{\prime }\)/S\(^{{\prime }\wedge }\) in any pipe section was interpreted to be the undegraded fraction of original polymer therein. Semi-log plots of (S\(^{\prime }\)/S\(^{{\prime }\wedge }\)) at a section vs transit time from pipe entrance thereto revealed first order kinetics, from which apparent degradation rate constants k_deg s^??1 and entrance severities ?ln(S\(^{\prime }\)/S\(^{{\prime }\wedge }\))₀ were extracted. At constant C, k_deg increased linearly with increasing wall shear stress Tw, and at constant Tw, k_deg was independent of C, providing a B1120 degradation modulus (k_deg/Tw) \(=\) (0.012 \(\pm \) 0.001) (Pa s)^??1 for 8 \(<\) Tw Pa \(<\) 600, 0.30 \(<\) C wppm \(<\) 10. Entrance severities were negligible below a threshold Twe \(\sim \) 30 Pa and increased linearly with increasing Tw for Tw \(>\) Twe. The foregoing methods were applied to Type A drag reduction by 0.10 to 10 wppm solutions of a polyethyleneoxide PEO P309, MW \(=\) 11x10⁶, in a smooth pipe of ID \(=\) 7.77 mm and L/D \(=\) 220 at Re from 4000 to 115000. P309 solutions that induced S\(^{\prime }\)/S\(^{\prime }_{\mathrm {m}} <\) 0.4 developed akin to solvent, with L_e,p/D \(=\) L_e,n/D \(<\) 23, while those at MDR had entrance lengths L_e,m/D \(\sim \) 93, roughly 4 times the solvent L_e,n/D. P309 solutions described a Type A fan distorted by polymer degradation. A typical trajectory departed the P-K line at an onset point Re√ f* followed by ascending and descending polymeric regime segments separated by a falloff point Re√f^∧, of maximum flow enhancement; for all P309 solutions, onset Re√ f* = 550 \(\pm \) 100 and falloff Re√f^∧ = 2550 \(\pm \) 250, the interval between them delineating Type A drag reduction unaffected by degradation. A plot of falloff S\(^{{\prime }\wedge }\) vs C for PEO P309 solutions bore a striking resemblance to the analogous S\(^{{\prime }\wedge }\) vs C plot for solutions of PAMH B1120, indicating that the initial Type A drag reduction by P309 after onset at Re√ f* had evolved to Type B drag reduction by falloff at Re√f^∧. Presuming that Type B behaviour persisted past falloff permitted inference of P309 degradation kinetics; k_deg was found to increase linearly with increasing Tw at constant C and was independent of C at constant Tw, providing a P309 degradation modulus (k_deg/Tw) \(=\) (0.011 \(\pm \) 0.002) (Pa s)^??1 for 4 \(<\) Tw Pa \(<\) 400, 0.10 \(<\) C wppm < 5.0. Comparisons between the present degradation kinetics and previous literature showed (k_deg/Tw) data from laboratory pipes of D \(\sim \) 0.01 m to lie on a simple extension of (k_deg/Tw) data from pipelines of D \(\sim \) 0.1 m and 1.0 m, along a power-law relation (k_deg/Tw) \(=\) 10^??5.4.D^??1.6. Intrinsic slips derived from PAMH B1120 and PEO P309-at-falloff experiments were compared with previous examples from Type B drag reduction by polymers with vinylic and glycosidic backbones, showing: (i) For a given polymer, [S\(^{\prime }\)] was independent of Re√ f and pipe ID, implying insensitivity to both micro- and macro-scales of turbulence; and (ii) [S\(^{\prime }\)] increased linearly with increasing polymer chain contour length Lc, the proportionality constant \(\beta =\) 0.053 \(\pm \) 0.036 enabling estimation of flow enhancement S\(^{\prime } =\) C.Lc.β for all Type B drag reduction by polymers.     

Temporal Homogenization of Linear ODEs,with Applications to Parametric Super-Resonance and Energy Harvest

We consider the temporal homogenization of linear ODEs of the form \({\dot{x}=Ax+\epsilon P(t)x+f(t)}\), where P(t) is periodic and \({\epsilon}\) is small. Using a 2-scale expansion approach, we obtain the long-time approximation \({x(t)\approx {\rm exp}(At) \left( \Omega(t)+\int_0^t {\rm exp}(-A \tau) f(\tau) {\rm d}\tau \right)}\), where \({\Omega}\) solves the cell problem \({\dot{\Omega}=\epsilon B \Omega + \epsilon F(t)}\) with an effective matrix B and an explicitly-known F(t). We provide necessary and sufficient conditions for the accuracy of the approximation (over a \({{\mathcal{O}}(\epsilon^{-1})}\) time-scale), and show how B can be computed (at a cost independent of \({\epsilon}\)). As a direct application, we investigate the possibility of using RLC circuits to harvest the energy contained in small scale oscillations of ambient electromagnetic fields (such as Schumann resonances). Although a RLC circuit parametrically coupled to the field may achieve such energy extraction via parametric resonance, its resistance R needs to be smaller than a threshold \({\kappa}\) proportional to the fluctuations of the field, thereby limiting practical applications. We show that if n RLC circuits are appropriately coupled via mutual capacitances or inductances, then energy extraction can be achieved when the resistance of each circuit is smaller than \({n\kappa}\). Hence, if the resistance of each circuit has a non-zero fixed value, energy extraction can be made possible through the coupling of a sufficiently large number n of circuits (\({n\approx 1000}\) for the first mode of Schumann resonances and contemporary values of capacitances, inductances and resistances). The theory is also applied to the control of the oscillation amplitude of a (damped) oscillator.     

Smooth positon solutions of the focusing modified Korteweg–de Vries equation

The n-fold Darboux transformation \(T_{n}\) of the focusing real modified Korteweg–de Vries (mKdV) equation is expressed in terms of the determinant representation. Using this representation, the n-soliton solutions of the mKdV equation are also expressed by determinants whose elements consist of the eigenvalues \(\lambda _{j}\) and the corresponding eigenfunctions of the associated Lax equation. The nonsingular n-positon solutions of the focusing mKdV equation are obtained in the special limit \(\lambda _{j}\rightarrow \lambda _{1}\), from the corresponding n-soliton solutions and by using the associated higher-order Taylor expansion. Furthermore, the decomposition method of the n-positon solution into n single-soliton solutions, the trajectories, and the corresponding “phase shifts” of the multi-positons are also investigated.     

Exponential Decay of the Vorticity in the Steady-State Flow of a Viscous Liquid Past a Rotating Body

We construct a Sobolev homeomorphism in dimension \({n \geqq 4,\,f \in W^{1,1}((0, 1)^n,\mathbb{R}^n)}\) such that \({J_f = {\rm det} Df > 0}\) on a set of positive measure and J _f  < 0 on a set of positive measure. It follows that there are no diffeomorphisms (or piecewise affine homeomorphisms) f _k such that \({f_k\to f}\) in \({W^{1,1}_{\rm loc}}\).     

Topological Pressure of Generic Points for Amenable Group Actions

Let (X, G) be a G-action topological dynamical system (t.d.s. for short), where G is a countably infinite discrete amenable group. In this paper, we study the topological pressure of the sets of generic points. We show that when the system satisfies the almost specification property, for any G-invariant measure \(\mu \) and any continuous map \(\varphi \),

$$\begin{aligned} P\left( X_{\mu },\varphi ,\{F_n\}\right) = h_{\mu }(X)+\int \varphi d\mu , \end{aligned}$$

where \(\{F_n\}\) is a Følner sequence, \(X_{\mu }\) is the set of generic points of \(\mu \) with respect to (w.r.t. for short) \(\{F_n\}\), \(P(X_{\mu },\varphi ,\{F_n\})\) is the topological pressure of \(X_{\mu }\) for \(\varphi \) w.r.t. \(\{F_n\}\) and \(h_{\mu }(X)\) is the measure-theoretic entropy.

     

Steady Nearly Incompressible Vector Fields in Two-Dimension: Chain Rule and Renormalization

Given bounded vector field \({b : {\mathbb{R}^{d}} \to {\mathbb{R}^{d}}}\), scalar field \({u : {\mathbb{R}^{d}} \to {\mathbb{R}}}\), and a smooth function \({\beta : {\mathbb{R}} \to {\mathbb{R}}}\), we study the characterization of the distribution \({{\rm div}(\beta(u)b)}\) in terms of div b and div(ub). In the case of BV vector fields b (and under some further assumptions), such characterization was obtained by L. Ambrosio, C. De Lellis and J. Malý, up to an error term which is a measure concentrated on the so-called tangential set of b. We answer some questions posed in their paper concerning the properties of this term. In particular, we construct a nearly incompressible BV vector field b and a bounded function u for which this term is nonzero. For steady nearly incompressible vector fields b (and under some further assumptions), in the case when d = 2, we provide complete characterization of div(\({\beta(u)b}\)) in terms of div b and div(ub). Our approach relies on the structure of level sets of Lipschitz functions on \({{\mathbb{R}^{2}}}\) obtained by G. Alberti, S. Bianchini and G. Crippa. Extending our technique, we obtain new sufficient conditions when any bounded weak solution u of \({\partial_t u + b \cdot \nabla u=0}\) is renormalized, that is when it also solves \({\partial_t \beta(u) + b \cdot \nabla \beta(u)=0}\) for any smooth function \({\beta \colon{\mathbb{R}} \to {\mathbb{R}}}\). As a consequence, we obtain new a uniqueness result for this equation.     

Existence of $$C^{k}$$-Invariant Foliations for Lorenz-Type Maps

Under conditions similar to those in Shashkov and Shil’nikov (Differ Uravn 30(4):586–595, 732, 1994) we show that a \(C^{k+1}\) Lorenz-type map T has a \(C^{k}\) codimension one foliation which is invariant under the action of T. This allows us to associate T to a \(C^{k}\) one-dimensional transformation.     

Global Attractors of Sixth Order PDEs Describing the Faceting of Growing Surfaces

A spatially two-dimensional sixth order PDE describing the evolution of a growing crystalline surface h(x, y, t) that undergoes faceting is considered with periodic boundary conditions, as well as its reduced one-dimensional version. These equations are expressed in terms of the slopes \(u_1=h_{x}\) and \(u_2=h_y\) to establish the existence of global, connected attractors for both equations. Since unique solutions are guaranteed for initial conditions in \(\dot{H}^2_{per}\), we consider the solution operator \(S(t): \dot{H}^2_{per} \rightarrow \dot{H}^2_{per}\), to gain our results. We prove the necessary continuity, dissipation and compactness properties.     

Stability of Traveling Waves of Nonlinear Schrödinger Equation with Nonzero Condition at Infinity

We study effective elastic behavior of the incompatibly prestrained thin plates, where the prestrain is independent of thickness and uniform through the plate’s thickness h. We model such plates as three-dimensional elastic bodies with a prescribed pointwise stress-free state characterized by a Riemannian metric G, and seek the limiting behavior as \({h \to 0}\). We first establish that when the energy per volume scales as the second power of h, the resulting \({\Gamma}\) -limit is a Kirchhoff-type bending theory. We then show the somewhat surprising result that there exist non-immersible metrics G for whom the infimum energy (per volume) scales smaller than h². This implies that the minimizing sequence of deformations carries nontrivial residual three-dimensional energy but it has zero bending energy as seen from the limit Kirchhoff theory perspective. Another implication is that other asymptotic scenarios are valid in appropriate smaller scaling regimes of energy. We characterize the metrics G with the above property, showing that the zero bending energy in the Kirchhoff limit occurs if and only if the Riemann curvatures R₁₂₁₃, R₁₂₂₃ and R₁₂₁₂ of G vanish identically. We illustrate our findings with examples; of particular interest is an example where \({G_{2 \times 2}}\), the two-dimensional restriction of G, is flat but the plate still exhibits the energy scaling of the Föppl–von Kármán type. Finally, we apply these results to a model of nematic glass, including a characterization of the condition when the metric is immersible, for \({G = Id_{3} + \gamma n \otimes n}\) given in terms of the inhomogeneous unit director field distribution \({ n \in \mathbb{R}^3}\).     

On the interplay between fast reaction and slow diffusion in the concrete carbonation process: a matched-asymptotics approach

A matched-asymptotics approach is proposed to show the occurrence of two distinct characteristic length scales in the carbonation process. The separation of these scales arises due to the strong competition between reaction and diffusion effects. We show that for sufficiently large times τ the width of the carbonated region is proportional to \(\sqrt{\tau}\), while the width of the reaction front is proportional to \(\tau^{\frac{p-1}{2(p+1)}}\) for carbonation-reaction rates with a power law structure like k[CO₂] ^p [Ca(OH)₂] ^q , where k>0 and p,q>1 and identify the proportionality coefficient asymptotically. We emphasize the occurrence of a water barrier in the reaction zone which may hinder the penetration of CO₂ by locally filling with water air parts of the pores. This non-linear effect may be one of the causes why a purely linear extrapolation of accelerated carbonation test results to natural carbonation settings is (even theoretically) not reasonable. Finally, we compare our asymptotic penetration law against measured penetration depths from Bune (Zum Karbonatisierungsbedingten Verlust der Dauerhaftigkeit von Außenbauteilen aus Stahlbeton, 1994). The novelty consists in the fact that the factor multiplying \(\sqrt{\tau}\) is now identified asymptotically by solving a non-linear system of ordinary differential equations, and hence, fitting arguments are not necessary to estimate its size. We offer an alternative to the (asymptotic) \(\sqrt{\tau}\) expression of the carbonation-front position obtained in Papadakis et al. (AIChE J. 35:1639, 1989).     

Local Fixed Point Indices of Iterations of Planar Maps

Let \({f: U\rightarrow {\mathbb R}^2}\) be a continuous map, where U is an open subset of \({{\mathbb R}^2}\). We consider a fixed point p of f which is neither a sink nor a source and such that {p} is an isolated invariant set. Under these assumption we prove, using Conley index methods and Nielsen theory, that the sequence of fixed point indices of iterations \({\{{\rm ind}(f^n,p)\}_{n=1}^\infty}\) is periodic, bounded from above by 1, and has infinitely many non-positive terms, which is a generalization of Le Calvez and Yoccoz theorem (Annals of Math., 146, 241–293 (1997)) onto the class of non-injective maps. We apply our result to study the dynamics of continuous maps on 2-dimensional sphere.     

Influence of hydrodynamic slip on convective transport in flow past a circular cylinder

The presence of a finite tangential velocity on a hydrodynamically slipping surface is known to reduce vorticity production in bluff body flows substantially while at the same time enhancing its convection downstream and into the wake. Here, we investigate the effect of hydrodynamic slippage on the convective heat transfer (scalar transport) from a heated isothermal circular cylinder placed in a uniform cross-flow of an incompressible fluid through analytical and simulation techniques. At low Reynolds (\({\textit{Re}}\ll 1\)) and high Péclet (\({\textit{Pe}}\gg 1\)) numbers, our theoretical analysis based on Oseen and thermal boundary layer equations allows for an explicit determination of the dependence of the thermal transport on the non-dimensional slip length \(l_s\). In this case, the surface-averaged Nusselt number, Nu transitions gradually between the asymptotic limits of \(Nu \sim {\textit{Pe}}^{1/3}\) and \(Nu \sim {\textit{Pe}}^{1/2}\) for no-slip (\(l_s \rightarrow 0\)) and shear-free (\(l_s \rightarrow \infty \)) boundaries, respectively. Boundary layer analysis also shows that the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) holds for a shear-free cylinder surface in the asymptotic limit of \({\textit{Re}}\gg 1\) so that the corresponding heat transfer rate becomes independent of the fluid viscosity. At finite \({\textit{Re}}\), results from our two-dimensional simulations confirm the scaling \(Nu \sim {\textit{Pe}}^{1/2}\) for a shear-free boundary over the range \(0.1 \le {\textit{Re}}\le 10^3\) and \(0.1\le {\textit{Pr}}\le 10\). A gradual transition from the lower asymptotic limit corresponding to a no-slip surface, to the upper limit for a shear-free boundary, with \(l_s\), is observed in both the maximum slip velocity and the Nu. The local time-averaged Nusselt number \(Nu_{\theta }\) for a shear-free surface exceeds the one for a no-slip surface all along the cylinder boundary except over the downstream portion where unsteady separation and flow reversal lead to an appreciable rise in the local heat transfer rates, especially at high \({\textit{Re}}\) and Pr. At a Reynolds number of \(10^3\), the formation of secondary recirculating eddy pairs results in appearance of additional local maxima in \(Nu_{\theta }\) at locations that are in close proximity to the mean secondary stagnation points. As a consequence, Nu exhibits a non-monotonic variation with \(l_s\) increasing initially from its lowermost value for a no-slip surface and then decreasing before rising gradually toward the upper asymptotic limit for a shear-free cylinder. A non-monotonic dependence of the spanwise-averaged Nu on \(l_s\) is observed in three dimensions as well with the three-dimensional wake instabilities that appear at sufficiently low \(l_s\), strongly influencing the convective thermal transport from the cylinder. The analogy between heat transfer and single-component mass transfer implies that our results can directly be applied to determine the dependency of convective mass transfer of a single solute on hydrodynamic slip length in similar configurations through straightforward replacement of Nu and \({\textit{Pr}}\) with Sherwood and Schmidt numbers, respectively.     

Vortex solitons of the (3+1)-dimensional spatially modulated cubic–quintic nonlinear Schrödinger equation with the transverse modulation

Vortex solitons in the spatially modulated cubic–quintic nonlinear media are governed by a (3+1)-dimensional cubic–quintic nonlinear Schrödinger equation with spatially modulated nonlinearity and transverse modulation. Via the variable separation principle with the similarity transformation, we derive two families of vortex soliton solutions in the spatially modulated cubic–quintic nonlinear media. For the disappearing and parabolic transverse modulation, vortex solitons with different configurations are constructed. The similar configurations of vortex solitons exist for the same value of \(l-k\) with the topological charge k and degree number l. Moreover, the number of the inner layer structure of vortex solitons getting rid of the package covering layer is related to \((n-1)/2+1\) with the soliton order number n. For the disappearing transverse modulation, there exist phase azimuthal jumps around their cores of vortex solitons with \(2\pi \) phase change in every jump, and any two jumps one after another realize the change in \(\pi \). For the parabolic transverse modulation, all phases of vortex soliton exist k-jump, and every jump realizes the change in \(2\pi /k\); thus, k-jumps totally realize the azimuthal change in \(2\pi \) around their cores.     

Fine Level Set Structure of Flat Isometric Immersions

A result by Pogorelov asserts that C ¹ isometric immersions u of a bounded domain \({S \subset \mathbb R^2}\) into \({\mathbb {R}^3}\) whose normal takes values in a set of zero area enjoy the following regularity property: the gradient \({f := \nabla u}\) is ‘developable’ in the sense that the nondegenerate level sets of f consist of straight line segments intersecting the boundary of S at both endpoints. Motivated by applications in nonlinear elasticity, we study the level set structure of such f when S is an arbitrary bounded Lipschitz domain. We show that f can be approximated by uniformly bounded maps with a simplified level set structure. We also show that the domain S can be decomposed (up to a controlled remainder) into finitely many subdomains, each of which admits a global line of curvature parametrization.     

Focusing Quantum Many-body Dynamics: The Rigorous Derivation of the 1D Focusing Cubic Nonlinear Schrödinger Equation

We consider the dynamics of N bosons in 1D. We assume that the pair interaction is attractive and given by \({N^{\beta-1}V(N^{\beta}.) where }\) where \({\int V \leqslant 0}\). We develop new techniques in treating the N-body Hamiltonian so that we overcome the difficulties generated by the attractive interaction and establish new energy estimates. We also prove the optimal 1D collapsing estimate which reduces the regularity requirement in the uniqueness argument by half a derivative. We derive rigorously the 1D focusing cubic NLS with a quadratic trap as the \({N \rightarrow \infty}\) limit of the N-body dynamic and hence justify the mean-field limit and prove the propagation of chaos for the focusing quantum many-body system.     

Minimizers of the Landau–de Gennes Energy Around a Spherical Colloid Particle

We consider energy minimizing configurations of a nematic liquid crystal around a spherical colloid particle, in the context of the Landau–de Gennes model. The nematic is assumed to occupy the exterior of a ball B_r0, and satisfy homeotropic weak anchoring at the surface of the colloid and approach a uniform uniaxial state as \({|x|\to\infty}\). We study the minimizers in two different limiting regimes: for balls which are small \({r_0\ll L^{\frac12}}\) compared to the characteristic length scale \({L^{\frac 12}}\), and for large balls, \({r_0\gg L^{\frac12}}\). The relationship between the radius and the anchoring strength W is also relevant. For small balls we obtain a limiting quadrupolar configuration, with a “Saturn ring” defect for relatively strong anchoring, corresponding to an exchange of eigenvalues of the Q-tensor. In the limit of very large balls we obtain an axisymmetric minimizer of the Oseen–Frank energy, and a dipole configuration with exactly one point defect is obtained.     

An Eulerian–Lagrangian Form for the Euler Equations in Sobolev Spaces

In 2000 Constantin showed that the incompressible Euler equations can be written in an “Eulerian–Lagrangian” form which involves the back-to-labels map (the inverse of the trajectory map for each fixed time). In the same paper a local existence result is proved in certain Hölder spaces \({C^{1,\mu}}\). We review the Eulerian–Lagrangian formulation of the equations and prove that given initial data in H ^s for \({n \geq 2}\) and \({s > \frac{n}{2}+1}\), a unique local-in-time solution exists on the n-torus that is continuous into H ^s and C ¹ into H ^s-1. These solutions automatically have C ¹ trajectories. The proof here is direct and does not appeal to results already known about the classical formulation. Moreover, these solutions are regular enough that the classical and Eulerian–Lagrangian formulations are equivalent, therefore what we present amounts to an alternative approach to some of the standard theory.