Particle-Image Velocimetry Measurements of Separation and Re-attachment of Airflow over Two-Dimensional Hills with Various Slope Angles and Approach-Flow Characteristics

Airflow over two-dimensional hills was investigated in a wind tunnel using particle image velocimetry. We focus on the flow separation behaviour. A trapezoidal hill shape was used in most of the experimental runs, but the critical slope angle for flow separation was approximately the same as that established for smooth hill shapes. The re-attachment point of the separated flow became farther from the hill as the slope angle $\theta $ increased, reaching a saturation of about seven times the hill height for $\theta \gtrsim 60^\circ $ . Increasing the upwind surface roughness length was found to suppress flow separation. This tendency is analogous to the previous experimental results for turbulent boundary layers on flat plates. The boundary-layer thickness varied by the presence or absence of Counihan-type spires and a castellated fence at the test-section entrance had negligible effect on the flow separation.     

Large-Eddy Simulation of Turbulent Organized Structures within and above Explicitly Resolved Cube Arrays

Large-eddy simulations have been performed for fully developed turbulent flow within and above explicitly resolved simple cube arrays. The results from our model, hereafter LES-CITY, are shown to agree with laboratory experiments. We investigated the systematic influence of cube density on turbulent flow characteristics by performing numerical experiments for cube areal densities from 0 to 44%. The following results were obtained: (1) The dispersive momentum flux was quite large within the canopy layer due to a mean stream re-circulation, whereas it was smaller above the canopy. The spatial variation of temporally averaged momentum in the roughness sub-layer was 20% or less of the total kinematic surface drag. (2) The temporally and spatially-averaged flow structure confirmed the existence of conventionally described canyon flow regimes; isolated, interfacial, and wake. However, the intermittency of the canyon flow for all cube densities was quite large and the stream patterns were never persistent. (3) Turbulent organized structures (TOS) similar to those observed in turbulent surface-layer flows were simulated, which are characterized by longitudinally-elongated low speed streaks and the corresponding shorter streamwise vortices. The streaks in sparse and dense canopy flows were likely to be aligned to the street line and to the roof lines, respectively. Such heterogeneity of TOS partially accounts for the large spatial variation of momentum flux. (4) In contrast to the mixing layer analogy of vegetation flows, the TOS and the resulting turbulent statistics of urban flow above the canopy resembled those in surface layers. The recirculation within the canopy significantly influenced the turbulent statistical properties.     

A preparation zone for volcanic explosions beneath Naka-dake crater,Aso volcano,as inferred from magnetotelluric surveys

The 1st crater of Naka-dake, Aso volcano, is one of the most active craters in Japan, and known to have a characteristic cycle of activity that consists of the formation of a crater lake, drying-up of the lake water, and finally a Strombolian-type eruption. Recent observations indicate an increase in eruptive activity including a decrease in the level of the lake water, mud eruptions, and red hot glows on the crater wall. Temporal variations in the geomagnetic field observed around the craters of Naka-dake also indicate that thermal demagnetization of the subsurface rocks has been occurring in shallow subsurface areas around the 1st crater. Volcanic explosions act to release the energy transferred from magma or volcanic fluids. Measurement of the subsurface electrical resistivity is a promising method in investigating the shallow structure of the volcanic edifices, where energy from various sources accumulates, and in investigating the behaviors of magma and volcanic fluids. We carried out audio-frequency magnetotelluric surveys around the craters of Naka-dake in 2004 and 2005 to determine the detailed electrical structure down to a depth of around 1 km. The main objective of this study is to identify the specific subsurface structure that acts to store energy as a preparation zone for volcanic eruption. Two-dimensional inversions were applied to four profiles across the craters, revealing a strongly conductive zone at several hundred meters depth beneath the 1st crater and surrounding area. In contrast, we found no such remarkable conductor at shallow depths beneath the 4th crater, which has been inactive for 70 years, finding instead a relatively resistive body. The distribution of the rotational invariant of the magnetotelluric impedance tensor is consistent with the inversion results. This unusual shallow structure probably reflects the existence of a supply path of high-temperature volcanic gases to the crater bottom. We propose that the upper part of the conductor identified beneath the 1st crater is mainly composed of hydrothermally altered zone that acts both as a cap to upwelling fluids supplied from deep-level magma and as a floor to infiltrating fluid from the crater lake. The relatively resistive body found beneath the 4th crater represents consolidated magma. These results suggest that the shallow conductor beneath the active crater is closely related to a component of the mechanism that controls volcanic activity within Naka-dake.     

An overview of the geology and major ore deposits of Central Africa: Explanatory note for the 1:4,000,000 map “Geology and major ore deposits of Central Africa”

This paper is prepared within the frameworks of IGCP Project 470 and the associated BRGM scientific project “Africa 1999–2004” to accompany the 1:4,000,000 scale map “Geology and major ore deposits of Central Africa, presented at the 20th Colloquium of African Geology in Orleans in June 2004. It incorporates geological and metallogenic data from eight countries in Central Africa (Angola, Cameroon, Chad, Central African Republic, Congo Brazzaville, Democratic Republic of Congo (DRC), Equatorial Guinea and Zambia). The map is a harmonised and geo-referenced preliminary map, based on a GIS at 1:2,000,000 scale, and focusses on the spatial and temporal distribution of selected major deposits.     

Liquid waveguide spectrophotometric measurements of arsenate and particulate arsenic,as well as phosphate and particulate phosphorus,in seawater

Sensitive methods for the determination of arsenate and particulate arsenic (PAs), as well as phosphate and particulate phosphorus (PP), in seawater are described. The method for arsenate and phosphate was established by applying automated liquid waveguide spectrophotometry. Because the reaction time for the formation of the arsenate-molybdate complex is longer than that for the phosphate-molybdate complex, a long Teflon tube submerged in a heating bath was installed in the conventional phosphate flow system. The arsenate was quantified as the difference between absorbances of molybdenum blue dyes with (only phosphate) and without (phosphate + arsenate) arsenate reduction treatment. Contamination was observed in the reagent for arsenate reduction and must be corrected. Linear dynamic ranges up to 1000 nM were confirmed for arsenate and phosphate. The detection limits for arsenate and phosphate were 5 and 4 nM, respectively. Freezing was a reliable sample preservation technique for both arsenate and phosphate. The method for PAs and PP was established by combining conventional persulfate oxidation of PP and the automated liquid waveguide spectrophotometry of arsenate and phosphate. The digestion efficiencies of organic As analogs were >93%. Contamination in the glass fiber filter was negligible. Field tests confirmed that the coefficients of variation (CVs) of 10–19 nM arsenate and 4–151 nM phosphate were 7–20 and 1–25%, respectively, while the CVs of 0.9 nM PAs and 10.2 nM PP were 11 and 4%, respectively.     

Thermal Image Velocimetry

A method for measuring the two-dimensional distribution of wind velocity vectors near a surface exposed to solar radiation, by tracking brightness temperature images instead of particle images, is proposed. It is based on time-sequential thermography with the algorithm used for particle image velocimetry. This thermal image velocimetry (TIV) was tested on a full-scale building wall covered by polystyrene boards attached side-by-side over a vertically elongated area measuring 22.2 m by 2.73 m. A thermal infrared camera was installed 8 m from the test wall to capture the wall-surface temperature at 30 Hz frequency. A sonic anemometer was also installed 35 mm from the surface used for validation of the TIV. The advection velocity estimated from thermal infrared imagery had a linear relationship with the wind velocity measured by the sonic anemometer, irrespective of the wind speed and direction. This linear slope was multiplied by the advection velocity of the thermal infrared image to rescale it to the wind velocity, and the term ‘TIV velocity’ was then used. A histogram and power spectra of the TIV velocity showed quantitatively good agreement with the velocity measured by the sonic anemometer, except for the high-frequency region of the spectra, where the TIV velocity was overestimated compared with that of the sonic anemometer. The method was also tested on ground covered by artificial turf to demonstrate its application to a horizontal plane with a wider area, extending for more than 80 m by 60 m.