Direct Numerical Simulation of a Spatially Evolving Supersonic Turbulent Boundary Layer at Ma=6

Direct numerical simulation is carried out for a spatially evolving supersonic turbulent boundary layer at free-stream Mach number 6. To overcome numerical instability, the seventh-order WENO scheme is used for the convection terms of Navier-Stokes equations, and fine mesh is adopted to minimize numerical dissipation. Compressibilty effects on the near-wall turbulent kinetic energy budget are studied. The cross-stream extended self-similarity and scaling exponents including the near-wall region are studied. In high Mach number flows, the coherence vortex structures are arranged to be smoother and streamwised, and the hair-pin vortices are less likely to occur.

Numerical solution of the incompressible Navier-Stokes equations with an upwind compact difference scheme

A new finite difference method for the discretization of the incompressible Navier-Stokes equations is presented. The scheme is constructed on a staggered-mesh grid system. The convection terms are discretized with a fifth-order-accurate upwind compact difference approximation, the viscous terms are discretized with a sixth-order symmetrical compact difference approximation, the continuity equation and the pressure gradient in the momentum equations are discretized with a fourth-order difference approximation on a cell-centered mesh. Time advancement uses a three-stage Runge-Kutta method. The Poisson equation for computing the pressure is solved with preconditioning. Accuracy analysis shows that the new method has high resolving efficiency. Validation of the method by computation of Taylor's vortex array is presented.

DPIV and Interferometry Technique for Velocity Field and Concentration Fields Measurement in Crystal Growth

The property of crystal depends seriously on the solution concentration distribution near the growth surface of a crystal. However, the concentration distributions are affected by the diffusion and convection of the solution. In the present experiment, the two methods of optical measurement are used to obtained velocity field and concentration field of NaClO₃ solution. The convection patterns in sodium chlorate (NaClO₃) crystal growth are measured by Digital Particle image Velocimetry (DPIV) technology. The 2-dimentional velocity distributions in the solution of NaClO3 are obtained from experiments. And concentration field are obtained by a Mach-Zehnder interferometer with a phase shift servo system. Interference patterns were recorded directly by a computer via a CCD camera. The evolution of velocity field and concentration field from dissolution to crystallization are visualized clearly. The structures of velocity fields were compared with that of concentration field.

对流扩散方程的变步长摄动有限差分格式

摄动有限差分(PFD)方法是构造高精度差分格式的一种新方法。变步长摄动有限差分方法是等步长摄动有限差分方法的发展和推广。对需要局部加密网格的计算问题，变步长PFD格式不需要对自变量进行数学变换，且和等步长PFD格式一样，具有如下的共同特点：从变步长一阶迎风格式出发，通过把非微商项(对流系数和源项)作变步长摄动展开，展开幂级数系数通过消去摄动格式修正微分方程的截断误差项求出，由此获得高精度变步长PFD格式。该格式在一、二和三维情况下分别仅使用三、五和七个基点，且具有迎风性。文中利用变步长PFD格式对对流扩散反应模型方程，变系数方程及Burgers方程等进行了数值模拟，并与一阶迎风和二阶中心格式及其问题的精确解作了比较。数值试验表明，与一阶迎风和二阶中心格式相比，变步长PFD格式具有精度高，稳定性与收敛性好的特点。变步长PFD格式与等步长PFD格式相比，变步长PFD解在薄边界层型区域的分辨率得到了明显的提高。

Effects of the orientation of microgravity on the concentration distribution in crystallization from solution

A quasi-steady state growth and dissolution in a 2-D rectangular enclosure is numerically investigated. This paper is an extension to indicate the effects of the orientation of gravity on the concentration field in crystallization from solution under microgravity, especially on the lateral non-uniformity of concentration distribution at the growth surface. The thermal and solute convection are included in this model.

Mixing process in estuaries

A theoretical expression for vertical profile-of horizontal velocity in terms of its depth-average is derived based on oscillatory boundary layer theory and estuarine flow characteristics. The derived theoretical profile is then incorporated into a vertical quasi-two-dimensional model, which is proved advantageous in more physical implications and less CPU time demand. To validate the proposed model, the calculated results are compared to the field data in the Yangtze River Estuary, exhibiting good agreement with observations. The proposed quasi-two-dimensional vertical model is used to study mixing process, especially dependence of salinity distribution and salt front strength on runoff and tides in estuaries.

Multiple Axial Cracks in A Coated Hollow Cylinder due to Thermal Shock

The transient thermal stress problem of an inner-surface-coated hollow cylinder with multiple pre-existing surface cracks contained in the coating is considered. The transient temperature, induced thermal stress, and the crack tip stress intensity factor (SIF) are calculated for the cylinder via finite element method (FEM), which is exposed to convective cooling from the inner surface. As an example, the material pair of a chromium coating and an underlying steel substrate 30CrNi2MoVA is particularly evaluated. Numerical results are obtained for the stress intensity factors as a function of normalized quantities such as time, crack length, convection severity, material constants and crack spacing. (c) 2005 Elsevier Ltd. All rights reserved.

Tailoring of dendritic microstructure in solidification processing by crucible vibration

We analyzed the effects of both natural convection and forced flows on solid–liquid interface morphology during upward Bridgman solidification of metallic alloys. Experiments were carried out on Al–3.5wt% Ni alloy, for a cylindrical sample. The influence of natural convection induced by radial thermal gradient on solidified microstructure was first analyzed as a function of the pulling rate. Then, the influence of axial vibration on solidification microstructure was experimentally investigated by varying vibration parameters (frequency and amplitude). Experimental results demonstrated that vibrations could be used to either attenuate fluid flow in the melt and obtain a uniform dendritic pattern or to promote a fragmented dendritic microstructure. However, no marked effect was observed for cellular growth. This pointed out the critical role of the mushy zone in the interaction between fluid flow and solidification microstructure.

Scour of the seabed under a pipeline in oscillating flow

The scour of the seabed under a pipeline is studied experimentally in this paper. Tests are carried out in a U-shaped oscillatory water tunnel with a box imbedded in the bottom of the test section. By use of the standard sand, clay and plastic grain as the seabed material, the influence of the bed material on the scour is studied. The relationship between the critical initial gap-to-diameter ratio above which no scour occurs and the parameters of the oscillating flow is obtained. The self-burial phenomenon. which occurs for the pipeline not fixed to two sidewalls of the test section, is not observed for the Bred pipeline. The effect of the pipe on sand wave formation is discussed. The maximum equilibrium scour depths For different initial gap-to-diameter ratios, different Kc numbers and different bed sands are also given in this paper.

利用PIV方法测量半浮区液桥热毛细对流的速度场

利用PIV法对硅油液桥热毛细对流的定常速度场进行了实时测量。为了便于测量，液桥上桥端面采取了铜环中嵌透明材料的方法，从液桥的顶部进行观测。当液桥上下桥有温差时，热毛细对流出现；本实验对于不同上下桥的温差，对液桥横剖面内的速度场分别进行了测量，研究外加温差对于流场速度分布的影响；并且在液桥中取了几个典型横截面进行测量，以期对大Pr数液桥的定常速度场有比较全面的定量测量。此外，实验结果也可作为数值模拟计算结果的验证。

Wave-soil-pipe coupling effect on submarine pipeline on-bottom stability

Wave-soil-pipe coupling effect on the untrenched pipeline stability on sands is for the first time investigated experimentally. Tests are conducted in the U-shaped water tunnel, which generates an oscillatory how, simulating the water particle movements with periodically changing direction under the wave action. Characteristic times and phases during the instability process are revealed. Linear relationship between Froude number and non-dimensional pipe weight is obtained. Effects of initial embedment and loading history are observed. Test results between the wavesoil-pipe interaction and pipe-soil interaction under cyclic mechanical loading are compared. The mechanism is briefly discussed. For applying in the practical design, more extensive and systematic investigations are needed.

Effect of liquid bridge volume on the instability in small-Prandtl-number half zones

A perturbation method is used to examine the linear instability of thermocapillary convection in a liquid bridge of floating half-zone filled with a small Prandtl number fluid. The influence of liquid bridge volume on critical Marangoni number and flow features is analyzed. The neutral modes show that the instability is mainly caused by the bulk flow that is driven by the nonuniform thermocapillary forces acting on the free surface. The hydrodynamic instability is dominant in the case of small Prandtl number fluid and the first instability mode is a stationary bifurcation. The azimuthal wave number for the most dangerous mode depends on the liquid bridge volume, and is not always two as in the case of a cylindrical liquid bridge with aspect ratio near 0.6. Its value may be equal to unity when the liquid bridge is relatively slender.

Simulating shock to detonation transition: Algorithm and results

An algorithm based on flux-corrected transport and the Lagrangian finite element method is presented for solving the problem of shock dynamics. It is verified through the model problem of one-dimensional strain elastoplastic shock wave propagation that the algorithm leads to stable, non-oscillatory results. Shock initiation and detonation wave propagation is simulated using the algorithm, and some interesting results are obtained. (C) 1999 Academic Press.

Three-Dimensional Modeling of Heat Transfer and Fluid Flow in Laminar-Plasma Material Re-Melting Processing

Modeling study is performed concerning the heat transfer and fluid flow for a laminar argon plasma jet impinging normally upon a flat workpiece exposed to the ambient air. The diffusion of the air into the plasma jet is handled by using the combined-diffusion-coefficient approach. The heat flux density and jet shear stress distributions at the workpiece surface obtained from the plasma jet modeling are then used to study the re-melting process of a carbon steel workpiece. Besides the heat conduction within the workpiece, the effects of the plasma-jet inlet parameters (temperature and velocity), workpiece moving speed, Marangoni convection, natural convection etc. on the re-melting process are considered. The modeling results demonstrate that the shapes and sizes of the molten pool in the workpiece are influenced appreciably by the plasma-jet inlet parameters, workpiece moving speed and Marangoni convection. The jet shear stress manifests its effect at higher plasma-jet inlet velocities, while the natural convection effect can be ignored. The modeling results of the molten pool sizes agree reasonably with available experimental data.

液层厚度对浮力-热毛细对流面型的影响

将Michelson光学干涉测量系统与图像处理技术相结合,发展形成一种实时诊断热毛细对流和浮力对流流体表面形貌的实验测量系统.采用光学干涉测量方法研究了两端带有温差的矩形池内薄层流体的对流、表面变形、以及表面波的基本问题.应用Fourier变换方法对实验结果进行计算和分析,得到了流体表面变形和表面波的定量的实验结果.实验结果表明了在浮力-热毛细对流的发展过程中,首先出现流体的表面变形,之后在该变形的基础上,叠加了一个表面波的信息,该表面变形和表面波与流体的温度梯度、表面张力、以及浮力有直接的关系;表面波隐藏在表面变形内.