Existence of transverse waves in anisotropic poroelastic media

Out of the four waves in an anisotropic poroelastic medium, two are termed as quasi-transverse waves. The prefix 'quasi' refers to their polarizations being nearly, but not exactly, perpendicular to direction of propagation. In this composite medium, unlike perfectly elastic medium, the propagation of a longitudinal wave along a phase direction may not be accompanied by transverse waves. The existence of a transverse wave in anisotropic poroelastic media is ensured by the two equations restricting the choice of elastic coefficients of porous aggregate as well as fluid–solid coupling. Necessary and sufficient conditions for the existence of transverse waves along the coordinate axes and in the coordinate planes for general anisotropy are discussed. The discussion is extended to the case of orthotropic materials and existence for few specific phase directions is also explored. The conditions for the transverse waves decided on the basis of their apparent polarizations, that is, particle motion being perpendicular to ray direction, are also discussed. For a particular numerical model, the existence of these apparent transverse waves is solved numerically for phase directions in coordinate planes. For general directions of phase propagation, the existence of these transverse waves is checked graphically for the chosen numerical model.     

Wedge impact with the influence of ice

This work presents a canonical study on a wedge entering water near a single piece of ice using computational-fluid-dynamics (CFD) and a Wagner-type theoretical model with corrections for non-linear effects. Calculations for a series of conditions with ice of different sizes and locations relative to the wedge are conducted. The hydrodynamic force due to impact, the pressure distribution on the wedge surface, and the pile-up phenomenon are examined to study the role of ice in the impact process. The theoretical model is shown to be accurate and can serve as a useful method to assess slamming loads under the influence of ice. It is shown that even for the case of a small piece of ice, the slamming force on the wedge can increase by 30%.     

Validation of a finite difference method for the simulation of vortex-induced vibrations on a circular cylinder

The Karman Vortex Street generated by a circular cylinder is investigated by the numerical solution of the compressible Navier–Stokes equations in the incompressible Mach number range (Mach<0.3) using the Beam and Warming implicit scheme. The agreement with the fully incompressible projection method (Chorin, 1968) is fairly good while convergence time is very much better. The investigation suggests that the compressible Navier–Stokes equations may be used as an efficient alternative to study incompressible flows as well. Mach numbers just below 0.3 are enough to simulate incompressible flow behavior and at the same time do not cause numerical ill-conditioning in the solution. In addition, some relevant features of the vortices generated and carried by the wake of the cylinder could be fairly well captured.     

构造应力环境对浅成岩体成矿系统的制约:从安庆月山岩体冷却过程动力学计算模拟结果分析

为探讨构造应力环境对浅成岩体成矿系统的制约,以安庆地区的月山岩体及其所形成的安庆铜矿床为例,在传统地质调查研究的基础上,通过假设不同的力学边界条件进行岩体冷却过程动力学的计算模拟实验研究。计算结果显示,在没有构造应力作用时,在岩体的接触带上不可能形成有利成矿定位的局部汇流扩容空间,挤压和拉伸作用都能够在岩体接触带形成局部汇流扩容空间,但二者的空间分布规律有着极大的差别,只有拉伸构造应力作用下的动力学计算模拟实验能在已发现矿体部位形成汇流扩容空间,表明构造应力环境对成矿定位的汇流空间的制约以及拉伸构造应力作用对月山岩体成矿的重要贡献。在拉伸构造应力和无构造应力作用时,接触带上的流体压力都是随着变形的进程而趋降低的,只有在挤压构造应力作用下,流体压力才会有一个超压的过程,流体超压是形成斑岩型矿床的重要机制。可见挤压环境才是一种有利于斑岩型矿床成矿的构造环境,铜陵-安庆地区缺乏与矽卡岩矿床相对应的斑岩型矿床主要原因之一是成矿时处于一种拉伸的构造应力作用环境。     

车排子凸起与沙湾凹陷过渡带的构造模型

车排子凸起位于准噶尔盆地扎伊尔山前,是夹持在南部四棵树凹陷和东部沙湾凹陷之间的三角构造带。然而对这个三角构造带是原地掩伏构造,山前推覆构造,还是走滑构造控制形成的认识存在争议。目前对于车排子凸起与沙湾凹陷过渡带的构造变形特征以及变形机制的认识尚不清楚。而过渡带是研究该区构造变形与演化的良好区域,精细研究该区域构造几何学与运动学对于揭示其构造变形机制以及认识盆山之间的关系具有积极作用,并对于该区油气勘探具有重要意义。本文以研究区的连片三维地震资料为基础,运用不整合面识别、轴面分析及平衡剖面等方法,对车排子凸起与沙湾凹陷过渡带进行精细的构造解析、建模与构造复原。研究表明,过渡带发育下二叠统（P₁j-P₁f）,中二叠统（P₂x-P₂w）,上二叠统（P₃w）,三叠系（T）,侏罗系（J）,白垩系（K）与古近系-第四系（E-Q）7个构造变形层序。发育佳木河组底部及石炭系内部滑脱层。深部（前二叠系）构造楔发育于整个过渡带,其个数和几何学形态的变化直接影响了上覆背斜形态。构造楔在早二叠世开始楔入,挤压作用持续至三叠纪末期。晚三叠世发育生长地层,显示了P-T构造层褶皱变形的时间,随后三叠纪末期,过渡带受到来自西北缘造山带方向的逆冲推覆构造的叠加作用,表现在后缘发育一系列叠瓦构造,推覆距离约为8 km,南部逆冲断层上盘二叠系-三叠系被剥蚀完全,北部残余少量下二叠统佳木河组。     

基于CFD的单个建筑物风场模拟参数选择

运用流体力学软件FloEFD对沽源单个建筑物周围的风场进行数值模拟。通过不断改变模型中计算参数的设置进行一系列模拟试验,对比模拟试验结果,并与观测资料进行比较,分析不同计算参数对模拟结果的影响,并获得适用于该模型的最佳参数。主要研究的计算参数包括计算域高度,初始网格等级,局部初始网格等级和不同平均风速剖面形式。结果表明:计算域高度从3倍建筑物高度开始,空腔区的长度、漩涡中心位置以及再发展区的边界位置基本保持稳定。随着初始网格等级的增加,空腔区的长度、再发展区的边界位置及计算时间逐渐增大。局部初始网格等级对模拟结果影响不显著。以两种不同平均风速剖面形式进行模拟,迎风漩涡长度不同,背风面影响不大。与观测资料比较显示,最优参数组合为:计算域高度为3倍建筑物高度,初始网格等级为4、局部初始网格等级为4、平均风速剖面形式为指数律。     

Application of a semi-spectral cloud water parameterization to cooling tower plumes simulations

In order to simulate the plume produced by large natural draft cooling towers, a semi-spectral warm cloud parameterization has been implemented in an anelastic and non-hydrostatic 3D micro-scale meteorological code. The model results are compared to observations from a detailed field experiment carried out in 1980 at Bugey (location of an electrical nuclear power plant in the Rhône valley in East Central France) including airborne dynamical and microphysical measurements. Although we observe a slight overestimation of the liquid-water content, the results are satisfactory for all the 15 different cases simulated, which include different meteorological conditions ranging from low wind speed and convective conditions in clear sky to high wind and very cloudy. Such parameterization, which includes semi-spectral determination for droplet spectra, seems to be promising to describe plume interaction with atmosphere especially for aerosols and cloud droplets.     

Influence of Thermal Effects on the Wind Field Within the Urban Environment

Micrometeorological conditions in the vicinity of urban buildings strongly influence the requirements that are imposed on building heating and cooling. The goal of the present study, carried out within the Advance Tools for Rational Energy Use towards Sustainability (ATREUS) European research network, is the evaluation of the wind field around buildings with walls heated by solar radiation. Two computational fluid dynamics (CFD) codes were validated against extensive wind-tunnel observations to assess the influence of thermal effects on model performance. The code selected from this validation was used to simulate the wind and temperature fields for a summer day in a specific region of the city of Lisbon. For this study, the meteorological data produced by a non-hydrostatic mesoscale atmospheric model (MM5) were used as boundary conditions for a CFD code, which was further applied to analyze the effects of local roughness elements and thermodynamic conditions on the air flow around buildings. The CFD modelling can also provide the inflow parameters for a Heating, Ventilation and Air Conditioning (HVAC) system, used to evaluate the building energy budgets and to predict performance of the air-conditioning system. The main finding of the present three-dimensional analyses is that thermal forcing associated with the heating of buildings can significantly modify local properties of the air flow.     

Methods for locating the proper position of a planned fishway entrance near a hydropower tailrace

The spatial distribution of upstream migrating Atlantic salmon (Salmo salar) spawners was studied in 2008 and 2009 in the surroundings of the tailrace from a hydropower station in the River Umeälven. This area is problematic because the fish have difficulties finding their way into the original riverbed, which prompted an investigation on the feasibility of adding a fishway in the area. Echo sounding was used in 2008 to investigate the spatial use of fish in the tailrace channel during the time of salmon migration. Presence of other fish species was so low that all echo sounding detections were assumed to be salmon. In 2009, data on wild radio-tagged salmon (n = 94) was collected in the same tailrace by an array of underwater antennas with a detection radius of approximately 10 m, to validate the results from the echo sounding. Both types of surveys showed aggregations of fish in one part of the tailrace. As a final step, Computational Fluid Dynamics (CFD) modeling was performed to analyze hydrodynamics. This CFD modeling showed a coincidence of fish detections in areas with turbulence intensities between 0.6 and 0.8, which may be an indication that the fish are holding in these areas to reduce energy expenditure during migration. A high proportion of the radio-tagged salmon were observed in the tailrace area (a median of 21 days between their first and last detection), indicating that the salmon are delayed in their upstream migration by the attraction to the tailrace in this area. The overall data on fish behavior in the tailrace led us to suggest a location for a new fishway where the fish aggregation was most pronounced. The number of detections from either technology had similar distributions over the tailrace, suggesting that the results are comparable. Thus, the split beam echo sounder can preferably be used to investigate which location is most appropriate for locations of entrances to new fishways since many more individual fishes are covered by this method compared to telemetry. Furthermore, there is no need to handle the fish in the echo sounding studies as is required in telemetry studies.     

地震波动方程的局部间断有限元方法数值模拟

近年来,随着地震波数值模拟对计算精度和效率的要求越来越高,间断有限元方法开始受到越来越多的关注.本文中,针对具有吸收边界条件的二维地震声波波动方程,作者提出了一种基于局部间断有限元方法的数值模拟算法.该算法在空间上使用局部间断有限元方法进行离散,在时间上采用了显式蛙跳格式.在这种时空离散的组合方式下,每个时间步上,此算法在空间剖分的每个单元上的求解计算是相互独立的,因而具有极高的并行性.通过数值算例,我们将该算法与连续有限元方法进行了比较.结果表明,本算法不仅具有对起伏构造的良好适应性,而且在计算效率和计算精度等方面,都具有优越性.