可压缩粘性流体力学变分原理的建立

给出了可压缩粘性流动问题变分泛函的代入构成法，即直接从场的方程出发，应用代入原则，构成相应的变分泛函。应用此法给出了最大功率消耗原理、广义变分原理等，同时认为目前采用的最大功率消耗原理值得商榷。

平行平板间流体的流动与传热的计算

在基于动坐标系的库特剪切流的物理模型上 ,提出了两平板间可压缩粘性流体的温度 ,流速与热流速率间的关联式。

刚 - 流 - 弹耦合系统动力方程及其动力边界条件的建立

应用Ｊｏｕｒｄａｉｎ变分原理，建立了由刚体、两种互不相溶可压缩粘性流体以及弹性体组成系统的动力方程及其动力边界条件，完善和推广了Ｒｕｍｙｊａｎｔｓｅｖ［４，５］、Ｌｉｕ［６］的工作。

冲击波作用下Air/SF_6斜界面不稳定性实验和数值模拟研究

开展了Mach数为1.23和1.41的冲击波作用下的Air/SF6斜界面不稳定性激波管实验,并利用王涛等人发展的可压缩多介质粘性流体和湍流大涡模拟程序MVFT(multi-viscous-fluid and turbulence),对该激波管实验进行了数值模拟,二者相比较一致性较好,包括界面图像、湍流混合区TMZ(turbulent mixing zone)宽度、气泡和尖钉位移,确认了该计算代码对界面不稳定性问题模拟的可靠性和有效性.数值模拟再现了冲击波作用下,Air/SF6斜界面的演化过程及流动中复杂波系结构的发展如冲击波的传播、折射和反射.结果还显示冲击波Mach数较大时,冲击波和界面相互作用时混合区获得的能量也较大,扰动界面发展的也更快.

冲击波三次加载下RM不稳定性数值模拟

采用Smagorinsky和Vreman亚格子尺度应力模型，发展了适用于可压缩多介质粘性流体和湍流大涡模拟的数值编码MVFT（Multi-Viscous—Fluid and Turbulence）。利用MVFT编码对某冲击波二次加载下的流体动力学界面不稳定性及其引起湍流混合进行了数值模拟，并和实验结果进行了比较，二者相吻合。Sma—gorinsky和Vreman亚格子尺度应力模型均是基于涡粘性假设，因此都是绝对耗散的模型，但是Vreman模型的耗散相对较小。从本文计算结果的对比可以看出，Vre—man模型的计算结果比Smagorinsky模型的结果要好。

Numerical Simulation of Interaction Between Laminar Flow and Elastic Sheet

A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.

Multigrid fictitious boundary method for incompressible viscous flows around moving airfoils

In this paper, an efficient multigrid fictitious boundary method （MFBM） coupled with the FEM solver package FEATFLOW was used for the detailed simulation of incompressible viscous flows around one or more moving NACA0012 airfoils. The calculations were carded on a fixed multigrid finite element mesh on which fluid equations were satisfied everywhere, and the airfoils were allowed to move freely through the mesh. The MFBM was employed to treat interactions between the fluid and the airfoils The motion of the airfoils was modeled by Newton-Euler equations. Numerical results of experiments verify that this method provides an efficient way to simulate incompressible viscous flows around moving airfoils.

Recent development of vortex method in incompressible viscous bluff body flows

Vortex methods have been alternative tools of finite element and finite difference methods for several decades. This paper presents a brief review of vortex method development in the last decades and introduces efficient vortex methods developed for high Reynolds number bluff body flows and suitable for running on parallel computer architectures. Included in this study are particle strength exchange methods, core-spreading method, deterministic particle method and hybrid vortex methods. Combined with conservative methods, vortex methods can comprise the most available tools for simulations of three-dimensional complex bluff body flows at high Reynolds numbers.

Modeling of Bubble Motion in Two Dimensional Space

The equations of motion of a bubble, expanding adiabatically through an incompressible viscous fluid, are deduced when the centre of the bubble moves in a vertical plane in the presence of gravitational acceleration, acting vertically downwards. The non-linear equations of motion obtained are solved numerically for different values of the various parameters of the problem. The path traced by the centre of the bubble and velocity of the centre, the change of radius R with time, and the influence of the buoyancy force, which is experienced by the expanding bubble for different values of the gravitational acceleration on these quantities, are investigated. The radius R（t） of the bubble is found to vary periodically with time when the acceleration due to gravity is small. But when the acceleration due to gravity increases, this periodicity in the value of R（t） with t is lost. The influence of viscosity in determining the periodicity of the bubble motion is also investigated.

Nonlinear instability for nonhomogeneous incompressible viscous fluids

We investigate the nonlinear instability of a smooth steady density profile solution to the threedimensional nonhomogeneous incompressible Navier-Stokes equations in the presence of a uniform gravitational field,including a Rayleigh-Taylor steady-state solution with heavier density with increasing height(referred to the Rayleigh-Taylor instability).We first analyze the equations obtained from linearization around the steady density profile solution.Then we construct solutions to the linearized problem that grow in time in the Sobolev space H k,thus leading to a global instability result for the linearized problem.With the help of the constructed unstable solutions and an existence theorem of classical solutions to the original nonlinear equations,we can then demonstrate the instability of the nonlinear problem in some sense.Our analysis shows that the third component of the velocity already induces the instability,which is different from the previous known results.

A semi-implicit three-step method based on SUPG finite element formulation for flow in lid driven cavities with different geometries

A numerical algorithm using a bilinear or linear finite element and semi-implicit three-step method is presented for the analysis of incompressible viscous fluid problems. The streamline upwind/Petrov-Galerkin (SUPG) stabilization scheme is used for the formulation of the Navier-Stokes equations. For the spatial discretization, the convection term is treated explicitly, while the viscous term is treated implicitly, and for the temporal discretization, a three-step method is employed. The present method is applied to simulate the lid driven cavity problems with different geometries at low and high Reynolds numbers. The results compared with other numerical experiments are found to be feasible and satisfactory.

Blowup mechanism for viscous compressible heat-conductive magnetohydrodynamic flows in three dimensions

We investigate initial-boundary-value problem for three-dimensional magnetohydrodynamic （MHD） system of compressible viscous heat-conductive flows and the three-dimensional full compressible Navier-Stokes equations. We establish a blowup criterion only in terms of the derivative of velocity field, similar to the Beale^Kato-Majda type criterion for compressible viscous barotropic flows by Huang et al. （2011）. The results indicate that the nature of the blowup for compressible MHD models of viscous media is similar to the barotropic compressible Navier-Stokes equations and does not depend on further sophistication of the MHD model, in particular, it is independent of the temperature and magnetic field. It also reveals that the deformation tensor of the velocity field plays a more dominant role than the electromagnetic field and the temperature in regularity theory. Especially, the similar results also hold for compressible viscous heat-conductive Navier-Stokes flows, which extend the results established by Fan et al. （2010）, and I-Iuang and Li （2009）. In addition, the viscous coefficients are only restricted by the physical conditions in this paper.

Pressure Oscillating Flow in Corrugated Parallel Channel

The approximate analytical solution of velocity is presented for incompressible and viscous fluid driven by the oscillation of the periodic pressure, between two slit parallel plates with corrugated walls by employing perturbation method. The corrugations of the two walls are described as periodic sinusoidal waves with small amplitude either in phase or half-period out of phase. Based on the analysis, we discuss the influence of the dimensionless parameters on velocity u±and mean velocity parameter φ±numerically, such as Reynolds number Re, nondimensional amplitude A of pressure gradient and wave number k.

On Regularity of Incompressible Fluid with Shear Dependent Viscosity

The authors consider a non-Newtonian fluid governed by equations with p-structure in a cubic domain.A fluid is said to be shear thinning(or pseudo-plastic) if 1 < p < 2,and shear thickening(or dilatant) if p > 2.The case p > 2 is considered in this paper.To improve the regularity results obtained by Crispo,it is shown that the secondorder derivatives of the velocity and the first-order derivative of the pressure belong to suitable spaces,by appealing to anisotropic Sobolev embeddings.

Weak Continuity and Compactness for Nonlinear Partial Differential Equations

This paper presents several examples of fundamental problems involving weak continuity and compactness for nonlinear partial differential equations, in which compensated compactness and related ideas have played a significant role. The compactness and convergence of vanishing viscosity solutions for nonlinear hyperbolic conservation laws are first analyzed, including the inviscid limit from the Navier-Stokes equations to the Euler equations for homentropic flow, the vanishing viscosity method to construct the global spherically symmetric solutions to the multidimensional compressible Euler equations, and the sonic-subsonic limit of solutions of the full Euler equations for multi-dimensional steady compressible fluids. Then the weak continuity and rigidity of the Gauss-Codazzi-Ricci system and corresponding isometric embeddings in differential geometry are revealed. Further references are also provided for some recent developments on the weak continuity and compactness for nonlinear partial differential equations.