Error analysis of Trefftz-discontinuous Galerkin methods for the time-harmonic Maxwell equations

In this paper, we extend to the time-harmonic Maxwell equations the p-version analysis technique developed in [R. Hiptmair, A. Moiola and I. Perugia, Plane wave discontinuous Galerkin methods for the 2D Helmholtz equation: analysis of the p-version, SIAM J. Numer. Anal., 49 (2011), 264-284] for Trefftz-discontinuous Galerkin approximations of the Helmholtz problem. While error estimates in a mesh-skeleton norm are derived parallel to the Helmholtz case, the derivation of estimates in a mesh-independent norm requires new twists in the duality argument. The particular case where the local Trefftz approximation spaces are built of vector-valued plane wave functions is considered, and convergence rates are derived.

Dissipation of shear-free turbulence near boundaries

The rapid-distortion model of Hunt & Graham (1978) for the initial distortion of turbulence by a flat boundary is extended to account fully for viscous processes. Two types of boundary are considered: a solid wall and a free surface. The model is shown to be formally valid provided two conditions are satisfied. The first condition is that time is short compared with the decorrelation time of the energy-containing eddies, so that nonlinear processes can be neglected. The second condition is that the viscous layer near the boundary, where tangential motions adjust to the boundary condition, is thin compared with the scales of the smallest eddies. The viscous layer can then be treated using thin-boundary-layer methods. Given these conditions, the distorted turbulence near the boundary is related to the undistorted turbulence, and thence profiles of turbulence dissipation rate near the two types of boundary are calculated and shown to agree extremely well with profiles obtained by Perot & Moin (1993) by direct numerical simulation. The dissipation rates are higher near a solid wall than in the bulk of the flow because the no-slip boundary condition leads to large velocity gradients across the viscous layer. In contrast, the weaker constraint of no stress at a free surface leads to the dissipation rate close to a free surface actually being smaller than in the bulk of the flow. This explains why tangential velocity fluctuations parallel to a free surface are so large. In addition we show that it is the adjustment of the large energy-containing eddies across the viscous layer that controls the dissipation rate, which explains why rapid-distortion theory can give quantitatively accurate values for the dissipation rate. We also find that the dissipation rate obtained from the model evaluated at the time when the model is expected to fail actually yields useful estimates of the dissipation obtained from the direct numerical simulation at times when the nonlinear processes are significant. We conclude that the main role of nonlinear processes is to arrest growth by linear processes of the viscous layer after about one large-eddy turnover time.

Homoleptic copper(II) and copper(I) complexes of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline. Six-coordinate copper(I) in solution

Reaction of 5,6-dihydro-5,6-epoxy-1,10-phenanthroline (L) with Cu(ClO(4))(2)center dot 6H(2)O in methanol in 3:1 M ratio at room temperature yields light green [CuL(3)](ClO(4))(2)center dot H(2)O (1). The X-ray crystal structure of the hemi acetonitrile solvate [CuL(3)](ClO(4))(2)center dot 0.5CH(3)CN has been determined which shows Jahn-Teller distortion in the CuN(6) core present in the cation [CuL(3)](2+). Complex 1 gives an axial EPR spectrum in acetonitrile-toluene glass with g(parallel to) = 2.262 (A(parallel to) = 169 x 10 (4) cm (1)) and g(perpendicular to) = 2.069. The Cu(II/I) potential in 1 in CH(2)Cl(2) at a glassy carbon electrode is 0.32 V versus NHE. This potential does not change with the addition of extra L in the medium implicating generation of a six-coordinate copper(I) species [CuL(3)](+) in solution. B3LYP/LanL2DZ calculations show that the six Cu-N bond distances in [CuL(3)](+) are 2.33, 2.25, 2.32, 2.25, 2.28 and 2.25 angstrom while the ideal Cu(I)-N bond length in a symmetric Cu(I)N(6) moiety is estimated as 2.25 angstrom. Reaction of L with Cu(CH(3)CN)(4)ClO(4) in dehydrated methanol at room temperature even in 4:1 M proportion yields [CuL(2)]ClO(4) (2). Its (1)H NMR spectrum indicates that the metal in [CuL(2)](+) is tetrahedral. The Cu(II/I) potential in 2 is found to be 0.68 V versus NHE in CH(2)Cl(2) at a glassy carbon electrode. In presence of excess L, 2 yields the cyclic voltammogram of 1. From (1)H NMR titration, the free energy of binding of L to [CuL(2)](+) to produce [CuL(3)](+) in CD(2)Cl(2) at 298 K is estimated as -11.7 (+/-0.2) kJ mol (1).

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.

Effects of non-uniform crosswind fields on scintillometry measurements

The effects of a non-uniform wind field along the path of a scintillometer are investigated. Theoretical spectra are calculated for a range of scenarios where the crosswind varies in space or time and compared to the ‘ideal’ spectrum based on a constant uniform crosswind. It is verified that the refractive-index structure parameter relation with the scintillometer signal remains valid and invariant for both spatially and temporally-varying crosswinds. However, the spectral shape may change significantly preventing accurate estimation of the crosswind speed from the peak of the frequency spectrum and retrieval of the structure parameter from the plateau of the power spectrum. On comparison with experimental data, non-uniform crosswind conditions could be responsible for previously unexplained features sometimes seen in observed spectra. By accounting for the distribution of crosswind, theoretical spectra can be generated that closely replicate the observations, leading to a better understanding of the measurements. Spatial variability of wind speeds should be expected for paths other than those that are parallel to the surface and over flat, homogenous areas, whilst fluctuations in time are important for all sites.

Granular flow in the marginal ice zone

The region of sea ice near the edge of the sea ice pack is known as the marginal ice zone (MIZ), and its dynamics are complicated by ocean wave interaction with the ice cover, strong gradients in the atmosphere and ocean and variations in sea ice rheology. This paper focuses on the role of sea ice rheology in determining the dynamics of the MIZ. Here, sea ice is treated as a granular material with a composite rheology describing collisional ice floe interaction and plastic interaction. The collisional component of sea ice rheology depends upon the granular temperature, a measure of the kinetic energy of flow fluctuations. A simplified model of the MIZ is introduced consisting of the along and across momentum balance of the sea ice and the balance equation of fluctuation kinetic energy. The steady solution of these equations is found to leading order using elementary methods. This reveals a concentrated region of rapid ice flow parallel to the ice edge, which is in accordance with field observations, and previously called the ice jet. Previous explanations of the ice jet relied upon the existence of ocean currents beneath the ice cover. We show that an ice jet results as a natural consequence of the granular nature of sea ice.

Experimental confirmation of transformation pathways between inverse double diamond and gyroid cubic phases

A macroscopically oriented double diamond inverse bicontinuous cubic phase (QIID) of the lipid glycerol monooleate is reversibly converted into a gyroid phase (QIIG). The initial QIID phase is prepared in the form of a film coating the inside of a capillary, deposited under flow, which produces a sample uniaxially oriented with a ⟨110⟩ axis parallel to the symmetry axis of the sample. A transformation is induced by replacing the water within the capillary tube with a solution of poly(ethylene glycol), which draws water out of the QIID sample by osmotic stress. This converts the QIID phase into a QIIG phase with two coexisting orientations, with the ⟨100⟩ and ⟨111⟩ axes parallel to the symmetry axis, as demonstrated by small-angle X-ray scattering. The process can then be reversed, to recover the initial orientation of QIID phase. The epitaxial relation between the two oriented mesophases is consistent with topologypreserving geometric pathways that have previously been hypothesized for the transformation. Furthermore, this has implications for the production of macroscopically oriented QIIG phases, in particular with applications as nanomaterial templates.

Dynamic simulation studies of structural formation and transition in electro-magneto-rheological fluids

The computer simulation method has been used to study the structural formation and transition of electro-magneto-rheological (EMR) fluids under compatible electric and magnetic fields. When the fields are applied simultaneously and perpendicularly to each other, the particles rapidly arrange into two-dimensional close-packed layer structures parallel to both fields. The layers then combine together to form thicker sheet-like structures, which finally relax into three-dimensional close-packed structures with the help of the thermal fluctuations. On the other hand, if the electric field is applied firstly to induce the body-centered tetragonal (BCT) columns in the system, and then the magnetic field is applied in the perpendicular direction. the BCT to face-centered cubic (FCC) structure transition is observed in very short time. Following that. the structure keeps on evolving due to the demagnetization effect and finally form the three-dimensional close-packed structures.

Squeeze flow viscoelasticity of electrorheological fluids based on microcrystalline cellulose

Dynamic viscoelasticity of electrorheological fluids based on microcrystalline cellulose/castor oil suspensions was experimentally investigated in squeeze flow. The dependence of storage modulus G' and loss modulus G" parallel to external electric field on electric fields and strain amplitudes is presented. The experiments show that, when external electric field is higher than the critical field, the viscoelasticity of the ER fluids converts from linear to nonlinear, and the ER fluids transfer from solid-like state to fluid state with the growth of strain amplitude. The influences of strain amplitude and oscillatory frequency on the nonlinearity of viscoelasticity were also studied.

Simulation of field-induced structural formation and transition in electromagnetorheological suspensions

A computer simulation method has been used to study the three-dimensional structural formation and transition of eleetromagnetorheological (EMR) suspensions under compatible electric and magnetic fields. When the fields are applied simultaneously and perpendicularly to each other, the particles rapidly arrange into single layer structures parallel to both fields. In each layer, there is a two-dimensional hexagonal lattice. The single layers then combine together to form thicker sheetlike structures. With the help of the thermal fluctuations, the thicker structures relax into three-dimensional close-packed structures, which may be face-centered cubic (fcc), hexagonal close-packed (hup) lattices, or, more probably, the mixture of them, depending on the initial configurations and the thermal fluctuations. On the other hand, if the electric field is applied first to induce the body-centered tetragonal (bct) columns in the system, and then the magnetic field is applied in the perpendicular direction, the bet to fee structure transition is observed in a very short time. Following that, the structure keeps on evolving due to the demagnetization effect and finally forms close-packed structures with fee and hcp lattice character. The simulation results are in agreement with the theoretical and experimental 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.

Evidence of component merging equatorward of the cusp

The Polar spacecraft passed through a region near the dayside magnetopause on May 29, 1996, at a geocentric distance of similar to 8 R-E and high, northern magnetic latitudes. The interplanetary magnetic field (IMF) was northward during the pass. Data from the Thermal Ion Dynamics Experiment revealed the existence of low-speed (similar to 50 km s(-1)) ion D-shaped distributions mixed with cold ions (similar to 2 eV) over a period of 2.5 hours. These ions were traveling parallel to the magnetic field toward the Northern Hemisphere ionosphere and were convecting primarily eastward. The D-shaped distributions are distinct from a convecting Maxwellian and, along with the magnetic field direction, are taken as evidence that the spacecraft was inside the magnetosphere and not in the magnetosheath. Furthermore, the absence of ions in the antiparallel direction is taken as evidence that low-shear merging was occurring at a location southward of the spacecraft and equatorward of the Southern Hemisphere cusp. The cold ions were of ionospheric origin, with initially slow field-aligned speeds, which were accelerated upon reflection from the magnetopause. These observations provide significant new evidence consistent with component magnetic merging sites equatorward of the cusp for northward IMF.

Multi-instrument ground-based observations of a travelling convection vortices event

A coordinated ground-based observational campaign using the IMAGE magnetometer network, EISCAT radars and optical instruments on Svalbard has made possible detailed studies of a travelling convection vortices (TCV) event on 6 January 1992. Combining the data from these facilities allows us to draw a very detailed picture of the features and dynamics of this TCV event. On the way from the noon to the drawn meridian, the vortices went through a remarkable development. The propagation velocity in the ionosphere increased from 2.5 to 7.4 km s−1, and the orientation of the major axes of the vortices rotated from being almost parallel to the magnetic meridian near noon to essentially perpendicular at dawn. By combining electric fields obtained by EISCAT and ionospheric currents deduced from magnetic field recordings, conductivities associated with the vortices could be estimated. Contrary to expectations we found higher conductivities below the downward field aligned current (FAC) filament than below the upward directed. Unexpected results also emerged from the optical observations. For most of the time there were no discrete aurora at 557.7 nm associated with the TCVs. Only once did a discrete form appear at the foot of the upward FAC. This aurora subsequently expanded eastward and westward leaving its centre at the same longitude while the TCV continued to travel westward. Also we try to identify the source regions of TCVs in the magnetosphere and discuss possible generation mechanisms.

Analysis of incoherent scatter radar data from non-thermal F-region plasma

A procedure is presented for fitting incoherent scatter radar data from non-thermal F-region ionospheric plasma, using theoretical spectra previously predicted. It is found that values of the shape distortion factor D∗, associated with deviations of the ion velocity distribution from a Maxwellian distribution, and ion temperatures can be deduced (the results being independent of the path of iteration) if the angle between the line-of-sight and the geomagnetic field is larger than about 15–20°. The procedure can be used with one or both of two sets of assumptions. These concern the validity of the adopted model for the line-of-sight ion velocity distribution in the one case or for the full three-dimensional ion velocity distribution function in the other. The distribution function employed was developed to describe the line-of-sight velocity distribution for large aspect angles, but both experimental data and Monte Carlo simulations indicate that the form of the field-perpendicular distribution can also describe the distribution at more general aspect angles. The assumption of this form for the line-of-sight velocity distribution at a general aspect angle enables rigorous derivation of values of the one-dimensional, line-of-sight ion temperature. With some additional assumptions (principally that the field-parallel distribution is always Maxwellian and there is a simple relationship between the ion temperature anisotropy and the distortion of the field-perpendicular distribution from a Maxwellian), fits to data for large aspect angles enable determination of line-of-sight temperatures at all aspect angles and hence, of the average ion temperature and the ion temperature anisotropy. For small aspect angles, the analysis is restricted to the determination of the line-of-sight ion temperature because the theoretical spectrum is insensitive to non-thermal effects when the plasma is viewed along directions almost parallel to the magnetic field. This limitation is expected to apply to any realistic model of the ion velocity distribution function and its consequences are discussed. Fit strategies which allow for mixed ion composition are also considered. Examples of fits to data from various EISCAT observing programmes are presented.

Analysis of incoherent scatter spectra from non-Maxwellian plasma

Incoherent scatter data from non-thermal F-region ionospheric plasma are analysed, using theoretical spectra predicted by Raman et al. It is found that values of the semi-empirical drift parameter D∗, associated with deviations of the ion velocity distribution from a Maxwellian, and the plasma temperatures can be rigorously deduced (the results being independent of the path of iteration) if the angle between the line-of-sight and the geomagnetic field is larger than about 15–20 degrees. For small aspect angles, the deduced value of the average (or 3-D) ion temperature remains ambiguous and the analysis is restricted to the determination of the line-of-sight temperature because the theoretical spectrum is insensitive to non-thermal effects when the plasma is viewed along directions almost parallel to the magnetic field. This limitation is expected to apply to any realistic model of the ion velocity distribution, and its consequences are discussed. Fit strategies which allow for mixed ion composition are also considered. Examples of fits to data from various EISCAT observing programmes are presented.