Effect of nonlinear wave-current interaction on flow fields and hydrodynamic forces

A fifth-order theory for solving the problem of interaction between Stokes waves and exponential profile currents is proposed. The calculated flow fields are compared with measurements. Then the errors caused by the linear superposition method and approximate theory are discussed. It is found that the total wave-current field consists of pure wave, pure current and interaction components. The shear current not only directly changes the flow field, but also indirectly does sx, by changing the wave parameters due to wave-current interaction. The present theory can predict the wave kinematics on shear currents satisfactorily. The linear superposition method may give rise to more than 40% loading error in extreme conditions. When the apparent wave period is used and the Wheeler stretching method is adopted to extrapolate the current, application of the approximate theory is the best.

Numerical study of bluff body flow structures

Our recent progress in numerical studies of bluff body flow structures and a new method for the numerical analysis of near wake flow field for high Reynolds number flow are introduced. The paper consists of three parts. In part one, the evolution of wake vortex structure and variation of forces on a flat plate in harmonic oscillatory flows and in in-line steady-harmonic combined flows are presented by an improved discrete vortex method, as the Keulegan-Carpenter number (KC) varies from 2 to 40 and ratios of U-m to U-0 are of O(10(-1)), O(10) and O(10), respectively. In part 2, a domain decomposition hybrid method, combining the finite-difference and vortex methods for numerical simulation of unsteady viscous separated flow around a bluff body, is introduced. By the new method, some high resolution numerical visualization on near wake evolution behind a circular cylinder at Re = 10(2), 10(3) and 3 x 10(3) are shown. In part 3, the mechanism and the dynamic process for the three-dimensional evolution of the Karman vortex and vortex filaments in braid regions as well as the early features of turbulent structure in the wake behind a circular cylinder are presented numerically by the vortex dynamics method.

The mechanism of energy storage by shearing the solar magnetic field

In this paper, a complete set of MHD equations have been solved by numerical calculations in an attempt to study the dynamical evolutionary processes of the initial equilibrium configuration and to discuss the energy storage mechanism of the solar atmosphere by shearing the magnetic field. The initial equilibrium configuration with an arch bipolar potential field obtained from the numerical solution is similar to the configuration in the vicinity of typical solar flare before its eruption. From the magnetic induction equation in the set of MHD equations and dealing with the non-linear coupling effects between the flow field and magnetic field, the quantitative relationship has been derived for their dynamical evolution. Results show that plasma shear motion at the bottom of the solar atmosphere causes the magnetic field to shear; meanwhile the magnetic field energy is stored in local regions. With the increase of time the local magnetic energy increases and it may reach an order of 4×10^25 J during a day. Thus the local storage of magnetic energy is large enough to trigger a big solar flare and can be considered as the energy source of solar flares. The energy storage mechanism by shearing the magnetic field can well explain the slow changes in solar active regions.

Vortex motion in the early stages of unsteady flow around a circular cylinder

In this paper, processes in the early stages of vortex motion and the development of flow structure behind an impulsively-started circular cylinder at high Reynolds number are investigated by combining the discrete vortex model with boundary layer theory, considering the separation of incoming flow boundary layer and rear shear layer in the recirculating flow region. The development of flow structure and vortex motion, particularly the formation and development of secondary vortex and a pair of secondary vortices and their effect on the flow field are calculated. The results clearly show that the flow structure and vortices motion went through a series of complicated processes before the symmetric main vortices change into asymmetric: development of main vortices induces secondary vortices; growth of the secondary vortices causes the main vortex sheets to break off and causes the symmetric main vortices to become “free” vortices, while a pair of secondary vortices is formed; then the vortex sheets, after breaking off, gradually extend downstream and the structure of a pair of secondary vortices becomes relaxed. These features of vortex motion look very much like the observed features in some available flow field visualizations. The action of the secondary vortices causes the main vortex sheets to break off and converts the main vortices into free vortices. This should be the immediate cause leading to the instability of the motion of the symmetric main vortices. The flow field structure such as the separation position of boundary layer and rear shear layer, the unsteady pressure distributions and the drag coefficient are calculated. Comparison with other results or experiments is also made.

Flow pattern in chandler loop - the artificial thrombus formation

In this paper a hydrodynamic approach is used to analyse carefully the flow field inChandler loop--the artificial thrombus formation. The results obtained show that near thelower meniscus where the thrombus is formed, there is a back flow accompanied with asecondary flow and its mainflow is toward the meniscus, thus providing a favourable condi-tion for corpuscle aggregation. Our finding is valuable for studying the mechanism ofthrombus formation in artificial organ and in vivo.

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.

Modelling Study On The Plasma Flow and Heat Transfer In A Laminar Arc Plasma Torch Operating At Atmospheric and Reduced Pressure

A modelling study is performed to investigate the characteristics of both plasma flow and heat transfer of a laminar non-transferred arc argon plasma torch operated at atmospheric and reduced pressure. It is found that the calculated flow fields and temperature distributions are quite similar for both cases at a chamber pressure of 1.0 atm and 0.1 atm. A fully developed flow regime could be achieved in the arc constrictor-tube between the cathode and the anode of the plasma torch at 1.0 atm for all the flow rates covered in this study. However the flow field could not reach the fully developed regime at 0.1 atm with a higher flow rate. The arc-root is always attached to the torch anode surface near the upstream end of the anode, i.e. the abruptly expanded part of the torch channel, which is in consistence with experimental observation. The surrounding gas would be entrained from the torch exit into the torch interior due to a comparatively large inner diameter of the anode channel compared to that of the arc constrictor-tube.

FLOW-PIPE-SEEPAGE COUPLING ANALYSIS OF SPANNING INITIATION OF A PARTIALLY-EMBEDDED PIPELINE

The initiation of pipeline spanning involves the coupling between the flow over the pipeline and the seepage-flow in the soil underneath the pipeline. The pipeline spanning initiation is experimentally observed and discussed in this article. It is qualitatively indicated that the pressure-drop induced soil seepage failure is the predominant cause for pipeline spanning initiation. A flow-pipe-seepage sequential coupling Finite Element Method (FEM) model is proposed to simulate the coupling between the water flow-field and the soil seepage-field. A critical hydraulic gradient is obtained for oblique seepage failure of the sand in the direction tangent to the pipe. Parametric study is performed to investigate the effects of inflow velocity, pipe embedment on the pressure-drop, and the effects of soil internal friction angle and pipe embedment-to-diameter ratio on the critical flow velocity for pipeline spanning initiation. It is indicated that the dimensionless critical flow velocity changes approximately linearly with the soil internal friction angle for the submarine pipeline partially-embedded in a sandy seabed.

Multi-mode Vortex-induced Vibration of Slender Cable Experiencing Shear Flow

The dynamic characteristics of slender cable often present serried modes with low frequencies due to large structure flexibility resulted from high aspect ratio (ratio of length to diameter of cable), while the flow velocity distributes non-uniformly along the cable span actually in practical engineering. Therefore, the prediction of the vertex-induce vibration of slender cable suffered from multi-mode and high-mode motions becomes a challenging problem. In this paper a prediction approach based on modal energy is developed to deal with multi-mode lock-in. Then it is applied to the modified wake-oscillator model to predict the VIV displacement and stress responses of cable in non-uniform flow field. At last, illustrative examples are given of which the VIV response of flexible cable in nonlinear shear flow field is analyzed. The effects of flow velocity on VIV are explored. Our results show that both displacement and stress responses become larger as the flow velocity increasing; especially higher stress response companied with higher frequency vibration should be paid enough attention in practical design of SFT because of its remarkable influence on structure fatigue life.

Environmental Mechanics Research in China

By recalling mankind's path during past 50 years in the present article, we mainly highlight the significance of environmental issues today. In particular, two major factors leading to environment deterioration in China such as water resources and coal burning are stressed on. Present-day environmental issues are obviously interdisciplinary, of multiple scales and multi-composition in nature. Therefore, a process-based approach for environment research is absolutely necessarily. A series of sub-processes, either physical, chemical or biological, are subsequently analyzed in order to established reasonable parameterization scheme and credible comprehensive model. And we are now in a position to answer questions still open to us, improve existing somewhat empirical engineering approaches and enhance quantitative accuracy in prediction. To illustrate this process-based research approach, three typical examples associated with the Yangtze River Estuary, Loess Plateau and Tenggeli Desert environments have been dealt with respectively. A theoretical model of vertical flow field accounting for runoff and tide interaction has been established to delineate salinity and sediment motion which are responsible for the formation of mouth bar at the outlet and the ecological evolution there. A kinematic wave theory combined with the revised Green-Ampt infiltration formula is applied to the prediction of runoff generation and erosion in three types of erosion region on the Loess Plateau. Three approaches describing water motion in SPAC system in arid areas at different levels have been improved by introducing vegetation sub-models. However, we have found that the formation of a dry sandy layer and biological crust skin are additional primary causes leading to deterioration of water supply and succession of ecological system.

Forebody Compressibility Research of Hypersonic Vehicle

Three kinds of forebody model of hypersonic vehicles were studied with numerical simulation method. It shows that the two- order compressive ramp model is the best selection among the three for its good evaluative parameters value at the cowl of the inlet . This model can provide higher value of flux coefficient and total pressure recovery coefficient and lower average Mach number compared with those of the other two models . Simultaneously different compressive angles may have different effects . The configuration which the firstorder of compressive angle is 4°and the second 5°is the optimum combination. Furthermore factors such as attack angle were concerned. Better result may be obtained with a range of attack angles . Based on the work above the integrated design for forebodyPinlet of a hypersonic vehicle was performed. The numerical result shows that this integrated model provides good flow field quality for inlet and engine work.

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

Direct numerical simulation of transition and turbulence in compressible mixing layer

The three-dimensional compressible Navier-Stokes equations are approximated by a fifth order upwind compact and a sixth order symmetrical compact difference relations combined with three-stage Ronge-Kutta method. The computed results are presented for convective Mach number Mc = 0.8 and Re = 200 with initial data which have equal and opposite oblique waves. From the computed results we can see the variation of coherent structures with time integration and full process of instability, formation of Lambda-vortices, double horseshoe vortices and mushroom structures. The large structures break into small and smaller vortex structures. Finally, the movement of small structure becomes dominant, and flow field turns into turbulence. It is noted that production of small vortex structures is combined with turning of symmetrical structures to unsymmetrical ones. It is shown in the present computation that the flow field turns into turbulence directly from initial instability and there is not vortex pairing in process of transition. It means that for large convective Mach number the transition mechanism for compressible mixing layer differs from that in incompressible mixing layer.

Fourth order accurate compact scheme with group velocity control (GVC)

For solving complex flow field with multi-scale structure higher order accurate schemes are preferred. Among high order schemes the compact schemes have higher resolving efficiency. When the compact and upwind compact schemes are used to solve aerodynamic problems there are numerical oscillations near the shocks. The reason of oscillation production is because of non-uniform group velocity of wave packets in numerical solutions. For improvement of resolution of the shock a parameter function is introduced in compact scheme to control the group velocity. The newly developed method is simple. It has higher accuracy and less stencil of grid points.

COIL亚声速段横向喷流混合流场数值分析

COIL是一个气体动力学、化学反应动力学以及光学相互耦合的复杂过程.作为高总压COIL系统研究的第一步,利用三维CFD技术对传统的COIL亚声速段横向射流混合过程进行了数值分析,讨论了包括压力梯度驱动项的分子扩散机制,得到了横向射流的较精细的结构,如马蹄形射流界面、逆旋涡对以及射流剪切层.结果表明,压力梯度驱动项对重、轻组分的扩散作用相反,重组分沿压力梯度正向扩散,轻组分逆压力梯度方向扩散.在射流穿透不足的情况下,仍然在喷管出口得到了约为0.01 cm-1的小信号增益系数.