He II常流体粘度实验测量技术研究进展

目前还没有现成的理论可以对He II粘度进行合理的说明和计算,在大多数情况下还只能依靠实验技术来满足其工程需要。本文对传统实验技术,目前出现的各种改进型实验技术进行了概述,事实表明不同的实验技术所得的结果之间偶有符合,但不能从根本上解决He II粘度对实验技术的依赖性,目前的实验技术还无法提供He II粘度数值的一致性结论。但是通过对He II粘度与温度之间的变化趋势、数值的数量级等方面的分析,证明了He II粘度在1.6 K^1.8 K之间形成了最小值区域。当温度继续低于1.6 K时,He II粘度会迅速增加,同样在温度高于1.8 K时,随着温度的增加,He II粘度也会迅速增加。

常密度与变密度流体双组分溶质运移实验研究

为准确预测地下水系统中污染物的运移路径及浓度分布,需深入理解不同密度条件下多组分溶质的稀释及运移机理。本文通过室内实验,分别在常密度流体与变密度流体中,系统观测了荧光素钠与氯化钠两种化合物在不同混合比时的运移现象,量化了出口处的浓度分布、稀释指数及两种溶质的浓度标准差。研究发现:在常密度流体中,两种化合物的运移相互独立,其浓度分布取决于化合物各自的水动力扩散与弥散;而在变密度流体中,两种溶质的运移相互影响,其运移轨迹与浓度分布均相同,取决于占比较大化合物的运移结果;变密度对流导致溶质羽和周围流体的接触面积增大,相较于常密度条件,溶质的稀释程度得到大幅提升。

不可压缩流体的数学和物理解析

不可压缩流体的概念和内涵在冶金传输原理基本理论和计算冶金学的教学中具有重要地位。在以往的教材和教学中,不可压缩流体等同于常密度流体。为了探究这一思想的由来,回顾了不可压缩流体的数学定义,分析了不可压缩流体、均质流体和常密度流体之间的区别,给出了工程中不可压缩流体的判断标准。基于Boussinesq近似解释了将不可压缩流体视作常密度流体的理论基础,根据不可压缩流体和定常流动的定义式推导出作定常流动的不可压缩流体的密度是常数的结论。

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.

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.

Numerical computation and analysis of unsteady viscous flow around autonomous underwater vehicle with propellers based on sliding mesh

The flexible transmission shaft and wheel propeller are combined as the kinetic source equipment, which realizes the nmlti-motion modes of the autonomous underwater vehicle （AUV） such as vectored thruster and wheeled movement. In order to study the interactional principle between the hull and the wheel propellers while the AUV navigating in water, the computational fluid dynamics （CFD） method is used to simulate numerically the unsteady viscous flow around AUV with propellers by using the Reynolds-averaged Navier-Stokes （RANS） equations, shear-stress transport （SST） k-w model and pressure with splitting of operators （PISO） algorithm based on sliding mesh. The hydrodynamic parameters of AUV with propellers such as resistance, pressure and velocity are got, which reflect well the real ambient flow field of AUV with propellers. Then, the semi-implicit method for pressure-linked equations （SIMPLE） algorithm is used to compute the steady viscous flow field of AUV hull and propellers, respectively. The computational results agree well with the experimental data, which shows that the numerical method has good accuracy in the prediction of hydrodynamic performance. The interaction between AUV hull and wheel propellers is predicted qualitatively and quantitatively by comparing the hydrodynamic parameters such as resistance, pressure and velocity with those from integral computation and partial computation of the viscous flow around AUV with propellers, which provides an effective reference to the shady on noise and vibration of AUV hull and propellers in real environment. It also provides technical support for the design of new AUVs.

Analytical Study of Magnetohydrodynamic Propulsion Stability

In this paper an analytical solution for the stability of the fully developed flow drive in a magneto-hydro-dynamic pump with pulsating transverse Eletro-magnetic fields is presented. To do this, a theoretical model of the flow is developed and the analytical results are obtained for both the cylindrical and Cartesian configurations that are proper to use in the propulsion of marine vessels. The governing parabolic momentum PDEs are transformed into an ordinary differential equation using approximate velocity distribution. The numerical results are obtained and asymptotic analyses are built to discover the mathematical behavior of the solutions. The maximum velocity in a magneto-hydro-dynamic pump versus time for various values of the Stuart number, electro-magnetic interaction number, Reynolds number, aspect ratio, as well as the magnetic and electrical angular frequency and the shift of the phase angle is presented. Results show that for a high Stuart number there is a frequency limit for stability of the fluid flow in a certain direction of the flow. This stability frequency is dependent on the geometric parameters of a channel.

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.

Anomalous Transport Characteristics of High Temperature Superconductors and Josephson Currents

The transport characteristics of high temperature superconductor current and Josephson current is inves-tigated in the framework of the modified time-dependent Ginzburg-Landau model and the Lawrence-Doniach model.We evaluated the vortex equation and found that the signs of the high temperature superconductor current and theJosephson current can reverse. Some explicit expressions for different cases are derived, which accord with experimentaldata.

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.

A New Unsteady Fluid Network Approach to Simulate the Characteristics of the Air System of a Gas Turbine System

In this paper, a novel unsteady fluid network simulation method to compute the air system of jet engine was coded to predict the characteristics of pressure, temperature and mass flow rate of the flow and the temperature of the solid in the gas turbine engine. The fluid and solid areas are divided into the network comprised of branches and nodes, and the method solves transient mass, energy conservation equations at each node and momentum conservation equation at each branch by a newly deduced numerical method. With this method, to simulate complicated fluid and solid system in short time becomes possible. To verify the code developed, it has been applied to simulate a gas turbine model against the widely used commercial software Flowmaster. And the comparisons show that the two are in good agreement. Then the verified program is applied to the prediction of the characteristics of a designed turbine disk and air-cooling system associated to it, and useful information is obtained.

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.

Simulation of Hydrodynamic Performance of Drag and Double Reverse Propeller Podded Propulsors

The unsteady performance of drag and double reverse propeller podded propulsors in open water was numerically simulated using a computational fluid dynamics （CFD） method. A moving mesh method was used to more realistically simulate propulsor working conditions, and the thrust, torque, and lateral force coefficients of both propulsors were compared and analyzed. Forces acting on different parts of the propulsors along with the flow field distribution of steady and unsteady results at different advance coefficients were compared. Moreover, the change of the lateral force and the difference between the abovementioned two methods were mainly analyzed. It was shown that the thrust and torque results of both methods were similar, with the lateral force results having the highest deviation

A New Differential Quadrature Methodology Based on Bernstein Polynomials for Solving the Equations Governing the Unsteady Flow of a Polytropic Gas

Aims of this paper are to improve ADI differential quadrature method （ADI-DQM） based on Bernstein polynomials and add a new application to the differential quadrature method. By using the new methodology, the numerical solutions of the governing equations of unsteady two-dimensional flow of a polytropic gas are investigated. The numerical results reveal that the new technique is very effective and gives high accuracy, good convergence and reasonable stability.

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.

Developing Spline Based Overset Grid Assembling Approach and Application to Unsteady Flow Around a Moving Body

Overset or Chimera grid approach is one of methods to cope with complex geometries. A spline based overset grid assembling system has been developed. The system is based on structured grid approach and covers comprehensive features for overset assembling, i.e., grid generation, grid modification, and computing Domain Connectivity Information （DCI） for overset interpolation. Ferguson spline curve is used to compute curves through grid points and cell centers for trimming grids and computing DCI robustly and accurately. Flow simulation around Kriso Container Ship （KCS） with jointed grids shows good continuity of flow field between the grids. The overset grid assembling is enhanced to unsteady problem as dynamic overset approach coupled with a solver which also has been developed in National Maritime Research Institute, Japan. Computed results for pitchup spheroid are compared with measured data and show good agreement in unsteady force acting on the spheroid. It is confirmed that the system has capability to simulate flow field around jointed grids and unsteady flow with dynamic overset assembling approach practically.

Numerical Study of Unsteady Behavior of Partial Cavitation on Two Dimensional Hydrofoils

This paper deals with time dependent performance characteristics of cavitating hydrofoils, the flow around which has been simulated using pressure-based finite volume method. A bubble dynamics cavitation model was used to investigate the unsteady behavior of cavitating flow and describe the generation and evaporation of vapor phase. For choosing the turbulence model and mesh size a non cavitating study was conducted. Three turbulence models such as Spalart-Allmaras, Shear Stress Turbulence （SST） κ-ω model, Re-Normalization Group （RNG） κ-ε model with enhanced wall treatment are used to capture the turbulent boundary layer along the hydrofoil surface. The cavitating study presents an unsteady behavior of the partial cavity attached to the foil at different time steps for σ = 0.8 and σ = 0.4. Moreover, this study is focused on cavitation inception, the shape and general behavior of sheet cavitation, lift and drag forces for different cavitation numbers.

PIV measurement on internal instantaneous flows of a centrifugal pump

In this paper,the internal flow field in a centrifugal pump working at the design flow rate operating condition has been measured using the particle image velocimetry(PIV)technique with the fluorescent particles and the index-matched fluid technology.The index-marching fluid with the same refractive index as the transparent material has been prepared and applied in the present test of pump with geometrical complex walls.The comparison between velocity distributions of PIV results both with and without fluorescent particles,as well as with and without index-marching fluid are conducted to find the differences.The DES(Detached Eddy Simulation)has been employed to calculate the three-dimensional unsteady turbulent flow in the pump to examine and to certify the PIV measurement results.The DES results of instantaneous flow velocity fields agree with PIV test results with fluorescent particles and index-marching fluid.It is necessary to perform the PIV measurement of pumps with fluorescent seeds and index-marching fluid in order to get exact results.The experimental results show the distributions of velocity,steamlines,and the principal Reynolds normal stress(PRNS)and the principal Reynolds shear stress(PRSS).

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.

Numerical Simulation of Flow Characteristics in Micro Shock Tubes

Recently micro shock tubes have been widely used in many engineering and industrial fields, but the characteristics of unsteady flow are not well known to date in micro shock tubes. Compared to conventional shock tubes with macro scales, flows related to shock waves in micro shock tubes are highly complicated. Stronger viscous and dissipative interactions make shock wave dynamic behaviors significantly different from theoretical predictions. In the present study, a CFD work was applied to the unsteady compressible Navier-Stokes equations which were solved using a fully implicit finite volume scheme. The diaphragm pressure ratio and shock tube diameter were varied to investigate their effects on micro shock tube flows. Different wall boundary conditions were also performed to observe shock wave and contact surface propagation with no slip and slip walls. Detailed flow characteristics at the foot of shock wave and contact surface propagation were known from the present numerical simulations.