Large-Eddy-Simulation of 3-Dimensional Rayleigh-Taylor Instability in Incompressible Fluids

The 3-dimensiqnal incompressible Rayleigh-Taylor instability is numerically studied through the large-eddy-simulation (LES) approach based on the passive scalar transport model. Both the instantaneous velocity and the passive scalar fields excited by sinu

Probe for Measurements of Density/Conductivity in Flows of Conducting Fluids

A probe utilizing the bipolar pulse method to measure the density of a conducting fluid has been developed. The probe is specially designed such that the concentration of a stream tube can be sampled continuously. The density was determined indirectly from the measurement of solution conductivity. The probe was calibrated using standard NaCl solutions of varying molarity and was able to rapidly determine the density of a fluid with continuously varying conductance. Measurements of the conductivity profiles, corresponding density profiles, and their fluctuation levels are demonstrated in a channel flow with an electrolyte injected from a slot in one wall.

Rossby Waves With Linear Topography In Barotropic Fluids

Rossby waves are the most important waves in the atmosphere and ocean, and are parts of a large-scale system in fluid. The theory and observation show that, they satisfy quasi-geostrophic and quasi-static equilibrium approximations. In this paper, solitary Rossby waves induced by linear topography in barotropic fluids with a shear flow are studied. In order to simplify the problem, the topography is taken as a linear function of latitude variable y, then employing a weakly nonlinear method and a perturbation method, a KdV (Korteweg-de Vries) equation describing evolution of the amplitude of solitary Rossby waves induced by linear topography is derived. The results show that the variation of linear topography can induce the solitary Rossby waves in barotropic fluids with a shear flow, and extend the classical geophysical theory of fluid dynamics.

Vertically forced surface wave in weakly viscous fluids bounded in a circular cylindrical vessel

Two-time scale perturbation expansions were developed in weakly viscous fluids to investigate surface wave motions by linearizing the Navier-Stokes equation in a circular cylindrical vessel which is subject to a vertical oscillation. The fluid field was divided into an outer potential flow region and an inner boundary layer region. A linear amplitude equation of slowly varying complex amplitude, which incorporates a damping term and external excitation, was derived for the weakly viscid fluids. The condition for the appearance of stable surface waves was obtained and the critical curve was determined. In addition, an analytical expression for the damping coefficient was determined and the relationship between damping and other related parameters (such as viscosity, forced amplitude, forced frequency and the depth of fluid, etc.) was presented. Finally, the influence both of the surface tension and the weak viscosity on the mode formation was described by comparing theoretical and experimental results. The results show that when the forcing frequency is low, the viscosity of the fluid is prominent for the mode selection. However, when the forcing frequency is high, the surface tension of the fluid is prominent.

Proper orthogonal decomposition of oscillatory Marangoni flow in half-zone liquid bridges of low-Pr fluids

Proper orthogonal decomposition (POD) using method of snapshots was performed on three different types of oscillatory Marangoni flows in half-zone liquid bridges of low-Pr fluid (Pr = 0.01). For each oscillation type, a series of characteristic modes (eigenfunctions) have been extracted from the velocity and temperature disturbances, and the POD provided spatial structures of the eigenfunctions, their oscillation frequencies, amplitudes, and phase shifts between them. The present analyses revealed the common features of the characteristic modes for different oscillation modes: four major velocity eigenfunctions captured more than 99% of the velocity fluctuation energy form two pairs, one of which is the most energetic. Different from the velocity disturbance, one of the major temperature eigenfunctions makes the dominant contribution to the temperature fluctuation energy. On the other hand, within the most energetic velocity eigenfuction pair, the two eigenfunctions have similar spatial structures and were tightly coupled to oscillate with the same frequency, and it was determined that the spatial structures and phase shifts of the eigenfunctions produced the different oscillatory disturbances. The interaction of other major modes only enriches the secondary spatio-temporal structures of the oscillatory disturbances. Moreover, the present analyses imply that the oscillatory disturbance, which is hydrodynamic in nature, primarily originates from the interior of the liquid bridge. (C) 2007 Elsevier B.V. All rights reserved.

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.

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%.

Electronic structure of paramagnetic In1-xMnx as nanowires

The electronic structure, spin splitting energies, and g factors of paramagnetic In1-xMnxAs nanowires under magnetic and electric fields are investigated theoretically including the sp-d exchange interaction between the carriers and the magnetic ion. We find that the effective g factor changes dramatically with the magnetic field. The spin splitting due to the sp-d exchange interaction counteracts the Zeeman spin splitting. The effective g factor can be tuned to zero by the external magnetic field. There is also spin splitting under an electric field due to the Rashba spin-orbit coupling which is a relativistic effect. The spin-degenerated bands split at nonzero k(z) (k(z) is the wave vector in the wire direction), and the spin-splitting bands cross at k(z) = 0, whose k(z)-positive part and negative part are symmetrical. A proper magnetic field makes the k(z)-positive part and negative part of the bands asymmetrical, and the bands cross at nonzero k(z). In the absence of magnetic field, the electron Rashba coefficient increases almost linearly with the electric field, while the hole Rashba coefficient increases at first and then decreases as the electric field increases. The hole Rashba coefficient can be tuned to zero by the electric field.

Phase equilibria and plate-fluid interfacial tensions for associating hard sphere fluids confined in slit pores

The excess Helmholtz free energy functional for associating hard sphere fluid is formulated by using a modified fundamental measure theory [Y. X. Yu and J. Z. Wu, J. Chem. Phys. 117, 10156 (2002)]. Within the framework of density functional theory, the thermodynamic properties including phase equilibria for both molecules and monomers, equilibrium plate-fluid interfacial tensions and isotherms of excess adsorption, average molecule density, average monomer density, and plate-fluid interfacial tension for four-site associating hard sphere fluids confined in slit pores are investigated. The phase equilibria inside the hard slit pores and attractive slit pores are determined according to the requirement that temperature, chemical potential, and grand potential in coexistence phases should be equal and the plate-fluid interfacial tensions at equilibrium states are predicted consequently. The influences of association energy, fluid-solid interaction, and pore width on phase equilibria and equilibrium plate-fluid interfacial tensions are discussed.

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.

Effect of volume ratio on thermocapillary flow in liquid bridges of high-Prandtl-number fluids

In present study, the transition of thermocapillary convection from the axisymmetric stationary flow to oscillatory flow in liquid bridges of 5cst silicon oil (aspect ratio 1.0 and 1.6) is investigated in microgravity conditions by the linear instability analysis. The corresponding marginal instability boundary is closely related to the gas/liquid configuration of the liquid bridge noted as volume ratio. With the increasing volume ratio, the marginal instability boundary consists of the increasing branch and the decreasing branch. A gap region exists between the branches where the critical Marangoni number of the corresponding axisymmetric stationary flow increases drastically. Particularly, a unique axisymmetric oscillatory flow (the critical azimuthal wave number is m=0) in the gap region is reported for the liquid bridge of aspect ratio 1.6. Moreover, the energy transfer between the basic state and the disturbance fields of the thermocapillary convection is analyzed at the corresponding critical Marangoni number, which reveals different major sources of the energy transfer for the development of the disturbances in regimes of the increasing branch, the gap region and the decreasing branch, respectively.

A novel structure analysis method for lanthanide complex - None paramagnetic shift tri-lanthanide mixture method

A new structure analysis method for lanthanide complexes was proposed, that is, none paramagnetic shift tri-lanthanide mixture method, It was found that the paramagnetic induced shift could be cancelled by mixing three kinds of paramagnetic lanthanide ions in appropriate proportion. As a result, the chelating sites would he seen simplely from the half widths and the relative distances between lanthanide ion and the ligand nucleus could be calculated from the relaxation time (T-1) or the half width. Care should be addressed that the analysis method is suitable for the systems in which intramolecular arrangements and intermolecular ligand exchanges are relatively fasten NMR time scale used.