A Raviart–Thomas mixed finite element scheme for the two‐dimensional three‐temperature heat conduction problems

For the two‐dimensional three‐temperature radiative heat conduction problem appearing in the inertial confinement numerical stimulations, we choose the Freezing coefficient method to linearize the nonlinear equations, and initially apply the well‐known mixed finite element scheme with the lowest order Raviart–Thomas element associated with the triangulation to the linearized equations, and obtain the convergence with one order with respect to the space direction for the temperature and flux function approximations, and design a simple but efficient algorithm for the discrete system. Three numerical examples are displayed. The former two verify theoretical results and show the super‐convergence for temperature and flux functions at the barycenter of the element, which is helpful for solving the radiative heat conduction problems. The third validates the robustness of this scheme with small energy conservative error and one order convergence for the time discretization. Copyright © 2016 John Wiley & Sons, Ltd.     

An adaptive non‐conforming finite‐element method for Reissner–Mindlin plates

Adaptive algorithms are important tools for efficient finite‐element mesh design. In this paper, an error controlled adaptive mesh‐refining algorithm is proposed for a non‐conforming low‐order finite‐element method for the Reissner–Mindlin plate model. The algorithm is controlled by a reliable and efficient residual‐based a posteriori error estimate, which is robust with respect to the plate's thickness. Numerical evidence for this and the efficiency of the new algorithm is provided in the sense that non‐optimal convergence rates are optimally improved in our numerical experiments. Copyright © 2003 John Wiley & Sons, Ltd.     

A face‐centred finite volume method for second‐order elliptic problems

This work proposes a novel finite volume paradigm, ie, the face‐centred finite volume (FCFV) method. Contrary to the popular vertex and cell‐centred finite volume methods, the novel FCFV defines the solution on the mesh faces (edges in two dimensions) to construct locally conservative numerical schemes. The idea of the FCFV method stems from a hybridisable discontinuous Galerkin formulation with constant degree of approximation, and thus inheriting the convergence properties of the classical hybridisable discontinuous Galerkin. The resulting FCFV features a global problem in terms of a piecewise constant function defined on the faces of the mesh. The solution and its gradient in each element are then recovered by solving a set of independent element‐by‐element problems. The mathematical formulation of FCFV for Poisson and Stokes equation is derived, and numerical evidence of optimal convergence in two dimensions and three dimensions is provided. Numerical examples are presented to illustrate the accuracy, efficiency, and robustness of the proposed methodology. The results show that, contrary to other finite volume methods, the accuracy of the FCFV method is not sensitive to mesh distortion and stretching. In addition, the FCFV method shows its better performance, accuracy, and robustness using simplicial elements, facilitating its application to problems involving complex geometries in three dimensions.     

High‐order space‐time finite element schemes for acoustic and viscodynamic wave equations with temporal decoupling

We revisit a method originally introduced by Werder et al. (in Comput. Methods Appl. Mech. Engrg., 190:6685–6708, 2001) for temporally discontinuous Galerkin FEMs applied to a parabolic partial differential equation. In that approach, block systems arise because of the coupling of the spatial systems through inner products of the temporal basis functions. If the spatial finite element space is of dimension D and polynomials of degree r are used in time, the block system has dimension (r + 1)D and is usually regarded as being too large when r > 1. Werder et al. found that the space‐time coupling matrices are diagonalizable over for r ?100, and this means that the time‐coupled computations within a time step can actually be decoupled. By using either continuous Galerkin or spectral element methods in space, we apply this DG‐in‐time methodology, for the first time, to second‐order wave equations including elastodynamics with and without Kelvin–Voigt and Maxwell–Zener viscoelasticity. An example set of numerical results is given to demonstrate the favourable effect on error and computational work of the moderately high‐order (up to degree 7) temporal and spatio‐temporal approximations, and we also touch on an application of this method to an ambitious problem related to the diagnosis of coronary artery disease. Copyright © 2014 The Authors. International Journal for Numerical Methods in Engineering published by John Wiley & Sons Ltd.     

An anisotropic Zienkiewicz–Zhu‐type error estimator for 3D applications

We extend the anisotropic Zienkiewicz–Zhu a posteriori error estimator of (Proceedings of the ENUMATH‐2009, Uppsala, Sweden, 29 June–3 July 2009) to three dimensions. Like the standard Zienkiewicz–Zhu estimator, the proposed estimator is designed to be independent of the problem at hand, is cheap to compute and easy to implement. In contrast to the standard Zienkiewicz–Zhu estimator, the elementwise counterpart of the proposed estimator explicitly takes into account the geometrical properties of the actual tetrahedron. Thus, in a wide variety of applications, the estimator is able to detect the anisotropic features exhibited by the solution of the governing equations. A metric‐based optimization procedure, rigorously addressed, drives the adaptation of the mesh. It is shown numerically to yield quasi‐optimal triangulations, dictating the accuracy‐vs‐number of elements behaviour. Despite being heuristic to some extent, in practice the overall anisotropic adaptation procedure turns out to be effective. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical simulation of flapping wings using a panel method and a high‐order Navier–Stokes solver

The design of efficient flapping wings for human engineered micro aerial vehicles (MAVs) has long been an elusive goal, in part because of the large size of the design space. One strategy for overcoming this difficulty is to use a multifidelity simulation strategy that appropriately balances computation time and accuracy. We compare two models with different geometric and physical fidelity. The low‐fidelity model is an inviscid doublet lattice method with infinitely thin lifting surfaces. The high‐fidelity model is a high‐order accurate discontinuous Galerkin Navier–Stokes solver, which uses an accurate representation of the flapping wing geometry. To compare the performance of the two methods, we consider a model flapping wing with an elliptical planform and an analytically prescribed spanwise wing twist, at size scales relevant to MAVs. Our results show that in many cases, including those with mild separation, low‐fidelity simulations can accurately predict integrated forces, provide insight into the flow structure, indicate regions of likely separation, and shed light on design–relevant quantities. But for problems with significant levels of separation, higher‐fidelity methods are required to capture the details of the flow field. Inevitably high‐fidelity simulations are needed to establish the limits of validity of the lower fidelity simulations.Copyright © 2012 John Wiley & Sons, Ltd.     

A hybrid discontinuous in space and time Galerkin method for wave propagation problems

Medium‐frequency regime and multi‐scale wave propagation problems have been a subject of active research in computational acoustics recently. New techniques have attempted to overcome the limitations of existing discretization methods that tend to suffer from dispersion. One such technique, the discontinuous enrichment method, incorporates features of the governing partial differential equation in the approximation, in particular, the solutions of the homogeneous form of the equation. Here, based on this concept and by extension of a conventional space–time finite element method, a hybrid discontinuous Galerkin method (DGM) for the numerical solution of transient problems governed by the wave equation in two and three spatial dimensions is described. The discontinuous formulation in both space and time enables the use of solutions to the homogeneous wave equation in the approximation. In this contribution, within each finite element, the solutions in the form of polynomial waves are employed. The continuity of these polynomial waves is weakly enforced through suitably chosen Lagrange multipliers. Results for two‐dimensional and three‐dimensional problems, in both low‐frequency and medium‐frequency regimes, show that the proposed DGM outperforms the conventional space–time finite element method. Copyright © 2014 John Wiley & Sons, Ltd.     

A fracture framework for Euler–Bernoulli beams based on a full discontinuous Galerkin formulation/extrinsic cohesive law combination

A new full Discontinuous Galerkin discretization of Euler–Bernoulli beam is presented. The main interest of this framework is its ability to simulate fracture problems by inserting a cohesive zone model in the formulation. With a classical Continuous Galerkin method, the use of the cohesive zone model is difficult because inserting a cohesive element between bulk elements is not straightforward. On one hand if the cohesive element is inserted at the beginning of the simulation, there is a modification of the structure stiffness and on the other hand inserting the cohesive element during the simulation requires modification of the mesh during computation. These drawbacks are avoided with the presented formulation as the structure is discretized in a stable and consistent way with full discontinuous elements and inserting cohesive elements during the simulation becomes straightforward. A new cohesive law based on the resultant stresses (bending moment and membrane) of the thin structure discretization is also presented. This model allows propagating fracture while avoiding through‐the‐thickness integration of the cohesive law. Tests are performed to show that the proposed model releases, during the fracture process, an energy quantity equal to the fracture energy for any combination of tension‐bending loadings. Copyright © 2010 John Wiley & Sons, Ltd.     

A high‐order immersed boundary discontinuous‐Galerkin method for Poisson's equation with discontinuous coefficients and singular sources

We adopt a numerical method to solve Poisson's equation on a fixed grid with embedded boundary conditions, where we put a special focus on the accurate representation of the normal gradient on the boundary. The lack of accuracy in the gradient evaluation on the boundary is a common issue with low‐order embedded boundary methods. Whereas a direct evaluation of the gradient is preferable, one typically uses post‐processing techniques to improve the quality of the gradient. Here, we adopt a new method based on the discontinuous‐Galerkin (DG) finite element method, inspired by the recent work of [A.J. Lew and G.C. Buscaglia. A discontinuous‐Galerkin‐based immersed boundary method. International Journal for Numerical Methods in Engineering, 76:427‐454, 2008]. The method has been enhanced in two aspects: firstly, we approximate the boundary shape locally by higher‐order geometric primitives. Secondly, we employ higher‐order shape functions within intersected elements. These are derived for the various geometric features of the boundary based on analytical solutions of the underlying partial differential equation. The development includes three basic geometric features in two dimensions for the solution of Poisson's equation: a straight boundary, a circular boundary, and a boundary with a discontinuity. We demonstrate the performance of the method via analytical benchmark examples with a smooth circular boundary as well as in the presence of a singularity due to a re‐entrant corner. Results are compared to a low‐order extended finite element method as well as the DG method of [1]. We report improved accuracy of the gradient on the boundary by one order of magnitude, as well as improved convergence rates in the presence of a singular source. In principle, the method can be extended to three dimensions, more complicated boundary shapes, and other partial differential equations. Copyright © 2014 John Wiley & Sons, Ltd.     

Implicit characteristic Galerkin method for convection–diffusion equations

This paper presents a characteristic Galerkin finite element method with an implicit algorithm for solving multidimensional, time‐dependent convection–diffusion equations. The method is formulated on the basis of the combination of both the precise and the implicit numerical integration procedures aiming to reference particles. The precise integration procedure with a 2^N algorithm is taken as a tool to determine the material (Lagrangian) derivative of the convective function in the operator splitting procedure. The stability analysis of the algorithm and numerical results illustrate good performance of the present method in stability and accuracy. Copyright © 2000 John Wiley & Sons, Ltd.     

A modification of the Petrov–Galerkin method for the transient convection–diffusion equation

A variation of the Petrov–Galerkin method of solution of a partial differential equation is presented in which the weight function applied to the time derivative term of the transient convection–diffusion equation is different from the weight function applied to the special derivatives. This allows for the formulation of fourth-order explicit and centred difference schemes. Comparison with analytic solutions show that these methods are able to capture steep wave fronts. The ability of the explicit method to capture wave fronts increases as the amount of convective transport increases.     

An enhanced cell‐based smoothed finite element method for the analysis of Reissner–Mindlin plate bending problems involving distorted mesh

In this work, an enhanced cell‐based smoothed finite element method (FEM) is presented for the Reissner–Mindlin plate bending analysis. The smoothed curvature computed by a boundary integral along the boundaries of smoothing cells in original smoothed FEM is reformulated, and the relationship between the original approach and the present method in curvature smoothing is established. To improve the accuracy of shear strain in a distorted mesh, we span the shear strain space over the adjacent element. This is performed by employing an edge‐based smoothing technique through a simple area‐weighted smoothing procedure on MITC4 assumed shear strain field. A three‐field variational principle is utilized to develop the mixed formulation. The resultant element formulation is further reduced to a displacement‐based formulation via an assumed strain method defined by the edge‐smoothing technique. As the result, a new formulation consisting of smoothed curvature and smoothed shear strain interpolated by the standard transverse displacement/rotation fields and smoothing operators can be shown to improve the solution accuracy in cell‐based smoothed FEM for Reissner–Mindlin plate bending analysis. Several numerical examples are presented to demonstrate the accuracy of the proposed formulation.Copyright © 2013 John Wiley & Sons, Ltd.     

A rheological interface model and its space–time finite element formulation for fluid–structure interaction

This contribution discusses extended physical interface models for fluid–structure interaction problems and investigates their phenomenological effects on the behavior of coupled systems by numerical simulation. Besides the various types of friction at the fluid–structure interface the most interesting phenomena are related to effects due to additional interface stiffness and damping. The paper introduces extended models at the fluid–structure interface on the basis of rheological devices (Hooke, Newton, Kelvin, Maxwell, Zener). The interface is decomposed into a Lagrangian layer for the solid‐like part and an Eulerian layer for the fluid‐like part. The mechanical model for fluid–structure interaction is based on the equations of rigid body dynamics for the structural part and the incompressible Navier–Stokes equations for viscous flow. The resulting weighted residual form uses the interface velocity and interface tractions in both layers in addition to the field variables for fluid and structure. The weak formulation of the whole coupled system is discretized using space–time finite elements with a discontinuous Galerkin method for time‐integration leading to a monolithic algebraic system. The deforming fluid domain is taken into account by deformable space–time finite elements and a pseudo‐structure approach for mesh motion. The sensitivity of coupled systems to modification of the interface model and its parameters is investigated by numerical simulation of flow induced vibrations of a spring supported fluid‐immersed cylinder. It is shown that the presented rheological interface model allows to influence flow‐induced vibrations. Copyright © 2010 John Wiley & Sons, Ltd.     

Tensor‐based Gauss–Jacobi numerical integration for high‐order mass and stiffness matrices

In this work, we choose the points and weights of the Gauss–Jacobi, Gauss–Radau–Jacobi and Gauss–Lobatto–Jacobi quadrature rules that optimize the number of operations for the mass and stiffness matrices of the high‐order finite element method. The procedure is particularly applied to the mass and stiffness matrices using the tensor‐based nodal and modal shape functions given in (Int. J. Numer. Meth. Engng 2007; 71 (5):529–563). For square and hexahedron elements, we show that it is possible to use tensor product of the 1D mass and stiffness matrices for the Poisson and elasticity problem. For the triangular and tetrahedron elements, an analogous analysis given in (Int. J. Numer. Meth. Engng 2005; 63 (2):1530–1558) was considered for the selection of the optimal points and weights for the stiffness matrix coefficients for triangles and mass and stiffness matrices for tetrahedra. Copyright © 2009 John Wiley & Sons, Ltd.     

A space–time discontinuous Galerkin method for the solution of the wave equation in the time domain

In recent years, the focus of research in the field of computational acoustics has shifted to the medium frequency regime and multiscale wave propagation. This has led to the development of new concepts including the discontinuous enrichment method. Its basic principle is the incorporation of features of the governing partial differential equation in the approximation. In this contribution, this concept is adapted for the simulation of transient problems governed by the wave equation. We present a space–time discontinuous Galerkin method with Lagrange multipliers, where the shape approximation in space and time is based on solutions of the homogeneous wave equation. The use of hierarchical wave‐like basis functions is enabled by means of a variational formulation that allows for discontinuities in both the spatial and the temporal discretizations. Numerical examples in one space dimension demonstrate the outstanding performance of the proposed method compared with conventional space–time finite element methods. Copyright © 2008 John Wiley & Sons, Ltd.     

Interface handling for three‐dimensional higher‐order XFEM‐computations in fluid–structure interaction

Three‐dimensional higher‐order eXtended finite element method (XFEM)‐computations still pose challenging computational geometry problems especially for moving interfaces. This paper provides a method for the localization of a higher‐order interface finite element (FE) mesh in an underlying three‐dimensional higher‐order FE mesh. Additionally, it demonstrates, how a subtetrahedralization of an intersected element can be obtained, which preserves the possibly curved interface and allows therefore exact numerical integration. The proposed interface algorithm collects initially a set of possibly intersecting elements by comparing their ‘eXtended axis‐aligned bounding boxes’. The intersection method is applied to a highly reduced number of intersection candidates. The resulting linearized interface is used as input for an elementwise constrained Delaunay tetrahedralization, which computes an appropriate subdivision for each intersected element. The curved interface is recovered from the linearized interface in the last step. The output comprises triangular integration cells representing the interface and tetrahedral integration cells for each intersected element. Application of the interface algorithm currently concentrates on fluid–structure interaction problems on low‐order and higher‐order FE meshes, which may be composed of any arbitrary element types such as hexahedra, tetrahedra, wedges, etc. Nevertheless, other XFEM‐problems with explicitly given interfaces or discontinuities may be tackled in addition. Multiple structures and interfaces per intersected element can be handled without any additional difficulties. Several parallelization strategies exist depending on the desired domain decomposition approach. Numerical test cases including various geometrical exceptions demonstrate the accuracy, robustness and efficiency of the interface handling. Copyright © 2009 John Wiley & Sons, Ltd.     

FE/BE coupling for an acoustic fluid–structure interaction problem. Residual a posteriori error estimates

In this paper, we developed an a posteriori error analysis of a coupling of finite elements and boundary elements for a fluid–structure interaction problem in two and three dimensions. This problem is governed by the acoustic and the elastodynamic equations in time‐harmonic vibration. Our methods combined integral equations for the exterior fluid and FEMs for the elastic structure. It is well‐known that because of the reduction of the boundary value problem to boundary integral equations, the solution is not unique in general. However, because of superposition of various potentials, we consider a boundary integral equation that is uniquely solvable and avoids the irregular frequencies of the negative Laplacian operator of the interior domain. In this paper, two stable procedures were considered; one is based on the nonsymmetric formulation and the other is based on a symmetric formulation. For both formulations, we derived reliable residual a posteriori error estimates. From the estimators we computed local error indicators that allowed us to develop an adaptive mesh refinement strategy. For the two‐dimensional case we performed an adaptive algorithm on triangles, and for the three‐dimensional case we used hanging nodes on hexahedrons. Numerical experiments underline our theoretical results. Copyright © 2011 John Wiley & Sons, Ltd.     

A second‐order decoupled implicit/explicit method of the 3D primitive equations of ocean II: finite element spatial discretization

A fully discrete second‐order decoupled implicit/explicit method is proposed for solving 3D primitive equations of ocean in the case of Dirichlet boundary conditions on the side, where a second‐order decoupled implicit/explicit scheme is used for time discretization, and a finite element method based on the P₁(P₁) ? P₁?P₁(P₁) elements for velocity, pressure and density is used for spatial discretization of these primitive equations. Optimal H¹?L²?H¹ error estimates for numerical solution and an optimal L² error estimate for are established under the convergence condition of 0 < h≤β₁,0 < τ≤β₂, and τ≤β₃h for some positive constants β₁,β₂, and β₃. Furthermore, numerical computations show that the H¹?L²?H¹ convergence rate for numerical solution is of O(h + τ²) and an L² convergence rate for is O(h²+τ²) with the assumed convergence condition, where h is a mesh size and τ is a time step size. More practical calculations are performed as a further validation of the numerical method. Copyright © 2016 John Wiley & Sons, Ltd.