Numerical Analysis on Intermittent Flow in a Longitudinal Section of Bingham Fluid along a Slope

The transition process of intermittent flow in a longitudinal section of Bingham fluid from initial distribution to fully developed state was numerically investigated in this paper. The influences of slope dimensionless runoff Q* and viscosity μ0* on the dimensionless surge speed U* were also presented in a wide range of parameters. By one typical example, the intermittent flow possessed wave characteristics and showed a supercritical flow conformation for a fully developed flow. The distributions of gravity and bed drag along the flow path and the velocity distribution of flow field were also analyzed.

Influence of insoluble surfactant on the deformation and breakup of a bubble or thread in a viscous fluid

The influence of surfactant on the breakup of a prestretched bubble in a quiescent viscous surrounding is studied by a combination of direct numerical simulation and the solution of a long-wave asymptotic model. The direct numerical simulations describe the evolution toward breakup of an inviscid bubble, while the effects of small but non-zero interior viscosity are readily included in the long-wave model for a fluid thread in the Stokes flow limit. The direct numerical simulations use a specific but realizable and representative initial bubble shape to compare the evolution toward breakup of a clean or surfactant-free bubble and a bubble that is coated with insoluble surfactant. A distinguishing feature of the evolution in the presence of surfactant is the interruption of bubble breakup by formation of a slender quasi-steady thread of the interior fluid. This forms because the decrease in surface area causes a decrease in the surface tension and capillary pressure, until at a small but non-zero radius, equilibrium occurs between the capillary pressure and interior fluid pressure. The long-wave asymptotic model, for a thread with periodic boundary conditions, explains the principal mechanism of the slender thread's formation and confirms, for example, the relatively minor role played by the Marangoni stress. The large-time evolution of the slender thread and the precise location of its breakup are, however, influenced by effects such as the Marangoni stress and surface diffusion of surfactant. © 2008 Cambridge University Press.

Turbulence modeling for flows around convex features giving rapid eddy distortion

Reynolds averaged Navier-Stokes model performances in the stagnation and wake regions for turbulent flows with relatively large Lagrangian length scales (generally larger than the scale of geometrical features) approaching small cylinders (both square and circular) is explored. The effective cylinder (or wire) diameter based Reynolds number, ReW ≤ 2.5 × 103. The following turbulence models are considered: a mixing-length; standard Spalart and Allmaras (SA) and streamline curvature (and rotation) corrected SA (SARC); Secundov's νt-92; Secundov et al.'s two equation νt-L; Wolfshtein's k-l model; the Explicit Algebraic Stress Model (EASM) of Abid et al.; the cubic model of Craft et al.; various linear k-ε models including those with wall distance based damping functions; Menter SST, k-ω and Spalding's LVEL model. The use of differential equation distance functions (Poisson and Hamilton-Jacobi equation based) for palliative turbulence modeling purposes is explored. The performance of SA with these distance functions is also considered in the sharp convex geometry region of an airfoil trailing edge. For the cylinder, with ReW ≈ 2.5 × 103 the mixing length and k-l models give strong turbulence production in the wake region. However, in agreement with eddy viscosity estimates, the LVEL and Secundov νt-92 models show relatively little cylinder influence on turbulence. On the other hand, two equation models (as does the one equation SA) suggest the cylinder gives a strong turbulence deficit in the wake region. Also, for SA, an order or magnitude cylinder diameter decrease from ReW = 2500 to 250 surprisingly strengthens the cylinder's disruptive influence. Importantly, results for ReW ≪ 250 are virtually identical to those for ReW = 250 i.e. no matter how small the cylinder/wire its influence does not, as it should, vanish. Similar tests for the Launder-Sharma k-ε, Menter SST and k-ω show, in accordance with physical reality, the cylinder's influence diminishing albeit slowly with size. Results suggest distance functions palliate the SA model's erroneous trait and improve its predictive performance in wire wake regions. Also, results suggest that, along the stagnation line, such functions improve the SA, mixing length, k-l and LVEL results. For the airfoil, with SA, the larger Poisson distance function increases the wake region turbulence levels by just under 5%. © 2007 Elsevier Inc. All rights reserved.

Large eddy simulations for fan-OGV broadband noise prediction

The work in this paper forms part of a project on the use of large eddy simulation (LES) for broadband rotor-stator interaction noise prediction. Here we focus on LES of the flow field near a fan blade trailing edge. The first part of the paper aims to evaluate LES suitability for predicting the near-field velocity field for a blunt NACA-0012 airfoil at moderate Reynolds numbers (2× 10 5 and 4× 10 5). Preliminary computations of turbulent mean and root-mean-square velocities, as well as energy spectra at the trailing edge, are compared with those from a recent experiment.1 The second part of the paper describes preliminary progress on an LES calculation of the fan wakes on a fan rig. 2 The CFD code uses a mixed element unstructured mesh with a median dual control volume. A wall-adapting local eddy-viscosity sub-grid scale model is employed. A very small amount of numerical dissipation is added in the numerical scheme to keep the compressible solver stable. Further results for the fan turbulentmean and RMS velocity, and especially the aeroacoustics field will be presented at a later stage. Copyright © 2008 by Qinling LI, Nigel Peake & Mark Savill.

The interaction between poly(vinylpyrrolidone) and reversed micelles of water/AOT/n-heptane

The interactions between poly(vinylpyrrolidone) (PVP) and the reversed micelles composed of water, AOT, and n-heptane are investigated with the aid of phase diagram, measurements of conductivity and viscosity, Fourier transform infrared (FTIR) spectrum, and dynamic light scattering (DLS). The phase diagrams of water/AOT/heptane in the presence of and absence of PVP are given. The conductivity of the water/AOT/heptane reversed micelle without PVP initially increases and then decreases with the increase of water content, ω0 (the molar ratio of water to AOT), while the plots of conductivity (K) versus ω0 of the reversed micelle in the presence of PVP depend on the PVP concentrations. The plot of K versus ω0 with 2.0%wt PVP is similar to that without PVP. Only the ω0,max (the water content that the maximum conductivity corresponds to) is larger than that without PVP. Nevertheless, the conductivity of the reversed micelle containing more than 4%wt PVP always rises with the increase of the water content in the measured range. The DLS results indicate that the hydrodynamic radius (Rh) in the presence and absence of PVP rises with the increase of ω0. The plots with PVP and without PVP have almost the same value when ω0<17; and after that, it quickly increases with the increase of ω0. It is interesting to find that there is almost no effect of the PVP concentration on the viscosity and Rh of the reversed micelle at ω0 = 15. The FTIR results suggest that the contents of SO3--bound water and Na+-bound water both decrease with PVP added, while the content of the bulky-like water increases. However, the trapped water in the hydrophobic chain of the surfactant is nearly unaffected by PVP. It is also found from the FTIR that the carbonyl group stretching vibration of AOT is fitted into two sub-peaks, which center at 1740 and 1729 cm-1, corresponding to the trans and cis conformations of AOT, respectively.

A Microtube Viscometer with a Thermostat

The viscometer presented in this paper is suit-able for measuring the viscosity of liquids in micro-litre quantities. It consists of a micro-flow experimental system with a thermostat. Using the measurements of the flow rates and pressure drops of a liquid passing through a microtube, the liquid's viscosity can be calculated from on Hagen-Poiseuille theory. After calibration, the viscometer was used to measure viscosities of deionized water and ethyl alcohol at temperatures ranging from 0 to 40 "C. For both test liquids, the relative deviation of the measured values from those quoted in the literature (obtained using other viscometers) was less than 2.6o/o. The relative uncertainty of the experimental system was reduced to +-l.8% using the relative measuring method. Due to the micro-scale of the test section, only a micro-litre quantity of liquid is needed for a test, this is a potential advantage for measurement of bio-liquid viscosities.

Physical parameters affecting sonoluminescence: A self-consistent hydrodynamic study

We studied the dependence of thermodynamic variables in a sonoluminescing ~SL! bubble on various physical factors, which include viscosity, thermal conductivity, surface tension, the equation of state of the gas inside the bubble, as well as the compressibility of the surrounding liquid. The numerical solutions show that the existence of shock waves in the SL parameter regime is very sensitive to these factors. Furthermore, we show that even without shock waves, the reflection of continuous compressional waves at the bubble center can produce the high temperature and picosecond time scale light pulse of the SL bubble, which implies that SL may not necessarily be due to shock waves.

Numerical Study Of Heat Transfer Mechanism In Turbulent Supercritical Co2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Molecular Dynamics Simulation Of Adhesion Forces In A Dipalmitoylphosphatidylcholine Membrane

Adhesion forces of Dipalmitoylphosphatidylcholine ( DPPC) membrane in the gel phase are investigated by molecular dynamics ( MD) simulation. In the simulations, individual DPPC molecules are pulled out of DPPC membranes with different rates and we get the maximum adhesion forces of DPPC membrane. We find that the maximum adhesion forces increase with pull rate, from about 400 to 700 pN when pull rates are from 0.001 to 0.03 nm/ps. We analyze the relationship between pull rate and adhesion forces of different origins using Brownian dynamics and notice that viscosity of solvent plays an important role in adhesion forces. Then we simulate the motion of a single DPPC molecule in solvent and it elucidates that the maximum drag force is almost linear with respect to the pull rate. We use Stokes' relation to describe the motion of a single DPPC molecule and deduce the effective length of a DPPC molecule. Conformational analyses indicate that the free energy variation of a DPPC molecule inside and outside of the DPPC membrane is an essential part of adhesion energy.

On-Board Experimental Study Of Bubble Thermocapillary Migration In A Recoverable Satellite

An on-board space experiment of bubble thermocapillary migration was performed in the Chinese 22nd recoverable satellite in 2005. Silicone oil of nominal viscosity 5cSt was used as the continuous phase in the experiment. Air bubbles were injected into the liquid in the same direction as the constant temperature gradient in the liquid. The velocities of bubbles were obtained by recording the paths of the bubbles. The results indicate that the scaled velocity of bubbles decreases with an increase of the Marangoni number extended to 9288, which agrees with the results of previous space experiments and numerical simulation. In addition, the interaction between two bubbles was also observed in the space experiment. The trajectories and the velocities of the bubbles were obtained. The two-bubble experiment results are also consistent with the theoretical analysis.

Experimental Investigation On The Holdup Distribution Of Oil-Water Two-Phase Flow In Horizontal Parallel Tubes

An experimental investigation was conducted to study the holdup distribution of oil and water two-phase flow in two parallel tubes with unequal tube diameter. Tests were performed using white oil (of viscosity 52 mPa s and density 860 kg/m(3)) and tap water as liquid phases at room temperature and atmospheric outlet pressure. Measurements were taken of water flow rates from 0.5 to 12.5 m(3)/h and input oil volume fractions from 3 to 94 %. Results showed that there were different flow pattern maps between the run and bypass tubes when oil-water two-phase flow is found in the parallel tubes. At low input fluid flow rates, a large deviation could be found on the average oil holdup between the bypass and the run tubes. However, with increased input oil fraction at constant water flow rate, the holdup at the bypass tube became close to that at the run tube. Furthermore, experimental data showed that there was no significant variation in flow pattern and holdup between the run and main tubes. In order to calculate the holdup in the form of segregated flow, the drift flux model has been used here.

Sensitivity Analysis Of The Dimensionless Parameters In Scaling A Polymer Flooding Reservoir

A set of scaling criteria of a polymer flooding reservoir is derived from the governing equations, which involve gravity and capillary force, compressibility of water, oil, and rock, non-Newtonian behavior of the polymer solution, absorption, dispersion, and diffusion, etc. A numerical approach to quantify the dominance degree of each dimensionless parameter is proposed. With this approach, the sensitivity factor of each dimensionless parameter is evaluated. The results show that in polymer flooding, the order of the sensitivity factor ranges from 10(-5) to 10(0) and the dominant dimensionless parameters are generally the ratio of the oil permeability under the condition of the irreducible water saturation to water permeability under the condition of residual oil saturation, density, and viscosity ratios between water and oil, the reduced initial oleic phase saturation and the shear rate exponent of the polymer solution. It is also revealed that the dominant dimensionless parameters may be different from case to case. The effect of some physical variables, such as oil viscosity, injection rate, and permeability, on the dominance degree of the dimensionless parameters is analyzed and the dominant ones are determined for different cases.

Effects Of Subgrid-Scale Modeling On Lagrangian Statistics In Large-Eddy Simulation

The application of large-eddy simulation (LES) to turbulent transport processes requires accurate prediction of the Lagrangian statistics of flow fields. However, in most existing SGS models, no explicit consideration is given to Lagrangian statistics. In this paper, we focus on the effects of SGS modeling on Lagrangian statistics in LES ranging from statistics determining single-particle dispersion to those of pair dispersion and multiparticle dispersion. Lagrangian statistics in homogeneous isotropic turbulence are extracted from direct numerical simulation (DNS) and the LES with a spectral eddy-viscosity model. For the case of longtime single-particle dispersion, it is shown that, compared to DNS, LES overpredicts the time scale of the Lagrangian velocity correlation but underpredicts the Lagrangian velocity fluctuation. These two effects tend to cancel one another leading to an accurate prediction of the longtime turbulent dispersion coefficient. Unlike the single-particle dispersion, LES tends to underestimate significantly the rate of relative dispersion of particle pairs and multiple-particles, when initial separation distances are less than the minimum resolved scale due to the lack of subgrid fluctuations. The overprediction of LES on the time scale of the Lagrangian velocity correlation is further confirmed by a theoretical analysis using a turbulence closure theory.

Derivation Of A Nonlinear Reynolds Stress Model Using Renormalization Group Analysis And Two-Scale Expansion Technique

Adopting Yoshizawa's two-scale expansion technique, the fluctuating field is expanded around the isotropic field. The renormalization group method is applied for calculating the covariance of the fluctuating field at the lower order expansion. A nonlinear Reynolds stress model is derived and the turbulent constants inside are evaluated analytically. Compared with the two-scale direct interaction approximation analysis for turbulent shear flows proposed by Yoshizawa, the calculation is much more simple. The analytical model presented here is close to the Speziale model, which is widely applied in the numerical simulations for the complex turbulent flows.

Effect of Fe content on the thermal stability and dynamic mechanical behavior of NdAlNiCu bulk metallic glasses

The dependence of microstructure and thermal stability on Fe content of bulk Nd60Al10Ni10Cu20-xFex (0 less than or equal to x less than or equal to 20) metallic glasses is investigated by means of differential scanning calorimetry (DSC), X-ray diffraction (XRD) and high-resolution transmission electron micrograph (HRTEM). All samples exhibit typical amorphous feature under the detect limit of XRD, however, HRTEM results show that the microstructure of Nd60Al10Ni10Cu20-xFex alloys changes from a homogeneous amorphous phase to a composite structure consisting of clusters dispersed in amorphous matrix by increasing Fe content. Dynamic mechanical properties of these alloys with controllable microstructure are studied, expressed via storage modulus, the loss modulus and the mechanical damping. The results reveal that the storage modulus of the alloy without Fe added shows a distinct decrease due to the main a relaxation. This decrease weakens and begins at a higher temperature with increasing Fe content. The mechanism of the effect of Fe addition on the microstructure and thermal stability in this system is discussed in terms of thermodynamics viewpoints. (C) 2004 Elsevier B.V. All rights reserved.