Application of a Hierarchical Chromosome Based Genetic Algorithm to the Problem of Finding Optimal Initial Meshes for the Self-Adaptive hp-FEM

The paper presents an algorithm for finding the optimal initial mesh for the self-adaptive hp Finite Element Method (hp-FEM) calculations. We propose the application of the hierarchical chromosome based genetic algorithm for optimal selection of the initial mesh. The selection of the optimal initial mesh will optimize the convergence rate of the numerical error of the solution over the sequence of meshes generated by the self-adaptive hp-FEM. This is especially true in the case when material data are selected as a result of some stochastic algorithm and it is not possible to design optimal initial mesh by hand. The algorithm has been tested on the non-stationary mass transport problem modeling phase transition phenomenon

Hypergraph Grammars in hp-adaptive Finite Element Method

AbstractThe paper presents the hypergraph grammar for modelling the hp-adaptive finite element method algorithm with rectangular elements. The finite element mesh is represented by a hypergraph. All mesh transformations are modelled by means of hypergraph grammar rules. These rules allow to generate the initial mesh, to assign values of polynomial order to the element nodes, to generate the matrix for each element, to solve the problem and to perform the hp-adaptation

Hypergrammar-based parallel multi-frontal solver for grids with point singularities

This paper describes the application of hypergraph grammars to drive linear computationalcost solver for grids with point singularities. Such graph grammar productions are the rstmathematical formalism used to describe solver algorithm and each of them indicates thesmallest atomic task that can be executed in parallel, which is very useful in case of parallelexecution. In particular the partial order of execution of graph grammar productions can befound, and the sets of independent graph grammar productions can be localized. They canbe scheduled set by set into shared memory parallel machine. The graph grammar basedsolver has been implemented with NIVIDIA CUDA for GPU. Graph grammar productionsare accompanied by numerical results for 2D case. We show that our graph grammar basedsolver with GPU accelerator is order of magnitude faster than state of the art MUMPSsolver

Petri Nets Modeling of Dead-End Refinement Problems in a 3D Anisotropic hp-Adaptive Finite Element Method

We consider two graph grammar based Petri nets models for anisotropic refinements of three dimensional hexahedral grids. The first one detects possible dead-end problems during the graph grammar based anisotropic refinements of the mesh. The second one employs an enhanced graph grammar model that is actually dead-end free. We apply the resulting algorithm to the simulation of resistivity logging measurements for estimating the location of underground oil and/or gas formations. The graph grammar based Petri net models allow to fix the self-adaptive mesh refinement algorithm and finish the adaptive computations with the required accuracy needed by the numerical solution

Applications of a hyper-graph grammar system in adaptive finite-element computations

This paper describes application of a hyper-graph grammar system for modeling a three-dimensional adaptive finite element method. The hyper-graph grammar approach allows obtaining a linear computational cost of adaptive mesh transformations and computations performed over refined meshes. The computations are done by a hyper-graph grammar driven algorithm applicable to three-dimensional problems. For the case of typical refinements performed towards a point or an edge, the algorithm yields linear computational cost with respect to the mesh nodes for its sequential execution and logarithmic cost for its parallel execution. Such hyper-graph grammar productions are the mathematical formalism used to describe the computational algorithm implementing the finite element method. Each production indicates the smallest atomic task that can be executed concurrently. The mesh transformations and computations by using the hyper-graph grammar-based approach have been tested in the GALOIS environment. We conclude the paper with some numerical results performed on a shared-memory Linux cluster node, for the case of three-dimensional computational meshes refined towards a point, an edge and a face

Non-carious tooth loss in terms of erosion – a literature review

The loss of tooth structure caused by the process of chemical dissolution and mechanical abrasion is becoming a great problem in modern dentistry. Patient awareness as well as diagnostic alertness of a dentist are crucial for proper prevention and treatment. Late diagnosis and treatment of extensive loss of tooth structure caused by erosion, abrasion or abfraction poses many difficulties for clinicians. Chronically progressive loss of tooth structure can cause loss of vertical dimension of the bite, secondary orthodontic complications and also aesthetic and functional problems. The treatment is often associated with long-term diagnosis and interdisciplinary, expensive rehabilitation of the entire stomatognathic system. The aim of this study was to review current publications on the prevalence of tooth wear in developmental age, with focus on etiology, risk factors, prevention and treatment methods

Petri Nets Modeling of Dead-End Refinement Problems in a 3D Anisotropic hp-Adaptive Finite Element Method

We consider two graph grammar based Petri nets models for anisotropic refinements of three dimensional hexahedral grids. The first one detects possible dead-end problems during the graph grammar based anisotropic refinements of the mesh. The second one employs an enhanced graph grammar model that is actually dead-end free. We apply the resulting algorithm to the simulation of resistivity logging measurements for estimating the location of underground oil and/or gas formations. The graph grammar based Petri net models allow to fix the self-adaptive mesh refinement algorithm and finish the adaptive computations with the required accuracy needed by the numerical solution

Graph-grammar based algorithm for asteroid tsunami simulations

On January 18, 2022, around 1 million kilometers from Earth, five times the distance from Earth to the Moon, a large asteroid passed without harm to the Earth. Theoretically, however, the event of the asteroid falling into Earth, causing the tsunami, is possible since there are over 27,000 near-Earth asteroids [1], and the Earth’s surface is covered in 71 percent by water. We introduce a novel graph-grammar-based framework for asteroid tsunami simulations. Our framework adaptively generates the computational mesh of the Earth model. It is built from triangular elements representing the seashore and the seabed. The computational mesh is represented as a graph, with graph vertices representing the computational mesh element’s interiors and edges. Mesh refinements are often performed by the longest-edge refinement algorithm. We have expressed this algorithm by only two graph-grammar productions. The resulting graph represents the terrain approximating the topography with a prescribed accuracy. We generalize the graph-grammar mesh refinement algorithm to work on the entire Earth model, allowing the generation of the terrain topography, including the seabed. Having the seashore and the seabed represented by a graph, we introduce the finite element method simulations of the tsunami wave propagation. We illustrate the framework with simulations of the disastrous asteroid falling into the Baltic sea

Automatic stabilization of finite-element simulations using neural networks and hierarchical matrices

Petrov–Galerkin formulations with optimal test functions allow for the stabilization of finite element simulations. In particular, given a discrete trial space, the optimal test space induces a numerical scheme delivering the best approximation in terms of a problem-dependent energy norm. This ideal approach has two shortcomings: first, we need to explicitly know the set of optimal test functions; and second, the optimal test functions may have large supports inducing expensive dense linear systems. A concise proposal on how to overcome these shortcomings has been raised during the last decade by the Discontinuous Petrov–Galerkin (DPG) methodology. However, DPG has also some limitations and difficulties: the method requires ultraweak variational formulations, obtained through a hybridization process, which is not trivial to implement at the discrete level. Our motivation is to offer a simpler alternative for the case of parametric PDEs, which can be used with any variational formulation. Indeed, parametric families of PDEs are an example where it is worth investing some (offline) computational effort to obtain stabilized linear systems that can be solved efficiently in an online stage, for a given range of parameters. Therefore, as a remedy for the first shortcoming, we explicitly compute (offline) a function mapping any PDE parameter, to the matrix of coefficients of optimal test functions (in some basis expansion) associated with that PDE parameter. Next, as a remedy for the second shortcoming, we use the low-rank approximation to hierarchically compress the (non-square) matrix of coefficients of optimal test functions. In order to accelerate this process, we train a neural network to learn a critical bottleneck of the compression algorithm (for a given set of PDE parameters). When solving online the resulting (compressed) Petrov–Galerkin formulation, we employ a GMRES iterative solver with inexpensive matrix–vector multiplications thanks to the low-rank features of the compressed matrix. We perform experiments showing that the full online procedure is as fast as an (unstable) Galerkin approach. We illustrate our findings by means of 2D–3D Eriksson–Johnson problems, together with 2D Helmholtz equation