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.     

Prediction of hydrodynamic loading on a mini TLP with free surface effects

This paper describes the extension of a fluid-flow simulations method to capture the free surface evolution around a full-scale Tension Leg Platform (TLP). The focus is on the prediction of the resulting hydrodynamic loading on the various elements of the TLP in turbulent flow conditions and, in particular, on quantifying the effects of the free surface distortion on this loading. The basic method uses finite-volume techniques to discretize the differential equations governing conservation of mass and momentum in three dimensions. The time-averaged forms of the equations are used, and the effects of turbulence are accounted for by using a two-equation, eddy-viscosity closure. The method is extended here via the incorporation of surface-tracking algorithm on a moving grid to predict the free-surface shape. The algorithm was checked against experimental measurements from two benchmark flows: the flow over a submerged semi-circular cylinder and the flow around a floating parabolic hull. Predictions of forces on a model TLP were then obtained both with and without allowing for the deformation of the free surface. The results suggest that the free surface effects on the hydrodynamic loads are small for the values of Froude number typically encountered in offshore engineering practice.     

Applied hydrodynamic wave-resistance computation by Fourier transform

An applied Fourier transform computation for the hydrodynamic wave-resistance coefficient is shown, oriented to potential flows with a free surface and infinity depth. The presence of a ship-like body is simulated by its equivalent pressure disturbance imposed on the un-perturbed free surface, where a linearized free surface condition is used. The wave-resistance coefficient is obtained from the wave-height downstream. Two examples with closed solutions are considered: a submerged dipole, as a test-case, and a parabolic pressure distribution of compact support. In the three dimensional case, a dispersion relation is included which is a key resource for an inexpensive computation of the wave pattern far downstream like fifteen ship-lengths.     

Micro-fracture instabilities in granular solids

We study the aspect of unstable behavior (like strain localization bands) in elastic solids as a consequence of micro-fracturing. A two-scale approach of computational homogenization is considered. The macroscopic behavior is investigated by finite element computations on a unit cell. At the micro-level, we consider a granular structure with elastic grains. The inter-granular boundaries are modeled with cohesive laws, friction and unilateral contact. We show that decohesion between grains gives rise to macro-instabilities, indicated by the loss of ellipticity, typical for deformation localization bands. The relation between the microscopic softening on inter-granular boundaries and the onset of macro-instabilities is studied through numerical examples. The influence of the cohesive law and friction parameters is analyzed. For periodic distributions of granular structures, we prove the loss of periodicity by failure and the corresponding size dependence effect in the homogenized response. We present numerical examples of bifurcation of solutions for granular cell structures and of particular solutions specific to elementary volumes with periodic cell distribution. Size dependence appears in the unstable regime and is strongly influenced by cohesion and friction parameters.     

大别山-苏鲁碰撞造山带构造几何学、运动学和岩石变形分析

按照构造几何学特点和运动学特征我们把大别山-苏鲁造山带的分为三个构造单元:南部,中部和北部。造山带南部为一套构造堆叠体系;中部为一个混合岩穹窿,浅变质的砂岩、板岩和片岩构成了大别山-苏鲁造山带的北部构造单元。造山带南部的构造堆叠体系主要由前陆褶皱带构成:未变质的新元古代-早三叠世的沉积地层;由“宿松群”北部和苏北地区的“海州群”构成的高压变质岩石单元及含柯石英和金刚石的超高压变质岩石单元。造山带中部的混合岩穹隆由大别山地区的罗田穹隆和苏鲁地区的莱西-栖霞穹隆构成。同样大别山北部的浅变质“佛子岭-卢镇关群”和胶东地区浅变质的“蓬莱群”构成了造山带北部的构造堆叠体系。同时大别山和苏鲁两个构造地体均经历了相似的多期构造变形:沿 NW-SE 向矿物拉伸线理发育的上部指北的剪切变形代表着造山带主变形期的变形;早期向南逆冲的韧性剪切变形和沿中部混合岩穹隆边缘发育的重力滑脱变形体系,后者代表了混合岩穹隆形成时的垂向缩短作用。正是由于构造几何学和多期变形的可对比性决定了这两个变质地体具有相同的地球动力学背景。     

Geostatistical inversion of coupled problems: dealing with computational burden and different types of data

Practical application of geostatistical inversion to coupled problems is hampered by a number of difficulties. In this paper, we address two of them: first, the computational cost of sensitivity (Jacobian) matrices and, second, the evaluation of the relative weights of different types of data. Regarding the first, we revise the adjoint state equations to propose a form whose cost is independent of the number of unknown parameters and only grows with the number of observation wells. Regarding the second, we derive expressions for the relative weights of different types of data. These expressions are based on minimizing the expected likelihood, rather than the likelihood itself. The efficiency of both improvements is tested on a synthetic example. The example analyzes a wide range of groundwater flow and solute transport conditions. Yet, the expected likelihood consistently yields the optimal weights. The proposed form of the adjoint state equations leads to one order of magnitude reduction in CPU time with respect to the conventional sensitivity equations.     

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

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

A critical CAE analysis of the bottom shape of a multi stepped air cavity planing hull

One of the most important aspects, in the ACS and multi stepped hull design, is the choice of the geometrical shape of the cavity and the steps. In this article a complete experimental and numerical campaign on a multi stepped ACS has been carried out, varying the velocity and the air flow rate under the hull. The experimental tests have been conducted in an ITTC Towing Tank allowing to validate the numerical simulations obtained by means of a CFD U-RANSe (Unesteady Reynolds Averaged Navier-Stokes equations) code. The CFD setup is described in detail. From this campaign a critical analysis of the Froude number influence on the air cushion shape has been argued. The authors identified four different behaviours, from low to very high Froude numbers. The use of CFD has allowed to observe quantities of difficult evaluation by means of traditional experimental test, as e.g. the frictional component of the resistance, the airflow path lines and the volume of fraction in transversal and longitudinal sections. The results have been discussed.     

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%.     

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

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