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.

A perturbational h4 exponential finite-difference scheme for the convective diffusion equation

A perturbational h4 compact exponential finite difference scheme with diagonally dominant coefficient matrix and upwind effect is developed for the convective diffusion equation. Perturbations of second order are exerted on the convective coefficients and source term of an h2 exponential finite difference scheme proposed in this paper based on a transformation to eliminate the upwind effect of the convective diffusion equation. Four numerical examples including one- to three-dimensional model equations of fluid flow and a problem of natural convective heat transfer are given to illustrate the excellent behavior of the present exponential schemes, the h4 accuracy of the perturbational scheme is verified using double precision arithmetic.

Crystal-growth in floating zone with phase-change and thermosolutal convections

Floating zone crystal growth in microgravity environment is investigated numerically by a finite element method for semiconductor growth processing, which involves thermocapillary convection, phase change convection, thermal diffusion and solutal diffusion. The configurations of phase change interfaces and distributions of velocity, temperature and concentration fields are analyzed for typical conditions of pulling rates and segregation coefficients. The influence of phase change convection on the distribution of concentration is studied in detail. The results show that the thermocapillary convection plays an important role in mixing up the melt with dopant. The deformations of phase change interfaces by thermal convection-diffusion and pulling rods make larger variation of concentration field in comparison with the case of plane interfaces.

Coupling of melting-solidification phase-change convection with floating-zone convection under microgravity

A finite element algorithm is used to analyze the process of floating zone crystal growth under microgravity. The effect of phase change convection coupled with surface tension convection is considered. The results show that the rate of crystal growth is very important. The single-crystal-melt interface is steeper than the feed-melt interface during the process of crystal growth. When the rate exceeds a critical value, the Marangoni vortex near the feed-melt interface will become so large that a secondary vortex will exist.

A finite-element analysis of viscoelastically damped sandwich plates

A finite element analysis associated with an asymptotic solution method for the harmonic flexural vibration of viscoelastically damped unsymmetrical sandwich plates is given. The element formulation is based on generalization of the discrete Kirchhoff theory (DKT) element formulation. The results obtained with the first order approximation of the asymptotic solution presented here are the same as those obtained by means of the modal strain energy (MSE) method. By taking more terms of the asymptotic solution, with successive calculations and use of the Padé approximants method, accuracy can be improved. The finite element computation has been verified by comparison with an analytical exact solution for rectangular plates with simply supported edges. Results for the same plates with clamped edges are also presented.

An asymptotic analysis for one dimensional stress wave propagation in damaged viscoelastic media

A regular perturbation technique is suggested to deal with the problem of one dimensional stress wave propagation in viscoelastic media with damage. Based upon the first order asymptotic solution obtained, the characteristics of wave attenuation are studied. In fact, there exist three different time-dependent phenomena featuring the dynamic response of the materials, the first expressing the characteristics of wave propagation, the second indicating the innate effect of visco-elastic matrix and the third coming from the time dependent damage. The comparision of first order asymptotic solution with the numerical results calculated by a finite difference procedure shows that the perturbation expansion technique may offer a useful approach to the problem concerned.

Standing wave of finite-amplitude in an circular basin with uneven bottom

The effect of the particle cover over the density interface between two layers of fluids and of the suspended solid particles in the upper turbulcnt layer on the turbulent entrainment has been studied experimentally. The entrainment distance D is a function of the time of power: D=kt, where =0.200-0.130p. For suspended particles in the upper layer and pure 2-layer fluid is equal to 0.200, but the value of k for the suspended particles is smaller than that for the pure 2-layer fluid. The non-dimensional entrainment velocity is E=KRiln, where n=1.50+0.93 p. It is shown that the particle cover over the interface changes the power of Ril in the entrainment and hinders the turbulent entrainment. The variation rule of E for the suspended particles is the same as that for the pure 2-layer fluid, but the K value of the former is smaller than that of the latter. The turbulent mixing mechanism has been discussed.

Study on cracked plates, shells and 3 dimensional bodies

This paper presents a summary of the authors' recent work in following areas: (1) The stress-strain fields at crack tip in Reissner's plate. (2) The calculations of the stress intensity factors in finite size plates. (3) The stress-strain fields at crack tip in Reissner's shell. (4) The calculations of the stress intensity factors and bulging coefficients in finite size spherical shells. (5) The stress-strain fields along crack tip in three dimensional body with surface crack. (6) The calculation of stress intensity factors in a plate with surface crack.

Characteristics of a CW HF chemical-laser calculated from a simplified two-dimensional model

A two-dimensional simplified model of an HF chemical laser is introduced. Using an implicit finite difference scheme, the solution of two adjacent parallel streams with diffusion mixing and chemical reaction is generated. A contour of mixing and reaction boundary is obtained without presupposition. The distribution of the HF(v) concentrations, gas temperature and the optical small signal gain (alpha sub V, J) on the flowing plane (X, Y) are presented. Compared with the solution solved directly from a set of Navier-Stokes equations, the results of these two methods agree with each other qualitatively. The influences of the different velocity, temperature (T sub 0) and composition of the two streams on the small signal gain after the nozzle exit are investigated. It is interesting that for larger J with a fixed v, the peaks of alpha sub v-T sub 0 profiles move towards higher T sub 0. The computing method is simple and only a short computing time is needed.

Short-surface waves radiated by a partially-immersed 3-dimensional oscillating body

The short-surface waves generated by a 3-D arbitrarily oscillating body floating onwater are discussed. In the far-field off the body, the phase and the amplitude functions ofthe radiated waves are determined by the ray method. An undetermined constant is includ-ed in the amplitude function. From the result of Ref. [1], the near-field boundary layersolution near the body waterline is obtained. The amplitude of this solution depends on thewhole wall shape of the body and the slope at the body waterline on the cross-sections per-pendicular to the waterline. By matching the far-field solution with the near-field bound-ary layer solution, the undetermined constant in the amplitude function of the far-fieldradiated waves is determined. For the special case of a half-submerged sphere which per-forms vertical oscillating motion, the result obtained in this paper is in agreement withthat of Ref. [ 2 ].

Nonequilibrium statistical mechanics of two-dimensional turbulence

The vorticity dynamics of two-dimensional turbulence are investigated analytically, applying the method of Qian (1983). The vorticity equation and its Fourier transform are presented; a set of modal parameters and a modal dynamic equation are derived; and the corresponding Liouville equation for the probability distribution in phase space is solved using a Langevin/Fokker-Planck approach to obtain integral equations for the enstrophy and for the dynamic damping coefficient eta. The equilibrium spectrum for inviscid flow is found to be a stationary solution of the enstrophy equation, and the inertial-range spectrum is determined by introducing a localization factor in the two integral equations and evaluating the localized versions numerically.

3D Finite Volume Scheme for Czochralski Crystal Growth

A general three-dimensional model is developed for simulation of the growth process of silicon single crystals by Czochralski technique. The numerical scheme is based on the curvilinear non-orthogonal finite volume discretization. Numerical solutions show that the flow and temperature fields in the melt are asymmetric and unsteady for 8’’ silicon growth. The effects of rotation of crystal on the flow structure are studied. The rotation of crystal forms the Ekman layer in which the temperature gradient along solid/melt surface is small.

A smoothing technique for discrete delta functions with application to immersed boundary method in moving boundary simulations.

The effects of complex boundary conditions on flows are represented by a volume force in the immersed boundary methods. The problem with this representation is that the volume force exhibits non-physical oscillations in moving boundary simulations. A smoothing technique for discrete delta functions has been developed in this paper to suppress the non-physical oscillations in the volume forces. We have found that the non-physical oscillations are mainly due to the fact that the derivatives of the regular discrete delta functions do not satisfy certain moment conditions. It has been shown that the smoothed discrete delta functions constructed in this paper have one-order higher derivative than the regular ones. Moreover, not only the smoothed discrete delta functions satisfy the first two discrete moment conditions, but also their derivatives satisfy one-order higher moment condition than the regular ones. The smoothed discrete delta functions are tested by three test cases: a one-dimensional heat equation with a moving singular force, a two-dimensional flow past an oscillating cylinder, and the vortex-induced vibration of a cylinder. The numerical examples in these cases demonstrate that the smoothed discrete delta functions can effectively suppress the non-physical oscillations in the volume forces and improve the accuracy of the immersed boundary method with direct forcing in moving boundary simulations.