Abstract:Highlights • An efficient high-order immersed finite element method for the automated analysis of CAD models is proposed. • Discretisation relies on three distinct structured grids for implicit geometry, finite elements, and quadrature. • All three grids are adaptive and distributed for parallel computation. • Several auxiliary techniques for high-order, robust treatment of boundaries are presented. • Analysis of CAD models with tens of millions of finite element DOFs within a few minutes is demonstrated. The automated finite element analysis of complex CAD models using boundary-fitted meshes is rife with difficulties. Immersed finite element methods are intrinsically more robust but usually less accurate. In this work, we introduce an efficient, robust, high-order immersed finite element method for complex CAD models. Our approach relies on three adaptive structured grids: a geometry grid for representing the implicit geometry, a finite element grid for discretising physical fields and a quadrature grid for evaluating the finite element integrals. The geometry grid is a sparse VDB (Volumetric Dynamic B+ tree) grid that is highly refined close to physical domain boundaries. The finite element grid consists of a forest of octree grids distributed over several processors, and the quadrature grid in each finite element cell is an octree grid constructed in a bottom-up fashion. The resolution of the quadrature grid ensures that finite element integrals are evaluated with sufficient accuracy and that any sub-grid geometric features, like small holes or corners, are resolved up to a desired resolution. The conceptual simplicity and modularity of our approach make it possible to reuse open-source libraries, i.e. openVDB and p4est for implementing the geometry and finite element grids, respectively, and BDDCML for iteratively solving the discrete systems of equations in parallel using domain decomposition. We demonstrate the efficiency and robustness of the proposed approach by solving the Poisson equation on domains described by complex CAD models and discretised with tens of millions of degrees of freedom. The solution field is discretised using linear and quadratic Lagrange basis functions.

Software and Algorithmic Aspects of Automating Finite-element Discretization

Automatic Generation and Fast Implement of Finite Element Mesh

Automation of Numerical Simulation on One-way Fluid-structure Interaction of Fluid-filled Pipeline

Discrete Model Based Analysis Algorithm for Surface Smoothness

Automating the Finite Element Method

Analysis of structures with arbitrary kinematic boundary conditions by the semi-analytical finite element method

A Three-Grid High-Order Immersed Finite Element Method for the Analysis of CAD Models

IMPLEMENTATION OF A FINITE ELEMENT CLASS LIBRARY USING GENERALIZED PROGRAMMING

Analysis of the application of various numerical integration methods and their practical implementation in MATLAB CAD

An Optimization Model and Algorithm for the Smoothing of Sizing Functions of Three-Dimensional Finite Element Meshes

High order unfitted finite element discretizations for explicit boundary representations

Mathematical and programming toolkit of the computer aided design of the axonometric piping diagrams

CURRENT TRENDS IN THE DEVELOPMENT OF AUTOMATED COMPUTER MODELING SYSTEMS

Integrating CAD and Numerical Analysis: 'Dirty Geometry' handling using the Finite Cell Method

Software Engineering Development of Finite Element Method Programming Applications in 2D Frame Structures Using Python Programs

DOLFIN: Automated Finite Element Computing

Development of a Software Package of Smoothed Finite Element Method (S-FEM) for Solid Mechanics Problems

From CAD to representations suitable for isogeometric analysis: a complete pipeline

An unstructured immersed finite element method for nonlinear solid mechanics

Finite Volume Discretizations for Eigenvalue Problems with Applications to Electronic Structure Calculations.