A framework for developing a mimetic tensor artificial viscosity for Lagrangian hydrocodes on arbitrary polygonal meshes

We construct a new mimetic tensor artificial viscosity on general polygonal meshes. The tensor artificial viscosity is based on discretization of coordinate invariant operators, divergence of a tensor and gradient of a vector. The focus of this paper is on the non-symmetric form, div(μ∇u), of the tensor artificial viscosity. The discretizations of this operator is derived for the case of a full tensor coefficient μ. However, in the numerical experiments, we only use scalar μ. We prove that the new tensor viscosity preserves spatial symmetry on special meshes. We demonstrate performance of the new viscosity for the Noh implosion, Sedov explosion and Saltzman piston problems on a set of various polygonal meshes in both Cartesian and axisymmetric coordinate systems.     

Shock capturing with PDE-based artificial viscosity for DGFEM: Part I. Formulation

Artificial viscosity can be combined with a higher-order discontinuous Galerkin finite element discretization to resolve a shock layer within a single cell. However, when a non-smooth artificial viscosity model is employed with an otherwise higher-order approximation, element-to-element variations induce oscillations in state gradients and pollute the downstream flow. To alleviate these difficulties, this work proposes a higher-order, state-based artificial viscosity with an associated governing partial differential equation (PDE). In the governing PDE, a shock indicator acts as a forcing term while grid-based diffusion is added to smooth the resulting artificial viscosity. When applied to heat transfer prediction on unstructured meshes in hypersonic flows, the PDE-based artificial viscosity is less susceptible to errors introduced by grid edges oblique to captured shocks and boundary layers, thereby enabling accurate heat transfer predictions.     

Bad behavior of Godunov mixed methods for strongly anisotropic advection–dispersion equations

We study the performance of Godunov mixed methods, which combine a mixed-hybrid finite element solver and a Godunov-like shock-capturing solver, for the numerical treatment of the advection–dispersion equation with strong anisotropic tensor coefficients. It turns out that a mesh locking phenomenon may cause ill-conditioning and reduce the accuracy of the numerical approximation especially on coarse meshes. This problem may be partially alleviated by substituting the mixed-hybrid finite element solver used in the discretization of the dispersive (diffusive) term with a linear Galerkin finite element solver, which does not display such a strong ill conditioning. To illustrate the different mechanisms that come into play, we investigate the spectral properties of such numerical discretizations when applied to a strongly anisotropic diffusive term on a small regular mesh. A thorough comparison of the stiffness matrix eigenvalues reveals that the accuracy loss of the Godunov mixed method is a structural feature of the mixed-hybrid method. In fact, the varied response of the two methods is due to the different way the smallest and largest eigenvalues of the dispersion (diffusion) tensor influence the diagonal and off-diagonal terms of the final stiffness matrix. One and two dimensional test cases support our findings.     

拉氏数值模拟中的人工粘性方法

较好的人工粘性需要满足较小的计算开销、不能去除真实具有的涡运动等条件.提出一种应用于拉氏数值模拟中基于Lew人工粘性,同时增加了限制器的人工粘性方法.可以有效减少数值模拟结果对网格的依赖;采用特征值限制器控制施加的人工粘性大小,通过限制器能够区分激波压缩和等熵压缩;方便应用在二维、三维,结构网格或者非结构网格上.     

三维笛卡儿坐标系中Lagrange流体力学的显式相容有限元方法(英文)

将Caramana等人提出的相容算法思想和有限元方法相结合,提出三维笛卡儿坐标系中Lagrange流体力学的显式相容有限元方法.采用三线性六面体单元和交错网格进行空间离散,利用质量集中进行显式求解,无需求解线性代数方程组.时间离散可采用两步显式Runge-Kutta格式.用边人工粘性消除激波振荡,用子网格扰动压力抑制网格的非物理变形.给出若干标准算例.数值算例表明,该方法具有较高的计算精度和计算效率,同时具有很好的对称性和总能量守恒性,总能量计算误差为计算机浮点计算截断误差.     

Behavior of viscous solutions in Lagrangian formulation

In this paper, the behavior of shock-capturing methods in Lagrangian coordinate is investigated. The relation between viscous shock and inviscid one is analyzed quantitatively, and the procedure of a viscous shock formation and propagation with a jump type initial data is described. In general, a viscous shock profile and a discontinuous one include different energy and momentum, and these discrepancies result in the generation of waves in all families when a single wave Riemann problem (shock or rarefaction) is solved. Employing this method, some anomalous behavior, such as, viscous shock interaction, shock passing through ununiform grids, postshock oscillations and lower density phenomenon is explained well. Using some classical schemes to solve the inviscid flow in Lagrangian coordinate may be not adequate enough to correctly describe flow motion in the discretized space. Partial discrepancies between von Neumann artificial viscosity method and Godunov method are exhibited. Some reviews are given to those methods which can ameliorate even eliminate entropy errors. A hybrid scheme based on the understanding to the behavior of viscous solution is proposed to suppress the overheating error.     

基于近似Riemann解的有限体积ALE方法

研究二维平面坐标系和二维轴对称坐标系中四边形网格上可压缩流体力学的有限体积ALE(Arbitrary Lagrangian Eulerian)方法.数值方法采用节点中心有限体积法,数值通量采用适用于任意状态方程的HLLC(Harten-Lax-Van Leer-Collela)通量.空间二阶精度通过用WENO(weighted essentially non-oscillatory)方法对原始变量进行重构获得,时间离散采用两步显式Runge-Kutta格式.数值例子显示,方法具有良好的激波分辨能力和高精度的数值逼近能力.     

High order well-balanced finite volume schemes for simulating wave propagation in stratified magnetic atmospheres

Wave propagation in idealized stellar atmospheres is modeled by the equations of ideal MHD, together with the gravity source term. The waves are modeled as small perturbations of isothermal steady states of the system. We consider a formulation of ideal MHD based on the Godunov–Powell form, with an embedded potential magnetic field appearing as a parameter. The equations are discretized by finite volume schemes based on approximate Riemann solvers of the HLL type and upwind discretizations of the Godunov–Powell source terms. Local hydrostatic reconstructions and suitable discretization of the gravity source term lead to a well-balanced scheme, i.e., a scheme which exactly preserves a discrete version of the relevant steady states. Higher order of accuracy is obtained by employing suitable minmod, ENO and WENO reconstructions, based on the equilibrium variables, to construct a well-balanced scheme. The resulting high order well-balanced schemes are validated on a suite of numerical experiments involving complex magnetic fields. The schemes are observed to be robust and resolve the complex physics well.     

Numerical simulation of premixed combustion using an enriched finite element method

In this paper we present a novel discretization technique for the simulation of premixed combustion based on a locally enriched finite element method (FEM). Use is made of the G-function approach to premixed combustion in which the domain is divided into two parts, one part containing the burned and another containing the unburned gases. A level-set or G-function is used to define the flame interface separating burned from unburned gases. The eXtended finite element method (X-FEM) is employed, which allows for velocity and pressure fields that are discontinuous across the flame interface. Lagrange multipliers are used to enforce the correct essential interface conditions in the form of jump conditions across the embedded flame interface. A persisting problem with the use of Lagrange multipliers in X-FEM has been the discretization of the Lagrange multipliers. In this paper the distributed Lagrange multiplier technique is adopted. We will provide results from a spatial convergence analysis showing good convergence. However, a small modification of the interface is required to ensure a unique solution. Finally, results are presented from the application of the method to the problems of moving flame fronts, the Darrieus–Landau instability and a piloted Bunsen burner flame.     

An evaluation of the FCT method for high-speed flows on structured overlapping grids

This study considers the development and assessment of a flux-corrected transport (FCT) algorithm for simulating high-speed flows on structured overlapping grids. This class of algorithm shows promise for solving some difficult highly-nonlinear problems where robustness and control of certain features, such as maintaining positive densities, is important. Complex, possibly moving, geometry is treated through the use of structured overlapping grids. Adaptive mesh refinement (AMR) is employed to ensure sharp resolution of discontinuities in an efficient manner. Improvements to the FCT algorithm are proposed for the treatment of strong rarefaction waves as well as rarefaction waves containing a sonic point. Simulation results are obtained for a set of test problems and the convergence characteristics are demonstrated and compared to a high-resolution Godunov method. The problems considered are an isolated shock, an isolated contact, a modified Sod shock tube problem, a two-shock Riemann problem, the Shu–Osher test problem, shock impingement on single cylinder, and irregular Mach reflection of a strong shock striking an inclined plane.     

Hydrodynamics of resonances

We derive system of equations describing fluidity of the medium consisting of non-relativistic particles with finite mass-widths. For that we use expressions for the kinetic Noether 4-current and the Noether energy–momentum tensor being conserved provided one uses self-consistent approximations to the gradient expanded Kadanoff–Baym equations. Kinetic coefficients entering equations of non-ideal hydrodynamics of resonances are obtained in terms of the real and imaginary parts of the self-energies within a relaxation time approximation.     

Provably unconditionally stable,second-order time-accurate,mixed variational methods for phase-field models

We introduce provably unconditionally stable mixed variational methods for phase-field models. Our formulation is based on a mixed finite element method for space discretization and a new second-order accurate time integration algorithm. The fully-discrete formulation inherits the main characteristics of conserved phase dynamics, namely, mass conservation and nonlinear stability with respect to the free energy. We illustrate the theory with the Cahn–Hilliard equation, but our method may be applied to other phase-field models. We also propose an adaptive time-stepping version of the new time integration method. We present some numerical examples that show the accuracy, stability and robustness of the new method.     

Numerical simulation of bubble dynamics in a Phan-Thien–Tanner liquid: Non-linear shape and size oscillatory response under periodic pressure

We study non-linear bubble oscillations driven by an acoustic pressure with the bubble being immersed in a viscoelastic, Phan-Thien–Tanner liquid. Solution is provided numerically through a method which is based on a finite element discretization of the Navier–Stokes flow equations. The proposed computational approach does not rely on the solution of the simplified Rayleigh–Plesset equation, is not limited in studying only spherically symmetric bubbles and provides coupled solutions for the velocity, stress fields and bubble interface. We present solutions for non-spherical bubbles, with asphericity being addressed by means of Legendre polynomials or associated Legendre functions. A parametric investigation of the bubble dynamical oscillatory response as a function of the fluid rheological properties shows that the amplitude of bubble oscillations drastically increases as liquid elasticity (quantified by the Deborah number) increases or as liquid viscosity decreases (quantified by the Reynolds number). Extensive numerical calculations demonstrate that increasing elasticity and/or viscosity of the surrounding liquid tend to stabilize the shape anisotropy of an initially non-spherical bubble. Results are shown for pressure amplitudes 0.2–2 MPa and Deborah, Reynolds numbers in the intervals of 1–8 and 0.094–1.256, respectively.     

Lovelock tensor as generalized Einstein tensor

We show that the splitting feature of the Einstein tensor, as the first term of the Lovelock tensor, into two parts, namely the Ricci tensor and the term proportional to the curvature scalar, with the trace relation between them is a common feature of any other homogeneous terms in the Lovelock tensor. Motivated by the principle of general invariance, we find that this property can be generalized, with the aid of a generalized trace operator which we define, for any inhomogeneous Euler–Lagrange expression that can be spanned linearly in terms of homogeneous tensors. Then, through an application of this generalized trace operator, we demonstrate that the Lovelock tensor analogizes the mathematical form of the Einstein tensor, hence, it represents a generalized Einstein tensor. Finally, we apply this technique to the scalar Gauss–Bonnet gravity as an another version of string–inspired gravity. This work was partially supported by a grant from the MSRT/Iran.     

非结构网格边粘性公式

利用Kuropatenko提出的人为粘性形式和Christiansen建议的TVD(Total Variation Diminishing)输运限制子构造非结构网格边粘性,用最大角原则选择邻边构造限制子,研制包括边粘性的相容流体程序,给出边粘性与体粘性数值模拟比较,对Saltzman活塞问题、Noh问题进行数值实验.数值实验结果比较说明,非结构网格边粘性准确检出强速度梯度不存在的情形,减弱人为粘性力的作用,有效捕捉流体的激波间断,提高模拟精度.     

A least-squares/finite element method for the numerical solution of the Navier–Stokes-Cahn–Hilliard system modeling the motion of the contact line

In this article we discuss the numerical solution of the Navier–Stokes-Cahn–Hilliard system modeling the motion of the contact line separating two immiscible incompressible viscous fluids near a solid wall. The method we employ combines a finite element space approximation with a time discretization by operator-splitting. To solve the Cahn–Hilliard part of the problem, we use a least-squares/conjugate gradient method. We also show that the scheme has the total energy decaying in time property under certain conditions. Our numerical experiments indicate that the method discussed here is accurate, stable and efficient.     

二维三温辐射(磁)流体力学拉氏计算中的人工粘性分裂技术

针对使用边粘性和张量粘性进行二维三温辐射(磁)流体力学拉氏计算时激波处离子温度出现非物理尖峰的问题,提出一种人工粘性分裂技术,将人工粘性耗散能量按照电子压和离子压的比分成两部分,并分别加在电子内能和离子内能上.在二维三温辐射磁流体力学程序MARED上实现该技术,并对Z-pinch的内爆动力学过程进行数值模拟.消除了模拟结果中激波处离子温度的尖峰.     

A piecewise linear finite element discretization of the diffusion equation for arbitrary polyhedral grids

We develop a piecewise linear (PWL) Galerkin finite element spatial discretization for the multi-dimensional radiation diffusion equation. It uses recently introduced piecewise linear weight and basis functions in the finite element approximation and it can be applied on arbitrary polygonal (2D) or polyhedral (3D) grids. We first demonstrate some analytical properties of the PWL method and perform a simple mode analysis to compare the PWL method with Palmer’s vertex-centered finite-volume method and with a bilinear continuous finite element method. We then show that this new PWL method gives solutions comparable to those from Palmer’s. However, since the PWL method produces a symmetric positive-definite coefficient matrix, it should be substantially more computationally efficient than Palmer’s method, which produces an asymmetric matrix. We conclude that the Galerkin PWL method is an attractive option for solving diffusion equations on unstructured grids.     

Multi-domain Fourier-continuation/WENO hybrid solver for conservation laws

We introduce a multi-domain Fourier-continuation/WENO hybrid method (FC–WENO) that enables high-order and non-oscillatory solution of systems of nonlinear conservation laws, and which enjoys essentially dispersionless, spectral character away from discontinuities, as well as mild CFL constraints (comparable to those of finite difference methods). The hybrid scheme employs the expensive, shock-capturing WENO method in small regions containing discontinuities and the efficient FC method in the rest of the computational domain, yielding a highly effective overall scheme for applications with a mix of discontinuities and complex smooth structures. The smooth and discontinuous solution regions are distinguished using the multi-resolution procedure of Harten [J. Comput. Phys. 115 (1994) 319–338]. We consider WENO schemes of formal orders five and nine and a FC method of order five. The accuracy, stability and efficiency of the new hybrid method for conservation laws is investigated for problems with both smooth and non-smooth solutions. In the latter case, we solve the Euler equations for gas dynamics for the standard test case of a Mach three shock wave interacting with an entropy wave, as well as a shock wave (with Mach 1.25, three or six) interacting with a very small entropy wave and evaluate the efficiency of the hybrid FC–WENO method as compared to a purely WENO-based approach as well as alternative hybrid based techniques. We demonstrate considerable computational advantages of the new FC-based method, suggesting a potential of an order of magnitude acceleration over alternatives when extended to fully three-dimensional problems.