双三角翼拉升流场特性数值模拟研究

发展了适用于双三角翼大攻角非定常分离流场模拟的数值方法,开展双三角翼拉升运动的动态流场特性研究。通过减缩频率、转轴位置和起始攻角的变化,仔细分析了这些运动参数对动态流场施加影响的物理机制,有利于提高对双三角翼在拉升条件下的非定常特性和流场滞后效应等非线性现象的认识。     

高压储气罐泄漏喷射时冲击载荷数值计算

采用作者研发的用于含运动可变形物体流场计算的自适应非结构网格生成技术、动边界条件及数值处理方法，对地下高压储气罐发生泄漏喷射时罐盖（或破片）运动过程以及周围环境目标点上的冲击载荷进行了数值计算，讨论了在不同的罐盖质量、喷射孔径条件下目标点上冲击载荷的变化情况。研究结果表明,罐盖（或破片）的运动对爆炸场冲击载荷的影响是显著的。     

基于多孔介质模型的养殖网箱周围流场特性研究

为探讨鱼类养殖网箱对周围海域流场的影响,基于雷诺时均方程和Realizable k-ε湍流模型,耦合多孔介质模型,建立了研究鱼类养殖网箱周围流场的三维水动力模型.根据典型养殖场的实际数据,分别对仅网衣作用、仅鱼类作用及网衣和鱼类共同作用下网箱周围流场进行数值模拟.结果表明综合考虑网衣和鱼类的作用时,网箱下游流速衰减最大...     

水平布置等径气泡上升运动数值模拟分析

采用FTM(Front Tracking Method)方法,使用Fortran语言自编程序,对气泡上升运动进行数值模拟,并且选用小气泡作为研究对象,忽略气泡变形对气泡运动的影响。通过对比分析1~4个气泡的流场与运动速度变化,发现气泡间相互作用主要受流场变化的影响,由于双气泡的漩涡抵消和尾迹合并作用,双气泡的终极速度要比单气泡小,而三个气泡和四个气泡向上运动时中间气泡超越两侧气泡之后逐渐表现为单气泡和双气泡的运动特性,但发展到终极速度的时间更长。     

环形激波绕射、反射和聚焦流场的间断有限元模拟

采用间断有限元方法对环形激波在圆柱形激波管内绕射、反射和聚焦流场进行了数值模拟。将二维守恒方程的间断有限元方法发展到轴对称Euler方程,并对环形激波绕后台阶流动进行了数值计算。计算结果表明,采用间断有限元方法能够有效地捕捉运动激波在圆柱形激波管内传播的复杂流场结构;在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力,在聚焦点附近不会产生振荡或抹平间断现象。     

自由场水中爆炸气泡的物理特性

将水中爆炸气泡运动阶段周围流场假设为无粘、无旋、不可压缩的理想流体,运用边界元法模拟自由场中气泡的运动,在气泡运动模拟过程中引入数值光顺技术及弹性网格技术,避免因网格扭曲而导致的数值发散,并开发计算程序。计算值与实验值吻合良好,误差小于10%。从自由场水中爆炸气泡的基本现象入手,基于本文中开发的程序系统地研究了自由场中气泡的动力学特性。对流场中不同方位的压力进行分析,得出气泡中心的迁移方向及射流的攻击方向压力载荷比其他方向均大,说明气泡射流的攻击方向压力载荷最大,对水中结构造成严重毁伤,表明了气泡载荷的不对称性。计算了流场中不同位置的速度变化曲线,结果表明随着距气泡中心距离的增大,气泡运动引起的滞后流的速度迅速减小,且随着气泡的膨胀和坍塌,滞后流的方向逆转,总结了滞后流的衰减及变化规律。     

NUMERICAL ANALYSIS ON THE VELOCITY AND PRESSURE FIELDS INDUCED BYMULTI-OSCILLATIONS OF AN UNDERWATER EXPLOSION BUBBLE

假设水下爆炸气泡的内部气体在膨胀收缩过程中满足绝热条件,周围流体无黏无旋不可压缩. 基于势流理论,采用边界元法研究气泡动力学行为,重点关注气泡引起的流场脉动载荷以及滞后流特性,给出了相关的理论推导和数值计算方法. 通过将数值结果与解析解、实验值进行对比,数值模型的收敛性和有效性能够得到保证. 利用编写的程序进行计算和分析,发现在气泡加速膨胀阶段,流场压力在气泡径向不一定是逐渐衰减,还有可能以先增后减的规律变化;气泡射流后,为了能够继续描述环状气泡的运动以及流场特性,将此时的流场分为无旋场和一个布置在气泡内部涡环的叠加,计算过程中采用了一些数值技巧处理气泡的拓扑结构,得以连续模拟多个周期的气泡运动. 环状气泡具有相对较高的上浮迁移速度,而且在其顶部和底部附近分别形成两个高压区,顶部的高压区峰值相对较大,底部的高压区范围相对较大. 环状气泡中心轴上的流场速度会在气泡中心有一个加速过程,在气泡顶部附近又迅速减小.     

气泡多周期运动时引起的流场压力与速度

假设水下爆炸气泡的内部气体在膨胀收缩过程中满足绝热条件,周围流体无黏无旋不可压缩. 基于势流理论,采用边界元法研究气泡动力学行为,重点关注气泡引起的流场脉动载荷以及滞后流特性,给出了相关的理论推导和数值计算方法. 通过将数值结果与解析解、实验值进行对比,数值模型的收敛性和有效性能够得到保证. 利用编写的程序进行计算和分析,发现在气泡加速膨胀阶段,流场压力在气泡径向不一定是逐渐衰减,还有可能以先增后减的规律变化;气泡射流后,为了能够继续描述环状气泡的运动以及流场特性,将此时的流场分为无旋场和一个布置在气泡内部涡环的叠加,计算过程中采用了一些数值技巧处理气泡的拓扑结构,得以连续模拟多个周期的气泡运动. 环状气泡具有相对较高的上浮迁移速度,而且在其顶部和底部附近分别形成两个高压区,顶部的高压区峰值相对较大,底部的高压区范围相对较大. 环状气泡中心轴上的流场速度会在气泡中心有一个加速过程,在气泡顶部附近又迅速减小.      

基于动态面搭接技术的直升机旋翼流场分析

为了进一步提高直升机旋翼流场计算的精度与效率,本文根据直升机旋翼的运动特点,发展了一种基于结构化网格动态面搭接技术及低速预处理技术的直升机旋翼流场的非定常数值分析方法。采用该方法分别对Caradonna-Tung、UH-60A旋翼模型的悬停流场及NASA某旋翼无升力前飞流场进行了数值分析;将计算所得的旋翼表面压力分布与风洞试验结果进行了对比,对旋翼的尾涡结构进行了相关分析。分析结果表明:本文所发展的方法具有较高的精度且计算速度快,并能够比较准确地捕捉直升机旋翼的尾涡结构,可以用于直升机旋翼流场的分析研究,从而为开展考虑挥舞及变距运动的直升机旋翼前飞状态非定常气动特性的研究奠定了基础。     

淹没磨料射流的空泡运动分析

通过对淹没条件下磨料射流的空泡运动分析研究,建立了淹没磨料射流的空泡运动方程,揭示了淹没磨料射流中空泡的溃灭特性,数值模拟了淹没磨料射流的磨料体积浓度以及空泡所处流场压力对空泡运动及溃灭的影响规律.分析表明:淹没磨料射流中磨料的存在增大流体的粘性系数,增大空泡溃灭历时,减弱射流的空蚀破坏能力;流场压力的改变对空泡溃灭过程影响显著,压力越高,空泡溃灭历时越短.     

Development of a parallel storm surge model

A new parallel storm surge model, the Parallel Environmental Model (PEM), is developed and tested by comparisons with analytic solutions. The PEM is a 2‐D vertically averaged, wetting and drying numerical model and can be operated in explicit, semi‐implicit and fully implicit modes. In the implicit mode, the propagation, Coriolis and bottom friction terms can all be treated implicitly. The advection and diffusion terms are solved with a parallel Eulerian–Lagrangian scheme developed for this study. The model is developed specifically for use on parallel computer systems and will function accordingly in either explicit of implicit modes. Storm boundary conditions are based on a simple exponential decay of pressure from the centre of a storm. The simulated flooding caused by a major Category 5 hurricane making landfall in the Indian River Lagoon, Florida is then presented as an example application of the PEM. Copyright © 2003 John Wiley & Sons, Ltd.     

Calculation of wind currents in the edge zone of a sea ice sheet

The steady wind currents in the ocean near the edge of an ice sheet simulated by a rigid cover are studied. The characteristic horizontal scales of the problem correspond to several tens of kilometers. In this case, depending on the depth and viscosity of the fluid, the effect of the Coriolis force on its motion may be comparable with the effect of the wind shear stresses applied to the ice-free surface. Taking the Coriolis force into account leads to the formation of a spiral flow structure in the vertical direction.Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.6, pp. 132–141, November–December, 1992.     

波浪载荷导致黄河口潮坪沉积物垂向运移现场观测研究

2005年8月7~8日,现代黄河三角洲刁口地区经历了一次由台风"麦莎"引起的风暴潮过程。通过对比分析风暴潮前后受保护潮坪滩面沉积物的粒度特征变化情况,发现了风暴潮期间在强烈的波浪载荷导致的渗流作用下,粉质土海床出现了极细粉粒由下向上运移并输出滩面的现象。结合现场试验期间采集到的孔隙水压力监测数据,本文基于海床动力响应的观点分析了其成因机制。     

波-流相互作用的缓坡方程及其波作用量守恒

当表面波从开阔海域传播至近岸水域时，普遍的波一流相互作用经受着海底的强烈影响．运用水波Hamilton变分原理，建立了近岸水域任意水深变化海底上波一流相互作用的缓坡方程．它可包含波、流和水深一般变化的二阶效应，约化为某些典型的缓坡型方程．据此得出广义程函方程，并且证明该缓坡方程的波作用量守恒．     

The effects of mean sea level rise and strengthened winds on extreme sea levels in the Baltic Sea

Mean sea level rise and climatological wind speed changes occur as part of the ongoing climate change and future projections of both variables are still highly uncertain. Here the Baltic Sea's response in extreme sea levels to perturbations in mean sea level and wind speeds is investigated in a series of simulations with a newly developed storm surge model based on the nucleus for European modeling of the ocean(NEMO)-Nordic. A simple linear model with only two tunable parameters is found to capture the changes in the return levels extremely well. The response to mean sea level rise is linear and nearly spatially uniform, meaning that a mean sea level rise of 1 m increases the return levels by a equal amount everywhere. The response to wind speed perturbations is more complicated and return levels are found to increase more where they are already high. This behaviour is alarming as it suggests that already flooding prone regions like the Gulf of Finland will be disproportionally adversely affected in a future windier climate.     

Modelling the non-linear interaction of wind and tide: Its influence on current profiles

A single-point model in the vertical is used to examine the coupling between tidal currents and wind-driven flows in shallow near-coastal regions. Calculations using both a linear slip and a no-slip condition at the sea bed clearly show that coupling between tidal and wind-driven currents cannot occur in a linear model with a time-independent eddy viscosity. However with a physically more realistic time-varying viscosity related to the flow field, coupling does occur, the magnitude of this non-linear interaction depending upon the change in eddy viscosity over a tidal cycle and the intensity of shear in the vertical. A point model in the vertical with flow induced by an oscillatory pressure gradient and an additional constant wind stress is used to examine the influence of viscosity parametrization and water depth upon this coupling. The solution in the vertical is accomplished using both a functional approach and a finite difference method. Some conclusions as to the relative merits of these approaches, particularly the use of a transformed grid in the case of high-shear surface and bed boundary layers, are made in the paper.     

Normal and typhoon wind loadings on a large cooling tower: A comparative study

Large cooling towers are sensitive to wind effects with their increasing heights and flexibilities. Unlike traditional approaches, which employed Code-defined normal winds to check the loading characteristics, this paper developed a framework for checking the typhoon-induced wind loading on a large cooling tower using Monte Carlo simulations and multi-fan wind tunnel tests. Some distinct characteristics of typhoon winds were compared with those of Code-defined normal winds. Furthermore, wind characteristics of incoming normal and typhoon winds in terms of vertical profile of mean wind speed, turbulence intensity, turbulence integral length scale and power spectrum density of fluctuating winds were well reproduced by a feedback control process of a multi-fan actively controlled wind tunnel. The surface wind pressure distributions of a large cooling tower under these conditions were then investigated by testing a 1:600 reduced scale model. Mean and fluctuating external wind pressures along the circumferential direction under various incoming winds were discussed and quantitatively formulated with eight-term trigonometric equations. Moreover, the cross correlations of wind pressures in the circumferential and meridian directions and correlations with structure base forces, i.e. integral drag and lift forces, were investigated. Non-Gaussian characteristics in terms of skewness and kurtosis of fluctuating wind pressures were also analyzed under two wind climates. Peak factors for modeling extreme wind pressures were examined and compared with those of various models. Finally, the extreme wind loads on a large cooling tower obtained from different wind pressure combinations were compared with peak-factor-theory-based results to identify an appropriate combination for structural design.     

On the Effect of Friction on Wind Driven Shallow Lakes

We investigate the steady motion of a liquid in a lake, modeled as a thin domain. We assume the motion is governed by Navier--Stokes equations, while a Robin-type traction condition, and a friction condition is prescribed at the surface and at the bottom, respectively. We also take into account Coriolis forces. We derive an asymptotic model as the aspect ratio d \delta = depth/width of the domain goes to 0. When the Reynolds number is not too large, this is mathematically justified and the three-dimensional limit velocity is given in terms of wind, bathymetry, depth and of a two-dimensional potential. Numerical simulation is carried out and the influence of traction condition reading is experienced.     

Impact of Atmospheric Front Parameters on Free and Forced Oscillations of Level and Current in the Sea of Azov

Level and current oscillations in the basin of the Sea of Azov have been studied by hydrodynamic modeling using the Princeton ocean model (POM). The hypothesis on the role of the resonance mechanism in the occurrence of extremely large amplitudes of storm surge and seiche oscillations depending on the velocity and time of motion of atmospheric fronts of the Sea of Azov has been tested. It is found that at the same wind, pressure perturbations moving over the Sea of Azov induce forced oscillations, and after the perturbations cease, free oscillations with amplitudes that are 14% higher than those obtained at constant atmospheric pressure. It is shown that the motion of the atmospheric front (whose velocity and time are selected under the assumption that waves with maximum amplitudes are generated) plays an important but not decisive role in the formation of the structure of currents and level oscillations in the Sea of Azov.     

阻抗式和简单调压室甩荷时水位波动的显式计算方法

目前对调压室水位波动的研究主要有图解法和数值计算。尚未看到考虑阻力项后调压室水位波动随时间变化的显函数计算方法。本文从调压室的基本微分方程出发，采用泰勒级数展开式推导了阻抗式和简单式调压室甩荷时水位波动的显式计算方法。运用这种方法不仅可以计算调压室水位波动的过程，而且可以方便的计算调压室的最高涌浪和第二振幅。