In this paper, a conservative parallel iteration scheme is constructed to solve nonlinear diffusion equations on unstructured polygonal meshes. The design is based on two main ingredients: the first is that the parallelized domain decomposition is embedded into the nonlinear iteration; the second is that prediction and correction steps are applied at subdomain interfaces in the parallelized domain decomposition method. A new prediction approach is proposed to obtain an efficient conservative parallel finite volume scheme. The numerical experiments show that our parallel scheme is second-order accurate, unconditionally stable, conservative and has linear parallel speed-up.

We study the gradient superconvergence of bilinear finite volume element (FVE) solving the elliptic problems. First, a superclose weak estimate is established for the bilinear form of the FVE method. Then, we prove that the gradient approximation of the FVE solution has the superconvergence property:

where denotes the average gradient on elements containing point P and S is the set of optimal stress points composed of the mesh points, the midpoints of edges and the centers of elements.

This paper presents a new approach to verify the accuracy of computational simulations. We develop mathematical theorems which can serve as robust a posteriori error estimation techniques to identify numerical pollution, check the performance of adaptive meshes, and verify numerical solutions. We demonstrate performance of this methodology on problems from flow thorough porous media. However, one can extend it to other models. We construct mathematical properties such that the solutions to Darcy and Darcy-Brinkman equations satisfy them. The mathematical properties include the total minimum mechanical power, minimum dissipation theorem, reciprocal relation, and maximum principle for the vorticity. All the developed theorems have firm mechanical bases and are independent of numerical methods. So, these can be utilized for solution verification of finite element, finite volume, finite difference, lattice Boltzmann methods and so forth. In particular, we show that, for a given set of boundary conditions, Darcy velocity has the minimum total mechanical power of all the kinematically admissible vector fields. We also show that a similar result holds for Darcy-Brinkman velocity. We then show for a conservative body force, the Darcy and Darcy-Brinkman velocities have the minimum total dissipation among their respective kinematically admissible vector fields. Using numerical examples, we show that the minimum dissipation and total mechanical power theorems can be utilized to identify pollution errors in numerical solutions. The solutions to Darcy and Darcy-Brinkman equations are shown to satisfy a reciprocal relation, which has the potential to identify errors in the numerical implementation of boundary conditions. It is also shown that the vorticity under both steady and transient Darcy-Brinkman equations satisfy maximum principles if the body force is conservative and the permeability is homogeneous and isotropic. A discussion on the nature of vorticity under steady and transient Darcy equations is also presented. Using several numerical examples, we will demonstrate the predictive capabilities of the proposed a posteriori techniques in assessing the accuracy of numerical solutions for a general class of problems, which could involve complex domains and general computational grids.

A moving mesh method is proposed for solving reaction-diffusion equations. The finite element method is used to solving the partial different equation system, and an efficient numerical scheme is applied to implement mesh moving. In the practical calculations, the moving mesh step and the problem equation solver are performed alternatively. Several numerical examples are presented, including the Gray-Scott, the Activator-Inhibitor and a case with a growing domain. It is illustrated numerically that the moving mesh methods costs much lower, compared with the numerical schemes on a fixed mesh. Even in the case of complex pattern dynamics described by the reaction-diffusion systems, the adapted meshes can capture the details successfully.

In this paper, we develop a novel energy-preserving wavelet collocation method for solving general multi-symplectic formulations of Hamiltonian PDEs. Based on the autocorrelation functions of Daubechies compactly supported scaling functions, the wavelet collocation method is conducted for spatial discretization. The obtained semi-discrete system is shown to be a finite-dimensional Hamiltonian system, which has an energy conservation law. Then, the average vector field method is used for time integration, which leads to an energy-preserving method for multi-symplectic Hamiltonian PDEs. The proposed method is illustrated by the nonlinear Schrödinger equation and the Camassa-Holm equation. Since differentiation matrix obtained by the wavelet collocation method is a cyclic matrix, we can apply Fast Fourier transform to solve equations in numerical calculation. Numerical experiments show the high accuracy, effectiveness and conservation properties of the proposed method.

We prove that a class of A-stable symplectic Runge–Kutta time semi-discretizations (including the Gauss–Legendre methods) applied to a class of semilinear Hamiltonian partial differential equations (PDEs) that are well posed on spaces of analytic functions with analytic initial data can be embedded into a modified Hamiltonian flow up to an exponentially small error. Consequently, such time semi-discretizations conserve the modified Hamiltonian up to an exponentially small error. The modified Hamiltonian is O(hp )-close to the original energy, where p is the order of the method and h is the time-step size. Examples of such systems are the semilinear wave equation, and the nonlinear Schrödinger equation with analytic nonlinearity and periodic boundary conditions. Standard Hamiltonian interpolation results do not apply here because of the occurrence of unbounded operators in the construction of the modified vector field. This loss of regularity in the construction can be taken care of by projecting the PDE to a subspace in which the operators occurring in the evolution equation are bounded, and by coupling the number of excited modes and the number of terms in the expansion of the modified vector field with the step size. This way we obtain exponential estimates of the form O(exp(–c/h 1/(1+q))) with c > 0 and q ⩾ 0; for the semilinear wave equation, q = 1, and for the nonlinear Schrödinger equation, q = 2. We give an example which shows that analyticity of the initial data is necessary to obtain exponential estimates.

AMG preconditioners are typically designed for partial differential equation solvers and divergence-interpolation in a moving mesh strategy. Here we introduce an AMG preconditioner to solve the unsteady Navier-Stokes equations by a moving mesh finite element method. A 4P1 – P1 element pair is selected based on the data structure of the hierarchy geometry tree and two-layer nested meshes in the velocity and pressure. Numerical experiments show the efficiency of our approach.

A generalised Hermite spectral method for Fisher's equation in genetics with different asymptotic solution behaviour at infinities is proposed, involving a fully discrete scheme using a second order finite difference approximation in the time. The convergence and stability of the scheme are analysed, and some numerical results demonstrate its efficiency and substantiate our theoretical analysis.

In this paper, we analyze a nonconforming finite element method for the computation of transmission eigenvalues and the corresponding eigenfunctions. The error estimates of the eigenvalue and eigenfunction approximation are given, respectively. Finally, some numerical examples are provided to validate the theoretical results.

In this paper, we extend using the Runge-Kutta discontinuous Galerkin method together with the front tracking method to simulate the compressible two-medium flow on unstructured meshes. A Riemann problem is constructed in the normal direction in the material interfacial region, with the goal of obtaining a compact, robust and efficient procedure to track the explicit sharp interface precisely. Extensive numerical tests including the gas-gas and gas-liquid flows are provided to show the proposed methodologies possess the capability of enhancing the resolutions nearby the discontinuities inside of the single medium flow and the interfacial vicinities of the two-medium flow in many occasions.

Semi-Lagrangian (S-L) methods have no CFL stability constraint, and are more stable than the Eulerian methods. In the literature, the S-L method for the level-set re-initialization equation was complicated, which may be unnecessary. Since the re-initialization procedure is auxiliary, we propose to use the first-order S-L scheme coupled with a projection technique to improve the accuracy at the grid points just adjacent to the interface. Standard second-order S-L method is used for evolving the level-set convection equation. The implementation is simple, including on the block-structured adaptive mesh. The efficiency of the S-L method is demonstrated by extensive numerical examples including passive convection of interfaces with corners/kinks/large deformation under given velocity fields, a geometrical flow with topological changes, simulations of bubble/ droplet dynamics in incompressible two-phase flows. In terms of accuracy it is comparable to the other existing methods.

In this paper, a method is proposed for extracting fracture parameters in anisotropic thermoelasticity cracking via interaction integral method within the framework of extended finite element method (XFEM). The proposed method is applied to linear thermoelastic crack problems. The numerical results of the stress intensity factors (SIFs) are presented and compared with those reported in related references. The good agreement of the results obtained by the developed method with those obtained by other numerical solutions proves the applicability of the proposed approach and confirms its capability of efficiently extracting thermoelasticity fracture parameters in anisotropic materials.

In this paper we develop explicit fast exponential Runge-Kutta methods for the numerical solutions of a class of parabolic equations. By incorporating the linear splitting technique into the explicit exponential Runge-Kutta schemes, we are able to greatly improve the numerical stability. The proposed numerical methods could be fast implemented through use of decompositions of compact spatial difference operators on a regular mesh together with discrete fast Fourier transform techniques. The exponential Runge-Kutta schemes are easy to be adopted in adaptive temporal approximations with variable time step sizes, as well as applied to stiff nonlinearity and boundary conditions of different types. Linear stabilities of the proposed schemes and their comparison with other schemes are presented. We also numerically demonstrate accuracy, stability and robustness of the proposed method through some typical model problems.

Fixed-point iterative sweeping methods were developed in the literature to efficiently solve static Hamilton-Jacobi equations. This class of methods utilizes the Gauss-Seidel iterations and alternating sweeping strategy to achieve fast convergence rate. They take advantage of the properties of hyperbolic partial differential equations (PDEs) and try to cover a family of characteristics of the corresponding Hamilton-Jacobi equation in a certain direction simultaneously in each sweeping order. Different from other fast sweeping methods, fixed-point iterative sweeping methods have the advantages such as that they have explicit forms and do not involve inverse operation of nonlinear local systems. In principle, it can be applied in solving very general equations using any monotone numerical fluxes and high order approximations easily. In this paper, based on the recently developed fifth order WENO schemes which improve the convergence of the classical WENO schemes by removing slight post-shock oscillations, we design fifth order fixed-point sweeping WENO methods for efficient computation of steady state solution of hyperbolic conservation laws. Especially, we show that although the methods do not have linear computational complexity, they converge to steady state solutions much faster than regular time-marching approach by stability improvement for high order schemes with a forward Euler time-marching.

High fidelity modeling and simulation of moderately dense sprays at relatively low cost is still a major challenge for many applications. For that purpose, we introduce a new multi-fluid model based on a two-size moment formalism in sections, which are size intervals of discretization. It is derived from a Boltzmann type equation taking into account drag, evaporation and coalescence, which are representative of the complex terms that arise in multi-physics environments. The closure of the model comes from a reconstruction of the distribution. A piecewise affine reconstruction in size is thoroughly analyzed in terms of stability and accuracy, a key point for a high-fidelity and reliable description of the spray. Robust and accurate numerical methods are then developed, ensuring the realizability of the moments. The model and method are proven to describe the spray with a high accuracy in size and size-conditioned variables, resorting to a lower number of sections compared to one size moment methods. Moreover, robustness is ensured with efficient and tractable algorithms despite the numerous couplings and various algebra thanks to a tailored overall strategy. This strategy is successfully tested on a difficult 2D unsteady case, which proves the efficiency of the modeling and numerical choices.

The goal of this work is to construct and study hybrid and multiplicative two-level overlapping Schwarz algorithms with standard coarse spaces for the almost incompressible linear elasticity and Stokes systems, discretized by mixed finite and spectral element methods with discontinuous pressures. Two different approaches are considered to solve the resulting saddle point systems: a) a preconditioned conjugate gradient (PCG) method applied to the symmetric positive definite reformulation of the almost incompressible linear elasticity system obtained by eliminating the pressure unknowns; b) a GMRES method with indefinite overlapping Schwarz preconditioner applied directly to the saddle point formulation of both the elasticity and Stokes systems. Condition number estimates and convergence properties of the proposed hybrid and multiplicative overlapping Schwarz algorithms are proven for the positive definite reformulation of almost incompressible elasticity. These results are based on our previous study [8] where only additive Schwarz preconditioners were considered for almost incompressible elasticity. Extensive numerical experiments with both finite and spectral elements show that the proposed overlapping Schwarz preconditioners are scalable, quasi-optimal in the number of unknowns across individual subdomains and robust with respect to discontinuities of the material parameters across subdomains interfaces. The results indicate that the proposed preconditioners retain a good performance also when the quasi-monotonicity assumption, required by the available theory, does not hold.

A comparative study of two classes of third-order implicit time integration schemes is presented for a third-order hierarchical WENO reconstructed discontinuous Galerkin (rDG) method to solve the 3D unsteady compressible Navier-Stokes equations: — 1) the explicit first stage, single diagonally implicit Runge-Kutta (ESDIRK3) scheme, and 2) the Rosenbrock-Wanner (ROW) schemes based on the differential algebraic equations (DAEs) of Index-2. Compared with the ESDIRK3 scheme, a remarkable feature of the ROW schemes is that, they only require one approximate Jacobian matrix calculation every time step, thus considerably reducing the overall computational cost. A variety of test cases, ranging from inviscid flows to DNS of turbulent flows, are presented to assess the performance of these schemes. Numerical experiments demonstrate that the third-order ROW scheme for the DAEs of index-2 can not only achieve the designed formal order of temporal convergence accuracy in a benchmark test, but also require significantly less computing time than its ESDIRK3 counterpart to converge to the same level of discretization errors in all of the flow simulations in this study, indicating that the ROW methods provide an attractive alternative for the higher-order time-accurate integration of the unsteady compressible Navier-Stokes equations.

We propose and analyse finite volume Godunov type methods based on discontinuous flux for a 2×2 system of non-linear partial differential equations proposed by Hadeler and Kuttler to model the dynamics of growing sandpiles generated by a vertical source on a flat bounded rectangular table. The problem considered here is the so-called partially open table problem where sand is blocked by a wall (of infinite height) on some part of the boundary of the table. The novelty here is the corresponding modification of boundary conditions for the standing and the rolling layers and generalization of the techniques of the well-balancedness proposed in [1]. Presence of walls may lead to unbounded or discontinuous surface flow density at equilibrium resulting in solutions with singularities propagating from the extreme points of the walls. A scheme has been proposed to approximate efficiently the Hamiltonians with the coefficients which can be unbounded and discontinuous. Numerical experiments are presented to illustrate that the proposed schemes detect these singularities in the equilibrium solutions efficiently and comparisons are made with the previously studied finite difference and Semi-Lagrangian approaches by Finzi Vita et al.

Multiscale modelling is a powerful technique, which allows for computational efficiency while retaining small-scale details when they are essential for understanding a finer behaviour of the studied system. In the case of materials modelling, one of the effective multiscaling concepts is domain partitioning, which implies the existence of an explicit interface between various material descriptions, for instance atomistic and continuum regions. When dynamic material behaviour is considered, the major problem for this technique is dealing with reflections of high frequency waves from the interface separating two scales. In this article, a new method is suggested, which overcomes this problem for the case of magnetisation dynamics. The introduction of a damping band at the interface between scales, which absorbs high frequency waves, is suggested. The idea is verified using a number of one-dimensional examples with fine/coarse scale discretisation of a continuum problem of spin wave propagation. This work is the first step towards establishing a reliable atomistic/continuum multiscale transition for the description of the evolution of magnetic properties of ferromagnets.

An optimal control problem is considered to find a stable surface traction, which minimizes the discrepancy between a given displacement field and its estimation. Firstly, the inverse elastic problem is constructed by variational inequalities, and a stable approximation of surface traction is obtained with Tikhonov regularization. Then a finite element discretization of the inverse elastic problem is analyzed. Moreover, the error estimation of the numerical solutions is deduced. Finally, a numerical algorithm is detailed and three examples in two-dimensional case illustrate the efficiency of the algorithm.