Cut Finite Elements for Convection in Fractured Domains

We develop a cut finite element method (CutFEM) for the convection problem in a so called fractured domain which is a union of manifolds of different dimensions such that a $d$ dimensional component always resides on the boundary of a $d+1$ dimensional component. This type of domain can for instance be used to model porous media with embedded fractures that may intersect. The convection problem can be formulated in a compact form suitable for analysis using natural abstract directional derivative and divergence operators. The cut finite element method is based on using a fixed background mesh that covers the domain and the manifolds are allowed to cut through a fixed background mesh in an arbitrary way. We consider a simple method based on continuous piecewise linear elements together with weak enforcement of the coupling conditions and stabilization. We prove a priori error estimates and present illustrating numerical examples

Design sensitivity analysis using EAL. Part 1: Conventional design parameters

A numerical implementation of design sensitivity analysis of builtup structures is presented, using the versatility and convenience of an existing finite element structural analysis code and its database management system. The finite element code used in the implemenatation presented is the Engineering Analysis Language (EAL), which is based on a hybrid method of analysis. It was shown that design sensitivity computations can be carried out using the database management system of EAL, without writing a separate program and a separate database. Conventional (sizing) design parameters such as cross-sectional area of beams or thickness of plates and plane elastic solid components are considered. Compliance, displacement, and stress functionals are considered as performance criteria. The method presented is being extended to implement shape design sensitivity analysis using a domain method and a design component method

Real time plasma equilibrium reconstruction in a Tokamak

The problem of equilibrium of a plasma in a Tokamak is a free boundary problemdescribed by the Grad-Shafranov equation in axisymmetric configurations. The right hand side of this equation is a non linear source, which represents the toroidal component of the plasma current density. This paper deals with the real time identification of this non linear source from experimental measurements. The proposed method is based on a fixed point algorithm, a finite element resolution, a reduced basis method and a least-square optimization formulation

Exhibiting cross-diffusion-induced patterns for reaction-diffusion systems on evolving domains and surfaces

The aim of this manuscript is to present for the first time the application of the finite element method for solving reaction-diffusion systems with cross-diffusion on continuously evolving domains and surfaces. Furthermore we present pattern formation generated by the reaction-diffusion systemwith cross-diffusion on evolving domains and surfaces. A two-component reaction-diffusion system with linear cross-diffusion in both u and v is presented. The finite element method is based on the approximation of the domain or surface by a triangulated domain or surface consisting of a union of triangles. For surfaces, the vertices of the triangulation lie on the continuous surface. A finite element space of functions is then defined by taking the continuous functions which are linear affine on each simplex of the triangulated domain or surface. To demonstrate the role of cross-diffusion to the theory of pattern formation, we compute patterns with model kinetic parameter values that belong only to the cross-diffusion parameter space; these do not belong to the standard parameter space for classical reaction-diffusion systems. Numerical results exhibited show the robustness, flexibility, versatility, and generality of our methodology; the methodology can deal with complicated evolution laws of the domain and surface, and these include uniform isotropic and anisotropic growth profiles as well as those profiles driven by chemical concentrations residing in the domain or on the surface

Solving the incompressible surface Navier-Stokes equation by surface finite elements

We consider a numerical approach for the incompressible surface Navier-Stokes equation on surfaces with arbitrary genus $g(\mathcal{S})$. The approach is based on a reformulation of the equation in Cartesian coordinates of the embedding $\mathbb{R}^3$, penalization of the normal component, a Chorin projection method and discretization in space by surface finite elements for each component. The approach thus requires only standard ingredients which most finite element implementations can offer. We compare computational results with discrete exterior calculus (DEC) simulations on a torus and demonstrate the interplay of the flow field with the topology by showing realizations of the Poincar\'e-Hopf theorem on $n$-tori

A fully implicit multi-axial solution strategy for direct ratchet boundary evaluation : theoretical development

Ensuring sufficient safety against ratchet is a fundamental requirement in pressure vessel design. Determining the ratchet boundary can prove difficult and computationally expensive when using a full elastic-plastic finite element analysis and a number of direct methods have been proposed that overcome the difficulties associated with ratchet boundary evaluation. Here, a new approach based on fully implicit Finite Element methods, similar to conventional elastic-plastic methods, is presented. The method utilizes a two-stage procedure. The first stage determines the cyclic stress state, which can include a varying residual stress component, by repeatedly converging on the solution for the different loads by superposition of elastic stress solutions using a modified elastic-plastic solution. The second stage calculates the constant loads which can be added to the steady cycle whilst ensuring the equivalent stresses remain below a modified yield strength. During stage 2 the modified yield strength is updated throughout the analysis, thus satisfying Melan’s Lower bound ratchet theorem. This is achieved utilizing the same elastic plastic model as the first stage, and a modified radial return method. The proposed methods are shown to provide better agreement with upper bound ratchet methods than other lower bound ratchet methods, however limitations in these are identified and discussed

Design sensitivity analysis using EAL. Part 2: Shape design parameters

A numerical implementation of shape design sensitivity analysis of built-up structures is presented, using the versatility and convenience of an existing finite element structural analysis code and its data base management system. This report is a continuation of a previous report on conventional design parameters. The finite element code used in the implementation presented is the Engineering Analysis Language (EAL), which is based on a hybrid analysis method. It has been shown that shape design sensitivity computations can be carried out using the database management system of EAL, without writing a separate program and a separate data base. The material derivative concept of continuum mechanics and an adjoint variable method of design sensitivity analysis are used to derive shape design sensitivity information of structural performances. A domain method of shape design sensitivity analysis and a design component method are used. Displacement and stress functionals are considered as performance criteria

A Study on Using Hierarchical Basis Error Estimates in Anisotropic Mesh Adaptation for the Finite Element Method

A common approach for generating an anisotropic mesh is the M-uniform mesh approach where an adaptive mesh is generated as a uniform one in the metric specified by a given tensor M. A key component is the determination of an appropriate metric which is often based on some type of Hessian recovery. Recently, the use of a global hierarchical basis error estimator was proposed for the development of an anisotropic metric tensor for the adaptive finite element solution. This study discusses the use of this method for a selection of different applications. Numerical results show that the method performs well and is comparable with existing metric tensors based on Hessian recovery. Also, it can provide even better adaptation to the solution if applied to problems with gradient jumps and steep boundary layers. For the Poisson problem in a domain with a corner singularity, the new method provides meshes that are fully comparable to the theoretically optimal meshes.Comment: The final publication is available at http://www.springerlink.co