MILU-CG method and the numerical study on the flow around a rotating circular cylinder

A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.

Vortex dislocation in wake-type flow caused by local spanwise nonuniformity

DNS of spatiotemporal evolution of a wake-type flow is performed. In the incoming flow, a local spanwise nonuniformity in momentum defect is initially imposed. Results show that the spanwise nonuniformity leads to a series of symmetric twist vortex dislocation in downstream of the flow. Vortex line variations and substantial transition of vorticity from spanwise to the streamwise and vertical directions clearly feature the generation of a vortex dislocation and the real vortex linking in the dislocation. Dynamical process and the mechanism responsible for the vortex dislocation are described.

Vortex dislocations in wake-type flow induced by spanwise disturbances

Vortex dislocations in wake-type flow induced by three types of spanwise disturbances superimposed on an upstream velocity profile are investigated by direct numerical simulations. Three distinct modes of vortex dislocations and flow transitions have been found. A local spanwise exponential decay disturbance leads to the appearance of a twisted chainlike mode of vortex dislocation. A stepped spanwise disturbance causes a streamwise periodic spotlike mode of vortex dislocation. A spanwise sinusoidal wavy disturbance with a moderate waviness causes a strong unsteadiness of wake behavior. This unsteadiness starts with a systematic periodic mode of vortex dislocation in the spanwise direction followed by the spanwise vortex shedding suppressed completely with increased time and the near wake becoming a steady shear flow. Characteristics of these modes of vortex dislocation and complex vortex linkages over the dislocation, as well as the corresponding dynamic processes related to the appearance of dislocations, are described by examining the variations of vortex lines and vorticity distribution. The nature of the vortex dislocation is demonstrated by the substantial vorticity modification of the spanwise vortex from the original spanwise direction to streamwise and vertical directions, accompanied by the appearance of noticeable vortex branching and complex vortex linking, all of which are produced at the locations with the biggest phase difference or with a frequency discontinuity between shedding cells. The effect of vortex dislocation on flow transition, either to an unsteady irregular vortex flow or suppression of the Kaacutermaacuten vortex shedding making the wake flow steady state, is analyzed. Distinct similarities are found in the mechanism and main flow phenomena between the present numerical results obtained in wake-type flows and the experimental-numerical results of cylinder wakes reported in previous studies.

Modeling study on the flow, heat transfer and energy conversion characteristics of low-power arc-heated hydrogen/nitrogen thrusters

A modeling study is conducted to investigate the effect of hydrogen content in propellants on the plasma flow, heat transfer and energy conversion characteristics of low-power (kW class) arc-heated hydrogen/nitrogen thrusters (arcjets). 1:0 (pure hydrogen), 3:1 (to simulate decomposed ammonia), 2:1 (to simulate decomposed hydrazine) and 0:1 (pure nitrogen) hydrogen/nitrogen mixtures are chosen as the propellants. Both the gas flow region inside the thruster nozzle and the anode-nozzle wall are included in the computational domain in order to better treat the conjugate heat transfer between the gas flow region and the solid wall region. The axial variations of the enthalpy flux, kinetic energy flux, directed kinetic-energy flux, and momentum flux, all normalized to the mass flow rate of the propellant, are used to investigate the energy conversion process inside the thruster nozzle. The modeling results show that the values of the arc voltage, the gas axial-velocity at the thruster exit, and the specific impulse of the arcjet thruster all increase with increasing hydrogen content in the propellant, but the gas temperature at the nitrogen thruster exit is significantly higher than that for other three propellants. The flow, heat transfer, and energy conversion processes taking place in the thruster nozzle have some common features for all the four propellants. The propellant is heated mainly in the near-cathode and constrictor region, accompanied with a rapid increase of the enthalpy flux, and after achieving its maximum value, the enthalpy flux decreases appreciably due to the conversion of gas internal energy into its kinetic energy in the divergent segment of the thruster nozzle. The kinetic energy flux, directed kinetic energy flux and momentum flux also increase at first due to the arc heating and the thermodynamic expansion, assume their maximum inside the nozzle and then decrease gradually as the propellant flows toward the thruster exit. It is found that a large energy loss (31-52%) occurs in the thruster nozzle due to the heat transfer to the nozzle wall and too long nozzle is not necessary. Modeling results for the NASA 1-kW class arcjet thruster with hydrogen or decomposed hydrazine as the propellant are found to compare favorably with available experimental data.

Effect of variations in acid properties of HZSM-5 on the coking behavior and reaction stability in butene aromatization

In order to investigate the effect of acid properties on the coke behavior and stability of butene aromatization, we prepared the AHZSM-5 samples with various acid properties by the methods of hydrothernial treatment and K addition. The reaction of butene aromatization was carried out at 350 degrees C and 0.5 MPa in a continuous flow fixed bed. The characterization of the fresh/coked catalysts with NH3-TPD, N-2 adsorption-desorption measurement, and TG techniques has shown that a large amount of acid sites (high acid density) of the AHZMS-5 catalyst can cause a large quantity of coke deposit and serious channel blockage, and so result in a rapid loss of aromatization activity. On the contrary, after a great reduction in strong acid sites of AHZSM-5 catalyst resulting from some K-modification, the presence of only many weak acid sites also could not lessen the formation of coke nor improve the reaction stability of butene aromatization. Interestingly, the simultaneous reduction in the strong and weak acid sites to a desirable level by hydrothermal treating the AHZSM-5 catalyst at a proper temperature can effectively suppress the coke formation and channel blockage, and thus improve its olefin aromatization stability. (c) 2005 Elsevier B.V. All rights reserved.