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.     

Relaxation of the Navier–Stokes–Korteweg equations for compressible two‐phase flow with phase transition

The Navier–Stokes–Korteweg (NSK) system is a classical diffuse‐interface model for compressible two‐phase flow. However, the direct numerical simulation based on the NSK system is quite expensive and in some cases even not possible. We propose a lower‐order relaxation of the NSK system with hyperbolic first‐order part. This allows applying numerical methods for hyperbolic conservation laws and removing some of the difficulties of the original NSK system. To illustrate the new ansatz, we first present a local discontinuous Galerkin method in one and two spatial dimensions. It is shown that we can compute initial boundary value problems with realistic density ratios and perform stable computations for small interfacial widths. Second, we show that it is possible to construct a semi‐discrete finite‐volume scheme that satisfies a discrete entropy inequality. Copyright © 2015 John Wiley & Sons, Ltd.     

塑料SHPB实验中的端面摩擦效应

基于H.Meng的动摩擦实验结果,建立了一个端面动摩擦模型,以描述端面动摩擦因数随SHPB 塑料试样与钢压杆接触界面间最大径向相对滑动速度的关系。在SHPB实验的数值模拟中,分别考虑了该 端面动摩擦模型和以往被广泛采用的常摩擦因数模型,结果发现,选取不同的端面摩擦模型不影响通过塑料 SHPB实验的数值模拟获得的转折应变率(量级为102s-1),然而当应变率超出转折应变率以上时,不同的端 面摩擦模型会对塑料SHPB实验的数值模拟结果产生显著影响,使用提出的端面动摩擦模型可得到更准确 的塑料SHPB实验的数值模拟结果。     

Application of Hénon method in numerical estimation of the stick–slip transitions existing in Filippov-type discontinuous dynamical systems with dry friction

Accuracy of numerical modeling of any discontinuous dynamical system plays an important role in the proper use of the tools applied for its analysis. No less important in this matter is the numerical estimation of the phase trajectories, bifurcation diagrams, and Lyapunov exponents. This paper meets these expectations presenting application of Hénon’s method to obtain good numerical estimations of the stick–slip transitions existing in the Filippov-type discontinuous dynamical systems with dry friction. Subsequent sections are focused on the problem definition, block diagrams of the numerical procedure coded in Python, application of the method that was originally proposed for Poincaré maps, its use in estimation of phase trajectories, bifurcation diagrams of tangent points, and on estimation of the Lyapunov exponents for a selected two-dimensional system with stick–slip effect.     

discontinuous galerkin method for discontinuous temperature field problems

针对不连续温度场问题建立了一种间断Galerkin有限元方法,该方法的主要特点是允 许插值函数在单元边界上存在跳变. 在建立有限元方程时,通过在单元边界上引入数值通量 项和稳定性项来处理间断效应,并且数值通量可以直接由接触热阻的定义式导出. 数值算例 表明该方法可以很方便且准确地捕捉到结构内部由于接触热阻而引起的温度跳变,同时在局 部高梯度温度场的模拟方面也比常规连续Galerkin有限元方法效率明显要高. 该方法也为研 究由接触热阻引起的温度场与应力场之间的耦合问题提供了一种新的数值模拟手段.     

Interfacial stress analysis and prediction of debonding for a thin plate bonded to a curved substrate

This paper focuses on the analytical and numerical modeling of the interface between a rigid substrate with simple constant curvature and a thin bonded plate. The interfacial behavior is modeled by independent cohesive laws in the normal and tangential directions, coupled with a mixed-mode fracture criterion. The newly developed analytical model determines the interfacial shear and normal stress distributions as functions of the substrate curvature, during the various behavioral stages of the interface prior to the initiation of debonding. The model is also able to predict the debonding load and the effective bond length. In the numerical model the interface is modeled by zero-thickness node-to-segment contact elements, in which both the geometrical relationships between the nodes of the discretized problem and the interface constitutive laws are suitably defined. Numerical results and comparisons between the predictions of the two models are presented.     

基于高斯过程分类的结构贝叶斯可靠性分析

贝叶斯可靠性方法是处理不完备信息条件下结构可靠性问题的有效途径之一。在实际应用中,由于可靠性分析的计算量较大,常须采用各种近似替代模型以提高计算效率。传统的替代模型方法是对结构的功能函数予以近似建模。这种方法不易定量考虑模型误差对可靠性分析的影响,且难以应用于诸如功能函数不连续和失效域不连通等情况。为此,本文提出一种基于高斯过程分类的替代模型,直接辨识结构的极限状态曲面,并将其应用于结构贝叶斯可靠性分析之中。分析了替代模型不确定性对可靠性预测结果的影响,给出了失效概率分布参数的方差算式,进而提出了改善模型精度的补充采样准则。通过算例验证了方法的适用性和有被性．     

不连续温度场问题的间断Galerkin方法

针对不连续温度场问题建立了一种间断Galerkin有限元方法,该方法的主要特点是允许插值函数在单元边界上存在跳变.在建立有限元方程时,通过在单元边界上引入数值通量项和稳定性项来处理间断效应,并且数值通量可以直接由接触热阻的定义式导出.数值算例表明该方法可以很方便且准确地捕捉到结构内部由于接触热阻而引起的温度跳变,同时在局部高梯度温度场的模拟方面也比常规连续Galerkin有限元方法效率明显要高.该方法也为研究由接触热阻引起的温度场与应力场之间的耦合问题提供了一种新的数值模拟手段.     

A design of residual error estimates for a high order BDF‐DGFE method applied to compressible flows

We deal with the numerical solution of the non‐stationary compressible Navier–Stokes equations with the aid of the backward difference formula – discontinuous Galerkin finite element method. This scheme is sufficiently stable, efficient and accurate with respect to the space as well as time coordinates. The nonlinear algebraic systems arising from the backward difference formula – discontinuous Galerkin finite element discretization are solved by an iterative Newton‐like method. The main benefit of this paper are residual error estimates that are able to identify the computational errors following from the space and time discretizations and from the inexact solution of the nonlinear algebraic systems. Thus, we propose an efficient algorithm where the algebraic, spatial and temporal errors are balanced. The computational performance of the proposed method is demonstrated by a list of numerical experiments. Copyright © 2013 John Wiley & Sons, Ltd.     

Nonlinear parameter identification of models for masonry

This paper investigates the parameter estimation problem for brick masonry models. An identification procedure is proposed in which the uncertainties of known parameters and/or errors of measurements are its elements of distinction. The minimization process of the discrepancies between experimental data and theoretical measurements takes place by means of a first order iterative method. The identification procedure is applied to two different problems: the calibration of an interface model for brick–mortar joint in its functional form through monotonic experimental tests; to evaluate the unknown parameters of a continuum model for brick masonry walls in its non-holonomic form by means of in-plane cyclic shear–compression test of masonry panels. The general framework of the non-linear estimate methodology, the parameter identification problems and the numerical results are presented.     

两介质流界面-激波相互作用RKDG 方法应用分析

为精确模拟多介质流界面运动现象,采用RKDG方法结合虚拟流体方法对气-气、气-液和液-气等多种界面-激波相互作用问题展开研究。数值结果表明,RKDG方法的时空高精度特征使其能够精确、稳健地求解各种复杂界面运动问题。最后,对水下激波自由面折射问题用多种DG格式限制器进行了计算,对比了它们的间断捕捉能力。     

A discontinuous Galerkin‐based sharp‐interface method to simulate three‐dimensional compressible two‐phase flow

A numerical method for the simulation of compressible two‐phase flows is presented in this paper. The sharp‐interface approach consists of several components: a discontinuous Galerkin solver for compressible fluid flow, a level‐set tracking algorithm to follow the movement of the interface and a coupling of both by a ghost‐fluid approach with use of a local Riemann solver at the interface. There are several novel techniques used: the discontinuous Galerkin scheme allows locally a subcell resolution to enhance the interface resolution and an interior finite volume Total Variation Diminishing (TVD) approximation at the interface. The level‐set equation is solved by the same discontinuous Galerkin scheme. To obtain a very good approximation of the interface curvature, the accuracy of the level‐set field is improved and smoothed by an additional P_NP_M‐reconstruction. The capabilities of the method for the simulation of compressible two‐phase flow are demonstrated for a droplet at equilibrium, an oscillating ellipsoidal droplet, and a shock‐droplet interaction problem at Mach 3. Copyright © 2015 John Wiley & Sons, Ltd.     

Wall modeling via function enrichment within a high‐order DG method for RANS simulations of incompressible flow

We present a novel approach to wall modeling for the Reynolds‐averaged Navier‐Stokes equations within the discontinuous Galerkin method. Wall functions are not used to prescribe boundary conditions as usual, but they are built into the function space of the numerical method as a local enrichment, in addition to the standard polynomial component. The Galerkin method then automatically finds the optimal solution among all shape functions available. This idea is fully consistent and gives the wall model vast flexibility in separated boundary layers or high adverse pressure gradients. The wall model is implemented in a high‐order discontinuous Galerkin solver for incompressible flow complemented by the Spalart‐Allmaras closure model. As benchmark examples, we present turbulent channel flow starting from Re_τ=180 and up to Re_τ=100000 as well as flow past periodic hills at Reynolds numbers based on the hill height of Re_H=10595 and Re_H=19000.     

Identification of shallow sea models

In this paper we consider a parameter estimation procedure for shallow sea models. The method is formulated as a minimization problem. An adjoint model is used to calculate the gradient of the criterion which is to be minimized. In order to obtain a robust estimation method, the uncertainty of the open boundary conditions can be taken into acoount by allowing random noise inputs to act on the open boundaries. This method avoids the possibility that boundary errors are interpreted by the estimation procedure as parameter fluctuations. We apply the parameter estimation method to identify a shallow sea model of the entire European continental shelf. First, a space-varying bottom friction coefficient is estimated simultaneously with the depth. The second application is the estimation of the parameterization of the wind stress coefficient as a function of the wind velocity. Finally, an uncertain open boundary condition is included. It is shown that in this case the parameter estimation procedure does become more robust and produces more realistic estimates. Furthermore, an estimate of the open boundary conditions is also obtained.     

热冲击作用下层合皮肤组织热传导过程数值模拟

在急剧温度变化等强间断温度冲击作用下的生物层合组织非傅里叶热传导分析中,经典时域连续有限元方法(如Newmark等方法)会在波阵面以后的和层合组织界面附近的区域表现出强烈的数值振荡。这类数值振荡会影响问题求解精度,并带来较大不确定性。针对这类现象,本文发展了改进时域间断Galerkin有限元方法,进一步开展了相关问题的数值模拟。其控制方程的基本未知数(温度)及其时间导数在指定时间间隔内假设存在间断且独立插值。在有限元离散列式中引入比例刚度阵人工阻尼,以成功消除波前位置的虚假数值振荡行为。通过算例对比分析,相比Newmark方法和传统间断Galerkin方法,所提出的改进时域间断Galerkin有限元方法较好消除了波前、波后以及组织界面处的数值振荡,有效捕捉了波阵面的间断行为,提高了计算的精度。     

Local discontinuous galerkin method for radial porous flow with dispersion and adsorption

Based on the local discontinuous Galerkin methods for time-dependent convection-diffusion systems newly developed by Corkburn and Shu, according to the form of the generalized convection-diffusion equations which model the radial porous flow with dispersion and adsorption, a local discontinuous Galerkin method for radial porous flow with dispersion and adsorption was developed, a high order accurary new scheme for radial porous flow is obtained. The presented method was applied to the numerical tests of two cases of radial porous, i. e. , the convection-dispersion flow and the convection-dispersion-adsorption flow, the corresponding parts of the numerical results are in good agreement with the published solutions, so the presented method is reliable. Reckoning of the computational cost also shows that the method is practicable.     

An improved three-dimensional model for interface pressure calculations in free-surface flows

A three-dimensional method for the calculation of interface pressure in the computational modeling of free surfaces and interfaces is developed. The methodology is based on the calculation of the pressure force at the interfacial cell faces and is mainly designed for volume of fluid (VOF) interface capturing approach. The pressure forces at the interfacial cell faces are calculated according to the pressure imposed by each fluid on the portion of the cell face that is occupied by that fluid. Special formulations for the pressure in the interfacial cells are derived for different orientations of an interface. The present method, referred to as pressure calculation based on the interface location (PCIL), is applied to both static and dynamic cases. First, a three-dimensional motionless drop of liquid in an initially stagnant fluid with no gravity force is simulated as the static case and then two different small air bubbles in water are simulated as dynamic cases. A two-fluid, piecewise linear interface calculation VOF method is used for numerical simulation of the interfacial flow. For the static case, both the continuum surface force (CSF) and the continuum surface stress (CSS) methods are used for surface tension calculations. A wide range of Ohnesorge numbers and density and viscosity ratios of the two fluids are tested. It is shown that the presence of spurious currents (artificial velocities present in case of considerable capillary forces) is mainly due to the inaccurate calculation of pressure forces in the interfacial computational cells. The PCIL model reduces the spurious currents up to more than two orders of magnitude for the cases tested.

Also for the dynamic bubble rise case, it is shown that using the numerical solver employed here, without PCIL, the magnitude of spurious currents is so high that it is not possible to simulate this type of surface tension dominated flows, while using PCIL, we are able to simulate bubble rise and obtain results in close agreement with the experimental data.     

含摩擦与碰撞平面多刚体系统动力学线性互补算法

基于接触力学理论和线性互补问题的算法, 给出了一种含接触、碰撞以及库伦干摩擦, 同时具有理想定常约束(铰链约束) 和非定常约束(驱动约束) 的平面多刚体系统动力学的建模与数值计算方法. 将系统中的每个物体视为刚体, 但考虑物体接触点的局部变形, 将物体间的法向接触力表示成嵌入量与嵌入速度的非线性函数,其切向摩擦力采用库伦干摩擦模型. 利用摩擦余量和接触点的切向加速度等概念, 给出了摩擦定律的互补关系式; 并利用事件驱动法, 将接触点的黏滞-滑移状态切换的判断及黏滞状态下摩擦力的计算问题转化成线性互补问题的求解. 利用第一类拉格朗日方程和鲍姆加藤约束稳定化方法建立了系统的动力学方程, 由此可降低约束的漂移, 并可求解该系统的运动、法向接触力和切向摩擦力, 还可以求解理想铰链约束力和驱动约束力. 最后以一个类似夯机的平面多刚体系统为例, 分析了其动力学特性, 并说明了相关算法的有效性.      

On the Influence of Polynomial De-aliasing on Subgrid Scale Models

In this work we investigate the interplay of polynomial de-aliasing and sub-grid scale models for large eddy simulations based on discontinuous Galerkin discretizations. It is known that stability is a major concern when simulating underresolved turbulent flows with high order nodal collocation type discretizations. By changing the interpolatory character of the nodal collocation type discretization to a projection based discretization by increasing the number of quadrature points (polynomial de-aliasing), one is able to remove the aliasing induced stability problems. We focus on this effect and on the consequence for large eddy simulations with explicit subgrid scale models. Often, subgrid scale models have to achieve two possibly conflicting tasks in a single simulation: firstly stabilizing the numerics and secondly modeling the physical effect of the missing scales. Within a discontinuous Galerkin approach, it is possible to use either a fast (but potentially aliasing afflicted) nodal collocation discretization or a projection-based (but computationally costly) variant in combination with an explicit subgrid scale model. We use this framework to investigate the effect on the appropriate model parameter of a standard Smagorinsky subgrid scale model and of a Variational Multiscale Smagorinsky formulation. For this we first consider the 3-D viscous Taylor-Green vortex example to investigate the impact on the stability of the method and second the turbulent flow past a circular cylinder to investigate and compare the accuracy of the results. We show that the aliasing instabilities of collocative discretizations severely limit the choice of the model constant, in particular for high order schemes, while for de-aliased DG schemes, the closure model parameters can be chosen independently from the numerical scheme. For the cylinder flow, we also find that for the same model settings, the projection-based results are in better agreement with the reference DNS than those of the collocative scheme.