A High-Order Discontinuous Galerkin Method for 2D Incompressible Flows

In this paper we introduce a high-order discontinuous Galerkin method for two-dimensional incompressible flow in the vorticity stream-function formulation. The momentum equation is treated explicitly, utilizing the efficiency of the discontinuous Galerkin method. The stream function is obtained by a standard Poisson solver using continuous finite elements. There is a natural matching between these two finite element spaces, since the normal component of the velocity field is continuous across element boundaries. This allows for a correct upwinding gluing in the discontinuous Galerkin framework, while still maintaining total energy conservation with no numerical dissipation and total enstrophy stability. The method is efficient for inviscid or high Reynolds number flows. Optimal error estimates are proved and verified by numerical experiments.     

任意不可压多介质流的粒子有限元方法模拟

基于粒子有限元方法(particle finite element method,PFEM),利用细分混合单元的界面识别思想,模拟种类任意多的不可压多介质流问题.对分步算法采用基于有限增量微积分理论的稳定措施,以适应流体特性差异;将混合单元细分为代表单一流体的小单元,进而得到流体间的边界;通过加密边界、控制粒子速度、自动检查穿透来防止粒子穿透外部边界.瑞利-泰勒不稳定性和水柱在空气中倒塌的模拟与已有结果的对比验证了PFEM及界面识别方法的可靠性和准确性.七种流体混合的模拟结果表明PFEM可有效处理任意多种类不相溶流体的混合流动问题.     

Incompressible Flows of an Ideal Fluid¶with Unbounded Vorticity

In this paper we study solutions to the Euler equations of an ideal incompressible fluid in R ⁿ singular at the origin with a finite symmetry group. For an “admissible” class of finite groups we prove a local existence and uniqueness theorem. In even dimensions this theorem covers some symmetric flows with essentially unbounded vorticity. In arbitrary dimension (including n=3) we construct local in time solutions with vorticity that behaves, e.g., like a function of homogeneous degree zero near the origin. The symmetry condition provides necessary additional cancellations and is preserved by the evolution due to uniqueness. Received: 31 March 1999 / Accepted: 10 July 2000     

A Volume of Fluid Based Method for Fluid Flows with Phase Change

This paper presents a numerical method directed towards the simulation of flows with mass transfer due to changes of phase. We use a volume of fluid (VOF) based interface tracking method in conjunction with a mass transfer model and a model for surface tension. The bulk fluids are viscous, conducting, and incompressible. A one-dimensional test problem is developed with the feature that a thin thermal layer propagates with the moving phase interface. This test problem isolates the ability of a method to accurately calculate the thermal layers responsible for driving the mass transfer in boiling flows. The numerical method is tested on this problem and then is used in simulations of horizontal film boiling.     

一种模拟大密度比多相流的混合算法

在扩散界面法(diffuse interface method,DIM)的基础上提出一种能够处理大密度比(large density ratio)的多相流混合算法.流场信息通过格子波尔兹曼方法(lattice Boltzmann method,LBM)获得;相界面通过直接求解Cahn-Hilliard(C-H)方程确定.为保证在大密度比情况下求解界面方程的稳定性,采用二阶迎风格式来离散方程的对流项.通过对Rayleigh-Taylor(R-T)不稳定、液体中的气泡上升及液滴撞击干燥壁面的数值模拟,验证了方法的可行性.     

不可压缩流的并行两水平稳定有限元算法

提出一种数值求解定常不可压缩Stokes方程的并行两水平Grad-div稳定有限元算法。首先在粗网格中求解Grad-div稳定化的全局解, 再在相互重叠的细网格子区域上并行纠正。通过对稳定化参数、粗细网格尺寸恰当的选取, 该方法可得到最优收敛率, 数值结果验证了算法的高效性。     

Simulation of Incompressible Viscous Flows by Local DFD-Immersed Boundary Method

A local domain-free discretization-immersed boundary method (DFD-IBM) is presented in this paper to solve incompressible Navier-Stokes equations in the primitive variable form. Like the conventional immersed boundary method (IBM), the local DFD-IBM solves the governing equations in the whole domain including exterior and interior of the immersed object. The effect of immersed boundary to the surrounding fluids is through the evaluation of velocity at interior and exterior dependent points. To be specific, the velocity at interior dependent points is computed by approximate forms of solution and the velocity at exterior dependent points is set to the wall velocity. As compared to the conventional IBM, the present approach accurately implements the non-slip boundary condition. As a result, there is no flow penetration, which is often appeared in the conventional IBM results. The present approach is validated by its application to simulate incompressible viscous flows around a circular cylinder. The obtained numerical results agree very well with the data in the literature.     

Statistics of Advective Stretching in Three-dimensional Incompressible Flows

We present a method to quantify kinematic stretching in incompressible, unsteady, isoviscous, three-dimensional flows. We extend the method of Kellogg and Turcotte (J. Geophys. Res. 95:421–432, 1990) to compute the axial stretching/thinning experienced by infinitesimal ellipsoidal strain markers in arbitrary three-dimensional incompressible flows and discuss the differences between our method and the computation of Finite Time Lyapunov Exponent (FTLE). We use the cellular flow model developed in Solomon and Mezic (Nature 425:376–380, 2003) to study the statistics of stretching in a three-dimensional unsteady cellular flow. We find that the probability density function of the logarithm of normalised cumulative stretching (log S) for a globally chaotic flow, with spatially heterogeneous stretching behavior, is not Gaussian and that the coefficient of variation of the Gaussian distribution does not decrease with time as t^-\frac12t^{-\frac{1}{2}} . However, it is observed that stretching becomes exponential log S∼t and the probability density function of log S becomes Gaussian when the time dependence of the flow and its three-dimensionality are increased to make the stretching behaviour of the flow more spatially uniform. We term these behaviors weak and strong chaotic mixing respectively. We find that for strongly chaotic mixing, the coefficient of variation of the Gaussian distribution decreases with time as t^-\frac12t^{-\frac{1}{2}} . This behavior is consistent with a random multiplicative stretching process.     

Multimode Coupling Theory for Kelvin-Helmholtz Instability in Incompressible Fluid

A weakly nonlinear model is proposed for multimode Kelvin-Helmholtz instability. The second-order mode coupling formula for Kelvin-Helmholtz instability in two-dimensional incompressible fluid is presented by expanding the perturbation velocity potential to second order. It is found that there is an important resonance in the course of the sum frequency mode coupling but the difference frequency mode coupling does not have. This resonance makes the sum frequency mode coupling process relatively complex. The sum frequency mode coupling is strongly dependent on time especially when the density of the two fluids is adjacent and the difference frequency mode coupling is not.     

An Incompressible Three-Dimensional Multiphase Particle-in-Cell Model for Dense Particle Flows

A three-dimensional, incompressible, multiphase particle-in-cell method is presented for dense particle flows. The numerical technique solves the governing equations of the fluid phase using a continuum model and those of the particle phase using a Lagrangian model. Difficulties associated with calculating interparticle interactions for dense particle flows with volume fractions above 5% have been eliminated by mapping particle properties to an Eulerian grid and then mapping back computed stress tensors to particle positions. A subgrid particle, normal stress model for discrete particles which is robust and eliminates the need for an implicit calculation of the particle normal stress on the grid is presented. Interpolation operators and their properties are defined which provide compact support, are conservative, and provide fast solution for a large particle population. The solution scheme allows for distributions of types, sizes, and density of particles, with no numerical diffusion from the Lagrangian particle calculations. Particles are implicitly coupled to the fluid phase, and the fluid momentum and pressure equations are implicitly solved, which gives a robust solution.     

粘性不可压流体流动问题用直角坐标网格的贴体解法

研究一种新的全贴体的求解粘性不可压流体流动问题的非结构化直角坐标网格方法.该方法在于利用直角坐标网格但通过在边界附近保留不规则控制体,使得算法是完全贴体的.这有别于目前流行的各种非结构化直角坐标网格方法.通过对两个典型流动问题的计算对该数值方法进行验证.对比结果表明,本方法计算的结果与精确解和STAR-CD的结果在一定Re数和网格数时是很接近的,可以满足一定的精度要求,说明该数值计算方法是可行的.还对二维钝头体周围的流场进行了计算,计算的流场与STAR-CD的结果相当吻和,说明该算法还可计算较复杂的流动现象.     

Ghost Fluid方法与双介质可压缩流动计算

应用带有Isobaric修正的GhostFluid方法配合LevelSet方法计算可压缩双介质无粘流动.该方法可以消除计算流体界面时所产生的数值跳动和耗散,且编程上比界面跟踪法简单.应用WENO格式数值求解欧拉方程和LevelSet方程,对由刚性气体状态方程所支配的一二维双介质流动进行数值计算,得到了分辨率较高的计算结果.     

Linear Stability in an Ideal Incompressible Fluid

 We give an explicit construction of approximate eigenfunctions for a linearized Euler operator in dimensions two and three with periodic boundary conditions, and an estimate from below for its spectral bound in terms of an appropriate Lyapunov exponent. As a consequence, we prove that in dimension 2 the spectral and growth bounds for the corresponding group are equal. Therefore, the linear hydrodynamic stability of a steady state for the Euler equations in dimension 2 is equivalent to the fact that the spectrum of the linearized operator is pure imaginary. In dimension 3 we prove the estimate from below for the spectral bound that implies the same equality for every example where the relevant Lyapunov exponents could be effectively computed. For the kinematic dynamo operator describing the evolution of a magnetic field in an ideally conducting incompressible fluid we prove that the growth bound equals the spectral bound in dimensions 2 and 3. Received: 20 May 2002 / Accepted: 5 September 2002 Published online: 10 January 2003 RID="*" ID="*" The first author was partially supported by the Twinning Program of the National Academy of Sciences and National Science Foundation, and by the Research Council and Research Board of the University of Missouri. RID="**" ID="**" The second author was partially supported by the National Science Foundation grant DMS 9876947 and CRDF grant RM1-2084. Acknowledgements. The authors thank Susan Friedlander for useful discussions. Communicated by P. Constantin     

Development of Lattice Boltzmann Flux Solver for Simulation of Incompressible Flows

A lattice Boltzmann flux solver (LBFS) is presented in this work for simulation of incompressible viscous and inviscid flows. The new solver is based on Chapman-Enskog expansion analysis, which is the bridge to link Navier-Stokes (N-S) equations and lattice Boltzmann equation (LBE). The macroscopic differential equations are discretized by the finite volume method, where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers. The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. LBFS is validated by its application to simulate the viscous decaying vortex flow, the driven cavity flow, the viscous flow past a circular cylinder, and the inviscid flow past a circular cylinder. The obtained numerical results compare very well with available data in the literature, which show that LBFS has the second order of accuracy in space, and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.     

界面捕捉Level Set方法的(AMR)数值模拟

在流体力学方程的计算中采用高精度WENO格式,用AMR(adaptive mesh refinement)方法提高流场局部分辨率,在采用Level Set函数标定物质界面的计算中用GFM(ghost fluid method)方法进行界面处理,尝试将AMR技术与界面追踪技术相互融合并应用于数值模拟,对不同的模拟结果进行了比较.     

基于高阶残差量化的光谱二值编码新方法

光谱二值和多值编码技术能够实现目标光谱的快速匹配、识别和分类等应用,但这类量化编码方法会损失大量的光谱细节信息,且不能解码出与原始光谱近似的重构光谱,应用有限。为了解决上述问题,提出一种高阶残差量化的光谱编码新方法HOBC(high-order binary coding)。首先,对光谱向量进行去均值的规范化处理,得到值域为(-1, 1)的光谱序列;然后,求解规范化光谱的±1编码、编码系数和残差（即一阶残差）;基于一阶残差,逐阶解算2至K阶残差的±1编码及其系数;最后得到K个编码序列及其系数,即为HOBC的编码结果。选择典型波谱库数据集,对比光谱0/1二值编码BC01(binary coding with 0 and 1)、光谱分析编码SPAM(spectral analysis manager)、二值/四值混合编码SDFC(spectral derivative feature coding)和DNA四值编码等4种方法,进行了光谱量化编码和解码重构实验,分别统计了光谱形状特征和斜率特征编码的信息熵和存储量、光谱形状特征编码与原始光谱之间的光谱矢量距离SVD (spectral vector distance)、谱间Pearson相关系数SCC (spectral correlation coefficient)和光谱角SAM (spectral angle mapping)。结果表明,在编码存储量上,HOBC的1～4阶编码分别与以上4种编码相等;在编码信息熵上,HOBC的1～2阶编码分别与BC01和SPAM相等,而HOBC的3～4阶编码分别高于SDFC和DNA编码;在SCC上,HOBC1阶编码与BC01相等,而2～4阶编码均分别优于SPAM,SDFC和DNA编码;在SAM方面,HOBC 1～4阶编码均分别明显优于4种对比方法;4种对比方法不能明确解码重构,而HOBC可简便重构出与原始光谱近似的解码序列,且SVD逐阶递减。进一步,基于临泽草地试验站公开光谱数据集,进行了10类地物目标的光谱编码和监督分类实验,实验结果表明,在Kappa系数,总体分类精度和平均分类精度等3种性能评价指标上,HOBC均明显优于4种对比方法,尤其是,HOBC 4阶编码优于原始光谱的分类性能;对样本数量较少且类间相似性较高的难分类地物,HOBC亦具有优于其他算法的鲁棒性。说明HOBC编码在大幅压缩数据量的同时,其编码序列能保留较高的信息量,且具有较高的光谱可分性,可用于光谱高精度快速识别和分类;其解码重构序列与原始光谱序列具有较高的相似性,理论上可适用于目标识别和分类等应用。     

一种投影式激光方向测量方法

为了在激光告警中实现全方位角度探测,且达到一定的探测精度,提出了一种投影式激光方向测量方法.通过计算激光照射在参量给定的遮光板上形成的投影量.求解在探测器上由不同探测象限产生的不同的探测电流的比值.再根据投影面积的比值与角度一一对应的关系可事先编码的原理.由对应程序判断激光入射方向.在软件Matlab仿真中计箅得到产生最小分辨率的四个极限角度分别为(0°,90°)、(33°,21°)、(33°,69°)、(45°,90°),并由实验得到最小分辨电压为0.05 V.通过仿真计算与实验数据证明该设计可达到最小识别角度1°的探测要求.     

Lattice-Boltzmann Simulations of Fluid Flows in MEMS

The lattice Boltzmann model is a simplified kinetic method based on the particle distribution function. We use this method to simulate problems in MEMS, in which the velocity slip near the wall plays an important role. It is demonstrated that the lattice Boltzmann method can capture the fundamental behaviors in micro-channel flow, including velocity slip, nonlinear pressure drop along the channel and mass flow rate variation with Knudsen number. The Knudsen number dependence of the position of the vortex center and the pressure contour in micro-cavity flows is also demonstrated.     

一种新颖的变步长光伏系统MPPT控制方法

太阳能光伏阵列的输出功率随外界环境因素的变化而变化,为了能高效地利用太阳能电池,需对光伏阵列进行最大功率点跟踪(简称MPPT);针对定步长电导增量控制法存在的缺点,提出一种新颖的变步长电导增量法;该方法通过瞬时功率相对电流的变化值确定MPPT步长大小,使光伏系统MPPT快速准确且无震荡;经Matlab仿真及实验证明:该方法能快速准确地跟踪外部环境变化,并能保证系统的稳定性。