Viscosity-induced mode splitting and potential energy criterion for mode stability

This paper points out that viscosity can induce mode splitting in a uniform infinite cylinder of an incompressible fluid with self-gravitation, and that the potential energy criterion cannot be appropriate to all normal modes obtained, i.e., there will be stable modes with negative potential energy (<0). Therefore the condition >0 is not necessary, although sufficient, for the stability of a mode in an incompressible static fluid or magnetohydrodynamics (MHD) system, which is a correction of both Hare's [Philos. Mag. 8, 1305 (1959)] and Chandrasekhar's [Hydrodynamic and Hydromagnetic Stability (Oxford U.P., Oxford, 1961), p. 604] stability criterion for a mode. These results can also be extended to compressible systems with a polytropic exponent.

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

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.

LOCAL BENDING OF THIN FILM ON VISCOUS LAYER

Effects of deposition layer position and number/density on local bending of a thin film are systematically investigated. Because the deposition layer interacts with the thin film at the interface and there is an offset between the thin film neutral surface and the interface, the deposition layer generates not only axial stress but also bending moment. The bending moment induces an instant out-of-plane deflection of the thin film, which may or may not cause the so-called local bending. The deposition layer is modeled as a local stressor, whose location and density are demonstrated to be vital to the occurrence of local bending. The thin film rests on a viscous layer, which is governed by the Navier-Stokes equation and behaves like an elastic foundation to exert transverse forces on the thin film. The unknown feature of the axial constraint force makes the governing equation highly nonlinear even for the small deflection case. The constraint force and film transverse deflection are solved iteratively through the governing equation and the displacement constraint equation of immovable edges. This research shows that in some special cases, the deposition density increase does not necessarily reduce the local bending. By comparing the thin film deflections of different deposition numbers and positions, we also present the guideline of strengthening or suppressing the local bending.

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.

一种基于粒子的牛顿流体与粘弹性流体统一模拟方法

在越来越受到人们关注的基于物理流体动画领域,目前分别模拟牛顿流体或粘弹性流体的方法很多,但很少有统一模拟两者的方法.文中基于光滑粒子流体动力学方法,通过对传统纳维-斯托克斯方程添加弹性应力项,提出了一种新的统一模拟牛顿流体和粘弹性流体的方法.通过实验说明该方法不仅有效,易于实现,而且具有良好的可控性,仅仅通过调节参数就可以模拟不同粘弹性、不同类型的流体现象.

Viscous-force-dominated tensile deformation behavior of oriented polyethylene

High-density polyethylene with shish-kebab structure, prepared by a melt extrusion drawing, was employed to investigate the effect of the well-defined lamellar orientation on the deformation characteristics under uniaxial tensile deformation along the drawing direction. This was done by investigating the true stress-true strain dependencies at different strain rates, recovery properties, and stress relaxation measurements. Measurements were complemented by recording in-situ wide-angle X-ray scattering patterns during the deformation process. The oriented samples showed not only a higher modulus, but different from analogous isotropic samples, a homogeneous deformation without necking. The true strain associated with the onset of fibrillation was determined. Because of the preorientation, it is shifted to 0.3, which is below the value 0.6 of the isotropic counterpart. The main finding is a strong enhancement of the Viscous force, as was revealed by stress relaxation experiments; the viscous force takes up 70% of the total stress. The presence of shish-kebabs, i.e., interconnected lamellae in a stack, seems to be responsible for the high viscous force in the oriented samples. The absence of necking has to be ascribed to the high viscous force.

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.

吐合盆地鲁克沁构造带断块油藏形成机制及油藏预测

Lukeqin arc belt is a compound structure generated by multi-movements and composed of 6 sub-structural zones, which are connected by Huoyanshan Mountain. General characteristics of the arc belt are multi-patterns of structure, multi-phases for petroleum, multi-types of trap and multi-layers for reservoirs. As a part of the eastern Lukeqin arc belt located on the south of Taibei depression, Lukeqin structural zone behaves as a complex faulted-fold zone, in which the formation and distribution of hydrocarbons are controlled by structures. As the dominant source of dynamics for the second migration of hydrocarbon, structure stress field is closely related with the potentials of hydrodynamics. Results derived from the simulations of stress field by finite element method indicate that the northwest tending faults prefer seal to the northeast tending ones. The reason is that the northwest tending faults were squeezed more strongly than the northeast tending ones. Therefor, the northeast tending faults become always the paths for oil to migrate southeastward. Lukeqin structural zone is the main site for oil to concentration because it is surrounded by high stress. Situated on the front of the foreland basin of Turpan-Hami, Lukeqing arc belt is a dam to hold back the southward migrating oil from Shengbei depression. The axis line of Shenquan-Shengnan-Yanmuxi, Lukeqin and Yubei controls the migrating paths and concentrating process of oil and gas. Results derived from stress simulation and structure analyses indicate consistently that both Yubei and Lukeqin structural zones are the favorite areas for oil to migrate. The generally southward paths for oil to migrate out of Taibei depression can be two ways. One of them is from Taibei depression to Yubei structural zone and the other is from Taibei depression to Lukeqin structural zone. By the both ways, oil migrated upward along the faults and southeastward along the structural axis to concentrate in either Permian or Triassic system. The newly ascertained path for oil migration, which is accurately southeastward instead of coarsely southward, indicates the directions for further explorations on the compound Lukeqin block zone. Five kinds of seal models of fault are all found in Lukeqin block zone by studying the seal features of faults occurred in the zone. Having studied the fault seal and their controlling factors by fuzzy set method, the paper deems that the northwest tended faults are better than the northeast tended ones for oil to concentrate. The most important factors to decide the seal extent of faults in this zone are the characteristics of main stress and fluids instead of capillary pressure differences between the two sides of fault and smear mud factors. There exist seal differences not only between the faults of different time but also between the sections within a fault due to the variation of depths, strata and positions. The general distribution rules of reservoirs were dominated by the seal characteristics of a fault during the time reservoirs formed. While the current features of fault seal decide the conservation of reservoirs and heights of oil accumulations. Seal or not of a fault is not absolute because the essential for fault to seal is the distribution of permeability of fault zone. Therefor, the multi cyclical activities of faults create the space-time variation of seal features of the fault. Totally, the seal extent of the faults within the area is not as perfect as to accumulate ordinary crude. Crude oil can only be sealed when it becomes viscous. Process for crude oil to become viscous and viscous happened strongly because of the fault-fold movements. Shallowly burying and even revealing of the objective layers of the reservoirs made the crude oil to be thickened by water washing biologically degradation and oxidation degradation. The northwestward deepening during or after the reservoir formation of the structural zone provided the power for oil to migrate one or more times. The main reason for oil accumulation is the formation of Lukeqin block zone during Xishanyao stage, middle Jurassic Period, Early Yanshanian Movement. While the main reason for reservoir conservation is the placidity of Triassic blocks after the formation of reservoirs. Contrasting to former opinions, it is concluded that the reservoirs in Lukeqin zone, including viscous reservoirs, were formed by one time but not more times. So the author proposes the opinion that the reservoirs of viscous oil were formed by viscous oil migration under the conditions of aptitude sets of fault seals controlled by fluid and other factors. To grope the distribution rules outside Taibei depression and discuss the formation mechanism of Anjurassic reservoirs, it is necessary to study the dominate factors for the formation of reservoirs in Lukeqin structural zone such as structural stress, fault seals and thickening mechanism of crude oil. Also, the necessary studies are the key to break through the Taibei depression and Anjurassic systems. Therefor, they are significant for the future exploration and reserve increasing of hydrocarbon within the Turpan-Hami basin. The paper studied the distribution rules of block reservoirs and forecasted the favorable zones for further exploration in Turpan-Hami basin. Conclusions can be useful for not only the exploration in the area but also the theory consult in the adjacent areas.

Mantle-derived gaseous components in ore-forming fluids of the Xiangshan uranium

The Xiangshan U deposit, the largest hydrothermal U deposit in China, is hosted in late Jurassic felsic volcanic rocks although the U mineralization post dates the volcanics by at least 20 Ma. The mineralization coincides with intrusion of local mantle-derived mafic dykes formed during Cretaceous crustal extension in South China. Ore-forming fluids are rich in CO2, and U in the fluid is thought to have been dissolved in the form of UO2 (CO3)22− and UO2 (CO3) 34− complexes. This paper provides He and Ar isotope data of fluid inclusions in pyrites and C isotope data of calcites associated with U mineralization (pitchblende) in the Xiangshan U deposit. He isotopic compositions range between 0.1 and 2.0Ra (where Ra is the 3He/4He ratio of air=1.39×10−6) and correlates with 40Ar/36Ar; although there is potential for significant 3He production via 6Li(n,α)3H(β)3He reactions in a U deposit (due to abundant neutrons), nucleogenic production cannot account for either the 3He concentration in these fluids, nor the correlations between He and Ar isotopic compositions. It is more likely that the high 3He/4He ratios represent trapped mantle-derived gases. A mantle origin for the volatiles of Xiangshan is consistent with the δ13C values of calcites, which vary from −3.5‰ to −7.7‰, overlapping the range of mantle CO2. The He, Ar and CO2 characteristics of the ore-forming fluids responsible for the deposit are consistent with mixing between 3He- and CO2-rich mantle-derived fluids and CO2-poor meteoric fluids. These fluids were likely produced during Cretaceous extension and dyke intrusion which permitted mantle-derived CO2 to migrate upward and remobilize U from the acid volcanic source rocks, resulting in the formation of the U deposit. Subsequent decay of U within the fluid inclusions has reduced the 3He/4He ratio, and variations in U/3He result in the range in 3He/4He observed with U/3He ratios in the range 5–17×103 likely corresponding to U concentrations in the fluids b0.2 ppm.