A Fourier spectral method for the Navier–Stokes equations with volume penalization for moving solid obstacles

An efficient numerical scheme to compute flows past rigid solid bodies moving through viscous incompressible fluid is presented. Solid obstacles of arbitrary shape are taken into account using the volume penalization method to impose no-slip boundary condition. The 2D Navier–Stokes equations, written in the vorticity-streamfunction formulation, are discretized using a Fourier pseudo-spectral scheme. Four different time discretization schemes of the penalization term are proposed and compared. The originality of the present work lies in the implementation of time-dependent penalization, which makes the above method capable of solving problems where the obstacle follows an arbitrary motion. Fluid–solid coupling for freely falling bodies is also implemented. The numerical method is validated for different test cases: the flow past a cylinder, Couette flow between rotating cylinders, sedimentation of a cylinder and a falling leaf with elliptical shape.     

计算流体成像技术重构方法及其应用

建立系统的计算流体成像(CFI)技术,并论证其为流动数值模拟与风洞实验相互佐证的重要手段。采用纯数值模拟技术重构得到流场计算阴影、纹影以及干涉图,主要技术包括复杂湍流模拟技术、三线性插值技术和流动图像重构技术等。同时将CFI技术应用于激波管内绕障碍物三维复杂流动计算图像重构,并与双向全息实验结果进行了比较分析。研究证明该CFI技术可以用于复杂几何边界、分区域、非定常流动图像重构,不仅能够为流场结构特征分析提供重要参考依据,而且能够用于数值模拟与风洞实验结果的比较印证。     

Shock-Capturing and Front-Tracking Methods for Granular Avalanches

Shock formations are observed in granular avalanches when supercritical flow merges into a region of subcritical flow. In this paper we employ a shock-capturing numerical scheme for the one-dimensional Savage–Hutter theory of granular flow to describe this phenomenon. A Lagrangian moving mesh scheme applied to the nonconservative form of the equations reproduces smooth solutions of these free boundary problems very well, but fails when shocks are formed. A nonoscillatory central (NOC) difference scheme with TVD limiter or WENO cell reconstruction for the conservative equations is therefore introduced. For the avalanche free boundary problems it must be combined with a front-tracking method, developed here, to properly describe the margin evolution. It is found that this NOC scheme combined with the front-tracking module reproduces both the shock wave and the smooth solution accurately. A piecewise quadratic WENO reconstruction improves the smoothness of the solution near local extrema. The schemes are checked against exact solutions for (1) an upward moving shock wave, (2) the motion of a parabolic cap down an inclined plane, and (3) the motion of a parabolic cap down a curved slope ending in a flat run-out region, where a shock is formed as the avalanche comes to a halt.     

SVD–GFD scheme to simulate complex moving body problems in 3D space

The present paper presents a hybrid meshfree-and-Cartesian grid method for simulating moving body incompressible viscous flow problems in 3D space. The method combines the merits of cost-efficient and accurate conventional finite difference approximations on Cartesian grids with the geometric freedom of generalized finite difference (GFD) approximations on meshfree grids. Error minimization in GFD is carried out by singular value decomposition (SVD). The Arbitrary Lagrangian–Eulerian (ALE) form of the Navier–Stokes equations on convecting nodes is integrated by a fractional-step projection method. The present hybrid grid method employs a relatively simple mode of nodal administration. Nevertheless, it has the geometrical flexibility of unstructured mesh-based finite-volume and finite element methods. Boundary conditions are precisely implemented on boundary nodes without interpolation. The present scheme is validated by a moving patch consistency test as well as against published results for 3D moving body problems. Finally, the method is applied on low-Reynolds number flapping wing applications, where large boundary motions are involved. The present study demonstrates the potential of the present hybrid meshfree-and-Cartesian grid scheme for solving complex moving body problems in 3D.     

A high-order WENO-Z finite difference based particle-source-in-cell method for computation of particle-laden flows with shocks

A high-order particle-source-in-cell (PSIC) algorithm is presented for the computation of the interaction between shocks, small scale structures, and liquid and/or solid particles in high-speed engineering applications. The improved high-order finite difference weighted essentially non-oscillatory (WENO-Z) method for solution of the hyperbolic conservation laws that govern the shocked carrier gas flow, lies at the heart of the algorithm. Finite sized particles are modeled as points and are traced in the Lagrangian frame. The physical coupling of particles in the Lagrangian frame and the gas in the Eulerian frame through momentum and energy exchange, is numerically treated through high-order interpolation and weighing. The centered high-order interpolation of the fluid properties to the particle location is shown to lead to numerical instability in shocked flow. An essentially non-oscillatory interpolation (ENO) scheme is devised for the coupling that improves stability. The ENO based algorithm is shown to be numerically stable and to accurately capture shocks, small flow features and particle dispersion. Both the carrier gas and the particles are updated in time without splitting with a third-order Runge–Kutta TVD method. One and two-dimensional computations of a shock moving into a particle cloud demonstrates the characteristics of the WENO-Z based PSIC method (PSIC/WENO-Z). The PSIC/WENO-Z computations are not only in excellent agreement with the numerical simulations with a third-order Rusanov based PSIC and physical experiments in [V. Boiko, V.P. Kiselev, S.P. Kiselev, A. Papyrin, S. Poplavsky, V. Fomin, Shock wave interaction with a cloud of particles, Shock Waves, 7 (1997) 275–285], but also show a significant improvement in the resolution of small scale structures. In two-dimensional simulations of the Mach 3 shock moving into forty thousand bronze particles arranged in the shape of a rectangle, the long time accuracy of the high-order method is demonstrated. The fifth-order PSIC/WENO-Z method with the fifth-order ENO interpolation scheme improves the small scale structure resolution over the third-order PSIC/WENO-Z method with a second-order central interpolation scheme. Preliminary analysis of the particle interaction with the flow structures shows that sharp particle material arms form on the side of the rectangular shape. The arms initially shield the particles from the accelerated flow behind the shock. A reflected compression wave, however, reshocks the particle arm from the shielded area and mixes the particles.     

A New Method for Analysis of the Fluid Interaction with a Deformable Membrane

The lattice Boltzmann cellular automaton method has been successfully extended for analysis of fluid interactions with a deformable membrane or web. The hydrodynamic forces on the solid web are obtained through computation of the fluid flow stress at the moving boundary using the lattice Boltzmann method. Analysis of solid boundary deformation or vibration due to hydrodynamic force is based on Newtonian dynamics and a molecular dynamic type approach.     

中厚板弯曲问题的Kriging插值无网格法

基于Reissner-Mindlin板弯曲理论,将Kriging插值无网格法应用于中厚板弯曲问题,推导相应的离散方程.该方法可以只依赖于一组离散的节点建立试函数,有效地避免了复杂的网格划分和网格畸变的影响.相对于无网格法中常用的移动最小二乘近似而言,滑动Kriging插值法的形函数满足Kronecker delta函数性质,可以直接施加本质边界条件.算例分析表明,用Kriging插值无网格法分析中厚板弯曲问题,具有效率高,精度高和易于实现等优点.     

Solving Navier–Stokes equation for flow past cylinders using single-block structured and overset grids

Many complex fluid motions are driven by physical processes of instability, transition and turbulence dependent upon nonlinear mechanisms. Here, we solve the flow past cylinder(s) using single-block structured and overset grids by computing Navier–Stokes equation in two-dimensions. The suitability of a compact scheme in discretizing convection and diffusion terms are investigated first by looking at relevant numerical properties. Also, for the overset grid method, one of the methods is identified that shows the best results in minimizing interpolation error at sub-domain boundaries for an analytical test function. We provide extensive comparisons with experimental and other computational results for flow past a single cylinder, utilizing both single-block structured and Chimera or overset grids. Apart from showing instability of this flow calculated by these methods, we also compare the computed vorticity and velocity data using these two grids by employing the proper orthogonal decomposition (POD). We have analyzed and developed an overset grid method with compact scheme that does not need any filtering to control error. This has been ascertained by performing POD analysis. To show that the developed method is capable of handling complex geometries, we have computed flow past two cylinders in side-by-side arrangement. Results obtained capture the known flow characteristics for this arrangement well using relatively fewer number of grid points.     

A simple and efficient direct forcing immersed boundary framework for fluid–structure interactions

A direct forcing immersed boundary framework is presented for the simple and efficient simulation of strongly coupled fluid–structure interactions. The immersed boundary method developed by Yang and Balaras [J. Yang, E. Balaras, An embedded-boundary formulation for large-eddy simulation of turbulent flows interacting with moving boundaries, J. Comput. Phys. 215 (1) (2006) 12–40] is greatly simplified by eliminating several complicated geometric procedures without sacrificing the overall accuracy. The fluid–structure coupling scheme of Yang et al. [J. Yang, S. Preidikman, E. Balaras, A strongly-coupled, embedded-boundary method for fluid–structure interactions of elastically mounted rigid bodies, J. Fluids Struct. 24 (2008) 167–182] is also significantly expedited by moving the fluid solver out of the predictor–corrector iterative loop without altering the strong coupling property. Central to these improvements are the reformulation of the field extension strategy and the evaluation of fluid force and moment exerted on the immersed bodies, by taking advantage of the direct forcing idea in a fractional-step method. Several cases with prescribed motions are examined first to validate the simplified field extension approach. Then, a variety of strongly coupled fluid–structure interaction problems, including vortex-induced vibrations of a circular cylinder, transverse and rotational galloping of rectangular bodies, and fluttering and tumbling of rectangular plates, are computed. The excellent agreement between the present results and the reference data from experiments and other simulations demonstrates the accuracy, simplicity, and efficiency of the new method and its applicability in a wide range of complicated fluid–structure interaction problems.     

振动NACA0012翼型的移动网格数值模拟

利用流体力学方程的积分形式给出非结构移动网格上离散格式,利用自适应移动网格方法移动网格,进而得到网格速度.对振动Naca0012翼型问题,分三种类型确定网格速度,再结合Riemann问题的解法器构造数值通量,得到移动网格单元上新的物理量.数值实验表明这种格式同时具有高效、高分辨率的特点.     

运动激波和气泡串相互作用的初步数值模拟

通过对激波和流体界面相互作用诱导的大变形界面演化的数值模拟,验证Level set方法精确模拟多个流体界面的有效性.采用2阶迎风TVD求解欧拉方程得到流场解,采用5阶WENO求解Level set方程追踪多流体界面,采用GFM方法处理流体内界面.利用文[1]的计算结果校核本文程序.在此基础上,对运动激波和气泡串相互作用过程进行了初步数值模拟,得到了不同时刻运动激波和圆管内的两个气泡作用后的演化图象,包括压力和密度等值线分布.计算结果表明:针对推广后的多界面Level set方程,该方法仍可高质量地捕捉多个流体界面.     

A differentially interpolated direct forcing immersed boundary method for predicting incompressible Navier–Stokes equations in time-varying complex geometries

A dispersion-relation-preserving dual-compact scheme developed in Cartesian grids is applied together with the immersed boundary method to solve the flow equations in irregular and time-varying domains. The artificial momentum forcing term applied at certain points in cells containing fluid and solid allows an imposition of velocity condition to account for the motion of solid body. We develop in this study a differential-based interpolation scheme which can be easily extended to three-dimensional simulation. The results simulated from the proposed immersed boundary method agree well with other numerical and experimental results for the chosen benchmark problems. The accuracy and fidelity of the IB flow solver developed to predict flows with irregular boundaries are therefore demonstrated.     

A lattice Boltzmann method for immiscible multiphase flow simulations using the level set method

We consider the lattice Boltzmann method for immiscible multiphase flow simulations. Classical lattice Boltzmann methods for this problem, e.g. the colour gradient method or the free energy approach, can only be applied when density and viscosity ratios are small. Moreover, they use additional fields defined on the whole domain to describe the different phases and model phase separation by special interactions at each node. In contrast, our approach simulates the flow using a single field and separates the fluid phases by a free moving interface. The scheme is based on the lattice Boltzmann method and uses the level set method to compute the evolution of the interface. To couple the fluid phases, we develop new boundary conditions which realise the macroscopic jump conditions at the interface and incorporate surface tension in the lattice Boltzmann framework. Various simulations are presented to validate the numerical scheme, e.g. two-phase channel flows, the Young–Laplace law for a bubble and viscous fingering in a Hele-Shaw cell. The results show that the method is feasible over a wide range of density and viscosity differences.     

Lattice Boltzmann method for three-dimensional moving particles in a Newtonian fluid

A lattice Boltzmann method is developed to simulate three-dimensional solid particle motions in fluids. In the present model, a uniform grid is used and the exact spatial location of the physical boundary of the suspended particles is determined using an interpolation scheme. The numerical accuracy and efficiency of the proposed lattice Boltzmann method is demonstrated by simulating the sedimentation of a single sphere in a square cylinder. Highly accurate simulation results can be achieved with few meshes, compared with the previous lattice Boltzmann methods. The present method is expected to find applications on the flow systems with moving boundaries, such as the blood flow in distensible vessels, the particle－flow interaction and the solidification of alloys.     

激波和单列水柱、气水固多界面相互作用的数值模拟

研究可压缩多介质流场的激波和多介质界面相互作用问题.在Descartes固定网格采用level-set方法追踪界面,气/气界面边界条件处理采用OGFM方法,采用修正的rGFM方法提高气/水和气/固界面处构造Riemann问题精度,将Riemann近似解得到的界面参数外推到两侧真实和虚拟流体,采用五阶WENO方法求解流场Euler方程和界面level-set方程,给出不同时刻流场数值纹影图像.结果表明：在可压缩流场嵌入固体和水、气体等目标,本文方法可较精确地分辨平面运动激波和单列水柱及包含气/气、气/水和气/固等界面作用后产生的复杂激波结构.和传统的分区与贴体变换方法不同,为Descartes网格包含多介质界面复杂流场计算提供新途径.     

有激波的一维不定常气流的新的数值解法

解有激波的气体力学问题的数值解法,主要有特征线法和有限差分法两类。特征线法一般能给出高的精度,但当激波很弱、与特征线几乎平行时就需特殊处理^[1]。有限差分法的研究和发展更广,有人工粘性法^[2]、激波捕捉法^[8]、分离奇性法^[1]等处理激波的方案,它们又各有其特殊的技巧和问题。     

Remapping-free ALE-type kinetic method for flow computations

Based on the integral form of the fluid dynamic equations, a finite volume kinetic scheme with arbitrary control volume and mesh velocity is developed. Different from the earlier unified moving mesh gas-kinetic method [C.Q. Jin, K. Xu, An unified moving grid gas-kinetic method in Eulerian space for viscous flow computation, J. Comput. Phys. 222 (2007) 155–175], the coupling of the fluid equations and geometrical conservation laws has been removed in order to make the scheme applicable for any quadrilateral or unstructured mesh rather than parallelogram in 2D case. Since a purely Lagrangian method is always associated with mesh entangling, in order to avoid computational collapsing in multidimensional flow simulation, the mesh velocity is constructed by considering both fluid velocity (Lagrangian methodology) and diffusive velocity (Regenerating Eulerian mesh function). Therefore, we obtain a generalized Arbitrary-Lagrangian–Eulerian (ALE) method by properly designing a mesh velocity instead of re-generating a new mesh after distortion. As a result, the remapping step to interpolate flow variables from old mesh to new mesh is avoided. The current method provides a general framework, which can be considered as a remapping-free ALE-type method. Since there is great freedom in choosing mesh velocity, in order to improve the accuracy and robustness of the method, the adaptive moving mesh method [H.Z. Tang, T. Tang, Adaptive mesh methods for one-and two-dimensional hyperbolic conservation laws, SIAM J. Numer. Anal. 41 (2003) 487–515] can be also used to construct a mesh velocity to concentrate mesh to regions with high flow gradients.     

An optimized spectral difference scheme for CAA problems

In the implementation of spectral difference (SD) method, the conserved variables at the flux points are calculated from the solution points using extrapolation or interpolation schemes. The errors incurred in using extrapolation and interpolation would result in instability. On the other hand, the difference between the left and right conserved variables at the edge interface will introduce dissipation to the SD method when applying a Riemann solver to compute the flux at the element interface. In this paper, an optimization of the extrapolation and interpolation schemes for the fourth order SD method on quadrilateral element is carried out in the wavenumber space through minimizing their dispersion error over a selected band of wavenumbers. The optimized coefficients of the extrapolation and interpolation are presented. And the dispersion error of the original and optimized schemes is plotted and compared. An improvement of the dispersion error over the resolvable wavenumber range of SD method is obtained. The stability of the optimized fourth order SD scheme is analyzed. It is found that the stability of the 4th order scheme with Chebyshev–Gauss–Lobatto flux points, which is originally weakly unstable, has been improved through the optimization. The weak instability is eliminated completely if an additional second order filter is applied on selected flux points. One and two dimensional linear wave propagation analyses are carried out for the optimized scheme. It is found that in the resolvable wavenumber range the new SD scheme is less dispersive and less dissipative than the original scheme, and the new scheme is less anisotropic for 2D wave propagation. The optimized SD solver is validated with four computational aeroacoustics (CAA) workshop benchmark problems. The numerical results with optimized schemes agree much better with the analytical data than those with the original schemes.     

低渗透多孔介质非达西渗流动边界界面追踪

基于低渗透多孔介质非达西不稳定渗流的动边界数学模型,推导动边界移动速度的微分表达式,揭示动边界移动速度与动边界上地层压力关于径向距离的二次导数成正比;由此利用拉格朗日三点插值公式求得动边界附近控制方程的有限差分格式,并对下一时刻动边界的精确位置进行追踪.有限差分方法的数值结果表明界面追踪法可较好地反映低渗透多孔介质非达西不稳定渗流动边界的移动规律.     

湍流模式对转捩的初步研究

采用二阶时间精度的LU-SGS-τTS目的进行时间推进求解Navier-Stokes方程组,用基于Roe的六阶对称TVD格式离散对流项,并采用Baldwin-Lomax和Spalart-Allmaras模式对平板边界层、Aerospatial A-profile翼型低速大迎角及RAE-2822翼型跨声速激波/边界层干扰流动进行数值分析,自动捕捉流动的转捩过程,获取了有意义的结果.