Analysis of a finite element method for pressure/potential formulation of elastoacoustic spectral problems

A finite element method to approximate the vibration modes of a structure enclosing an acoustic fluid is analyzed. The fluid is described by using simultaneously pressure and displacement potential variables, whereas displacement variables are used for the solid. A mathematical analysis of the continuous spectral problem is given. The problem is discretized on a simplicial mesh by using piecewise constant elements for the pressure and continuous piecewise linear finite elements for the other fields. Error estimates are settled for approximate eigenvalues and eigenfrequencies. Finally, implementation issues are discussed.

     

Finite element analysis of a holonomic elastic-plastic problem

Summary A thorough analysis of the finite element method is given for a holonomic elastic-plastic problem. An inequality of the Cea's lemma type is proved, which is the basis of error estimates for various finite element solutions. Difficulty caused by a non-differentiable term in the problem can be overcome by using two convergent procedures, an iterative procedure and a regularization procedure. An a-posteriori quantitative error estimate is derived for the regularized solution.The work was done while the author was at the Department of Mathematics, University of Maryland, College Park. The research was partially supported by the National Science Foundation grant CCR-88-20279     

Finite element formulation of thermal boundary layers

In this article we discuss the finite element discretization of the two-dimensional, incompressible, and turbulent boundary layers. The formulation of the momentum equation is essentially due to Baker and Soliman [1] with some modifications.The versatility and the accuracy of the method is established by considering several test cases. The predictions are satisfactory and compare favorably with alternative numerical techniques.     

Finite element approximation of the volume-matching problem

Summary Optimal orderH ¹ andL error bounds are obtained for a continuous piecewise linear finite element approximation of the volume matching problem. This problem consists of minimising |v| _1, ² overvH ¹() subject to the inequality constraintv0 and a number of linear equality constraints. The presence of the equality constraints leads to Lagrange multipliers, which in turn lead to complications with the standard error analysis for variational inequalities. Finally we consider an algorithm for solving the resulting algebraic problem.Supported by a SERC research studentship     

Finite element approach to the Stokes problem

We approximate the Stokes problem by using a finite element method. This method utilizes the approach of Kleiser–Schumann, in which a boundary condition for the pressure is implicitly defined by a condition on the velocity. We consider a suitable splitting of the unknowns that allows one to reduce the Stokes problem to a cascade of classical Dirichlet problems and to a boundary integral equation.     

Finite element analysis of pressure vessel using beam on elastic foundation analysis

In this study, we first applied the variation principle to derive a new finite element method (FEM) based on the theory of beam on elastic foundation using line element. The derived FEM was then applied to solve, for the first time, the pressure vessel problems with uniform thickness. Our FEM results, obtained even by using only one line element, agreed exactly with the available closed-form solution, confirming the validity and computing efficiency of our finite element formulation. Moreover, we have applied our new FEM to solve pressure vessel problems with non-uniform thickness where no exact analytical solution is known to exist. The distributions of discontinuity stress in the cylindrical part were obtained. We found that shear force and bending moment were indeed discontinuous at the geometrically discontinuous juncture, due to the bending rigidity and elastic constant change by the non-uniform thickness. Finally, the case of discontinuity stresses in a bimetallic joint was also studied. The locations of maximum shear force and bending moment were found to be affected by the bending rigidity of the material.     

Finite element analysis of a coupling eigenvalue problem on overlapping domains

In this paper, we consider a nonstandard elliptic eigenvalue problem on a rectangular domain, consisting of two overlapping rectangles, where the interaction between the subdomains is expressed through an integral coupling condition on their intersection. For this problem we set up finite element (FE) approximations, without and with numerical quadrature. The involved error analysis is affected by the nonlocal coupling condition, which requires the introduction and error estimation of a suitably modified vector Lagrange interpolant on the overall FE mesh. As a consequence, the resulting error estimates are sub-optimal, as compared to the ones established, e.g., in Vanmaele and van Keer (RAIRO – Math. Mod. Num. Anal 29(3) (1995) 339–365) for classical eigenvalue problems with local boundary or transition conditions.     

Finite element analysis of a static contact problem with Coulomb friction

A unilateral contact problem with a variable coefficient of friction is solved by a simplest variant of the finite element technique. The coefficient of friction may depend on the magnitude of the tangential displacement. The existence of an approximate solution and some a priori estimates are proved.     

Finite element analysis of a current density-electric field formulation of Bean's model for superconductivity

We study a current density–electric field formulationof Bean's model for the experimental set-up of an infinitelylong cylindrical superconductor subject to a transverse magneticfield. We introduce a fully practical finite-element approximationof the model and prove an error between the exact solution andthe approximate solution for the current density of order (h+ t)^1/2. Numerical simulations for a variety of given sourcecurrents are presented.     

Finite element approximation of nonlocal parabolic problem

In this article, we present a finite element scheme combined with backward Euler method to solve a nonlocal parabolic problem. An important issue in the numerical solution of nonlocal problems while using Newton's method is related to its structure. In fact differently from the local case where the Jacobian matrix is sparse and banded, in the nonlocal case the Jacobian matrix is dense and computations are much more onerous compared to that for differential equations. In order to avoid this difficulty, we use the technique given by Gudi (SIAM J Numer Anal 50 (2012), 657–668) for elliptic nonlocal problem of Kirchhoff type. We discuss the well‐posedness of the weak formulation at continuous as well as at discrete levels. We also derive a priori error estimates for semidiscrete and fully discrete formulations in L² and H¹ norms. Results based on the usual finite element method are provided to confirm the theoretical estimates. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 786–813, 2017     

Finite element analysis of the vibration problem of a plate coupled with a fluid

Summary. We consider the approximation of the vibration modes of an elastic plate in contact with a compressible fluid. The plate is modelled by Reissner-Mindlin equations while the fluid is described in terms of displacement variables. This formulation leads to a symmetric eigenvalue problem. Reissner-Mindlin equations are discretized by a mixed method, the equations for the fluid with Raviart-Thomas elements and a non conforming coupling is used on the interface. In order to prove that the method is locking free we consider a family of problems, one for each thickness , and introduce appropriate scalings for the physical parameters so that these problems attain a limit when . We prove that spurious eigenvalues do not arise with this discretization and we obtain optimal order error estimates for the eigenvalues and eigenvectors valid uniformly on the thickness parameter t. Finally we present numerical results confirming the good performance of the method. Received February 4, 1998 / Revised version received May 26, 1999 / Published online June 21, 2000     

A combined nodal continuous-discontinuous finite element formulation for the Maxwell problem

Continuous Galerkin formulations are appealing due to their low computational cost, whereas discontinuous Galerkin formulation facilitate adaptative mesh refinement and are more accurate in regions with jumps of physical parameters. Since many electromagnetic problems involve materials with different physical properties, this last point is very important. For this reason, in this article we have developed a combined cG-dG formulation for Maxwell’s problem that allows arbitrary finite element spaces with functions continuous in patches of finite elements and discontinuous on the interfaces of these patches. In particular, the second formulation we propose comes from a novel continuous Galerkin formulation that reduces the amount of stabilization introduced in the numerical system. In all cases, we have performed stability and convergence analyses of the methods. The outcome of this work is a new approach that keeps the low CPU cost of recent nodal continuous formulations with the ability to deal with coefficient jumps and adaptivity of discontinuous ones. All these methods have been tested using a problem with singular solution and another one with different materials, in order to prove that in fact the resulting formulations can properly deal with these problems.     

Finite element error analysis of a zeroth order approximate deconvolution model based on a mixed formulation

A suitable discretization for the Zeroth Order Model in Large Eddy Simulation of turbulent flows is sought. This is a low order model, but its importance lies in the insight that it provides for the analysis of the higher order models actually used in practice by the pioneers Stolz and Adams [N.A. Adams, S. Stolz, On the approximate deconvolution procedure for LES, Phys. Fluids 2 (1999) 1699-1701; N.A. Adams, S. Stolz, Deconvolution methods for subgrid-scale approximation in large eddy simulation, in: B.J. Geurts (Ed.), Modern Simul. Strategies for Turbulent Flow, Edwards, Philadelphia, 2001, pp. 21-44] and others. The higher order models have proven to be of high accuracy. However, stable discretizations of them have proven to be tricky and other stabilizations, such as time relaxation and eddy viscosity, are often added. We propose a discretization based on a mixed variational formulation that gives the correct energy balance. We show it to be unconditionally stable and prove convergence.     

Finite element approximation to nonlinear coupled thermal problem

A nonlinear coupled elliptic system modelling a large class of engineering problems was discussed in [A.F.D. Loula, J. Zhu, Finite element analysis of a coupled nonlinear system, Comp. Appl. Math. 20 (3) (2001) 321–339; J. Zhu, A.F.D. Loula, Mixed finite element analysis of a thermally nonlinear coupled problem, Numer. Methods Partial Differential Equations 22 (1) (2006) 180–196]. The convergence analysis of iterative finite element approximation to the solution was done under an assumption of ‘small’ solution or source data which guarantees the uniqueness of the nonlinear coupled system. Generally, a nonlinear system may have multiple solutions. In this work, the regularity of the weak solutions is further studied. The nonlinear finite element approximations to the nonsingular solutions are then proposed and analyzed. Finally, the optimal order error estimates in H¹$H^1$-norm and L²$L^2$-norm as well as in W^1,p$W^{1,p}$-norm and L^p$L^p$-norm are obtained.     

Optimal preconditioner for boundary element formulation of the Dirichlet problem in elasticity

An optimal preconditioning procedure for the numerical solution of two-dimensional Dirichlet problem for Lamé equations by boundary element method is constructed. An efficient algorithm for the above problem is also developed.     

Finite element analysis of varitional crimes for a quasilinear elliptic problem in 3D

Summary. We examine a finite element approximation of a quasilinear boundary value elliptic problem in a three-dimensional bounded convex domain with a smooth boundary. The domain is approximated by a polyhedron and a numerical integration is taken into account. We apply linear tetrahedral finite elements and prove the convergence of approximate solutions on polyhedral domains in the -norm to the true solution without any additional regularity assumptions. Received May 23, 1997 / Published online December 6, 1999     

Finite element method for the diffraction problem in a waveguide

The convergence of the finite element method (FEM) as applied to the diffraction problem in a waveguide in the case when there is no damping in the medium filling the waveguide is proved. A functional space that takes into account the partial radiation conditions is introduced to carry out the proof. A highly accurate approximation method for the partial radiation conditions is considered.     

Finite element approximation of a free boundary plasma problem

In this article, we study a finite element approximation for a model free boundary plasma problem. Using a mixed approach (which resembles an optimal control problem with control constraints), we formulate a weak formulation and study the existence and uniqueness of a solution to the continuous model problem. Using the same setting, we formulate and analyze the discrete problem. We derive optimal order energy norm a priori error estimates proving the convergence of the method. Further, we derive a reliable and efficient a posteriori error estimator for the adaptive mesh refinement algorithm. Finally, we illustrate the theoretical results by some numerical examples.