Four-phase patterns in forced oscillatory systems

We investigate pattern formation in self-oscillating systems forced by an external periodic perturbation. Experimental observations and numerical studies of reaction-diffusion systems and an analysis of an amplitude equation are presented. The oscillations in each of these systems entrain to rational multiples of the perturbation frequency for certain values of the forcing frequency and amplitude. We focus on the subharmonic resonant case where the system locks at one fourth the driving frequency, and four-phase rotating spiral patterns are observed at low forcing amplitudes. The spiral patterns are studied using an amplitude equation for periodically forced oscillating systems. The analysis predicts a bifurcation (with increasing forcing) from rotating four-phase spirals to standing two-phase patterns. This bifurcation is also found in periodically forced reaction-diffusion equations, the FitzHugh-Nagumo and Brusselator models, even far from the onset of oscillations where the amplitude equation analysis is not strictly valid. In a Belousov-Zhabotinsky chemical system periodically forced with light we also observe four-phase rotating spiral wave patterns. However, we have not observed the transition to standing two-phase patterns, possibly because with increasing light intensity the reaction kinetics become excitable rather than oscillatory.Comment: 11 page

Advanced Numerical Methods and Software Approaches for Semiconductor Device Simulation

In this article we concisely present several modern strategies that are applicable to driftdominated carrier transport in higher-order deterministic models such as the driftdiffusion, hydrodynamic, and quantum hydrodynamic systems. The approaches include extensions of “upwind” and artificial dissipation schemes, generalization of the traditional Scharfetter – Gummel approach, Petrov – Galerkin and streamline-upwind Petrov Galerkin (SUPG), “entropy” variables, transformations, least-squares mixed methods and other stabilized Galerkin schemes such as Galerkin least squares and discontinuous Galerkin schemes. The treatment is representative rather than an exhaustive review and several schemes are mentioned only briefly with appropriate reference to the literature. Some of the methods have been applied to the semiconductor device problem while others are still in the early stages of development for this class of applications. We have included numerical examples from our recent research tests with some of the methods. A second aspect of the work deals with algorithms that employ unstructured grids in conjunction with adaptive refinement strategies. The full benefits of such approaches have not yet been developed in this application area and we emphasize the need for further work on analysis, data structures and software to support adaptivity. Finally, we briefly consider some aspects of software frameworks. These include dial-an-operator approaches such as that used in the industrial simulator PROPHET, and object-oriented software support such as those in the SANDIA National Laboratory framework SIERRA