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

The gravity wave momentum flux in hydrostatic flow with directional shear over elliptical mountains

Semi-analytical expressions for the momentum flux associated with orographic internal gravity waves, and closed analytical expressions for its divergence, are derived for inviscid, stationary, hydrostatic, directionally-sheared flow over mountains with an elliptical horizontal cross-section. These calculations, obtained using linear theory conjugated with a third-order WKB approximation, are valid for relatively slowly-varying, but otherwise generic wind profiles, and given in a form that is straightforward to implement in drag parametrization schemes. When normalized by the surface drag in the absence of shear, a quantity that is calculated routinely in existing drag parametrizations, the momentum flux becomes independent of the detailed shape of the orography. Unlike linear theory in the Ri → ∞ limit, the present calculations account for shear-induced amplification or reduction of the surface drag, and partial absorption of the wave momentum flux at critical levels. Profiles of the normalized momentum fluxes obtained using this model and a linear numerical model without the WKB approximation are evaluated and compared for two idealized wind profiles with directional shear, for different Richardson numbers (Ri). Agreement is found to be excellent for the first wind profile (where one of the wind components varies linearly) down to Ri = 0.5, while not so satisfactory, but still showing a large improvement relative to the Ri → ∞ limit, for the second wind profile (where the wind turns with height at a constant rate keeping a constant magnitude). These results are complementary, in the Ri > O(1) parameter range, to Broad’s generalization of the Eliassen–Palm theorem to 3D flow. They should contribute to improve drag parametrizations used in global weather and climate prediction models.

Dynamic response times of electrorheological fluids in steady shear

Transient responses of electrorheological fluids to square-wave electric fields in steady shear are investigated by computational simulation method. The structure responses of the fluids to the field with high frequency are found to be very similar to that to the field with very low frequency or the sudden applied direct current field. The stress rise processes are also similar in both cases and can be described by an exponential expression. The characteristic time tau of the stress response is found to decrease with the increase of the shear rate (gamma) over dot and the area fraction of the particles phi(2). The relation between them can be roughly expressed as tau proportional to(gamma) over dot(-3/4)phi(2)(-3/2). The simulation results are compared with experimental measurements. The aggregation kinetics of the particles in steady shear is also discussed according to these results.

Structure and viscoelasticity of an electrorheological fluid in oscillatory shear: computer simulation investigation

The three-dimensional molecular dynamics simulation method has been used to study the dynamic responses of an electrorheological (ER) fluid in oscillatory shear. The structure and related viscoelastic behaviour of the fluid are found to be sensitive to the amplitude of the strain. With the increase of the strain amplitude, the structure formed by the particles changes from isolated columns to sheet-like structures which may be perpendicular or parallel to the oscillating direction. Along with the structure evolution, the field-induced moduli decrease significantly with an increase in strain amplitude. The viscoelastic behaviour of the structures obtained in the cases of different strain amplitudes was examined in the linear response regime and an evident structure dependence of the moduli was found. The reason for this lies in the anisotropy of the arrangement of the particles in these structures. Short-range interactions between the particles cannot be neglected in determining the viscoelastic behaviour of ER fluids at small strain amplitude, especially for parallel sheets. The simulation results were compared with available experimental data and good agreement was reached for most of them.

The effect of electric field direction on the-shear stress of ER fluids

The effect of the direction of external electric field on the shear stress of an ER fluid has been studied by molecular-dynamics simulation. Due to the formation of inclined chains, the shear stress strongly depends on the direction of the field, and it may be very large under some special field direction. And theoretical model of ideal microstructure of ER fluids has proved this result. Thus the ER effect may be greatly enhanced just by choosing an optimum direction for the field without any additional requirement, suggesting a promising way to the practical application of ER fluids.

Polymer brushes under shear: molecular dynamics simulations compared to experiments

Surfaces coated with polymer brushes in a good solvent are known to exhibit excellent tribological properties. We have performed coarse-grained equilibrium and nonequilibrium molecular dynamics (MD) simulations to investigate dextran polymer brushes in an aqueous environment in molecular detail. In a first step, we determined simulation parameters and units by matching experimental results for a single dextran chain. Analyzing this model when applied to a multichain system, density profiles of end-tethered polymer brushes obtained from equilibrium MD simulations compare very well with expectations based on self-consistent field theory. Simulation results were further validated against and correlated with available experimental results. The simulated compression curves (normal force as a function of surface separation) compare successfully with results obtained with a surface forces apparatus. Shear stress (friction) obtained via nonequilibrium MD is contrasted with nanoscale friction studies employing colloidal-probe lateral force microscopy. We find good agreement in the hydrodynamic regime and explain the observed leveling-off of the friction forces in the boundary regime by means of an effective polymer–wall attraction.

Influence of inherent structure shear stress of supercooled liquids on their shear moduli

Configurations of supercooled liquids residing in their local potential minimum (i.e. in their inherent structure, IS) were found to support a non-zero shear stress. This IS stress was attributed to the constraint to the energy minimization imposed by boundary conditions, which keep size and shape of the simulation cell fixed. In this paper we further investigate the influence of these boundary conditions on the IS stress. We investigate its importance for the computation of the low frequency shear modulus of a glass obtaining a consistent picture for the low- and high frequency shear moduli over the full temperature range. Hence, we find that the IS stress corresponds to a non-thermal contribution to the fluctuation term in the Born-Green expression. This leads to an unphysical divergence of the moduli in the low temperature limit if no proper correction for this term is applied. Furthermore, we clarify the IS stress dependence on the system size and put its origin on a more formal basis.

Science for Humanitarian Emergencies and Resilience (SHEAR) scoping study: Annex 3 - Early warning system and risk assessment case studies

This report provides case studies of Early Warning Systems (EWSs) and risk assessments encompassing three main hazard types: drought; flood and cyclone. The case studies are taken from ten countries across three continents (focusing on Africa, South Asia and the Caribbean). The case studies have been developed to assist the UK Department for International Development (DFID) to prioritise areas for Early Warning System (EWS) related research under their ‘Science for Humanitarian Emergencies and Resilience’ (SHEAR) programme. The aim of these case studies is to ensure that DFID SHEAR research is informed by the views of Non-Governmental Organisations (NGOs) and communities engaged with Early Warning Systems and risk assessments (including community-based Early Warning Systems). The case studies highlight a number of challenges facing Early Warning Systems (EWSs). These challenges relate to financing; integration; responsibilities; community interpretation; politics; dissemination; accuracy; capacity and focus. The case studies summarise a number of priority areas for EWS related research: • Priority 1: Contextualising and localising early warning information • Priority 2: Climate proofing current EWSs • Priority 3: How best to sustain effective EWSs between hazard events? • Priority 4: Optimising the dissemination of risk and warning information • Priority 5: Governance and financing of EWSs • Priority 6: How to support EWSs under challenging circumstances • Priority 7: Improving EWSs through monitoring and evaluating the impact and effectiveness of those systems

CLEC-2 is not required for platelet aggregation at arteriolar shear

The C-type lectin receptor CLEC-2 is expressed primarily on the surface of platelets, where it is present as a dimer, and is found at low level on a subpopulation of other hematopoietic cells, including mouse neutrophils [1–4] Clustering of CLEC-2 by the snake venom toxin rhodocytin, specific antibodies or its endogenous ligand, podoplanin, elicits powerful activation of platelets through a pathway that is similar to that used by the collagen receptor glycoprotein VI (GPVI) [4–6]. The cytosolic tail of CLEC-2 contains a conserved YxxL sequence preceded by three upstream acidic amino acid residues, which together form a novel motif known as a hemITAM. Ligand engagement induces tyrosine phosphorylation of the hemITAM sequence providing docking sites for the tandem-SH2 domains of the tyrosine kinase Syk across a CLEC-2 receptor dimer [3]. Tyrosine phosphorylation of Syk by Src family kinases and through autophosphorylation leads to stimulation of a downstream signaling cascade that culminates in activation of phospholipase C γ2 (PLCγ2) [4,6]. Recently, CLEC-2 has been proposed to play a major role in supporting activation of platelets at arteriolar rates of flow [1]. Injection of a CLEC-2 antibody into mice causes a sustained depletion of the C-type lectin receptor from the platelet surface [1]. The CLEC-2-depleted platelets were unresponsive to rhodocytin but underwent normal aggregation and secretion responses after stimulation of other platelet receptors, including GPVI [1]. In contrast, there was a marked decrease in aggregate formation relative to controls when CLEC-2-depleted blood was flowed at arteriolar rates of shear over collagen (1000 s−1 and 1700 s−1) [1]. Furthermore, antibody treatment significantly increased tail bleeding times and mice were unable to occlude their vessels after ferric chloride injury [1]. These data provide evidence for a critical role for CLEC-2 in supporting platelet aggregation at arteriolar rates of flow. The underlying mechanism is unclear as platelets do not express podoplanin, the only known endogenous ligand of CLEC-2. In the present study, we have investigated the role of CLEC-2 in platelet aggregation and thrombus formation using platelets from a novel mutant mouse model that lacks functional CLEC-2.

Simulating startup shear of entangled polymer melts

Start-up shear rheology is a standard experiment used for characterizing polymer flow, and to test various models of polymer dynamics. A rich phenomenology is developed for behavior of entangled monodisperse linear polymers in such tests, documenting shear stress overshoots as a function of shear rates and molecular weights. A tube theory does a reasonable qualitative job at describing these phenomena, although it involves several drastic approximations and the agreement can be fortuitous. Recently, Lu and coworkers published several papers [e.g. Lu {\it et al.} {\it ACS Macro Lett}. 2014, 3, 569-573] reporting results from molecular dynamics simulations of linear entangled polymers, which contradict both theory and experiment. Based on these observations, they made very serious conclusions about the tube theory, which seem to be premature. In this letter, we repeat simulations of Lu {\it et al.} and systematically show that neither their simulation results, nor their comparison with theory are confirmed.

The utility of convection-permitting ensembles for the prediction of stationary convective bands

This study examines convection-permitting numerical simulations of four cases of terrain-locked quasi-stationary convective bands over the UK. For each case, a 2.2-km grid-length 12-member ensemble and 1.5-km grid-length deterministic forecast are analyzed, each with two different initialization times. Object-based verification is applied to determine whether the simulations capture the structure, location, timing, intensity and duration of the observed precipitation. These verification diagnostics reveal that the forecast skill varies greatly between the four cases. Although the deterministic and ensemble simulations captured some aspects of the precipitation correctly in each case, they never simultaneously captured all of them satisfactorily. In general, the models predicted banded precipitation accumulations at approximately the correct time and location, but the precipitating structures were more cellular and less persistent than the coherent quasi-stationary bands that were observed. Ensemble simulations from the two different initialization times were not significantly different, which suggests a potential benefit of time-lagging subsequent ensembles to increase ensemble size. The predictive skill of the upstream larger-scale flow conditions and the simulated precipitation on the convection-permitting grids were strongly correlated, which suggests that more accurate forecasts from the parent ensemble should improve the performance of the convection-permitting ensemble nested within it.

Intensities and self-broadening coefficients of the strongest water vapour lines in the 2.7 and 6.25μm absorption bands

Intensities and self-broadening coefficients are presented for about 460 of the strongest water vapour lines in the spectral regions 1400–1840 cm−1 and 3440–3970 cm−1 at room temperature, obtained from rather unique measurements using a 5-mm-path-length cell. The retrieved spectral line parameters are compared with those in the HITRAN database ver. 2008 and 2012 and with recent ab-initio calculations. Both the retrieved intensities and half-widths are on average in reasonable agreement with those in HITRAN-2012. Maximum systematic differences do not exceed 4% for intensities (1600 cm−1 band) and 7% for self-broadening coefficients (3600 cm−1 band). For many lines however significant disagreements were detected with the HITRAN-2012 data, exceeding the average uncertainty of the retrieval. In addition, water vapour line parameters for 5300 cm−1 (1.9 μm) band reported by us in 2005 were also compared with HITRAN-2012, and show average differences of 4–5% for both intensities and half-widths.

Effect of blood contamination on the shear bond strength at resin/dentin interface in primary teeth

Purpose: To assess in vitro the shear bond strength at the resin/dentin interface in primary teeth after contamination with fresh human blood. Methods: 75 crowns of primary molars were embedded in acrylic resin and mechanically ground to expose a flat dentin surface. The specimens were randomly assigned to five groups (n=15), according to the surface treatment. Group I (control) had no blood contamination. The other groups were blood-contaminated and subjected to different post-contamination protocols: in Group 2, the surfaces were rinsed with water; in Group 3, the surfaces were air-dried; in Group 4, the surfaces were rinsed and air-dried; and in Group 5, no post-contamination treatment was done. In all groups, a 3-mm dentin bonding site was demarcated, Single Bond adhesive system was applied and resin composite cylinders were bonded. After 24 hours in distilled water, shear bond strength was tested at a crosshead speed of 0.5 mm/minute. Results: Means (in MPa) were: Group 1: 7.1 (+/- 4.2); Group 2: 4.0 (+/- 1.8); Group 3: 0.9 (+/- 0.7); Group 4: 3.9 (+/- 2.2) and Group 5: 1.3 (+/- 1.5). Data were analyzed statistically by the Kruskal-Wallis test at 5% significance level. Groups 2 and 4 were similar to each other (P > 0.05) and both ware similar to Group 1 (P > 0.05). These groups (2, 3 and 4) had statistically significantly higher bond strengths than Groups 3 and 5 (P < 0.05). Blood contamination negatively affected the shear bond strength to primary tooth dentin. Among the blood-contaminated groups, water-rinsed specimens had higher bond strengths than those that were exclusively air-dried or not submitted to any post-contamination protocol before adhesive application.