Theoretical Derivation of a Peak Particle Velocity–Distance Law for the Prediction of Vibrations from Blasting

I would like to suggest a theoretical justification for the mathematical structure of some laws for predicting the maximum particle velocity vibration from blasting operations in the light of some basic notions of elastic and anelastic wave theory. Within this point of view, in dimensionally correct expressions, the terms pertaining to the rock, to the blast and to the seismic wave become evident and recognisable. A law is presented that can be used to forecast the maximum particle velocity on the basis of some blast design and rock parameters. Four tests of the proposed law performed with real data sets seem to confirm fairly well its reliability.     

Limits to scale invariance in alluvial rivers

Assumptions about fluvial processes and process–form relations are made in general models and in many site-specific applications. Many standard assumptions about reach-scale flow resistance, bed-material entrainment thresholds and transport rates, and downstream hydraulic geometry involve one or other of two types of scale invariance: a parameter (e.g. critical Shields number) has the same value in all rivers, or doubling one variable causes a fixed proportional change in another variable in all circumstances (e.g. power-law hydraulic geometry). However, rivers vary greatly in size, gradient, and bed material, and many geomorphologists regard particular types of river as distinctive. This review examines the tension between universal scaling assumptions and perceived distinctions between different types of river. It identifies limits to scale invariance and departures from simple scaling, and illustrates them using large data sets spanning a wide range of conditions. Scaling considerations and data analysis support the commonly made distinction between coarse-bed and fine-bed reaches, whose different transport regimes can be traced to the different settling-velocity scalings for coarse and fine grains. They also help identify two end-member sub-types: steep shallow coarse-bed ‘torrents’ with distinctive flow-resistance scaling and increased entrainment threshold, and very large, low-gradient ‘mega rivers’ with predominantly suspended load, subdued secondary circulation, and extensive backwater conditions. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

A modified Raupach's model applicable for shear-stress partitioning on surfaces covered with dense and flat-shaped gravel roughness elements

A commonly used measure to prevent soil wind erosion is to cover the surface with gravel. Gravel can inhibit soil erosion by covering the surface directly, changing the airflow field near the surface and sharing the shear stress of wind. Similar to other roughness elements, the protective effect of gravel on soil is usually expressed in terms of the ratio of the shear stress on the exposed soil surface to the total shear stress on the rough surface due to wind, i.e. through a shear-stress partitioning model. However, the existing shear-stress partitioning models, represented by Raupach's model (RM93), are only applicable when the lateral coverage of the roughness elements, λ < 0.10, and the applicability of the models to flat-shaped roughness elements is unclear. The purpose of this study is to verify the applicability of RM93 for dense and flat-shaped gravel roughness elements by using shear-stress data from wind-tunnel measurements pertaining to roughness elements with different densities (0.013 ≤ λ ≤ 0.318) and flat shapes (height-to-width ratios in the range 0.20 ≤ H/W ≤ 0.63), and to modify RM93 to enhance its predictive ability. The results indicate that RM93 cannot accurately predict the shear-stress partitioning for surfaces covered by densely distributed and flat-shaped gravel roughness elements. This phenomenon occurs because, when roughness elements are distributed densely or are flat-shaped, the proportion of the shear stress on the top surface of the roughness elements (τ_c) to the total shear stress (τ) is large; in this case, τ_c plays a dominant role and serves as an essential component in the shear-stress partitioning model. Consequently, RM93 is modified by incorporating τ_c into the calculation of τ. Under conditions of λ < 0.32 and H/W > 0.2, the modified RM93 can yield satisfactory predictions regarding the shear-stress partitioning.     

Roughness effect on the correction factor of surface velocity for rill flows

Flow velocity is one of the most important hydrodynamic variables for both channelized (rill and gullies) and interrill erosive phenomena. The dye tracer technique to measure surface flow velocity V_s is based on the measurement of the travel time of a tracer needed to cover a known distance. The measured V_s must be corrected to obtain the mean flow velocity V using a factor α_v = V/V_s which is generally empirically deduced. The V_s measurement can be influenced by the method applied to time the travel of the dye-tracer and α_v can vary in different flow conditions. Experiments were performed by a fixed bed small flume simulating a rill channel for two roughness conditions (sieved soil, gravel). The comparison between a chronometer-based (CB) and video-based (VB) technique to measure V_s was carried out. For each slope-discharge combination, 20 measurements of V_s, characterized by a sample mean V_m, were carried out. For both techniques, the frequency distributions of V_s/V_m resulted independent of slope and discharge. For a given technique, all measurements resulted normally distributed, with a mean equal to one, and featured by a low variability. Therefore, V_m was considered representative of surface flow velocity. Regardless of roughness, the V_m values obtained by the two techniques were very close and characterized by a good measurement precision. The developed analysis on α_v highlighted that it is not correlated with Reynolds number for turbulent flow regime. Moreover, α_v is correlated neither with the Froude number nor with channel slope. However, the analysis of the empirical frequency distributions of the correction factor demonstrated a slope effect. For each technique (CB, VB)-roughness (soil, gravel) combination, a constant correction factor was statistically representative even if resulted in less accurate V estimations compared to those yielded by the slope-specific correction factor.     

依兰—伊通断裂带黑龙江段构造运动特征

为研究依兰—伊通断裂带黑龙江段构造运动特征,基于2016—2019年GPS和地质资料,解算了该断裂的三维速度场,通过构建断层模型反演了滑动速率。结果显示:依兰—伊通断裂带黑龙江段总体呈下沉趋势,沉降速率在1~2 mm/a,断裂呈右旋走滑态势,闭锁层15 km以下走滑速率为(1.7±0.4)mm/a。佳木斯—萝北段以右旋走滑为主、拉张为辅;五常—佳木斯段以拉张为主、右旋走滑为辅。     

某海底沉管隧道风塔及下部结构抗震性能分析

为研究某海底隧道风塔及下部结构体系在塑性阶段的损伤破坏形态及特点,探讨在多重荷载作用下弹塑性静动力响应对结构的影响,建立了沉管隧道风塔及下部结构的大型三维有限元模型。采用基于能量原理的混凝土塑性及损伤本构模型,借助大型有限元软件ANSYS及ABAQUS分别对结构进行不同地震条件下的动力时程分析,对比分析振型、层间位移角及较为完整的塑性损伤破坏系数曲线。结果表明:不同设防地震下,结构整体性良好,振型及层间位移角满足规范要求;不同罕遇地震下,该结构的混凝土塑性拉压损伤最大时刻均发生在20 s,主要破坏区域在风塔与人防井及下部立柱的接触位置,且拉伸破坏系数明显高于压缩破坏。本文的研究成果可以为类似近海沉管隧道工程抗震设计提供一定的依据与指导。     

The “slip law” of the free surface

A “slip law” connects the excess velocity or “slip” of a wind-blown water surface, relative to the motion in the middle of the mixed layer, to the wind stress, the wind-wave field, and buoyancy flux. An inner layer-outer layer model of the turbulent shear flow in the mixed layer is appropriate, as for a turbulent boundary layer or Ekman layer over a solid surface, allowing, however, for turbulent kinetic energy transfer from the air-side via breaking waves, and for Stokes drift. Asymptotic matching of the velocity distributions in inner and outer portions of the mixed layer yields a slip law of logarithmic form, akin to the drag law of a turbulent boundary layer. The dominant independent variable is the ratio of water-side roughness length to mixed layer depth or turbulent Ekman depth. Convection due to surface cooling is also an important influence, reducing surface slip. Water-side roughness length is a wind-wave property, varying with wind speed similarly to air-side roughness. Slip velocity is typically 20 times water-side friction velocity or 3% of wind speed, varying within a range of about 2 to 4.5%. A linearized model of turbulent kinetic energy distribution shows much higher values near the surface than in a wall layer. Nondimensional dissipation peaks at a value of about eight, a short distance below the surface.     

The Venusian induced magnetosphere: A case study of plasma and magnetic field measurements on the Venus Express mission

Plasma and magnetic field measurements made onboard the Venus Express on June 1, 2006, are analyzed and compared with predictions of a global model. It is shown that in the orbit studied, the plasma and magnetic field observations obtained near the North Pole under solar minimum conditions were qualitatively and, in many cases also, quantitatively in agreement with the general picture obtained using a global numerical quasi-neutral hybrid model of the solar wind interaction (HYB-Venus). In instances where the orbit of Venus Express crossed a boundary referred to as the magnetic pileup boundary (MPB), field line tracing supports the suggestion that the MPB separates the region that is magnetically connected to the fluctuating magnetosheath field from a region that is magnetically connected to the induced magnetotail lobes.     

The magnetic field near Mars: A comparison between a hybrid model, Mars Global Surveyor and Mars Express observations

Mars Express (MEX) does not carry its own magnetometer which complicates interpretation of ASPERA-3/MEX ion measurements. The direction of the interplanetary magnetic field (IMF) is especially important because it, among other things, determines the direction of the convective electric field and orientation of the cross tail current sheet and tail lobes. In this paper we present a case study to show the properties of the magnetic field near Mars in a quasi-neutral hybrid (QNH) model at the orbits where the Mars Global Surveyor (MGS) has made measurements, present a method to derive the IMF clock angle by comparing fields in a hybrid model and the direction of the magnetic field measured by MGS by deriving the IMF clock angle. We also use H⁺ ring velocity distribution observations upstream of the bow shock measured by the IMA/ASPERA-3 instrument on board MEX spacecraft. These observations are used to indirectly provide the orientation of the IMF. We use a QNH model (HYB-Mars) where ions are modeled as particles while electrons form a mass-less charge neutralizing fluid. We found that the direct MGS and non-direct IMA observations of the orientation magnetic field vectors in non-crustal magnetic field regions are consistent with the global magnetic field draping pattern predicted by the global model.     

Variations of the magnetic field near Mars caused by magnetic crustal anomalies

The possible avenues for photoelectron transport were determined during southern hemisphere winter at Mars by using a mapping analysis of the theoretical magnetic field. Magnetic field line tracing was performed by superposing two magnetic field models: (1) magnetic field derived from a three-dimensional (3D) self-consistent quasi-neutral hybrid model which does not contain the Martian crustal magnetic anomalies and (2) a 3D map of the magnetic field associated with the magnetic anomalies based on Mars Global Surveyor magnetic field measurements. It was found that magnetic field lines connected to the nightside of the planet are mainly channeled within the optical shadow of the magnetotail whereas magnetic field lines connected to the dayside of the planet are observed to form the remainder of the magnetosphere. The simulation suggests that the crustal anomalies create “a magnetic shield” by decreasing the region near Mars which is magnetically connected to the Martian magnetosphere. The rotation of Mars causes periodic changes in magnetic connectivity, but not to qualitative changes in the overall magnetic field draping around Mars.