最优同伦渐近法在带滑移的6常数Oldroyd流体模型中的应用

应用最优同伦渐近法(OHAM)研究了6常数Oldroydl流体的平面Couette流、平面Poiseuille流和平面Couette Poiseuille流这三个非线性问题.介绍了6常数Oldroyd流体的平面Couette流、平面Poiseuille流和平面Couette Poiseuille流的控制方程和边界条件,应用最优同伦渐近法对上述问题进行求解,并给出模拟仿真,根据仿真结果讨论了滑移参数对速度分布的影响.

钕系顺丁橡胶-甲苯溶液的表观黏度及流动参数

以中国石油独山子石化公司研究院中试装置生产的钕系顺丁橡胶(NdBR)为胶样,以甲苯为溶剂,制备了NdBR-甲苯溶液,研究了胶液的质量分数(C)、温度(T)和剪切速率(γ)3个因素对胶液表观黏度(η)的影响,并采用外推法考察了胶液的流体常数(a)和零剪切黏度(η0)。结果表明,3个因素对η的影响顺序为:C≥γ≈T;在C分别为8.00%～10.00%的低浓度区和10.00%～12.00%高浓度区,其对应的γ分别为4.40～25.00 s-1和0.88～15.00 s-1时,η的回归方程分别为:lnη=0.815 2+3.533 5 ln C+13.405 4/T-0.210 8 lnγ和lnη=-6.593 7+6.751 9 ln C+13.466 9/T-0.212 0 lnγ;2个回归方程的适用T均为30～70℃。a=0.000 635;η0的经验方程为lnη0=-5.674 2+6.464 4 ln C+12.289 2/T,影响胶液η0因素的顺序为:C≥T。

国内外介电测井技术综述

以SHAD-2000高频阵列介电测井仪为例,对国内外介电测井技术进行综述,此项技术是利用烃和水的相对介电常数值较大来区分油水层。当常规电测井在某些地层中无法进行区分油水层时,介电仪器与对应解释成果可以更好区分油层和水层。对XX地区,白加XX井,王xx井,坪XX井介电解释成果表进行了分析研究,具有一定实用推广价值。

Modeling of fine coal flotation separation based on particle characteristics and hydrodynamic conditions

Flotation is a complex multifaceted process that is widely used for the separation of finely ground minerals. The theory of froth flotation is complex and is not completely understood. This fact has been brought many monitoring challenges in a coal processing plant. To solve those challenges, it is important to understand the effect of different parameters on the fine particle separation, and control flotation performance for a particular system. This study is going to indicate the effect of various parameters （particle Characteristics and hydrodynamic conditions） on coal flotation responses （flotation rate constant and recovery） by different modeling techniques. A comprehensive coal flotation database was prepared for the statistical and soft computing methods. Statistical factors were used for variable selections. Results were in a good agreement with recent theoretical flotation investigations. Computational models accurately can estimate flotation rate constant and coal recovery （correlation coefficient 0.85, and 0.99, respectively）. According to the results, it can be concluded that the soft computing models can overcome the complexity of process and be used as an expert system to control, and optimize parameters of coal flotation process.

Kerosene-alumina nanofluid flow and heat transfer for cooling application

Kerosene-alumina nanofluid flow and heat transfer in the presence of magnetic field are studied. The basic partial differential equations are reduced to ordinary differential equations which are solved semi analytically using differential transformation method. Velocity and temperature profiles as well as the skin friction coefficient and the Nusselt number are determined analytically. The influence of pertinent parameters such as magnetic parameter, nanofluid volume fraction, viscosity parameter and Eckert number on the flow and heat transfer characteristics is discussed. Results indicate that skin friction coefficient decreases with increase of magnetic parameter, nanofluid volume fraction and viscosity parameter. Nusselt number increases with increase of magnetic parameter and nanofluid volume fraction while it decreases with increase of Eckert number and viscosity parameter.

Calculation of Metzner Constant for Double Helical Ribbon Impeller by Computational Fluid Dynamic Method

Using the multiple reference frames （MRF） impeller method, the three-dimensional non-Newtonian flow field generated by a double helical ribbon （DHR） impeller has been simulated. The velocity field calculated by the numerical simulation was similar to the previous studies and the power constant agreed well with the experimental data. Three computational fluid dynamic （CFD） methods, labeled Ⅰ, Ⅱ and Ⅲ, were used to compute the Metzuer constant k5. The results showed that the calculated value from the slop method （method Ⅰ） was consistent with the experimental data. Method Ⅱ, which took the maximal circumference-average shear rate around the impeller as the effective shear rate to compute ks, also showed good agreement with the experiment. However, both methods suffer from the complexity of calculation procedures. A new method （method Ⅲ） was devised in this paper to use the area-weighted average viscosity around the impeller as the effective viscosity for calculating k5. Method Ⅲ showed both good accuracy and ease of use.

A CFD Based Investigation of the Unsteady Hydrodynamic Coefficients of 3-D Fins in Viscous Flow

The motion of the fins and control surfaces of underwater vehicles in a fluid is an interesting and challenging research subject.Typically the effect of fin oscillations on the fluid flow around such a body is highly unsteady, generating vortices and requiring detailed analysis of fluid-structure interactions.An understanding of the complexities of such flows is of interest to engineers developing vehicles capable of high dynamic performance in their propulsion and maneuvering.In the present study, a CFD based RANS simulation of a 3-D fin body moving in a viscous fluid was developed.It investigated hydrodynamic performance by evaluating the hydrodynamic coefficients (lift, drag and moment) at two different oscillating frequencies.A parametric analysis of the factors that affect the hydrodynamic performance of the fin body was done, along with a comparison of results from experiments.The results of the simulation were found in close agreement with experimental results and this validated the simulation as an effective tool for evaluation of the unsteady hydrodynamic coefficients of 3-D fins.This work can be further be used for analysis of the stability and maneuverability of fin actuated underwater vehicles.

Dependence of Reaction Rate Constants on Density in Supercritical Fluids

A new method,which correlates rate constants of chemical reactions and density or pressure in supercritical fluids,was developed.Based on the transition state theory and thermodynamic principles, the rate constant can be reasonably correlated with the density of the supercritical fluid,and a correlation equation was obtained. Coupled with the equation of state (EOS) of a supercritical solvent,the effect of pressure on reaction rate constant could be represented.Two typical systems were used to test this method.The result indicates that this method is suitable for dilute supercritical fluid solutions.

Response Difference and Related Discussion of Digital Fluid Records of the Beijing Wuliying Well before Strong Earthquakes

In this paper, the precursor response characteristics of digital fluid caused by the Wenchuan M8.0 and Yushu M7.1 earthquakes are studied, and the response difference of the observations of Wuliying well to the two strong earthquakes is compared. The result shows that the abnormal fluid response has a certain relationship with earthquake size and epicenter distance. The greater the earthquake, and the closer it is to the epicenter from the observatory, the more sensitive the response will be to fluid anomalies. Abnormal Helium release was first observed before both strong earthquakes in the fluid precursor observation. The release intensity is related to earthquake magnitude; the larger the magnitude, the stronger the abnormal changes. The large change in He release in a short period after the Wenchuan earthquake may be related to the continuous activity of strong aftershocks and the structural adjustment after the earthquake. Helium release increased significantly after both earthquakes, as contributed by the abnormal deep-sited Helium release. However, this process may be connected to the annual change from July to September or the hot-reservoir type gas release in the Wuliying well. In the earthquake preparation process, a large number of deep-derived Helium is released into the geothermal system and the entire composition is changed. Temporally, this gas release appears later than fault-type gas release, and the disappearance time of this anomaly is also much later. The response difference of the fluid precursor to the two strong earthquakes may be related to differences of deep structural setting and dynamic mechanism. It also shows there is a stronger correlation between Wenchuan M8.0 earthquake in the North-South Seismic Belt and North China region than the Yushu M7.1 earthquake in the internal Qinghai-Tibetan plateau. Helium gas can be a sensitive indicator for monitoring abnormal deep-gas activity of the region where the observation station is located. Hence, observation and research should be strengthened in the future.

Stability and self-adaption character of turbulence coherent structure in narrow-deep river bend

In a meandering fiver, a certain scale of turbulent vortex dominates the development of fiver morphology, making the river bend with s particular curvature. This kind of vortex is denoted as ＂bend-forming vortex＂. The coordinated relationship of bend-forming vortex and meandering fiver channel is then known as ＂self-adaption feature＂ of rivers. With these two concepts, this paper investigated the stability and self-adaption character of coherent vortex in the U-shape river bend with a constant curvature. On the basis of fluid mechanics theory and in consideration of turbulent coherent vortex as disturbance, the growth rate and the wave number response range of coherent vortex in meandering rivers with different curvatures were calculated in this paper. Moreover, the responses of different scales of coherent turbulence structure to river bend parameters were analyzed to explain the mechanism of fiver bend maintenance. These methods could provide a theoretical basis for further investigation on fiver meandering.

Numerical simulation of unsteady cavitating flows around a transient pitching hydrofoil

The objective of this paper is to improve the understanding of the influence of multiphase flow on the turbulent closure model, the interplay between vorticity fields and cavity dynamics around a pitching hydrofoil. The effects of pitching rate on the sub- cavitating and cavitating response of the pitching hydrofoil are also investigated. In particular, we focus on the interactions between cavity inception, growth, and shedding and the vortex flow structures, and their impacts on the hydrofoil performance. The calculations are 2-D and performed by solving the incompressible, multiphase Unsteady Reynolds Averaged Navier Stokes （URANS） equations via the commercial CFD code CFX. The k-co SST （Shear Stress Transport） turbulence model is used along with the transport equation-based cavitation models. The density correction function is considered to reduce the eddy viscosity according to the computed local fluid mixture density. The calculation results are validated with experiments conducted by Ducoin et al. （see Computational and experimental investigation of flow over a transient pitching hydrofoil, Eur J Mech/B Fluids, 2009, 28：728-743 and An experimental analysis of fluid structure interaction of a flexible hydrofoil in vari- ous flow regimes including cavitating flow, Eur J Mech B/fluids, 2012, 36： 63-74）. Results are shown for a NACA66 hydro- foil subject to slow （quasi static, t2=6~/s, ＆＊ =0.18） and fast （dynamic, ＆=63~/s, dr＂ =1.89） pitching motions from a =0~ to a =15~. Both subcavitaing （or =8.0） and cavitating （cr=3.0） flows are considered. For subcavitating flow （or=8.0）, low frequency fluctuations have been observed when the leading edge vortex shedding occurs during stall, and delay of stall is ob- served with increasing pitching velocity. For cavitating flow （tr=3.0）, small leading edge cavities are observed with the slow pitching case, which significantly modified the vortex dynamics at high angles of attack, leading to high frequency fluctuations of the hydrodynamic coefficients and different stall behaviors compared to the subcavitating flow at the same pitching rate. On the other hand, for the fast pitching case at or=3.0, large-scale sheet/cloud cavitation is observed, the cavity behavior is un- steady and has a strong impact on the hydrodynamic response, which leads to high amplitude fluctuations of the hydrodynamic coefficients, as well as significant changes in the stall and post-stall behavior. The numerical results also show that the local density modification helps to reduce turbulent eddy viscosity in the cavitating region, which significantly modifies the cavity lengths and shedding frequencies, particularly for the fast pitching case. In general, compared with the experimental visualiza- tions, the numerical results with local density correction have been found to agree well with experimental measurements and observations for both slow and fast transient pitching cases.

Convergence of some finite element iterative methods related to different Reynolds numbers for the 2D/3D stationary incompressible magnetohydrodynamics

Based on the finite element method(FEM), some iterative methods related to different Reynolds numbers are designed and analyzed for solving the 2D/3D stationary incompressible magnetohydrodynamics(MHD) numerically. Two-level finite element iterative methods, consisting of the classical m-iteration methods on a coarse grid and corrections on a fine grid, are designed to solve the system at low Reynolds numbers under the strong uniqueness condition. One-level Oseen-type iterative method is investigated on a fine mesh at high Reynolds numbers under the weak uniqueness condition. Furthermore, the uniform stability and convergence of these methods with respect to equation parameters R_e, R_m, S_c, mesh sizes h, H and iterative step m are provided. Finally, the efficiency of the proposed methods is confirmed by numerical investigations.