Shear banding in molecular dynamics of polymer melts

In order to establish constitutive equations for a viscoelastic fluid uniform shear flow is usually required. However, in the last 10 years S. Q. Wang and co-workers have demonstrated that some entangled polymers do not flow with the uniform shear rate as usually assumed, but instead choose to separate into fast and slow flowing regions. This phenomenon, known as shear banding, causes flow instabilities and in principle invalidates all rheological measurements when it occurs. In this Letter we report the first observation of shear banding in molecular dynamics simulations of entangled polymer melts. We show that our observations are in a very good agreement with the phenomenology developed by Fielding and Olmsted. Our findings provide a simple way of validating the empirical macroscopic phenomenology of shear banding. © 2012 American Physical Society

Specific detection and localization of microsporidian parasites in invertebrate hosts using in situ hybridization

We designed FISH-probes for two distinct microsporidian clades and demonstrated their application in detecting respectively Nosema/Vairimorpha and Dictyoceola species. We applied them to study the vertical transmission of two microsporidia infecting the amphipod Gammarus duebeni

Localization of calcitonin receptor-like receptor (CLR) and receptor activity-modifying protein 1 (RAMP1) in human gastrointestinal tract

Calcitonin gene-related peptide (CGRP) exerts its diverse effects on vasodilation, nociception, secretion, and motor function through a heterodimeric receptor comprising of calcitonin receptor-like receptor (CLR) and receptor activity-modifying protein 1 (RAMP1). Despite the importance of CLR.RAMP1 in human disease, little is known about its distribution in the human gastrointestinal (GI) tract, where it participates in inflammation and pain. In this study, we determined that CLR and RAMP1 mRNAs are expressed in normal human stomach, ileum and colon by RT-PCR. We next characterized antibodies that we generated to rat CLR and RAMP1 in transfected HEK cells. Having characterized these antibodies in vitro, we then localized CLR-, RAMP1-, CGRP- and intermedin-immunoreactivity (IMD-IR) in various human GI segments. In the stomach, nerve bundles in the myenteric plexus and nerve fibers throughout the circular and longitudinal muscle had prominent CLR-IR. In the proximal colon and ileum, CLR was found in nerve varicosities of the myenteric plexus and surrounding submucosal neurons. Interestingly, CGRP expressing fibers did not co-localize, but were in close proximity to CLR. However, CLR and RAMP1, the two subunits of a functional CGRP receptor were clearly localized in myenteric plexus, where they may form functional cell-surface receptors. IMD, another member of calcitonin peptide family was also found in close proximity to CLR, and like CGRP, did not co-localize with either CLR or RAMP1 receptors. Thus, CGRP and IMD appear to be released locally, where they can mediate their effect on their receptors regulating diverse functions such as inflammation, pain and motility.

Localization of calcitonin receptor-like receptor and receptor activity modifying protein 1 in enteric neurons, dorsal root ganglia, and the spinal cord of the rat

Calcitonin receptor-like receptor (CLR) and receptor activity modifying protein 1 (RAMP1) comprise a receptor for calcitonin gene related peptide (CGRP) and intermedin. Although CGRP is widely expressed in the nervous system, less is known about the localization of CLR and RAMP1. To localize these proteins, we raised antibodies to CLR and RAMP1. Antibodies specifically interacted with CLR and RAMP1 in HEK cells coexpressing rat CLR and RAMP1, determined by Western blotting and immunofluorescence. Fluorescent CGRP specifically bound to the surface of these cells and CGRP, CLR, and RAMP1 internalized into the same endosomes. CLR was prominently localized in nerve fibers of the myenteric and submucosal plexuses, muscularis externa and lamina propria of the gastrointestinal tract, and in the dorsal horn of the spinal cord of rats. CLR was detected at low levels in the soma of enteric, dorsal root ganglia (DRG), and spinal neurons. RAMP1 was also localized to enteric and DRG neurons and the dorsal horn. CLR and RAMP1 were detected in perivascular nerves and arterial smooth muscle. Nerve fibers containing CGRP and intermedin were closely associated with CLR fibers in the gastrointestinal tract and dorsal horn, and CGRP and CLR colocalized in DRG neurons. Thus, CLR and RAMP1 may mediate the effects of CGRP and intermedin in the nervous system. However, mRNA encoding RAMP2 and RAMP3 was also detected in the gastrointestinal tract, DRG, and dorsal horn, suggesting that CLR may associate with other RAMPs in these tissues to form a receptor for additional peptides such as adrenomedullin.

Localization and function of the Kv3.1b subunit in the rat medulla oblongata: focus on the nucleus tractus solitarii

The voltage-gated potassium channel subunit Kv3.1 confers fast firing characteristics to neurones. Kv3.1b subunit immunoreactivity (Kv3.1b-IR) was widespread throughout the medulla oblongata, with labelled neurones in the gracile, cuneate and spinal trigeminal nuclei. In the nucleus of the solitary tract (NTS), Kv3.1b-IR neurones were predominantly located close to the tractus solitarius (TS) and could be GABAergic or glutamatergic. Ultrastructurally, Kv3.1b-IR was detected in NTS terminals, some of which were vagal afferents. Whole-cell current-clamp recordings from neurones near the TS revealed electrophysiological characteristics consistent with the presence of Kv3.1b subunits: short duration action potentials (4.2 +/- 1.4 ms) and high firing frequencies (68.9 +/- 5.3 Hz), both sensitive to application of TEA (0.5 mm) and 4-aminopyridine (4-AP; 30 mum). Intracellular dialysis of an anti-Kv3.1b antibody mimicked and occluded the effects of TEA and 4-AP in NTS and dorsal column nuclei neurones, but not in dorsal vagal nucleus or cerebellar Purkinje cells (which express other Kv3 subunits, but not Kv3.1b). Voltage-clamp recordings from outside-out patches from NTS neurones revealed an outward K(+) current with the basic characteristics of that carried by Kv3 channels. In NTS neurones, electrical stimulation of the TS evoked EPSPs and IPSPs, and TEA and 4-AP increased the average amplitude and decreased the paired pulse ratio, consistent with a presynaptic site of action. Synaptic inputs evoked by stimulation of a region lacking Kv3.1b-IR neurones were not affected, correlating the presence of Kv3.1b in the TS with the pharmacological effects.

Effect of shear rupture on aggregate scale formation in sea ice

A discrete element model is used to study shear rupture of sea ice under convergent wind stresses. The model includes compressive, tensile, and shear rupture of viscous elastic joints connecting floes that move under the action of the wind stresses. The adopted shear rupture is governed by Coulomb’s criterion. The ice pack is a 400 km long square domain consisting of 4 km size floes. In the standard case with tensile strength 10 times smaller than the compressive strength, under uniaxial compression the failure regime is mainly shear rupture with the most probable scenario corresponding to that with the minimum failure work. The orientation of cracks delineating formed aggregates is bimodal with the peaks around the angles given by the wing crack theory determining diamond-shaped blocks. The ice block (floe aggregate) size decreases as the wind stress gradient increases since the elastic strain energy grows faster leading to a higher speed of crack propagation. As the tensile strength grows, shear rupture becomes harder to attain and compressive failure becomes equally important leading to elongation of blocks perpendicular to the compression direction and the blocks grow larger. In the standard case, as the wind stress confinement ratio increases the failure mode changes at a confinement ratio within 0.2–0.4, which corresponds to the analytical critical confinement ratio of 0.32. Below this value, the cracks are bimodal delineating diamond shape aggregates, while above this value failure becomes isotropic and is determined by small-scale stress anomalies due to irregularities in floe shape.

Do magnetospheric shear Alfvén waves generate sufficient electron energy flux to power the aurora?

Using a self-consistent drift-kinetic simulation code, we investigate whether electron acceleration owing to shear Alfvén waves in the plasma sheet boundary layer is sufficient to cause auroral brightening in the ionosphere. The free parameters used in the simulation code are guided by in situ observations of wave and plasma parameters in the magnetosphere at distances >4 RE from the Earth. For the perpendicular wavelength used in the study, which maps to ∼4 km at 110 km altitude, there is a clear amplitude threshold which determines whether magnetospheric shear Alfvén waves above the classical auroral acceleration region can excite sufficient electrons to create the aurora. Previous studies reported wave amplitudes that easily exceed this threshold; hence, the results reported in this paper demonstrate that auroral acceleration owing to shear Alfvén waves can occur in the magnetosphere at distances >4 RE from the Earth.

Electron trapping in shear Alfvén waves that power the aurora

Results from 1D Vlasov drift-kinetic plasma simulations reveal how and where auroral electrons are accelerated along Earth’s geomagnetic field. In the warm plasma sheet, electrons become trapped in shear Alfven waves, preventing immediate wave damping. As waves move to regions with larger vTe=vA, their parallel electric field decreases, and the trapped electrons escape their influence. The resulting electron distribution functions compare favorably with in situ observations, demonstrating for the first time a self-consistent link between Alfven waves and electrons that form aurora.

Timing and localization of ionospheric signatures associated with substorm expansion phase onset

In this paper, we present case studies of the optical and magnetic signatures of the characteristics of the first minute of substorm expansion phase onset observed in the ionosphere. We find that for two isolated substorms, the onset of magnetic pulsations in the 24–96 s period wavelet band are colocated in time and space with the formation and development of small-scale optical undulations along the most equatorward preexisting auroral arc prior to auroral breakup. These undulations undergo an inverse spatial cascade into vortices prior to the release of the westward traveling surge. We also present a case study of a multiple activation substorm, whereby discrete onsets of ULF wave power above a predetermined quiet time threshold are shown to be associated with specific optical intensifications and brightenings. Moreover, in the multiple activation substorm event, we show that neither the formation of the small-scale undulations nor the formation of similar structures along a north–south aligned arc is sufficient to produce auroral breakup associated with expansion phase onset. It is only ∼10 min after these two disparate activation regions initiate that auroral breakup and the subsequent formation of a westward traveling surge occur. We discuss the implications of these results in terms of the triggering mechanisms likely to be occurring during these specific events.

Rigorous bounds on the nonlinear saturation of instabilities to parallel shear flows

A novel method is presented for obtaining rigorous upper bounds on the finite-amplitude growth of instabilities to parallel shear flows on the beta-plane. The method relies on the existence of finite-amplitude Liapunov (normed) stability theorems, due to Arnol'd, which are nonlinear generalizations of the classical stability theorems of Rayleigh and Fjørtoft. Briefly, the idea is to use the finite-amplitude stability theorems to constrain the evolution of unstable flows in terms of their proximity to a stable flow. Two classes of general bounds are derived, and various examples are considered. It is also shown that, for a certain kind of forced-dissipative problem with dissipation proportional to vorticity, the finite-amplitude stability theorems (which were originally derived for inviscid, unforced flow) remain valid (though they are no longer strictly Liapunov); the saturation bounds therefore continue to hold under these conditions.

DK-1D: a drift-kinetic simulation tool for modelling the shear Alfvén wave and its interaction with collisionless plasma

We present a highly accurate tool for the simulation of shear Alfven waves (SAW) in collisionless plasma. SAW are important in space plasma environments because for small perpendicular scale lengths they can support an electric field parallel to the ambient magnetic field. Electrons can be accelerated by the parallel electric field and these waves have been implicated as the source of vibrant auroral displays. However, the parallel electric field carried by SAW is small in comparison to the perpendicular electric field of the wave, making it difficult to measure directly in the laboratory, or by satellites in the near-Earth plasma environment. In this paper, we present a simulation code that provides a means to study in detail the SAW-particle interaction in both space and laboratory plasma. Using idealised, small-amplitude propagating waves with a single perpendicular wavenumber, the simulation code accurately reproduces the damping rates and parallel electric field amplitudes predicted by linear theory for varying temperatures and perpendicular scale lengths. We present a rigorous kinetic derivation of the parallel electric field strength for small-amplitude SAW and show that commonly-used inertial and kinetic approximations are valid except for where the ratio of thermal to Alfv\'{e}n speed is between 0.7 and 1.0. We also present nonlinear simulations of large-amplitude waves and show that in cases of strong damping, the damping rates and parallel electric field strength deviate from linear predictions when wave energies are greater than only a few percent of the plasma kinetic energy, a situation which is often observed in the magnetosphere. The drift-kinetic code provides reliable, testable predictions of the parallel electric field strength which can be investigated directly in the laboratory, and will help to bridge the gap between studies of SAW in man-made and naturally occuring plasma.

On Rossby waves modified by weak sinusoidal shear

The spatial structure of beta-plane Rossby waves in a sinusoidal basic zonal flow U 0cos(γ,y) is determined analytically in the (stable) asymptotic limit of weak shear, U 0γ2 0/β≈1. The propagating neutral normal modes are found to take their greatest amplitude in the region of maximum westerly flow, while their most rapid phase variation is achieved in the region of maximum easterly flow. These results are shown to be consistent with what is obtained by ray-tracing methods in the limit of small meridional disturbance wavelength.

An exact local conservation theorem for finite-amplitude disturbances to non-parallel shear flows, with remarks on Hamiltonian structure and on Arnol'd's stability theorems

Disturbances of arbitrary amplitude are superposed on a basic flow which is assumed to be steady and either (a) two-dimensional, homogeneous, and incompressible (rotating or non-rotating) or (b) stably stratified and quasi-geostrophic. Flow over shallow topography is allowed in either case. The basic flow, as well as the disturbance, is assumed to be subject neither to external forcing nor to dissipative processes like viscosity. An exact, local ‘wave-activity conservation theorem’ is derived in which the density A and flux F are second-order ‘wave properties’ or ‘disturbance properties’, meaning that they are O(a2) in magnitude as disturbance amplitude a [rightward arrow] 0, and that they are evaluable correct to O(a2) from linear theory, to O(a3) from second-order theory, and so on to higher orders in a. For a disturbance in the form of a single, slowly varying, non-stationary Rossby wavetrain, $\overline{F}/\overline{A}$ reduces approximately to the Rossby-wave group velocity, where (${}^{-}$) is an appropriate averaging operator. F and A have the formal appearance of Eulerian quantities, but generally involve a multivalued function the correct branch of which requires a certain amount of Lagrangian information for its determination. It is shown that, in a certain sense, the construction of conservable, quasi-Eulerian wave properties like A is unique and that the multivaluedness is inescapable in general. The connection with the concepts of pseudoenergy (quasi-energy), pseudomomentum (quasi-momentum), and ‘Eliassen-Palm wave activity’ is noted. The relationship of this and similar conservation theorems to dynamical fundamentals and to Arnol'd's nonlinear stability theorems is discussed in the light of recent advances in Hamiltonian dynamics. These show where such conservation theorems come from and how to construct them in other cases. An elementary proof of the Hamiltonian structure of two-dimensional Eulerian vortex dynamics is put on record, with explicit attention to the boundary conditions. The connection between Arnol'd's second stability theorem and the suppression of shear and self-tuning resonant instabilities by boundary constraints is discussed, and a finite-amplitude counterpart to Rayleigh's inflection-point theorem noted

Localization of deep water formation: role of atmospheric moisture transport and geometrical constraints on ocean circulation

A series of coupled atmosphere–ocean–ice aquaplanet experiments is described in which topological constraints on ocean circulation are introduced to study the role of ocean circulation on the mean climate of the coupled system. It is imagined that the earth is completely covered by an ocean of uniform depth except for the presence or absence of narrow barriers that extend from the bottom of the ocean to the sea surface. The following four configurations are described: Aqua (no land), Ridge (one barrier extends from pole to pole), Drake (one barrier extends from the North Pole to 35°S), and DDrake (two such barriers are set 90° apart and join at the North Pole, separating the ocean into a large basin and a small basin, connected to the south). On moving from Aqua to Ridge to Drake to DDrake, the energy transports in the equilibrium solutions become increasingly “realistic,” culminating in DDrake, which has an uncanny resemblance to the present climate. Remarkably, the zonal-average climates of Drake and DDrake are strikingly similar, exhibiting almost identical heat and freshwater transports, and meridional overturning circulations. However, Drake and DDrake differ dramatically in their regional climates. The small and large basins of DDrake exhibit distinctive Atlantic-like and Pacific-like characteristics, respectively: the small basin is warmer, saltier, and denser at the surface than the large basin, and is the main site of deep water formation with a deep overturning circulation and strong northward ocean heat transport. A sensitivity experiment with DDrake demonstrates that the salinity contrast between the two basins, and hence the localization of deep convection, results from a deficit of precipitation, rather than an excess of evaporation, over the small basin. It is argued that the width of the small basin relative to the zonal fetch of atmospheric precipitation is the key to understanding this salinity contrast. Finally, it is argued that many gross features of the present climate are consequences of two topological asymmetries that have profound effects on ocean circulation: a meridional asymmetry (circumpolar flow in the Southern Hemisphere; blocked flow in the Northern Hemisphere) and a zonal asymmetry (a small basin and a large basin).

Chaotic destruction of Anderson localization in a nonlinear lattice

We consider a scattering problem for a nonlinear disordered lattice layer governed by the discrete nonlinear Schrodinger equation. The linear state with exponentially small transparency, due to the Anderson localization, is followed for an increasing nonlinearity, until it is destroyed via a bifurcation. The critical nonlinearity is shown to decay with the lattice length as a power law. We demonstrate that in the chaotic regimes beyond the bifurcation the field is delocalized and this leads to a drastic increase of transparency. Copyright (C) EPLA, 2008