Towards a compact high-order method for non-linear hyperbolic systems. I: The Hermite Least-Square Monotone (HLSM) reconstruction

A new Hermite Least-Square Monotone (HLSM) reconstruction to calculate accurately complex flows on non-uniform meshes is presented.     

L-stability of vertex-based MUSCL finite volume schemes on unstructured grids: Simulation of incompressible flows with high density ratios

This work is devoted to the design of multi-dimensional finite volume schemes for solving transport equations on unstructured grids. In the framework of MUSCL vertex-based methods we construct numerical fluxes such that the local maximum property is guaranteed under an explicit Courant–Friedrichs–Levy condition. The method can be naturally completed by adaptive local mesh refinements and it turns out that the mesh generation is less constrained than when using the competitive cell-centered methods. We illustrate the effectiveness of the scheme by simulating variable density incompressible viscous flows. Numerical simulations underline the theoretical predictions and succeed in the computation of high density ratio phenomena such as a water bubble falling in air.     

A new point-locating algorithm under three-dimensional hybrid meshes

When particles are traced individually in a numerical simulation of particle-fluid flow, a point-locating scheme is necessary to relate the coordinates of a given point to the grid cell containing it. A new point-locating scheme combining the advantages of the existing ones is presented in this paper. The new method is first explained in detail. Then its performance under three-dimensional hybrid meshes is compared with the existing approaches in the literature. The results show that the proposed method is efficient, multifunctional and easy to implement in a structured/unstructured fluid flow solver.     

Fluctuation Splitting Schemes for the Compressible and Incompressible Euler and Navier-Stokes Equations

Introduced in the late eighties by Roe, fluctuation splitting (or residual distribution) schemes have recently emerged as a viable alternative to Finite Volume and Finite Element methods for PDE based, fluid dynamics simulations using unstructured meshes. Their application to the numerical approximation of the compressible and incompressible Euler and Navier-Stokes equations is described, emphasizing low Mach number and incompressible applications. The advantages provided by time-preconditioning techniques are discussed and details of the implementation are given.     

CONVERGENCE AND SUPERCONVERGENCE ANALYSIS OF LAGRANGE RECTANGULAR ELEMENTS WITH ANY ORDER ON ARBITRARY RECTANGULAR MESHES

This paper is to study the convergence and superconvergence of rectangular finite elements under anisotropic meshes. By using of the orthogonal expansion method, an anisotropic Lagrange interpolation is presented. The family of Lagrange rectangular elements with all the possible shape function spaces are considered, which cover the Intermediate families, Tensor-product families and Serendipity families. It is shown that the anisotropic interpolation error estimates hold for any order Sobolev norm. We extend the convergence and superconvergence result of rectangular finite elements to arbitrary rectangular meshes in a unified way.     

Sobolev型方程各向异性网格下Wilson元的高精度分析

1引言 Sobolev方程在流体穿过裂缝岩石的渗透理论,土壤中湿气的转移问题,不同介质的热传导问题等许多物理问题中都有广泛的应用,因此已有许多文献研究此方程[2,6,8,9,11]但这些研究都是基于对剖分的正则性条件或拟一致假设[7],即满足hk/pK≤C或h/h≤C,vK∈Th,其中hk,PK分别是一般单元K的最大直径和最大内切圆直径,h=maxhk,h=minhk,C是一个与h无关的正常数.     

THE APPROXIMATE COMPUTATION OF HYPERSINGULAR INTEGRALS ON INTERVAL

In this paper some numerical methods for computing hypersingular integrals on interval are given. Using geometric meshes, these methods lead 10 an exponential convergence in the range of engineering compulation. A numerical example shows their effeclivity and accuracy.     

多尺度有限元法结合分层网格模拟二维奇异摄动的两端边界层问题

通过摄动系数建立分层网格,用多尺度有限元法捕捉对流扩散方程的两端边界层,研究二维奇异摄动模型。基于分层网格并利用多尺度基函数刻画了边界层的微观信息,用有限的计算资源、较短的计算时间,得到了不依赖于摄动系数、一致稳定的模拟结果。     

研究了雷诺数Re=200, 1000, 线速度比$alpha =0.5$,

研究了雷诺数Re=200, 1000, 线速度比$alpha =0.5$,2.0, 4.0, 强迫振荡频率$f_{s}=0.1sim 2.0$情况下的旋转振荡圆柱绕流问题. 通过基于非结构同位网格有限体积法对Navier-Stokes方程进行数值求解. 对流项、扩散项和非恒定项的离散格式均具有二阶精度,利用SIMPLE算法处理压力-速度耦合.计算得到了作用力系数随不同控制参数的变化规律. 通过对升力系数的频谱分析得到自然脱落频率和强迫振荡频率下的作用力振幅. 通过对不同频率作用力幅值的分析,得到频率之间的竞争关系,进而定量地给出了不同尾迹涡脱落模式的分区图.