Local refinement techniques for elliptic problems on cell-centered grids. III. Algebraic multilevel BEPS preconditioners.
Mesh r-adaptation for unilateral contact problems
We present a mesh adaptation method by node movement for two-dimensional linear elasticity problems with unilateral contact. The adaptation is based on a hierarchical estimator on finite element edges and the node displacement techniques use an analogy of the mesh topology with a spring network. We show, through numerical examples, the efficiency of the present adaptation method.
Mesh Refinement For Stabilized Convection Diffusion Equations
We derive a residual a posteriori error estimates for the subscales stabilization of convection diffusion equation. The estimator yields upper bound on the error which is global and lower bound that is local
Mixed-hybrid model of the fracture flow
Moving mesh for the axisymmetric harmonic map flow
We build corotational symmetric solutions to the harmonic map flow from the unit disc into the unit sphere which have constant degree. First, we prove the existence of such solutions, using a time semi-discrete scheme based on the idea that the harmonic map flow is the -gradient of the relaxed Dirichlet energy. We prove a partial uniqueness result concerning these solutions. Then, we compute numerically these solutions by a moving-mesh method which allows us to deal with the singularity at the...
Moving mesh for the axisymmetric harmonic map flow
We build corotational symmetric solutions to the harmonic map flow from the unit disc into the unit sphere which have constant degree. First, we prove the existence of such solutions, using a time semi-discrete scheme based on the idea that the harmonic map flow is the L2-gradient of the relaxed Dirichlet energy. We prove a partial uniqueness result concerning these solutions. Then, we compute numerically these solutions by a moving-mesh method which allows us to deal with the singularity at the...
Multilevel iterative methods for mixed finite element discretizations of elliptic problems.
Nonobtuse tetrahedral partitions that refine locally towards Fichera-like corners
Linear tetrahedral finite elements whose dihedral angles are all nonobtuse guarantee the validity of the discrete maximum principle for a wide class of second order elliptic and parabolic problems. In this paper we present an algorithm which generates nonobtuse face-to-face tetrahedral partitions that refine locally towards a given Fichera-like corner of a particular polyhedral domain.
Norm estimates for a maximal right inverse of the divergence operator in spaces of piecewise polynomials
Numerical integration for high order pyramidal finite elements
We examine the effect of numerical integration on the accuracy of high order conforming pyramidal finite element methods. Non-smooth shape functions are indispensable to the construction of pyramidal elements, and this means the conventional treatment of numerical integration, which requires that the finite element approximation space is piecewise polynomial, cannot be applied. We develop an analysis that allows the finite element approximation space to include non-smooth functions and show that,...
Numerical integration for high order pyramidal finite elements
We examine the effect of numerical integration on the accuracy of high order conforming pyramidal finite element methods. Non-smooth shape functions are indispensable to the construction of pyramidal elements, and this means the conventional treatment of numerical integration, which requires that the finite element approximation space is piecewise polynomial, cannot be applied. We develop an analysis that allows the finite element approximation space to include non-smooth functions and show that,...
Numerical solution of a scattering problem in a wave guide using a grid generation technique.
Object oriented design philosophy for scientific computing
This contribution gives an overview of current research in applying object oriented programming to scientific computing at the computational mechanics laboratory (LABMEC) at the school of civil engineering – UNICAMP. The main goal of applying object oriented programming to scientific computing is to implement increasingly complex algorithms in a structured manner and to hide the complexity behind a simple user interface. The following areas are current topics of research and documented within the...
Object oriented design philosophy for scientific computing
This contribution gives an overview of current research in applying object oriented programming to scientific computing at the computational mechanics laboratory (LABMEC) at the school of civil engineering – UNICAMP. The main goal of applying object oriented programming to scientific computing is to implement increasingly complex algorithms in a structured manner and to hide the complexity behind a simple user interface. The following areas are current topics of research and documented within the...
Observations Regarding Algorithms Required for Robust CFD Codes
Over the last three decades Computational Fluid Dynamics (CFD) has gradually joined the wind tunnel and flight test as a primary flow analysis tool for aerodynamic designers. CFD has had its most favorable impact on the aerodynamic design of the high-speed cruise configuration of a transport. This success has raised expectations among aerodynamicists that the applicability of CFD can be extended to the full flight envelope. However, the complex nature...
On a posteriori error estimators in the infinite element method on anisotropic meshes.
On FE-grid relocation in solving unilateral boundary value problems by FEM
We consider FE-grid optimization in elliptic unilateral boundary value problems. The criterion used in grid optimization is the total potential energy of the system. It is shown that minimization of this cost functional means a decrease of the discretization error or a better approximation of the unilateral boundary conditions. Design sensitivity analysis is given with respect to the movement of nodal points. Numerical results for the Dirichlet-Signorini problem for the Laplace equation and the...
On interpolation error on degenerating prismatic elements
We propose an analogue of the maximum angle condition (commonly used in finite element analysis for triangular and tetrahedral meshes) for the case of prismatic elements. Under this condition, prisms in the meshes may degenerate in certain ways, violating the so-called inscribed ball condition presented by P. G. Ciarlet (1978), but the interpolation error remains of the order in the -norm for sufficiently smooth functions.
On quadrature methods for the double layer potential equation over the boundary of a polyhedron.
On Synge-type angle condition for -simplices
The maximum angle condition of J. L. Synge was originally introduced in interpolation theory and further used in finite element analysis and applications for triangular and later also for tetrahedral finite element meshes. In this paper we present some of its generalizations to higher-dimensional simplicial elements. In particular, we prove optimal interpolation properties of linear simplicial elements in that degenerate in some way.