Features of the supercritical transition in the wake of a circular cylinder

The features of the wake behind a uniform circular cylinder at Re = 200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations by solving 3-D incompressible Navier-Stokes equations using mixed spectral-spectral-element method. The high-order splitting algorithm based on the mixed stiffly stable scheme is employed in the time discretization. Due to the nonlinear evolution of the secondary instability of the wake, the spanwise modes with different wavelengths emerge. The spanwise characteristic length determines the transition features and global properties of the wake. The existence of the spanwise phase difference of the primary vortices shedding is confirmed by Fourier analysis of the time series of the spanwise vorticity and attributed. to the dominant spanwise mode. The spatial energy distributions of various modes and the velocity profiles in the near wake are obtained. The numerical results indicate that the near wake is in 3-D quasi-periodic laminar state with transitional behaviors at this supercritical Reynolds number.

Semi-gas kinetics model for performance modeling of flowing chemical oxygen-iodine lasers (COIL)

A semi-gas kinetics (SGK) model for performance analyses of flowing chemical oxygen-iodine laser (COIL) is presented. In this model, the oxygen-iodine reaction gas flow is treated as a continuous medium, and the effect of thermal motions of particles of different laser energy levels on the performances of the COIL is included and the velocity distribution function equations are solved by using the double-parameter perturbational method. For a premixed flow, effects of different chemical reaction systems, different gain saturation models and temperature, pressure, yield of excited oxygen, iodine concentration and frequency-shift on the performances of the COIL are computed, and the calculated output power agrees well with the experimental data. The results indicate that the power extraction of the SGK model considering 21 reactions is close to those when only the reversible pumping reaction is considered, while different gain saturation models and adjustable parameters greatly affect the output power, the optimal threshold gain range, and the length of power extraction.

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise characteristic length linear stable mode may dominate the wake in place of mode A and determine the spanwise phase difference of the primary vortices shedding. At Re = 250 and 300 it is found that the streamwise vortices evolve into a new type of mode - "dual vortex pair mode" downstream. The streamwise vortex structures switch among mode A, mode B and dual vortex pair mode from near wake to downstream wake. At Re = 250, an independent low frequency f(m) in addition to the vortex shedding frequency f(s) is identified. Frequency coupling between f(m) and f(s) occurs. These result in the irregularity of the temporal signals and become a key feature in the transition of the wake. Based on the formation analysis of the streamwise vorticity in the vicinity of cylinder, it is suggested that mode A is caused by the emergence of the spanwise velocity due to three dimensionality of the incoming flow past the cylinder. Energy distribution on various wave numbers and the frequency variation in the wake are also described.

Compact finite difference - Fourier spectral method for three-dimensional incompressible Navier-Stokes equations

A new compact finite difference-Fourier spectral hybrid method for solving the three dimensional incompressible Navier-Stokes equations is developed in the present paper. The fifth-order upwind compact finite difference schemes for the nonlinear convection terms in the physical space, and the sixth-order center compact schemes for the derivatives in spectral space are described, respectively. The fourth-order compact schemes in a single nine-point cell for solving the Helmholtz equations satisfied by the velocities and pressure in spectral space is derived and its preconditioned conjugate gradient iteration method is studied. The treatment of pressure boundary conditions and the three dimensional non-reflecting outflow boundary conditions are presented. Application to the vortex dislocation evolution in a three dimensional wake is also reported.

Formation, microstructure and development of the localized shear deformation in low-carbon steels

An investigation has been made into the effect of microstructural parameters on the propensity for forming shear localization produced during high speed torsional testing by split Hopkinson bar with different average rates of 610, 650 and 1500 s(-1) in low carbon steels. These steels received the quenched, quenched and tempered as well as normalized treatments that provide wide microstructural parameters and mechanical properties. The results indicate that the occurrence of the shear localization is susceptible to the strength of the steels. In other words, the tendency of the quenched steel to form a shear band is higher than that of the other two steels. It is also found that there is a critical strain at which the shear localization occurs in the steels. The critical strain value is strongly dependent on the strength of the steels. Before arriving at this point, the material undergoes a slow work-hardening. After this point, the material suffers work-softening, corresponding to a process during which the deformation is gradually localized and eventually becomes spatially correlated to form a macroscopic shear band. Examinations by SEM reveal that the shear localization within the band involves a series of sequential crystallographic and non-crystallographic events including the change in crystal orientation, misorientation, generation and even perhaps damage in microstructures such as the initiation, growth and coalescence of the microcracks. It is expected that the sharp drop in the load-carrying capacity is associated with the growth and coalescence of the microcracks rather than the occurrence of the shear localization, but the shear localization is seen to accelerate the growth and coalescence of the microcracks. The thin foil observations by TEM reveal that the density of dislocations in the band is extremely high and the tangled arrangement and cell structure of dislocations tends to align along the shear direction. The multiplication and interaction of dislocations seems to be responsible for work-hardening of the steels. The avalanche of the dislocation cells corresponds to the sharp drop in shear stress at which the deformed specimen is broken. Double shear bands and kink bands are also observed in the present study. The principal band develops first and its width is narrower than that of the secondary band.

Predicting chaotic time-series with wavelet networks

A new technique, wavelet network, is introduced to predict chaotic time series. By using this technique, firstly, we make accurate short-term predictions of the time series from chaotic attractors. Secondly, we make accurate predictions of the values and bifurcation structures of the time series from dynamical systems whose parameter values are changing with time. Finally we predict chaotic attractors by making long-term predictions based on remarkably few data points, where the correlation dimensions of predicted attractors are calculated and are found to be almost identical to those of actual attractors.

Direct dynamical test for deterministic chaos and optimal embedding of a chaotic time-series

We propose here a local exponential divergence plot which is capable of providing an alternative means of characterizing a complex time series. The suggested plot defines a time-dependent exponent and a ''plus'' exponent. Based on their changes with the embedding dimension and delay time, a criterion for estimating simultaneously the minimal acceptable embedding dimension, the proper delay time, and the largest Lyapunov exponent has been obtained. When redefining the time-dependent exponent LAMBDA(k) curves on a series of shells, we have found that whether a linear envelope to the LAMBDA(k) curves exists can serve as a direct dynamical method of distinguishing chaos from noise.

Local exponential divergence plot and optimal embedding of a chaotic time-series

we propose here a local exponential divergence plot which is capable of providing a new means of characterizing chaotic time series. The suggested plot defines a time dependent exponent LAMBDA and a ''plus'' exponent LAMBDA+ which serves as a criterion for estimating simultaneously the minimal acceptable embedding dimension, the proper delay time and the largest Lyapunov exponent.

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.

Numerical experiments on one-dimensional model of turbulence

The initial-value problem of a forced Burgers equation is numerically solved by the Fourier expansion method. It is found that its solutions finally reach a steady state of 'laminar flow' which has no randomness and is stable to disturbances. Hence, strictly speaking, the so-called Burgers turbulence is not a turbulence. A new one-dimensional model is proposed to simulate the Navier-Stokes turbulence. A series of numerical experiments on this one-dimensional turbulence is made and is successful in obtaining Kolmogorov's (1941) k exp(-5/3) inertial-range spectrum. The (one-dimensional) Kolmogorov constant ranges from 0.5 to 0.65.

Singularity criteria for perturbation series

This paper is concerned with some mathematical aspects of the Van Dyke method inperturbation theory, i.e. the singularity criteria of perturbation series. The author suggestsa sign criterion and a Domb-syke plot for the cases with complex conjugate singularities, thussucceeding in extending the conclusions of Van Dyke's. Subsequently. effects of singularitiesof the lower order upon the criteria are taken into account. In addition, a method of locat-ing singular points is developed by analysing the new perturbation series derived by the Eulertransformation.

用Fourier变换莫尔法测量流体自由表面微幅波的振幅

采用一种非接触的光学方法——傅立叶变换莫尔法（Fourier transform method），结合数字图像处理技术，对微幅振荡的水表面波的振幅进行测量。它是对全场中每一个像素点进行测量，比接触测量法具有更高的灵敏度。它为微幅水表面波振幅的测量提供了一种手段。通过将计算机生成的周期性光栅图像经投影机直接投影到被测物体的参考平面，经CCD摄像头、图像板捕捉存储形成数字化的光栅图像，利用傅立叶变换莫尔法处理光栅图像，从而获得包含有水表面波的振幅的相位信息，再经适当的几何变换获得振幅信息。我们在生趣振荡装置上进行了不同激励频率和不同振幅的表面波的振幅测量。

Experimental Study on Bénard-Marangoni Convection in Double Immiscible Liquid Layer

An experimental investigation of Bénard-Marangoni convection has been performed in double immiscible liquid layers of rectangular configuration on the ground. The two kinds of liquid are 10cst silicon oil and FC-70 respectively. The size of rectangular chamber is 100mm×40mm in horizontal cross-section. The evolution processes of convection are observed in the differential thickness ratio of two liquid layers. The critical temperature difference was measured via the detections of fluid convection by a particle image velocimetry (PIV) in the vertical cross-section of the liquid layer. The critical temperature difference or the critical Marangoni number was given. And the influence of the thickness ratio of two liquid layers on the convection instability was discussed. The evolution processes of patterns and temperature distributions on the interface are displayed by using thermal liquid crystal. The velocity distributions on the interface were also obtained. In comparison with the thermocapillary effect, the effect of buoyancy convection will relatively increase when the depth of the liquid layer increases. Because of the coupling of buoyancy and thermocapillary effect, the convection instability is much more complex than that in the microgravity environment. And the critical convection depends on the change of the thickness of liquid layers and also the change of thickness ratio of two liquid layers.

Numerical Partition of Energy Absorption of Foam-filled Top-Hat and Double-Hat Sections

The interaction effect, i.e., the contribution of each component to the total energy absorption of an axially crushed foam-filled hat section was investigated quantitatively via numerical simulation. The FE results were first verified by experimental work of aluminum foam-filled top-hat and double-hat sections, then the contribution of foam-fillers and that of hat sections to the overall energy absorption were quantitatively obtained, respectively. When foam-filled, increase in energy absorption was found both in hat section component and foam-filler component, whereas the latter contributes predominantly to the interaction effect.