Fully implicit 1D radiation hydrodynamics: Validation and verification

A fully implicit finite difference scheme has been developed to solve the hydrodynamic equations coupled with radiation transport. Solution of the time-dependent radiation transport equation is obtained using the discrete ordinates method and the energy flow into the Lagrangian meshes as a result of radiation interaction is fully accounted for. A tridiagonal matrix system is solved at each time step to determine the hydrodynamic variables implicitly. The results obtained from this fully implicit radiation hydrodynamics code in the planar geometry agrees well with the scaling law for radiation driven strong shock propagation in aluminium. For the point explosion problem the self similar solutions are compared with results for pure hydrodynamic case in spherical geometry. Results obtained when radiation interaction is also accounted agree with those of point explosion with heat conduction for lower input energies. Having, thus, benchmarked the code, self convergence of the method w.r.t. time step is studied in detail for both the planar and spherical problems. Spatial as well as temporal convergence rates are ?1 as expected from the difference forms of mass, momentum and energy conservation equations. This shows that the asymptotic convergence rate of the code is realized properly.     

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.     

各向异性网格下Wilson元的超收敛性分析

在各向异性网格下研究了二阶椭圆边值问题的Wilson有限元方法,利用单元构造的特殊性和一些新的技巧得到相应的超逼近和超收敛结果.数值算例的结果与理论分析是相吻合的.     

A unified framework of multi-objective cost functions for partitioning unstructured finite element meshes

Good performance of parallel finite element computations on unstructured meshes requires high-quality mesh partitioning. Such a decomposition task is normally done by a graph-based partitioning approach. However, the main shortcoming of graph partitioning algorithms is that minimizing the so-called edge cut is not entirely the same as minimizing the communication overhead. This paper thus proposes a unified framework of multi-objective cost functions, which take into account several factors that are not captured by the graph-based partitioning approach. Freely adjustable weighting parameters in the framework also promote a flexible treatment of different optimization objectives. A greedy-style post-improvement procedure is designed to use these cost functions to improve the quality of subdomain meshes arising from the graph-based partitioning approach. Both serial and parallel implementation of the post-improvement procedure have been done. Numerical experiments show that communication overhead can indeed be reduced by this improvement procedure, thereby increasing the performance of parallel finite element computations.     

Parameter uniform approximations for time‐dependent reaction‐diffusion problems

Using discrete Green's functions techniques, we present a classification of fitted mesh methods for time‐dependent reaction diffusion problems, based on the analyses of Linß (Linß, Numer Algor 40 (2005), 23–32) for the analogous steady‐state problem and of Kopteva (Kopteva, Computing 66 (2001), 179–197) for time‐dependent convection‐diffusion problems. As examples of how to apply the analysis, we derive error estimates for the fitted meshes of Shishkin and Bakhvalov, and provide supporting numerical results. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

一类非正则网格上的非协调Mortar元的高精度分析

In this paper, the nonconforming mortar finite element with a class of meshes is studied without considering the global regularity condition or quasi-uniformly assumption. Meanwhile, the superclose result coincides with conventional methods is obtained by means of integral identities techniques.     

Mixed hp-DGFEM for incompressible flows II: Geometric edge meshes

We consider the Stokes problem of incompressible fluid flowin three-dimensional polyhedral domains discretized on hexahedralmeshes with hp-discontinuous Galerkin finite elements of typeQ_k for the velocity and Q_k–1 for the pressure. We provethat these elements are inf-sup stable on geometric edge meshesthat are refined anisotropically and non-quasiuniformly towardsedges and corners. The discrete inf-sup constant is shown tobe independent of the aspect ratio of the anisotropic elementsand is of O(k^–3/2) in the polynomial degree k, as in thecase of conforming Q_k–Q_k–2 approximations on thesame meshes.     

Exponential meshes and three‐dimensional computation of a magnetic field

We describe the simulation of an exterior problem using a magnetic field deriving from magnetostatics, with a numerical method mixing the approaches of C. I. Goldstein and L.‐A. Ying. This method is based on the use of a graded mesh obtained by gluing homothetic layers in the exterior domain. On this mesh, we use an edge elements discretization and a recently proposed mixed formulation. In this paper, we provide both a theoretical and numerical study of the method. We establish an error estimate, describe the implementation, propose some preconditioning techniques and show the numerical results. We also compare these results with those obtained from an equivalent boundary elements approach. In this way, we retain that our method leads to a practical numerical implementation, a saving of storage, and turns out to be an alternative to the classical boundary elements method. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 595–637, 2003     

多边形网格上扩散方程新的单调格式

 本文在星形多边形网格上, 构造了扩散方程新的单调有限体积格式.该格式与现有的基于非线性两点流的单调格式的主要区别是, 在网格边的法向流离散模板中包含当前边上的点, 在推导离散法向流的表达式时采用了定义于当前边上的辅助未知量, 这样既可适应网格几何大变形, 同时又兼顾了当前网格边上物理量的变化. 在光滑解情形证明了离散法向流的相容性.对于具有强各向异性、非均匀张量扩散系数的扩散方程, 证明了新格式是单调的, 即格式可以保持解析解的正性. 数值结果表明在扭曲网格上, 所构造的格式是局部守恒和保正的, 对光滑解有高于一阶的精度, 并且, 针对非平衡辐射限流扩散问题, 数值结果验证了新格式在计算效率和守恒精度上优于九点格式.     

A simple weighted essentially non-oscillatory limiter for the correction procedure via reconstruction (CPR) framework on unstructured meshes

In this paper, we adapt a simple weighted essentially non-oscillatory (WENO) limiter, originally designed for discontinuous Galerkin (DG) schemes on two-dimensional unstructured triangular meshes [39], to the correction procedure via reconstruction (CPR) framework for solving nonlinear hyperbolic conservation laws on two-dimensional unstructured triangular meshes with straight or curved edges. This is an extension of our earlier work [4] in which the WENO limiter was designed for the CPR framework on regular meshes. The objective of this simple WENO limiter is to simultaneously maintain uniform high order accuracy of the CPR framework in smooth regions and control spurious numerical oscillations near discontinuities. The WENO limiter we adopt in this paper uses information only from the target cell and its immediate neighbors. Hence, it is particularly simple to implement and will not harm the conservativeness and compactness of the CPR framework. Since the CPR framework with this WENO limiter does not in general satisfy the positivity preserving property, we also extend the positivity-preserving limiters [36], [33] to the CPR framework. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good behavior of this procedure.