Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

Direct Numerical Simulation of Shock/Turbulent Boundary Layer Interaction in a Supersonic Compression Ramp

A direct numerical simulation of the shock/turbulent boundary layer interaction flow in a supersonic 24-degree compression ramp is conducted with the free stream Mach number 2.9. The blow-and-suction disturbance in the upstream wall boundary is used to trigger the transition. Both the mean wall pressure and the velocity profiles agree with those of the experimental data, which validates the simulation. The turbulent kinetic energy budget in the separation region is analyzed. Results show that the turbulent production term increases fast in the separation region, while the turbulent dissipation term reaches its peak in the near-wall region. The turbulent transport term contributes to the balance of the turbulent conduction and turbulent dissipation. Based on the analysis of instantaneous pressure in the downstream region of the mean shock and that in the separation bubble, the authors suggest that the low frequency oscillation of the shock is not caused by the upstream turbulent disturbance, but rather the instability of separation bubble.

Numerical study on onset of oscillatory thermocapillary flow in rectangular liquid pool

Thermocapillary flow in a rectangular liquid pool of large Prandtl fluid (Pr = 105.6) is numerically studied in microgravity. Oscillatory thermocapillary flow arises when the imposed temperature difference between the sidewalls exceeds a critical value. The fluctuations of the oscillatory flow, accompanied by the propagation of the hydrothermal wave from the cold sidewall to the hot one, are much smaller than the time-averaged velocity and temperature fields. The corresponding disturbance cells arise in the centre of the liquid pool initially, and extend to the whole region with the increasing imposed temperature difference. The present study reveals the different characteristics of the oscillatory themocapillary flow in the rectangular liquid pool as compared to the cases in other configurations.

Numerical Study on the Distribution Characteristics of TP in Taihu Lake, China

A depth-averaged two-dimensional hydrodynamics model and a two-dimensional water quality model, bases on the alternating direction implicit (ADI) method, is developed to study the distribution characteristics of total phosphorus (TP) in Taihu Lake. Wind stress, the pollution source from the inflow rivers, releasing rate of TP from the bottom sediment and water diversion from Yangtze River are the effecting factors of TP distribution. By using the model proposed in this paper, the concentration field of the total phosphorus was simulated, which leads to the conclusion that the flow field has a great influence on the spatial and temporal distribution of TP in the Taihu Lake.

Numerical investigation of particle velocity distributions in aeolian sand transport

Particle velocity distribution in a blowing sand cloud is a reflection of saltation movement of many particles. Numerical analysis is performed for particle velocity distribution with a discrete particle model. The probability distributions of resultant particle velocity in the impact-entrainment process, particle horizontal and vertical velocities at different heights and the vertical velocity of ascending particles are analyzed. The probability distributions of resultant impact and lift-off velocities of saltating particles can be expressed by a log-normal function, and that of impact angle comply with an exponential function. The probability distribution of particle horizontal and vertical velocities at different heights shows a typical single-peak pattern. In the lower part of saltation layer, the particle horizontal velocity distribution is positively skewed. Further analysis shows that the probability density function of the vertical velocity of ascending particles is similar to the right-hand part of a normal distribution function, and a general equation is acquired for the probability density function of non-dimensional vertical velocity of ascending particles which is independent of diameter of saltating particles, wind strength and height. These distributions in the present numerical analysis are consistent with reported experimental results. The present investigation is important for understanding the saltation state in wind-blown sand movement. (C) 2009 Elsevier B.V. All rights reserved.

Numerical Study on the Spatial and Temporal Characteristics of Water Quality in the Taihu Lake of China

Taihu Lake is the third largest fresh water lake in China. With the fast economic development, abundant industrial and agricultural waste water has been discharged into Taihu Lake, causing the eutrophication of the water quality, which greatly affected the water utility. In the past decades, the treatment of Taihu Lake has witnessed limited success. Therefore, it is practically and theoretically significant to study the eutrophication of Taihu Lake. This research has focused on the issue of water quality including the characteristics of spatial and temporal distributions, and the rules of nutrient diffusion in the Taihu lake area. Based on the monitoring data, the basis distribution characteristics of water quality in Taihu Lake are analyzed. Comparing Taihu Lake with other Lakes shows that one important reason for Taihu eutrophication is the long period of water retention. A transporting and diffusing model of Taihu nutrient is developed by combining with the hydrodynamics model. Using the model, the concentration field of the total phosphorus (TP) and the influence of wind-driven current are numerically investigated, which leads to the conclusion that the flow field has a great influence on the spatial and temporal distributions of TP in Taihu Lake. Furthermore, the effect for improving the water quality by the project of water diversion from the Yangtze River to Taihu Lake was analyzed by simulation. The results demonstrate that short-term water diversion cannot improve the water quality of the heavily-polluted Meiliang Bay and the western bank areas of Taihu Lake.

Numerical simulation of thermal loading produced by shaped high power laser onto engine parts

Recently a new method for simulating the thermal loading on pistons of diesel engines was reported. The spatially shaped high power laser is employed as the heat source, and some preliminary experimental and numerical work was carried out. In this paper, a further effort was made to extend this simulation method to some other important engine parts such as cylinder heads. The incident Gaussian beam was transformed into concentric multi-circular patterns of specific intensity distributions, with the aid of diffractive optical elements (DOEs). By incorporating the appropriate repetitive laser pulses, the designed transient temperature fields and thermal loadings in the engine parts could be simulated. Thermal-structural numerical models for pistons and cylinder heads were built to predict the transient temperature and thermal stress. The models were also employed to find the optimal intensity distributions of the transformed laser beam that could produce the target transient temperature fields. Comparison of experimental and numerical results demonstrated that this systematic approach is effective in simulating the thermal loading on the engine parts. (C) 2009 Elsevier Ltd. All rights reserved.

Numerical modeling of heat transfer and plasma flow in a low power arcjet thruster

Modeling studies are preformed to investigate the plasma and heat transfer characteristics of a low power argon arcjet thruster. Computed temperature, velocity, static pressure, and Mach number distribution in arcjet thruster under typical operating condition are presented in this paper. It shows that the performance data from numerical modeling results are basically consistent with the experimental measured values.

Experimental and numerical study on power consumptions in a double-tube-socket pneumatic conveying system

A new methodology is proposed in this paper to predict the lowest power consumption for a double-tube-socket (DTS) pneumatic conveying system. This methodology is established on both experimental work and numerical simulation. After parametric studies by numerical simulation, the desired conveying cases which have the lowest power consumption were obtained. Finally those cases were carried out in our experimental system. The measured power consumption was close to that predicted. In this paper the experimental work is discussed and the numerical simulation introduced. (c) 2010 Elsevier B.V. All rights reserved.

Numerical study on NOx/CO emissions in the diffusion flames of high-temperature off-gas of steelmaking converter

The combustion of high-temperature off-gas of steelmaking converter with periodical change of temperature and CO concentration always leads to CO and NOx over-standard emissions. In the paper, high-temperature off-gas combustion is simulated by adopting counterflow diffusion flame model, and some influencing factors of CO and NOx emissions are investigated by adopting a detailed chemistry GRI 3.0 mechanism. The emission index of NOx (EINOx) decreases 1.7–4.6% when air stoichiometric ratio (SR) increase from 0.6 to 1.4, and it dramatically increases with off-gas temperature at a given SR when the off-gas temperature is above 1500 K. High-concentration CO in off-gas can result in high NOx emissions, and NOx levels increase dramatically with CO concentration when off-gas temperature is above 1700 K. Both SR and off-gas temperature are important for the increase of CO burnout index (BICO) when SR is less than 1.0, but BICO increase about 1% when off-gas temperature increases from 1100 K to 1900 K at SR > 1.0. BICO increases with CO concentration in off-gas, and the influence of off-gas temperature on BICO is marginal. BICO increases with the relative humidity (RH) in air supplied, but it increases about 0.5% when RH is larger than 30%.

Numerical Simulation of a Low-power Hydrazine Arcjet Thruster

A modeling study is conducted to investigate the plasma flow and heat transfer characteristics of low-power (kW class) arc-heated thrusters (arcjets) with 2:1 hydrogen/nitrogen to simulate decomposed hydrazine as the propellant. The all-speed SIMPLE algorithm is employed to solve the governing equations, which take into account the effects of compressibility, the Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. Typical computed results about the temperature, velocity and Mach number distributions within arcjet thruster are presented for the case with arc current of 9 A and inlet stagnant pressure of 3.3×105 Pa to show the flow and heat transfer characteristics. It is found that the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appearing near the cathode tip, and the flow transition from the subsonic to supersonic regime occurs within the constrictor region. The effect of gas viscosity on the plasma flow within arcjet thruster is examined by an additional numerical test using artificially reduced values of gas viscosity. The test results show that the gas viscosity appreciably affects the plasma flow and the performance of the arcjet thruster for the cases with the hydrazine or hydrogen as the propellant. The integrated axial Lorentz force in the thruster nozzle is also calculated and compared with the thrust force of the arcjet thruster. It is found that the integrated axial Lorentz force is much smaller than the thrust force for the low-power arcjet thruster. Modeling results for the NASA 1-kW class arcjet thruster with simulated hydrazine as the propellant are found to be reasonably consistent with available experimental data.