Search

We explore the chaotic dynamics of a large one-dimensional lattice of coupled maps with diffusive coupling of varying strength using the covariant Lyapunov vectors (CLVs). Using a lattice of diffusively coupled quadratic maps we quantify the growth of spatial structures in the chaotic dynamics as the strength of diffusion is increased. When the diffusion strength is increased from zero, we find that the leading Lyapunov exponent decreases rapidly from a positive value to zero to yield a small window of periodic dynamics which is then followed by chaotic dynamics. For values of the diffusion strength beyond the window of periodic dynamics, the leading Lyapunov exponent does not vary significantly with the strength of diffusion with the exception of a small variation for the largest diffusion strengths we explore. The Lyapunov spectrum and fractal dimension are described analytically as a function of the diffusion strength using the eigenvalues of the coupling operator. The spatial features of the CLVs are quantified and compared with the eigenvectors of the coupling operator. The chaotic dynamics are composed entirely of physical modes for all of the conditions we explore. The leading CLV is highly localized and the localization decreases with increasing strength of the spatial coupling. The violation of the dominance of Oseledets splitting indicates that the entanglement of pairs of CLVs becomes more significant between neighboring CLVs as the strength of the diffusion is increased., Comment: 33 pages, 15 figures

Even a relatively weak drive force is enough to push a typical nanomechanical resonator into the nonlinear regime. Consequently, nonlinearities are widespread in nanomechanics and determine the critical characteristics of nanoelectromechanical systems (NEMS) resonators. A thorough understanding of the nonlinear dynamics of higher eigenmodes of NEMS resonators would be beneficial for progress, given their use in applications and fundamental studies. Here, we characterize the nonlinearity and the linear dynamic range (LDR) of each eigenmode of two nanomechanical beam resonators with different intrinsic tension values up to eigenmode $n=11$. We find that the modal Duffing constant increases as $n^4$, while the critical amplitude for the onset of nonlinearity decreases as $1/n$. The LDR, determined from the ratio of the critical amplitude to the thermal noise amplitude, increases weakly with $n$. Our findings are consistent with our theory treating the beam as a string, with the nonlinearity emerging from stretching at high amplitudes. These scaling laws, observed in experiments and validated theoretically, can be leveraged for pushing the limits of NEMS-based sensing even further., Comment: 7 pages, 3 figures

Brownian motion imposes a hard limit on the overall precision of a nanomechanical measurement. Here, we present a combined experimental and theoretical study of the Brownian dynamics of a quintessential nanomechanical system, a doubly-clamped nanomechanical beam resonator, in a viscous fluid. Our theoretical approach is based on the fluctuation-dissipation theorem of statistical mechanics: We determine the dissipation from fluid dynamics; we incorporate this dissipation into the proper elastic equation to obtain the equation of motion; the fluctuation-dissipation theorem then directly provides an analytical expression for the position-dependent power spectral density (PSD) of the displacement fluctuations of the beam. We compare our theory to experiments on nanomechanical beams immersed in air and water, and obtain excellent agreement. Within our experimental parameter range, the Brownian force noise driving the nanomechanical beam has a colored PSD due to the ``memory" of the fluid; the force noise remains mode-independent and uncorrelated in space. These conclusions are not only important for nanomechanical sensing but also provide insight into the fluctuations of elastic systems at any length scale.

We numerically study the propagation of reacting fronts in a shallow and horizontal layer of fluid with solutal feedback and in the presence of a thermally driven flow field of counter-rotating convection rolls. We solve the Boussinesq equations along with a reaction-convection-diffusion equation for the concentration field where the products of the nonlinear autocatalytic reaction are less dense than the reactants. For small values of the solutal Rayleigh number the characteristic fluid velocity scales linearly, and the front velocity and mixing length scale quadratically, with increasing solutal Rayleigh number. For small solutal Rayleigh numbers the front geometry is described by a curve that is nearly antisymmetric about the horizontal midplane. For large values of the solutal Rayleigh number the characteristic fluid velocity, the front velocity, and the mixing length exhibit square-root scaling and the front shape collapses onto an asymmetric self-similar curve. In the presence of counter-rotating convection rolls, the mixing length decreases while the front velocity increases. The complexity of the front geometry increases when both the solutal and convective contributions are significant and the dynamics can exhibit chemical oscillations in time for certain parameter values. Lastly, we discuss the spatiotemporal features of the complex fronts that form over a range of solutal and thermal driving.

We explore the high-dimensional chaotic dynamics of the Lorenz-96 model by computing the variation of the fractal dimension with system parameters. The Lorenz-96 model is a continuous in time and discrete in space model first proposed by Edward Lorenz to study fundamental issues regarding the forecasting of spatially extended chaotic systems such as the atmosphere. First, we explore the spatiotemporal chaos limit by increasing the system size while holding the magnitude of the external forcing constant. Second, we explore the strong driving limit by increasing the external forcing while holding the system size fixed. As the system size is increased for small values of the forcing we find dynamical states that alternate between periodic and chaotic dynamics. The windows of chaos are extensive, on average, with relative deviations from extensivity on the order of 20%. For intermediate values of the forcing we find chaotic dynamics for all system sizes past a critical value. The fractal dimension exhibits a maximum deviation from extensivity on the order of 5% for small changes in system size and decreases non-monotonically with increasing system size. The length scale describing the deviations from extensivity and the natural chaotic length scale are approximately equal in support of the suggestion that deviations from extensivity are due to the addition of chaotic degrees of freedom as the system size is increased. As the forcing is increased at constant system size the fractal dimension exhibits a power-law dependence. The power-law behavior is independent of the system size and quantifies the decreasing size of chaotic degrees of freedom with increased forcing which we compare with spatial features of the patterns., Comment: 12 pages, 20 figures

Using a thermodynamic approach based upon the fluctuation-dissipation theorem we quantify the stochastic dynamics of rectangular and V-shaped microscale cantilevers immersed in a viscous fluid. We show that the stochastic cantilever dynamics as measured by the displacement of the cantilever tip or by the angle of the cantilever tip are different. We trace this difference to contributions from the higher modes of the cantilever. We find that contributions from the higher modes are significant in the dynamics of the cantilever tip-angle. For the V-shaped cantilever the resulting flow field is three-dimensional and complex in contrast to what is found for a long and slender rectangular cantilever. Despite this complexity the stochastic dynamics can be predicted using a two-dimensional model with an appropriately chosen length scale. We also quantify the increased fluid dissipation that results as a V-shaped cantilever is brought near a solid planar boundary., Comment: 10 pages, 15 images, corrected equation (8)

We consider the stochastic dynamics of an array of two closely spaced atomic force microscope cantilevers in a viscous fluid for use as a possible biomolecule sensor. The cantilevers are not driven externally, as is common in applications of atomic force microscopy, and we explore the stochastic cantilever dynamics due to the constant buffeting of fluid particles by Brownian motion. The stochastic dynamics of two adjacent cantilevers are correlated due to long range effects of the viscous fluid. Using a recently proposed thermodynamic approach the hydrodynamic correlations are quantified for precise experimental conditions through deterministic numerical simulations. Results are presented for an array of two readily available atomic force microscope cantilevers. It is shown that the force on a cantilever due to the fluid correlations with an adjacent cantilever is more than 3 times smaller than the Brownian force on an individual cantilever. Our results indicate that measurements of the correlations in the displacement of an array of atomic force microscopes can detect piconewton forces with microsecond time resolution., Comment: 7 page article with 11 images submitted to the International Journal of Nonlinear Mechanics

The dynamical behavior of propagating structures, determined from a Karhunen-Lo`eve decomposition, in turbulent pipe flow undergoing reverse transition to laminar flow is investigated. The turbulent flow data is generated by a direct numerical simulation started at a fully turbulent Reynolds number of Re_\tau=150, which is slowly decreased until Re_\tau=95. At this low Reynolds number the high frequency modes decay first, leaving only the decaying streamwise vortices. The flow undergoes a chugging phenomena, where it begins to relaminarize and the mean velocity increases. The remaining propagating modes then destabilize the streamwise vortices, rebuild the energy spectra, and eventually the flow regains its turbulent state. Our results capture three chugging cycles before the flow completely relaminarizes. The high frequency modes present in the outer layer decay first, establishing the importance of the outer region in the self-sustaining mechanism of wall bound turbulence., Comment: 8 pages, 20 figures. First post-review update

The results of a comparative analysis based upon a Karhunen-Lo\`{e}ve expansion of turbulent pipe flow and drag reduced turbulent pipe flow by spanwise wall oscillation are presented. The turbulent flow is generated by a direct numerical simulation at a Reynolds number $Re_\tau = 150$. The spanwise wall oscillation is imposed as a velocity boundary condition with an amplitude of $A^+ = 20$ and a period of $T^+ = 50$. The wall oscillation results in a 27% mean velocity increase when the flow is driven by a constant pressure gradient. The peaks of the Reynolds stress and root-mean-squared velocities shift away from the wall and the Karhunen-Lo\`{e}ve dimension of the turbulent attractor is reduced from 2453 to 102. The coherent vorticity structures are pushed away from the wall into higher speed flow, causing an increase of their advection speed of 34% as determined by a normal speed locus. This increase in advection speed gives the propagating waves less time to interact with the roll modes. This leads to less energy transfer and a shorter lifespan of the propagating structures, and thus less Reynolds stress production which results in drag reduction., Comment: 11 pages, 21 figures. Figure added

The results of an analysis of turbulent pipe flow based on a Karhunen-Lo`eve decomposition are presented. The turbulent flow is generated by a direct numerical simulation of the Navier-Stokes equations using a spectral element algorithm at a Reynolds number Re_\tau=150. This simulation yields a set of basis functions that captures 90% of the energy after 2,453 modes. The eigenfunctions are categorised into two classes and six subclasses based on their wavenumber and coherent vorticity structure. Of the total energy, 81% is in the propagating class, characterised by constant phase speeds; the remaining energy is found in the non propagating subclasses, the shear and roll modes. The four subclasses of the propagating modes are the wall, lift, asymmetric, and ring modes. The wall modes display coherent vorticity structures near the wall, the lift modes display coherent vorticity structures that lift away from the wall, the asymmetric modes break the symmetry about the axis, and the ring modes display rings of coherent vorticity. Together, the propagating modes form a wave packet, as found from a circular normal speed locus. The energy transfer mechanism in the flow is a four step process. The process begins with energy being transferred from mean flow to the shear modes, then to the roll modes. Energy is then transfer ed from the roll modes to the wall modes, and then eventually to the lift modes. The ring and asymmetric modes act as catalysts that aid in this four step energy transfer. Physically, this mechanism shows how the energy in the flow starts at the wall and then propagates into the outer layer., Comment: 28 pages, 20 figures. Updated with reviewer's comments / suggestions

The stochastic dynamics of micron and nanoscale cantilevers immersed in a viscous fluid are quantified. Analytical results are presented for long slender cantilevers driven by Brownian noise. The spectral density of the noise force is not assumed to be white and the frequency dependence is determined from the fluctuation-dissipation theorem. The analytical results are shown to be useful for the micron scale cantilevers that are commonly used in atomic force microscopy. A general thermodynamic approach is developed that is valid for cantilevers of arbitrary geometry as well as for arrays of multiple cantilevers whose stochastic motion is coupled through the fluid. It is shown that the fluctuation-dissipation theorem permits the calculation of stochastic quantities via straightforward deterministic methods. The thermodynamic approach is used with deterministic finite element numerical simulations to quantify the autocorrelation and noise spectrum of cantilever fluctuations for a single micron scale cantilever and the cross-correlations and noise spectra of fluctuations for an array of two experimentally motivated nanoscale cantilevers as a function of cantilever separation. The results are used to quantify the noise reduction possible using correlated measurements with two closely spaced nanoscale cantilevers., Comment: Submitted to Nanotechnology April 26, 2006

The coarsening and wavenumber selection of striped states growing from random initial conditions are studied in a non-relaxational, spatially extended, and far-from-equilibrium system by performing large-scale numerical simulations of Rayleigh-B\'{e}nard convection in a large-aspect-ratio cylindrical domain with experimentally realistic boundaries. We find evidence that various measures of the coarsening dynamics scale in time with different power-law exponents, indicating that multiple length scales are required in describing the time dependent pattern evolution. The translational correlation length scales with time as $t^{0.12}$, the orientational correlation length scales as $t^{0.54}$, and the density of defects scale as $t^{-0.45}$. The final pattern evolves toward the wavenumber where isolated dislocations become motionless, suggesting a possible wavenumber selection mechanism for large-aspect-ratio convection., Comment: 5 pages, 6 figures

The stochastic response of nanoscale oscillators of arbitrary geometry immersed in a viscous fluid is studied. Using the fluctuation-dissipation theorem it is shown that deterministic calculations of the governing fluid and solid equations can be used in a straightforward manner to directly calculate the stochastic response that would be measured in experiment. We use this approach to investigate the fluid coupled motion of single and multiple cantilevers with experimentally motivated geometries., Comment: 5 pages, 5 figures

We describe a numerical procedure to construct a modified velocity field that does not have any mean flow. Using this procedure, we present two results. Firstly, we show that, in the absence of mean flow, spiral defect chaos collapses to a stationary pattern comprising textures of stripes with angular bends. The quenched patterns are characterized by mean wavenumbers that approach those uniquely selected by focus-type singularities, which, in the absence of mean flow, lie at the zig-zag instability boundary. The quenched patterns also have larger correlation lengths and are comprised of rolls with less curvature. Secondly, we describe how mean flow can contribute to the commonly observed phenomenon of rolls terminating perpendicularly into lateral walls. We show that, in the absence of mean flow, rolls begin to terminate into lateral walls at an oblique angle. This obliqueness increases with Rayleigh number., Comment: 14 pages, 19 figures

Rayleigh-B\'{e}nard convection is studied and quantitative comparisons are made, where possible, between theory and experiment by performing numerical simulations of the Boussinesq equations for a variety of experimentally realistic situations. Rectangular and cylindrical geometries of varying aspect ratios for experimental boundary conditions, including fins and spatial ramps in plate separation, are examined with particular attention paid to the role of the mean flow. A small cylindrical convection layer bounded laterally either by a rigid wall, fin, or a ramp is investigated and our results suggest that the mean flow plays an important role in the observed wavenumber. Analytical results are developed quantifying the mean flow sources, generated by amplitude gradients, and its effect on the pattern wavenumber for a large-aspect-ratio cylinder with a ramped boundary. Numerical results are found to agree well with these analytical predictions. We gain further insight into the role of mean flow in pattern dynamics by employing a novel method of quenching the mean flow numerically. Simulations of a spiral defect chaos state where the mean flow is suddenly quenched is found to remove the time dependence, increase the wavenumber and make the pattern more angular in nature., Comment: 9 pages, 10 figures

Pattern formation in Rayleigh-Benard convection in a large-aspect-ratio cylinder with a radial ramp in the plate separation is studied analytically and numerically by performing numerical simulations of the Boussinesq equations. A horizontal mean flow and a vertical large scale counterflow are quantified and used to understand the pattern wavenumber. Our results suggest that the mean flow, generated by amplitude gradients, plays an important role in the roll compression observed as the control parameter is increased. Near threshold the mean flow has a quadrupole dependence with a single vortex in each quadrant while away from threshold the mean flow exhibits an octupole dependence with a counter-rotating pair of vortices in each quadrant. This is confirmed analytically using the amplitude equation and Cross-Newell mean flow equation. By performing numerical experiments the large scale counterflow is also found to aid in the roll compression away from threshold but to a much lesser degree. Our results yield an understanding of the pattern wavenumbers observed in experiment away from threshold and suggest that near threshold the mean flow and large scale counterflow are not responsible for the observed shift to smaller than critical wavenumbers., Comment: 10 pages, 13 figures

The origin of the power-law decay measured in the power spectra of low Prandtl number Rayleigh-Benard convection near the onset of chaos is addressed using long time numerical simulations of the three-dimensional Boussinesq equations in cylindrical domains. The power-law is found to arise from quasi-discontinuous changes in the slope of the time series of the heat transport associated with the nucleation of dislocation pairs and roll pinch-off events. For larger frequencies, the power spectra decay exponentially as expected for time continuous deterministic dynamics., Comment: (10 pages, 6 figures)