Variational bounds on the energy dissipation rate in body-forced shear flow

A new variational problem for upper bounds on the rate of energy dissipation in body-forced shear flows is formulated by including a balance parameter in the derivation from the Navier-Stokes equations. The resulting min-max problem is investigated computationally, producing new estimates that quantitatively improve previously obtained rigorous bounds. The results are compared with data from direct numerical simulations.Comment: 15 pages, 7 figure

Variational bound on energy dissipation in turbulent shear flow

We present numerical solutions to the extended Doering-Constantin variational principle for upper bounds on the energy dissipation rate in plane Couette flow, bridging the entire range from low to asymptotically high Reynolds numbers. Our variational bound exhibits structure, namely a pronounced minimum at intermediate Reynolds numbers, and recovers the Busse bound in the asymptotic regime. The most notable feature is a bifurcation of the minimizing wavenumbers, giving rise to simple scaling of the optimized variational parameters, and of the upper bound, with the Reynolds number.Comment: 4 pages, RevTeX, 5 postscript figures are available as one .tar.gz file from [email protected]

Transcription factor search for a DNA promoter in a three-states model

To ensure fast gene activation, Transcription Factors (TF) use a mechanism known as facilitated diffusion to find their DNA promoter site. Here we analyze such a process where a TF alternates between 3D and 1D diffusion. In the latter (TF bound to the DNA), the TF further switches between a fast translocation state dominated by interaction with the DNA backbone, and a slow examination state where interaction with DNA base pairs is predominant. We derive a new formula for the mean search time, and show that it is faster and less sensitive to the binding energy fluctuations compared to the case of a single sliding state. We find that for an optimal search, the time spent bound to the DNA is larger compared to the 3D time in the nucleus, in agreement with recent experimental data. Our results further suggest that modifying switching via phosphorylation or methylation of the TF or the DNA can efficiently regulate transcription.Comment: 4 pages, 3 figure

Variational bound on energy dissipation in plane Couette flow

We present numerical solutions to the extended Doering-Constantin variational principle for upper bounds on the energy dissipation rate in turbulent plane Couette flow. Using the compound matrix technique in order to reformulate this principle's spectral constraint, we derive a system of equations that is amenable to numerical treatment in the entire range from low to asymptotically high Reynolds numbers. Our variational bound exhibits a minimum at intermediate Reynolds numbers, and reproduces the Busse bound in the asymptotic regime. As a consequence of a bifurcation of the minimizing wavenumbers, there exist two length scales that determine the optimal upper bound: the effective width of the variational profile's boundary segments, and the extension of their flat interior part.Comment: 22 pages, RevTeX, 11 postscript figures are available as one uuencoded .tar.gz file from [email protected]

Rocket measurement of the secondary electron spectrum in an aurora

Secondary electron emission measurement in aurora using Aerobee rocket vehicl

Time-stepping approach for solving upper-bound problems: Application to two-dimensional Rayleigh-Benard convection

An alternative computational procedure for numerically solving a class of variational problems arising from rigorous upper-bound analysis of forced-dissipative infinite-dimensional nonlinear dynamical systems, including the Navier-Stokes and Oberbeck-Boussinesq equations, is analyzed and applied to Rayleigh-Benard convection. A proof that the only steady state to which this numerical algorithm can converge is the required global optimal of the relevant variational problem is given for three canonical flow configurations. In contrast with most other numerical schemes for computing the optimal bounds on transported quantities (e.g., heat or momentum) within the "background field" variational framework, which employ variants of Newton's method and hence require very accurate initial iterates, the new computational method is easy to implement and, crucially, does not require numerical continuation. The algorithm is used to determine the optimal background-method bound on the heat transport enhancement factor, i.e., the Nusselt number (Nu), as a function of the Rayleigh number (Ra), Prandtl number (Pr), and domain aspect ratio L in two-dimensional Rayleigh-Benard convection between stress-free isothermal boundaries (Rayleigh's original 1916 model of convection). The result of the computation is significant because analyses, laboratory experiments, and numerical simulations have suggested a range of exponents alpha and beta in the presumed Nu similar to (PrRa beta)-Ra-alpha scaling relation. The computations clearly show that for Ra <= 10(10) at fixed L = 2 root 2, Nu <= 0.106Pr(0)Ra(5/12), which indicates that molecular transport cannot generally be neglected in the "ultimate" high-Ra regime.NSF DMS-0928098 DMS-1515161 DMS-0927587 PHY-1205219Simons FoundationNSFONRInstitute for Computational Engineering and Sciences (ICES

Diffusion-Limited Aggregation Processes with 3-Particle Elementary Reactions

A diffusion-limited aggregation process, in which clusters coalesce by means of 3-particle reaction, A+A+A->A, is investigated. In one dimension we give a heuristic argument that predicts logarithmic corrections to the mean-field asymptotic behavior for the concentration of clusters of mass $m$ at time $t$, $c(m,t)~m^{-1/2}(log(t)/t)^{3/4}$, for $1 << m << \sqrt{t/log(t)}$. The total concentration of clusters, $c(t)$, decays as $c(t)~\sqrt{log(t)/t}$ at $t --> \infty$. We also investigate the problem with a localized steady source of monomers and find that the steady-state concentration $c(r)$ scales as $r^{-1}(log(r))^{1/2}$, $r^{-1}$, and $r^{-1}(log(r))^{-1/2}$, respectively, for the spatial dimension $d$ equal to 1, 2, and 3. The total number of clusters, $N(t)$, grows with time as $(log(t))^{3/2}$, $t^{1/2}$, and $t(log(t))^{-1/2}$ for $d$ = 1, 2, and 3. Furthermore, in three dimensions we obtain an asymptotic solution for the steady state cluster-mass distribution: $c(m,r) \sim r^{-1}(log(r))^{-1}\Phi(z)$, with the scaling function $\Phi(z)=z^{-1/2}\exp(-z)$ and the scaling variable $z ~ m/\sqrt{log(r)}$.Comment: 12 pages, plain Te

Temperature control of a liquid helium cooled Eigler-style scanning tunneling microscope

A procedure for operating an Eigler-style, low temperature scanning tunneling microscope (STM) at variable temperatures has been developed. A critical exchange gas pressure regime was found to allow for controlled variation of the STM temperature while it is encapsulated in a liquid helium Dewar. The sensitivity of various parameters to the ability to generate stable variable temperatures above 4 K is discussed

Dynamics and thermodynamics of axisymmetric flows: I. Theory

We develop new variational principles to study stability and equilibrium of axisymmetric flows. We show that there is an infinite number of steady state solutions. We show that these steady states maximize a (non-universal) $H$-function. We derive relaxation equations which can be used as numerical algorithm to construct stable stationary solutions of axisymmetric flows. In a second part, we develop a thermodynamical approach to the equilibrium states at some fixed coarse-grained scale. We show that the resulting distribution can be divided in a universal part coming from the conservation of robust invariants and one non-universal determined by the initial conditions through the fragile invariants (for freely evolving systems) or by a prior distribution encoding non-ideal effects such as viscosity, small-scale forcing and dissipation (for forced systems). Finally, we derive a parameterization of inviscid mixing to describe the dynamics of the system at the coarse-grained scale

Strong friction limit in quantum mechanics: the Quantum Smoluchowski equation

For a quantum system coupled to a heat bath environment the strong friction limit is studied starting from the exact path integral formulation. Generalizing the classical Smoluchowski limit to low temperatures a time evolution equation for the position distribution is derived and the strong role of quantum fluctuations in this limit is revealed.Comment: 4 pages, PRL in pres