BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

J. D. Gibbon

J. D. Gibbon contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
nlin.CDmath-phmath.MPphysics.flu-dynmath.APphysics.ao-phnlin.PSnlin.SIphysics.hist-ph

Trust snapshot

Quick read

Trust 21 - Emerging
20works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2021arXiv

A correspondence between the multifractal model of turbulence and the Navier-Stokes equations

We study a correspondence between the multifractal model of turbulence and the Navier-Stokes equations in $d$ spatial dimensions by comparing their respective dissipation length scales. In Kolmogorov's 1941 theory the key parameter $h$, which is an exponent in the Navier-Stokes invariance scaling, is fixed at $h=1/3$ but is allowed a spectrum of values in multifractal theory. Taking into account all derivatives of the Navier-Stokes equations, it is found that for this correspondence to hold the multifractal spectrum $C(h)$ must be bounded from below such that $C(h) \geq 1-3h$, which is consistent with the four-fifths law. Moreover, $h$ must also be bounded from below such that $h \geq (1-d)/3$. When $d=3$ the allowed range of $h$ is given by $h \geq -2/3$ thereby bounding $h$ away from $h=-1$. The implications of this are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2017arXiv

Martin David Kruskal: a biographical memoir

Martin David Kruskal was one of the most versatile theoretical physicists of his generation and is distinguished for his enduring work in several different areas, most notably plasma physics, a memorable detour into relativity, and his pioneering work in nonlinear waves. In the latter, together with Norman Zabusky, he invented the concept of the soliton and, with others, developed its application to classes of partial differential equations of physical significance.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Depletion of Nonlinearity in Magnetohydrodynamic Turbulence: Insights from Analysis and Simulations

We build on recent developments in the study of fluid turbulence [Gibbon \textit{et al.} Nonlinearity 27, 2605 (2014)] to define suitably scaled, order-$m$ moments, $D_m^{\pm}$, of $ω^\pm= ω\pm j$, where $ω$ and $j$ are, respectively, the vorticity and current density in three-dimensional magnetohydrodynamics (MHD). We show by mathematical analysis, for unit magnetic Prandtl number $P_M$, how these moments can be used to identify three possible regimes for solutions of the MHD equations; these regimes are specified by inequalities for $D_m^{\pm}$ and $D_1^{\pm}$. We then compare our mathematical results with those from our direct numerical simulations (DNSs) and thus demonstrate that 3D MHD turbulence is like its fluid-turbulence counterpart insofar as all solutions, which we have investigated, remain in \textit{only one of these regimes}; this regime has depleted nonlinearity. We examine the implications of our results for the exponents $q^{\pm}$ that characterize the power-law dependences of the energy spectra $\mathcal{E}^{\pm}(k)$ on the wave number $k$, in the inertial range of scales. We also comment on (a) the generalization of our results to the case $P_M \neq 1$ and (b) the relation between $D_m^{\pm}$ and the order-$m$ moments of gradients of hydrodynamic fields, which are used in characterizing intermittency in turbulent flows.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Nonlinear effects in buoyancy-driven variable density turbulence

We consider the time-dependence of a hierarchy of scaled $L^{2m}$-norms $D_{m,ω}$ and $D_{m,θ}$ of the vorticity $\boldsymbol ω = \boldsymbol{\nabla} \times {\mathbf u}$ and the density gradient $\boldsymbol{\nabla} θ$, where $θ=\log (ρ^*/ρ^*_0)$, in a buoyancy-driven turbulent flow as simulated by \cite{LR2007}. $ρ^*({\mathbf x},\,t) $ is the composition density of a mixture of two incompressible miscible fluids with fluid densities $ρ^*_2 > ρ^*_1$ and $ρ^*_{0}$ is a reference normalisation density. Using data from the publicly available Johns Hopkins Turbulence Database we present evidence that the $L^{2}$-spatial average of the density gradient $\boldsymbol{\nabla} θ$ can reach extremely large values, even in flows with low Atwood number $At = (ρ^*_{2} - ρ^*_{1})/(ρ^*_{2} + ρ^*_{1}) = 0.05$, implying that very strong mixing of the density field at small scales can arise in buoyancy-driven turbulence. This large growth raises the possibility that the density gradient $\boldsymbol{\nabla} θ$ might blow up in a finite time.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

High-low frequency slaving and regularity issues in the 3D Navier-Stokes equations

The old idea that an infinite dimensional dynamical system may have its high modes or frequencies slaved to low modes or frequencies is re-visited in the context of the $3D$ Navier-Stokes equations. A set of dimensionless frequencies $\{\tildeΩ_{m}(t)\}$ are used which are based on $L^{2m}$-norms of the vorticity. To avoid using derivatives a closure is assumed that suggests that the $\tildeΩ_{m}$ ($m>1$) are slaved to $\tildeΩ_{1}$ (the global enstrophy) in the form $\tildeΩ_{m} = \tildeΩ_{1}\mathcal{F}_{m}(\tildeΩ_{1})$. This is shaped by the constraint of two Hölder inequalities and a time average from which emerges a form for $\mathcal{F}_{m}$ which has been observed in previous numerical Navier-Stokes and MHD simulations. When written as a phase plane in a scaled form, this relation is parametrized by a set of functions $1 \leq λ_{m}(τ) \leq 4$, where curves of constant $λ_{m}$ form the boundaries between tongue-shaped regions. In regions where $2.5 \leq λ_{m} \leq 4$ and $1 \leq λ_{m} \leq 2$ the Navier-Stokes equations are shown to be regular\,: numerical simulations appear to lie in the latter region. Only in the central region $2 < λ_{m} < 2.5$ has no proof of regularity been found.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Stretching and folding processes in the 3D Euler and Navier-Stokes equations

Stretching and folding dynamics in the incompressible, stratified 3D Euler and Navier-Stokes equations are reviewed in the context of the vector $\bdB = \nabla q\times\nablaθ$ where $q=\bom\cdot\nablaθ$. The variable $θ$ is the temperature and $\bdB$ satisfies $\partial_{t}\bdB = \mbox{curl}\,(\bu\times\bdB)$. These ideas are then discussed in the context of the full compressible Navier-Stokes equations where $q$ takes the two forms $q = \bom\cdot\nablaρ$ and $q = \bom\cdot\nabla(\lnρ)$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Vorticity moments in four numerical simulations of the 3D Navier-Stokes equations

The issue of intermittency in numerical solutions of the 3D Navier-Stokes equations on a periodic box $[0,\,L]^{3}$ is addressed through four sets of numerical simulations that calculate a new set of variables defined by $D_{m}(t) = \left(\varpi_{0}^{-1}Ω_{m}\right)^{α_{m}}$ for $1 \leq m \leq \infty$ where $α_{m}= \frac{2m}{4m-3}$ and $\left[Ω_{m}(t)\right]^{2m} = L^{-3}\I |\bom|^{2m}dV$ with $\varpi_{0} = νL^{-2}$. All four simulations unexpectedly show that the $D_{m}$ are ordered for $m = 1\,,...,\,9$ such that $D_{m+1} < D_{m}$. Moreover, the $D_{m}$ squeeze together such that $D_{m+1}/D_{m}\nearrow 1$ as $m$ increases. The first simulation is of very anisotropic decaying turbulence\,; the second and third are of decaying isotropic turbulence from random initial conditions and forced isotropic turbulence at constant Grashof number respectively\,; the fourth is of very high Reynolds number forced, stationary, isotropic turbulence at up to resolutions of $4096^{3}$.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Conditional regularity of solutions of the three dimensional Navier-Stokes equations and implications for intermittency

Two unusual time-integral conditional regularity results are presented for the three-dimensional Navier-Stokes equations. The ideas are based on $L^{2m}$-norms of the vorticity, denoted by $Ω_{m}(t)$, and particularly on $D_{m} = Ω_{m}^{α_{m}}$, where $α_{m} = 2m/(4m-3)$ for $m\geq 1$. The first result, more appropriate for the unforced case, can be stated simply : if there exists an $1\leq m < \infty$ for which the integral condition is satisfied ($Z_{m}=D_{m+1}/D_{m}$) $$ \int_{0}^{t}\ln (\frac{1 + Z_{m}}{c_{4,m}}) dτ\geq 0$$ then no singularity can occur on $[0, t]$. The constant $c_{4,m} \searrow 2$ for large $m$. Secondly, for the forced case, by imposing a critical \textit{lower} bound on $\int_{0}^{t}D_{m} dτ$, no singularity can occur in $D_{m}(t)$ for \textit{large} initial data. Movement across this critical lower bound shows how solutions can behave intermittently, in analogy with a relaxation oscillator. Potential singularities that drive $\int_{0}^{t}D_{m} dτ$ over this critical value can be ruled out whereas other types cannot.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamics of scaled norms of vorticity for the three-dimensional Navier-Stokes and Euler equations

A series of numerical experiments is suggested for the three-dimensional Navier-Stokes and Euler equations on a periodic domain based on a set of $L^{2m}$-norms of vorticity $Ω_{m}$ for $m\geq 1$. These are scaled to form the dimensionless sequence $D_{m}= (\varpi_{0}^{-1}Ω_{m})^{α_{m}}$ where $\varpi_{0}$ is a constant frequency and $α_{m} = 2m/(4m-3)$. A numerically testable Navier-Stokes regularity criterion comes from comparing the relative magnitudes of $D_{m}$ and $D_{m+1}$ while another is furnished by imposing a critical lower bound on $\int_{0}^{t}D_{m}\,dτ$. The behaviour of the $D_{m}$ is also important in the Euler case in suggesting a method by which possible singular behaviour might also be tested.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Quasi-conservation laws for compressible 3D Navier-Stokes flow

We formulate the quasi-Lagrangian fluid transport dynamics of mass density $ρ$ and the projection $q=\bom\cdot\nablaρ$ of the vorticity $\bom$ onto the density gradient, as determined by the 3D compressible Navier-Stokes equations for an ideal gas, although the results apply for an arbitrary equation of state. It turns out that the quasi-Lagrangian transport of $q$ cannot cross a level set of $ρ$. That is, in this formulation, level sets of $ρ$ (isopychnals) are impermeable to the transport of the projection $q$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A hierarchy of length scales for weak solutions of the three-dimensional Navier-Stokes equations

Moments of the vorticity are used to define and estimate a hierarchy of time-averaged inverse length scales for weak solutions of the three-dimensional, incompressible Navier-Stokes equations on a periodic box. The estimate for the smallest of these inverse scales coincides with the inverse Kolmogorov length but thereafter the exponents of the Reynolds number rise rapidly for correspondingly higher moments. The implications of these results for the computational resolution of small scale vortical structures are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Enstrophy bounds and the range of space-time scales in the hydrostatic primitive equations

The hydrostatic primitive equations (HPE) form the basis of most numerical weather, climate and global ocean circulation models. Analytical (not statistical) methods are used to find a scaling proportional to $(Nu\,Ra\,Re)^{1/4}$ for the range of horizontal spatial sizes in HPE solutions, which is much broader than currently achievable computationally. The range of scales for the HPE is determined from an analytical bound on the time-averaged enstrophy of the horizontal circulation. This bound allows the formation of very small spatial scales, whose existence would excite unphysically large linear oscillation frequencies and gravity wave speeds.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Extreme events in solutions of hydrostatic and non-hydrostatic climate models

Initially this paper reviews the mathematical issues surrounding the hydrostatic (HPE) and non-hydrostatic (NPE) primitive equations that have been used extensively in numerical weather prediction and climate modelling. Cao and Titi (2005, 2007) have provided a new impetus to this by proving existence and uniqueness of solutions of viscous HPE on a cylinder with Neumann-like boundary conditions on the top and bottom. In contrast, the regularity of solutions of NPE remains an open question. With this HPE regularity result in mind, the second issue examined in this paper is whether extreme events are allowed to arise spontaneously in their solutions. Such events could include, for example, the sudden appearance and disappearance of locally intense fronts that do not involve deep convection. Analytical methods are used to show that for viscous HPE, the creation of small-scale structures is allowed locally in space and time at sizes that scale inversely with the Reynolds number.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stretching and folding diagnostics in solutions of the three-dimensional Euler and Navier-Stokes equations

Two possible diagnostics of stretching and folding (S&F) in fluid flows are discussed, based on the dynamics of the gradient of potential vorticity ($q = \bom\cdot\nablaθ$) associated with solutions of the three-dimensional Euler and Navier-Stokes equations. The vector $\bdB = \nabla q \times \nablaθ$ satisfies the same type of stretching and folding equation as that for the vorticity field $\bom $ in the incompressible Euler equations (Gibbon & Holm, 2010). The quantity $θ$ may be chosen as the potential temperature for the stratified, rotating Euler/Navier-Stokes equations, or it may play the role of a seeded passive scalar for the Euler equations alone. The first discussion of these S&F-flow diagnostics concerns a numerical test for Euler codes and also includes a connection with the two-dimensional surface quasi-geostrophic equations. The second S&F-flow diagnostic concerns the evolution of the Lamb vector $\bsD = \bom\times\bu$, which is the nonlinearity for Euler&#39;s equations apart from the pressure. The curl of the Lamb vector ($\boldsymbol{\varpi} := \bsD$) turns out to possess similar stretching and folding properties to that of the $\bdB$-vector.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The dynamics of the gradient of potential vorticity

The transport of the potential vorticity gradient $\bnabla{q}$ along surfaces of constant temperature $θ$ is investigated for the stratified Euler, Navier-Stokes and hydrostatic primitive equations of the oceans and atmosphere using the divergenceless flux vector $\bdB = \bnabla Q(q)\times\bnablaθ$, for any smooth function $Q(q)$. The flux $\bdB$ is shown to satisfy $$ \partial_t\bdB - {curl} (\bU\times\bdB) = - \bnabla\big[qQ&#39;(q) {div} \bU\big]\times\bnablaθ, $$ where $\bU$ is a formal transport velocity of PV flux. While the left hand side of this expression is reminiscent of the frozen-in magnetic field flux in magnetohydrodynamics, the non-zero right hand side means that $\bdB$ is not frozen into the flow of $\bU$ when ${div} \bU \neq 0$. The result may apply to measurements of potential vorticity and potential temperature at the tropopause.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The gradient of potential vorticity, quaternions and an orthonormal frame for fluid particles

The gradient of potential vorticity (PV) is an important quantity because of the way PV (denoted as $q$) tends to accumulate locally in the oceans and atmospheres. Recent analysis by the authors has shown that the vector quantity $\bdB = \bnabla q\times \bnablaθ$ for the three-dimensional incompressible rotating Euler equations evolves according to the same stretching equation as for $\bom$ the vorticity and $\bB$, the magnetic field in magnetohydrodynamics (MHD). The $\bdB$-vector therefore acts like the vorticity $\bom$ in Euler&#39;s equations and the $\bB$-field in MHD. For example, it allows various analogies, such as stretching dynamics, helicity, superhelicity and cross helicity. In addition, using quaternionic analysis, the dynamics of the $\bdB$-vector naturally allow the construction of an orthonormal frame attached to fluid particles\,; this is designated as a quaternion frame. The alignment dynamics of this frame are particularly relevant to the three-axis rotations that particles undergo as they traverse regions of a flow when the PV gradient $\bnabla q$ is large.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Regularity and singularity in solutions of the three-dimensional Navier-Stokes equations

Higher moments of the vorticity field $Ω_{m}(t)$ in the form of $L^{2m}$-norms ($1 \leq m < \infty$) are used to explore the regularity problem for solutions of the three-dimensional incompressible Navier-Stokes equations on the domain $[0, L]^{3}_{per}$. It is found that the set of quantities $$ D_{m}(t) = Ω_{m}^{α_{m}} ,\qquad\qquadα_{m} = \frac{2m}{4m-3}, $$ provide a natural scaling in the problem resulting in a bounded set of time averages $<D_{m}>_{T}$ on a finite interval of time $[0, T]$. The behaviour of $D_{m+1}/D_{m}$ is studied on what are called `good&#39; and `bad&#39; intervals of $[0, T]$ which are interspersed with junction points (neutral) $τ_{i}$. For large but finite values of $m$ with large initial data \big($Ω_{m}(0) \leq \varpi_{0}O(\Gr^{4})$\big), it is found that there is an upper bound $$ Ω_{m} \leq c_{av}^{2}\varpi_{0}\Gr^{4} ,\qquad\varpi_{0} = νL^{-2}, $$ which is punctured by infinitesimal gaps or windows in the vertical walls between the good/bad intervals through which solutions may escape. While this result is consistent with that of Leray \cite{Leray} and Scheffer \cite{Scheff76}, this estimate for $Ω_{m}$ corresponds to a length scale well below the validity of the Navier-Stokes equations.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Estimating intermittency in three-dimensional Navier-Stokes turbulence

The issue of why computational resolution in Navier-Stokes turbulence is so hard to achieve is addressed. It is shown that Navier-Stokes solutions can potentially behave differently in two distinct regions of space-time $\mathbb{R}^{\pm}$ where $\mathbb{R}^{-}$ is comprised of a union of disjoint space-time `anomalies&#39;. Large values of $|\nabla\bom|$ dominate $\mathbb{R}^{-}$, which is consistent with the formation of vortex sheets or tightly-coiled filaments. The local number of degrees of freedom $\mathcal{N}^{\pm}$ needed to resolve the regions in $\mathbb{R}^{\pm}$ satisfies $$\mathcal{N}^{\pm}(\bx, t)\lessgtr c_{\pm}\mathcal{R}_{u}^{3}$$ where $\mathcal{R}_{u} = uL/ν$ is a Reynolds number dependent on the local velocity field $u(\bx, t)$.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Ortho-normal quaternion frames, Lagrangian evolution equations and the three-dimensional Euler equations

More than 150 years after their invention by Hamilton, quaternions are now widely used in the aerospace and computer animation industries to track the paths of moving objects undergoing three-axis rotations. It is shown here that they provide a natural way of selecting an appropriate ortho-normal frame -- designated the quaternion-frame -- for a particle in a Lagrangian flow, and of obtaining the equations for its dynamics. How these ideas can be applied to the three-dimensional Euler fluid equations is then considered. This work has a bearing on the issue of whether the Euler equations develop a singularity in a finite time. Some of the literature on this topic is reviewed, which includes both the Beale-Kato-Majda theorem and associated work on the direction of vorticity by both Constantin, Fefferman & Majda and Deng, Hou and Yu. It is then shown how the quaternion formulation provides a further direction of vorticity result using the Hessian of the pressure.

0 citations0 reviewsNo code linkedOpen access