Solving the advection-dispersion equation using the discrete puff particle method

A flexible, mass-conservative numerical technique for solving the advection-dispersion equation for miscible contaminant transport is presented. The method combines features of puff transport models from air pollution studies with features from the random walk particle method used in water resources studies, providing a deterministic time-marching algorithm which is independent of the grid Peclet number and scales from one to higher dimensions simply. The concentration field is discretised into a number of particles, each of which is treated as a point release which advects and disperses over the time interval. The dispersed puff is itself discretised into a spatial distribution of particles whose masses can be pre-calculated. Concentration within the simulation domain is then calculated from the mass distribution as an average over some small volume. Comparison with analytical solutions for a one-dimensional fixed-duration concentration pulse and for two-dimensional transport in an axisymmetric flow field indicate that the algorithm performs well. For a given level of accuracy the new method has lower computation times than the random walk particle method.

Surface plasmons in layered structures with semiconductor and metallic films

The complete spectrum of eigenwaves including surface plasmon polaritons (SPP), dynamic (bulk) and complex waves in the layered structures containing semiconductor and metallic films has been explored. The effects of loss, geometry and the parameters of dielectric layers on the eigenmode spectrum and, particularly, on the SPP modes have been analysed using both the asymptotic and rigorous numerical solutions of the full-wave dispersion equation. The field and Poynting vector distributions have been examined to identify the modes and elucidate their properties. It has been shown that losses and dispersion of permittivity qualitatively alter the spectral content and the eigenwave properties. The SPP counter-directional power fluxes in the film and surrounding dielectrics have been attributed to vortices of power flow, which are responsible for the distinctive features of SPP modes. It has been demonstrated for the first time that the maximal attainable slow-wave factor of the SPP modes guided by thin Au films at optical frequencies is capped not by losses but the frequency dispersion of the actual Au permittivity. © 2009 EDP Sciences.

A Modified Pole-Zero Technique for the Synthesis of Waveguide Leaky-Wave Antennas Loaded With Dipole-Based FSS

An extension of the pole-zero matching method proposed by Stefano Maci et al. for the analysis of electromagnetic bandgap (EBG) structures composed by lossless dipole-based frequency selective surfaces (FSS) printed on stratified dielectric media, is presented in this paper. With this novel expansion, the dipoles length appears as a variable in the analytical dispersion equation. Thus, modal dispersion curves as a function of the dipoles length can be easily obtained with the only restriction of single Floquet mode propagation. These geometry-dispersion curves are essential for the efficient analysis and design of practical EBG structures, such as waveguides loaded with artificial magnetic conductors (AMC) for miniaturization, or leaky-wave antennas (LWA) using partially reflective surfaces (PRS). These two practical examples are examined in this paper. Results are compared with full-wave 2D and 3D simulations showing excellent agreement, thus validating the proposed technique and illustrating its utility for practical designs.

Determining the Reaeration Coefficient and Hydrodynamic Properties of Rivers Using Inert Gas Tracers

Various contaminants which can be aerobically degraded find their way directly or indirectly into surface water bodies. The reaeration coefficient (K2) characterises the rate at which oxygen can transfer from the atmosphere across the air-water interface following oxygen depletion in a water body. Other mechanisms (like advection, dispersion and transient storage) determine how quickly the contaminants can spread in the water, affecting their spatial and temporal concentrations. Tracer methods involving injection of a gas into the water body have traditionally been used for direct (in-situ) measurement of K2 in a given reach. This paper shows how additional modelling of tracer test results can be used to quantify also hydrodynamic mechanisms (e.g. dispersion and storage exchange coefficients, etc.). Data from three tracer tests conducted in the River Lagan (Northern Ireland) using an inert gas (krypton, Kr) are re-analysed using two solute transport models (ADM, TSM) and an inverse-modelling framework (OTIS-P). Results for K2 are consistent with previously published values for this reach (K2(20)~10-40 d-1). The storage area constituted 30-60% of the main cross-section area and the storage exchange rate was between 2.5×10-3-3.2×10-3s-1. The additional hydrodynamic parameters obtained give insight into transport and dispersion mechanisms within the reach.

Kelvin-Helmholtz instability in solar chromospheric jets: theory and observation

Using data obtained by the high-resolution CRisp Imaging SpectroPolarimeter instrument on the Swedish 1 m Solar Telescope, we investigate the dynamics and stability of quiet-Sun chromospheric jets observed at the disk center. Small-scale features, such as rapid redshifted and blueshifted excursions, appearing as high-peed jets in the wings of the Hα line, are characterized by short lifetimes and rapid fading without any descending behavior. To study the theoretical aspects of their stability without considering their formation mechanism, we model chromospheric jets as twisted magnetic flux tubes moving along their axis, and use the ideal linear incompressible magnetohydrodynamic approximation to derive the governing dispersion equation. Analytical solutions of the dispersion equation indicate that this type of jet is unstable to Kelvin–Helmholtz instability (KHI), with a very short (few seconds) instability growth time at high upflow speeds. The generated vortices and unresolved turbulent flows associated with the KHI could be observed as a broadening of chromospheric spectral lines. Analysis of the Hα line profiles shows that the detected structures have enhanced line widths with respect to the background. We also investigate the stability of a larger-scale Hα jet that was ejected along the line of sight. Vortex-like features, rapidly developing around the jet’s boundary, are considered as evidence of the KHI. The analysis of the energy equation in the partially ionized plasma shows that ion–neutral collisions may lead to fast heating of the KH vortices over timescales comparable to the lifetime of chromospheric jets.

A Correlation between the Dynamic Viscosity and Crystallisation of Felodipine from its Amorphous Polymer Solid Dispersion

Purpose Previously, it has been reported that molecular mobility determines the rate of molecular approach to crystal surfaces, while entropy relates to the probability of that approaching molecule having the desirable configuration for further growth of the existing crystal; and the free energy dictates the probability of that molecule not returning to the liquid phase1. If we plot the crystal growth rate and viscosity of a supercooled liquid in a log-log format, the relationship between the two is linear, indicating the influence viscosity has upon crystal growth rate. However, such approximation has been derived from pure drug compounds and it is apparent that further understanding of crystallization from drug-polymer solid dispersion is required in order to stabilise drugs embedded within amorphous polymeric solid dispersions. Methods Mixtures of felodipine and polymer (HPMCAS-HF, PVPK15 and Soluplus®) at specified compositions were prepared using a Restch MM200 ball mill. To examine crystal growth within amorphous solid dispersions, samples were prepared by melting 5-10 mg of ball milled mixture at 150°C for 3-5 minutes on a glass slip pre-cleaned with methanol and acetone. All prepared samples were confirmed to be crystal free by visual observation using a polarised light microscope (Olympus BX50). Prepared samples were stored at 0% RH (P2O5), inside desiccators, maintained in ovens at 80°C. For the dynamic viscosity measurement, approximately 100-200mg ball milled mixture was heated on the base plate of a rotational rheometer at 150°C for 5 minutes and the top plate was lowered to a defined gap to form a good contact with the material. The sandwiched amorphous material was heated to 80°C and the viscosity was measured. Results The equation was used to probe the correlation of viscosity to crystal growth rate. In comparison to the value of xi in log-log equation reported from pure drug compound, a value of 1.63 was obtained for FD-polymer solid dispersions irrespective of the polymer involved. &#8733 Conclusion The high xi value suggests stronger viscosity dependence may exist for amorphous FD once incorporated with amorphous polymer.