水下机器人粘性类水动力数值计算方法研究

通过对通用流体动力学仿真软件CFX的研究,提出了一套水下机器人粘性类水动力的数值计算方法.该方法采用标准k-ε湍流模型计算位置力系数,采用标准k-ω湍流模型计算旋转力系数及其它耦合水动力系数.对“CR-02”6000 m自治水下机器人的计算表明,通过这种方法获得的水动力系数具有较高的精度,可以满足水下机器人方案设计阶段的操纵性设计、运动预报和仿真等需求.

Computation of Viscous Incompressible Flow using Pressure Correction Method on Unstructured Chimera Grid

In this paper, a pressure correction algorithm for computing incompressible flows is modified and implemented on unstructured Chimera grid. Schwarz method is used to couple the solutions of different sub-domains. A new interpolation to ensure consistency between primary variables and auxiliary variables is proposed. Other important issues such as global mass conservation and order of accuracy in the interpolations are also discussed. Two numerical simulations are successfully performed. They include one steady case, the lid-driven cavity and one unsteady case, the flow around a circular cylinder. The results demonstrate a very good performance of the proposed scheme on unstructured Chimera grids. It prevents the decoupling of pressure field in the overlapping region and requires only little modification to the existing unstructured Navier–Stokes (NS) solver. The numerical experiments show the reliability and potential of this method in applying to practical problems.

Damping of Vertically Excited Surface Wave in Weakly Viscous Fluid

In a vertically oscillating circular cylindrical container, singular perturbation theory of two-time scale expansions is developed in weakly viscous fluids to investigate the motion of single free surface standing wave by linearizing the Navier-Stokes equation. The fluid field is divided into an outer potential flow region and an inner boundary layer region. The solutions of both two regions are obtained and a linear amplitude equation incorporating damping term and external excitation is derived. The condition to appear stable surface wave is obtained and the critical curve is determined. In addition, an analytical expression of damping coefficient is determined. Finally, the dispersion relation, which has been derived from the inviscid fluid approximation, is modified by adding linear damping. It is found that the modified results are reasonably closer to experimental results than former theory. Result shows that when forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when forcing frequency is high, the surface tension of the fluid is prominent.

Three-dimensional instability of an oscillating viscous flow past a circular cylinder

A systematically numerical study of the sinusoidally oscillating viscous flow around a circular cylinder was performed to investigate vortical instability by solving the three-dimensional incompressible Navier-Stokes equations. The transition from two- to three-dimensional flow structures along the axial direction due to the vortical instability appears, and the three-dimensional structures lie alternatively on the two sides of the cylinder. Numerical study has been taken for the Keulegan-Carpenter( KC) numbers from 1 to 3.2 and frequency parameters from 100 to 600. The force behaviors are also studied by solving the Morison equation. Calculated results agree well with experimental data and theoretical prediction.

Vertically forced surface wave in weakly viscous fluids bounded in a circular cylindrical vessel

Two-time scale perturbation expansions were developed in weakly viscous fluids to investigate surface wave motions by linearizing the Navier-Stokes equation in a circular cylindrical vessel which is subject to a vertical oscillation. The fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates a damping term and external excitation, was derived for the weakly viscid fluids. The condition for the appearance of stable surface waves was obtained and the critical curve was determined. In addition, an analytical expression for the damping coefficient was determined and the relationship between damping and other related parameters (such as viscosity, forced amplitude, forced frequency and the depth of fluid, etc.) was presented. Finally, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent.

Numerical simulation of hypersonic viscous flow around space shuttle with method of diffusion analogy

Hypersonic viscous flow around a space shuttle with M(infinity) = 7, Re = 148000 and angle of attack alpha = 5-degrees is simulated numerically with the special Jacobian matrix splitting technique and simplified diffusion analogy method. With the simplified diffusion analogy method the efficiency of computation and resolution of the shock can be improved.

Distribution, stability, bifurcations and catastrophe of steady rotation of a symmetrical heavy gyroscope with viscous-liquid-filled cavity

The number, the angles of orientation and the stability in Rumyantsev Movchan's sense of oblique steady rotations of a symmetric heavy gyroscope with a cavity completely filled with a uniform viscous liquid, possessing a fixed point 0 on its symmetric axis. are given for various values of the parameters. By taking the square of the upright component of the angular momentum M2 as a control parameter, three types of bifurcation diagrams of the steady rotations, two types of jumps and two kinds of local catastrophes, one being the symmetric reduced cusp type and the other being of the symmetric reduced butterfly type, are obtained. By taking account of the M2-damping owing to the moment of unavoidable faint friction, two different modes for the gyroscope, initially in a stable quasi-steady upright rotation with a nutation angle theta(s) equal to zero, to topple over are found.

Temperature transformation in relativistic thermodynamics

This paper reviews conflicting results on relativistic transformation formula for tem-perature obtained by different authors in the last half century, discusses the proper expres-sion for elementary work done under reversible processes, and presents a correct derivationof the transformation formula, rather similar in spirit to that of Einstein. It is pointed outthat the point of view adopted by Eddington, Ott, Mφller and Landsberg are erroneous,and that by correctly carrying out detailed analysis for Mφller's working model, which wasoriginally proposed to disprove Planck-Einstein result, we have arrived instead at a result incomplete agreement with that of Planck-Einstein. Thercupon this long standing controversyover the temperature transformation dilemma for relativistic thermodynamics is clarified.

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.

Smoothed Particle Hydrodynamics (SPH): an Overview and Recent Developments

Smoothed particle hydrodynamics (SPH) is a meshfree particle method based on Lagrangian formulation, and has been widely applied to different areas in engineering and science. This paper presents an overview on the SPH method and its recent developments, including (1) the need for meshfree particle methods, and advantages of SPH, (2) approximation schemes of the conventional SPH method and numerical techniques for deriving SPH formulations for partial differential equations such as the Navier-Stokes (N-S) equations, (3) the role of the smoothing kernel functions and a general approach to construct smoothing kernel functions, (4) kernel and particle consistency for the SPH method, and approaches for restoring particle consistency, (5) several important numerical aspects, and (6) some recent applications of SPH. The paper ends with some concluding remarks.

Few-cycle relativistic solitons in plasmas

A set of exact one-dimensional solutions to coupled nonlinear equations describing the propagation of a relativistic ultrashort circularly polarized laser pulse in a cold collisionless and bounded plasma where electrons have an initial velocity in the laser propagating direction is presented. The solutions investigated here are in the form of quickly moving envelop solitons at a propagation velocity comparable to the light speed. The features of solitons in both underdense and overdense plasmas with electrons having different given initial velocities in the laser propagating direction are described. It is found that the amplitude of solitons is larger and soliton width shorter in plasmas where electrons have a larger initial velocity. In overdense plasmas, soliton duration is shorter, the amplitude higher than that in underdense plasmas where electrons have the same initial velocity.

Ion cascade acceleration from the interaction of a relativistic femtosecond laser pulse with a narrow thin target

Particle-in-cell simulations are performed to study the acceleration of ions due to the interaction of a relativistic femtosecond laser pulse with a narrow thin target. The numerical results show that ions can be accelerated in a cascade by two electrostatic fields if the width of the target is smaller than the laser beam waist. The first field is formed in front of the target by the central part of the laser beam, which pushes the electron layer inward. The major part of the abaxial laser energy propagates along the edges to the rear side of the target and pulls out some hot electrons from the edges of the target, which form another electrostatic field at the rear side of the target. The ions from the front surface are accelerated stepwise by these two electrostatic fields to high energies at the rear side of the target. The simulations show that the largest ion energy gain for a narrow target is about four times higher than in the case of a wide target. (c) 2006 American Institute of Physics.

Control of single attosecond pulse generation from the reflection of a synthesized relativistic laser pulse on a solid surface

The influence of the carrier-envelope phase (CEP) of the driving laser pulse on the generation of single attosecond (as) pulses from surface harmonics by using the polarization gating technique is investigated in detail. It is found that the modulation depth of the high-order harmonic spectrum depends on the CEP, and a strong single 68 as pulse can be generated when the CEP is stable and has the proper value. The physical origin of the influence of the CEP is explained in terms of the oscillating mirror model. (C) 2008 American Institute of Physics. [DOI: 10.1063/1.2997342]

Relativistic channel-coupling focusing enhanced nonlinear guiding of an intense laser beam in a plasma channel

We employ the variational method to study the optical guiding of an intense laser beam in a preformed plasma channel without using the weakly relativistic approximation. Apart from the dependence on the laser power and the nonlinear channel strength parameter, the beam focusing properties is shown also to be governed by the laser intensity. Relativistic channel-coupling focusing, arising from the coupling between relativistic self-focusing and linear channel focusing, can enhance relativistic self-focusing but its strength is weaker than that of linear channel focusing. (C) 2008 Elsevier B.V. All rights reserved.