The aim of this paper is to develop a CFD solver that used to simulate the transient flow phenomena in a reflected shock tunnel. The transient flow phenomena in a shock tunnel include the reflected shock/boundary layer interaction and the starting process of nozzle flow that can affect the duration of test flow in actual conditions. To numerically simulate these transient flow features, a full-size, axisymmetric reflected shock tunnel model is used. The governing equations are a full Navier-Stokes equation, a species equation and a simplified polynomial correlation to simulate the real gas effects. The numerical code is developed based on the finite volume method coupled with the upwind Roe's scheme and the total variation diminishing (TVD) method. To increase the calculation efficiency, a multi-block and multi-mesh grid generation technique is employed in a huge computational domain. The present computational results have not only confirmed the theoretical characteristics of a shock tube, but have also qualitatively presented the phenomena of reflected shock/boundary layer interaction and the starting process of nozzle flow. This numerical code is a useful tool to demonstrate the actual flow phenomena and to assist the design of experiments.

Direct numerical simulations of 2D turbulent flows, freely decaying as well as forced, are performed to examine the mechanism of the enstrophy cascade and serve as a template of developing LES models. The stretching effect on the 2D vorticity gradients is emphasized on the analogy of the stretching effect on 3D vorticity. The enstrophy cascade rate, the Reynolds stresses and the associated eddy viscosity for 2D turbulence are correspondingly derived and investigated. Proposed herein is that the enstrophy cascade rate to be modeled in a large-eddy simulation can be and should be calculated using the only available large-eddy information, especially when the Reynolds number is not very large or when the flow is not stationary.

The simulation results suggest all Kolmogorov's, Kraichnan's, and Saffman's similarity spectra. The Kolmogorov's spectrum appears in front of forced wave numbers and creates a subrange of a zero enstrophy cascade rate and a constant energy cascade rate. The Saffman's spectrum is the dissipation spectrum at large wave numbers. Kraichnan's spectrum shows up at intermediate wave numbers when the Reynolds number is sufficiently high. When the Smagorinsky model is employed for a large eddy simulation, its inability of capturing the significant reverse cascade phenomenon as observed in the DNS data becomes a fatal defect. Nonetheless, if only the mean cascade rate is concerned, the required Smagorinsky constant is evaluated using the DNS data and compared with the theoretical prediction of the Kraichnan's spectrum.

The effects of blockage and swirl on the macro flow structures of the annular jet past a circular disc are experimentally studied through the time-averaged streamline patterns. In the blockage-effect regime, the flows present multiple modes, single bubble, dual rings, vortex breakdown, and triple rings, in different regimes of blockage ratio and swirl number. The topological models of the flow structures are proposed and discussed according to the measured flow fields to manifest the complex flow structures. The single bubble is a closed recirculation bubble with a stagnation point on the central axis. The dual-ring flow is an open-top recirculsation zone, in which a pair of counter-rotating vortex rings exists in the near wake. The fluids in the dual rings are expelled downstream through a central jet-like swirling flow. A vortex breakdown may occur in the central jet-like swirling flow if the exit swirl number exceeds critical values. When the vortex breakdown interacts with the dual rings, a complex triple-ring flow structure forms. Axial distributions of the local swirl number are presented and discussed. The local swirl number increases with the increase of the exit swirl number and attains the maximum in the dual-ring mode. At large exit swirl numbers where the vortex breakdown occurs, the local swirl number decreases drastically to a low value.

The flow-induced vibration of tubes in a rotated triangular array subject to cross flow is simulated numerically. In the study, the flow field around the tube bundle is computed by solving the continuity and Navier-Stokes equations with assumption of constant fluid properties, and the kε-model for turbulent Reynolds stress. With the flow field known, the fluid forces on the tube surfaces can be calculated, and then the displacement of each tube due to the fluid force can be evaluated. Iteration is needed to obtain the dynamic response of the tube structure in the fluid flow. The parameters in the study are inlet velocity of the cross flow and properties of the tube bundle including natural frequency, damping factor, and mass. Based on the tube response, the critical flow conditions of tube vibration are determined for varying mass damping. Once tube vibrations occur, it is shown that the vibrations of the tubes in the second and fourth tube rows are significant as compared to other tubes. The orbits of the tube vibration look like an ellipse with major axis in the cross-stream direction, implying large lift force on the tubes. The dominant frequency in the spectrum of lift coefficients of the tubes is the same as the natural frequency, and the corresponding amplitude is increased with increasing the inlet velocity. The calculated data predicted for the critical reduced velocity agrees well with the data by Kassera and Strohmeier [17].

Exact electroelastic analysis of hollow and solid circular cylinders subjected to a uniform temperature change and electromechanical loads that do not vary axially are presented. Generalized plane strain and torsion of a piezoelectric cylindrical body is formulated in a state space framework. Uniformly valid exact solutions are obtained for hollow and solid circular cylinders subjected to extension, torsion, pressuring, shearing, electric loadings and a uniform temperature change.

This study develops a neural network system to monitor the safety of a bridge structure. A truck of constant mass is driven at constant speed through the target bridge. Then, the maximal and minimal values of the bridge elongations are processed by a monitoring system to evaluate the current condition of the bridge. The monitoring system is composed of parallel backpropagation neural networks. Each neural network monitors a part of the bridge. The neural networks are trained using simulation data. The numerical example shows that the monitoring system is effective in the damage detection of the bridge.