Nonplanar multi-modal vibrations of fluid-conveying risers under shear cross flows

Risers/pipes conveying fluid are a typical kind of slender structures commonly used in marine engineering. It is of great academic significance and application value for us to evaluate and understand the vibration characteristics and nonlinear responses of these risers under the combined action of internal and external fluid flows. In this paper, the nonplanar vibrations and multi-modal responses of pinned-pinned risers in shear cross flow are numerically studied. With this objective in mind, the van der Pol wake oscillators are used to simulate the dynamical behavior of the vortex shedding in the wake. Two nonlinear equations of motion of the riser are proposed to govern the lateral responses of the riser structure. The nonplanar nonlinear equations for the riser and wake are then discretized by employing Galerkin's method and solved by using a fourth-order Runge–Kutta integration algorithm. Theoretical results show that the coupled frequencies for cross-flow (CF) and in-line (IL) motions and the corresponding coupled damping ratio could be influenced by the external and/or internal fluid velocities. Based on extensive calculations, the dynamical behavior of the riser with various internal and external flow velocities are presented in the form of bifurcation diagrams, time traces, phase portraits, oscillation trajectories and response spectrum curves. It is shown that some interesting dynamical phenomena, such as ‘lock-in’ state, ‘figure-of-eight’ trajectory and quasi-periodic oscillation, could occur in such a fluid-structure interaction system. Our results also demonstrate that the shear parameter can significantly affect the dynamic responses of the riser. When the shear parameter of the cross flow is large, multi-modal quasi-periodic responses of the riser can be excited, showing some new features undetected in the system of fluid-conveying risers in uniform cross flow.     

Preliminary seismic anisotropy in the upper crust of the south segment of Xiaojiang faults and its tectonic implications

The Xiaojiang faults,striking north-to-south(NS),and the Honghe faults,striking north-to-west(NW),are first-order block boundaries that intersect to form a concentrated stress zone at an acute angle in the southern part of the Sichuan-Yunnan rhombic block(SYB).It is also a crucial zone for material escaping from the Tibetan Plateau(TP)due to the collision between the Indian Plate and the Eurasian Plate.In December 2017,the Institute of Earthquake Forecasting of the China Earthquake Administration(CEA)deployed a linear temporary seismic broadband array,the Honghe-Xiaojiang temporary Seismic Array(HX Array),across first-order block boundaries in the southern SYB.By using the waveform data of small earthquakes recorded by stations in the HX Array across Xiaojiang faults from 2017 to 2019,and by permanent seismic stations of the China National Earthquake Networks from 2012 to 2019,this paper adopts the systematic analysis method of shear-wave splitting(SWS),SAM method,to obtain preliminary results for seismic anisotropy in the upper crust.The study area can be divided into two subzones according to the spatial distribution of the directions of polarization of the fast shear-wave(PFS)at the stations:the northern zone(zone A,where the HX Array is located)and the southern zone(zone B,to the south of the HX Array).The results show that the directions of the PFS at stations in zone A were highly consistent,dominant in the NE direction,correlated with the in-situ principal compressive stress,and were seemingly unaffected by the Xiaojiang faults.The directions of the PFS as recorded at stations in zone B were more complicated,and were dominant in the NS direction parallel to that of the regional principal compressive stress.This suggests the joint influence of complex tectonics and regional stress in this narrow wedge area.By referring to the azimuthal anisotropy derived from seismic ambient noise in the southeast margin of the TP,the NS direction of the PFS in the middle and lower crust,and its EW direction in the upper mantle,this paper concludes that azimuthal anisotropy in the upper crust differed from that in the lower crust in the south segment of Xiaojiang faults,at least beneath the observation area,and azimuthal anisotropy in the crust was different from that in the upper mantle.The results support the pattern of deformation of ductile flow in the lower crust,and the decoupling between the upper and lower crusts as well as that between the crust and the mantle in the study area.The crustal directions of the PFS appeared to be independent of the Xiaojiang faults,suggesting that the influence of the South China block on the SYB passed through the Xiaojiang faults to the Yimen region.The results of this study indicate that anisotropic studies based on data on the dense temporary seismic array can yield clearer tectonic information,and reveal the complex spatial distribution of stress and deformation in the upper crust of the south segment of Xiaojiang faults.     

Genesis of lithium pegmatites from the leinster granite margin,southeast ireland: Geochemical constraints

Extensive and economically-important lithium pegmatites have been described from the eastern flank of the large S-type peraluminous Leinster Granite batholith in southeast Ireland. This flank area is also the locus of a major dip-slip ductile shear zone which was active during emplacement of the Tullow Lowlands pluton of the batholith. All the known bedrock occurrences of lithium pegmatite are contained within this marginal deformation zone. Most previous studies have suggested derivation of the lithium pegmatites by fractionation from the Leinster Granite itself, though none are conclusive. More recently, an anatectic model has been proposed whereby the pegmatites were derived by melting of Li-rich sediments adjacent to the batholithic contact. Our investigation has been concerned with a detailed study of the Rb-Sr isotope systematics and trace element geochemistry of the pegmatites and their host granitoids in an attempt to resolve the problem of pegmatite genesis. The results suggest a strong geochemical and temporal coherence between the Leinster Granite and the lithium pegmatites. A model is proposed in which Li enrichment in residual magmas may be linked to the development of muscovite instead of biotite in the Leinster Granite itself. The Li excluded from entry to biotite then accumulates in highly mobile, low-viscosity pegmatitic fluids which migrate into the marginal shear zone of the Tullow Lowlands pluton where they crystallize at considerably lower temperatures.     

2003年云南大姚6.2级和6.1级地震余震S波分裂研究

利用2003年云南大姚6.2、 6.1级两次地震的余震资料对该区进行了S波分裂研究。 结果表明, 云南大姚地区的快剪切波偏振的优势方向约为N33°W, 与该区域的区域主压应力的方向一致。 通过对比6.2级余震和6.1级余震的偏振方向和延迟时间, 可以看到6.1级地震后明显的应力释放现象, 在几次强余震前, 均观测到延迟时间呈现增加的特征, 并且几次地震在即将发生之前的短时间内还观测到延迟时间有减小的现象, 这一现象符合震前应力长时间积累和短时间应力释放的结论, 时间延迟在强震前会出现下降具有重要的地震短临预测意义。 研究认为S波分裂参数可以反映区域应力场的动态变化信息, 为应力场研究与地震预测提供有用的信息。     

各向异性介质中扭转波分裂的实验观测

实验室可以产生两种振动模式的横波,一种是剪切振动,另一种是扭转振动,在各向同性介质中两模式的横波速度是相同的,但它们的振动特性不一样,前者表现出很强的偏振特性,后者为无偏振特性.实验测试表明无偏振特性的扭转波在各向异性介质中传播时也会出现两种速度不同的扭转波,速度值与剪切横波的快慢横波速度值一致.用扭转波换能器接收时,这快慢扭转波的波形振幅不受各向异性方位影响.通过两块均匀的各向异性样品,用实验观测揭示了各向异性介质中扭转波的一些传播特征.     

辽东的主要剪切带及其金矿化特征

辽东地区构造上位于辽中—吉南元古宙活动带西南部,为胶辽克拉通金矿床分布比较集中的地区。金矿化作用发生于古元古宙变质沉积岩和花岗片麻岩的韧性及脆-韧性剪切带之中,剪切活动与成金事件具同时性关系。矿化剪切带分为3大系统,其中近东—西向和北东—南西向两大系统发育于辽河群分布区,其受控于岩石非均匀性,分别借助于早(吕梁)期S₁(≈S₀)和S₂或D₃剪切带而形成,剪切活动(构造事件)与成金事件主要发生于印支运动;丹东韧性剪切带及其相关构造为本区最重要的含矿系统,矿床均顺主构造线分布,并规律性地受控于其次级构造,剪切带的形成和成金事件均在燕山运动时期。本区提供了一个特殊的克拉通成矿作用格局:容矿岩石时代老,成金事件时代新。     

显式完全平方守恒差分格式在近岸海流数值模拟中的应用

用显式完全平方守恒差分格式及其改进分解算法对南海月平均流和海面起伏进行了数值模拟，与隐式完全平方守恒差分格式相比，计算时间可省３—５倍，具有良好的时间效益，而且，其计算效果不比隐式完全平方守恒差分格式差。因此，显式完全平方守恒差分格式及其改进分解算法具有良好的实用价值。     

Numerical simulation of nearshore circulation on field topography under random wave environment

In order to investigate surf zone hydrodynamics through two-dimensional numerical simulations of nearshore circulation under random wave environment, a nearshore circulation model, SHORECIRC, and a random wave model, SWAN, were combined and utilized. Using this combined model, a numerical simulation of the October 2, 1997 SandyDuck field experiment was performed. For this simulation, field topography and an input offshore spectrum were constructed using observed data sets synchronized with the experiment. The wave-breaking model in SWAN was modified by using breaker parameters varied according to bottom slope. The simulation results were compared with the experimental data, which revealed a well-developed longshore current, as well as with results using other combinations which were SHORECIRC and its original monochromatic wave-driver, and SHORECIRC and the default of SWAN. The results from the novel combined model agreed well with the experimental data. The results of the present simulation also indicate that alongshore field topography influences shear fluctuation of longshore currents.     

新型珊瑚海水混凝土力学性能试验研究

为了解决全珊瑚海水混凝土(coral aggregate seawater concrete,简称CASC)的脆性问题,采用碱式硫酸镁水泥(basic magnesium sulfate cement,简称BMSC)为胶凝材料和剑麻纤维为增韧材料制备了一种新型全珊瑚海水混凝土(NCASC),系统测定了基本力学性能,并与普通硅酸盐水泥制备的CASC、碱式硫酸镁水泥砂石混凝土(basic magnesium sulfate cement concrete,简称BMSCC)、普通砂石混凝土(ordinary portland cement concrete,简称OPC)以及普通轻骨料混凝土(lightweight aggregate concrete,简称LAC)的力学性能进行对比。结果表明:立方体抗压强度(f_(cu),m)在15~35 MPa范围内的NCASC的轴心抗压强度(f_c,m)、劈裂抗拉强度(f_(sp),m)与f_(cu),m之间具有较显著的线性关系;对于f_(cu),m为15~35 MPa的混凝土,相同强度NCASC的f_c,m高于OPC、LAC及BMSCC,其f_c,m与f_(cu),m比值最大,大致在0.85~0.96之间,NCASC的f_(sp),m比OPC和LAC高,但比BMSCC低,且随着强度的提高,OPC与NCASC的f_(sp),m间的差距在增大,而LAC和BMSCC与NCASC的f_(sp),m间的差距在缩小;在强度等级C15~C55的范围中,NCASC和CASC的f_(sp),m、f_c,m与f_(cu),m之间的关系基本可以采用相同的线性回归方程来描述,利用这两个回归线性方程得到对应于标准强度等级的CASC和NCASC的f_(sp),m和f_c,m数据,为CASC力学性能研究与结构设计提供参考依据。     

Implementation of baroclinic terms in a three-dimensional hydrodynamic model and its application to Anzali Port, Iran

Baroclinic terms have been implemented in a three-dimensional fully hydrodynamic model developed by Badiei et al. [2008. A three-dimensional non-hydrostatic boundary fitted model for free surface flows. International Journal for Numerical Methods in Fluids, 56(6), 607-627] modifying its momentum equations to account for density gradients and utilizing the scalar (salinity, temperature, etc.) conservation equation (SCE) and a state equation for the calculation of density. In the solution of advection-diffusion terms of the governing Navier-Stokes equations (NSE) and SCE, a symmetric splitting method was applied to ensure the long-term stability of simulations. Correction terms proposed by Ruddic et al. (1995) were applied to SCE to ensure the conservation of the scalar quantity. In the presence of baroclinic terms, the zero gradient pressure in the vertical direction in the vicinity of surface and bottom boundaries assumed by Badiei et al. [2008. A three-dimensional non-hydrostatic boundary fitted model for free surface flows. International Journal for Numerical Methods in Fluids, 56(6), 607-627] created spurious currents. This problem was solved by assuming a hydrostatic pressure variation at those boundaries. The ability of extended model was validated by comparing its results with an experimental test case. The simulation of hydrodynamic and salt intrusion at Anzali Port located at the southern coasts of Caspian Sea in Iran was carried out by the model with both barotropic and baroclinic modes. The simulated results with baroclinic mode show a better agreement with measured data as compared to the results of barotropic mode that clearly demonstrate the significance of baroclinic terms in the simulation of cyclic intrusion of salt wedge into the Port Basin.