高雷诺数流动的控制方程体系和扩散抛物化Navier-Stokes方程组的意义和用途

在计算机发达的时代,高雷诺(Re)数绕流计算中有无必要使用简化NS方程组,本文讨论这个问题.主要内容如下:(1)高Re数绕流包含3种基本流动:所有方向对流占优流动、所有方向对流扩散竞争流动和部分方向对流占优部分方向对流扩散竞争流动(简称干扰剪切流动),3个基本流动的特征彼此不同且在流场中所占领域大小彼此相差悬殊,NS方程区域很小,它们的最简单控制方程组Euler、Navier-Stokes(NS)和扩散抛物化(DP)NS方程组的数学性质彼此不同,因此利用Euler-DPNS-NS方程组体系分析计算高Re数绕流流动就是一个合乎逻辑的选择,该法与利用单一NS方程组的常用方法可以彼此检验和补充.(2)流体之间以及流体与外界的动量、能量和质量交换,流态从层流到湍流的演化主要发生在干扰剪切流动中,干扰剪切流及其最简单控制方程--DPNS方程组具有基础意义;DPNS方程组笔者在1967年已提出.(3)诸简化NS方程组:DPNS、抛物化(P)NS、薄层(TL)NS、黏性层(VL)NS方程组的发展、相互关系,它们的历史贡献和今后的用途;它们的数学性质均为扩散抛物型,但它们包含的黏性项彼此有所不同;从流体力学角度来看,它们中只有DPNS方程组能够准确描述干扰剪切流动.提出把诸简化NS方程组统一为DPNS方程组的建议.(4)干扰剪切流--DPNS方程组与无干扰剪切流--边界层方程组之间的关系以及进一步研究干扰剪切流的意义.

Simulations of oil spread on a water surface

Numerical simulation of an oil slick spreading on still and wavy surfaces is described in this paper. The so-called sigma transformation is used to transform the time-varying physical domain into a fixed calculation domain for the water wave motions and, at the same time, the continuity equation is changed into an advection equation of wave elevation. This evolution equation is discretized by the forward time and central space scheme, and the momentum equations by the projection method. A damping zone is set up in front of the outlet boundary coupled with a Sommerfeld-Orlanski condition at that boundary to minimize the wave reflection. The equations for the oil slick are depth-averaged and coupled with the water motions when solving numerically. As examples, sinusoidal and solitary water waves, the oil spread on a smooth plane and on still and wavy water surfaces are calculated to examine the accuracy of simulating water waves by Navier-Stokes equations, the effect of damping zone on wave reflection and the precise structures of oil spread on waves.

Diffusive Wave Solutions for Open Channel Flows with Uniform and Concentrated Lateral Inflow

The diffusive wave equation with inhomogeneous terms representing hydraulics with uniform or concentrated lateral inflow intoa river is theoretically investigated in the current paper. All the solutions have been systematically expressed in a unified form interms of response function or so called K-function. The integration of K-function obtained by using Laplace transform becomesS-function, which is examined in detail to improve the understanding of flood routing characters. The backwater effects usuallyresulting in the discharge reductions and water surface elevations upstream due to both the downstream boundary and lateral infloware analyzed. With a pulse discharge in upstream boundary inflow, downstream boundary outflow and lateral inflow respectively,hydrographs of a channel are routed by using the S-functions. Moreover, the comparisons of hydrographs in infinite, semi-infiniteand finite channels are pursued to exhibit the different backwater effects due to a concentrated lateral inflow for various channeltypes.

The scattering of general SH plane wave by interface crack between two dissimilar viscoelastic bodies

The scattering of general SH plane wave by an interface crack between two dissimilar viscoelastic bodies is studied and the dynamic stress,intensity factor at the crack-tip is computed. The scattering problem can be decomposed into two problems: one is the reflection and refraction problem of general SH plane waves at perfect interface (with no crack); another is the scattering problem due to the existence of crack. For the first problem, the viscoelastic wave equation, displacement and stress continuity conditions across the interface are used to obtain the shear stress distribution at the interface. For the second problem, the integral transformation method is used to reduce the scattering problem into dual integral equations. Then, the dual integral equations are transformed into the Cauchy singular integral equation of first kind by introduction of the crack dislocation density function. Finally, the singular integral equation is solved by Kurtz's piecewise continuous function method. As a consequence, the crack opening displacement and dynamic stress intensity factor are obtained. At the end of the paper, a numerical example is given. The effects of incident angle, incident frequency and viscoelastic material parameters are analyzed. It is found that there is a frequency region for viscoelastic material within which the viscoelastic effects cannot be ignored.

Propagation Properties of Electromagnetic Wave in a Plasma Slab

In this paper, the calculated results about the propagation properties of electromagnetic wave in a plasma slab are described. The relationship of the propagation properties with frequencies of electromagnetic wave, and parameters of plasma (electron temperature, electron density, dimensionless collision frequency and the size of the plasma slab) is analyzed.

Hydrothermal Wave in a Shallow Liquid Layer

The oscillatory thermocapillary convection and hydrothermal wave in a shallow liquid layer, where a temperature difference is applied between two parallel sidewalls, have been numerically investigated in a two-dimensional model. The oscillatory thermocapillary convection and hydrothermal wave appear if the Marangoni number is larger than a critical value. The critical phase speed and critical wave number of the hydrothermal wave agree with the ones given analytically by Smith and Davis in the microgravity environment, and it travels in the direction opposed to the surface flow. Another wave traveled downstream in addition to the hydrothermal wave traveled upstream was observed in the case of earth gravity condition.

A New Version of Smoothed Point Method to Simulate Linear Elastic Wave Propagation

By using the kernel function of the smoothed particle hydrodynamics (SPH) and modification of statistical volumes of the boundary points and their kernel functions, a new version of smoothed point method is established for simulating elastic waves in solid. With the simplicity of SPH kept, the method is easy to handle stress boundary conditions, especially for the transmitting boundary condition. A result improving by de-convolution is also proposed to achieve high accuracy under a relatively large smooth length. A numerical example is given and compared favorably with the analytical solution.

The Optical Methods on Measuring Surface Deformation and Surface Wave in The Thermal Capillary Convection

An optical diagnostic system consisting of the Michelson interferometer with the image processor has been developed for the study of the kinetics of the thermal capillary convection. The capillary convection, surface deformation, surface wave and the velocity field in a rectangular cavity with different temperature's sidewalls have been investigated by optical interference method and PIV technique. In order to calculate the surface deformation from the interference fringe, Fourier transformation is used to grating analysis. The quantitative results of the surface deformation and surface wave have been calculated from the interference fringe pattern.

Influence of insoluble surfactant on the deformation and breakup of a bubble or thread in a viscous fluid

The influence of surfactant on the breakup of a prestretched bubble in a quiescent viscous surrounding is studied by a combination of direct numerical simulation and the solution of a long-wave asymptotic model. The direct numerical simulations describe the evolution toward breakup of an inviscid bubble, while the effects of small but non-zero interior viscosity are readily included in the long-wave model for a fluid thread in the Stokes flow limit. The direct numerical simulations use a specific but realizable and representative initial bubble shape to compare the evolution toward breakup of a clean or surfactant-free bubble and a bubble that is coated with insoluble surfactant. A distinguishing feature of the evolution in the presence of surfactant is the interruption of bubble breakup by formation of a slender quasi-steady thread of the interior fluid. This forms because the decrease in surface area causes a decrease in the surface tension and capillary pressure, until at a small but non-zero radius, equilibrium occurs between the capillary pressure and interior fluid pressure. The long-wave asymptotic model, for a thread with periodic boundary conditions, explains the principal mechanism of the slender thread's formation and confirms, for example, the relatively minor role played by the Marangoni stress. The large-time evolution of the slender thread and the precise location of its breakup are, however, influenced by effects such as the Marangoni stress and surface diffusion of surfactant. © 2008 Cambridge University Press.

Influence of surfactant solubility on the deformation and breakup of a bubble or capillary jet in a viscous fluid

In a previous study [M. Hameed, J. Fluid Mech. 594, 307 (2008)] the authors investigated the influence of insoluble surfactant on the evolution of a stretched, inviscid bubble surrounded by a viscous fluid via direct numerical simulation of the Navier-Stokes equations, and showed that the presence of surfactant can cause the bubble to contract and form a quasisteady slender thread connecting parent bubbles, instead of proceeding directly toward pinch-off as occurs for a surfactant-free bubble. Insoluble surfactant significantly retards pinch-off and the thread is stabilized by a balance between internal pressure and reduced capillary pressure due to a high concentration of surfactant that develops during the initial stage of contraction. In the present study we investigate the influence of surfactant solubility on thread formation. The adsorption-desorption kinetics for solubility is in the diffusion controlled regime. A long-wave model for the evolution of a capillary jet is also studied in the Stokes flow limit, and shows dynamics that are similar to those of the evolving bubble. With soluble surfactant, depending on parameter values, a slender thread forms but can pinch-off later due to exchange of surfactant between the interface and exterior bulk flow. © 2009 American Institute of Physics.

Shock wave diffraction and reflection around a dusty square cavity

The diffraction and reflection of planar shock wave around a dusty square cavity is investigated numerically, which is embedded in the net bottom surface of a two-dimensional channel, and the induced gas-particle two-phase now. The wave patterns at different times are obtained for three different values of the particle diameter. The computational results show that the existence of particles affects appreciably the shock wave diffraction and cavity flow.

Variation of drag coefficient with wind wave development

Turbulent air flows over developing wind waves in the air-sea boundary layer are numerically simulated without considering wave breaking. Influences of wind waves on air flows are considered using a model of significant wave and surface roughness, with a formula proposed for calculating the surface roughness, k - epsilon model is adopted to simulate turbulent flows. The results of the drag coefficient and turbulence characteristics agree well with the observations.

Numerical Investigation of Toroidal Shock Wave Focusing in a Cylindrical Chamber

In this paper, focusing of a toroidal shock wave propagating from an annular shock tube into a cylindrical chamber was investigated numerically with the dispersion controlled dissipation (DCD) scheme. The first case for an incident Mach number of 1.5 was conducted and compared with experiments for validation. Then, several cases were calculated for higher incident Mach numbers varying from 2.0 to 5.0, and complicated flow structures were observed. The numerical study was mainly focused on two aspects: focusing process and flow structures. The process, including diffraction, focusing, and reflection, is displayed to reveal the focusing mechanism, and the flow structures at different incident. Mach numbers are used to demonstrate shock reflection styles and focusing characteristics.

Improved Analysis Method for Wave-Induced Pipeline Stability on Sandy Seabed

The existing Det Norske Veritas DNV Recommended Practice RP E305 for pipeline on-bottom stability is mainly based on the pipe–soil interaction model reported by Wagner et al. in 1987, and the wake model reported by Lambrakos et al. in 1987, to calculate the soil resistance and the hydrodynamic forces upon pipeline, respectively. Unlike the methods in the DNV Practice, in this paper, an improved analysis method is proposed for the on-bottom stability of a submarine pipeline, which is based on the relationships between Um/ gD 0.5 and Ws / D2 for various restraint conditions obtained by the hydrodynamic loading experiments, taking into account the coupling effects between wave, pipeline, and sandy seabed. The analysis procedure is illustrated with a detailed flow chart. A comparison is made between the submerged weights of pipeline predicted with the DNV Practice and those with the new method. The proposed analysis method may provide a helpful tool for the engineering practice of pipeline on-bottom stability design.

GSPS photonic analogue to digital conversion system using broadband continuous wave source

An 80 GSPS photonic ADC system is demonstrated, using broadband MLL and dispersive fibre to form a continuous waveform with time-wavelength mapping, and AWG to channelise. Tests are carried out for RF signals up to 10GHz. © 2005 Optical Society of America.