数值模拟二元喷管的三维可压粘性流场

发展了一种用扫描有限元素法求解二元喷管内部湍流流场的数值模拟方法，并提出了ＮＳ方程中二阶导数项新的离散方法，推导了更完善、更简洁的扫描有限元方程的表达式。用本文发展的方法对二维平面收扩喷管、两个大宽高比的二元喷管（宽高比分别为７和４）及其圆转方转接段的三维粘性流场等算例的计算结果表明，该方法具有较宽广的适用范围，有较高的应用和研究价值。

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.

A Numerical Modeling of the Vortex-induced Vibration of Cascade in Turbomachinery using Immersed Boundary Method

Based on the immersed boundary method,a fast simulation for solving unsteady,incompressible,viscous flow associated with the oscillating cascade is established on a quasi-three-dimensional coordinate system.The numerical method is applied to the simulation of the flow passing an oscillating circular cylinder which is forced to move in X direction under prescribed motions in water at rest at low Keulegan-Carpenter numbers.Then vor-tex-induced vibration of a cylinder with two degrees of freedom which oscillates in in-line direction and transverse direction is simulated using this method.The results are in good agreement with the previous research.Then the method is extended to the oscillating cascade simulation of making various comparisons.It is found that the IBPA(inter blade phase angle) will change as the time goes on,because of the non-uniformity of the flow in the circumferential direction,until the oscillating cascade goes to a stable situation.The reduced velocity and the number of blades are chosen to investigate the effects of them on IBPA.The results indicate that both the reduced velocity and the number of blades are the main factors which influence IBPA.It is worth noting that the coupling process is not necessary to generate any body-fitting grids,which makes it much faster in computational process for such a complicated fluid-structure interaction problem.

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.

Suitability of artificial viscosity discontinuous Galerkin method for compressible turbulence

A discontinuous Galerkin method based on an artificial viscosity model is investigated in the context of the simulation of compressible turbulence. The effects of artificial viscosity on shock capturing ability, broadband accuracy and under-resolved instability are examined combined with various orders and mesh resolutions. For shock-dominated flows, the superior accuracy of high order methods in terms of discontinuity resolution are well retained compared with lower ones. For under-resolved simulations, the artificial viscosity model is able to enhance stability of the eighth order discontinuous Galerkin method despite of detrimental influence for accuracy. For multi-scale flows, the artificial viscosity model demonstrates biased numerical dissipation towards higher wavenumbers. Capability in terms of boundary layer flows and hybrid meshes is also demonstrated.It is concluded that the fourth order artificial viscosity discontinuous Galerkin method is comparable to typical high order finite difference methods in the literature in terms of accuracy for identical number of degrees of freedom, while the eighth order is significantly better unless the under-resolved instability issue is raised. Furthermore, the artificial viscosity discontinuous Galerkin method is shown to provide appropriate numerical dissipation as compensation for turbulent kinetic energy decaying on moderately coarse meshes, indicating good potentiality for implicit large eddy simulation.

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.