Behaviour of a puff of resuspended sediment: a conceptual model

The particulate matter concentration above the seabed is usually assumed to decrease with height, following an exponential or Rouse profile. Many particulate matter concentration profiles with a peak were found on the North Mediterranean bottom water at a few tens of metres above the bottom. A particle size signal at the same altitude was found in this area and on the New York Eight shelf. It is assumed that this unexpected shape is due to a cloud of resuspended cohesive sediments originating from an impulse resuspension process. A simplified three-dimensional numerical model is proposed to describe the behaviour of resuspended particulate matter that originates from a sediment impulse vertically injected in the bottom water. This model reproduces the concentration profile shape observed, and it gives indications concerning the length and time characteristics of such a cloud, depending on the water velocity and bottom boundary layer properties.

Experimental Investigation Of Particle Velocity Distributions In Windblown Sand Movement

With the PDPA (Phase Doppler Particle Analyzer) measurement technology, the probability distributions of particle impact and lift-off velocities on bed surface and the particle velocity distributions at different heights are detected in a wind tunnel. The results show that the probability distribution of impact and lift-off velocities of sand grains can be expressed by a log-normal function, and that of impact and lift-off angles complies with an exponential function. The mean impact angle is between 28 degrees and 39 degrees, and the mean lift-off angle ranges from 30 degrees to 44 degrees. The mean lift-off velocity is 0.81-0.9 times the mean impact velocity. The proportion of backward-impacting particles is 0.05-0.11, and that of backward-entrained particles ranges from 0.04 to 0.13. The probability distribution of particle horizontal velocity at 4 mm height is positive skew, the horizontal velocity of particles at 20 mm height varies widely, and the variation of the particle horizontal velocity at 80 mm height is less than that at 20 mm height. The probability distribution of particle vertical velocity at different heights can be described as a normal function.

Measurement of velocity profiles in a rectangular microchannel with aspect ratio α = 0.35

In this work, we measured 14 horizontal velocity profiles along the vertical direction of a rectangular microchannel with aspect ratio alpha = h/w = 0.35 (h is the height of the channel and w is the width of the channel) using microPIV at Re = 1.8 and 3.6. The experimental velocity profiles are compared with the full 3D theoretical solution, and also with a Poiseuille parabolic profile. It is shown that the experimental velocity profiles in the horizontal and vertical planes are in agreement with the theoretical profiles, except for the planes close to the wall. The discrepancies between the experimental data and 3D theoretical results in the center vertical plane are less than 3.6%. But the deviations between experimental data and Poiseuille's results approaches 5%. It indicates that 2D Poiseuille profile is no longer a perfect theoretical approximation since a = 0.35. The experiments also reveal that, very near the hydrophilic wall (z = 0.5-1 mu m), the measured velocities are significantly larger than the theoretical velocity based on the no-slip assumption. A proper discussion on some physical effects influencing the near wall velocity measurement is given.

Complex Transition of Double-Diffusive Convection in a Rectangular Enclosure with Height-to-Length Ratio Equal to 4: Part I

This is the first part of direct numerical simulation (DNS) of double-diffusive convection in a slim rectangular enclosure with horizontal temperature and concentration gradients. We consider the case with the thermal Rayleigh number of 10^5, the Pradtle number of 1, the Lewis number of 2, the buoyancy ratio of composition to temperature being in the range of [0,1], and height-to-width aspect ration of 4. A new 7th order upwind compact scheme was developed for approximation of convective terms, and a three-stage third-order Runge-Kutta method was employed for time advancement. Our DNS suggests that with the buoyancy ratio increasing form 0 to 1, the flow of transition is a complex series changing fromthe steady to periodic, chaotic, periodic, quasi-periodic, and finally back to periodic. There are two types of periodic flow, one is simple periodic flow with single fundamental frequency (FF), and another is complex periodic flow with multiple FFs. This process is illustrated by using time-velocity histories, Fourier frequency spectrum analysis and the phase-space rajectories.

Computational study of the vortex path variation with the VG height

An extensive range of conventional, vane-type, passive vortex generators (VGs) are in use for successful applications of flow separation control. In most cases, the VG height is designed with the same thickness as the local boundary layer at the VG position. However, in some applications, these conventional VGs may produce excess residual drag. The so-called low-profile VGs can reduce the parasitic drag associated to this kind of passive control devices. As suggested by many authors, low-profile VGs can provide enough momentum transfer over a region several times their own height for effective flow-separation control with much lower drag. The main objective of this work is to study the variation of the path and the development of the primary vortex generated by a rectangular VG mounted on a flat plate with five different device heights h = delta, h(1) = 0.8 delta, h(2) = 0.6 delta, h(3) = 0.4 delta and h(4) = 0.2 delta, where delta is the local boundary layer thickness. For this purpose, computational simulations have been carried out at Reynolds number Re = 1350 based on the height of the conventional VG h = 0.25m with the angle of attack of the vane to the oncoming flow beta = 18.5 degrees. The results show that the VG scaling significantly affects the vortex trajectory and the peak vorticity generated by the primary vortex.

Estudo dos componentes metabólicos, hormonais, marcadores inflamatórios e função microvascular em crianças com excesso de peso

Estudos da microcirculação cutânea demonstram que a disfunção microvascular neste sítio está relacionada a diversos fatores de risco cardiovascular. Existem poucos estudos avaliando a reatividade microvascular em crianças e a interferência da puberdade está presente na maioria deles. O objetivo deste estudo foi avaliar se a disfunção microvascular está presente em crianças pré-púberes com excesso de peso através da técnica de videocapilaroscopia de leito periungueal. Realizou-se um estudo transversal com 52 obesos, 18 sobrepesos e 28 eutróficos, com idade de 7,44 1,22 anos. Avaliou-se o comportamento dos fatores de risco e a função microvascular. A reatividade microvascular foi testada através da avaliação da densidade capilar funcional, da velocidade de deslocamento das hemácias em repouso e após uma isquemia de 1 min, e do tempo de reperfusão durante a hiperemia reativa. Análise de função disciminante canônica foi utilizada de forma multivariada para testar a possibilidade de separação dos grupos conforme o grau de adiposidade. Nos pacientes estudados não observamos diferença na reatividade microvascular, em nenhuma da variáveis testadas. Conforme esperado, os grupos obeso e sobrepeso apresentavam maiores valores para a circunferência da cintura (p<0,001), a relação cintura/altura (p<0,001), a pressão arterial média (p<0,001), o homeostasis model assessment for insulin resistance (HOMA-IR) (p<0,001) e os níveis de insulina (p<0,001), leptina (p<0,0001), glicose (p=0,02), triglicerídeos (p<0,05), colesterol total (p=0,004), ácido úrico (p=0,007) e proteína C reativa (p<0,0001) do que os eutróficos. A análise multivariada demonstrou a associação de variáveis metabólicas, antropométricas e microvasculares, sendo que estas foram separadas pelo grau de adiposidade corporal. Concluímos que nessa população estudada, apesar das diferenças nos perfis metabólico, inflamatório e hormonal, não houve diferença na reatividade microvascular. Entretanto, a associação entre variáveis clínico-antropométricas com aquelas relacionadas com a reatividade microvascular esteve presente nestas crianças pré-púberes e o grau de adiposidade corporal foi capaz de influenciar estas associações.

Velocity distributions for low pressure turbines

A parametric set of velocity distributions has been investigated using a flat plate experiment. Three different diffusion factors and peak velocity locations were tested. These were designed to mimic the suction surfaces of Low Pressure (LP) turbine blades. Unsteady wakes, inherent in real turbomachinery flows, were generated using a moving bar mechanism. A turbulence grid generated a freestream turbulence level that is believed to be typical of LP turbines. Measurements were taken across a Reynolds number range of 50,000-220,000 at three reduced frequencies (0.314, 0.628, 0.942). Boundary layer traverses were performed at the nominal trailing edge using a Laser Doppler Anemometry system and hot-films were used to examine the boundary layer behaviour along the surface. For every velocity distribution tested, the boundary layer separated in the diffusing flow downstream of the peak velocity. The loss production is dominated by the mixing in the reattachment process, mixing in the turbulent boundary layer downstream of reattachment and the effects of the unsteady interaction between the wakes and the boundary layer. A sensitive balance governs the optimal location of peak velocity on the surface. Moving the velocity peak forwards on the blade was found to be increasingly beneficial when bubblegenerated losses are high, i.e. at low Reynolds number, at low reduced frequency and at high levels of diffusion. Copyright © 2008 by ASME.

The rhythm of fountains: The length and time scales of rise height fluctuations at low and high Froude numbers

The magnitude and frequency of vertical fluctuations of the top of an axisymmetric miscible Boussinesq fountain forms the focus of this work. We present measurements of these quantities for saline-aqueous fountains in uniform quiescent surroundings. Our results span source Froude numbers 0.3 ≤ Fr 0 ≤ 40 and, thereby, encompass very weak, weak, intermediate and forced classes of fountain. We identify distinct scalings, based on known quantities at the fountain source, for the frequency of fountain height fluctuations which collapse our data within bands of Fr0. Notably, our scalings reveal that the (dimensionless) frequency takes a constant value within each band. These results highlight characteristic time scales for the fluctuations which we decompose into a single, physically apparent, length scale and velocity scale within each band. Moreover, within one particular band, spanning source Froude numbers towards the lower end of the full range considered, we identify unexpectedly long-period fluctuations indicating a near balance of inertia and (opposing) buoyancy at the source. Our analysis identifies four distinct classes of fluctuation behaviour (four bands of Fr 0) and this classification matches well with existing classifications of fountains based on rise heights. As such, we show that an analysis of the behaviour of the fountain top alone, rather than the entire fountain, provides an alternative approach to classifying fountains. The similarity of classifications based on the two different methods confirms that the boundaries between classes mark tangible changes in the physics of fountains. For high Fr0 we show that the dominant fluctuations occur at the scale of the largest eddies which can be contained within the fountain near its top. Extending this, we develop a Strouhal number, Strtop, based on experimental measures of the fountain top, defined such that Strtop = 1 would suggest the dominant fluctuations are caused by a continual cycle of eddies forming and collapsing at this largest physical scale. For high- Fr 0 fountains we find Strtop ≈ 0. 9. © 2013 Cambridge University Press.

Dynamic compression of foam supported plates impacted by high velocity soil

The response of back-supported buffer plates comprising a solid face sheet and foam core backing impacted by a column of high velocity particles (sand slug) is investigated via a lumped parameter model and coupled discrete/continuum simulations. The buffer plate is either resting on (unattached) or attached to a rigid stationary foundation. The lumped parameter model is used to construct maps of the regimes of behaviour with axes of the ratio of the height of the sand slug to core thickness and the normalised core strength. Four regimes of behaviour are identified based on whether the core compression ends prior to the densification of the sand slug or vice versa. Coupled discrete/continuum simulations are also reported and compared with the lumped parameter model. While the model predicted regimes of behaviour are in excellent agreement with numerical simulations, the lumped parameter model is unable to predict the momentum transmitted to the supports as it neglects the role of elasticity in both the buffer plate and the sand slug. The numerical calculations show that the momentum transfer is minimised for intermediate values of the core strength when the so-called "soft-catch" mechanism is in play. In this regime the bounce-back of the sand slug is minimised which reduces the momentum transfer. However, in this regime, the impulse reduction is small (less than 10% of that transferred to a rigid structure). For high values of the core strength, the response of the buffer plate resembles a rigid plate with nearly no impulse mitigation while at low values of core strength, a slap event occurs when the face sheet impinges against the foundation due to full densification of the foam core. This slap event results in a significant enhancement of the momentum transfer to the foundation. The results demonstrate that appropriately designed buffer plates have potential as impulse mitigators in landmine loading situations. © 2013 Elsevier Ltd. All rights reserved.

Comparison of acoustic velocity perturbation measurements using PIV vs. two-microphone technique

Understanding combustion instabilities requires accurate measurements of the acoustic velocity perturbation into injectors. This is often accomplished via the use of the two microphone technique, as this only requires two pressure transducers. However, measurements of the actual velocities emerging from the injectors are not often taken, leaving questions regarding the assumptions about the acoustic velocity. A comparison of velocity measured at downstream of the injector with that of two-microphone technique can show the accuracy and limitations of two-microphone technique. In this paper, velocity measurements are taken using both particle image velocimetry (PIV) and the two-microphone technique in a high pressure facility designed for aeroengine injector measurements. The flow is excited using an area modulation device installed on the choked end of the combustion chamber, with PIV measurements enabled by optical access downstream of the injector through a quartz tube and windows. Acoustic velocity perturbations at the injector are determined by considering the Fourier transformed pressure fluctuations for two microphones installed at a known distance upstream of the injector. PIV measurements are realized by seeding the air flow with micrometric water particles under 2.5 bar pressure at ambient temperature. Phase locked velocity fields are realized by synchronizing the acquisition of PIV images with the revolution of the acoustic modulator using the pressure signal measured at the face of injector. The mean velocity fluctuation is calculated as the difference between maximum and minimum velocities, normalized by the mean velocity of the unforced case. Those values are compared with the peak-to-peak velocity fluctuation amplitude calculated by the two-microphone technique. Although the ranges of velocity fluctuations for both techniques are similar, the variation of fluctuation with forcing frequencies diverges significantly with frequency. The differences can be attributed to several limitations associated with of both techniques, such as the quality of the signal, the signal/noise ratio, the accuracy of PIV measurements and the assumption of isentropic flow of the particle velocity from the plenum through the injector. We conclude that two-microphone methods can be used as a reference value for the velocity fluctuation in low order applications such as flame transfer functions, but not for drawing conclusions regarding the absolute velocity fluctuations in the injector. Copyright © 2013 by ASME.

Response of clamped sandwich beams subjected to high-velocity impact by sand slugs

The dynamic response of end-clamped sandwich and monolithic beams of equal areal mass subjected to loading via high-velocity slugs of dry and water-saturated sand is measured using a novel laboratory-based method. The sandwich beams comprise aluminium face sheets and an aluminium honeycomb core: the effect of sandwich core strength and beam thickness on the dynamic beam deflection is investigated by varying the orientation and height of the anisotropic aluminium honeycomb core material. High-speed imaging is used to measure the transient transverse deflection of the beams and to record the dynamic modes of deformation. The measurements show that sandwich beams with thick, strong cores are optimal and that these beams significantly outperform monolithic beams of equal mass. The water-saturated sand slugs cause significantly higher deflections compared to the dry sand slugs having the same mean slug velocity and we demonstrate that this enhanced deflection is due to the larger mass of the water-saturated slugs. Finally, we show that the impact of sand slugs is equivalent to the impact of a crushable foam projectile. The experiments using foam projectiles are significantly simpler to perform and thus represent a more convenient laboratory technique. © 2014 Elsevier Ltd. All rights reserved.

边坡动力响应分析及应用研究

As a marginal subject, dynamic responses of slopes is not only an important problem of engineering geology (Geotechnical problem), but also of other subjects such as seismology, geophysics, seismic engineering and engineering seismic and so on. Owning to the gulf between different subjects, it is arduous to study dynamic responses of slopes and the study is far from ripeness. Studying on the dynamic responses of slopes is very important in theories as well as practices. Supported by hundreds of bibliographies, this paper systemically details the development process of this subject, introduces main means to analyze this subject, and then gives brief remarks to each means respectively. Engineering geology qualitative analysis is the base of slopes dynamic responses study. Because of complexity of geological conditions, engineering geology qualitative analysis is very important in slopes stability study, especially to rock slopes with complex engineering geology conditions. Based on research fruits of forerunners, this paper summarizes factors influencing slopes dynamic stability into five aspects as geology background, stratums, rock mass structure, and topography as well as hydrogeology condition. Based on rock mass structure controlling theory, engineering geology model of the slope is grouped into two typical classes, one is model with obvious controlling discontinuities, which includes horizontal bedded slope, bedding slope, anti-dip slope, slide as well as slope with base rock and weathered crust; the other is model without obvious controlling discontinuities, which includes homogeneous soil slope, joint rock mass slope. Study on slope failure mechanism under dynamic force, the paper concludes that there are two effects will appear in slope during strong earthquake, one is earthquake inertia force, the other is ultra pore pressure buildup. The two effects lead to failure of the slope. To different types of slope failure, the intensity of two effects acting on the slope is different too. To plastic flow failure, pore pressure buildup is dominant; to falling rock failure and toppling failure, earthquake inertia force is dominant in general. This paper briefly introduces the principle of Lagrangian element method. Through a lot of numerical simulations with FLAC3D, the paper comprehensively studies dynamic responses of slopes, and finds that: if the slope is low, displacement, velocity and acceleration are linear enlarging with elevation increasing in vertical direction; if the slope is high enough, displacement, velocity and acceleration are not linear with elevation any more, on the other hand, they fluctuate with certain rhythm. At the same time, the rhythm appears in the horizontal direction in the certain area near surface of the slope. The distribution form of isoline of displacement, velocity and acceleration in the section of the slope is remarkably affected by the slope angle. In the certain area near the slope surface, isoline of displacement, velocity and acceleration is parallel to the surface of the slope, in the mean time, the strike direction of the extreraum area is parallel to the surface of the slope too. Beyond this area, the isoline direction and the strike direction of the extremum area turn to horizontal with invariable distance. But the rhythm appearing or not has nothing to with the slope angle. The paper defines the high slope effect and the low slope effect of slopes dynamic responses, discusses the threshold height H^t of the dynamic high slope effect, and finds that AW is proportional to square root of the dynamic elastic moduli El P , at the same time, it is proportional to period Tof the dynamic input. Thus, the discriminant of H^t is achieved. The discriminant can tell us that to a slope, if its height is larger than one fifth of the wavelength, its response regular will be the dynamic high slope effect; on the other hand, its response regular will be the dynamic low slope effect. Based on these, the discriminant of different slopes taking on same response under the same dynamic input is put forward in this paper. At the same time, the paper studies distribution law of the rhythm extremum point of displacement, velocity and acceleration, and finds that there exists relationship of N = int among the slope height H, the number of the rhythm extremum VHlhro) point N and ffthre- Furthermore, the paper points out that if N^l, the response of the slope will be dynamic high slope effect; \fNpeak ground acceleration) of the site is decided (about 197.1cm/s2) by engineering seismology, and then the seismic input used in the slope dynamic analysis is determined. Comprehensive studies are carried out on the slope, especially to its deep fractures, and then the paper concludes that the deep fractures of slope are the result of the combination unloading effect of gravity and tectonic stress. At the same time, the failure model of the slope under dynamic input is attained. Based on these, the stability is comprehensively studied for section IV-IV with numerical simulation method as well as method putting forward in chapter 5. At last, the paper concludes that under dynamic input, the section will slide along fault f9, some deep fractures and fault/5 with certain permanent displacement, and this must be taken into consideration in the engineering.

CaII K interstellar observations towards early-type disc and halo stars - distances to intermediate- and high-velocity clouds

We compare existing high spectral resolution (R = lambda/Deltalambda similar to 40 000) Ca II Kobservations (lambda(air) = 3933.66 Angstrom) towards 88 mainly B-type stars, and new observations taken using the Intermediate dispersion Spectrograph and Imaging System (ISIS) on the William Herschel Telescope at R similar to 10 000 towards three stars taken from the Palomar-Green Survey, with 21-cm HI emission-line profiles, in order to search for optical absorption towards known intermediate- and high-velocity cloud complexes. Given certain assumptions, limits to the gas phase abundance of Ca II are estimated for the cloud components. We use the data to derive the following distances from the Galactic plane (z). (i) Tentative lower z-height limits of 2800 and 4100 pc towards complex C using lack of absorption in the spectra of HD341617 and PG 0855 + 294, respectively. (ii) A weak lower z-height of 1400 pc towards complex WA-WB using lack of absorption in EC 09470-1433 and a weak lower limit of 2470 pc using lack of absorption in EC 09452-1403. (iii) An upper z- height of 2470 pc towards a southern intermediate- velocity cloud (IVC) with v(LSR) = -55 km s(-1) using PG 2351 + 198. (iv) Detection of a possible IVC in Ca II absorption at v(LSR) = +52 km s(-1) using EC 20104-2944. No associated HI in emission is detected. At this position, normal Galactic rotation predicts velocities of up to similar to+ 25 km s(-1). The detection puts an upper z-height of 1860 pc to the cloud. (v) Tentative HI and Ca II K detections towards an IVC at similar to+70 km s(-1) in the direction of high-velocity cloud (HVC) complex WE, sightline EC 06387-8045, indicating that the IVC may be at a z-height lower than 1770 pc. (vi) Detection of Ca II K absorption in the spectrum of PG 0855 + 294 in the direction of IV20, indicating that this IVC has a z-height smaller than 4100 pc. (vii) A weak lower z-height of 4300 pc towards a small HVC with v(LSR) = +115 km s(-1) at l, b = 200degrees, + 52degrees, using lack of absorption in the Ca II K spectrum of PG 0955 + 291.

Multiwavelength observations of the M15 intermediate-velocity cloud

We present Westerbork Synthesis Radio Telescope HI images, Lovell telescope multibeam H I wide-field mapping, William Herschel Telescope long-slit echelle Ca II observations, Wisconsin Halpha Mapper (WHAM) facility images, and IRAS ISSA 60- and 100-mum co-added images towards the intermediate- velocity cloud (IVC) at + 70 km s(-1), located in the general direction of the M15 globular cluster. When combined with previously published Arecibo data, the H I gas in the IVC is found to be clumpy, with a peak H I column density of similar to1.5 x 10(20) cm(-2), inferred volume density (assuming spherical symmetry) of similar to24 cm(-3)/D (kpc) and a maximum brightness temperature at a resolution of 81 x 14 arcsec(2) of 14 K. The major axis of this part of the IVC lies approximately parallel to the Galactic plane, as does the low- velocity H I gas and IRAS emission. The H I gas in the cloud is warm, with a minimum value of the full width at half-maximum velocity width of 5 km s(-1) corresponding to a kinetic temperature, in the absence of turbulence, of similar to540 K. From the H I data, there are indications of two-component velocity structure. Similarly, the Ca II spectra, of resolution 7 km s(-1), also show tentative evidence of velocity structure, perhaps indicative of cloudlets. Assuming that there are no unresolved narrow-velocity components, the mean values of log(10)[N(Ca II K) cm(2)] similar to 12.0 and Ca II/H I similar to2 5 x 10(-8) are typical of observations of high Galactic latitude clouds. This compares with a value of Ca II/H I>10(-6) for IVC absorption towards HD 203664, a halo star of distance 3 kpc, some 3.degrees1 from the main M15 IVC condensation. The main IVC condensation is detected by WHAM in Halpha with central local-standard-of-rest velocities of similar to60-70 km s(-1), and intensities uncorrected for Galactic extinction of up to 1.3 R, indicating that the gas is partially ionized. The FWHM values of the Halpha IVC component, at a resolution of 1degrees, exceed 30 km s(-1). This is some 10 km s(-1) larger than the corresponding H I value at a similar resolution, and indicates that the two components may not be mixed. However, the spatial and velocity coincidence of the Halpha and H I peaks in emission towards the main IVC component is qualitatively good. If the Halpha emission is caused solely by photoionization, the Lyman continuum flux towards the main IVC condensation is similar to2.7 x 10(6) photon cm(-2) s(-1). There is not a corresponding IVC Halpha detection towards the halo star HD 203664 at velocities exceeding similar to60 km s(- 1). Finally, both the 60- and 100-mum IRAS images show spatial coincidence, over a 0.675 x 0 625 deg(2) field, with both low- and intermediate-velocity H I gas (previously observed with the Arecibo telescope), indicating that the IVC may contain dust. Both the Halpha and tentative IRAS detections discriminate this IVC from high-velocity clouds, although the H I properties do not. When combined with the H I and optical results, these data point to a Galactic origin for at least parts of this IVC.

Optical spectroscopy of the M 15 intermediate velocity cloud

We present echelle spectrograph observations in the Na D lines, at resolutions of 6.2-8.5 km s(-1), for 11 stars located in the line-of-sight to the M15 intermediate velocity cloud (IVC), which has a radial velocity of similar to +70 km s(-1) in the Local Standard of Rest. This cloud is a part of IVC Complex gp. The targets range in magnitude from m(V) = 13.3-14.8. Seven of the observed stars are in the M15 globular cluster, the remaining four being field stars. Three of the observed cluster stars are located near a peak in intensity of the IVC Hi column density as observed at a resolution of similar to 1 arcmin. Intermediate velocity gas is detected in absorption towards 7 stars, with equivalent widths in NaD2 ranging from similar to0.09-0.20 Angstrom, corresponding to log(10)(N-Na cm(-2)) similar to 11.8-12.5, and Na I/H I column density ratios (neglecting the HII component) ranging from similar to(1-3) x 10(-8). Over scales ranging from 30 arcsec to 1 arcmin, the Na i column density and the Na i/H i ratio varies by upto 70 per cent and a factor of similar to 2, respectively. Combining the current sightlines with previously obtained Nai data from Kennedy et al. (1998b), the Na i/H i column density ratio over cluster sightlines varies by upto a factor of similar to 25, when using Hi data of resolution similar to 2 x 1 arcmin. One cluster star, M15 ZNG-1, was also observed in the Ca i (lambda(air) = 4226.728 Angstrom) and Ca ii (lambda(air) = 3933.663 Angstrom) lines. A column density ratio N(Ca i)/N(Ca ii) <0.03 was found, typical of values seen in the warm ionised interstellar medium. Towards this sightline, the IVC has a Nai/Ca ii column density ratio of &SIM; 0.25, similar to that observed in the local interstellar medium. Finally, we detect tentative evidence for IV absorption in Ki (?(air) = 7698:974 &ANGS) towards 3 cluster stars, which have N(K i)/N(H i) ratios of &SIM;0.5-3 x 10(-9).