利用Ghost Fluid方法模拟激波与柱形界面相互作用

利用Ghost Fluid方法(后面简称Ghost方法)和γ-model方法，在同样的时空离散精度条件下，对激波与柱形界面相互作用的二维可压缩流场进行了直接模拟，并与实验结果相比较．从模拟结果看，在短时间内， Ghost方法和γ-model方法模拟的结果与实验结果基本相同，两种方法均正确地模拟出界面的位置、激波的强度和速度．但随着时间的发展，具有较大数值耗散的γ-model方法的计算结果与实验差别越来越大；而数值耗散较小的Ghost方法能较为正确地模拟界面的运动．

Two Bubbles Thermocapillary Migration Experimental System

An experimental investigation will be performed on the thermocapillary motion of two bubbles in Chinese return-satellite. The experiment will study the migration process of bubble caused by thermocapillary effect in microgravity environment, and their interaction between two bubbles. The bubble is driven by the thermocapillary stress on the surface on account on the variation of the surface tension with temperature. The interaction between two bubbles becomes significant as the separation distance between them is reduced drastically so that the bubble interaction has to be considered. Recently, the problem has been discussed on the method of successive reflections, and accurate migration velocities of two arbitrarily oriented bubbles were derived for the limit of small Marangoni and Reynolds numbers. Numerical results for the migration of the two bubbles show that the interaction between two bubbles has significant influence on their thermocapillary migration velocities with a bubble approaching another. However, there is a lack of experimental validate for the theoretic results. Now the experimental facility is designed for experimenting time after time. A cone-shaped top cover is used to expel bubble from the cell after experiment. But, the cone-shaped top cover can cause temperature uniformity on horizontal plane in whole cell. Therefore, a metal board with multi-holes is fixed under the top cover. The board is able to let the temperature distribution on the board uniform because of their high heat conductivity, and the bubble can pass through it. In the system two bubbles are injected into the test cell respectively by two sets of cylinder. And the bubbles sizes are controlled by two sets of step-by-step motor. It is very important problem that bubble can be divorced from the injecting mouth in microgravity environment. Thus, other two sets of device for injecting mother liquid were used to push bubble. The working principle of injecting mother liquid is to utilize pressure difference directly between test cell and reservoir

Computation of flow over a rotating body on unstructured Chimera mesh

Flow around moving boundary is ubiquitous in engineering applications. To increse the efficienly of the algorithm to handle moving boundaries is still a major challenge in Computational Fluid Dynamics (CFD). The Chimera grid method is one type of method to handle moving boundaries. A concept of domain de-composition has been proposed in this paper. In this method, sub-domains are meshed independently and governing equations are also solved separately on them. The Chimera grid method was originally used only on structured (curvilinear) meshes. However, in a problem which involves both moving boundary and complex geometry, the number of sub-domains required in a traditional (structured) Chimera method becomes fairly large. Thus the time required in the interior boundary locating, link-building and data exchanging also increases. The use of unstructured Chimera grid can reduce the time consumption significantly by the reduction of domain(block) number. Generally speaking, unstructured Chimera grid method has not been developed. In this paper, a well-known pressure correction scheme - SIMPLEC is modified and implemented on unstructured Chimera mesh. A new interpolation scheme regarding the pressure correction is proposed to prevent the possible decoupling of pressure. A moving-mesh finite volume approach is implemented in an inertial reference frame. This approach is then used to compute incompressible flow around a rotating circular and elliptic cylinder. These numerical examples demonstrate the capability of the proposed scheme in handling moving boundaries. The numerical results are in good agreement with other experimental and computational data in literature. The method proposed in this paper can be efficiently applied to more challenge cases such as free-falling objects or heavy particles in fluid.

Adiabatic shear localization: occurrence, theories, and applications

谈镐生文集

本文集收入了著名力学家、应用数学家、中国科学院院士谈镐生先生在流体力学、稀薄气体动力学和应用数学研究领域的论文和研究报告26篇，谈镐生先生倡导和支持力学基础研究的有关文章和论述21篇，谈先生的学术活动和生活图片多幅，以及谈镐生先生生平，最后附有谈先生生平年表。

 

 

A stationary visual census technique for quantitatively assessing community structure of coral reef fishes

A new method is described and evaluated for visually sampling reef fish community structure in environments with highly diverse and abundant reef fish populations. The method is based on censuses of reef fishes taken within a cylinder of 7.5 m radius by a diver at randomly selected, stationary points. The method provides quantitative data on frequency of occnrrence, fish length, abundance, and community composition, and is simple, fast, objective, and repeatable. Species are accumulated rapidly for listing purposes, and large numbers of samples are easily obtained for statistical treatment. The method provides an alternative to traditional visual sampling methods. Observations showed that there were no significant differences in total numbers of species or individuals censused when visibility ranged between 8 and 30 m. The reefs and habitats sampled were significant sources of variation in number of species and individuals censused, but the diver was not a significant influence. Community similarity indices were influenced significantly by the specific sampling site and the reef sampled, but were not significantly affected by the habitat or diver (PDF file contains 21 pages.)

再入尾迹电磁特性的湍流效应研究

本文旨意在于通过探讨高超声速再入尾迹中的湍流等离子体与电磁波相互作用的机理，以及建立能反映此机理的应用性理论模型，从而提供一套可进行目标特性分析的方法，以便为工程部门的突防技术服务。本题目在再入气动物理现象研究中具有重要意义。综合分析指出，地面雷达观测到的非相干散射信号主要来源于再入尾迹的亚密湍流区产生的体积散射。因此，电磁散射特性分析主要针对尾迹亚密湍流等离子体。并且，这里所有的分析都是根据在工程应用中最成熟的一阶畸变波Born近似理论模型。再入尾迹电磁特性的湍流效应研究，着眼点就在于湍流等离子体场的研究。对湍流等离子体场理论模型，本文试图通过模式理论来表达，即求解平均化的全Navier-Stokes方程及其封闭方程k-ε-g模型，从而准确获得流动平均场和脉动场信息。这种表达方式较以前有了较大改进。注意到高超声速流动具有强烈可压缩性的特点，故使用的N-S平均方程由质量加权平均过程产生，湍流模型方程也经过可压缩性修正。方程的离散求解方法，都是运用带矢通量分裂的二阶TVD格式的有限体积法。再入尾迹湍流场的初始条件由近尾迹（底部）流动经N-S方程求解给定，初始值更加准确可靠。尾迹从层流到湍流的转捩过程采用相对成熟的半经验公式确定。飞行器的高超声速再入过程必然导致它周围的空气温度升高，使得流动表现出真实气体效应。对重点考察的湍流流动而言，真实气体效应主要表现为气体处于热化学平衡状态。就工程部门面临的实际问题，把一阶畸变波Born近似的解算方法做些改进，使其能够处理的范围从轴对称尾迹扩展到三维湍流等离子体场是必要的。这为深入的理论分析提供了有力的保障。在能够准确模拟湍流流动的刻划雷达散射截面的基础上，考察亚密湍流等离子体对电磁散射的影响。通过选择的几个有代表性的因素进行讨论，初步结果表明：湍流转捩方式、湍流尺度对尾迹雷达散射截面值计算影响不大，而电子组份脉动能初始值影响较明显，且在特定条件下湍流模型的影响亦不大。但由于湍流模型涉及脉动初始值，其影响需进一步确定。同时，一些今后开展继续此项研究工作的有益建议也提了出来。

均匀来流绕旋转及旋转振荡圆柱绕流研究

该文通过数值方法求解二维不可压Navier-Stokes方程，对均匀来流中静止、旋转和旋转振荡圆柱绕流进行了系统的数值模拟.该文采用有限体积法对控制方程进行离散，选用元结构化四边形网格剖分计算区域，关于速度-压力耦合的处理使用了SIMPLEC方法.经过了大量的数值模拟，分盺166L鸬玫搅苏饧钢秩屏鞯氖的Ｄ饨峁?该文重点是用快速傅里叶变换（FFT）方法对旋转振荡圆柱绕流中的频率耦合现象进行研究，并分析在不同频率耦合作用下涡形成、发展和脱落的规律.

A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.

Effects of geometric shape on the hydrodynamics of a self-propelled flapping foil

The hydrodynamics of a free flapping foil is studied numerically. The foil undergoes a forced vertical oscillation and is free to move horizontally. The effect of chord-thickness ratio is investigated by varying this parameter while fixing other ones such as the Reynolds number, the density ratio, and the flapping amplitude. Three different flow regimes have been identified when we increase the chord-thickness ratio, i.e., left-right symmetry, back-and-forth chaotic motion, and unidirectional motion with staggered vortex street. It is observed that the chord-thickness ratio can affect the symmetry-breaking bifurcation, the arrangement of vortices in the wake, and the terminal velocity of the foil. The similarity in the symmetry-breaking bifurcation of the present problem to that of a flapping body under constraint is discussed. A comparison between the dynamic behaviors of an elliptic foil and a rectangular foil at various chord-thickness ratios is also presented.