A cell-centered Lagrangian hydrodynamics scheme on general unstructured meshes in arbitrary dimension

We describe a cell-centered Godunov scheme for Lagrangian gas dynamics on general unstructured meshes in arbitrary dimension. The construction of the scheme is based upon the definition of some geometric vectors which are defined on a moving mesh. The finite volume solver is node based and compatible with the mesh displacement. We also discuss boundary conditions. Numerical results on basic 3D tests problems show the efficiency of this approach. We also consider a quasi-incompressible test problem for which our nodal solver gives very good results if compared with other Godunov solvers. We briefly discuss the compatibility with ALE and/or AMR techniques at the end of this work. We detail the coefficients of the isoparametric element in the appendix.     

The influence of cell geometry on the Godunov scheme applied to the linear wave equation

By studying the structure of the discrete kernel of the linear acoustic operator discretized with a Godunov scheme, we clearly explain why the behaviour of the Godunov scheme applied to the linear wave equation deeply depends on the space dimension and, especially, on the type of mesh. This approach allows us to explain why, in the periodic case, the Godunov scheme applied to the resolution of the compressible Euler or Navier–Stokes system is accurate at low Mach number when the mesh is triangular or tetrahedral and is not accurate when the mesh is a 2D (or 3D) cartesian mesh. This approach confirms also the fact that a Godunov scheme remains accurate when it is modified by simply centering the discretization of the pressure gradient.     

Lagrange方法基于保正性的时间步长

对交错网格上Lagrange预估-校正显示格式的时间步长选取提出新的方法.与经典CFL稳定性理论时间步长选取方法不同,新方法考虑了原始微分方程组的非线性效应,并基于物理量保正性给出了自适应的时间步长选取方法.数值实验验证了该方法有效.     

类锂铝复合机制产生X光激光的等离子体状态

本文从速率方程出发,讨论了类锂铝复合等离子体的激发态结构,衰减常数,反转率和小信号增益等表征介质增益特性的物理量以及它们随电子温度,电子密度和光子逃逸几率的变化。找到了进行类锂铝离子通过复合机制产生X光激光设计应创造的等离子体状态目标区域。还讨论了这些物理量随原子序数变化的定标律。     

Comparison of the generalized Riemann solver and the gas-kinetic scheme for inviscid compressible flow simulations

The generalized Riemann problem (GRP) scheme for the Euler equations and gas-kinetic scheme (GKS) for the Boltzmann equation are two high resolution shock capturing schemes for fluid simulations. The difference is that one is based on the characteristics of the inviscid Euler equations and their wave interactions, and the other is based on the particle transport and collisions. The similarity between them is that both methods can use identical MUSCL-type initial reconstructions around a cell interface, and the spatial slopes on both sides of a cell interface involve in the gas evolution process and the construction of a time-dependent flux function. Although both methods have been applied successfully to the inviscid compressible flow computations, their performances have never been compared. Since both methods use the same initial reconstruction, any difference is solely coming from different underlying mechanism in their flux evaluation. Therefore, such a comparison is important to help us to understand the correspondence between physical modeling and numerical performances. Since GRP is so faithfully solving the inviscid Euler equations, the comparison can be also used to show the validity of solving the Euler equations itself. The numerical comparison shows that the GRP exhibits a slightly better computational efficiency, and has comparable accuracy with GKS for the Euler solutions in 1D case, but the GKS is more robust than GRP. For the 2D high Mach number flow simulations, the GKS is absent from the shock instability and converges to the steady state solutions faster than the GRP. The GRP has carbuncle phenomena, likes a cloud hanging over exact Riemann solvers. The GRP and GKS use different physical processes to describe the flow motion starting from a discontinuity. One is based on the assumption of equilibrium state with infinite number of particle collisions, and the other starts from the non-equilibrium free transport process to evolve into an equilibrium one through particle collisions. The different mechanism in the flux evaluation deviates their numerical performance. Through this study, we may conclude scientifically that it may NOT be valid to use the Euler equations as governing equations to construct numerical fluxes in a discretized space with limited cell resolution. To adapt the Navier–Stokes (NS) equations is NOT valid either because the NS equations describe the flow behavior on the hydrodynamic scale and have no any corresponding physics starting from a discontinuity. This fact alludes to the consistency of the Euler and Navier–Stokes equations with the continuum assumption and the necessity of a direct modeling of the physical process in the discretized space in the construction of numerical scheme when modeling very high Mach number flows. The development of numerical algorithm is similar to the modeling process in deriving the governing equations, but the control volume here cannot be shrunk to zero.     

Alternating direction implicit schemes for the two-dimensional fractional sub-diffusion equation

New numerical techniques are presented for the solution of a two-dimensional anomalous sub-diffusion equation with time fractional derivative. In these methods, standard central difference approximation is used for the spatial discretization, and, for the time stepping, two new alternating direction implicit (ADI) schemes based on the L₁ approximation and backward Euler method are considered. The two ADI schemes are constructed by adding two different small terms, which are different from standard ADI methods. The solvability, unconditional stability and H¹ norm convergence are proved. Numerical results are presented to support our theoretical analysis and indicate the efficiency of both methods.     

Unconditionally Stable Z-Transform ADI-PML Algorithm for Lorentzian Double Negative Meta-Materials

Unconditionally stable formulations of the anisotropic perfectly matched layer (APML) are presented for truncating double negative (DNG) meta-material finite difference time domain (FDTD) grids. In the proposed formulations, the Z-transform theory is employed in the alternating direction implicit FDTD (ADI-FDTD) scheme to obtain update equations for the field components in the DNG meta-material domains. Numerical examples carried out in one dimensional Lorentzian type DNG meta-material domains are included to show the validity of the proposed formulations.     

轴对称射流气动声场的数值模拟

构造了Kirchhoff积分和CFD计算相结合的算法,采用高精度差分格式对不同喷口马赫数的轴对称射流进行数值模拟,并以此作为近场声源,运用Kirchhoff积分方法研究远场气动噪声.数值结果表明,远场噪声具有方向性,噪声声压在离开对称轴20°处达到最大值.随着传播距离增大,噪声方向性逐渐减弱.     

分辨率优化的混合WENO格式

为提高有限差分格式的分辨率,利用傅里叶分析对WENO格式进行色散及耗散优化,并给出优化的线性权重.用优化后的WENO格式与保单调格式（MP）进行加权混合,得到新的加权混合WENO格式（H-WENO）.通过一维激波管问题、Shu-Osher问题及二维双Mach反射问题及R-T不稳定性问题对格式进行数值测试.结果显示,新格式具有强健的激波捕捉能力和对小尺度波结构的高分辨率,与原WENO格式相比改进明显.     

A Finite Volume Method Preserving Maximum Principle for the Conjugate Heat Transfer Problems with General Interface Conditions

In this paper, we present a unified finite volume method preserving discrete maximum principle (DMP) for the conjugate heat transfer problems with general interface conditions. We prove the existence of the numerical solution and the DMP-preserving property. Numerical experiments show that the nonlinear iteration numbers of the scheme in [24] increase rapidly when the interfacial coefficients decrease to zero. In contrast, the nonlinear iteration numbers of the unified scheme do not increase when the interfacial coefficients decrease to zero, which reveals that the unified scheme is more robust than the scheme in [24]. The accuracy and DMP-preserving property of the scheme are also veried in the numerical experiments.