Large-Eddy Simulation of Coherent Flow Structures within a Cubical Canopy

Instantaneous flow structures “within” a cubical canopy are investigated via large-eddy simulation. The main topics of interest are, (1) large-scale coherent flow structures within a cubical canopy, (2) how the structures are coupled with the turbulent organized structures (TOS) above them, and (3) the classification and quantification of representative instantaneous flow patterns within a street canyon in relation to the coherent structures. We use a large numerical domain (2,560 m × 2,560 m × 1,710 m) with a fine spatial resolution (2.5 m), thereby simulating a complete daytime atmospheric boundary layer (ABL), as well as explicitly resolving a regular array of cubes (40 m in height) at the surface. A typical urban ABL is numerically modelled. In this situation, the constant heat supply from roof and floor surfaces sustains a convective mixed layer as a whole, but strong wind shear near the canopy top maintains the surface layer nearly neutral. The results reveal large coherent structures in both the velocity and temperature fields “within” the canopy layer. These structures are much larger than the cubes, and their shapes and locations are shown to be closely related to the TOS above them. We classify the instantaneous flow patterns in a cavity, specifically focusing on two characteristic flow patterns: flushing and cavity-eddy events. Flushing indicates a strong upward motion, while a cavity eddy is characterized by a dominant vortical motion within a single cavity. Flushing is clearly correlated with the TOS above, occurring frequently beneath low-momentum streaks. The instantaneous momentum and heat transport within and above a cavity due to flushing and cavity-eddy events are also quantified.     

Experimental study on the collisional disruption of porous gypsum spheres

Abstract— In order to study the catastrophic disruption of porous bodies such as asteroids and planetesimals, we conducted several impact experiments using porous gypsum spheres (porosity: 50%). We investigated the fragment mass and velocity of disrupted gypsum spheres over a wide range of specific energies from 3 times 10³ J/kg to 5 times 10⁴ J/kg. We compared the largest fragment mass (m₁/M_t) and the antipodal velocity (V_a) of gypsum with those of non‐porous materials such as basalt and ice. The results showed that the impact strength of gypsum was notably higher than that of the non‐porous bodies; however, the fragment velocity of gypsum was slower than that of the non‐porous bodies. This was because the micro‐pores dispersed in the gypsum spheres caused a rapid attenuation of shock pressure in them. From these results, we expect that the collisional disruption of porous bodies could be significantly different from that of non‐porous bodies.     

Impacts of explosive volcanism on distal alluvial sedimentation: Examples from the Pliocene–Holocene volcaniclastic successions of Japan

Volcanic activities can create cataclysmic hazards to surrounding environments and human life not only during the eruption but also by hydrologic remobilisation (lahar) processes after the cessation of eruptive activity. Although there are many studies dealing with the assessment and mitigation of volcanic hazards, these are mostly concentrated on primary eruptive processes in areas proximal to active volcanoes. However, the influence of volcaniclastic resedimentation may surpass the impacts of primary eruptive activity in terms of both extent and persistence, and can ultimately result in severe hazards in downstream areas.Examination of the volcaniclastic successions of non-marine Pliocene–Holocene sedimentary basins in Japan has revealed hydrological volcaniclastic sedimentation in fluvial and lacustrine environments hundreds of kilometres from the inferred source volcano. Impacts on these distal and often spatially separated basins included drastic changes in depositional systems caused by sudden massive influxes of remobilised pyroclastic material. Typical volcaniclastic beds comprise centimetre- to decimetre-thick primary pyroclastic fall deposits overlain by metre- to 10s of metres-thick resedimented volcaniclastic deposits, intercalated in sedimentary successions of non-volcanic provenance. The relatively low component of primary pyroclastic fall deposits in the volcaniclastic beds suggests that: 1) potential volcanic hazards would be underestimated on the basis of primary pyroclastic fall events alone; and 2) the majority of resedimented material was likely derived from erosion of non-welded pyroclastic flow deposits in catchment areas rather than remobilisation of local fallout deposits from surrounding hillslopes.The nature, distribution and sequence of facies developed by distal volcaniclastic sediments reflect the influence of: 1) proximity to ignimbrite, but not directly with the distance to the eruptive centre; 2) ignimbrite nature (non-welded or welded) and volume; 3) temporal changes in sediment flux from the source area; 4) the physiography and drainage patterns of the source area and the receiving basin, and any intervening areas; and 5) the formation of ephemeral dam-lakes and intra-caldera lakes whose potential catastrophic failure can impact distal areas. Models of the styles and timing of distal volcaniclastic resedimentation are thus more complicated than those developed for proximal settings of stratovolcanoes and their volcaniclastic aprons and hence present different challenges for hazard assessment and mitigation.     

Organized Structure of Active Turbulence Over an Array of Cubes within the Logarithmic Layer of Atmospheric Flow

We investigate the coherent structure of atmosphere turbulence over very large roughness within a fully rough, high Reynolds number turbulent flow. The horizontal distributions of coherent turbulence were determined by multipoint measurements of velocity fluctuations using sonic anemometers in a comprehensive outdoor scale model experiment for urban climate (COSMO). COSMO is made up of 512 cubical obstacles, each 1.5 m on a side, arranged in a rectangular pattern on a flat 50 m × 100 m concrete plate. A total of 15 sets of sonic anemometers were aligned horizontally within the logarithmic layer above this site. The velocity fluctuations observed in COSMO were decomposed into active and inactive contributions by applying a spatial-filtering method, and which used a simple moving average along the spanwise direction of the predominant flow as a filter function. The size of the filter should be between the sizes of the active and inactive fluctuations. This method potentially eliminates the considerable portion of low frequency modes included in the horizontal velocity fluctuation, while preserving well the Reynolds stress. The structural characteristics of the active turbulence were qualitatively similar to those measured over various surface configurations. Overall, the observed structures of the active turbulence are composed of very large streaks of low momentum fluid elongated in the streamwise direction with some sub-structures included in the streaks. The sub-structures were the main cause of the ejections, which accompany horizontal vortices. The active motion, including the streaky structures, did not reproduce the lower frequency peak of the bi-modal distribution of the horizontal velocity spectra, but reproduced the higher frequency mode that robustly follows inner-layer similarity (i.e. Monin–Obukhov similarity).     

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.     

Theoretical Evaluation of the Mechanical Behavior of Unsaturated Soils

An evaluation method for the mechanical behavior of unsaturated soils is studied in this paper. Although the mechanical behavior of unsaturated soils is complicated, a simple modeling is preferable in practice. This is because the soil properties are not homogeneous and ground data is limited when structures are being designed. In addition, in order to evaluate the reliability of the design, the physical meanings of the parameters applied in the prediction model should be clear. Firstly, the authors study the relationship between compaction curves and compression indexes in the unsaturated state that is used in the proposed constitutive model. Based on the constitutive model, the stress paths for constant volume shear tests are formulated under a constant void ratio condition and the stress paths for undrained shear tests are calculated under a constant water content condition. In the case of unsaturated specimens, the volume of these specimens changes with the shear deformation and the stress paths depend on the initial degree of saturation. The results of the calculation qualitatively describe the test results by considering the changes in effective confining pressure in the undrained condition and the water retention curves.     

Thermal equation of state of CaIrO<Subscript>3</Subscript> post-perovskite

The pressure–volume–temperature (P–V–T) relation of CaIrO₃ post-perovskite (ppv) was measured at pressures and temperatures up to 8.6 GPa and 1,273 K, respectively, with energy-dispersive synchrotron X-ray diffraction using a DIA-type, cubic-anvil apparatus (SAM85). Unit-cell dimensions were derived from the Le Bail full profile refinement technique, and the results were fitted using the third-order Birth-Murnaghan equation of state. The derived bulk modulus \( K_{T0} \) at ambient pressure and temperature is 168.3 ± 7.1 GPa with a pressure derivative \( K_{T0}^{\prime } \) = 5.4 ± 0.7. All of the high temperature data, combined with previous experimental data, are fitted using the high-temperature Birch-Murnaghan equation of state, the thermal pressure approach, and the Mie-Grüneisen-Debye formalism. The refined thermoelastic parameters for CaIrO₃ ppv are: temperature derivative of bulk modulus \( (\partial K_{T} /\partial T)_{P} \) = ?0.038 ± 0.011 GPa K^?1, \( \alpha K_{T} \) = 0.0039 ± 0.0001 GPa K^?1, \( \left( {\partial K_{T} /\partial T} \right)_{V} \) = ?0.012 ± 0.002 GPa K^?1, and \( \left( {\partial^{2} P/\partial T^{2} } \right)_{V} \) = 1.9 ± 0.3 × 10^?6 GPa² K^?2. Using the Mie-Grüneisen-Debye formalism, we obtain Grüneisen parameter \( \gamma_{0} \) = 0.92 ± 0.01 and its volume dependence q = 3.4 ± 0.6. The systematic variation of bulk moduli for several oxide post-perovskites can be described approximately by the relationship K _T0  = 5406.0/V(molar) + 5.9 GPa.     

Measurements in a heated turbulent wall jet on a concave wall

The paper contains the results of extensive single-point hot-wire and resistance-thermometer measurements in a wall jet on a heated concave wall with an external free stream. It is found that the turbulence in the inner (wall) and outer layers is sensitive to the distortion produced by the curved wall, broadly confirming current views of these effects. The effect on the turbulence of streamline curvature is stabilising in the outer layer, destabilising in the inner. Consequently the point of zero shear stress is closer to the point of maximum mean velocity than in flat and convex wall jets with which these new results are compared. The rate of growth is about 80% of that of the equivalent flat wall flow, and about half that of the convex wall flow. Changes in the wall shear stress and heat flux, which are increased relative to the flat wall flow, are significant but less than the changes in the corresponding convex wall flow. The greatest changes occur in the triple products.List of Symbols c _f skin friction coefficient, 2 _w / U _M ² - q ² u ² +v ²+w ² - R radius of curvature of wall - r radius of curvature - St Stanton number,-q _w /c _p U _M (T –T _w ) - T mean temperature - U streamwise mean velocity - U _J initial wall-jet velocity - U ₁ initial free-stream velocity - U ₀ U _M –U _PW - u, v, w fluctuating velocity components - x, y distances measured along and normal to the wall - y _0.5,y _M wall-jet thicknesses defined in the text - y _S ,y _T wall-jet thicknesses defined in the text - T T _W –T - shear-layer thickness - ₂ momentum thickness - temperature fluctuation - M point of maximum velocity - P potential flow - s point of zero shear stress - W wall value - free stream value