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 developed 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.

Non-modal instabilities of two-dimensional disturbances in plane Couette flow of a power-law fluid

Instabilities of fluid flows have traditionally been investigated by normal mode analysis, i.e. by linearizing the equations of flow and testing for unstable eigenvalues of the linearized problem. However, the results of eigenvalue analysis agree poorly in many cases with experiments, especially for shear flows. In this paper we study the instabilities of two-dimensional Couette flow of a polymeric fluid in the framework of non-modal stability theory rather than normal mode analysis. A power-law model is used to describe the polymeric liquid. We focus on the response to external excitations and initial conditions by examining the pseudospectra structures and the transient energy growths. For both Newtonian and non-Newtonian flows, the results show that there can be a rather large transient growth even though the linear operator of Couette flow has no unstable eigenvalue. The effects of non-Newtonian viscosity on the transient behaviors are examined in this study. The results show that the "shear-thinning/shear-thickening" effect increases/decreases the amplitude of responses to external excitations and initial conditions. (C) 2010 Elsevier B.V. All rights reserved.

Underwater Acoustics and Cavitating Flow of Water Entry

The fluid mechanics of water entry is studied through investigating the underwater acoustics and the supercavitation. Underwater acoustic signals in water entry are extensively measured at about 30 different positions by using a PVDF needle hydrophone. From the measurements we obtain (1) the primary shock wave caused by the impact of the blunt body on free surface; (2) the vapor pressure inside the cavity; (3) the secondary shock wave caused by pulling away of the cavity from free surface; and so on. The supercavitation induced by the blunt body is observed by using a digital high-speed video camera as well as the single shot photography. The periodic and 3 dimensional motion of the supercavitation is revealed. The experiment is carried out at room temperature.

固井工程中的流动问题

同井是油气在建过程中的重要工程，一般分下套管和注水泥两个过程．本文结合作者近年来的研究，着重阐述了钻井液、前置液、水泥浆在环形空间中流动的流体力学问题，包括：高温高压下的水泥浆流变性、偏心环形空间中浆液的流动特性、钻井液-前置液和前置液-水泥浆界面的稳定性等，介绍了国内外学者对这些问题的研究进展，提出了今后应深入研究的问题.

Characters of Surface Deformation and Surface Wave in Thermal Capillary Convection

In the field of fluid mechanics, free surface phenomena is one of the most important physical processes. In the present research work, the surface deformation and surface wave caused by temperature difference of sidewalls in a rectangular cavity have been investigated. The horizontal cross-section of the container is 52 mmx42 mm, and there is a silicon oil layer of height 3.5 mm in the experimental cavity. Temperature difference between the two side walls of the cavity is increased gradually, and the flow on the liquid layer will develop from stable convection to un-stable convection. An optical diagnostic system consisting of a modified Michelson interferometer and image processor has been developed for study of the surface deformation and surface wave of thermal capillary convection. The Fourier transformation method is used to interferometer fringe analysis. The quantitative results of surface deformation and surface wave have been calculated from a serial of the interference fringe patterns.The characters of surface deformation and surface wave have been obtained. They are related with temperature gradient and surface tension. Surface deformation is fluctuant with time, which shows the character of surface wave. The cycle period of the wave is 4.8 s, and the amplitudes are from 0 to 0.55 mu m. The phase of the wave near the cool side of the cavity is opposite and correlative to that near the hot side. The present experiment proves that the surface wave of thermal capillary convection exists on liquid free surface, and it is wrapped in surface deformation.

Pulverized-coal combustion in Wall-Protecting-Jets combustor (WPJC)

A new type of pulverized-coal combustor, called "Wall-Protecting-Jets Combustor" (hereafter, WPJC has been proposed, designed and studied with both CFD (Computational Fluid Dynamics) and experimental methods. The WPJC is based on a novel concept in which all inlet jets are along the combustor wall. Pilot combustion experiments were conducted to investigate the combustion performance of WPJC. Two-phase flows and pulverized-coal combustion were simulated to study the mechanism of),WPJC using the commercial software FLUENT. The results show that the WPJC has many remarkable advantages: wall-protection by the cold jets without the use of refractory materials; low-temperature and three-stage combustion with low NOx emission; negligible ash/slag-deposition; multiple functions with convenient switching between them; effective adjustment of the combustion intensity and the ignition position.

High Aspect Ratio (L/D) Riser Viv Prediction Using Wake Oscillator Model

A two-dimensional (2-D) vortex-induced vibration (VIV) prediction model for high aspect ratio (LID) riser subjected to uniform and sheared flow is studied in this paper. The nonlinear structure equations are considered. The near wake dynamics describing the fluctuating nature of vortex shedding is modeled using classical van der Pol equation. A new approach was applied to calibrate the empirical parameters in the wake oscillator model. Compared the predicted results with the experimental data and computational fluid dynamic (CFD) results. Good agreements are observed. It can be concluded that the present model can be used as simple computational tool in predicting some aspects of VIV of long flexible structures. (C) 2008 Elsevier Ltd. All rights reserved.

The Real-Time Mach-Zehnder Interferometer Used In Space Experiment

The optical interference method is a promising technique for measuring temperature, density, and concentration in fluids. The non-intrusive and non-invasive nature of its optical techniques to the measured section are its most outstanding features. However, the adverse experiment environment, especially regarding shaking and vibrating, greatly restricts the application of the interferometer. In the present work, an optical diagnostic system consisting of a Mach-Zehnder interferometer (named after physicists Ludwig Mach) and an image processor has been developed that increases the measuring sensitivity compared to conventional experimental methods in fluid mechanics. An image processor has also been developed for obtaining quantitative results by using Fourier transformation. The present facility has been used in observing and measuring the mass transfer process of a water droplet in EAFP protein solution under microgravity condition provided by the satellite Shi Jian No. 8.

Compact difference approximation with consistent boundary condition

For simulating multi-scale complex flow fields it should be noted that all the physical quantities we are interested in must be simulated well. With limitation of the computer resources it is preferred to use high order accurate difference schemes. Because of their high accuracy and small stencil of grid points computational fluid dynamics (CFD) workers pay more attention to compact schemes recently. For simulating the complex flow fields the treatment of boundary conditions at the far field boundary points and near far field boundary points is very important. According to authors' experience and published results some aspects of boundary condition treatment for far field boundary are presented, and the emphasis is on treatment of boundary conditions for the upwind compact schemes. The consistent treatment of boundary conditions at the near boundary points is also discussed. At the end of the paper are given some numerical examples. The computed results with presented method are satisfactory.

Difference schemes on non-uniform mesh and their application

High order accurate schemes are needed to simulate the multi-scale complex flow fields to get fine structures in simulation of the complex flows with large gradient of fluid parameters near the wall, and schemes on non-uniform mesh are desirable for many CFD (computational fluid dynamics) workers. The construction methods of difference approximations and several difference approximations on non-uniform mesh are presented. The accuracy of the methods and the influence of stretch ratio of the neighbor mesh increment on accuracy are discussed. Some comments on these methods are given, and comparison of the accuracy of the results obtained by schemes based on both non-uniform mesh and coordinate transformation is made, and some numerical examples with non-uniform mesh are presented.

Waverider configurations derived from general conical flowfields

A method based on the computational fluid dynamics (CFD) is presented for a flexible waverider's design. The generating bodies of this method could be any cones. In addition, either the leading edge or the profile of the scramjet's inlet is used as the waverider's definition curve, parameterized by the quadric function, the sigmoid function or the B-spline function. Furthermore, several numerical examples are carried out to validate the method and the relevant codes. The CFD results of the configurations show that all the designs are successful. Moreover, primary suggestions are proposed for practical design by comparing the geometrical and aerodynamic performances of the cone-derived waveriders at Mach 6.

Numerical analysis of spatial evolution of the small signal gain in a chemical oxygen-iodine laser operating without primary buffer gas

A chemical oxygen iodine laser (COIL) that operates without primary buffer gas has become a new way of facilitating the compact integration of laser systems. To clarify the properties of spatial gain distribution, three-dimensional (3-D) computational fluid dynamics (CFD) technology was used to study the mixing and reactive flow in a COIL nozzle with an interleaving jet configuration in the supersonic section. The results show that the molecular iodine fraction in the secondary flow has a notable effect on the spatial distribution of the small signal gain. The rich iodine condition produces some negative gain regions along the jet trajectory, while the lean iodine condition slows down the development of the gain in the streamwise direction. It is also found that the new configuration of an interleaving jet helps form a reasonable gain field under appropriate operation conditions. (c) 2007 Elsevier Ltd. All rights reserved.

Free surface vibration in oscillatory convection of half floating zone

Projecting an orthographical grating mask (20pl/mm) on the surface of a small liquid bridge and receiving the reflected distortion image, one can calculate out reversely the shape of free surface of a liquid bridge. In this way we measured the surface shape of a small floating zone and the two-dimensional deformation of its vibration. The mechanism of thermocapillary oscillatory convection and the three-dimensional variation of the free surface are revealed experimentally. The principle for space experiment has been studied in our laboratory.

H4 exponential finite-difference scheme for convective diffusion-equations

Perturbations are applied to the convective coefficients and source term of a convection-diffusion equation so that second-order corrections may be applied to a second-order exponential scheme. The basic Structure of the equations in the resulting fourth-order scheme is identical to that for the second order. Furthermore, the calculations are quite simple as the second-order corrections may be obtained in a single pass using a second-order scheme. For one to three dimensions, the fourth-order exponential scheme is unconditionally stable. As examples, the method is applied to Burgers' and other fluid mechanics problems. Compared with schemes normally used, the accuracies are found to be good and the method is applicable to regions with large gradients.