Numerical simulation of a single fluid particle rising in non-Newtonian fluids

In this work, a level set method is developed for simulating the motion of a fluid particle rising in non-Newtonian fluids described by generalized Newtonian as well as viscoelastic model fluids. As the shear-thinning model we use a Carreau-Yasuda model, and the viscoelastic effect can be modeled with Oldroyd-B constitutive equations. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The level set method is implemented to compute the motion of a bubble in a Newtonian fluid as one of typical examples for validation, and the computational results are in good agreement with the reported experimental data．The level set method is also applied for simulating a Newtonian drop rising in Carreau-Yasuda and Oldroyd-B fluids．Numerical results including noticeably negative wake behind the drop and viscosity field are obtained, and compare satisfactorily with the known literature data.

An Improved Ce/Se Scheme For Multi-Material Elastic-Plastic Flows And Its Applications

In this paper, a new definition of SE and CE, which is based on the hexahedron mesh and simpler than Chang's original CE/SE method (the space-time Conservation Element and Solution Element method), is proposed and an improved CE/SE scheme is constructed. Furthermore, the improved CE/SE scheme is extended in order to solve the elastic-plastic flow problems. The hybrid particle level set method is used for tracing the interfaces of materials. Proper boundary conditions are presented in interface tracking. Two high-velocity impact problems are simulated numerically and the computational results are carefully compared with the experimental data, as well as the results from other literature and LS-DYNA software. The comparisons show that the computational scheme developed currently is clear in physical concept, easy to be implemented and high accurate and efficient for the problems considered. (C) 2008 Elsevier Ltd. All rights reserved.

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.

An Improved CE/SE Scheme for Numerical Simulation of Gaseous and Two-Phase Detonations

A new structure of solution elements and conservation elements based on rectangular mesh was pro- posed and an improved space-time conservation element and solution element (CE/SE) scheme with sec- ond-order accuracy was constructed. Furthermore, the application of improved CE/SE scheme was extended to detonation simulation. Three models were used for chemical reaction in gaseous detonation. And a two-fluid model was used for two-phase (gas–droplet) detonation. Shock reflections were simu- lated by the improved CE/SE scheme and the numerical results were compared with those obtained by other different numerical schemes. Gaseous and gas–droplet planar detonations were simulated and the numerical results were carefully compared with the experimental data and theoretical results based on C–J theory. Mach reflection of a cellular detonation was also simulated, and the numerical cellular pat- terns were compared with experimental ones. Comparisons show that the improved CE/SE scheme is clear in physical concept, easy to be implemented and high accurate for above-mentioned problems.

Microscopic Model Of Nano-Scale Particles Removal In High Pressure Co2-Based Solvents

A physical model is presented to describe the kinds of static forces responsible for adhesion of nano-scale copper metal particles to silicon surface with a fluid layer. To demonstrate the extent of particle cleaning, Received in revised form equilibrium separation distance (ESD) and net adhesion force (NAF) of a regulated metal particle with different radii (10-300 nm) on the silicon surface in CO2-based cleaning systems under different pressures were simulated. Generally, increasing the pressure of the cleaning system decreased the net adhesion force between spherical copper particle and silicon surface entrapped with medium. For CO2 + isopropanol cleaning system, the equilibrium separation distance exhibited a maximum at temperature 313.15 K in the Equilibrium separation distance regions of pressure space (1.84-8.02 MPa). When the dimension of copper particle was given, for example, High pressure 50 nm radius particles, the net adhesion force decreased and equilibrium separation distance increased with increased pressure in the CO2 + H2O cleaning system at temperature 348.15 K under 2.50-12.67 MPa pressure range. However, the net adhesion force and equilibrium separation distance both decreased with an increase in surfactant concentration at given pressure (27.6 or 27.5 MPa) and temperature (318 or 298 K) for CO2 + H2O with surfactant PFPE COO-NH4+ or DiF(8)-PO4-Na+. (C) 2008 Elsevier B.V. All rights reserved.

Influence Of Gravity On Structure Of Colloidal Crystal Using Simulated Microgravity

Liquid mixtures of water and deuterium oxide as the liquid phase, were used to match the density of charged colloidal particles. Kossel diffraction method was used to detect the crystal structures. The experiments under the density-matched (g=0) and unmatched (g=1) conditions are compared to examine the influence of gravity on the crystal structures formed by self-assembly of 110 nm (in diameter) polystyrene microspheres. The result shows that die gravity tends to make the lattice constants of colloidal crystals smaller at lower positions, which indicates that the effect of gravity should be taken into account in the study of the colloidal crystals.

Investigation on average void fraction for air/non-Newtonian power-law fluids two-phase flow in downward inclined pipes

The present work has been carried out to investigate on the average void fraction of gas/non-Newtonian fluids flow in downward inclined pipes. The influences of pipe inclination angle on the average void fraction were studied experimentally. A simple correlation, which incorporated the method of Vlachos et al. for gas/Newtonain fluid horizontal flow, the correction factor of Farooqi and Richardson and the pipe inclination angle, was proposed to predict the average void fraction of gas/non-Newtonian power-law stratified flow in downward inclined pipes. The correlation was based on 470 data points covering a wide range of flow rates for different systems at diverse angles. A good agreement was obtained between theory and data and the fitting results could describe the majority of the experimental data within ±20%.