The adaptive immersed interface finite element method for elliptic and Maxwell interface problems

We propose self-adaptive finite element methods with error control for solving elliptic and electromagnetic problems with discontinuous coefficients. The meshes in the methods do not need to fit the interfaces. New error indicators are introduced to control the error due to non-body-fitted meshes. Flexible h-adaptive strategies are developed, which can be systematically extended to a large class of interface problems. Extensive numerical experiments are performed to support the theoretical results and to show the competitive behavior of the adaptive algorithm even for interfaces involving corner or tip singularities.     

Morley type non‐C0 nonconforming rectangular plate finite elements on anisotropic meshes

In this article, two Morley type non‐C⁰ nonconforming rectangular finite elements are discussed to numerically solve the fourth order plate bending problem under anisotropic meshes. The optimal anisotropic interpolation error and consistency error estimates are obtained by using some novel approaches. Some numerical tests are given to confirm the theoretical analysis. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Unstructured finite volume method for matrix free explicit solution of stress-strain fields in two dimensional problems with curved boundaries in equilibrium condition

In this paper, a plane stress structural solver which uses a matrix free unstructured finite volume method based on Galerkin approach is introduced for solution of weak form of two dimensional Cauchy equations on linear triangular element meshes. The developed shape function free Galerkin finite volume structural solver explicitly computes stresses and displacements in cartesian coordinate directions for the two dimensional solid mechanic problems in equilibrium condition. The accuracy of the introduced algorithm is assessed by comparison of computed results of two plane-stress cases with curved boundaries under uniformly distributed loads with available analytical solutions. The results of the introduced method are presented in terms of stress and strain contours and its effective parameters on convergence behaviour to equilibrium condition are assessed.     

A higher-order macroscopic model for pedestrian flows

This paper develops a higher-order macroscopic model of pedestrian crowd dynamics derived from fluid dynamics that consists of two-dimensional Euler equations with relaxation. The desired directional motion of pedestrians is determined by an Eikonal-type equation, which describes a problem that minimizes the instantaneous total walking cost from origin to destination. A linear stability analysis of the model demonstrates its ability to describe traffic instability in crowd flows. The algorithm to solve the macroscopic model is composed of a splitting technique introduced to treat the relaxation terms, a second-order positivity-preserving central-upwind scheme for hyperbolic conservation laws, and a fast-sweeping method for the Eikonal-type equation on unstructured meshes. To test the applicability of the model, we study a challenging pedestrian crowd flow problem of the presence of an obstruction in a two-dimensional continuous walking facility. The numerical results indicate the rationality of the model and the effectiveness of the computational algorithm in predicting the flux or density distribution and the macroscopic behavior of the pedestrian crowd flow. The simulation results are compared with those obtained by the two-dimensional Lighthill-Whitham-Richards pedestrian flow model with various model parameters, which further shows that the macroscopic model is able to correctly describe complex phenomena such as “stop-and-go waves” observed in empirical pedestrian flows.     

A brief survey on numerical methods for solving singularly perturbed problems

In the present paper, a brief survey on computational techniques for the different classes of singularly perturbed problems is given. This survey is a continuation of work performed earlier by the first author and contains the literature of the work done by the researchers during the years 2000-2009. However some older important relevant papers are also included in this survey. We also mentioned those papers which are not surveyed in the previous survey papers by the first author of this paper, see [Appl. Math. Comput. 30 (1989) 223-259, 130 (2002) 457-510, 134 (2003) 371-429] for details. Thus this survey paper contains a surprisingly large amount of literature on singularly perturbed problems and indeed can serve as an introduction to some of the ideas and methods for the singular perturbation problems.     

Finite element simulation of compressible particle-laden gas flows

A macroscopic two-fluid model of compressible particle-laden gas flows is considered. The governing equations are discretized by a high-resolution finite element method based on algebraic flux correction. A multidimensional limiter of TVD type is employed to constrain the local characteristic variables for the continuous gas phase and conservative fluxes for a suspension of solid particles. Special emphasis is laid on the efficient computation of steady state solutions at arbitrary Mach numbers. To avoid stability restrictions and convergence problems, the characteristic boundary conditions are imposed weakly and treated in a fully implicit manner. A two-way coupling via the interphase drag force is implemented using operator splitting. The Douglas-Rachford scheme is found to provide a robust treatment of the interphase exchange terms within the framework of a fractional-step solution strategy. Two-dimensional simulation results are presented for a moving shock wave and for a steady nozzle flow.     

A high accuracy hybrid method for two-dimensional Navier–Stokes equations

A dual-mesh hybrid numerical method is proposed for high Reynolds and high Rayleigh number flows. The scheme is of high accuracy because of the use of a fourth-order finite-difference scheme for the time-dependent convection and diffusion equations on a non-uniform mesh and a fast Poisson solver DFPS2H based on the HODIE finite-difference scheme and algorithm HFFT [R.A. Boisvert, Fourth order accurate fast direct method for the Helmholtz equation, in: G. Birkhoff, A. Schoenstadt (Eds.), Elliptic Problem Solvers II, Academic Press, Orlando, FL, 1984, pp. 35–44] for the stream function equation on a uniform mesh. To combine the fast Poisson solver DFPS2H and the high-order upwind-biased finite-difference method on the two different meshes, Chebyshev polynomials have been used to transfer the data between the uniform and non-uniform meshes. Because of the adoption of a hybrid grid system, the proposed numerical model can handle the steep spatial gradients of the dependent variables by using very fine resolutions in the boundary layers at reasonable computational cost. The successful simulation of lid-driven cavity flows and differentially heated cavity flows demonstrates that the proposed numerical model is very stable and accurate within the range of applicability of the governing equations.     

An edge-based pressure stabilization technique for finite elements on arbitrarily anisotropic meshes

In this paper, we analyze a stabilized equal-order finite element approximation for the Stokes equations on anisotropic meshes. In particular, we allow arbitrary anisotropies in a subdomain, for example, along the boundary of the domain, with the only condition that a maximum angle is fulfilled in each element. This discretization is motivated by applications on moving domains as arising, for example, in fluid-structure interaction or multiphase-flow problems. To deal with the anisotropies, we define a modification of the original continuous interior penalty stabilization approach. We show analytically the discrete stability of the method and convergence of order in the energy norm and in the L²-norm of the velocities. We present numerical examples for a linear Stokes problem and for a nonlinear fluid-structure interaction problem, which substantiate the analytical results and show the capabilities of the approach.     

惯性约束聚变内爆中基于多块结构网格的高效辐射扩散并行算法

针对三维球形靶丸内爆高效模拟需求和传统笛卡尔正交网格上辐射加源困难的问题, 发展一种多块结构非正交网格生成方法, 并基于此种计算网格提出高效的三维扩散格式并行算法, 将其应用于辐射流体方程组的求解和三维内爆不对称性的数值模拟, 数值结果显示了算法的有效性。并行性能测试显示该算法可扩展到5400个核上, 并行效率达到69%。     

Parallel direct Poisson solver for discretisations with one Fourier diagonalisable direction

In the context of time-accurate numerical simulation of incompressible flows, a Poisson equation needs to be solved at least once per time-step to project the velocity field onto a divergence-free space. Due to the non-local nature of its solution, this elliptic system is one of the most time consuming and difficult to parallelise parts of the code.