A local discontinuous Galerkin method for pattern formation dynamical model in polymerizing action flocks

In this paper, we apply local discontinuous Galerkin methods to the pattern formation dynamical model in polymerizing action flocks. Optimal error estimates for the density and filament polarization in different norms are established. We use a semi-implicit spectral deferred correction time method for time discretization, which allows a relative large time step and avoids computation of a Jacobian matrix. Numerical experiments are presented to verify the theoretical analysis and to show the capability for simulations of action wave formation.

     

A local discontinuous Galerkin method for the pattern formation dynamical model in polymerizing action flocks

In this paper, we apply local discontinuous Galerkin methods to the pattern formation dynamical model in polymerizing action flocks. Optimal error estimates for the density and filament polarization in different norms are established. We use a semi-implicit spectral deferred correction time method for time discretization,which allows a relative large time step and avoids computation of a Jacobian matrix. Numerical experiments are presented to verify the theoretical analysis and to show the capab...     

Mixed finite element analysis for dissipative SRLW equations with damping term

In this paper a mixed finite element (MFE) formulation is proposed for the initial-boundary value problem of dissipative symmetric regularized long wave (SRLW) equations with damping. Existence and uniqueness of its generalized solution and of the fully discrete mixed finite element solution are proved. Error estimates based on energy methods are given. Numerical experiments verify the theoretical analysis.     

Nonconforming finite element approximation of crystalline microstructure

We consider a class of nonconforming finite element approximations of a simply laminated microstructure which minimizes the nonconvex variational problem for the deformation of martensitic crystals which can undergo either an orthorhombic to monoclinic (double well) or a cubic to tetragonal (triple well) transformation. We first establish a series of error bounds in terms of elastic energies for the approximation of derivatives of the deformation in the direction tangential to parallel layers of the laminate, for the approximation of the deformation, for the weak approximation of the deformation gradient, for the approximation of volume fractions of deformation gradients, and for the approximation of nonlinear integrals of the deformation gradient. We then use these bounds to give corresponding convergence rates for quasi-optimal finite element approximations.

     

Adaptive finite element analysis of structures under transient dynamic loading using modal superposition

In this paper, the error estimation and adaptive strategy developed for the linear elastodynamic problem under transient dynamic loading based on the Z–Z criterion is utilized for 2D and plate bending problems. An automatic mesh generator based on “growth meshing” is utilized effectively for adaptive mesh refinement. Optimal meshes are obtained iteratively corresponding to the prescribed domain discretization error limit and for a chosen number of basis modes satisfying modal truncation errors. Numerous examples show the effectiveness of the integrated approach in achieving the target accuracy in finite element transient dynamic analysis.     

An enhanced finite volume method to model 2D linear elastic structures

This paper details the evaluation and enhancement of the vertex-centred finite volume method for the purpose of modelling linear elastic structures undergoing bending. A matrix-free edge-based finite volume procedure is discussed and compared with the traditional isoparametric finite element method via application to a number of test-cases. It is demonstrated that the standard finite volume approach exhibits similar disadvantages to the linear Q4 finite element formulation when modelling bending. An enhanced finite volume approach is proposed to circumvent this and a rigorous error analysis conducted. It is demonstrated that the developed finite volume method is superior to both standard finite volume and Q4 finite element methods, and provides a practical alternative to the analysis of bending-dominated solid mechanics problems.     

Rate of weak convergence of the finite element method for the stochastic heat equation with additive noise

The stochastic heat equation driven by additive noise is discretized in the spatial variables by a standard finite element method. The weak convergence of the approximate solution is investigated and the rate of weak convergence is found to be twice that of strong convergence. M. Kovács and S. Larsson supported by the Swedish Research Council (VR). Part of this work was done at Institut Mittag-Leffler. S. Larsson supported by the Swedish Foundation for Strategic Research (SSF) through GMMC, the Gothenburg Mathematical Modelling Centre.     

Discrete finite element characterizations of source fields for volume and boundary constraints in electromagnetic problems

The consideration of electromagnetic field sources in potential formulations necessitates the definition of source fields. Such source fields are first defined for both volume and boundary constraints in static electromagnetic models. Then, automatic procedures are proposed to conveniently and efficiently characterize discrete source fields, with regard to their use in finite element formulations, their supports, their direct expression requiring no pre-computation, and their associated constraints. Two application examples are proposed to illustrate the approach.     

Two-grid methods for characteristic finite volume element solution of semilinear convection-diffusion equations

Two-grid methods for characteristic finite volume element solutions are presented for a kind of semilinear convection-dominated diffusion equations. The methods are based on the method of characteristics, two-grid method and the finite volume element method. The nonsymmetric and nonlinear iterations are only executed on the coarse grid (with grid size H). And the fine-grid solution (with grid size h) can be obtained by a single symmetric and linear step. It is proved that the coarse grid can be much coarser than the fine grid. The two-grid methods achieve asymptotically optimal approximation as long as the mesh sizes satisfy H = O(h^1/3).     

A fully discrete stabilized finite element method for the time-dependent Navier-Stokes equations

In this article, we consider a fully discrete stabilized finite element method based on two local Gauss integrations for the two-dimensional time-dependent Navier-Stokes equations. It focuses on the lowest equal-order velocity-pressure pairs. Unlike the other stabilized method, the present approach does not require specification of a stabilization parameter or calculation of higher-order derivatives, and always leads to a symmetric linear system. The Euler semi-implicit scheme is used for the time discretization. It is shown that the proposed fully discrete stabilized finite element method results in the optimal order bounds for the velocity and pressure.     

The MMOC and MMOCAA schemes for the finite volume element method of convection-diffusion problems

We present and analyze the modified method of characteristics (MMOC) and the modified method of characteristics with adjusted advection (MMOCAA) for the finite volume element (FVE) method of convection-diffusion problems. These two schemes maintain the advantages of both the MMOC and the FVE method. And the MMOCAA scheme discussed herein conserves the conservation law globally at a minor additional computational cost. Optimal-order error estimates in the H¹-norm are proved for these schemes. A numerical example is presented to confirm the estimates.     

A projection-based stabilized finite element method for steady-state natural convection problem

We formulate a projection-based stabilization finite element technique for solving steady-state natural convection problems. In particular, we consider heat transport through combined solid and fluid media. This stabilization does not act on the large flow structures. Based on the projection stabilization idea, finite element error analysis of the problem is investigated and optimal errors for the velocity, temperature and pressure are established. We also present some numerical tests which both verify the theoretical predictions and demonstrate the method?s promise.     

两相渗流驱动问题的体积有限元数值计算和分析

本文在一般的三角形剖分上对两相渗流驱动提出了全离散体积有限元 ,并分析了带有弥散项时格式的收敛性 ,得到H1 模的最优估计 .     

Robust error estimates for the finite element approximation of elliptic optimal control problems

In this paper, we consider the finite element approximation of an elliptic optimal control problem. Based on an assumption on the adjoint state of the continuous problem with a small parameter, which represents a regularization of the bang-bang type control problem, we derive robust a priori error estimates for optimal control and state and a posteriori error estimate is also presented. Numerical experiments confirm our theoretical results.     

On error estimation of finite element approximations to the elliptic equations in nonconvex polygonal domains

Numerical verification methods, so-called Nakao's methods, on existence or uniqueness of solutions to PDEs have been developed by Nakao and his group including the authors. They are based on the error estimation of approximate solutions which are mainly computed by FEM.     

A novel finite element model for helical springs

A general and accurate finite element model for helical springs subject to axial loads (extension or/and torsion) is developed in this paper. Due to the establishment of precise boundary conditions, only a slice of the wire cross-section needs to be modelled; hence, more accurate results can be achieved. An example application to a circular cross-sectional spring is analysed in detail.     

On the use of physical spline finite element method for acoustic scattering

This paper presents the comparison of physical spline finite element method (PSFEM), in which differential equations are incorporated into interpolations of basic elements, with least-squares finite element method (LSFEM) and mixed Galerkin finite element method (MGFEM) on the numerical solution of one dimensional Helmholtz equation applied to an acoustic scattering problem. Firstly, all three methods are explained in detail and then it is shown that PSFEM reaches higher precision in a shorter time with fewer nodes than the other methods. It is also observed that this method is well suited for high frequency acoustic problems. Consequently, the results of PSFEM point out better efficiency in terms of number of unknowns and accuracy level.     

Convergence analysis of an upwind mixed element method for advection diffusion problems

We consider a upwinding mixed element method for a system of first order partial differential equations resulting from the mixed formulation of a general advection diffusion problem. The system can be used to model the transport of a contaminant carried by a flow. We use the lowest order Raviart-Thomas mixed finite element space. We show the first order convergence both for concentration and concentration flux in L²(Ω).     

Inverse mode problems for the finite element model of a vibrating rod

The inverse mode problems for the finite element model of an axially vibrating rod are formulated and solved. It is known that for the finite element model, based on linear shape functions, of the rod, the mass and stiffness matrices are both tridiagonal. It is shown that the finite element model of the rod can be constructed from two eigenvalues, their corresponding eigenvectors and the total mass of the rod. The necessary and sufficient conditions for the construction of a physically realizable rod with positive mass and stiffness elements from two eigenpairs and the total mass of the rod are established. If these conditions are satisfied, then the construction of the model is unique.