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

Charles R. Doering

Charles R. Doering contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.flu-dynmath.DSmath.OCmath-phmath.APmath.FAmath.MPmath.NAnlin.CDphysics.comp-phPopulations and Evolution

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
10works
0followers
11topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

10 published item(s)

preprint2022arXiv

Heat transport in a hierarchy of reduced-order convection models

Reduced-order models (ROMs) are systems of ordinary differential equations (ODEs) designed to approximate the dynamics of partial differential equations (PDEs). In this work, a distinguished hierarchy of ROMs is constructed for Rayleigh's 1916 model of natural thermal convection. These models are distinguished in the sense that they preserve energy and vorticity balances derived from the governing equations, and each is capable of modeling zonal flow. Various models from the hierarchy are analyzed to determine the maximal heat transport in a given model, measured by the dimensionless Nusselt number, for a given Rayleigh number. Lower bounds on the maximal heat transport are ascertained by computing the Nusselt number among equilibria of the chosen model using numerical continuation. A method known as sum-of-squares optimization is applied to construct upper bounds on the time-averaged Nusselt number. In this case, the sum-of-squares approach involves constructing a polynomial quantity whose global nonnegativity implies the upper bound along all solutions to a chosen ROM. The minimum such bound is determined through a type of convex optimization called semidefinite programming. For the ROMs studied in this work, the Nusselt number is maximized by equilibria whenever the Rayleigh number is sufficiently small. In this range of Rayleigh number, the equilibria maximizing heat transport are those that bifurcate first from the zero state. Analyzing this primary equilibrium branch provides a possible mechanism for the increase in heat transport near the onset of convection.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Steady Rayleigh--Bénard convection between no-slip boundaries

The central open question about Rayleigh--Bénard convection -- buoyancy-driven flow in a fluid layer heated from below and cooled from above -- is how vertical heat flux depends on the imposed temperature gradient in the strongly nonlinear regime where the flows are typically turbulent. The quantitative challenge is to determine how the Nusselt number $Nu$ depends on the Rayleigh number $Ra$ in the $Ra\to\infty$ limit for fluids of fixed finite Prandtl number $Pr$ in fixed spatial domains. Laboratory experiments, numerical simulations, and analysis of Rayleigh's mathematical model have yet to rule out either of the proposed `classical' $Nu \sim Ra^{1/3}$ or `ultimate' $Nu \sim Ra^{1/2}$ asymptotic scaling theories. Among the many solutions of the equations of motion at high $Ra$ are steady convection rolls that are dynamically unstable but share features of the turbulent attractor. We have computed these steady solutions for $Ra$ up to $10^{14}$ with $Pr=1$ and various horizontal periods. By choosing the horizontal period of these rolls at each $Ra$ to maximize $Nu$, we find that steady convection rolls achieve classical asymptotic scaling. Moreover, they transport more heat than turbulent convection in experiments or simulations at comparable parameters. If heat transport in turbulent convection continues to be dominated by heat transport in steady rolls as $Ra\to\infty$, it cannot achieve the ultimate scaling.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Optimal time averages in non-autonomous nonlinear dynamical systems

The auxiliary function method allows computation of extremal long-time averages of functions of dynamical variables in autonomous nonlinear ordinary differential equations via convex optimization. For dynamical systems defined by autonomous polynomial vector fields, it is operationally realized as a semidefinite program utilizing sum of squares technology. In this contribution we review the method and extend it for application to periodically driven non-autonomous nonlinear vector fields involving trigonometric functions of the dynamical variables. The damped driven Duffing oscillator and periodically driven pendulum are presented as examples to illustrate the auxiliary function method's utility.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Wall-to-wall optimal transport in two dimensions

Gradient ascent methods are developed to compute incompressible flows that maximize heat transport between two isothermal no-slip parallel walls. Parameterizing the magnitude of velocity fields by a Péclet number $\text{Pe}$ proportional to their root-mean-square rate-of-strain, the schemes are applied to compute two-dimensional flows optimizing convective enhancement of diffusive heat transfer, i.e., the Nusselt number $\text{Nu}$ up to $\text{Pe} \approx 10^5$. The resulting transport exhibits a change of scaling from $\text{Nu}-1 \sim \text{Pe}^{2}$ for $\text{Pe} < 10$ in the linear regime to $\text{Nu} \sim \text{Pe}^{0.54}$ for $\text{Pe} > 10^3$. Optimal fields are observed to be approximately separable, i.e., products of functions of the wall-parallel and wall-normal coordinates. Analysis employing a separable ansatz yields a conditional upper bound $\lesssim \text{Pe}^{6/11} = \text{Pe}^{0.\overline{54}}$ as $\text{Pe} \rightarrow \infty$ similar to the computationally achieved scaling. Implications for heat transfer in buoyancy-driven Rayleigh-Bénard convection are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stirring up trouble: Multi-scale mixing measures for steady scalar sources

The mixing efficiency of a flow advecting a passive scalar sustained by steady sources and sinks is naturally defined in terms of the suppression of bulk scalar variance in the presence of stirring, relative to the variance in the absence of stirring. These variances can be weighted at various spatial scales, leading to a family of multi-scale mixing measures and efficiencies. We derive a priori estimates on these efficiencies from the advection--diffusion partial differential equation, focusing on a broad class of statistically homogeneous and isotropic incompressible flows. The analysis produces bounds on the mixing efficiencies in terms of the Peclet number, a measure the strength of the stirring relative to molecular diffusion. We show by example that the estimates are sharp for particular source, sink and flow combinations. In general the high-Peclet number behavior of the bounds (scaling exponents as well as prefactors) depends on the structure and smoothness properties of, and length scales in, the scalar source and sink distribution. The fundamental model of the stirring of a monochromatic source/sink combination by the random sine flow is investigated in detail via direct numerical simulation and analysis. The large-scale mixing efficiency follows the upper bound scaling (within a logarithm) at high Peclet number but the intermediate and small-scale efficiencies are qualitatively less than optimal. The Peclet number scaling exponents of the efficiencies observed in the simulations are deduced theoretically from the asymptotic solution of an internal layer problem arising in a quasi-static model.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Demographic Fluctuations versus Spatial Variation in the Competition between Fast and Slow Dispersers

Dispersal is an important strategy that allows organisms to locate and exploit favorable habitats. The question arises: given competition in a spatially heterogeneous landscape, what is the optimal rate of dispersal? Continuous population models predict that a species with a lower dispersal rate always drives a competing species to extinction in the presence of spatial variation of resources. However, the introduction of intrinsic demographic stochasticity can reverse this conclusion. We present a simple model in which competition between the exploitation of resources and stochastic fluctuations leads to victory by either the faster or slower of two species depending on the environmental parameters. A simplified limiting case of the model, analyzed by closing the moment and correlation hierarchy, quantitatively predicts which species will win in the complete model under given parameters of spatial variation and average carrying capacity.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Models and measures of mixing and effective diffusion

Mixing a passive scalar field by stirring can be measured in a variety of ways including tracer particle dispersion, via the flux-gradient relationship, or by suppression of scalar concentration variations in the presence of inhomogeneous sources and sinks. The mixing efficiency or efficacy of a particular flow is often expressed in terms of enhanced diffusivity and quantified as an effective diffusion coefficient. In this work we compare and contrast several notions of effective diffusivity. We thoroughly examine the fundamental case of a steady sinusoidal shear flow mixing a scalar sustained by a steady sinusoidal source-sink distribution to explore apparent quantitative inconsistencies among the measures. Ultimately the conflicts are attributed to the noncommutative asymptotic limits of large P$\acute{\text{e}}$clet number and large length-scale separation. We then propose another approach, a generalization of Batchelor&#39;s 1949 theory of diffusion in homogeneous turbulence, that helps unify the particle dispersion and concentration variance suppression measures.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Symmetric factorization of the conformation tensor in viscoelastic fluid models

The positive definite symmetric polymer conformation tensor possesses a unique symmetric square root that satisfies a closed evolution equation in the Oldroyd-B and FENE-P models of viscoelastic fluid flow. When expressed in terms of the velocity field and the symmetric square root of the conformation tensor, these models&#39; equations of motion formally constitute an evolution in a Hilbert space with a total energy functional that defines a norm. Moreover, this formulation is easily implemented in direct numerical simulations resulting in significant practical advantages in terms of both accuracy and stability.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Extreme vorticity growth in Navier-Stokes turbulence

According to statistical turbulence theory, the ensemble averaged squared vorticity rho_E is expected to grow not faster than drho_E/dt ~ rho_E^{3/2}. Solving a variational problem for maximal bulk enstrophy (E) growth, velocity fields were found for which the growth rate is as large as dE/dt ~ E^3. Using numerical simulations with well resolved small scales and a quasi-Lagrangian advection to track fluid subvolumes with rapidly growing vorticity, we study spatially resolved statistics of vorticity growth. We find that the volume ensemble averaged growth bound is satisfied locally to a remarkable degree of accuracy. Elements with dE/dt ~ E^3 can also be identified, but their growth tends to be replaced by the ensemble-averaged law when the intensities become too large.

0 citations0 reviewsNo code linkedOpen access