A frequency-independent boundary element method for scattering by two-dimensional screens and apertures

We propose and analyse a hybrid numerical–asymptotic hp boundary element method (BEM) for time-harmonic scattering of an incident plane wave by an arbitrary collinear array of sound-soft two-dimensional screens. Our method uses an approximation space enriched with oscillatory basis functions, chosen to capture the high-frequency asymptotics of the solution. We provide a rigorous frequency-explicit error analysis which proves that the method converges exponentially as the number of degrees of freedom N increases, and that to achieve any desired accuracy it is sufficient to increase N in proportion to the square of the logarithm of the frequency as the frequency increases (standard BEMs require N to increase at least linearly with frequency to retain accuracy). Our numerical results suggest that fixed accuracy can in fact be achieved at arbitrarily high frequencies with a frequency-independent computational cost, when the oscillatory integrals required for implementation are computed using Filon quadrature. We also show how our method can be applied to the complementary ‘breakwater’ problem of propagation through an aperture in an infinite sound-hard screen.

Ionospheric Signatures of Pulsed Magnetopause Reconnection

The concept of zero-flow equilibria of the magnetosphere-ionosphere system leads to a large number of predictions concerning the ionospheric signatures of pulsed magnetopause reconnection. These include: poleward-moving F-region electron temperature enhancements and associated transient 630nm emission; associated poleward plasma flow which, compared to the pulsed variation of the reconnection rate, is highly smoothed by induction effects; oscillatory latitudinal motion of the open/closed field line boundary; phase lag of plasma flow enhancements after equatorward motions of the boundary; azimuthal plasma flow bursts, coincident in time and space with the 630nm-dominant auroral transients, only when the magnitude of the By component of the interplanetary magnetic field (IMF) is large; azimuthal-then-poleward motion of 630nm-dominant transients at a velocity which at all times equals the internal plasma flow velocity; 557.7nm-dominant transients on one edge of the 630nm-dominant transient (initially, and for large |By|, on the poleward or equatorward edge depending on the polarity of IMF By); tailward expansion of the flow response at several km s-1; and discrete steps in the cusp ion dispersion signature between the polewardmoving structures. This paper discusses these predictions and how all have recently been confirmed by combinations of observations by optical instruments on the Svalbard Islands, the EISCAT radars and the DMSP and DE satellites.

Movement intention detection from autocorrelation of EEG for BCI

Movement intention detection is important for development of intuitive movement based Brain Computer Interfaces (BCI). Various complex oscillatory processes are involved in producing voluntary movement intention. In this paper, temporal dynamics of electroencephalography (EEG) associated with movement intention and execution were studied using autocorrelation. It was observed that the trend of decay of autocorrelation of EEG changes before and during the voluntary movement. A novel feature for movement intention detection was developed based on relaxation time of autocorrelation obtained by fitting exponential decay curve to the autocorrelation. This new single trial feature was used to classify voluntary finger tapping trials from resting state trials with peak accuracy of 76.7%. The performance of autocorrelation analysis was compared with Motor-Related Cortical Potentials (MRCP).

Subjective evaluation of experimental dyspnoea: effects of isocapnia and repeated exposure

Resistive respiratory loading is an established stimulus for the induction of experimental dyspnoea. In comparison to unloaded breathing, resistive loaded breathing alters end-tidal CO2 (PETCO2), which has independent physiological effects (e.g. upon cerebral blood flow). We investigated the subjective effects of resistive loaded breathing with stabilized PETCO2 (isocapnia) during manual control of inspired gases on varying baseline levels of mild hypercapnia increased PETCO2). Furthermore, to investigate whether perceptual habituation to dyspnoea stimuli occurs, the study was repeated over four experimental sessions. Isocapnic hypercapnia did not affect dyspnoea unpleasantness during resistive loading. A post hoc analysis revealed a small increase of respiratory unpleasantness during unloaded breathing at +0.6 kPa, the level that reliably induced isocapnia. We didnot observe perceptual habituation over the four sessions. We conclude that isocapnic respiratory loading allows stable induction of respiratory unpleasantness, making it a good stimulus for multi-session studies of dyspnoea.

EMD performance comparison: single vs double floating points

Empirical mode decomposition (EMD) is a data-driven method used to decompose data into oscillatory components. This paper examines to what extent the defined algorithm for EMD might be susceptible to data format. Two key issues with EMD are its stability and computational speed. This paper shows that for a given signal there is no significant difference between results obtained with single (binary32) and double (binary64) floating points precision. This implies that there is no benefit in increasing floating point precision when performing EMD on devices optimised for single floating point format, such as graphical processing units (GPUs).

Objective Empirical Mode Decomposition metric

Empirical Mode Decomposition (EMD) is a data driven technique for extraction of oscillatory components from data. Although it has been introduced over 15 years ago, its mathematical foundations are still missing which also implies lack of objective metrics for decomposed set evaluation. Most common technique for assessing results of EMD is their visual inspection, which is very subjective. This article provides objective measures for assessing EMD results based on the original definition of oscillatory components.

Why cancer patients choose in-patient complementary therapy in palliative care: a qualitative study at Arokhayasala Hospice in Thailand

Cancer patients often choose complementary and alternative medicine (CAM) in palliative care, often in addition to conventional treatment and without medical advice or approval. Herbal medicines (HM) are the most commonly used type of CAM, but rarely available on an in-patient basis for palliative care. The motivations which lead very ill patients to travel far to receive such therapies are not clear. A qualitative study was therefore carried out to investigate influences on choosing to attend a CAM herbal hospice, to identify cancer patients’ main concerns about end-of-life care. Semi-structured interviews with 32 patients were conducted and analysed using thematic analysis. Patients were recruited from Arokhayasala, a Buddhist cancer hospice in Thailand which provides CAM, in the form of HM, a restricted diet, Thai yoga, deep-breathing exercises, meditation, chanting, Dhamma, laughter and music therapy, free-of-charge. The main factors influencing decision-making were a positive attitude towards HMs and previous use of them, dissatisfaction with conventional treatment, the home environment and their relationships with hospital doctors. Patients’ own perceptions and experiences were more important in making the decision to use CAM, and especially HM, in palliative cancer care than referral by healthcare professionals or scientific evidence of efficacy. Patients were prepared to travel far and live away from home to receive such care, especially as it was cost-free. In view of patients’ previously stated satisfaction with the regime at the Arokhayasala, these findings may be relevant to the provision of in-patient cancer palliative care to other patients.

Film forming solutions based on gelatin and poly(vinyl alcohol) blends: Thermal and rheological characterizations

The objective of this work was to study the theological and thermal properties of film forming solutions (FFS) based on blends of gelatin and poly(vinyl alcohol) (PVA). The effect of the PVA concentration and plasticizer presence on the flow behavior, and viscoelastic and thermal properties of FFS was studied by steady-shear flow and oscillatory experiments, and also, by microcalorimetry. The FB presented Newtonian behavior at 30 degrees C, and the viscosity was not affected neither by the PVA concentration nor by the plasticizer. All FFS presented a phase transition during tests applying temperature scanning. It was verified that the PVA affected the viscoelastic properties of FFS by dilution of gelatin. This behavior was confirmed by microcalorimetric analysis. The behaviors of the storage (G`) and loss (G ``) moduli as a function of frequency of FFS obtained at 5 degrees C were typical of physical gels; with the G` higher than the G ``. The strength of the gels was affected by the PVA concentration. (C) 2009 Elsevier Ltd. All rights reserved.

SIMULATIONS OF WINDS OF WEAK-LINED T TAURI STARS. II. THE EFFECTS OF A TILTED MAGNETOSPHERE AND PLANETARY INTERACTIONS

Based on our previous work, we investigate here the effects on the wind and magnetospheric structures of weak-lined T Tauri stars due to a misalignment between the axis of rotation of the star and its magnetic dipole moment vector. In such a configuration, the system loses the axisymmetry presented in the aligned case, requiring a fully three-dimensional (3D) approach. We perform 3D numerical magnetohydrodynamic simulations of stellar winds and study the effects caused by different model parameters, namely the misalignment angle theta(t), the stellar period of rotation, the plasma-beta, and the heating index.. Our simulations take into account the interplay between the wind and the stellar magnetic field during the time evolution. The system reaches a periodic behavior with the same rotational period of the star. We show that the magnetic field lines present an oscillatory pattern. Furthermore, we obtain that by increasing theta(t), the wind velocity increases, especially in the case of strong magnetic field and relatively rapid stellar rotation. Our 3D, time-dependent wind models allow us to study the interaction of a magnetized wind with a magnetized extrasolar planet. Such interaction gives rise to reconnection, generating electrons that propagate along the planet`s magnetic field lines and produce electron cyclotron radiation at radio wavelengths. The power released in the interaction depends on the planet`s magnetic field intensity, its orbital radius, and on the stellar wind local characteristics. We find that a close-in Jupiter-like planet orbiting at 0.05 AU presents a radio power that is similar to 5 orders of magnitude larger than the one observed in Jupiter, which suggests that the stellar wind from a young star has the potential to generate strong planetary radio emission that could be detected in the near future with LOFAR. This radio power varies according to the phase of rotation of the star. For three selected simulations, we find a variation of the radio power of a factor 1.3-3.7, depending on theta(t). Moreover, we extend the investigation done in Vidotto et al. and analyze whether winds from misaligned stellar magnetospheres could cause a significant effect on planetary migration. Compared to the aligned case, we show that the timescale tau(w) for an appreciable radial motion of the planet is shorter for larger misalignment angles. While for the aligned case tau(w) similar or equal to 100 Myr, for a stellar magnetosphere tilted by theta(t) = 30 degrees, tau(w) ranges from similar to 40 to 70 Myr for a planet located at a radius of 0.05 AU. Further reduction on tau(w) might occur for even larger misalignment angles and/or different wind parameters.

Some Possible Dynamical Constraints for Life`s Origin

Oscillating biochemical reactions are common in cell dynamics and could be closely related to the emergence of the life phenomenon itself. In this work, we study the dynamical features of some classical chemical or biochemical oscillators where the effect of cell volume changes is explicitly considered. Such analysis enables us to find some general conditions about the cell membrane to preserve such oscillatory patterns, of possible relevance to hypothetical primitive cells in which these structures first appeared.

Melanopsin and Clock Genes: Regulation by Light and Endothelin in the Zebrafish ZEM-2S Cell Line

It is well known that clocks are present in brain regions other than the suprachiasmatic nucleus and in many peripheral tissues. In the teleost, Danio rerio, peripheral oscillators can be directly synchronized by light. Danio rerio ZEM-2S embryonic cells respond to light with differential growth: cells kept in constant light exhibited a strong inhibition of proliferation, whereas in cells kept in light:dark (LD) cycles (14L:10D and 10L:14D) or in constant darkness (DD), the doubling times were not statistically different. We demonstrated by RT-PCR followed by PCR that ZEM-2S cells express two melanopsins, Opn4x and Opn4m, and the six Cry genes. The presence of the protein OPN4x was demonstrated by immunocytochemistry. The pattern of temporal expression of the genes Opn4x, Per1, Cry1b, and Clock was studied in ZEM-2S cells kept for five days in 12L:12D or DD. In 12L:12D, the clock genes Per 1 and Cry1b exhibited robust circadian expression, while Opn4x and Clock expression seemed to vary in an ultradian pattern. Both Per1 and Cry1b genes had higher expression during the L phase; Clock gene had an increase in expression coincident with the D phase, and during the subjective night. In DD, the temporal variation of Per1 and Cry1b genes was greatly attenuated but not extinguished, and the higher expressions were shifted to the transition times between subjective day and night, demonstrating that Per and Cry1b were synchronized by the LD cycle. Clock and Opn4x kept the ultradian oscillation, but the rhythm was not statistically significant. As endothelins (ET) have been reported to be a potent stimulator of Per genes in rodents, we investigated the effect of endothelin on ZEM-2S cells, which express ETA receptors. Cells were kept in 12D:12L for five days, and then treated with 10-11 to 10-8M ET-1 for 24h. ET-1 exhibited a biphasic effect on Opn4x expression. At 10-11M, the hormone exerted a highly significant stimulation of Opn4x expression during the L phase and introduced a circadian oscillatory pattern. At 10-10M, a significant increase was seen at ZT21 and ZT0 (i.e., at the end of the D phase and beginning of the L phase), whereas 10-9 and 10-8M ET-1 inhibited the expression of Opn4x at most ZTs. Clock expression was unaffected by 10-8M ET-1; however, in the presence of lower concentrations, the expression was enhanced at some ZTs, strengthening the ultradian oscillation. ET-1 at 10-11 and 10-10M had no effect on Per1 circadian expression; however, 10-9 and 10-8M ET-1 reduced the amplitude of Per1 expression in the beginning of the L phase. ET-1 effects were less evident on Cry 1b. For both genes, the reduction in expression was not sufficient to abolish the circadian oscillatory pattern. Based on these results and data in the literature, a link between ET-1 stimulation of ETA receptors may be established by E4BP4 binding to the promoters and consequent inhibition of gene expression.

Ventrolateral medulla mechanisms involved in cardiorespiratory responses to central chemoreceptor activation in rats

A rise in arterial PCO(2) stimulates breathing and sympathetic activity to the heart and blood vessels. In the present study, we investigated the involvement of the retrotrapezoid nucleus (RTN) and glutamatergic mechanisms in the Botzinger/C1 region (Botz/C1) in these responses. Splanchnic sympathetic nerve discharge (sSND) and phrenic nerve discharge (PND) were recorded in urethane-anesthetized, sino-aortic-denervated, vagotomized, and artificially ventilated rats subjected to hypercapnia (end-expiratory CO(2) from 5% to 10%). Phrenic activity was absent at end-expiratory CO(2) of 4%, and strongly increased when end-expiratory CO(2) reached 10%. Hypercapnia also increased sSND by 103 +/- 7%. Bilateral injections of the GABA-A agonist muscimol (2 mM) into the RTN eliminated the PND and blunted the sSND activation (Delta = +56 +8%) elicited by hypercapnia. Injections of NMDA receptor antagonist AP-5 (100 mM), non-NMDA receptor antagonist 6,7-dinitro-quinoxaline-2,3-dione (DNQX; 100 mM) or metabotropic glutamate receptor antagonist (+/-)-alpha-methyl-4-carboxyphenylglycine (MCPG; 100 mM) bilaterally into the Botz/C1 reduced PND (Delta = +43 +/- 7%, +52 +/- 6% or +56 +/- 11%, respectively). MCPG also reduced sSND (Delta = +41 +/- 7%), whereas AP-5 and DNQX had no effect. In conclusion, the increase in sSND caused by hypercapnia depends on increased activity of the RTN and on metabotropic receptors in the Botz/C1, whereas PND depends on increased RTN activity and both ionotropic and metabotropic receptors in the Botz/C1.

Exact Search Directions for Optimization of Linear and Nonlinear Models Based on Generalized Orthonormal Functions

A novel technique for selecting the poles of orthonormal basis functions (OBF) in Volterra models of any order is presented. It is well-known that the usual large number of parameters required to describe the Volterra kernels can be significantly reduced by representing each kernel using an appropriate basis of orthonormal functions. Such a representation results in the so-called OBF Volterra model, which has a Wiener structure consisting of a linear dynamic generated by the orthonormal basis followed by a nonlinear static mapping given by the Volterra polynomial series. Aiming at optimizing the poles that fully parameterize the orthonormal bases, the exact gradients of the outputs of the orthonormal filters with respect to their poles are computed analytically by using a back-propagation-through-time technique. The expressions relative to the Kautz basis and to generalized orthonormal bases of functions (GOBF) are addressed; the ones related to the Laguerre basis follow straightforwardly as a particular case. The main innovation here is that the dynamic nature of the OBF filters is fully considered in the gradient computations. These gradients provide exact search directions for optimizing the poles of a given orthonormal basis. Such search directions can, in turn, be used as part of an optimization procedure to locate the minimum of a cost-function that takes into account the error of estimation of the system output. The Levenberg-Marquardt algorithm is adopted here as the optimization procedure. Unlike previous related work, the proposed approach relies solely on input-output data measured from the system to be modeled, i.e., no information about the Volterra kernels is required. Examples are presented to illustrate the application of this approach to the modeling of dynamic systems, including a real magnetic levitation system with nonlinear oscillatory behavior.

Oscillation for a second-order neutral differential equation with impulses

We consider a certain type of second-order neutral delay differential systems and we establish two results concerning the oscillation of solutions after the system undergoes controlled abrupt perturbations (called impulses). As a matter of fact, some particular non-impulsive cases of the system are oscillatory already. Thus, we are interested in finding adequate impulse controls under which our system remains oscillatory. (C) 2009 Elsevier Inc. All rights reserved.

A Mass-Conserving Algorithm for Dynamical Lubrication Problems With Cavitation

A numerical algorithm for fully dynamical lubrication problems based on the Elrod-Adams formulation of the Reynolds equation with mass-conserving boundary conditions is described. A simple but effective relaxation scheme is used to update the solution maintaining the complementarity conditions on the variables that represent the pressure and fluid fraction. The equations of motion are discretized in time using Newmark`s scheme, and the dynamical variables are updated within the same relaxation process just mentioned. The good behavior of the proposed algorithm is illustrated in two examples: an oscillatory squeeze flow (for which the exact solution is available) and a dynamically loaded journal bearing. This article is accompanied by the ready-to-compile source code with the implementation of the proposed algorithm. [DOI: 10.1115/1.3142903]