An h-adaptive solution of the spherical blast wave problem

Shock waves and contact discontinuities usually appear in compressible flows, requiring a fine mesh in order to achieve an acceptable accuracy of the numerical solution. The usage of a mesh adaptation strategy is convenient as uniform refinement of the whole mesh becomes prohibitive in three-dimensional (3D) problems. An unsteady h-adaptive strategy for unstructured finite element meshes is introduced. Non-conformity of the refined mesh and a bounded decrease in the geometrical quality of the elements are some features of the refinement algorithm. A 3D extension of the well-known refinement constraint for 2D meshes is used to enforce a smooth size transition among neighbour elements with different levels of refinement. A density-based gradient indicator is used to track discontinuities. The solution procedure is partially parallelised, i.e. the inviscid flow equations are solved in parallel with a finite element SUPG formulation with shock capturing terms while the adaptation of the mesh is sequentially performed. Results are presented for a spherical blast wave driven by a point-like explosion with an initial pressure jump of 10⁵ atmospheres. The adapted solution is compared to that computed on a fixed mesh. Also, the results provided by the theory of self-similar solutions are considered for the analysis. In this particular problem, adapting the mesh to the solution accounts for approximately 4% of the total simulation time and the refinement algorithm scales almost linearly with the size of the problem.     

Numerical simulation of compressible two-phase flow using a diffuse interface method

In this article, a high-resolution diffuse interface method is investigated for simulation of compressible two-phase gas–gas and gas–liquid flows, both in the presence of shock wave and in flows with strong rarefaction waves similar to cavitations. A Godunov method and HLLC Riemann solver is used for discretization of the Kapila five-equation model and a modified Schmidt equation of state (EOS) is used to simulate the cavitation regions. This method is applied successfully to some one- and two-dimensional compressible two-phase flows with interface conditions that contain shock wave and cavitations. The numerical results obtained in this attempt exhibit very good agreement with experimental results, as well as previous numerical results presented by other researchers based on other numerical methods. In particular, the algorithm can capture the complex flow features of transient shocks, such as the material discontinuities and interfacial instabilities, without any oscillation and additional diffusion. Numerical examples show that the results of the method presented here compare well with other sophisticated modeling methods like adaptive mesh refinement (AMR) and local mesh refinement (LMR) for one- and two-dimensional problems.     

一种新的求解对流占优问题的自适应网格细化方法

在均匀网格上求解对流占优问题时,往往会产生数值震荡现象,因此需要局部加密网格来提高解的精度。针对对流占优问题,设计了一种新的自适应网格细化算法。该方法采用流线迎风SUPG(Petrov-Galerkin)格式求解对流占优问题,定义了网格尺寸并通过后验误差估计子修正来指导自适应网格细化,以泡泡型局部网格生成算法BLMG为网格生成器,通过模拟泡泡在区域中的运动得到了高质量的点集。与其他自适应网格细化方法相比,该方法可在同一框架内实现网格的细化和粗化,同时在所有细化层得到了高质量的网格。数值算例结果表明,该方法在求解对流占优问题时具有更高的数值精度和更好的收敛性。     

双曲偏微分方程的局部伪弧长方法研究

重点研究了局部伪弧长方法在处理偏微分方程,尤其是双曲型偏微分方程出现激波间断的奇异性问题,对比分析了全局伪弧长方法空间转化的形式及其网格自适应的性质。为提高求解效率,提出了局部伪弧长方法,利用激波间断的性质,给出了判断奇异点位置以及模板选择的方法,涉及如何处理激波振荡,如何引入弧长参数,以及怎样求解间断等问题。通过数值算例验证了局部伪弧长在激波捕捉和追踪方面的可行性,通过比较局部伪弧长方法与Godunov方法处理不同初值条件的双曲问题,显示出局部伪弧长方法处理双曲偏微分方程的优越性,为伪弧长方法应用到物理问题奠定基础。     

A simple shock‐capturing technique for high‐order discontinuous Galerkin methods

This article presents a novel shock‐capturing technique for the discontinuous Galerkin (DG) method. The technique is designed for compressible flow problems, which are usually characterized by the presence of strong shocks and discontinuities. The inherent structure of standard DG methods seems to suggest that they are especially adapted to capture shocks because of the numerical fluxes based on suitable approximate Riemann solvers, which, in practice, introduces some stabilization. However, the usual numerical fluxes are not sufficient to stabilize the solution in the presence of shocks for large high‐order elements. Here, a new basis of shape functions is introduced. It has the ability to change locally between a continuous or discontinuous interpolation depending on the smoothness of the approximated function. In the presence of shocks, the new discontinuities inside an element introduce the required stabilization because of numerical fluxes. Large high‐order elements can therefore be used and shocks captured within a single element, avoiding adaptive mesh refinement and preserving the locality and compactness of the DG scheme. Several numerical examples for transonic and supersonic flows are studied to demonstrate the applicability of the proposed approach. Copyright © 2011 John Wiley & Sons, Ltd.     

Evaluation of adaptive mesh refinement and coarsening for the computation of compressible flows on unstructured meshes

The compressible gas flows of interest to aerospace applications often involve situations where shock and expansion waves are present. Decreasing the characteristic dimension of the computational cells in the vicinity of shock waves improves the quality of the computed flows. This reduction in size may be accomplished by the use of mesh adaption procedures. In this paper an analysis is presented of an adaptive mesh scheme developed for an unstructured mesh finite volume upwind computer code. This scheme is tailored to refine or coarsen the computational mesh where gradients of the flow properties are respectively high or low. The refinement and coarsening procedures are applied to the classical gas dynamic problems of the stabilization of shock waves by solid bodies. In particular, situations where oblique shock waves interact with an expansion fan and where bow shocks arise around solid bodies are considered. The effectiveness of the scheme in reducing the computational time, while increasing the solution accuracy, is assessed. It is shown that the refinement procedure alone leads to a number of computational cells which is 20% larger than when alternate passes of refinement and coarsening are used. Accordingly, a reduction of computational time of the same order of magnitude is obtained. Copyright © 2005 John Wiley & Sons, Ltd.     

基于适合分析T样条的高阶数值流形方法

数值流形方法是一种非常灵活的数值计算方法,连续体的有限单元方法和块体系统的非连续变形分析方法只是这一数值方法的特例.数值流形方法中高阶位移函数的构造可通过提高权函数的阶次来实现,这种方法往往需要沿单元边界配置适当的边内节点,这些结点的出现增加了前处理的复杂性,特别是对于大型复杂的空间问题.另一方面,在数值流形方法中可通过缩小单元尺寸(h加密)来提高求解精度.当模拟裂纹扩展时,这种细化策略可用来克服裂纹尖端的奇异性.一个传统的解决方案是细化整个网格,但这会导致计算效率的显著降低.将适合分析的T样条(analysis-suitable T-spline,AST)引入数值流形方法中来建立高阶数值流形方法的分析格式,有效的避免了该问题的出现.AST样条基函数具有线性无关,单位分解,局部加密等许多重要性质,使得其非常适合用于工程设计及分析.在引入AST样条后,可通过改变数学覆盖的构造形式建立不同阶次的数值流形方法分析格式;AST样条自身的局部加密性质也使得数值流形方法中的数学网格局部加密更容易实现.算例结果表明:随着AST样条基函数阶次的提高,数值流形方法的计算结果有了明显的改善;基于AST样条基函数的数值流形方法在保持计算精度的前提下降低了自由度的数量.     

求解弱不连续问题的p型自适应有限元方法

在实际工程计算中,存在大量的弱不连续问题,如含夹杂问题。利用通常的有限元方法,为确保界面上各点满足给定高精度,往往需要采用全域网格加密或全域提高单元阶次的方法,这将会导致计算机的物理内存和CPU时间的剧烈增长。p-型自适应有限元方法是一种能通过自适应分析逐步增加单元阶次以改善计算精度的数值方法。本文,我们针对弱不连续问题设计了相应的p-型自适应有限元方法,重点讨论了容许误差控制标准对界面上各点计算结果的影响,并对几类典型的弱不连续问题进行了数值计算与模拟。数值结果表明,本文设计的p-型自适应有限元方法对求解弱不连续问题是非常有效的,用较少的单元得到精度可靠的数值结果,可大大提高其有限元分析效率。     

Numerical simulation of compressible multifluid flows using an adaptive positivity-preserving RKDG-GFM approach

The Runge-Kutta discontinuous Galerkin method together with a refined real-ghost fluid method is incorporated into an adaptive mesh refinement environment for solving compressible multifluid flows, where the level set method is used to capture the moving material interface. To ensure that the Riemann problem is exactly along the normal direction of the material interface, a simple and efficient modification is introduced into the original real-ghost fluid method for constructing the interfacial Riemann problem, and the initial conditions of the Riemann problem are obtained directly from the solution polynomials of the discontinuous Galerkin finite element space. In addition, a positivity-preserving limiter is introduced into the Runge-Kutta discontinuous Galerkin method to suppress the failure of preserving positivity of density or pressure for the problems involving strong shock wave or shock interaction with material interface. For interfacial cells in adaptive mesh refinement, the data transfer between different grid levels is achieved by using a L² projection approach along with the least squares fitting. Various numerical cases, including multifluid shock tubes, underwater explosions, and shock-induced collapse of a underwater air bubble, are computed to assess the capability of the present adaptive positivity-preserving RKDG-GFM approach, and the simulated results show that the present approach is quite robust and can provide relatively reasonable results across a wide variety of flow regimes, even for problems involving strong shock wave or shock wave impacting high acoustic impedance mismatch material interface.     

A high-order Runge-Kutta discontinuous Galerkin method with a subcell limiter on adaptive unstructured grids for two-dimensional compressible inviscid flows

A robust, adaptive unstructured mesh refinement strategy for high-order Runge-Kutta discontinuous Galerkin method is proposed. The present work mainly focuses on accurate capturing of sharp gradient flow features like strong shocks in the simulations of two-dimensional inviscid compressible flows. A posteriori finite volume subcell limiter is employed in the shock-affected cells to control numerical spurious oscillations. An efficient cell-by-cell adaptive mesh refinement is implemented to increase the resolution of our simulations. This strategy enables to capture strong shocks without much numerical dissipation. A wide range of challenging test cases is considered to demonstrate the efficiency of the present adaptive numerical strategy for solving inviscid compressible flow problems having strong shocks.     

An adaptive enrichment algorithm for advection‐dominated problems

We are interested in developing a numerical framework well suited for advection–diffusion problems when the advection part is dominant. In that case, given Dirichlet type boundary condition, it is well known that a boundary layer develops. To resolve correctly this layer, standard methods consist in increasing the mesh resolution and possibly increasing the formal accuracy of the numerical method. In this paper, we follow another path: we do not seek to increase the formal accuracy of the scheme but, by a careful choice of finite element, to lower the mesh resolution in the layer. Indeed the finite element representation we choose is locally the sum of a standard one plus an enrichment. This paper proposes such a method and with several numerical examples, we show the potential of this approach. In particular, we show that the method is not very sensitive to the choice of the enrichment and develop an adaptive algorithm to automatically choose the enrichment functions.Copyright © 2012 John Wiley & Sons, Ltd.     

计算动力学中的伪弧长方法研究

动力学问题通常采用微分方程来描绘,但由于工程实际问题的复杂性,微分方程模型常伴随着解的不连续性、刚性或激波间断奇异性特点,传统方法很难求解,奇异性问题是计算动力学难点,同时也是国内外学者研究的热点.伪弧长数值算法是针对计算动力学中的奇异性问题所提出的,其基本思想为通过在解曲线上引入伪弧长参数,并增加一个约束方程,在伪弧长参数作用下,使得原始离散单元发生扭曲形变,从而达到消除或减弱奇异性的目的.本文首先介绍伪弧长方法求解定常对流-扩散方程的奇异性问题,并提出针对双曲守恒定律的局部伪弧长算法,其思想在于首先通过间断解的梯度变换来确定强间断所处位置,进而通过局部网格点重构以及数值修正来达到强间断处奇异性消除与降低的目的.针对高维问题,提出全局伪弧长方法,通过对整个计算区域内的网格点进行重构,使得所有网格点向奇异间断点处移动,从而降低间断点的影响域,达到降低奇异性的目的.重点讨论了三维全局伪弧长算法问题的计算难点,即三维空间网格扭曲大变形导致的数值算法不收敛,并提出在算法设计过程中采用分块重构与整体计算相结合的策略,实现了三维空间中的伪弧长数值算法,最后通过数值实验来验证伪弧长算法对于奇异性问题的有效性.     

Moving grid method for numerical simulation of stratified flows

We develop a second‐order accurate Navier–Stokes solver based on r‐adaptivity of the underlying numerical discretization. The motion of the mesh is based on the fluid velocity field; however, certain adjustments to the Lagrangian velocities are introduced to maintain quality of the mesh. The adjustments are based on the variational approach of energy minimization to redistribute grid points closer to the areas of rapid solution variation. To quantify the numerical diffusion inherent to each method, we monitor changes in the background potential energy, computation of which is based on the density field. We demonstrate on a standing interfacial gravity wave simulation how using our method of grid evolution decreases the rate of increase of the background potential energy compared with using the same advection scheme on the stationary grid. To further highlight the benefit of the proposed moving grid method, we apply it to the nonhydrostatic lock‐exchange flow where the evolution of the interface is more complex than in the standing wave test case. Naive grid evolution based on the fluid velocities in the lock‐exchange flow leads to grid tangling as Kelvin–Helmholtz billows develop at the interface. This is remedied by grid refinement using the variational approach. Copyright © 2012 John Wiley & Sons, Ltd.     

An adaptive numerical scheme for solving incompressible 2‐phase and free‐surface flows

In this paper, we present a numerical scheme for solving 2‐phase or free‐surface flows. Here, the interface/free surface is modeled using the level‐set formulation, and the underlying mesh is adapted at each iteration of the flow solver. This adaptation allows us to obtain a precise approximation for the interface/free‐surface location. In addition, it enables us to solve the time‐discretized fluid equation only in the fluid domain in the case of free‐surface problems. Fluids here are considered incompressible. Therefore, their motion is described by the incompressible Navier‐Stokes equation, which is temporally discretized using the method of characteristics and is solved at each time iteration by a first‐order Lagrange‐Galerkin method. The level‐set function representing the interface/free surface satisfies an advection equation that is also solved using the method of characteristics. The algorithm is completed by some intermediate steps like the construction of a convenient initial level‐set function (redistancing) as well as the construction of a convenient flow for the level‐set advection equation. Numerical results are presented for both bifluid and free‐surface problems.     

Adaptive h-refinement on 3D unstructured grids for transient problems

An adaptive finite element scheme for transient problems is presented. The classic h-enrichment/coarsening is employed in conjunction with a tetrahedral finite element discretization in three dimensions. A mesh change is performed every n time steps, depending on the Courant number employed and the number of ‘protective layers’ added ahead of the refined region. In order to simplify the refinement/coarsening logic and to be as fast as possible, only one level of refinement/coarsening is allowed per mesh change. A high degree of vectorizability has been achieved by pre-sorting the elements and then performing the refinement/coarsening groupwise according to the case at hand. Further reductions in CPU requirements arc realized by optimizing the identification and sorting of elements for refinement and deletion. The developed technology has been used extensively for shock-shock and shock-object interaction runs in a production mode. A typical example of this class of problems is given.     

Adaptive mesh refinement for simulation of thin film flows

We present a robust and accurate numerical method for simulating gravity-driven, thin-film flow problems. The convection term in the governing equation is treated by a semi-implicit, essentially non-oscillatory scheme. The resulting nonlinear discrete equation is solved using a nonlinear full approximation storage multigrid algorithm with adaptive mesh refinement techniques. A set of representative numerical experiments are presented. We show that the use of adaptive mesh refinement reduces computational time and memory compared to the equivalent uniform mesh results. Our simulation results are consistent with previous experimental observations.     

基于XFEM与自适应网格的非均质材料裂纹扩展模拟

针对含有间断的非均匀材料的断裂问题，本文将虚节点多边形单元的形函数引入到扩展有限元（XFEM）中，提出了一种基于四叉树结构的动态网格细化方法，该方法可对间断面附近单元实现可调控的多层级细化，特别是对于裂纹扩展问题，可实现裂尖附近单元的动态网格细化与粗化。基于以上网格细化方法，本文提出了针对非均匀材质裂纹扩展问题的计算方法VP-XFEM。为验证算法的准确性与计算效率，针对含有孔洞及材料界面的断裂问题，本文给出了相应的算例。结果显示，与传统的一致性网格的XFEM相比，VP-XFEM能够明显改善计算精度与计算效率。     

On the steady solution of linear advection problems in steady recirculating flows

Numerical modelling of non-Newtonian flows typically involves the coupling between equations of motion characterized by an elliptic behaviour, and the fluid constitutive equation, which is an advection equation linked to the fluid history. In this paper we prove that linear steady advection problems in steady recirculating flows have only one solution when the kinematics differs from a rigid motion. We also give a numerical procedure to determine this steady solution. We will describe this numerical procedure for two linear models the first will be the SFRT flow model and the second will be a simplified linear formulation of the Pom–Pom viscoelastic model.     

A high‐order element based adaptive mesh refinement strategy for three‐dimensional unstructured grid

Adaptive mesh refinement (AMR) shows attractive properties in automatically refining the flow region of interest, and with AMR, better prediction can be obtained with much less labor work and cost compared to manually remeshing or the global mesh refinement. Cartesian AMR is well established; however, AMR on hybrid unstructured mesh, which is heavily used in the high‐Reynolds number flow simulation, is less matured and existing methods may result in degraded mesh quality, which mostly happens in the boundary layer or near the sharp geometric features. User intervention or additional constraints, such as freezing all boundary layer elements or refining the whole boundary layer, are required to assist the refinement process. In this work, a novel AMR strategy is developed to handle existing difficulties. In the new method, high‐order unstructured elements are first generated based on the baseline mesh; then the refinement is conducted in the parametric space; at last, the mesh suitable for the solver is output. Generating refined elements in the parametric space with high‐order elements is the key of this method and this helps to guarantee both the accuracy and robustness. With the current method, 3‐dimensional hybrid unstructured mesh of huge size and complex geometry can be automatically refined, without user intervention nor additional constraints. With test cases including the 2‐dimensional airfoil and 3‐dimensional full aircraft, the current AMR method proves to be accurate, simple, and robust.