Design of Parametrically Forced Patterns and Quasipatterns

The Faraday wave experiment is a classic example of a system driven by parametric forcing, and it produces a wide range of complex patterns, including superlattice patterns and quasipatterns. Nonlinear three-wave interactions between driven and weakly damped modes play a key role in determining which patterns are favored. We use this idea to design single and multifrequency forcing functions that produce examples of superlattice patterns and quasipatterns in a new model PDE with parametric forcing. We make quantitative comparisons between the predicted patterns and the solutions of the PDE. Unexpectedly, the agreement is good only for parameter values very close to onset. The reason that the range of validity is limited is that the theory requires strong damping of all modes apart from the driven pattern-forming modes. This is in conflict with the requirement for weak damping if three-wave coupling is to influence pattern selection effectively. We distinguish the two different ways that three-wave interactions can be used to stabilize quasipatterns, and we present examples of 12-, 14-, and 20-fold approximate quasipatterns. We identify which computational domains provide the most accurate approximations to 12-fold quasipatterns and systematically investigate the Fourier spectra of the most accurate approximations

Numerical simulation of Faraday waves

We simulate numerically the full dynamics of Faraday waves in three dimensions for two incompressible and immiscible viscous fluids. The Navier-Stokes equations are solved using a finite-difference projection method coupled with a front-tracking method for the interface between the two fluids. The domain of calculation is periodic in the horizontal directions and bounded in the vertical direction by two rigid horizontal plates. The critical accelerations and wavenumbers, as well as the temporal behaviour at onset are compared with the results of the linear Floquet analysis of Kumar and Tuckerman [J. Fluid Mech. 279, 49-68 (1994)]. The finite amplitude results are compared with the experiments of Kityk et al. [Phys. Rev. E 72, 036209 (2005)]. In particular we reproduce the detailed spatiotemporal spectrum of both square and hexagonal patterns within experimental uncertainty

Coupling between plate vibration and acoustic radiation

A detailed numerical investigation of the coupling between the vibration of a flexible plate and the acoustic radiation is performed. The nonlinear Euler equations are used to describe the acoustic fluid while the nonlinear plate equation is used to describe the plate vibration. Linear, nonlinear, and quasi-periodic or chaotic vibrations and the resultant acoustic radiation are analyzed. We find that for the linear plate response, acoustic coupling is negligible. However, for the nonlinear and chaotic responses, acoustic coupling has a significant effect on the vibration level as the loading increases. The radiated pressure from a plate undergoing nonlinear or chaotic vibrations is found to propagate nonlinearly into the far-field. However, the nonlinearity due to wave propagation is much weaker than that due to the plate vibrations. As the acoustic wave propagates into the far-field, the relative difference in level between the fundamental and its harmonics and subharmonics decreases with distance

Optimal Subharmonic Entrainment

For many natural and engineered systems, a central function or design goal is the synchronization of one or more rhythmic or oscillating processes to an external forcing signal, which may be periodic on a different time-scale from the actuated process. Such subharmonic synchrony, which is dynamically established when N control cycles occur for every M cycles of a forced oscillator, is referred to as N:M entrainment. In many applications, entrainment must be established in an optimal manner, for example by minimizing control energy or the transient time to phase locking. We present a theory for deriving inputs that establish subharmonic N:M entrainment of general nonlinear oscillators, or of collections of rhythmic dynamical units, while optimizing such objectives. Ordinary differential equation models of oscillating systems are reduced to phase variable representations, each of which consists of a natural frequency and phase response curve. Formal averaging and the calculus of variations are then applied to such reduced models in order to derive optimal subharmonic entrainment waveforms. The optimal entrainment of a canonical model for a spiking neuron is used to illustrate this approach, which is readily extended to arbitrary oscillating systems

A new frequency domain representation and analysis for subharmonic oscillation

For a weakly nonlinear oscillator, the frequency domain Volterra kernels, often called the generalised frequency response functions, can provide accurate analysis of the response in terms of amplitudes and frequencies, in a transparent algebraic way. However a Volterra series representation based analysis will become void for nonlinear oscillators that exhibit subharmonics, and the problem of finding a solution in this situation has been mainly treated by the traditional analytical approximation methods. In this paper a novel method is developed, by extending the frequency domain Volterra representation to the subharmonic situation, to allow the advantages and the benefits associated with the traditional generalised frequency response functions to be applied to severely nonlinear systems that exhibit subharmonic behaviour

Squeezed-light generation in a nonlinear planar waveguide with a periodic corrugation

Two-mode nonlinear interaction (second-harmonic and second-subharmonic generation) in a planar waveguide with a small periodic corrugation at the surface is studied. Scattering of the interacting fields on the corrugation leads to constructive interference that enhances the nonlinear process provided that all the interactions are phase matched. Conditions for the overall phase matching are found. Compared with a perfectly quasi-phase-matched waveguide, better values of squeezing as well as higher intensities are reached under these conditions. Procedure for finding optimum values of parameters for squeezed-light generation is described.Comment: 14 pages, 14 figure

Resonances and superlattice pattern stabilization in two-frequency forced Faraday waves

We investigate the role weakly damped modes play in the selection of Faraday wave patterns forced with rationally-related frequency components m*omega and n*omega. We use symmetry considerations to argue for the special importance of the weakly damped modes oscillating with twice the frequency of the critical mode, and those oscillating primarily with the "difference frequency" |n-m|*omega and the "sum frequency" (n+m)*omega. We then perform a weakly nonlinear analysis using equations of Zhang and Vinals (1997, J. Fluid Mech. 336) which apply to small-amplitude waves on weakly inviscid, semi-infinite fluid layers. For weak damping and forcing and one-dimensional waves, we perform a perturbation expansion through fourth order which yields analytical expressions for onset parameters and the cubic bifurcation coefficient that determines wave amplitude as a function of forcing near onset. For stronger damping and forcing we numerically compute these same parameters, as well as the cubic cross-coupling coefficient for competing waves travelling at an angle theta relative to each other. The resonance effects predicted by symmetry are borne out in the perturbation results for one spatial dimension, and are supported by the numerical results in two dimensions. The difference frequency resonance plays a key role in stabilizing superlattice patterns of the SL-I type observed by Kudrolli, Pier and Gollub (1998, Physica D 123).Comment: 41 pages, 13 figures; corrected figure 1b and minor typos in tex