BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Charles R. Doering

Charles R. Doering appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynmath-phmath.MPmath.OCmath.APmath.DSnlin.CDastro-ph.SRcond-mat.stat-mechmath.FAmath.NAMolecular Networksphysics.ao-phphysics.comp-phphysics.soc-phPopulations and EvolutionQuantitative MethodsSocial and Information Networks

Catalog footprint

What is connected

20works
18topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

20 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
preprint2016arXiv

Bounds for convection between rough boundaries

We consider Rayleigh-Bénard convection in a layer of fluid between rough no-slip boundaries where the top and bottom boundary heights are functions of the horizontal coordinates with square-integrable gradients. We use the background method to derive an upper bound on mean heat flux across the layer for all admissible boundary geometries. This flux, normalized by the temperature difference between the boundaries, can grow with the Rayleigh number ($Ra$) no faster than ${\cal O}(Ra^{1/2})$ as $Ra \rightarrow \infty$. Our analysis yields a family of similar bounds, depending on how various estimates are tuned, but every version depends explicitly on the boundary geometry. In one version the coefficient of the ${\cal O}(Ra^{1/2})$ leading term is $0.242 + 2.925\Vert\nabla h\Vert^2$, where $\Vert\nabla h\Vert^2$ is the mean squared magnitude of the boundary height gradients. Application to a particular geometry is illustrated for sinusoidal boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Gene expression dynamics with stochastic bursts: exact results for a coarse-grained model

We present a theoretical framework to analyze the dynamics of gene expression with stochastic bursts. Beginning with an individual-based model which fully accounts for the messenger RNA (mRNA) and protein populations, we propose a novel expansion of the master equation for the joint process. The resulting coarse-grained model reduces the dimensionality of the system, describing only the protein population while fully accounting for the effects of discrete and fluctuating mRNA population. Closed form expressions for the stationary distribution of the protein population and mean first-passage times of the coarse-grained model are derived and large-scale Monte Carlo simulations show that the analysis accurately describes the individual-based process accounting for mRNA population, in contrast to the failure of commonly proposed diffusion-type models.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Wall to Wall Optimal Transport

The calculus of variations is employed to find steady divergence-free velocity fields that maximize transport of a tracer between two parallel walls held at fixed concentration for one of two constraints on flow strength: a fixed value of the kinetic energy or a fixed value of the enstrophy. The optimizing flows consist of an array of (convection) cells of a particular aspect ratio Gamma. We solve the nonlinear Euler-Lagrange equations analytically for weak flows and numerically (and via matched asymptotic analysis in the fixed energy case) for strong flows. We report the results in terms of the Nusselt number Nu, a dimensionless measure of the tracer transport, as a function of the Peclet number Pe, a dimensionless measure of the energy or enstrophy of the flow. For both constraints the maximum transport Nu_{MAX}(Pe) is realized in cells of decreasing aspect ratio Gamma_{opt}(Pe) as Pe increases. For the fixed energy problem, Nu_{MAX} \sim Pe and Gamma_{opt} \sim Pe^{-1/2}, while for the fixed enstrophy scenario, Nu_{MAX} \sim Pe^{10/17} and Gamma_{opt} \sim Pe^{-0.36}. We also interpret our results in the context of certain buoyancy-driven Rayleigh-Benard convection problems that satisfy one of the two intensity constraints, enabling us to investigate how the transport scalings compare with upper bounds on Nu expressed as a function of the Rayleigh number \Ra. For steady convection in porous media, corresponding to the fixed energy problem, we find Nu_{MAX} \sim \Ra and Gamma_{opt} \sim Ra^{-1/2}$, while for steady convection in a pure fluid layer between free-slip isothermal walls, corresponding to fixed enstrophy transport, Nu_{MAX} \sim Ra^{5/12} and Gamma_{opt} \sim Ra^{-1/4}.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bounds on Surface Stress Driven Shear Flow

The background method is adapted to derive rigorous limits on surface speeds and bulk energy dissipation for shear stress driven flow in two and three dimensional channels. By-products of the analysis are nonlinear energy stability results for plane Couette flow with a shear stress boundary condition: when the applied stress is gauged by a dimensionless Grashoff number $Gr$, the critical $Gr$ for energy stability is 139.5 in two dimensions, and 51.73 in three dimensions. We derive upper bounds on the friction (a.k.a. dissipation) coefficient $C_f = τ/\bar{u}^2$, where $τ$ is the applied shear stress and $\bar{u}$ is the mean velocity of the fluid at the surface, for flows at higher $Gr$ including developed turbulence: $C_f le 1/32$ in two dimensions and $C_f \le 1/8$ in three dimensions. This analysis rigorously justifies previously computed numerical estimates.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The small-world effect is a modern phenomenon

The "small-world effect" is the observation that one can find a short chain of acquaintances, often of no more than a handful of individuals, connecting almost any two people on the planet. It is often expressed in the language of networks, where it is equivalent to the statement that most pairs of individuals are connected by a short path through the acquaintance network. Although the small-world effect is well-established empirically for contemporary social networks, we argue here that it is a relatively recent phenomenon, arising only in the last few hundred years: for most of mankind's tenure on Earth the social world was large, with most pairs of individuals connected by relatively long chains of acquaintances, if at all. Our conclusions are based on observations about the spread of diseases, which travel over contact networks between individuals and whose dynamics can give us clues to the structure of those networks even when direct network measurements are not available. As an example we consider the spread of the Black Death in 14th-century Europe, which is known to have traveled across the continent in well-defined waves of infection over the course of several years. Using established epidemiological models, we show that such wave-like behavior can occur only if contacts between individuals living far apart are exponentially rare. We further show that if long-distance contacts are exponentially rare, then the shortest chain of contacts between distant individuals is on average a long one. The observation of the wave-like spread of a disease like the Black Death thus implies a network without the small-world effect.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

"Ultimate state" of two-dimensional Rayleigh-Benard convection between free-slip fixed temperature boundaries

Rigorous upper limits on the vertical heat transport in two dimensional Rayleigh-Benard convection between stress-free isothermal boundaries are derived from the Boussinesq approximation of the Navier-Stokes equations. The Nusselt number Nu is bounded in terms of the Rayleigh number Ra according to $Nu \leq 0.2295 Ra^{5/12}$ uniformly in the Prandtl number Pr. This Nusselt number scaling challenges some theoretical arguments regarding the asymptotic high Rayleigh number heat transport by turbulent convection.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Internal heating driven convection at infinite Prandtl number

We derive an improved rigorous bound on the space and time averaged temperature $<T>$ of an infinite Prandtl number Boussinesq fluid contained between isothermal no-slip boundaries thermally driven by uniform internal heating. A novel approach is used wherein a singular stable stratification is introduced as a perturbation to a non-singular background profile, yielding the estimate $<T>\geq 0.419[R\log(R)]^{-1/4}$ where $R$ is the heat Rayleigh number. The analysis relies on a generalized Hardy-Rellich inequality that is proved in the appendix.

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

Destabilizing Taylor-Couette flow with suction

We consider the effect of radial fluid injection and suction on Taylor-Couette flow. Injection at the outer cylinder and suction at the inner cylinder generally results in a linearly unstable steady spiralling flow, even for cylindrical shears that are linearly stable in the absence of a radial flux. We study nonlinear aspects of the unstable motions with the energy stability method. Our results, though specialized, may have implications for drag reduction by suction, accretion in astrophysical disks, and perhaps even in the flow in the earth's polar vortex.

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'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

Optimal stirring strategies for passive scalar mixing

We address the challenge of optimal incompressible stirring to mix an initially inhomogeneous distribution of passive tracers. As a quantitative measure of mixing we adopt the $H^{-1}$ norm of the scalar fluctuation field, equivalent to the (square-root of the) variance of a low-pass filtered image of the tracer concentration field. First we establish that this is a useful gauge even in the absence of molecular diffusion: its vanishing as $t --> \infty$ is evidence of the stirring flow's mixing properties in the sense of ergodic theory. Then we derive absolute limits on the total amount of mixing, as a function of time, on a periodic spatial domain with a prescribed instantaneous stirring energy or stirring power budget. We subsequently determine the flow field that instantaneously maximizes the decay of this mixing measure---when such a flow exists. When no such `steepest descent' flow exists (a possible but non-generic situation) we determine the flow that maximizes the growth rate of the $H^{-1}$ norm's decay rate. This local-in-time optimal stirring strategy is implemented numerically on a benchmark problem and compared to an optimal control approach using a restricted set of flows. Some significant challenges for analysis are outlined.

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' 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
preprint2010arXiv

The Mixing efficiency of open flows

Mixing in open incompressible flows is studied in a model problem with inhomogeneous passive scalar injection on an inlet boundary. As a measure of the efficiency of stirring, the bulk scalar concentration variance is bounded and the bound is shown to be sharp at low Peclet number. Although no specific flow saturating the bound at high Peclet number is produced here, the estimate is conjectured to be approached for flows possessing sufficiently sustained chaotic regions.

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