基于二重网格的动网格松弛法改进

相比传统的弹簧法等方法,基于球松弛算法的动网格松弛法在复杂边界大变形条件下可以得到质量更高的边界网格以及更大的极限变形量,但该方法在时间效率上还有提升的空间.引入二重网格,采用动网格松弛法进行稀疏网格的网格变形,将边界位移传递到整个网格计算域;再利用二重网格映射,将稀疏网格位移映射到原有计算网格的节点上.算例表明,改进...     

Upwind-biased FORCE schemes with applications to free-surface shallow flows

In this paper we develop numerical fluxes of the centred type for one-step schemes in conservative form for solving general systems of conservation laws in multiple-space dimensions on structured meshes. The proposed method is an extension of the multidimensional FORCE flux developed by Toro et al. (2009) [14]. Here we introduce upwind bias by modifying the shape of the staggered mesh of the original FORCE method. The upwind bias is evaluated using an estimate of the largest eigenvalue, which in any case is needed for selecting a time step. The resulting basic flux is first-order accurate and monotone. For the linear advection equation, the proposed UFORCE method reproduces exactly the upwind Godunov method. Extension to non-linear systems has been done empirically via the two-dimensional inviscid shallow water equations. Second order of accuracy in space and time on structured meshes is obtained in the framework of finite volume methods. The proposed method improves the accuracy of the solution for small Courant numbers and intermediate waves associated with linearly degenerate fields (contact discontinuities, shear waves and material interfaces). It achieves comparable accuracy to that of upwind methods with approximate Riemann solvers, though retaining the simplicity and efficiency of centred methods. The performance of the schemes is assessed on a suite of test problems for the two-dimensional shallow water equations.     

UNIFORM SUPERCONVERGENCE OF GALERKIN METHODS FOR SINGULARLY PERTURBED PROBLEMS

In this paper, we are concerned with uniform superconvergence of Galerkin methods for singularly perturbed reaction-diffusion problems by using two Shishkin-type meshes. Based on an estimate of the error between spline interpolation of the exact solution and its numerical approximation, an interpolation post-processing technique is applied to the original numerical solution. This results in approximation exhibit superconvergence which is uniform in the weighted energy norm. Numerical examples are presented to demonstrate the effectiveness of the interpolation post-processing technique and to verify the theoretical results obtained in this paper.     

Small perturbations of polynomial meshes

We show that the property of being a (weakly) admissible mesh for multivariate polynomials is preserved by small perturbations on real and complex Markov compacts. Applications are given to smooth transformations of polynomial meshes and to polynomial interpolation.     

各向异性网格下抛物方程一个新的非协调混合元收敛性分析

本文将 Crouzeix-Raviart 型非协调线性三角形元应用到抛物方程,建立了一个新的混合元格式.在抛弃传统有限元分析的必要工具 Ritz 投影算子的前提下,直接利用单元的插值性质和导数转移技巧, 分别得到了各向异性剖分下关于原始变量u 的H^-1-模和积分意义下L²-模以及通量p=-▽u 在L²-模下的最优阶误差估计.数值结果与我们的理论分析是相吻合的.     

A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

Virtual element methods for parabolic problems on polygonal meshes

The virtual element method (VEM) is a recent technology that can make use of very general polygonal/polyhedral meshes without the need to integrate complex nonpolynomial functions on the elements and preserving an optimal order of convergence. In this article, we develop for the first time, the VEM for parabolic problems on polygonal meshes, considering time‐dependent diffusion as our model problem. After presenting the scheme, we develop a theoretical analysis and show the practical behavior of the proposed method through a large array of numerical tests. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2110–2134, 2015     

A residual‐based Allen–Cahn phase field model for the mixture of incompressible fluid flows

The hydrodynamics of fluid mixtures is receiving more and more attention in many science and engineering applications. Within the techniques for dealing with front displacements and moving boundaries between different density and/or viscosity fluids, phase fields are a class of models in which a diffusive transition region is taken into account instead of a steep interface. Although these models have a physical motivation, they require the definition of extra parameters. In order to make it less parameter dependent, the classic Allen–Cahn phase field model is modified, exploring its similarities with residual‐based discontinuity‐capturing schemes, making the phase field equation dependent on its own residual. We solve the coupling between incompressible viscous fluid flow and the phase field advective–diffusive–reactive transport to simulate the main processes in interface tension and/or buoyancy driven problems. For the solution of the Navier–Stokes and transport equations, we use a stabilized finite element formulation. The implementation has been performed using the libMesh finite element library, written in C++ , which provides support for adaptive mesh refinement and coarsening. A chemical convection benchmark problem is used to validate the proposed model, and then we solve two bubble interaction problems. Copyright © 2014 John Wiley & Sons, Ltd.     

On the accuracy of a nonlinear finite volume method for the solution of diffusion problems using different interpolations strategies

In this paper, we consider a nonlinear finite volume method to solve the steady‐state diffusion equation in nonhomogeneous and non‐isotropic media. The method is nonlinear even if the original problem is linear. In its original form, the scheme is monotone, because the coefficient matrix is monotone under certain assumptions and, as a consequence, whenever the analytic operator demands, it preserves the positivity of numerical solutions. On the other hand, the scheme is unable to reproduce piecewise linear solutions exactly. In order to recover this interesting feature, we use two different interpolation strategies. In this case, even though we are unable to prove monotonicity, we show some numerical evidences that the combined method has an improved behavior, producing second order accurate solutions, even for nonhomogeneous and strongly anisotropic media. Copyright © 2013 John Wiley & Sons, Ltd.