A Neural Pattern Generator that Exhibits Arousal-Dependant Human Gait Transitions

A neural pattern generator based upon a non-linear cooperative-competitive feedback neural network is presented. It can generate the two standard human gaits: the walk and the run. A scalar arousal or GO signal causes a bifurcation from one gait to the next. Although these two gaits are qualitatively different, they both have the same limb order and may exhibit oscillation frequencies that overlap. The model simulates the walk and the run via qualitatively different waveform shapes. The fraction of cycle that activity is above threshold distinguishes the two gaits, much as the duty cycles of the feet are longer in the walk than in the run.

Statistical Properties of Single and Competing Non-Linear Fast-Slow Oscillators in Noise

Statistical properties offast-slow Ellias-Grossberg oscillators are studied in response to deterministic and noisy inputs. Oscillatory responses remain stable in noise due to the slow inhibitory variable, which establishes an adaptation level that centers the oscillatory responses of the fast excitatory variable to deterministic and noisy inputs. Competitive interactions between oscillators improve the stability in noise. Although individual oscillation amplitudes decrease with input amplitude, the average to'tal activity increases with input amplitude, thereby suggesting that oscillator output is evaluated by a slow process at downstream network sites.

A Neural Pattern Generator that Exhibits Frequency-Dependent In-Phase and Anti-Phase Oscillations

This article describes a. neural pattern generator based on a cooperative-competitive feedback neural network. The two-channel version of the generator supports both in-phase and anti-phase oscillations. A scalar arousal level controls both the oscillation phase and frequency. As arousal increases, oscillation frequency increases and bifurcations from in-phase to anti-phase, or anti-phase to in-phase oscillations can occur. Coupled versions of the model exhibit oscillatory patterns which correspond to the gaits used in locomotion and other oscillatory movements by various animals.

Monitoring the vegetation start of season (SOS) across the island of Ireland using the MERIS global vegetation index

The aim of this study was to develop a methodology, based on satellite remote sensing, to estimate the vegetation Start of Season (SOS) across the whole island of Ireland on an annual basis. This growing body of research is known as Land Surface Phenology (LSP) monitoring. The SOS was estimated for each year from a 7-year time series of 10-day composited, 1.2 km reduced resolution MERIS Global Vegetation Index (MGVI) data from 2003 to 2009, using the time series analysis software, TIMESAT. The selection of a 10-day composite period was guided by in-situ observations of leaf unfolding and cloud cover at representative point locations on the island. The MGVI time series was smoothed and the SOS metric extracted at a point corresponding to 20% of the seasonal MGVI amplitude. The SOS metric was extracted on a per pixel basis and gridded for national scale coverage. There were consistent spatial patterns in the SOS grids which were replicated on an annual basis and were qualitatively linked to variation in landcover. Analysis revealed that three statistically separable groups of CORINE Land Cover (CLC) classes could be derived from differences in the SOS, namely agricultural and forest land cover types, peat bogs, and natural and semi-natural vegetation types. These groups demonstrated that managed vegetation, e.g. pastures has a significantly earlier SOS than in unmanaged vegetation e.g. natural grasslands. There was also interannual spatio-temporal variability in the SOS. Such variability was highlighted in a series of anomaly grids showing variation from the 7-year mean SOS. An initial climate analysis indicated that an anomalously cold winter and spring in 2005/2006, linked to a negative North Atlantic Oscillation index value, delayed the 2006 SOS countrywide, while in other years the SOS anomalies showed more complex variation. A correlation study using air temperature as a climate variable revealed the spatial complexity of the air temperature-SOS relationship across the Republic of Ireland as the timing of maximum correlation varied from November to April depending on location. The SOS was found to occur earlier due to warmer winters in the Southeast while it was later with warmer winters in the Northwest. The inverse pattern emerged in the spatial patterns of the spring correlates. This contrasting pattern would appear to be linked to vegetation management as arable cropping is typically practiced in the southeast while there is mixed agriculture and mostly pastures to the west. Therefore, land use as well as air temperature appears to be an important determinant of national scale patterns in the SOS. The TIMESAT tool formed a crucial component of the estimation of SOS across the country in all seven years as it minimised the negative impact of noise and data dropouts in the MGVI time series by applying a smoothing algorithm. The extracted SOS metric was sensitive to temporal and spatial variation in land surface vegetation seasonality while the spatial patterns in the gridded SOS estimates aligned with those in landcover type. The methodology can be extended for a longer time series of FAPAR as MERIS will be replaced by the ESA Sentinel mission in 2013, while the availability of full resolution (300m) MERIS FAPAR and equivalent sensor products holds the possibility of monitoring finer scale seasonality variation. This study has shown the utility of the SOS metric as an indicator of spatiotemporal variability in vegetation phenology, as well as a correlate of other environmental variables such as air temperature. However, the satellite-based method is not seen as a replacement of ground-based observations, but rather as a complementary approach to studying vegetation phenology at the national scale. In future, the method can be extended to extract other metrics of the seasonal cycle in order to gain a more comprehensive view of seasonal vegetation development.

Two closely coupled lasers, dynamics & applications

The dynamics of two mutually coupled identical single-mode semi-conductor lasers are theoretically investigated. For small separation and large coupling between the lasers, symmetry-broken one-colour states are shown to be stable. In this case the light output of the lasers have significantly different intensities whilst at the same time the lasers are locked to a single common frequency. For intermediate coupling we observe stable two-colour states, where both single-mode lasers lase simultaneously at two optical frequencies which are separated by up to 150 GHz. For low coupling but possibly large separation, the frequency of the relaxation oscillations of the freerunning lasers defines the dynamics. Chaotic and quasi-periodic states are identified and shown to be stable. For weak coupling undamped relaxation oscillations dominate where each laser is locked to three or more odd number of colours spaced by the relaxation oscillation frequency. It is shown that the instabilities that lead to these states are directly connected to the two colour mechanism where the change in the number of optical colours due to a change in the plane of oscillation. At initial coupling, in-phase and anti-phase one colour states are shown to emerge from “on” uncoupled lasers using a perturbation method. Similarly symmetry-broken one-colour states come from considering one free-running laser initially “on” and the other laser initially “off”. The mechanism that leads to a bi-stability between in-phase and anti-phase one-colour states is understood. Due to an equivariant phase space symmetry of being able to exchange the identical lasers, a symmetric and symmetry-broken variant of all states mentioned above exists and is shown to be stable. Using a five dimensional model we identify the bifurcation structure which is responsible for the appearance of symmetric and symmetry-broken one-colour, symmetric and symmetry-broken two-colour, symmetric and symmetry-broken undamped relaxation oscillations, symmetric and symmetry-broken quasi-periodic, and symmetric and symmetry-broken chaotic states. As symmetry-broken states always exist in pairs, they naturally give rise to bi-stability. Several of these states show multistabilities between symmetric and symmetry-broken variants and among states. Three memory elements on the basis of bi-stabilities in one and two colour states for two coupled single-mode lasers are proposed. The switching performance of selected designs of optical memory elements is studied numerically.

Vulnerability and adaptation to climate change in Guam – a case-study of the perceptions, attitudes, behaviors, and knowledge of a small rural community towards their watershed

This thesis argues that examining the attitudes, perceptions, behaviors, and knowledge of a community towards their specific watershed can reveal their social vulnerability to climate change. Understanding and incorporating these elements of the human dimension in coastal zone management will lead to efficient and effective strategies that safeguard the natural resources for the benefit of the community. By having healthy natural resources, ecological and community resilience to climate change will increase, thus decreasing vulnerability. In the Pacific Ocean, climate and SLR are strongly modulated by the El Niño Southern Oscillation. SLR is three times the global average in the Western Pacific Ocean (Merrifield and Maltrud 2011; Merrifield 2011). Changes in annual rainfall in the Western North Pacific sub‐region from 1950-2010 show that islands in the east are getting much less than in the past, while the islands in the west are getting slightly more rainfall (Keener et al. 2013). For Guam, a small island owned by the United States and located in the Western Pacific Ocean, these factors mean that SLR is higher than any other place in the world and will most likely see increased precipitation. Knowing this, the social vulnerability may be examined. Thus, a case-study of the community residing in the Manell and Geus watersheds was conducted on the island of Guam. Measuring their perceptions, attitudes, knowledge, and behaviors should bring to light their vulnerability to climate change. In order to accomplish this, a household survey was administered from July through August 2010. Approximately 350 surveys were analysed using SPSS. To supplement this quantitative data, informal interviews were conducted with the elders of the community to glean traditional ecological knowledge about perceived climate change. A GIS analysis was conducted to understand the physical geography of the Manell and Geus watersheds. This information about the human dimension is valuable to CZM managers. It may be incorporated into strategic watershed plans, to better administer the natural resources within the coastal zone. The research conducted in this thesis is the basis of a recent watershed management plan for the Guam Coastal Management Program (see King 2014).

Nonlinear piezoelectricity in electroelastic energy harvesters: Modeling and experimental identification

We propose and experimentally validate a first-principles based model for the nonlinear piezoelectric response of an electroelastic energy harvester. The analysis herein highlights the importance of modeling inherent piezoelectric nonlinearities that are not limited to higher order elastic effects but also include nonlinear coupling to a power harvesting circuit. Furthermore, a nonlinear damping mechanism is shown to accurately restrict the amplitude and bandwidth of the frequency response. The linear piezoelectric modeling framework widely accepted for theoretical investigations is demonstrated to be a weak presumption for near-resonant excitation amplitudes as low as 0.5 g in a prefabricated bimorph whose oscillation amplitudes remain geometrically linear for the full range of experimental tests performed (never exceeding 0.25% of the cantilever overhang length). Nonlinear coefficients are identified via a nonlinear least-squares optimization algorithm that utilizes an approximate analytic solution obtained by the method of harmonic balance. For lead zirconate titanate (PZT-5H), we obtained a fourth order elastic tensor component of c1111p =-3.6673× 1017 N/m2 and a fourth order electroelastic tensor value of e3111 =1.7212× 108 m/V. © 2010 American Institute of Physics.

Surface plasmon resonance in nanostructured metal films under the Kretschmann configuration

We systematically investigated the surface plasmon resonance in one-dimensional (1D) subwavelength nanostructured metal films under the Kretschmann configuration. We calculated the reflectance, transmittance, and absorption for varying the dielectric fill factor, the period of the 1D nanostructure, and the metal film thickness. We have found that the small dielectric slits in the metal films reduce the surface plasmon resonance angle and move it toward the critical angle for total internal reflection. The reduction in surface plasmon resonance angle in nanostructured metal films is due to the increased intrinsic free electron oscillation frequency in metal nanostructures. Also we have found that the increasing the spatial frequency of the 1D nanograting reduces the surface plasmon resonance angle, which indicates that less momentum is needed to match the momentum of the surface plasmon-polariton. The variation in the nanostructured metal film thickness changes the resonance angle slightly, but mainly remains as a mean to adjust the coupling between the incident optical wave and the surface plasmon-polariton wave. © 2009 American Institute of Physics.

Pulsating tandem microbubble for localized and directional single-cell membrane poration.

The interaction of laser-generated tandem microbubble (maximum diameter of about 50  μm) with single (rat mammary carcinoma) cells is investigated in a 25-μm liquid layer. Antiphase and coupled oscillation of the tandem microbubble leads to the formation of alternating, directional microjets (with max microstreaming velocity of 10  m/s) and vortices (max vorticity of 350 000  s{-1}) in opposite directions. Localized and directional membrane poration (200 nm to 2  μm in pore size) can be produced by the tandem microbubble in an orientation and proximity-dependent manner, which is absent from a single oscillating microbubble of comparable size and at the same stand-off distance.

Indication of electron neutrino appearance from an accelerator-produced off-axis muon neutrino beam.

The T2K experiment observes indications of ν(μ) → ν(e) appearance in data accumulated with 1.43×10(20) protons on target. Six events pass all selection criteria at the far detector. In a three-flavor neutrino oscillation scenario with |Δm(23)(2)| = 2.4×10(-3)  eV(2), sin(2)2θ(23) = 1 and sin(2)2θ(13) = 0, the expected number of such events is 1.5±0.3(syst). Under this hypothesis, the probability to observe six or more candidate events is 7×10(-3), equivalent to 2.5σ significance. At 90% C.L., the data are consistent with 0.03(0.04) < sin(2)2θ(13) < 0.28(0.34) for δ(CP) = 0 and a normal (inverted) hierarchy.

A natural extension of the conventional finite volume method into polygonal unstructured meshes for CFD application.

A new general cell-centered solution procedure based upon the conventional control or finite volume (CV or FV) approach has been developed for numerical heat transfer and fluid flow which encompasses both structured and unstructured meshes for any kind of mixed polygon cell. Unlike conventional FV methods for structured and block structured meshes and both FV and FE methods for unstructured meshes, the irregular control volume (ICV) method does not require the shape of the element or cell to be predefined because it simply exploits the concept of fluxes across cell faces. That is, the ICV method enables meshes employing mixtures of triangular, quadrilateral, and any other higher order polygonal cells to be exploited using a single solution procedure. The ICV approach otherwise preserves all the desirable features of conventional FV procedures for a structured mesh; in the current implementation, collocation of variables at cell centers is used with a Rhie and Chow interpolation (to suppress pressure oscillation in the flow field) in the context of the SIMPLE pressure correction solution procedure. In fact all other FV structured mesh-based methods may be perceived as a subset of the ICV formulation. The new ICV formulation is benchmarked using two standard computational fluid dynamics (CFD) problems i.e., the moving lid cavity and the natural convection driven cavity. Both cases were solved with a variety of structured and unstructured meshes, the latter exploiting mixed polygonal cell meshes. The polygonal mesh experiments show a higher degree of accuracy for equivalent meshes (in nodal density terms) using triangular or quadrilateral cells; these results may be interpreted in a manner similar to the CUPID scheme used in structured meshes for reducing numerical diffusion for flows with changing direction.

A two-dimensional numerical investigation of the oscillatory flow behaviour in rectangular fire compartments with a single horizontal ceiling vent

This paper presents a comparison of fire field model predictions with experiment for the case of a fire within a compartment which is vented (buoyancydriven) to the outside by a single horizontal ceiling vent. Unlike previous work, the mathematical model does not employ a mixing ratio to represent vent temperatures but allows the model to predict vent temperatures a priori. The experiment suggests that the flow through the vent produces oscillatory behaviour in vent temperatures with puffs of smoke emerging from the fire compartment. This type of flow is also predicted by the fire field model. While the numerical predictions are in good qualitative agreement with observations, they overpredict the amplitudes of the temperature oscillations within the vent and also the compartment temperatures. The discrepancies are thought to be due to three-dimensional effects not accounted for in this model as well as using standard ‘practices’ normally used by the community with regards to discretization and turbulence models. Furthermore, it is important to note that the use of the k–ε turbulence model in a transient mode, as is used here, may have a significant effect on the results. The numerical results also suggest that a linear relationship exists between the frequency of vent temperature oscillation (n) and the heat release rate (Q0) of the type n∝Q0.290, similar to that observed for compartments with two horizontal vents. This relationship is predicted to occur only for heat release rates below a critical value. Furthermore, the vent discharge coefficient is found to vary in an oscillatory fashion with a mean value of 0.58. Below the critical heat release rate the mean discharge coefficient is found to be insensitive to fire size.

Modeling the dynamics of electromagnetically agitated metallurgical flows: levitation, cold crucible, aluminium electrolysis, etc.

We present practical modelling techniques for electromagnetically agitated liquid metal flows involving dynamic change of the fluid volume and shape during melting and the free surface oscillation. Typically the electromagnetic field is strongly coupled to the free surface dynamics and the heat-mass transfer. Accurate pseudo-spectral code and the k-omega turbulence model modified for complex and transitional flows with free surfaces are used for these simulations. The considered examples include magnetic suspension melting, induction scull remelting (cold crucible), levitation and aluminium electrolysis cells. The process control and the energy savings issues are analysed.

Dynamic solid mechanics using finite volume methods

A procedure for evaluating the dynamic structural response of elastic solid domains is presented. A prerequisite for the analysis of dynamic fluid–structure interaction is the use of a consistent set of finite volume (FV) methods on a single unstructured mesh. This paper describes a three-dimensional (3D) FV, vertex-based method for dynamic solid mechanics. A novel Newmark predictor–corrector implicit scheme was developed to provide time accurate solutions and the scheme was evaluated on a 3D cantilever problem. By employing a small amount of viscous damping, very accurate predictions of the fundamental natural frequency were obtained with respect to both the amplitude and period of oscillation. This scheme has been implemented into the multi-physics modelling software framework, PHYSICA, for later application to full dynamic fluid structure interaction.

Modelling electromagnetically levitated liquid droplet oscillations

This work comprises accurate computational analysis of levitated liquid droplet oscillations in AC and DC magnetic fields. The AC magnetic field interacting with the induced electric current within the liquid metal droplet generates intense fluid flow and the coupled free surface oscillations. The pseudo-spectral technique is used to solve the turbulent fluid flow equations for the continuously dynamically transformed axisymmetric fluid volume. The volume electromagnetic force distribution is updated with the shape and position change. We start with the ideal fluid test case for undamped Rayleigh frequency oscillations in the absence of gravity, and then add the viscous and the DC magnetic field damping. The oscillation frequency spectra are further analysed for droplets levitated against gravity in AC and DC magnetic fields at various combinations. In the extreme case electrically poorly conducting, diamagnetic droplet (water) levitation dynamics are simulated. Applications are aimed at pure electromagnetic material processing techniques and the material properties measurements in uncontaminated conditions.