A high order accurate difference scheme for complex flow fields

A high order accurate finite difference method for direct numerical simulation of coherent structure in the mixing layers is presented. The reason for oscillation production in numerical solutions is analyzed, It is caused by a nonuniform group velocity of wavepackets. A method of group velocity control for the improvement of the shock resolution is presented. In numerical simulation the fifth-order accurate upwind compact difference relation is used to approximate the derivatives in the convection terms of the compressible N-S equations, a sixth-order accurate symmetric compact difference relation is used to approximate the viscous terms, and a three-stage R-K method is used to advance in time. In order to improve the shock resolution the scheme is reconstructed with the method of diffusion analogy which is used to control the group velocity of wavepackets. (C) 1997 Academic Press.

Periodic vortex breakdown in wide spherical gaps

The flow field with vortex breakdown in wide spherical gaps was studied numerically by a finite difference method under the axisymmetric condition. The result shows that the flow bifurcates to periodic motion as the Reynolds number or the eccentricity of the spheres increases. (C) 1997 American Institute of Physics.

Dynamic analysis of arrest of buckle propagation on a beam on a nonlinear elastic foundation by FEM

Based on the dynamic governing equation of propagating buckle on a beam on a nonlinear elastic foundation, this paper deals with an important problem of buckle arrest by combining the FEM with a time integration technique. A new conclusion completely different from that by the quasi-static analysis about the buckle arrestor design is drawn. This shows that the inertia of the beam cannot be ignored in the analysis under consideration, especially when the buckle propagation is suddenly stopped by the arrestors.

Numerical study of bifurcation solutions of spherical Taylor-Couette flow

The steady bifurcation flows in a spherical gap (gap ratio sigma=0.18) with rotating inner and stationary outer spheres are simulated numerically for Re(c1)less than or equal to Re less than or equal to 1 500 by solving steady axisymmetric incompressible Navier-Stokes equations using a finite difference method. The simulation shows that there exist two steady stable flows with 1 or 2 vortices per hemisphere for 775 less than or equal to Re less than or equal to 1 220 and three steady stable flows with 0, 1, or 2 vortices for 1 220

Influence of low-gravity level on crystal-growth in floating-zone

The influence of low gravity level on crystal growth in the floating zone, which involves thermocapillary convection, phase change convection, thermal and solutal diffusion, is investigated numerically by a finite element method for the silicon crystal growth process. The velocity, temperature, concentration fields and phase change interfaces depending on heating temperature and growth rates are analyzed. The influence of low gravity level on the concentration is studied especially. The results show that the non-uniformities of concentration are about 10(-3) for growth rate nu(p) = 5.12 x 10(-8) m/s, 10(-2) for nu(p) = 5.12 x 10(-7) m/s and relatively larger for larger growth rate in the gravity level g = 0-9.8 m/s2. The thermocapillary effect is strong in comparison with the Bridgman system, and the level of low gravity is relatively insensitive for lower growth rates.

Equilibrium of isolated magnetostatic flux structure

Two-dimensional magnetostatic models of flux structure confined in stratified atmosphere are discussed in the present paper. The magnetic field in the flux structure is assumed to be force-free at the first step. Numerical solutions for this nonlinear free boundary problem are obtained by finite element method. Results show clearly the relation between the inside fields and outside pressure, especially the influence of atmospheric pressure distribution on the flux structure.

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.

Dynamic interaction of an elastic container with fluid

The vibration analysis of an elastic container with partially filled fluid was investigated in this paper. The container is made of a thin cylinder and two circular plates at the ends. The axis of the cylinder is in the horizontal direction. It is difficult to solve this problem because the complex system is not axially symmetric. The equations of motion for this system were derived. An incompressible and ideal fluid model is used in the present work. Solutions of the equations were obtained by the generalized variational method. The solution was expressed in a series of normalized generalized Fourier's functions. This series converged rapidly, and so its approximate solution was obtained with high precision. The agreement of the calculated values with the experimental result is good. It should be mentioned that with our method, the computer time is less than that with the finite-element method.

Excitation mechanism of thermocapillary oscillatory convection

It is suggested that the oscillation of thermocapillary convection may be excited by the buoyancy instability. By means of numerical simulation of the finite-element method, the temperature distributions in the liquid bridge are qualitatively analyzed. The temperature gradient in a certain flow region of liquid bridge may turn to be parallel to the direction of gravity when the temperature difference △T between two boundary rods of liquid bridge is larger than the critical value. The buoyancy instability may be excited, and then the thermocapillary oscillatory convection appears, as the temperature difference increases further. The distribution of the critical Marangoni number in the micro-gravity environment is derived from the data on the ground experiments. The results show that the onset of thermocapillary oscillatory convection is delayed in the case of smaller typical scale of liquid bridge and lower gravity environment.

Prediction of mixed-mode cracking direction in random, short-fibre composite-materials

The prediction of cracking direction in composite materials is of significance to the design of composite structures. This paper presents several methods for predicting the cracking direction in the double grooved tension-shear specimen which gives mixed-mode cracking. Five different criteria are used in this analysis: two of them have been used by other investigators and the others are proposed by the present authors. The strain energy density criterion proposed by G.C. Sih is modified to take account of the influence of the anisotropy of the strength on the direction of crack. The two failure criteria of Tsai-Hill and Norris are extended to predict the crack orientation. The stress distributions in the near-notch zone are calculated by using the 8-node quadrilateral isoparametric finite element method. The predictions of all the criteria except one are in good agreement with the experimental measurement. In addition, on the basis of the FEM results, the size of the zone in which the singular term is dominant is estimated.

永乐大钟及悬挂支撑系统的撞钟过程瞬态分析

采用ANSYS有限元分析软件对永乐大钟及其悬挂支撑系统建立了有限元模型,进行了撞钟过程有限元瞬态分析及动力强度校核。通过有限元计算,获得了大钟各局部考察对象及整体的应力、位移分布情况和各部分的动力响应和对强度的影响,为合理撞钟和加固提供了科学依据与技术指导。

Chemical Nonequilibrium Effects on Flow Field for Reusable Launch Vehicles

High temperature chemical non-equilibrium phenomena have a great effect on the flow field around a reentry vehicle. A set of three dimensional Navier-Stokes equations have been solved by implicit finite volume NND scheme. Both ideal gas viscous flow and chemical non-equilibrium flow are calculated for a spherical-cone at a small angle of attack. The results of the two flows have been compared and the effect of chemical non-equilibrium has been analyzed. The effect of wall material's properties, such as catalysis and radiation were studied. The results are in good agreement with the referenced paper.

Integration property of a wedge-elliptic cone waverider with scramjet engine

The hypersonic waverider forebody is designed in this paper. For the present waverider, the undersurface is carved out as a stream surface of a hypersonic inviscid flow field around wedge-elliptic cone, and the upper surface is assumed to be a freestream surface. A finite-volume code is used to generate the three-dimensional flow field. The leading edge is determined by satisfying the condition that the lip is situated at the intersection line of shocks.

Torsion fracture behavior of drawn pearlitic steel wires with different heat treatments

In this paper, torsion fracture behavior of drawn pearlitic steel wires with different heat treatments was investigated. Samples with different heat treatments was investigated. Samples with different heat treatment conditions were subjected to torsion and tensile tests. The shear strain along the torsion sample after fracture was measured. Fracture surface of wires was examined by Scanning Electron Microscopy. In addition, the method of Differential Scanning Calorimetry was used to characterize the thermodynamic process in the heat treatment. A numerical simulation via finite element method on temperature field evolution for the wire during heat treatment process was performed. The results show that both strain aging and recovery process occur in the material within the temperature range between room temperature and 435 degrees C. It was shown that the ductility measured by the number of twists drops at short heating times and recovers after further heating in the lead bath of 435 degrees C. On the other hand, the strenght of the wire increases at short heating times and decreases after further heating. The microstructure inhomogeneity due to short period of heat treatment, coupled with the gradient characteristics of shear deformation during torsion results in localized shear deformation of the wire. In this situation, shear cracks nucleate between lamella and the wire breaks with low number of twists.

Modelling rail corrugation with specific track parameters focusing on ballasted track and slab track

The objective of this paper is to compare 3 types of track (high performance ballasted track, STEDEF and AFTRAV) from the corrugation growth point of view. This work has considered different vehicle speeds and track radii, and the results have taken into account the four wheels of a bogie. These tracks have been studied using Finite Elements with Nastran-Patran and RACING, a tool developed in Matlab by the authors which estimates the corrugation growth tendency. The tracks are studied using the Finite Strip Method and the Periodic Structure Theory. Lateral and vertical receptances for track and vehicle have been obtained, as well as the corrugation growth functions. In the paper the tracks are ranked according to corrugation development.