skip to main content


Title: Circulation and Energy Theorem Preserving Stochastic Fluids
Abstract Smooth solutions of the incompressible Euler equations are characterized by the property that circulation around material loops is conserved. This is the Kelvin theorem. Likewise, smooth solutions of Navier–Stokes are characterized by a generalized Kelvin's theorem, introduced by Constantin–Iyer (2008). In this note, we introduce a class of stochastic fluid equations, whose smooth solutions are characterized by natural extensions of the Kelvin theorems of their deterministic counterparts, which hold along certain noisy flows. These equations are called the stochastic Euler–Poincaré and stochastic Navier–Stokes–Poincaré equations respectively. The stochastic Euler–Poincaré equations were previously derived from a stochastic variational principle by Holm (2015), which we briefly review. Solutions of these equations do not obey pathwise energy conservation/dissipation in general. In contrast, we also discuss a class of stochastic fluid models, solutions of which possess energy theorems but do not, in general, preserve circulation theorems.  more » « less
Award ID(s):
1703997
NSF-PAR ID:
10271649
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Proceedings of the Royal Society of Edinburgh: Section A Mathematics
Volume:
150
Issue:
6
ISSN:
0308-2105
Page Range / eLocation ID:
2776 to 2814
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract We formulate a class of stochastic partial differential equations based on Kelvin’s circulation theorem for ideal fluids. In these models, the velocity field is randomly transported by white-noise vector fields, as well as by its own average over realizations of this noise. We call these systems the Lagrangian averaged stochastic advection by Lie transport (LA SALT) equations. These equations are nonlinear and non-local, in both physical and probability space. Before taking this average, the equations recover the Stochastic Advection by Lie Transport (SALT) fluid equations introduced by Holm (Proc R Soc A 471(2176):20140963, 2015). Remarkably, the introduction of the non-locality in probability space in the form of momentum transported by its own mean velocity gives rise to a closed equation for the expectation field which comprises Navier–Stokes equations with Lie–Laplacian ‘dissipation’. As such, this form of non-locality provides a regularization mechanism. The formalism we develop is closely connected to the stochastic Weber velocity framework of Constantin and Iyer (Commun Pure Appl Math 61(3):330–345, 2008) in the case when the noise correlates are taken to be the constant basis vectors in $$\mathbb {R}^3$$ R 3 and, thus, the Lie–Laplacian reduces to the usual Laplacian. We extend this class of equations to allow for advected quantities to be present and affect the flow through exchange of kinetic and potential energies. The statistics of the solutions for the LA SALT fluid equations are found to be changing dynamically due to an array of intricate correlations among the physical variables. The statistical properties of the LA SALT physical variables propagate as local evolutionary equations which when spatially integrated become dynamical equations for the variances of the fluctuations. Essentially, the LA SALT theory is a non-equilibrium stochastic linear response theory for fluctuations in SALT fluids with advected quantities. 
    more » « less
  2. null (Ed.)
    Prior mathematical work of Constantin & Iyer ( Commun. Pure Appl. Maths , vol. 61, 2008, pp. 330–345; Ann. Appl. Probab. , vol. 21, 2011, pp. 1466–1492) has shown that incompressible Navier–Stokes solutions possess infinitely many stochastic Lagrangian conservation laws for vorticity, backward in time, which generalize the invariants of Cauchy ( Sciences mathématiques et physique , vol. I, 1815, pp. 33–73) for smooth Euler solutions. We reformulate this theory for the case of wall-bounded flows by appealing to the Kuz'min ( Phys. Lett. A , vol. 96, 1983, pp. 88–90)–Oseledets ( Russ. Math. Surv. , vol. 44, 1989, p. 210) representation of Navier–Stokes dynamics, in terms of the vortex-momentum density associated to a continuous distribution of infinitesimal vortex rings. The Constantin–Iyer theory provides an exact representation for vorticity at any interior point as an average over stochastic vorticity contributions transported from the wall. We point out relations of this Lagrangian formulation with the Eulerian theory of Lighthill (Boundary layer theory. In Laminar Boundary Layers (ed. L. Rosenhead), 1963, pp. 46–113)–Morton ( Geophys. Astrophys. Fluid Dyn. , vol. 28, 1984, pp. 277–308) for vorticity generation at solid walls, and also with a statistical result of Taylor ( Proc. R. Soc. Lond. A , vol. 135, 1932, pp. 685–702)–Huggins ( J. Low Temp. Phys. , vol. 96, 1994, pp. 317–346), which connects dissipative drag with organized cross-stream motion of vorticity and which is closely analogous to the ‘Josephson–Anderson relation’ for quantum superfluids. We elaborate a Monte Carlo numerical Lagrangian scheme to calculate the stochastic Cauchy invariants and their statistics, given the Eulerian space–time velocity field. The method is validated using an online database of a turbulent channel-flow simulation (Graham et al. , J. Turbul. , vol. 17, 2016, pp. 181–215), where conservation of the mean Cauchy invariant is verified for two selected buffer-layer events corresponding to an ‘ejection’ and a ‘sweep’. The variances of the stochastic Cauchy invariants grow exponentially backward in time, however, revealing Lagrangian chaos of the stochastic trajectories undergoing both fluid advection and viscous diffusion. 
    more » « less
  3. Abstract In Jayanti and Trivisa (2022 J. Math. Fluid Mech. 24 46), the authors proved the existence of local-in-time weak solutions to a model of superfluidity. The system of governing equations was derived in Pitaevskii (1959 Sov. Phys. JETP 8 282–287) and couples the nonlinear Schrödinger equation and the Navier–Stokes equations. In this article, we prove a weak–strong type uniqueness theorem for these weak solutions. Only some of their regularity properties are used, allowing room for improved existence theorems in the future, with compatible uniqueness results. 
    more » « less
  4. Whether the 3D incompressible Navier–Stokes equations can develop a finite time sin- gularity from smooth initial data is one of the most challenging problems in nonlinear PDEs. In this paper, we present some new numerical evidence that the incompress- ible axisymmetric Navier–Stokes equations with smooth initial data of finite energy seem to develop potentially singular behavior at the origin. This potentially singular behavior is induced by a potential finite time singularity of the 3D Euler equations that we reported in a companion paper published in the same issue, see also Hou (Poten- tial singularity of the 3D Euler equations in the interior domain. arXiv:2107.05870 [math.AP], 2021). We present numerical evidence that the 3D Navier–Stokes equa- tions develop nearly self-similar singular scaling properties with maximum vorticity increased by a factor of 107. We have applied several blow-up criteria to study the potentially singular behavior of the Navier–Stokes equations. The Beale–Kato–Majda blow-up criterion and the blow-up criteria based on the growth of enstrophy and neg- ative pressure seem to imply that the Navier–Stokes equations using our initial data develop a potential finite time singularity. We have also examined the Ladyzhenskaya– Prodi–Serrin regularity criteria (Kiselev and Ladyzhenskaya in Izv Akad Nauk SSSR Ser Mat 21(5):655–690, 1957; Prodi in Ann Math Pura Appl 4(48):173–182, 1959; Serrin in Arch Ration Mech Anal 9:187–191, 1962) that are based on the growth rate of Lqt Lxp norm of the velocity with 3/p + 2/q ≤ 1. Our numerical results for the cases of (p,q) = (4,8), (6,4), (9,3) and (p,q) = (∞,2) provide strong evidence for the potentially singular behavior of the Navier–Stokes equations. The critical case of (p,q) = (3,∞) is more difficult to verify numerically due to the extremely slow growth rate in the L3 norm of the velocity field and the significant contribution from the far field where we have a relatively coarse grid. Our numerical study shows that while the global L3 norm of the velocity grows very slowly, the localized version of the L 3 norm of the velocity experiences rapid dynamic growth relative to the localized L 3 norm of the initial velocity. This provides further evidence for the potentially singular behavior of the Navier–Stokes equations. 
    more » « less
  5. Reducing wall drag in turbulent pipe and channel flows is an issue of great practical importance. In engineering applications, end-functionalized polymer chains are often employed as agents to reduce drag. These are polymers which are floating in the solvent and attach (either by adsorption or through irreversible chemical binding) at one of their chain ends to the substrate (wall). We propose a PDE model to study this setup in the simple setting where the solvent is a viscous incompressible Navier–Stokes fluid occupying the bulk of a smooth domain Ω⊂ℝ𝑑, and the wall-grafted polymer is in the so-called mushroom regime (inter-polymer spacing on the order of the typical polymer length). The microscopic description of the polymer enters into the macroscopic description of the fluid motion through a dynamical boundary condition on the wall-tangential stress of the fluid, something akin to (but distinct from) a history-dependent slip-length. We establish the global well-posedness of strong solutions in two-spatial dimensions and prove that the inviscid limit to the strong Euler solution holds with a rate. Moreover, the wall-friction factor ⟨𝑓⟩ and the global energy dissipation ⟨𝜀⟩ vanish inversely proportional to the Reynolds number 𝐑𝐞. This scaling corresponds to Poiseuille’s law for the friction factor ⟨𝑓⟩∼1/𝐑𝐞 for laminar flow and thereby quantifies drag reduction in our setting. These results are in stark contrast to those available for physical boundaries without polymer additives modeled by, for example, no-slip conditions, where no such results are generally known even in two-dimensions. 
    more » « less