三维波流数值水槽的构建方法

为了模拟波浪与水流的联合作用,基于Reynolds平均的Navier-Stokes方法,构建了一个三维波流数值水槽. 使用该模型对不同出流边界条件进行了对比分析. 结果表明:引入合理的入流和出流边界条件,可以保证在计算稳定后,水槽平均水面基本保持不变. 垂向时均流速分布的计算结果与实验数据吻合良好. 这表明:该波流数值水槽具有较好的波流特性,可用于计算波浪与水流、以及与三维固定结构之间相互作用.

可压缩湍流直接数值模拟

本书是以作者及其研究小组多年来的研究工作为基础撰写的，主要是可压缩湍流的直接数值模拟。全书共分两部分。第一部分为数值模拟方法：针对可压缩湍流的多尺度非定常非稳定流动特征，所构造的高精度差分格式（高于三阶精度），分析了各种高精度差分格式的数值行为特征和模拟小尺度物理结构、分辨非定常间断面的能力，重点分析了正确模拟多尺度物理结构及其非线性干扰对数值方法的要求。第二部分是物理问题的直接数值模拟：根据可压缩湍流的流动物理特征，所建立的求解可压缩Navier-Stokes方程的高精度数值模拟方法，详细描述了不稳定流和可压缩湍流的直接数值模拟，以及在此基础上所发现的一些新的物理现象。直接数值模拟主要针对典型可压缩湍流，包括均匀各向同性湍流、槽道湍流、平面混合流、射流、平板边界层湍流和钝楔、钝锥边界层湍流等，重点是可压缩湍流从层流转捩到湍流的直接数值模拟及其流动机理的研究。书中还介绍了我课组所研制的高精度、高效并行算法的软件。

Enhancements of the simulation method on the edge effect in resin transfer molding processes

In resin transfer molding processes, small clearances exist between the fiber preform and the mold edges, which result in a preferential resin flow in the edge channel and then disrupt the flow patterns during the mold filling stage. A mathematical model including the effect of cavity thickness on resin flow was developed for flow behavior involving the interface between an edge channel and a porous medium. According to the mathematical analysis of momentum equations in a fully developed rectangular duct and formulations of the equivalent edge permeability, comparing with three-dimensional Navier-Stokes equations, the governing equations were modified in the edge channel. The volume of fluid (VOF) method was applied to track the flow front. A simple case is numerically simulated using the modified governing equations. The effects of edge channel width and cavity thickness on flow front and inlet pressure are analyzed, and the evolution characteristics of simulated results are in agreement with the experimental results. (c) 2007 Elsevier B.V. All rights reserved

Reactive Mold Filling in Resin Transfer Molding Processes with Edge Effects

Reactive mold filling is one of the important stages in resin transfer molding processes, in which resin curing and edge effects are important characteristics. On the basis of previous work, volume-averaging momentum equations involving viscous and inertia terms were adopted to describe the resin flow in fiber preform, and modified governing equations derived from the Navier-Stokes equations are introduced to describe the resin flow in the edge channel. A dual-Arrhenius viscosity model is newly introduced to describe the chemorheological behavior of a modified bismaleimide resin. The influence of the curing reaction and processing parameters on the resin flow patterns was investigated.

Numerical wave flume study on wave motion around submerged plates

Nonlinear interaction between surface waves and a submerged horizontal plate is investigated in the absorbed numerical wave flume developed based on the volume of fluid (VOF) method. The governing equations of the numerical model are the continuity equation and the Reynolds-Averaged Navier-Stokes (RANS) equations with the k-epsilon turbulence equations. Incident waves are generated by an absorbing wave-maker that eliminates the waves reflected from structures. Results are obtained for a range of parameters, with consideration of the condition under which the reflection coefficient becomes maximal and the transmission coefficient minimal. Wave breaking over the plate, vortex shedding downwave, and pulsating flow below the plate are observed. Time-averaged hydrodynamic force reveals a negative drift force. All these characteristics provide a reference for construction of submerged plate breakwaters.

Effects of buoyancy on Langmuir circulation

Based on the Navier-Stokes equation, an equation describing the Langmuir circulation is derived by a perturbation method when the influences of Coriolis force and buoyancy force are both considered. The approach used in the analysis is similar to the works carried out by Craik and Leibovich [J. Fluid Mech. 73 (1976) 401], Leibovich [J. Fluid Mech. 79 (1977) 715] and Huang [J. Fluid Mech. 91 (1979) 191]. Potential applications of the equation proposed are discussed in the area of Antarctic circumpolar current.

浮力对Langmuir环流的影响

Langmuir环流是发生在上层海洋中纵向的螺旋状涡旋运动，其轴向与风向基本一致。它是一常见的重要海洋现象，与海洋研究的许多重要问题，如海－气交换、海洋混合、海洋环境、海洋遥感、海洋生态和海洋灾害等有着十分密切的关系。 本文在同时考虑科氏力与浮力作用下，从Navier-Stokes方程出发，应用摄动理论，导出一个包含科氏力与浮力影响的Langmuir环流控制方程组。所采用的方法类似于Craik and Leibovich（1976）、Leibovich（1977）及Huang（1979）等人研究Langmuir环流时所采用的研究方法。此外，本文还讨论了所导出的方程组在南极绕极流区域的潜在应用。

Second-order Stokes wave solutions for intefacial waves in three-layer stratified fluid with background current

In this paper, interfacial waves in three-layer stratified fluid with background current are investigated using a perturbation method, and the second-order asymptotic solutions of the velocity potentials and the second-order Stokes wave solutions of the associated elevations of the interfacial waves are presented based on the small amplitude wave theory, and the Kelvin-Helmholtz instability of interfacial waves is studied. As expected, for three-layer stratified fluid with background current, the first-order asymptotic solutions (linear wave solutions), dispersion relation and the second-order asymptotic solutions derived depend on not only the depths and densities of the three-layer fluid but also the background current of the fluids, and the second-order Stokes wave solutions of the associated elevations of the interfacial waves describe not only the second-order nonlinear wave-wave interactions between the interfacial waves but also the second-order nonlinear interactions between the interfacial waves and currents. It is also noted that the solutions obtained from the present work include the theoretical results derived by Chen et al (2005) as a special case. It also shows that with the given wave number k (real number) the interfacial waves may show Kelvin-Helmholtz instability.

Stokes二阶波对圆柱的绕射

本文介绍绕圆柱的二阶Stokes波的简单计算方法,并将它与其他计算结果作比较。结果表明本方法有简单及足够的精度的优点。

Stokes波高阶谐波系数递推的数值计算

本文用W.H.Hui提出的方法,在半物理平面内重新表述了Stokes波的数学模型和边界条件,提出了两种更有效的数值计算方法来获得Stokes波高阶谐波系数,并可递推至无穷。通过小参数转换,重新得到了Cokelet(1977)的波速和半波高的摄动展开式。

Stokes流的积分方程法

<正> Stokes流,或称零雷诺数流,指的是尺寸微小、速度缓慢的流动。它的理论在化工、生物力学、物理化学、环境保护、选矿、地球物理和气象科学等各个领域都有重要的应用。 零雷诺数流可用Stokes方程来描述:式中μ,V和P分别是流体的粘度、速度向量和压力。直到本世纪60年代,只有数目非常有

非线性Kelvin-Helmholtz不稳定性

本文用导数展开法对液体薄层与亚音速气流接壤时的界面稳定性作非线性分析.文中考虑了液体的表面张力与体积力,故非线性的Rayleigh-Taylor不稳定性可作为特例而导出;液体与气体均不计粘性.虽然Nayfeh曾算过这一情况,但其三阶方程有遗漏(如213页的式(2.29)).同时解也不自洽(如其一阶解(2.31)并不满足他的初始条件(2.20)),此外,在截止波数附近,对行波他并未考虑.本文弥补了这些,并得出了新的结论.

Tokamak中的Kelvin-Helmholtz不稳定性

<正> Tokamak中的一个重要问题是加热。中性束注入加热是加热的一个有效手段,它使美国PLT上的离子温度达到7.1KeV.但PLT上的中性束注入的不对称性引起等离子体的快速环向旋转,转速可达1×10~7厘米/秒。1979年5月Suckewer等在PLT上测量了速度分布。 在具有速度剪切进行旋转的等离子体中,会不会形成新的磁流体力学不稳定性?1980年

星系螺旋结构三维密度波的流体动力学理论(Ⅰ)——局部解及局部Jeans稳定性的讨论

本文是作者的盘状星系三维密度波工作的第一部分.文中给出了在厚度参数ε及准单色波参数ε_