Source author record

C. P. Caulfield

C. P. Caulfield appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.flu-dyn physics.geo-ph cond-mat.soft nlin.CD physics.ao-ph

Catalog footprint

What is connected

7works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

preprint2022arXiv

Implications of inertial subrange scaling for stably stratified mixing

The effects of turbulent dynamic range on scalar mixing in stably stratified turbulence are investigated by an adaptation of the theoretical passive scalar modelling arguments of Beguier et al. (1978) and demonstrated statistically using direct numerical simulations of statistically stationary homogeneous stratified and sheared turbulence (SHSST). By analysis of inertial and inertial-convective subrange scaling, we show that the relationship between active scalar and turbulence time scales is predicted by the ratio of the Kolmogorov and Oboukhov-Corrsin constants provided there is sufficient scale separation for inertial and inertial-convective subrange scalings to be valid. With this analysis, we show that the turbulent mixing coefficient, $Γ\equiv χ/ε$, that is, within this context defined to be the ratio of available potential energy ($E_p$) and turbulent kinetic energy ($E_k$) dissipation rates, can be estimated by $E_p,E_k$ and a universal constant provided a Reynolds number is sufficiently high, observed here at $Re_b \equiv ε/ νN^2 \gtrapprox 300$ where $ν$ is the kinematic viscosity and $N$ is the characteristic buoyancy frequency. We propose a model for diapycnal mixing with robust theoretical parametrisation and asymptotic behaviour in this high-$Re_b$ limit.

preprint2021arXiv

Layering and vertical transport in sheared double diffusive convection in the diffusive regime

A sequence of two and three-dimensional simulations is conducted for the double diffusive convection (DDC) flows in the diffusive regime subjected to an imposed shear. The flow is confined between two horizontal plates which are maintained at different constant temperature, salinity, and different velocity, thus setting up a shear across the flow. The lower plate is fixed at higher temperature and salinity, while the overall (unperturbed) density gradient is statically stable. For a wide range of control parameters, and for sufficiently strong perturbation of the conductive initial state, we find that staircase-like structures spontaneously develop, with relatively well-mixed layers separated by sharp interfaces of enhanced scalar gradient. Such staircases appear to be robust even in the presence of strong shear over very long times, although we typically observe early time coarsening of the number of observed layers. For the same set of control parameters, different asymptotic layered states, with markedly different vertical scalar fluxes, can arise for different initial perturbation structures. The imposed shear does significantly spatio-temporally modify the vertical transport of the various scalars. The flux ratio (i.e., the ratio between the density fluxes due to the total (convective and diffusive) salt flux and the total heat flux) is found, at steady state, to be essentially equal to the square root of the ratio of the salt diffusivity to the thermal diffusivity, consistently with the physical model originally proposed by Linden and Shirtcliffe (1978) and the variational arguments presented by Stern (1982) for unsheared double diffusive convection.

preprint2020arXiv

Confronting Grand Challenges in Environmental Fluid Mechanics

Environmental fluid mechanics underlies a wealth of natural, industrial and, by extension, societal challenges. In the coming decades, as we strive towards a more sustainable planet, there are a wide range of grand challenge problems that need to be tackled, ranging from fundamental advances in understanding and modeling of stratified turbulence and consequent mixing, to applied studies of pollution transport in the ocean, atmosphere and urban environments. A workshop was organized in the Les Houches School of Physics in France in January 2019 with the objective of gathering leading figures in the field to produce a road map for the scientific community. Five subject areas were addressed: multiphase flow, stratified flow, ocean transport, atmospheric and urban transport, and weather and climate prediction. This article summarizes the discussions and outcomes of the meeting, with the intent of providing a resource for the community going forward.

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.

preprint2014arXiv

Bulldozing of granular material

We investigate the bulldozing motion of a granular sandpile driven forwards by a vertical plate. The problem is set up in the laboratory by emplacing the pile on a table rotating underneath a stationary plate; the continual circulation of the bulldozed material allows the dynamics to be explored over relatively long times, and the variation of the velocity with radius permits one to explore the dependence on bulldozing speed within a single experiment. We measure the time-dependent surface shape of the dune for a range of rotation rates, initial volumes and radial positions, for four granular materials, ranging from glass spheres to irregularly shaped sand. The evolution of the dune can be separated into two phases: a rapid initial adjustment to a state of quasi-steady avalanching perpendicular to the blade, followed by a much slower phase of lateral spreading and radial migration. The quasi-steady avalanching sets up a well-defined perpendicular profile with a nearly constant slope. This profile can be scaled by the depth against the bulldozer to collapse data from different times, radial positions and experiments onto common master curves that are characteristic of the granular material and depend on the local Froude number. The lateral profile of the dune along the face of the bulldozer varies more gradually with radial position, and evolves by slow lateral spreading. The spreading is asymmetrical, with the inward progress of the dune eventually arrested and its bulk migrating to larger radii. A one-dimensional depth-averaged model recovers the nearly linear perpendicular profile of the dune, but does not capture the finer nonlinear details of the master curves. A two-dimensional version of the model leads to an advection-diffusion equation that reproduces the lateral spreading and radial migration.

preprint2012arXiv

A variational framework for flow optimization using semi-norm constraints

When considering a general system of equations describing the space-time evolution (flow) of one or several variables, the problem of the optimization over a finite period of time of a measure of the state variable at the final time is a problem of great interest in many fields. Methods already exist in order to solve this kind of optimization problem, but sometimes fail when the constraint bounding the state vector at the initial time is not a norm, meaning that some part of the state vector remains unbounded and might cause the optimization procedure to diverge. In order to regularize this problem, we propose a general method which extends the existing optimization framework in a self-consistent manner. We first derive this framework extension, and then apply it to a problem of interest. Our demonstration problem considers the transient stability properties of a one-dimensional (in space) averaged turbulent model with a space- and time-dependent model "turbulent viscosity". We believe this work has a lot of potential applications in the fluid dynamics domain for problems in which we want to control the influence of separate components of the state vector in the optimization process.

preprint2011arXiv

Variational identification of minimal seeds to trigger transition in plane Couette flow

A variational formulation incorporating the full Navier-Stokes equations is used to identify initial perturbations with finite kinetic energy E_{0} which generate the largest gain in perturbation kinetic energy (across all possible time intervals) for plane Couette flow. Two different representative flow geometries are chosen corresponding to those used previously by Butler & Farrell (1992) and Monokrousos et al. (2011). In the former (smaller geometry) case as E_{0} increases from 0, we find an optimal which is a smooth nonlinear continuation of the well-known linear result at $E_{0} = 0$. At $E_{0} = E_{c}$, however, completely unrelated states are uncovered which trigger turbulence and our algorithm consequently fails to converge. As $E_{0} \rightarrow E^{+}_{c}, we find good evidence that the turbulence triggering initial conditions approach a 'minimal seed' which corresponds to the state of lowest energy on the laminar-turbulent basin boundary or 'edge'. This situation is repeated in the Monokrousos et al. (2011) (larger) geometry albeit with one notable new feature - the appearance of a nonlinear optimal (as found recently in pipe flow by Pringle & Kerswell (2010) and boundary layer flow by Cherubini et al. (2010)) at finite $E_{0} < E_{c}$ which has a very different structure to the linear optimal. Again the minimal seed at $E_{0} = E_{c}$ does not resemble the linear or now the nonlinear optimal. Our results support the first of two conjectures recently posed by Pringle et al. (2011) but contradict the second. Importantly, their prediction that the form of the functional optimised is not important for identifying E_{c} providing heightened values are produced by turbulent flows is confirmed: we find the what looks to be the same E_{c} and minimal seed using energy gain as opposed to total dissipation in the Monokrousos et al. (2011) geometry.