A compact high‐order HLL solver for nonlinear hyperbolic systems

We present a new finite‐volume method for calculating complex flows on non‐uniform meshes. This method is designed to be highly compact and to accurately capture all discontinuities that may arise within the solution of a nonlinear hyperbolic system. In the first step, we devise a fourth‐degree Hermite polynomial to interpolate the solution. The coefficients defining this polynomial are calculated by using a least‐square method. To introduce monotonicity conditions within the procedure, two constraints are added into the least‐square system. Those constraints are derived by locally matching the high‐order Hermite polynomial with a low‐order TVD polynomial. To emulate these constraints only in regions of discontinuities, data‐depending weights are defined; these weights are based upon normalized indicators of smoothness of the solution and are parameterized by an O(1) quantity. The reconstruction so generated is highly compact and is fifth‐order accurate when the solution is smooth; this reconstruction becomes first order in regions of discontinuities. In the second step, this reconstruction is inserted in an HLL approximate Riemann solver. This solver is designed to correctly capture all discontinuities that may arise into the solution. To this aim, we introduce the contribution of a possible contact discontinuity into the HLL Riemann solver. Thus, a spatially fifth‐order non‐oscillatory method is generated. This method evolves in time the solution and its first derivative. In a one‐dimensional context, a linear spectral analysis and extensive numerical experiments make it possible to assess the robustness and the advantages of the method in computing multi‐scale problems with embedded discontinuities. Copyright © 2008 John Wiley & Sons, Ltd.     

Stokes问题各向异性网格Q2-P1混合元超收敛分析

讨论Stokes问题在各向异性冈格下的Q2-P1混合有限元方法，利用积分恒等式技巧得到了与传统方法相同的超逼近性质，同时基于插值后处理的技巧，构造了速度和压力的一对插值后处理算子，并且前者具有备向异性特征，从而导出了整体超收敛结果．     

A modified Bakhvalov mesh

We consider a few numerical methods for solving a one-dimensional convection–diffusion singularly perturbed problem. More precisely, we introduce a modified Bakvalov mesh generated using some implicitly defined functions. Properties of this mesh and convergence results for several methods on it are given. Numerical results are presented in support of the theoretical considerations.     

Generating Layer-Adapted Meshes Using Mesh Partial Differential Equations

We present a new algorithm for generating layer-adapted meshes for the finite element solution of singularly perturbed problems based on mesh partial differential equations (MPDEs). The ultimate goal is to design meshes that are similar to the well-known Bakhvalov meshes, but can be used in more general settings: specifically two-dimensional problems for which the optimal mesh is not tensor-product in nature. Our focus is on the efficient implementation of these algorithms, and numerical verification of their properties in a variety of settings. The MPDE is a nonlinear problem, and the efficiency with which it can be solved depends adversely on the magnitude of the perturbation parameter and the number of mesh intervals. We resolve this by proposing a scheme based on $h$-refinement. We present fully working FEniCS codes [Alnaes et al., Arch. Numer. Softw., 3 (100) (2015)] that implement these methods, facilitating their extension to other problems and settings.     

二阶椭圆问题的各向异性非协调Crouzeix-Raviart型有限元方法

对满足最大角条件和坐标系条件的二维区域中的各向异性一般三角形网格,研究了二阶椭圆问题的非协调Crouzeix-Raviart型线性三角形有限元逼近,得到了最优的能量模和L2-模误差估计结果．     

Integrating Hydrogeochemical and Geophysical Data for Testing a Finite Volume Based Numerical Model for Saltwater Intrusion

In this paper, hydrogeological and geophysical data are used to validate a numerical model developed to predict seawater intrusion into coastal aquifers. The cell-centered finite volume method is adopted here to solve the set of coupled partial differential equations describing the motion of saltwater and freshwater separated by a sharp interface. These equations are based on the Dupuit approximation and are obtained from integration of 3D flow equations for fresh and salt water zones over the vertical dimension. In order to have flexibility upon complex configurations domain, non structured grid meshing is utilized. To approximate the diffusion fluxes, Green-Gauss type reconstruction, based on diamond-cell and least squares interpolation, is performed. The model is first validated using academic test case studies with known closed form solutions. The mathematical model has been calibrated using hydrogeochemical and geophysical data. The geophysical method applied in this study has been a frequency domain electromagnetic method. In this method the apparent electrical conductivity is measured by induction using two separate hand-held transmitter and receiver coils. During the operation the transmitter coil is energized by a low frequency alternating current that radiates an electromagnetic field and the receiver coil detects the resulting field. Taking into account the relationship between the bulk conductivity of the subsoil and the conductivity of groundwater, EM soundings have been interpreted to provide complementary information to hydrogeochemical data to outline the fresh–saltwater interface. This methodology has been applied to the case of saltwater intrusion into the Llobregat delta aquifer, near Barcelona, Spain.     

A posteriori error estimates for optimal control problems constrained by convection-diffusion equations

We propose a characteristic finite element discretization of evolutionary type convection-diffusion optimal control problems. Nondivergence-free velocity fields and bilateral inequality control constraints are handled. Then some residual type a posteriori error estimates are analyzed for the approximations of the control, the state, and the adjoint state. Based on the derived error estimators, we use them as error indicators in developing efficient multi-set adaptive meshes characteristic finite element algorithm for such optimal control problems. Finally, one numerical example is given to check the feasibility and validity of multi-set adaptive meshes refinements.     

A finite volume scheme preserving extremum principle for convection–diffusion equations on polygonal meshes

We propose a nonlinear finite volume scheme for convection–diffusion equation on polygonal meshes and prove that the discrete solution of the scheme satisfies the discrete extremum principle. The approximation of diffusive flux is based on an adaptive approach of choosing stencil in the construction of discrete normal flux, and the approximation of convection flux is based on the second‐order upwind method with proper slope limiter. Our scheme is locally conservative and has only cell‐centered unknowns. Numerical results show that our scheme can preserve discrete extremum principle and has almost second‐order accuracy. Copyright © 2017 John Wiley & Sons, Ltd.     

有限元计算中疏密网格过渡方法研究

工程计算中出于节省计算量的目的,往往需要在一个有限元模型中布置粗细不同的网格。为保证计算结果的准确性,必须保证网格突变情况下的位移协调问题。本文工作之一是在强天驰界面过渡单元的基础上,引入虚拟节点和子单元,在子单元中应用节理元思想,提出了基于最小势能原理的弹簧节理单元法。简化了积分运算,避免了精度要求极高的坐标转换,从而提高了方法的精度和实用性;二是提出了基于位移约束的主从自由度法,简便实用,只需简单的矩阵运算即可实现。两种方法均实现了不同尺寸网格间位移的协调性和刚度的匹配,从而使之满足有限元收敛准则,且生成的刚度阵具有对称性及带状性。算例证明两种方法精度良好,并可方便地应用于求解大规模工程问题。     

Adaptive unstructured finite element method for two-dimensional detonation simulations

A non-equilibrium reacting flow methodology has been added to a conservative, monotonic, compressible flow solver to allow numerical simulations of gas detonations. This flow solver incorporates unstructured dynamically adaptive meshes with the Finite Element Method – Flux Corrected Transport (FEM-FCT) scheme, which has shown excellent predictive capability of various non-reacting compressible flows. A two-step induction parameter model was used to model the combustion of the gas phase coupled with an energy release equation which was simulated with a point implicit finite element scheme. This combustion model was then applied to a two-dimensional detonation test case of a hypothetical fuel:oxygen mixture. The detonation simulation yielded two transverse waves which continued to propagate. This feature and the detonation shock speed mean and fluctuations were found to be grid-independent based on a resolution of about twenty elements within the average induction length. The resolution was efficiently achieved with the unstructured dynamically adaptive finite elements, which were three orders of magnitude less in number then required for uniform discretization. Received 26 August 1996 / Accepted 31 March 1997