FINITE ELEMENT ANALYSIS OF A VERTICAL RECTANGULAR PLATE COUPLED WITH AN UNBOUNDED FLUID DOMAIN ON ONE SIDE USING A TRUNCATED FAR BOUNDARY

The dynamic pressure distribution on a rectangular plate attached to a rigid wall and supporting an infinitely large extent of fluid subjected to a harmonic ground excitation is evaluated in the time domain. Governing equations for the fluid domain are set considering the compressibility of the fluid with negligibly small change in density and a linearized free surface. A far boundary condition for the three-dimensional fluid domain is developed so that the far boundary is truncated at a closer proximity to the structure. The coupled problem is solved independently for the structure and the fluid domain by transferring the acceleration of the plate to the fluid and pressure of the fluid to the plate in sequence. Helmholtz equation for the three-dimensional fluid domain and Mindlin's theory for the two-dimensional plate are used for the solution of the interacting domains. Finite element technique is adopted for the solution of this problem with pressure as nodal variable for the fluid domain and displacement for the plate. The time dependent equations are solved in each of the interacting domain using Newmark-b method. The effectiveness of the technique is demonstrated and the influences of surface wave, exciting frequency and flexibility of the plate on dynamic pressure are investigated.     

THREE-DIMENSIONAL LARGE EDDY SIMULATION OF FREE SURFACE TURBULENT FLOW IN OPEN CHANNEL WITHIN SUBMERGED VEGETATION DOMAIN

1 .　INTRODUCTIONInactualrivers ,theriverbediscommonlycoveredbyvegetation .Thevegetationinriversplaysanimportantroleontheenvironmentalfunc tionoftherivers,thatis,theamenitytohumanlife ,livingspaceformanylives ,purificationofwaterqualityandsoon .However ,thevegetationinriversalsoplaysakeyroleonthecarryingca pacityofachannelwithfloodplains .Therefore ,thehydraulicsoftheflowwithvegetationshouldbemoreimprovedorrefined .Logicalexplanationandaccurateevaluationaboutturbulentflowinopenchannelwith…     

均匀来流中大长径比深海立管涡激振动特性

提出了一种对大长径比深海立管涡激振动特性进行三维耦合数值模拟的方法。其中,流场采用大涡模拟方法进行三维数值模拟,结构振动采用基于薄壳模型的有限元方法进行数值计算,并采用一种将流场和结构响应数据进行实时传输的新方法,实现了流体与固体之间的三维耦合数值模拟。在雷诺数Re=3400和11000下,对某长径比为1200的深海立管模型在均匀来流中的横向涡激振动特性进行了数值模拟与分析。结果表明,大长径比柔性立管的横向振动响应对其尾涡动力特性会产生显著影响,包括会出现涡泄频率显著增大和多频涡泄等复杂现象。尽管来流是均匀的,但立管的横向振动响应会出现明显的多模态特征,除了包含有与涡泄频率一致的振动响应外,还包含有其它频率成分的振动响应。特别地,还会出现多种高阶模态振动共存的非对称复杂弯曲变形等现象。研究表明,该方法为研究深海立管涡激振动特性及其工程预报的相关问题提供了一种有效的途径。     

大流量进水工况下主隧系统水气运动特性

为给深层隧道排水系统工程的主隧系统结构设计和工况控制提供依据,采用数值模拟与物理模型试验相结合的方法,建立与物理模型体型结构相同,包含两竖井和相连主隧段的主隧系统三维数值模型,利用该数值模型对两种大流量极端工况下的进水过程进行三维动态数值模拟.结果表明:主隧系统在大流量下的进水过程可归纳为水流从两端竖井进入主隧,水流碰...     

基于CFX的离心泵内部流场三维数值模拟

基于三维不可压缩流体的N-S方程和标准k-ε湍流模型,采用隐式修正SIMPLE算法,利用CFX软件,对离心泵装置进行全流场三维数值模拟,共计算额定转速下30～80 m3/h流量范围内8个工况点。对比分析小流量工况、设计工况和大流量工况下泵装置流态和压力分布,并分析叶片表面静压,揭示内部流动规律,所得结果对预测水力性能,提高水泵效率及进一步结构优化具有重要的参考价值。     

三维裂隙网络非稳定渗流数值分析

研究裂隙岩体中渗透水流随时间、边界条件的变化规律,从而客观反映渗透水流情况,具有重要的意义。本文针对分布特别稀疏的三维主干岩体裂隙网络,以渗流区内的水量平衡原理为依据建立基本方程,推导三维岩体裂隙网络非稳定渗流的基本方程;并将三维裂隙网络渗流理论及有限元数值方法运用于主干岩体裂隙网络渗流分析中。基于裂隙岩体网络非稳定渗流数学模型,编制相应的程序。对于工程实例,根据实测情况利用Monte-Carlo模拟出岩体中二级裂隙网络分布情况进行分析,然后进行了三维非稳定渗流分析。计算结果充分体现了岩体中裂隙网络的储水性,当边界条件发生变化时,裂隙岩体中水头分布变化存在一定的滞后性;并表明主干裂隙的强导水性,使得裂隙岩体中的主干裂隙的渗流分析成为渗流分析中的重点。     

Time domain simulations of radiation and diffraction by a Rankine panel method

The radiation and the diffraction of a ship with a forward speed are studied by using a time domain Rankine panel method. The free surface conditions are linearized onto an undisturbed free surface based on the double body flow. The linearized body boundary condition is applied on the mean wetted hull surface. The fluid domain boundary is discretized by a collection of quadric panels. The unknown quantities, including the free surface elevation, the normal flux over the free surface and the potential on the fluid domain boundary, are determined at each time step. The numerical results are compared with experimental data and other numerical solutions, showing satisfactory agreements.     

A three-dimensional fictitious domain method for the simulation of fluid-structure interactions

The Distributed-Lagrange-Multiplier based Fictitious Domain (DLM/FD) method previously proposed by Yu [J. Comput. Phys. 207 (2005), 1–27] for the simulation of fluid/flexible-body interaction problems is extended from the 2D to the 3D case in this study. We first demonstrate that the Lagrange multiplier problem can be more efficiently solved with a direct-forcing scheme instead of the original Uzawa iterations without the sacrifice of the accuracy. The type of the interpolation function (i.e. smoothed delta function) for the transfer of the quantities between the Eulerian and Lagragian frames is shown not to affect the results significantly. Secondly, the fictitious domain method is implemented in the three-dimensional case. The new 3D codes are applied to the flapping of a 3D flexible plate in a uniform flow, and the deformation of tri-leaflets at opening stage.

     

泥石流柔性防护系统破坏原因分析与设计对策

泥石流柔性防护系统对于中小型泥石流有很好的防治效果,但在实际工程中很容易发生破坏。为研究泥石流柔性防护系统破坏原因,介绍了泥石流柔性防护系统的构成与优势、破坏原理、破坏形式,并针对四川省道217某段公路改建工程的泥石流柔性防护工程,分析了系统的设计、构件、基础等方面的问题,得出主要破坏原因是设计不合理,未完全按照定型化设计,基础发生破坏,构件发生破坏。同时对泥石流柔性防护工程的设计提出了一些建议,以期提高其拦截能力。     

溢流面板坝泄槽锚杆设计参数敏感性分析

针对设锚杆的溢流面板坝泄槽流激振动性能复杂的特点,建立拟三维仿真模型,采用时域分析法对锚杆设计参数的敏感性进行了研究。结果表明:在一定的水流限制流速条件下(v≤40 m/s),在影响溢流面板坝流激振动性能的3个锚杆设计参数中,锚杆布设间距对关键点的拉应力和接缝处法向变形响应幅值影响最大,锚杆布设长度及直径影响相对较小。大致规律为:关键点的应力随锚杆布设间距和锚杆直径的增加而增大,随布设长度的增大而减小;关键点的法向变形随板后布设长度和锚杆直径的增加而增大,随锚杆布设间距的增加而减小;3种因素存在显著性顺序:锚杆布设间距,布设长度,锚杆直径。研究成果可为同类坝型泄槽锚杆参数的设计提供指导。     

不同长细比圆柱绕流的大涡模拟

自从人们对层流圆柱绕流现象有了系列研究及清楚的认识后,逐渐把目光投向湍流的圆柱绕流问题,但高雷诺数下的圆柱绕流具有很强的三维特性,因此在数值计算时模型轴向长度的选取是大家关心的问题。该文基于开源软件OpenFOAM运用LES计算了4种不同长细比圆柱的静止圆柱绕流(包括一个二维算例),重点分析和对比了结果中的一些基础参数如St数、平均阻力系数等,与实验结果吻合良好。文中还进一步分析了时均顺流向速度U在流场中沿流向和横向的分布,发现长细比为π的圆柱算例中,计算域的轴向尺度已经能满足圆柱下游靠近圆柱处的流场中的三维结构的发展。     

浮体的非线性辐射效应

本文在空间上应用面源法，在时间域上应用步进法，研究了浮体的非线性辐射效应，本文探讨了：（1）不同形式的辐射条件对非线性辐射问题解的影响；（2）不同距离的辐射面对开边界条件的使用所产生的影响；（3）大振幅的浮体受近运动所产和的非线性效应。     

胥浦活水泵站肘形进水流道流态分析及优化

胥浦活水泵站为一座长江边低扬程引水泵站,设计流量5.0 m3/s,最大扬程3.3 m。为适应其低扬程、小流量的特点,选用了全贯流潜水泵,在原胥浦节制闸底板上改建安装。为了保证在长江低水位时,水泵进水口完全淹没,进水流道型线平顺,流道内无涡带或其他不良流态,机组启动正常和运行稳定,利用三维数值建模、边界条件设置与三维湍流流动模拟,运用计算流体动力学CFD方法,进行进水流道内部流动数值模拟,多方案分析和比较进水流道内部流态、水泵进水条件和流道的水力损失。计算结果表明,进口尺寸1.1 m×3.4 m(高×宽),出口直径1.2 m的型线B方案进水条件好,流道水头损失小,满足水泵高效运行的要求。优化设计方案为泵站安全运行提供最优的设计进水流道型线和设计参数。     

三维薄翼入水水动力的计算

本文提供了一个计算三维薄翼高速入水水动力性能的非定常数值升力面方法。假设薄翼高速入水时,翼的背面为全空泡状态,当翼的随边进入液面后,此空泡域按一定的方向拖至自由液面。在翼的拱弧面上布置线涡和线源,在翼后空泡面上布置线源,所有的奇点系在自由液面上方存在映象以满足自由液面条件。根据翼迎流面的法向不可穿透条件和翼背面及翼后空泡面上的空泡动力学条件确定所有待定线涡和线源的强度,并由此可以计算各时刻翼上压力分布。利用本文提供的方法首先计算了二维薄翼入水的水动力性能,将结果与理论解进行比较,两者结果较为一致。然后,研究了三维薄翼的入水问题。     

Turbulent structure and bursting process in multi-bend meander channel

The flow characteristics in a meander channel are fully three-dimensional. With the primary flow in the streamwise direction, the secondary flow in the transverse and vertical directions induced by the channel bends are significant in the analyses of the turbulent structure. Some of the analyses in the straight channel, for instance, the quadrant analysis for the bursting phenomena, are inadequate in investigating the meander channel, since the flow in the transverse direction is not taken into account. In order to reveal the flow characteristics in a multi-bend meander channel, especially, the bursting process, experiments are conducted in the present study. With the three-dimensional quadrant analysis, the influence of the transverse flow velocity during the bursting process could be correctly addressed. The analyses and discussions are presented in this paper.     

适合多目标种群的窄深型仿自然鱼道研究

传统的仿自然鱼道断面以宽浅型为宜,但对于某些水头较高、河道狭窄的水电站,为了降低工程造价,并增加过鱼目标种群数量,需要建成窄深型的仿自然鱼道。窄深型仿自然鱼道内各项水力参数受边壁及内部结构影响较大,水流结构复杂,易产生流速集中、旋涡等不利流场,为了形成满足多种群鱼类上溯的水流条件,需要开展全面的研究。按照将流场“垂向分层、横向分区”的研究思路,有效利用窄深型断面的特点,结合理论分析、三维数值模拟和物理模型试验的方法,分层分区地满足不同鱼类的上溯条件,在解决了工程问题的基础上,实现了对窄深型仿自然鱼道水力学特性的初步探索研究,增加了仿自然鱼道建设的选择多样性。     

THE SIMULATION AND ANALYSIS OF TANK SLOSHING WITH POROSITY GIRDER BY MULTI-DOMAIN BOUNDARY ELEMENT METHOD

The simulations of nonlinear sloshing in both two-dimensional and three-dimensional liquid tanks with porosity girder are carried out in time domain by using the Multi-domain Boundary Element Method (MBEM). The comparison of the present results and the references is conducted and shows they agree well with each other. The three-dimensional effect and the influence of the width of liquid tank and the height and porosity of the girder on the natural period of the model are also discussed.     

RECENT PROGRESS IN INCOMPRESSIBLE REYNOLDSAVERAGED NAVIER-STOKES SOLVERS

ＲＥＣＥＮＴＰＲＯＧＲＥＳＳＩＮＩＮＣＯＭＰＲＥＳＳＩＢＬＥＲＥＹＮＯＬＤＳＡＶＥＲＡＧＥＤＮＡＶＩＥＲ－ＳＴＯＫＥＳＳＯＬＶＥＲＳ￥Ｃ．Ｈ．Ｓｕｎｇ；Ｔ．Ｔ．Ｈｕａｎｇ（ＤａｖｉｄＴａｙｌｏｒＭｏｄｅｌＢａｓｉｎ，ＣＤ／ＮＳＷＣＢｅｔｈｅｓｄａ，Ｍ...     

NUMERICAL SIMULATION OF 3-D TURBULENT FLOW IN THE MULTI- INTAKES SUMP OF THE PUMP STATION

In this article, a numerical model for three-dimensional turbulent flow in the sump of the pump station was presented. A reasonable boundary condition for the flow in the sump with several water intakes at different flow rates was proposed. The finite volume method was employed to solve the governing equations with the body fitted grid generated by the multi-block grid technique. By using the Fluent software, the fluid flow in a model sump of the pump station was calculated. Compared with the experimental result, the numerical result of the example is fairly good.     

Method of coupling 1-D unsaturated flow with 3-D saturated flow on large scale

A coupled unsaturated-saturated water flow numerical model was developed. The water flow in the unsaturated zone is considered the one-dimensional vertical flow, which changes in the horizontal direction according to the groundwater table and the atmospheric boundary conditions. The groundwater flow is treated as the three-dimensional water flow. The recharge flux to groundwater from soil water is considered the bottom flux for the numerical simulation in the unsaturated zone, and the upper flux for the groundwater simulation. It connects and unites the two separated water flow systems. The soil water equation is solved based on the assumed groundwater table and the subsequent predicted recharge flux. Then, the groundwater equation is solved with the predicted recharge flux as the upper boundary condition. Iteration continues until the discrepancy between the assumed and calculated groundwater nodal heads have a certain accuracy. Illustrative examples with different water flow scenarios regarding the Dirichlet boundary condition, the Neumann boundary condition, the atmospheric boundary condition, and the source or sink term were calculated by the coupled model. The results are compared with those of other models, including Hydrus-1D, SWMS-2D, and FEFLOW, which demonstrate that the coupled model is effective and accurate and can significantly reduce the computational time for the large number of nodes in saturated-unsaturated water flow simulation.