Direct Numerical Simulation Of Three-Dimensional Richtmyer-Meshkov Instability

Direct numerical simulation (DNS) is used to study flow characteristics after interaction of a planar shock with a spherical media interface in each side of which the density is different. This interfacial instability is known as the Richtmyer-Meshkov (R-M) instability. The compressible Navier-Stoke equations are discretized with group velocity control (GVC) modified fourth order accurate compact difference scheme. Three-dimensional numerical simulations are performed for R-M instability installed passing a shock through a spherical interface. Based on numerical results the characteristics of 3D R-M instability are analysed. The evaluation for distortion of the interface, the deformation of the incident shock wave and effects of refraction, reflection and diffraction are presented. The effects of the interfacial instability on produced vorticity and mixing is discussed.

Study of drag reduction by gas injection for power-law fluid flow in horizontal stratified and slug flow regimes

In this work, the drag reduction by gas injection for power-law fluid flow in stratified and slug flow regimes has been studied. Experimentswere conducted to measure the pressure gradient within air/CMC solutions in a horizontal Plexiglas pipe that had a diameter of 50mm and a length of 30 m. The drag reduction ratio in stratified flow regime was predicted using the two-fluid model. The results showed that the drag reduction should occur over the large range of the liquid holdup when the flow behaviour index remained at the low value. Furthermore, for turbulent gas-laminar liquid stratified flow, the drag reduction by gas injection for Newtonian fluid was more effective than that for shear-shinning fluid, when the dimensionless liquid height remained in the area of high value. The pressure gradient model for a gas/Newtonian liquid slug flow was extended to liquids possessing the Ostwald–de Waele power law model. The proposed model was validated against 340 experimental data point over a wide range of operating conditions, fluid characteristics and pipe diameters. The dimensionless pressure drop predicted was well inside the 20% deviation region for most of the experimental data. These results substantiated the general validity of the model presented for gas/non-Newtonian two-phase slug flows.

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.

Direct dynamical test for deterministic chaos

We present a direct and dynamical method to distinguish low-dimensional deterministic chaos from noise. We define a series of time-dependent curves which are closely related to the largest Lyapunov exponent. For a chaotic time series, there exists an envelope to the time-dependent curves, while for a white noise or a noise with the same power spectrum as that of a chaotic time series, the envelope cannot be defined. When a noise is added to a chaotic time series, the envelope is eventually destroyed with the increasing of the amplitude of the noise.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

Optimization of 2-stage gas gun for fusion refueling pellet injection with consideration of real effects

Analytic expression of pellet acceleration by constant base pressure with consideration of gas-wall friction, heat transfer and viscous dissipation that important for high speed injection is obtained. The process of compression stage is formulated by a set of governing equations and can be numerically integrated. Excellent confirmation with experiments is obtained and the ways to optimum match the compression stage with the launch stage are suggested.

Investigation of Vaporized Kerosene Injection in a Supersonic Model Combustor

Characteristics of vaporized aviation kerosene injection in a supersonic model combustor were preliminarily investigated. The electrically storage type heater has a volume capacity of heating kerosene of 0.8 kg up to 670 K at a pressure of 5.5 Mpa. The temperature to cause pressurized kerosene jet being fully vaporized in Quiescent atmosphere was found to be 550 K at 4 Mpa however the pressurized hot kerosene remains in liquid state within the tube. The correspondent jet spray in Mach 2.5 vitiated air cross-flow were visualized by using stop schlieren photograph.It was found the penetration depth of the hot pressurized kerosene jet is approximately same with the temperature varied from 290 K to 550 k. at pressure of 4 Mpa. This results showed that the atomization process of hot kerosene jet spray in supersonic combustor could be bypassed and directly transferred to be gas state at temperature 550 K and pressure of 4 Mpa.

Numeric Calculation for Transient Thermal Stress Field in Wall-Shape Metal Part During Laser Direct Forming

In order to investigate the transient thermal stress field in wall-shape metal part during laser direct forming, a FEM model basing on ANSYS is established, and its algorithm is also dealt with. Calculation results show that while the wall-shape metal part is being deposited, in X direction, the thermal stress in the top layer of the wall-shape metal part is tensile stress and in the inner of the wall-shape metal part is compressive stress. The reason causing above-mentioned thermal stress status in the wall-shape metal part is illustrated, and the influence of the time and the processing parameters on the thermal stress field in wall-shape metal part is also studied. The calculation results are consistent with experimental results in tendency.

Performance of a Supersonic Model Combustor using Vaporized Kerosene Injection

Characteristics of supersonic combustion by injecting kerosene vapor into a Mach 2.5 crossflow at various preheat temperatures and pressures were investigated experimentally. A two-stage heating system has been designed and tested, which can prepare heated kerosene of 0.8 kg up to 820 K at pressure of 5.5 Mpa with minimum/negligible fuel coking. In order to simulate the thermophysical properties of kerosene over a wide range of thermodynamic conditions, a three-component surrogate that matches the compound class of the parent fuel was employed. The flow rate of kerosene vapor was calibrated using a sonic nozzle. Computed flow rates using the surrogate fuel are in agreement with the experimental data. Kerosene jets at various preheat temperatures injecting into both quiescent environment and Mach 2.5 crossflow were visualized. It was found that at injection pressure of 4 Mpa and preheat temperature of 550 K the kerosene jet was completely in vapor phase, while keeping almost the same penetration depth as compared to the liquid kerosene injection. Supersonic combustion tests were also carried out to compare the combustor performance for the cases of vaporized kerosene injection, liquid kerosene injection, and effervescent atomization with hydrogen barbotage, under the similar stagnation conditions. Experimental results demonstrated that the use of vaporized kerosene injection leads to better combustor performance. Further parametric study on vaporized kerosene injection in a supersonic model combustor is needed to assess the combustion efficiency as well as to identify the controlling mechanism for the overall combustion enhancement.

Investigation of Vaporized Kerosene Injection and Combustion in a Supersonic Model Combustor

Injection and combustion of vaporized kerosene was experimentally investigated in a Mach 2.5 model combustor at various fuel temperatures and injection pressures. A unique kerosene heating and delivery system, which can prepare heated kerosene up to 820 K at a pressure of 5.5 MPa with negligible fuel coking, was developed. A three-species surrogate was employed to simulate the thermophysical properties of kerosene. The calculated thermophysical properties of surrogate provided insight into the fuel flow control in experiments. Kerosene jet structures at various preheat temperatures injecting into both quiescent environment and a Mach 2.5 crossflow were characterized. It was shown that the use ofvaporized kerosene injection holds the potential of enhancing fuel-air mixing and promoting overall burning. Supersonic combustion tests further confirmed the preceding conjecture by comparing the combustor performances of supercritical kerosene with those of liquid kerosene and effervescent atomization with hydrogen barbotage. Under the similar flow conditions and overall kerosene equivalence ratios, experimental results illustrated that the combustion efficiency of supercritical kerosene increased approximately 10-15% over that of liquid kerosene, which was comparable to that of effervescent atomization.

Direct initiation of detonation with ignition tube

Experiments of direct initiation of hydrogen-oxygen by means of a hot turbulent jet were made. Results indicate that the length of ignition tube is the dominant factor in determining the ignition capability of hot turbulent jet, and that the ignition capability of turbulence jet increases with the length of ignition tube. Because this ignition capability can meet the demands of a gas-detonation-driver shock tunnel and it doesn't require additional facilities, the hot turbulent jet initiation method can be applied to large hydrogen-oxygen detonation-driver shock tunnels. Influences of obstacles on the ignition capability were also studied. It was found that the presence of obstacles weakens the ignition capability of a hot turbulent jet.

Development of a Cell Size Relaxed Scheme for the Direct Simulation Monte Carlo Method

The conventional direct simulation Monte Carlo (DSMC) method has a strong restriction on the cell size because simulated particles are selected randomly within the cell for collisions. Cells with size larger than the molecular mean free path are generally not allowed in correct DSMC simulations. However, the cell-size induced numerical error can be controlled if the gradients of flow properties are properly involved during collisions. In this study, a large cell DSMC scheme is proposed to relax the cell size restriction. The scheme is applied to simulate several test problems and promising results are obtained even when the cell size is greater than 10 mean free paths of gas molecules. However, it is still necessary, of course, that the cell size be small with respect to the flow field structures that must be resolved.

Improved sampling techniques for the direct simulation Monte Carlo method

The direct simulation Monte Carlo (DSMC) method is a widely used approach for flow simulations having rarefied or nonequilibrium effects. It involves heavily to sample instantaneous values from prescribed distributions using random numbers. In this note, we briefly review the sampling techniques typically employed in the DSMC method and present two techniques to speedup related sampling processes. One technique is very efficient for sampling geometric locations of new particles and the other is useful for the Larsen-Borgnakke energy distribution.

Direct numerical simulation of hypersonic boundary layer transition over a blunt cone with a small angle of attack

The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [Li et al., AIAA J. 46, 2899(2008)], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 0° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves VLFWs� are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.The direct numerical simulation of boundary layer transition over a 5° half-cone-angle blunt cone is performed. The free-stream Mach number is 6 and the angle of attack is 1°. Random wall blow-and-suction perturbations are used to trigger the transition. Different from the authors’ previous work [ Li et al., AIAA J. 46, 2899 (2008) ], the whole boundary layer flow over the cone is simulated (while in the author’s previous work, only two 45° regions around the leeward and the windward sections are simulated). The transition location on the cone surface is determined through the rapid increase in skin fraction coefficient (Cf). The transition line on the cone surface shows a nonmonotonic curve and the transition is delayed in the range of 20° ≤ θ ≤ 30° (θ = 0° is the leeward section). The mechanism of the delayed transition is studied by using joint frequency spectrum analysis and linear stability theory (LST). It is shown that the growth rates of unstable waves of the second mode are suppressed in the range of 20° ≤ θ ≤ 30°, which leads to the delayed transition location. Very low frequency waves (VLFWs) are found in the time series recorded just before the transition location, and the periodic times of VLFWs are about one order larger than those of ordinary Mack second mode waves. Band-pass filter is used to analyze the low frequency waves, and they are deemed as the effect of large scale nonlinear perturbations triggered by LST waves when they are strong enough.

Two-Phase Flow In Anode Flow Field Of A Small Direct Methanol Fuel Cell In Different Gravities

An in-situ visualization of two-phase flow inside anode flow bed of a small liquid fed direct methanol fuel cells in normal and reduced gravity has been conducted in a drop tower. The anode flow bed consists of 11 parallel straight channels. The length, width and depth of single channel, which had rectangular cross section, are 48.0, 2.5 and 2.0 mm, respectively. The rib width was 2.0 mm. The experimental results indicated that when the fuel cell orientation is vertical, two-phase flow pattern in anode channels can evolve from bubbly flow in normal gravity into slug flow in microgravity. The size of bubbles in the reduced gravity is also bigger. In microgravity, the bubbles rising speed in vertical channels is obviously slower than that in normal gravity. When the fuel cell orientation is horizontal, the slug flow in the reduced gravity has almost the same characteristic with that in normal gravity. It implies that the effect of gravity on two-phase flow is small and the bubbles removal is governed by viscous drag. When the gas slugs or gas columns occupy channels, the performance of liquid fed direct methanol fuel cells is failing rapidly. It infers that in long-term microgravity, flow bed and operating condition should be optimized to avoid concentration polarization of fuel cells.