近海湍流三维数值模拟中的重整化群模式

利用重整化群理论来确定近海湍流数值模式中的垂向涡动粘性系数，并对冬末春初琼州海峡的潮流场特征进行数值模拟。结果表明，重整化群模式的待定常数少，应用方便可行，其计算结果也比常用的湍流封闭模式的计算结果更加合理。     

湍流数值模拟中封闭模式应用的局限性

通过对琼州海峡的潮流场特征进行数值模拟，指出了选择不同的特征混合长度表达式对数值模拟结果的影响，表明了基于特征混合长度理论的流封闭模式在近海湍流数值模拟中应用的局限性。     

海洋湍流模式应用研究

海洋湍流模型研究自二十世纪 70年代中发展至今 ,在海洋动力学研究中 ,特别是关于混合 /层化研究中广泛应用 ,近年来由于认识到湍流对海洋生物过程的重要影响 ,对海洋湍流的客观描述更加关注。文中详细介绍了几个主要海洋动力学模型中的湍流封闭模式 ,如 HAMSOM中的Prandtl混合长模型、Johns模式中的 k-方程模型、POM中的 k- kl模型、水动力学中常用的 k-ε模型等等 ,介绍了海洋湍流模型的应用。对于湍流模型的使用提出针对具体问题选择的原则 ,复杂的并非最优的。     

基于大涡模拟与被动示踪物模型的人工鱼礁数值研究

投放人工鱼礁是解决海洋生态环境问题的重要措施之一。通过使用并行大涡模拟模式（the parallelizeda large-eddy simulation model,PALM）及被动示踪物模型模块,研究了不同流速条件（0.1,0.2,0.4和0.6 m/s）下方型人工鱼礁对流场形态、营养盐的抬升作用、和湍流动能收支的影响。研究表明,鱼礁的存在使得其附近垂向速度增大,产生上升流。受到上升流的抬升作用,鱼礁底部的示踪物迅速进入海洋上层,之后遇到鱼礁后方的背涡流,示踪物的抬升受阻,高度逐渐降低。在上升流区域以及背涡流区域的共同影响下,示踪物抬升区域的最大高度与来流流速无关。不同上升流定义对应的上升流区域的高度与来流流速均不相关;上升流区域的最大速度、平均速度与来流流速都成线性增加的关系;然而随着来流流速的增大,不同上升流定义对应的上升流区域面积的变化趋势却完全不同。投放鱼礁后,鱼礁区域底部的混合增强,区域底部的能量被输运至上层。这说明,鱼礁的存在不仅对营养物质具有抬升作用,还能将能量向上输运。     

斜压海洋动力学的一种三维数值模式：Ⅱ.温度,盐度和垂直涡动粘性系数的计？

根据海水温度和盐度平流扩散方程给出一种数值计算方案并采用混合长度理论给出生趣涡动粘性系数的计算方法。对于温度和盐度方程其平流过程采用了Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式,水平扩散彩 显著 式,生趣扩散采用隐格式。     

上混合中剪切湍流和朗缪尔环流动力特征差异

Large eddy simulation(LES) is used to investigate contrasting dynamic characteristics of shear turbulence(ST)and Langmuir circulation(LC) in the surface mixed layer(SML). ST is usually induced by wind forcing in SML. LC can be driven by wave-current interaction that includes the roles of wind, wave and vortex forcing. The LES results show that LC suppresses the horizontal velocity and greatly modifies the downwind velocity profile, but increases the vertical velocity. The strong downwelling jets of LC accelerate and increase the downward transport of energy as compared to ST. The vertical eddy viscosity Km of LC is much larger than that of ST. Strong mixing induced by LC has two locations. They are located in the 2ds–3ds(Stokes depth scale) and the lower layer of the SML,respectively. Its value and position change periodically with time. In contrast, maximum Km induced by ST is located in the middle depth of the SML. The turbulent kinetic energy(TKE) generated by LC is larger than that by ST. The differences in vertical distributions of TKE and Km are evident. Therefore, the parameterization of LC cannot be solely based on TKE. For deep SML, the convection of large-scale eddies in LC plays a main role in downward transport of energy and LC can induce stronger velocity shear(S2) near the SML base. In addition, the large-scale eddies and S2 induced by LC is changing all the time, which needs to be fully considered in the parameterization of LC.     

基于大涡模拟耦合模式的小尺度海气相互作用研究

大气和海洋是影响地球气候系统的两个重要因素,它们之间的相互作用是海洋和大气研究的重要课题,海气耦合模式则是研究海气相互作用的重要工具,而海气耦合模式重点考虑的参数是海气通量。针对传统的大尺度海气耦合模式缺少湍流尺度分析的问题,本文使用并行大涡模拟海气耦合模式（The Parallelized Large-Eddy Simulation Model,PALM）,在小尺度上探究风速对海气通量及湍流动能收支（Turbulence Kinetic Energy Budget,TKE Budget）的影响,设置了5、10和15m/s三种地转风速度对大气边界层（Atmospheric Boundary Layer,ABL）和海洋混合层（Oceanic Mixed Layer,OML）进行海气耦合模拟。研究表明：海气通量的分布与风速大小密切相关,风速越大,净热通量和浮力通量相对越大,由于温度上升导致海水蒸发加剧,使得大气的淡水通量增大;海洋湍流动能收支各项在近海面处受风速影响较大,且随着深度加深而逐渐减弱。本研究初步展示了小尺度海气耦合模式在海气通量研究中的应用,对进行小尺度海气相互作用研究具有一定的意义。     

基于大涡模拟的三角形地形之层结流体湍动能收支模型

为了探究海底地形对湍流动能收支的影响,本文使用并行大涡模拟模式(The Parallelized Large-Eddy Simulation Model,PALM),以坡陡作为无量纲地形参数(δ),设置了亚临界、临界和超临界三种地形状态(δ=0.5,1和2)进行数值模拟。文章计算并呈现了地形作用下的流体速度和湍流动能收支分布,以及湍流动能平衡方程各参量。通过回归分析和量纲分析重点讨论了地形顶点处耗散和海表处能量的耗散,与地形坡陡的关系,得出其关系均呈指数形式。     

地形和风速影响下的海气相互作用大涡模拟研究

当水流经过海洋地形时,水流的不稳定性会引起垂向混合并伴随大量湍流过程。针对传统海气耦合模式缺少在湍流尺度上讨论海洋地形与风速对海气相互作用影响的问题,使用并行大涡模拟海气耦合模式(the parallelized large eddy simulation model, PALM)在5 m/s的背景风场下,引入理想立方体地形,对比有无地形的影响;设置地形边长为L,高为3L (其中大气部分高L), L与水深H之比为L/H=1/2;然后保持地形条件不变。设置5、10和15 m/s三种风速,讨论风速对小尺度海气相互作用的影响。研究表明:地形在大气部分减弱顺风向速度,增强侧风向速度,影响0～5L的高度区域,而对垂向作用较小;无地形条件下湍流垂向涡黏系数K_m在-0.3L时,水深达到最大值0.024 m²/s,有地形条件下K_m在-0.8L时,达到最大值为0.16 m²/s,地形的存在使得上层海洋混合加强, K_m最大值增加1个数量级。随风速增大海洋和大气中的净热通量、淡水通量和浮力通量都相应...     

不同湍流模型在波浪水槽数值模拟中的 应用与比较

为了研究不同湍流模型在建立数值水槽中的适用性,选用计算流体力学软件,采用S-A、标准的k-着及RNG k-着 等8种不同的湍流模型对水槽中的波浪进行了数值模拟,通过对比得到：RNG k-着模型不仅具有较高的精度,而且又节省计 算资源,较其它模型更具优势。     

Numerical Simulation of Sloshing Using the MPS-FSI Method with Large Eddy Simulation

A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.     

伴随同化方法在中国近海海洋数值模拟中的应用

简述了几种常用的资料同化方法,并较为详细地阐述了伴随同化方法的基本原理。基于其在潮汐模式中、在海水温度场中和在海洋生态模型中的应用等三方面综述了伴随方法在中国近海海洋数值模拟中的应用研究现状和进展。并对伴随方法的应用前景进行了初步的评述,指出该方法可对海洋观测方案的优化发挥重要作用,其在生态学流动力学模式中的应用亦应得到充分重视。     

Large Eddy Simulation of Open Ocean Deep Convection with Application to the Deep Water Formation in the East Sea (Japan Sea)

Various important features could be found on the open ocean deep convection and the subsequent deep water formation from large eddy simulation (LES), and the results were applied to the East Sea (Japan Sea). It was found that under a strong cold wind outburst with the heat flux of 1000 Wm⁻² for 5 days generates a deep convection which can penetrate to the depth 1500 m, but under the continuous cooling with the heat flux of 250 Wm⁻² the growth of a mixed layer is suppressed at 700 m. The effects of the spatial and temporal variations of the surface forcing were investigated with regard to the penetrative depth of convection, the generation of baroclinic eddies, the volume of the water mass formation, and the intensity of the rim current. The deep water formations in the intermediate and deep layer of the East Sea were explained in terms of the simulation results, and the intensity of the consequent circulation and the volume of water mass formation were compared with the observation data. This revised version was published online in August 2006 with corrections to the Cover Date.     

Large Eddy Simulation for Wave Breaking in the Surf Zone

In this paper, (he large eddy simulation method is used combined with the marker and cell method to study the wave propagation or shoaling and breaking process. As wave propagates into shallow water, the shoaling leads lo the increase of wave height, and then at a certain position, the wave will be breaking. The breaking wave is a powerful agent for generating turbulence, which plays an important role in most of the fluid dynamic processes throughout the surf zone, such as transformation of wave energy, generation of near-shore current and diffusion of materials. So a proper numerical model for describing the turbulence effect is needed. In this paper, a revised Smagorinsky subgrid-scale mode! is used to describe the turbulence effect. The present study reveals that the coefficient of the Smagorinsky model for wave propagation or breaking simulation may be taken as a varying function of the water depth and distance away from the wave breaking point. The large eddy simulation model presented in this pape     

Large Eddy Simulation for Plunge Breaker and Sediment Suspension

Breaking waves are a powerful agent for generating turbulence that plays an important role in many fluid dynamical processes, particularly in the mixing of materials. Breaking waves can dislodge sediment and throw it into suspension, which will then be carried by wave-induced steady current and tidal flow. In order to investigate sediment suspension by breaking waves, a numerical model based on large-eddy-simulation (LES) is developed. This numerical model can be used to simulate wave breaking and sediment suspension. The model consists of a free-surface model using the surface marker method combined with a two-dimensional model that solves the flow equations. The turbulence and the turbulent diffusion are described by a large-eddy-simulation (LES) method where the large turbulence features are simulated by solving the flow equations, and a subgrid model represents the small-scale turbulence that is not resolved by the flow model. A dynamic eddy viscosity subgrid scale stress model has been used for the     

A Study of Eddy Viscosity Coefficient in Numerical Tidal Simulation

Based on the fluid motion equations, the physical meaning of eddy viscosity coefficient and the rationality of the Boussinesq hypothesis are discussed in this paper. The effect of the coefficient on numerical stability is analyzed briefly. A semi-enclosed rectangular sea area, with an orthogonal spur dike, is applied in a 2-D numerical model to study the effect of horizontal eddy viscosity coefficient (AH). The computed result shows that AH has little influence on the tidal level and averaged flow velocity, but has obvious influence on the intensity and the range of return flow around near the spur dike. Correspondingly, a wind-driven current pool and an annular current are applied in a 3-D numerical model respectively to study the effect of vertical eddy viscosity coefficient (Av). The computed result shows that the absolute value of Av is inversely proportional to that of horizontal velocity, and the vertical gradient value of Av determines the vertical distribution of horizontal velocity. The distrib     

联合海浪过程的数值模拟

早期,海浪被抽象成一些简单的模式(如正弦波、有限振幅波等等),而加以研究,也得到了许多能近似符合实际的结果。但是,对真实海浪而言,其内部结构和外观特征都是复杂的,且都受着随机性的制约。因此,自五十年代初期开始,人们已广泛地使用随机函数(随机过程)来作为描述海浪的模式。当用随机函数这种模式研究海浪时,海浪谱是一种有效的工具,而如何获得海浪谱便构成海浪研究中重要的问题。海浪谱可以通过多种途径获得,但利用海浪过程的现实(观测记录)求谱是最基本、最可靠的方法,这方面国内外已进行了许多的研究。     

大气湍流引起飞机颠簸的理论分析和数值试验

基于飞机载荷因数变量方程、无量纲的大气垂直运动方程和热力学方程,建立起描述大气湍流和飞机颠簸的湍流模型,并对该模型进行理论分析和数值试验.研究指出,负的理查森数(Ri)对应着静力不稳定和动力不稳定流动,在不稳定大气层结条件下,当Ri<-Ra/PrRe2时,对流状态的大气运动必将由对流转变成湍流,并对运行于其中的飞机产生飞机颠簸.但在稳定大气层结条件下,当理查森数大于临界理查森数时(Ri>Ric=-Ra/PrRe2),正的Richardson数是动力稳定的,非对流状态的大气运动表现为重力波.当理查森数小于临界理查森数(O相似文献   

北印度洋半日潮波的数值模拟研究

基于FVCOM(Finite Volume Coast and Ocean Model)模型,建立北印度洋海域(31^°~102^°E,16^°S~31^°N)的M₂和S₂分潮潮波数值模式,研究北印度洋半日潮潮汐、潮流分布特征。对底摩擦系数进行数值试验,利用代价函数梯度下降法,得到分潮调和常数向量均方根偏差(RMSE)的变化曲线,逼近并确定最优的底摩擦系数。将采用该系数的模拟结果与TOPEX/Poseidon卫星高度计交叉点的调和常数数据、国际海道测量组织(IHO)及部分文献中的验潮站数据进行比较与验证,一致性较好。其中对比卫星数据的振幅偏差为2~4 cm、迟角偏差为7^°~8^°,与验潮站数据的振幅偏差为3~6 cm、迟角偏差为8^°~9^°。根据模拟结果,分析了北印度洋海域M₂和S₂分潮潮波传播特征和潮流椭圆的空间分布特征等。M₂分潮潮波在阿拉伯海南部有1个无潮点,在波斯湾内有2个无潮点,最大振幅超过80 cm;潮流在西北印度洋和孟加拉湾中部大多为顺时针旋转,其余海域大多为逆时针旋转;流速在阿拉伯海东北部、安达曼海、波斯湾和孟加拉湾北部较大,最大流速为160 cm/s,其他海域较小。S₂分潮的潮波传播特征、无潮点的位置和潮流椭圆的空间分布特征等都与M₂分潮类似,但潮波振幅和潮流流速等都相对M₂分潮较小。研究完善了北印度洋海域2个主要半日分潮M₂和S₂的整体特征。