Marine Mammal Strandings in the United States: Proceedings of the Second Marine Mammal Stranding Workshop, Miami, Florida, December 3-5, 1987

Over a decade ago, in August 1977, the First Marine Mammal Stranding Workshop was convened in Athens, Georgia. That workshop, organized by j.R. Geraci and D.J. St. Aubin, not only considered biology and pathology of stranded marine mammals, but it also served as a springboard for the formation of regional marine mammal stranding networks in the United States. The ramifications have been extremely important to the field of marine mammalogy since, for some species, examination or rehabilitation of stranded specimens serves as virtually the only source of information on distribution, anatomy, physiology, reproduction, and pathology. The First Marine Mammal Stranding Workshop led to increased awareness of the marine mammals themselves, as well as the logistic and legal factors associated with effective handling of the animals. A number of individuals indicated that they felt that a Second Marine Mammal Stranding Workshop held prior to the Seventh Biennial Conference on the Biology of Marine Mammals (Miami, Florida; December 1987) would be both timely and productive. Accordingly, we organized the workshop and scheduled it to occur on 3-5 December. Our goals for the workshop were several, including 1) providing descriptions of some research, especially new techniques, regarding stranded marine mammals; 2) providing a forum where scientists could interact and possibly initiate cooperative research activities; 3) presenting information regarding procedures used effectively to handle stranded animals; 4) assessing ways to standardize data and specimen collection, archiving, and retrieval; and 5) providing a forum for assessing accomplishments and status of regional stranding networks to date, as well as for making recommendations regarding future activities of the networks. Nearly 100 individuals representing Federal and State governments, academic institutions, the oceanarium industry, consulting groups, conservation organizations, and the private sector attended the workshop (see Workshop Participants, this volume). (PDF file contains 166 pages.)

可压平面混合层稳定性分析及数值模拟

可压平面混合层是包含复杂多时空尺度运动的非定常流体力学部问题，具有深刻的理论意义和广泛的应用背景。针对该问题所涉及内容的多面性，本文的目的是，基于高精度、高分辨率数值算法的构造、发展和数值行为分析，采用线性稳定性分析和直接数值模拟方法。从理论和计算两方面集中研究压缩性效应、粘性效应、初值效应以及燃烧反应放热效应等对可压平面混合层早期稳定性行为和大尺度拟序涡结构非线性演化的影响。以混合层已有研究成果的分析和综述为开端，论文主体共包括四部分：第一部分是可压平面混合层时间/空间模式数值线性稳定性分析。实现了高精度对称紧致差分格式（SCD）对可压粘性扰动线性稳定性边值问题的求解，对导出的线性和非线性离散特征值问题，提出了两个高效局部解法。研究涉及二维/三维扰动波、无粘/粘性扰动波、特征函数和特征值谱、第一/第二模态、超声速快/慢模态、速度比和密度比等。验证了对流Mach数Mc为一个合理的压缩性参数。指出压缩性效应和粘性效应对最不稳定扰动波的波数（频率）和增长率呈相拟的抑制作用，且时间模式稳定性分析结果在许多方面是可信的。从随机和线性扰动场出发，采用高精度五阶迎风紧致和六阶对称紧致混合差分算法（UCD5/SCD6）对可压平面混合层的稳定性特征进行了直接数值模拟，揭示了初始主导线性扰动与一些实际涡结构非线性作用形态间的内在关联，印证了线性稳定性分析方法的合理性和有效性。第二部分是高精度迎风紧致差分格式（UCD）时空全离散数值行为分析。导出了其一维/二维一般色散表达式。研究表明，UCD格式在高波数区具有内在的全离散耗散和色散特性；其数值群速度的快/慢特征可因CFL数不同而改变；在稳定CFL数下简单附加人工粘性可强化UCD格式在高波数区的耗散量；提高时间精度可放宽稳定CFL数限制；UCD格式的二维全离散色散介质中存在三个不同性质的数值波，其全离散稳定性由数值声波主控。第三部分实现了高精度UCD5/SCD6差分算法对空间发展可压平面混合层的直接数值模拟。通过亚谐扰动波的个数和扰动频率的控制，捕捉到了基频涡的饱和、一次和二次对并等现象，显示了大尺度涡结构与入中初始扰动方式之间的内在联系。利用参数Mc观察了压缩性效应对大尺度涡空间演化及其相互作用的影响。第四部分实现了高精度UCD5/SCD6差分算法对非预混扩散火焰化学反应平面混合层的直接数值模拟。研究指出，放热效应可抑制和延迟涡的形成，使基频涡卷拉伸甚至丧失，混合层Reynolds 应力ρu'v'和流向速度波动关联项u'v'下降，以致涡结构与外流动量交换和标量输运减少，脉动输运能力被削弱，从而混合效率、产物生成率和混合层增长率下降，放热主要通过膨胀效应和斜压效应来抑制大尺度涡的演化。

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.

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4:2244-2251; Int. J. Multiphase Flow 2000; 26:1583-1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes.

Generalization of the JTZ model to open plane wakes

The JTZ model [C. Jung, T. T¶el and E. Ziemniak, Chaos 3, (1993) 555], as a theoretical model of a plane wake behind a circular cylinder in a narrow channel at a moderate Reynolds number, has previously been employed to analyze phenomena of chaotic scattering. It is ex- tended here to describe an open plane wake without the con¯ned nar- row channel by incorporating a double row of shedding vortices into the intermediate and far wake. The extended JTZ model is found in qualitative agreement with both direct numerical simulations and ex- perimental results in describing streamlines and vorticity contours. To further validate its applications to particle transport processes, the in- teraction between small spherical particles and vortices in an extended JTZ model °ow is studied. It is shown that the particle size has signif- icant in°uences on the features of particle trajectories, which have two characteristic patterns: one is rotating around the vortex centers and the other accumulating in the exterior of vortices. Numerical results based on the extended JTZ model are found in qualitative agreement with experimental ones in the normal range of particle sizes.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

AUTOROTATION OF TWO TANDEM TRIANGULAR CYLINDERS

The autorotation of two tandem triangular cylinders at different gap distances is investigated by numerical simulations. At the Reynolds number of 200， three distinct regimes are observed with the increase of gap distance: namely, angular oscillation, quasi-periodic autorotation and ‘chaotic’ autorotation. For various gap distances, the characteristic of vortex shedding and vortex interaction are discussed. The phase graphs (angular acceleration vs. angular velocity) and the power spectra of moment are analyzed to characterize the motion of the cylinder. The Lyapunov exponent is also calculated to identify the existence of chaos.

聚二甲基硅氧烷（PDMS）表面液滴动力学的实验研究

液滴是自然界中普遍存在的一种物质形态。非连续微流体（液滴）是近年来微流体技术重要发展方向之一。对液滴的产生、启动、移动、合并、分离和碰撞过程的研究对于航天、微纳系统、电子显示、计算机冷却、喷墨、生物医学等学科领域有着重要的应用价值。液滴属于软物质，其力学性质介于流体和固体之间，其类固体（solid-like）行为来自于曲率产生的Laplace压力和表面张力的约束。对液滴动力学行为的研究有着重要的学术价值。 本文的主要工作是针对生物微电子机械系统（Bio-MEMS）以及柔性微纳电子加工中常用的材料聚二甲基硅氧烷（Polydimethylsiloxane，PDMS）为基底的液滴动力学实验研究。 液滴是一个理想的微反应器，许多实验可以集成在一个液滴或多个液滴内完成。液滴本身的动力学特性对于实验的完成效率和质量有着重要的影响。液滴的微操控技术包括多相流法、电润湿法、热毛细管法、介电泳法等。液滴的动力学特性受到基底的影响非常大，包括基频、振动模态、运动过程等均随基底的润湿性、弹性模量的变化而有所变化。 在Bio-MEMS以及柔性微纳电子加工当中，PDMS扮演着越来越重要的角色，尤其是PDMS的润湿性和电润湿特性。目前的PDMS在Bio-MEMS当中主要是用于制备各种微流道。常见的问题主要是一方面PDMS是疏水材料，影响流体的输运。另一方面是液体在这种低Reynolds数情况下不易混合，反应效率低。本文提出了在PDMS表面溅射纳米厚度的金来减小PDMS表观接触角的方法。这种方式在特定喷金量的情况下可以在PDMS表面产生多层次的压应力波纹。这种压应力波纹对于柔性微纳电子加工，以及微流道中加速流体混合有着非常重要的作用。 电润湿是另一种可以使PDMS亲水化的方法。实验证明，PDMS具有较好的电润湿性质。此外电润湿也是目前操纵液滴的主要方式。目前一个常见的问题是电击穿现象阻碍了驱动电压的低压化，且低Reynolds数情况下液滴的混合效率偏低。此外电极还会由于少量电解的发生导致腐蚀及对液体样品的污染。本文提出了接触式的电润湿，在电极逐渐触碰液滴的过程中，液滴发生百Hz的失稳振动，稳定后接触角减小。这种电润湿模式可以有效的提高临界击穿电压，避免液滴被腐蚀后的电极污染，同时可以加快液滴的混合效率。其失稳特征时间在10 ms量级，这恰是所用液滴特征尺度在1 mm左右的电润湿器件的最快响应时间。并采用液滴振动的理论估算了液滴的失稳时间，同时还考虑了基底润湿性对液滴振动过程的影响。 液滴的启动是电润湿操控液滴过程中的重要环节。通常的液滴启动都是在非连续基底上依靠逻辑电路产生的电势变化来驱动液滴。无论是逻辑电路的设计还是驱动装置的加工都非常复杂。本文首次实现了在超疏水生物样品荷叶上的液滴启动，启动速度为数十毫米/秒，启动时间为10 ms量级。并利用PDMS成功的仿制了荷叶结构实现了超疏水的PDMS表面，荷叶同仿荷叶的PDMS超疏水表面具有相近的润湿性。 在数字微流体操控液滴的过程中，液滴的合并涉及液滴的碰撞，而且MEMS系统当中利用液滴撞击进行冷却的实验已经有所开展。同时理解液滴碰撞还对许多领域包括生物、化学、喷墨、大气物理等有着非常重要的作用。本文实验研究了Weber数和毛细数对液滴碰撞过程的影响，通过改变Weber数和毛细数得到了四种不同的响应模式。

Numerical Investigation On The Flowfield Of "Swallowtail" Cavity For Supersonic Mixing Enhancement

A "swallowtail" cavity for the supersonic combustor was proposed to serve as an efficient flame holder for scramjets by enhancing the mass exchange between the cavity and the main flow. A numerical study on the "swallowtail" cavity was conducted by solving the three-dimensional Reynolds-averaged Navier-Stokes equations implemented with a k-epsilon turbulence model in a multi-block mesh. Turbulence model and numerical algorithms were validated first, and then test cases were calculated to investigate into the mechanism of cavity flows. Numerical results demonstrated that the certain mass in the supersonic main flow was sucked into the cavity and moved spirally toward the combustor walls. After that, the flow went out of the cavity at its lateral end, and finally was efficiently mixed with the main flow. The comparison between the "swallowtail" cavity and the conventional one showed that the mass exchanged between the cavity and the main flow was enhanced by the lateral flow that was induced due to the pressure gradient inside the cavity and was driven by the three-dimensional vortex ring generated from the "swallowtail" cavity structure.

Experimental Study Of Multi-Hole Cooling For Integrally-Woven, Ceramic Matrix Composite Walls For Gas Turbine Applications

In this paper, multi-hole cooling is studied for an oxide/oxide ceramic specimen with normal injection holes and for a SiC/SiC ceramic specimen with oblique injection holes. A special purpose heat transfer tunnel was designed and built, which can provide a wide range of Reynolds numbers (10(5)similar to 10(7)) and a large temperature ratio of the primary flow to the coolant (up to 2.5). Cooling effectiveness determined by the measured surface temperature for the two types of ceramic specimens is investigated. It is found that the multi-hole cooling system for both specimens has a high cooling efficiency and it is higher for the SiC/SiC specimen than for the oxide/oxide specimen. Effects on the cooling effectiveness of parameters including blowing ratio, Reynolds number and temperature ratio, are studied. In addition, profiles of the mean velocity and temperature above the cooling surface are measured to provide further understanding of the cooling process. Duplication of the key parameters for multi-hole cooling, for a representative combustor flow condition (without radiation effects), is achieved with parameter scaling and the results show the high efficiency of multi-hole cooling for the oblique hole, SiC/SiC specimen. (C) 2008 Elsevier Ltd. All rights reserved.

A CE/SE Scheme for Flows in Porous Media and Its Applications

In this paper, a new computational scheme for solving flows in porous media was proposed. The scheme was based on an improved CE/SE method (the space-time Conservation Element and Solution Element method). We described porous flows by adopting DFB (Brinkman-Forchheimer extended Darcy) equation. The comparison between our computational results and Ghia's confirmed the high accuracy, resolution, and efficiency of our CE/SE scheme. The proposed first-order CE/SE scheme is a new reliable way for numerical simulations of flows in porous media. After investigation of effects of Darcy number on porous flow, it shows that Darcy number has dominant influence on porous flow for the Reynolds number and porosity considered.

Heat Transfer of Aviation Kerosene at Supercritical Conditions

The heat transfer characteristics of China no. 3 kerosene were investigated experimentally and analytically under conditions relevant to a regenerative cooling system for scramjet applications. A test facility developed for the present study can handle kerosene in a temperature range of 300-1000 K, a pressure range of 2.6-5 MPa, and a mass How rate range of 10-100 g/s. In addition, the test section was uniquely designed such that both the wall temperature and the bulk fuel temperature were measured at the same location along the flowpath. The measured temperature distributions were then used to analytically deduce the local heat transfer characteristics. A 10-component kerosene surrogate was proposed and employed to calculate the fuel thermodynamic and transport properties that were required in the heat transfer analysis. Results revealed drastic changes in the fuel flow properties and heat transfer characteristics when kerosene approached its critical state. Convective heat transfer enhancement was also found as kerosene became supercritical. The heat transfer correlation in the relatively low-fuel-temperature region yielded a similar result to other commonly used jet fuels, such as JP-7 and JP-8, at compressed liquid states. In the high-fuel-temperature region, near and beyond the critical temperature, heat transfer enhancement was observed; hence, the associated correlation showed a more significant Reynolds number dependency.

Simulation of gas flow over a micro cylinder

Gas flow over a micro cylinder is simulated using both a compressible Navier-Stokes solver and a hybrid continuum /particle approach. The micro cylinder flow has low Reynolds number because of the small length scale and the low speed, which also indicates that the rarefied gas effect exists in the flow. A cylinder having a diameter of 20 microns is simulated under several flow conditions where the Reynolds number ranges from 2 to 50 and the Mach number varies from 0.1 to 0.8. It is found that the low Reynolds number flow can be compressible even when the Mach number is less than 0.3, and the drag coefficient of the cylinder increases when the Reynolds number decreases. The compressible effect will increase the pressure drag coefficient although the friction coefficient remains nearly unchanged. The rarefied gas effect will reduce both the friction and pressure drag coefficients, and the vortex in the flow may be shrunk or even disappear.

IP方法模拟硬盘读写头气膜润滑问题

本论文中采用IP方法模拟硬盘读写头的气膜润滑问题，并选择耗时但精准的DSMC方法作为该问题中IP方法的检验标准，IP方法计算得到的读写头表面压力分布及读写头所受净升力为读写头的设计提供可靠的参考数据。 首先，将读写头复杂表面简化为平面并假定尺寸为微米（远小于真实的毫米），微米尺寸时可以得到DSMC的模拟结果用于比较。IP法对该问题的模拟结果与DSMC的模拟结果完全吻合。 其次，选定一个真实的读写头表面几何，但仍然假定尺寸为微米。IP法在该问题的应用过程中遇到四点困难，我们提出或采取了相应的解决办法。在复杂表面几何的读写头气膜润滑问题中，IP法的模拟结果与DSMC的模拟结果也相符一致。 最后完成真实几何、真实尺寸（毫米）的读写头气膜润滑问题的模拟，这也是本论文的最终研究目标。IP法模拟该问题时需要采用大网格，为此我们构造二维检验模型验证了大网格的合理性。 真实几何读写头气膜润滑问题的模拟结果表明，在流场的局部区域压力沿着高度方向变化显著，这意味着读写头问题的传统计算方法（即概括性雷诺方程）在该局部区域不成立，故用它计算得到的结果值得认真检验，但文献中尚未给出过这类检验。这正是该论文的研究背景和采用新方法（IP方法）的原因。 本论文研究的问题属于过渡领域、低速流动问题，故本论文还介绍了该领域另外一种粒子模拟方法（LVDSMC方法）以及我们对该方法提出的一些改进。 关键词：硬盘读写头，气膜润滑，概括性雷诺方程，IP方法，DSMC方法，过渡领域，LVDSMC方法

A Continuation Method Applied to the Study of Thermocapillary Instabilities in Liquid Bridges

A continuation method is applied to investigate the linear stability of the steady, axisymmetric thermocapillary flows in liquid bridges. The method is based upon an appropriate extended system of perturbation equations depending on the nature of transition of the basic flow. The dependence of the critical Reynolds number and corresponding azimuthal wavenumber on serval parameters is presented for both cylindrical and non-cylindrical liquid bridges.