On acoustic cavitation of slightly subcritical bubbles

The classical Blake threshold indicates the onset of quasistatic evolution leading to cavitation for gas bubbles in liquids. When the mean pressure in the liquid is reduced to a value below the vapor pressure, the Blake analysis identifies a critical radius which separates quasistatically stable bubbles from those which would cavitate. In this work, we analyze the cavitation threshold for radially symmetric bubbles whose radii are slightly less than the Blake critical radius, in the presence of time-periodic acoustic pressure fields. A distinguished limit equation is derived that predicts the threshold for cavitation for a wide range of liquid viscosities and forcing frequencies. This equation also yields frequency-amplitude response curves. Moreover, for fixed liquid viscosity, our study identifies the frequency that yields the minimal forcing amplitude sufficient to initiate cavitation. Numerical simulations of the full Rayleigh-Plesset equation confirm the accuracy of these predictions. Finally, the implications of these findings for acoustic pressure fields that consist of two frequencies will be discussed.Comment: 14 pages, Presented at APS/DFD conference in Philadelphia 199

Deformations of Gabor Frames

The quantum mechanical harmonic oscillator Hamiltonian generates a one-parameter unitary group W(\theta) in L^2(R) which rotates the time-frequency plane. In particular, W(\pi/2) is the Fourier transform. When W(\theta) is applied to any frame of Gabor wavelets, the result is another such frame with identical frame bounds. Thus each Gabor frame gives rise to a one-parameter family of frames, which we call a deformation of the original. For example, beginning with the usual tight frame F of Gabor wavelets generated by a compactly supported window g(t) and parameterized by a regular lattice in the time-frequency plane, one obtains a family of frames F_\theta generated by the non-compactly supported windows g_\theta=W(theta)g, parameterized by rotated versions of the original lattice. This gives a method for constructing tight frames of Gabor wavelets for which neither the window nor its Fourier transform have compact support. When \theta=\pi/2, we obtain the well-known Gabor frame generated by a window with compactly supported Fourier transform. The family F_\theta therefore interpolates these two familiar examples.Comment: 8 pages in Plain Te

Adiabatic stability under semi-strong interactions: The weakly damped regime

We rigorously derive multi-pulse interaction laws for the semi-strong interactions in a family of singularly-perturbed and weakly-damped reaction-diffusion systems in one space dimension. Most significantly, we show the existence of a manifold of quasi-steady N-pulse solutions and identify a "normal-hyperbolicity" condition which balances the asymptotic weakness of the linear damping against the algebraic evolution rate of the multi-pulses. Our main result is the adiabatic stability of the manifolds subject to this normal hyperbolicity condition. More specifically, the spectrum of the linearization about a fixed N-pulse configuration contains essential spectrum that is asymptotically close to the origin as well as semi-strong eigenvalues which move at leading order as the pulse positions evolve. We characterize the semi-strong eigenvalues in terms of the spectrum of an explicit N by N matrix, and rigorously bound the error between the N-pulse manifold and the evolution of the full system, in a polynomially weighted space, so long as the semi-strong spectrum remains strictly in the left-half complex plane, and the essential spectrum is not too close to the origin

Pitchfork bifurcation along a slow parameter ramp: coherent structures in the critical scaling

We investigate the slow passage through a pitchfork bifurcation in a spatially extended system, when the onset of instability is slowly varying in space. We focus here on the critical parameter scaling, when the instability locus propagates with speed $c\sim \varepsilon^{1/3}$, where $\varepsilon$ is a small parameter that measures the gradient of the parameter ramp. Our results establish how the instability is mediated by a front traveling with the speed of the parameter ramp, and demonstrate scalings for a delay or advance of the instability relative to the bifurcation locus depending on the sign of $c$, that is on the direction of propagation of the parameter ramp through the pitchfork bifurcation. The results also include a generalization of the classical Hastings-McLeod solution of the Painlev\'e-II equation to Painlev\'e-II equations with a drift term

The Dynamics of Hybrid Metabolic-Genetic Oscillators

The synthetic construction of intracellular circuits is frequently hindered by a poor knowledge of appropriate kinetics and precise rate parameters. Here, we use generalized modeling (GM) to study the dynamical behavior of topological models of a family of hybrid metabolic-genetic circuits known as "metabolators." Under mild assumptions on the kinetics, we use GM to analytically prove that all explicit kinetic models which are topologically analogous to one such circuit, the "core metabolator," cannot undergo Hopf bifurcations. Then, we examine more detailed models of the metabolator. Inspired by the experimental observation of a Hopf bifurcation in a synthetically constructed circuit related to the core metabolator, we apply GM to identify the critical components of the synthetically constructed metabolator which must be reintroduced in order to recover the Hopf bifurcation. Next, we study the dynamics of a re-wired version of the core metabolator, dubbed the "reverse" metabolator, and show that it exhibits a substantially richer set of dynamical behaviors, including both local and global oscillations. Prompted by the observation of relaxation oscillations in the reverse metabolator, we study the role that a separation of genetic and metabolic time scales may play in its dynamics, and find that widely separated time scales promote stability in the circuit. Our results illustrate a generic pipeline for vetting the potential success of a potential circuit design, simply by studying the dynamics of the corresponding generalized model