Non-relativistic conformal symmetries in fluid mechanics

The symmetries of a free incompressible fluid span the Galilei group, augmented with independent dilations of space and time. When the fluid is compressible, the symmetry is enlarged to the expanded Schrodinger group, which also involves, in addition, Schrodinger expansions. While incompressible fluid dynamics can be derived as an appropriate non-relativistic limit of a conformally invariant relativistic theory, the recently discussed conformal Galilei group, obtained by contraction from the relativistic conformal group, is not a symmetry. This is explained by the subtleties of the non-relativistic limit.

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

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

Generation of internal and surface waves by seafloor movement in a two layer fluid system

Internal and surface waves generated by the deformations of the solid bed in a two layer fluid system of infinite lateral extent and uniform depth are investigated. An integral solution is developed for an arbitrary bed displacement on the basis of a linear approximation of the complete description of wave motion using a transform method (Laplace in time and Fourier in space) analogous to that used to study the generation of tsunamis by many researchers. The theoretical solutions are presented for three interesting specific deformations of the seafloor; the spatial variation of each seafloor displacement consists of a block section of the seafloor moving vertically either up or down while the time-displacement history of the block section is varied. The generation process and the profiles of the internal and surface waves for the case of the exponential bed movement are numerically illustrated, and the effects of the deformation parameters, densities and depths of the two layers on the solutions are discussed. As expected, the solutions derived from the present work include as special cases that obtained by Kervella et al. [Theor Comput Fluid Dyn 21:245-269, 2007] for tsunamis cased by an instantaneous seabed deformation and those presented by Hammack [J Fluid Mech 60:769-799, 1973] for the exponential and the half-sine bed displacements when the density of the upper fluid is taken as zero.

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.

西峰油田特低渗透油藏渗流机理及开发对策研究

Micro-pore-throat, micro-fracture and low permeability are the most obvious characters of Xifeng ultra-low permeability reservoir, and threshold pressure gradient and medium deformation during the period of oilfield developing results non-linear seepage feature of the formation liquid flowing in the porous medium underground. It is impossible to solve some problems in the ultra-low permeability reservoir development by current Darcy filtration theory and development techniques. In the view of the characters of ultra-low permeability and powerful-diagenesis and fracture up-growth, the paper quantitatively characterizes of through-going scope for reservoir parameters together with some materials such as similarity field outcrop, rock core, drilling, well logging and production dynamic, which provides geological base for further development adjustment. Based on the displacement experiment of different kinds of seepage fluid and oil-water two phases, this paper proves the relationship between threshold pressure gradient and formation permeability in experiment and theory, which is power function and its index is about -1. The variation rule and the mechanism of oil-water two phases threshold pressure gradient are studied. At the same time, based on the experiment of medium deformation, the variation rule of formation physical property parameters and the deformation mechanism are researched, and the influential factors on the medium deformation are analyzed systematically. With elastic unsteady filtration theory, nonlinear mathematical models of the steady and unsteady flow of single phase as well as horizontal well flow and oil-water two phases flow are deduced with the influence of nonlinear factors including threshold pressure gradient and media deformation. The influences of nonlinear factors upon well deliverability and reservoir pressure distribution as well as the saturation variation pattern of oil-water front are analyzed. By means of the researches such as reasonable well pattern, reasonable well array ration, artificial fracture length optimization advisable water flood timing and feasibility of advanced water flooding, it is necessary to find out effective techniques in order to improve development result of this kind of reservoir. This research result develops and improves on low-velocity nonlinear seepage theory, and offers ways to study similar kind of reservoir; it is meaningful to the development of the ultra-low permeability oil and gas reservoir.

流体动力学的研究:面向国家需求

流体动力学这门科学描述液体和气体的运动及其和固体的相互作用.它几乎涉及日常生活的方方面面,并在大部分科学和工程领域中处于核心地位,所以,这是一门宽广的多学科的学科领域.流体力学影响着国防、国土安全、交通、制造业、医学、生物学、能源和环境科学.预测人体内血液的流动,微流动器件装置的行为,飞机、汽车、船舶的流体动力学性能,电子元件的冷却,或者灾害性天气和气候的预报等,都需要对流体动力学有细致的了解和深入的研究.

Abel Inversion of Axially-Symmetric Shock Wave Flows

Finite-fringe interferograms produced for axisymmetric shock wave flows are analyzed by Fourier transform fringe analysis and an Abel inversion method to produce density field data for the validation of numerical models. For the Abel inversion process, we use basis functions to model phase data from axially-symmetric shock wave structure. Steady and unsteady flow problems are studied, and compared with numerical simulations. Good agreement between theoretical and experimental results is obtained when one set of basis functions is used during the inversion process, but the shock front is smeared when another is used. This is because each function in the second set of basis functions is infinitely differentiable, making them poorly-suited to the modelling of a step function as is required in the representation of a shock wave.

Pulverized-coal combustion in Wall-Protecting-Jets combustor (WPJC)

A new type of pulverized-coal combustor, called "Wall-Protecting-Jets Combustor" (hereafter, WPJC has been proposed, designed and studied with both CFD (Computational Fluid Dynamics) and experimental methods. The WPJC is based on a novel concept in which all inlet jets are along the combustor wall. Pilot combustion experiments were conducted to investigate the combustion performance of WPJC. Two-phase flows and pulverized-coal combustion were simulated to study the mechanism of),WPJC using the commercial software FLUENT. The results show that the WPJC has many remarkable advantages: wall-protection by the cold jets without the use of refractory materials; low-temperature and three-stage combustion with low NOx emission; negligible ash/slag-deposition; multiple functions with convenient switching between them; effective adjustment of the combustion intensity and the ignition position.

Numerical exploration of new features on three-dimensional transition of a cylinder wake

The three-dimensional transition of the wake flow behind a circular cylinder is studied in detail by direct numerical simulations using 3D incompressible N-S equations for Reynolds number ranging from 200 to 300. New features and vortex dynamics of the 3D transition of the wake are found and investigated. At Re = 200, the flow pattern is characterized by mode A instability. However, the spanwise characteristic length of the cylinder determines the transition features. Particularly for the specific spanwise characteristic length linear stable mode may dominate the wake in place of mode A and determine the spanwise phase difference of the primary vortices shedding. At Re = 250 and 300 it is found that the streamwise vortices evolve into a new type of mode - "dual vortex pair mode" downstream. The streamwise vortex structures switch among mode A, mode B and dual vortex pair mode from near wake to downstream wake. At Re = 250, an independent low frequency f(m) in addition to the vortex shedding frequency f(s) is identified. Frequency coupling between f(m) and f(s) occurs. These result in the irregularity of the temporal signals and become a key feature in the transition of the wake. Based on the formation analysis of the streamwise vorticity in the vicinity of cylinder, it is suggested that mode A is caused by the emergence of the spanwise velocity due to three dimensionality of the incoming flow past the cylinder. Energy distribution on various wave numbers and the frequency variation in the wake are also described.

Waverider configurations derived from general conical flowfields

A method based on the computational fluid dynamics (CFD) is presented for a flexible waverider's design. The generating bodies of this method could be any cones. In addition, either the leading edge or the profile of the scramjet's inlet is used as the waverider's definition curve, parameterized by the quadric function, the sigmoid function or the B-spline function. Furthermore, several numerical examples are carried out to validate the method and the relevant codes. The CFD results of the configurations show that all the designs are successful. Moreover, primary suggestions are proposed for practical design by comparing the geometrical and aerodynamic performances of the cone-derived waveriders at Mach 6.

Numerical investigation of dynamic stall of an oscillating aerofoil

The flow structure around an NACA 0012 aerofoil oscillating in pitch around the quarter-chord is numerically investigated by solving the two-dimensional compressible N-S equations using a special matrix-splitting scheme. This scheme is of second-order accuracy in time and space and is computationally more efficient than the conventional flux-splitting scheme. A 'rigid' C-grid with 149 x 51 points is used for the computation of unsteady flow. The freestream Mach number varies from 0.2 to 0.6 and the Reynolds number from 5000 to 20,000. The reduced frequency equals 0.25-0.5. The basic flow structure of dynamic stall is described and the Reynolds number effect on dynamic stall is briefly discussed. The influence of the compressibility on dynamic stall is analysed in detail. Numerical results show that there is a significant influence of the compressibility on the formation and convection of the dynamic stall vortex. There is a certain influence of the Reynolds number on the flow structure. The average convection velocity of the dynamic stall vortex is approximately 0.348 times the freestream velocity.

Free surface vibration in oscillatory convection of half floating zone

Projecting an orthographical grating mask (20pl/mm) on the surface of a small liquid bridge and receiving the reflected distortion image, one can calculate out reversely the shape of free surface of a liquid bridge. In this way we measured the surface shape of a small floating zone and the two-dimensional deformation of its vibration. The mechanism of thermocapillary oscillatory convection and the three-dimensional variation of the free surface are revealed experimentally. The principle for space experiment has been studied in our laboratory.

High-speed streams from coronal holes and the accelerating mechanism of the solar-wind

In this paper, the general Mach number equation is derived, and the influence of typical energy forms in the solar wind is analysed in detail. It shows that the accelerating process of the solar wind is influenced critically by the form of heating in the corona, and that the transonic mechanism is mainly the result of the adjustment of the variation of the crosssection of flowing tubes and the heat source term.The accelerating mechanism for both the high-speed stream from the coronal hole and the normal solar wind is similar. But, the temperature is low in the lower level of the coronal hole and more heat energy supply in the outside is required, hence the high speed of the solar wind; while the case with the ordinary coronal region is just the opposite, and the velocity of the solar wind is therefore lower. The accelerating process for various typical parameters is calculated, and it is found that the high-speed stream may reach 800 km/sec.

封闭和半封闭通风空间的环境流体问题

半封闭与全封闭空间的环境流体力学包括系统内的空气流动、温度和湿度分布以及污染物的凝积等问题，其中涉及控制空间的空气质量、通风设计以及预防水汽凝结等多个方面的研究．本文着重介绍国内外近年来相关的工作，研究如何从理论分析、实验研究以及数值模拟等方面入手，解决通风问题，以达到获得高通风效率、低能耗和提高空气舒适度的目的．研究的对象包括房间，建筑物，城市等半封闭空间，以及汽车、地铁、飞机、潜艇、载人航天飞行器等不同程度的全封闭空间．现有的理论研究已经从解释单点源纯浮力对流的基础模型发展到能够描述具有体积流量和动量流量的扩散流的湍流喷泉模型．然而，理论研究上的进步还不足以使它能够处理实际问题中多样的通风情况和复杂的空间结构，工程应用中还需更多依靠实验研究和数值模拟的结论．对载人航天飞行器中的环境流体问题的研究虽然刚刚起步，但其研究手段和经验可以从相对成熟的半封闭与全封闭空间的相关研究中获得借鉴．

三维直叶栅非定常流动的并行计算研究

流动分离直接关系到压气机运行的安全性与效率，对分离流动的研究是叶轮机械真实流动研究中的一个重大课题。本文针对三维压气机单转子叶片中截面所构成的三维直叶栅跨音速分离流开发了通用数值计算程序。该程序基于B-L湍流模型及高精度差分方法。多种工况的数值计算显示本程序结果与实验值吻合比较理想，验证了程序的正确性。10°攻角下分离区脉动压力的频谱与实验结果的数量级吻合，说明本程序能够较好地模拟大攻角分离流这种非定常复杂流动，为了提高计算规模及计算速度，作者对程序进行了并行化并针对微机机群系统进行了并行优化，实际计算表明本程序具有较高的并行效率。