Combining Boundary-Conforming Finite Element Meshes on Moving Domains Using a Sliding Mesh Approach

For most finite element simulations, boundary-conforming meshes have significant advantages in terms of accuracy or efficiency. This is particularly true for complex domains. However, with increased complexity of the domain, generating a boundary-conforming mesh becomes more difficult and time consuming. One might therefore decide to resort to an approach where individual boundary-conforming meshes are pieced together in a modular fashion to form a larger domain. This paper presents a stabilized finite element formulation for fluid and temperature equations on sliding meshes. It couples the solution fields of multiple subdomains whose boundaries slide along each other on common interfaces. Thus, the method allows to use highly tuned boundary-conforming meshes for each subdomain that are only coupled at the overlapping boundary interfaces. In contrast to standard overlapping or fictitious domain methods the coupling is broken down to few interfaces with reduced geometric dimension. The formulation consists of the following key ingredients: the coupling of the solution fields on the overlapping surfaces is imposed weakly using a stabilized version of Nitsche's method. It ensures mass and energy conservation at the common interfaces. Additionally, we allow to impose weak Dirichlet boundary conditions at the non-overlapping parts of the interfaces. We present a detailed numerical study for the resulting stabilized formulation. It shows optimal convergence behavior for both Newtonian and generalized Newtonian material models. Simulations of flow of plastic melt inside single-screw as well as twin-screw extruders demonstrate the applicability of the method to complex and relevant industrial applications

Solving multi-physics problems using adaptive finite elements with independently refined meshes

In this thesis, we study a numerical tool named multi-mesh method within the framework of the adaptive finite element method. The aim of this method is to minimize the size of the linear system to get the optimal performance of simulations. Multi-mesh methods are typically used in multi-physics problems, where more than one component is involved in the system. During the discretization of the weak formulation of partial differential equations, a finite-dimensional space associated with an independently refined mesh is assigned to each component respectively. The usage of independently refined meshes leads less degrees of freedom from a global point of view. To our best knowledge, the first multi-mesh method was presented at the beginning of the 21st Century. Similar techniques were announced by different mathematics researchers afterwards. But, due to some common restrictions, this method is not widely used in the field of numerical simulations. On one hand, only the case of two-mesh is taken into scientists\' consideration. But more than two components are common in multi-physics problems. Each is, in principle, allowed to be defined on an independent mesh. Besides that, the multi-mesh methods presented so far omit the possibility that coefficient function spaces live on the different meshes from the trial and test function spaces. As a ubiquitous numerical tool, the multi-mesh method should comprise the above circumstances. On the other hand, users are accustomed to improving the performance by taking the advantage of parallel resources rather than running simulations with the multi-mesh approach on one single processor, so it would be a pity if such an efficient method was only available in sequential. The multi-mesh method is actually used within local assembling process, which should not be conflict with parallelization. In this thesis, we present a general multi-mesh method without the limitation of the number of meshes used in the system, and it can be applied to parallel environments as well. Chapter 1 introduces the background knowledge of the adaptive finite element method and the pioneering work, on which this thesis is based. Then, the main idea of the multi-mesh method is formally derived and the detailed implementation is discussed in Chapter 2 and 3. In Chapter 4, applications, e.g. the multi-phase flow problem and the dendritic growth, are shown to prove that our method is superior in contrast to the standard single-mesh finite element method in terms of performance, while accuracy is not reduced

Proceedings of the FEniCS Conference 2017

Proceedings of the FEniCS Conference 2017 that took place 12-14 June 2017 at the University of Luxembourg, Luxembourg

Concurrent multiscale analysis without meshing: Microscale representation with CutFEM and micro/macro model blending

In this paper, we develop a novel unfitted multiscale framework that combines two separate scales represented by only one single computational mesh. Our framework relies on a mixed zooming technique where we zoom at regions of interest to capture microscale properties and then mix the micro and macroscale properties in a transition region. Furthermore, we use homogenization techniques to derive macro model material properties. The microscale features are discretized using CutFEM. The transition region between the micro and macroscale is represented by a smooth blending function. To address the issues with ill-conditioning of the multiscale system matrix due to the arbitrary intersections in cut elements and the transition region, we add stabilization terms acting on the jumps of the normal gradient (ghost-penalty stabilization). We show that our multiscale framework is stable and is capable to reproduce mechanical responses for heterogeneous structures in a mesh-independent manner. The efficiency of our methodology is exemplified by 2D and 3D numerical simulations of linear elasticity problems

Aeronautical engineering: A continuing bibliography with indexes (supplement 227)

This bibliography lists 418 reports, articles, and other documents introduced into the NASA scientific and technical information system in May, 1988

Anisotropic Adaptivity and Subgrid Scale Modelling for the Solution of the Neutron Transport Equation with an Emphasis on Shielding Applications

This thesis demonstrates advanced new discretisation and adaptive meshing technologies that improve the accuracy and stability of using finite element discretisations applied to the Boltzmann transport equation (BTE). This equation describes the advective transport of neutral particles such as neutrons and photons within a domain. The BTE is difficult to solve, due to its large phase space (three dimensions of space, two of angle and one each of energy and time) and the presence of non-physical oscillations in many situations. This work explores the use of a finite element method that combines the advantages of the two schemes: the discontinuous and continuous Galerkin methods. The new discretisation uses multiscale (subgrid) finite elements that work locally within each element in the finite element mesh in addition to a global, continuous, formulation. The use of higher order functions that describe the variation of the angular flux over each element is also explored using these subgrid finite element schemes. In addition to the spatial discretisation, methods have also been developed to optimise the finite element mesh in order to reduce resulting errors in the solution over the domain, or locally in situations where there is a goal of specific interest (such as a dose in a detector region). The chapters of this thesis have been structured to be submitted individually for journal publication, and are arranged as follows. Chapter 1 introduces the reader to motivation behind the research contained within this thesis. Chapter 2 introduces the forms of the BTE that are used within this thesis. Chapter 3 provides the methods that are used, together with examples, of the validation and verification of the software that was developed as a result of this work, the transport code RADIANT. Chapter 4 introduces the inner element subgrid scale finite element discretisation of the BTE that forms the basis of the discretisations within RADIANT and explores its convergence and computational times on a set of benchmark problems. Chapter 5 develops the error metrics that are used to optimise the mesh in order to reduce the discretisation error within a finite element mesh using anisotropic adaptivity that can use elongated elements that accurately resolves computational demanding regions, such as in the presence of shocks. The work of this chapter is then extended in Chapter 6 that forms error metrics for goal based adaptivity to minimise the error in a detector response. Finally, conclusions from this thesis and suggestions for future work that may be explored are discussed in Chapter 7.Open Acces

A geometry preserving, conservative, mesh-to-mesh isogeometric interpolation algorithm for spatial adaptivity of the multigroup, second-order even-parity form of the neutron transport equation

In this paper a method is presented for the application of energy-dependent spatial meshes applied to the multigroup, second-order, even-parity form of the neutron transport equation using Isogeometric Analysis (IGA). The computation of the inter-group regenerative source terms is based on conservative interpolation by Galerkin projection. The use of Non-Uniform Rational B-splines (NURBS) from the original computer-aided design (CAD) model allows for efficient implementation and calculation of the spatial projection operations while avoiding the complications of matching different geometric approximations faced by traditional finite element methods (FEM). The rate-of-convergence was verified using the method of manufactured solutions (MMS) and found to preserve the theoretical rates when interpolating between spatial meshes of different refinements. The scheme’s numerical efficiency was then studied using a series of two-energy group pincell test cases where a significant saving in the number of degrees-of-freedom can be found if the energy group with a complex variation in the solution is refined more than an energy group with a simpler solution function. Finally, the method was applied to a heterogeneous, seven-group reactor pincell where the spatial meshes for each energy group were adaptively selected for refinement. It was observed that by refining selected energy groups a reduction in the total number of degrees-of-freedom for the same total L2 error can be obtained

Cut Finite Element Methods on Overlapping Meshes: Analysis and Applications

This thesis deals with both analysis and applications of cut finite element methods (CutFEMs) on overlapping meshes. By overlapping meshes we mean a mesh hierarchy with a background mesh at the bottom and one or more overlapping meshes that are stacked on top of it. Overlapping meshes can be used as an alternative to costly remeshing for problems with changing geometry. The main content of the thesis is the five appended papers. The thesis consists of an analysis part and an applications part.In the analysis part (Paper I and Paper II), we consider cut finite element methods on overlapping meshes for a time-dependent\ua0parabolic model problem: the heat equation on two overlapping meshes, where one mesh is allowed to move around on top of the other. In Paper I, the overlapping mesh is prescribed a cG(1) movement, meaning that its location as a function of time is continuous and piecewise linear. The cG(1) mesh movement results in a space-time discretization for which existing analysis methodologies either fail or are unsuitable. We therefore propose, to the best of our knowledge, a new energy analysis framework that is general enough to be applicable to the current setting. In Paper II, the overlapping mesh is prescribed a dG(0) movement, meaning that its location as a function of time is\ua0discontinuous\ua0and\ua0piecewise constant. The dG(0) mesh movement results in a space-time discretization for which existing analysis methodologies work with some modifications to handle the shift in the overlapping mesh\u27s location at discrete times.The applications part (Paper III, IV, and V) concerns cut finite element methods on overlapping meshes for\ua0stationary\ua0PDE-problems. We consider two potential applications for CutFEM on overlapping meshes. The first application, presented in Paper III, presents methodology for evaluating configurations of buildings based on wind and view. The wind model is based on a CutFEM on overlapping meshes for Stokes equations. The second application, presented in Paper IV and Paper V, concerns a software application (app). The app lets a user define and solve physical problems governed by PDEs in an immersive and interactive augmented reality environment