沿岸海域三维斜压场的数值模拟——Ⅰ、渤海潮流数值计算

本文建立了渤海海域三维斜压场潮流模型。该模型考虑了热、盐效应及密度的空间变化,引入了计算网格无法分辨的湍运动能量,并以M_2分潮为例进行了计算,与实测值比较,获得了比较满意的结果。并对潮致欧拉余流垂直分量的空间变化进行了探讨。     

湍致次生环流的数值研究

本文给出了一个数值模式,用于描述具有横向粗糙度变化的平面上方的湍流场内的次生环流。在此模式中,雷诺平均运动方程组的闭合,采用了代数化的雷诺应力输运方程和k-ε二方程湍流模式。计算结果和实验结果相吻合。本文还讨论了利用数值模式进行参数实验的一些结果。     

抛弃式探头下沉运动数值计算的网格划分方法研究

各类抛弃式探头下沉运动的数值计算研究都需要对探头表面和整个计算域进行网格划分的前处理过程,选择不同的网格划分方法和划分精度,会影响数值计算速度和计算精度,最终生成的网格质量决定了计算结果的收敛性和准确性。针对抛弃式探头的复杂结构,对比了两种网格划分方法、三种网格划分精度,完成探头计算区域的网格划分,采用k-ε湍流模型,进行抛弃式探头下沉运动的数值计算,并将计算结果与水箱和水库实验结果进行对比,验证了混合网格划分方法和普通精度网格对抛弃式探头下沉运动数值计算的适用性,研究结果对类似较为复杂结构的水下运动体的数值计算前处理过程具有一定的参考和借鉴意义。     

湖／陆风过程中湍流摩擦作用的数值分析

本文利用一个湍流能量闭合的二维原始方程模式分析讨论了湍流交换对湖／陆风过程的影南。文中较详细地介绍了湍能闭合方法．分析了在湖(陆)风发展的不同阶段湍流交换强度的分布特征并讨论了大尺度风场对湍流交换系数场的影响。     

三维数值模式用于海底摩擦系数的研究

本文讨论了在风暴潮和潮汐作用下海底摩擦力与平均流速的关系公式,即;利用湍流能量损耗理论得到了一个新的求拖曳系数C_b的公式;并且通过一个三维模式,在理想的矩形海区中计算了各种条件下的C_b分布。     

湖/陆风过程中湍流摩擦作用的数值分析(英文)

本文利用一个湍流能量闭合的二维原始方程模式分析讨论了湍流交换对湖/陆风过程的影响。文中较详细地介绍了湍能闭合方法,分析了在湖(陆)风发展的不同阶段湍流交换强度的分布特征并讨论了大尺度风场对湍流交换系数场的影响。     

XCP投弃式海洋探头阻力系数数值计算及试验研究

以XCP探头为模型,采用有限体积法结合κ-ε两方程湍流模型,对XCP探头阻力系数进行了数值计算,研究了XCP探头的阻力特性。获得了不同雷诺数下阻力系数的变化趋势和水下稳定运动时的极限速度。结果表明阻力系数随着雷诺数的增加而减小,其数值计算结果与海上实测数据基本吻合,验证了该数值计算方法的正确性,可为XCP探头结构设计提供参考。     

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

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

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

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

船上海气之间湍流通量的观测研究

用６个加速度计组阵系统来监测船体姿态和运动，校正船上观测到的风湍流资料，然后用涡动相关技术计算海气之间动量，热量和水汽通量，结果表明，本观测系统可精确地消除船体运动对风湍流的影响，并获得精确的海气湍流通量。     

A Three-Dimensional Prediction Method for Thermal Diffusion

A three-dimensional,first order turbulence closure,thermal diffusion model is described inthis paper.The governing equations consist of an equation of continuity,three components of momentum,conservation equations for salt,temperature and subgridscale energy,and an equation of state.In the mod-el,according to the hypothesis of Kolmogorov and Prandtl,the viscosity coefficient of turbulent flow ofhomogeneous fluid is related to the local turbulent energy,and the horizontal and vertical exchangecoefficients of mass,heat and momentum are computed with the introduction of subgridscale turbulenceenergy.The governing equations are solved by finite difference techniques.This model is applied to theJiaozhou bay to predict thermal pollution by the Huangdao power plant.An instantaneous tidal currentfield is computed,then the distribution of temperature increment is predicted,and finally the effect of windstress on thermal discharge is discussed.     

海浪破碎对海洋上混合层中湍能量收支的影响

海浪破碎产生一向下输入的湍动能通量,在近海表处形成一湍流生成明显增加的次层,加强了海洋上混合层中的湍流垂向混合。为了研究海浪破碎对混合层中湍能量收支的影响,文中分析了海浪破碎对海洋上混合层中湍流生成的影响机制,采用垂向一维湍封闭混合模式,通过改变湍动能方程的上边界条件,引入了海浪破碎产生的湍动能通量,并分别对不同风速下海浪破碎的影响进行了数值研究,分析了混合层中湍能量收支的变化。当考虑海浪破碎影响时,近海表次层中的垂直扩散项和耗散项都有显著的增加,该次层中被耗散的湍动能占整个混合层中耗散的总的湍能量的92.0%,比无海浪破碎影响的结果增加了近1倍;由于平均流场切变减小,混合层中的湍流剪切生成减小了3.5%,形成一种存在于湍动能的耗散和垂直扩散之间的局部平衡关系。在该次层以下,局部平衡关系与壁层定律的结论一致,即湍动能的剪切生成与耗散相平衡。研究结果表明,海浪破碎在海表产生的湍动能通量影响了海洋上混合层中的各项湍能量收支间的局部平衡关系。     

Yangtze Brackish water plume—circulation and diffusion

A system of 3-D linearized momentum equations, the continuity equation and a simplified equation of state for sea water are analytically solved. The solution comprises three parts. The first part is a wind-driven current, the second a jet-like current caused by Yangtze River outflow and the last a density flow. The computed current velocity is obtained by solving an advection-diffusion equation of salinity by a finite-difference method. The results show that the Yangtze brackish water plume spreads as a large low salinity tongue. At the surface, the turning of the tongue axis is mainly controlled by wind stress, Yangtze run-off and bottom topography. Near the bottom, the axis changes little all the year around. At about 10–15 m depth, there is a sharp halocline in summer, but at deeper layers the vertical distribution is almost always uniform. The computed salinity distribution agrees fairly well with the observed one.     

非线性流函数在实验风浪场计算中的应用

风作用于水面产生风浪, 其中由于波流紊动产生的动量和能量的交换机制是一个很复杂的过程。风应力一般用来描述这种能量交换, 可以分为3个部分: 水面的剪切力、波生应力以及紊动应力。采用一种有效的非线性波流分离方法——NSFM(Nonlinear Stream Function Method)对波流运动的动量和能量输移进行定性描述。构造能够有效表达非线性波浪的解析流函数, 摄动求解使其满足拉普拉斯方程、动力边界条件和运动边界条件, 结合实验室风浪数据, 分离出波生速度场。通过交叉谱分析, 得到波生雷诺应力在不同风速下对风应力的贡献。结果表明: NSFM对不同工况条件下的风浪的处理具有较高的精度, 模型适应性良好; 且风速越大, 波生应力沿着水深衰减得越快, 且自由面波生应力在动量输移中的比重会逐渐减弱。     

An Extended Mild-Slope Equation

On the assumption that the vortex and the vertical velocity component of the current aresmall,a mild-slope equation for wave propagation on non-uniform flows is deduced from the basichydrodynamic equations,with the terms of (V_hh)~2 and (V_h~2)h included in the equation.The terms of bot-tom friction,wind energy input and wave nonlinearity are also introduced into the equation.The wind en-ergy input functions for wind waves and swells are separately considered by adopting Wen′s(1989)empiri-cal formula for wind waves and Snyder′s observation results for swells.Thus,an extended mild-slope equa-tion is obtained,in which the effects of refraction,diffraction,reflection,current,bottom friction,wind en-ergy input and wave nonlinearity are considered synthetically.     

A pair of Langmuir cells in a laboratory tank (I) wind-only experiment

Three velocity components of subsurface flow, observed in a rectangular tank under the action of a constant wind speed, are measured systematically at mesh points distributed uniformly over a vertical cross-section of the tank. Measurements are carried out for two cases: 1) reference wind speedU _r =7.5 m/s and fetchF=10 m; and 2)U _r =10 m/s andF=25 m. A pair of Langmuir cells is observed for both cases; downwelling zones are found along both of the sidewalls and an upwelling zone in the centre of the tank. Near the water surface, the vertical momentum flux is dominated by the Reynolds stress resulting from small-scale turbulence, while over the entire cross-section except near the surface, the Reynolds stress due to the Langmuir cells dominates the vertical momentum flux. As the result of the occurrence of this Langmuir cells, the vertical momentum flux, which consists of both mean advection and small-scale turbulence, is markedly inhomogeneous in the spanwise direction; for example, the largest vertical flux of the order of the wind stress is observed in the downwelling zone near one sidewall, while at the centre of the tank, the vertical momentum flux occupies only 30% of the wind stress. This indicates that a pair of Langmuir cells plays more important role than small-scale turbulence in the mixing process in a greater part of the wind-wave tank.Address after April 1, 1992: Department of Civil Engineering, Hiroshima Institute of Technology, Miyake 2-1-1, Saeki-ku, Hiroshima 731-51, Japan.     

Three-parameter model of turbulence for the atmospheric boundary layer over an urbanized surface

A modified three-parameter model of turbulence for a thermally stratified atmospheric boundary layer (ABL) is presented. The model is based on tensor-invariant parametrizations for the pressure-strain and pressure-temperature correlations that are more complete than the parametrizations used in the Mellor-Yamada model of level 3.0. The turbulent momentum and heat fluxes are calculated with explicit algebraic models obtained with the aid of symbol algebra from the transport equations for momentum and heat fluxes in the approximation of weakly equilibrium turbulence. The turbulent transport of heat and momentum fluxes is assumed to be negligibly small in this approximation. The three-parameter E ? ε ? 2> model of thermally stratified turbulence is employed to obtain closed-form algebraic expressions for the fluxes. A computational test of a 24-h ABL evolution is implemented for an idealized two-dimensional region. Comparison of the computed results with the available observational data and other numerical models shows that the proposed model is able to reproduce both the most important structural features of the turbulence in an urban canopy layer near the urbanized ABL surface and the effect of urban roughness on a global structure of the fields of wind and temperature over a city. The results of the computational test for the new model indicate that the motion of air in the urban canopy layer is strongly influenced by mechanical factors (buildings) and thermal stratification.     

Modelling of the upper turbulent layer in the ocean

A differential model of the upper turbulent layer in the ocean is considered. A closed system of equations includes equations of motion, balance, and dissipation of kinetic turbulence energy. Boundary conditions at the surface are determined using a solution of the atmospheric problem taking into account the interaction between the two media. The formulated algorithm allows for a relationship between turbulent energy dissipation and flux and the parameters of wind disturbance. The vertical profiles of turbulence and drift current characteristics are presented as well as parameters of the ocean-atmosphere interaction for various values of impulse jump within the limits of the wave layer with waves collapsing and not collapsing.UDK 551.456.152     

On the application of a turbulence closure modified model to the description of the density jump evolution in a stably stratified medium

The self-similar turbulent density jump evolution has been studied in the scope of a turbulence closure modernized theory which takes into account the anisotropy and mutual transformation of the turbulent fluctuation kinetic and potential energy for a stably stratified fluid. The numerical calculation, performed using the equations for the average density and kinetic and potential energies of turbulent fluctuations, indicates that the vertical profiles of the buoyancy frequency, turbulence scale, and kinetic and potential energies drastically change when the turbulence anisotropy is strong. The vertical profiles of the corresponding energy and spatial discontinuity parameters, calculated at a weaker anisotropy, indicate that similar drastic changes are absent and a qualitative agreement exists with the known analytical solution, which describes the density jump evolution in a freshwater basin and was obtained previously [5, 8] in the scope of a turbulence local-similarity hypothesis applied in combination with the budget equation for the turbulent fluctuation kinetic energy.