BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Lennard Kamenski

Lennard Kamenski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.NANumerical Analysis

Catalog footprint

What is connected

10works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

10 published item(s)

preprint2019arXiv

Conditioning of implicit Runge-Kutta integration for finite element approximation of linear diffusion equations on anisotropic meshes

The conditioning of implicit Runge-Kutta (RK) integration for linear finite element approximation of diffusion equations on general anisotropic meshes is investigated. Bounds are established for the condition number of the resulting linear system with and without diagonal preconditioning for the implicit Euler and general implicit RK methods. Two solution strategies are considered for the linear system resulting from general implicit RK integration: the simultaneous solution (the system is solved as a whole) and a successive solution which follows the commonly used implementation of implicit RK methods to first transform the system into smaller systems using the Jordan normal form of the RK matrix and then solve them successively. For the simultaneous solution in case of a positive semidefinite symmetric part of the RK coefficient matrix and for the successive solution it is shown that . If the smallest eigenvalue of the symmetric part of the RK coefficient matrix is negative and the simultaneous solution strategy is used, an upper bound on the time step is given so that the system matrix is positive definite. The obtained bounds for the condition number have explicit geometric interpretations and take the interplay between the diffusion matrix and the mesh geometry into full consideration. They show that there are three mesh-dependent factors that can affect the conditioning: the number of elements, the mesh nonuniformity measured in the Euclidean metric, and the mesh nonuniformity with respect to the inverse of the diffusion matrix. They also reveal that the preconditioning using the diagonal of the system matrix, the mass matrix, or the lumped mass matrix can effectively eliminate the effects of the mesh nonuniformity measured in the Euclidean metric. Numerical examples are given.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Why do we need Voronoi cells and Delaunay meshes? Essential properties of the Voronoi finite volume method

Unlike other schemes that locally violate the essential stability properties of the analytic parabolic and elliptic problems, Voronoi finite volume methods (FVM) and boundary conforming Delaunay meshes provide good approximation of the geometry of a problem and are able to preserve the essential qualitative properties of the solution for any given resolution in space and time as well as changes in time scales of multiple orders of magnitude. This work provides a brief description of the essential and useful properties of the Voronoi FVM, application examples, and a motivation why Voronoi FVM deserve to be used more often in practice than they are currently.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stability of explicit one-step methods for P1-finite element approximation of linear diffusion equations on anisotropic meshes

We study the stability of explicit one-step integration schemes for the linear finite element approximation of linear parabolic equations. The derived bound on the largest permissible time step is tight for any mesh and any diffusion matrix within a factor of $2(d+1)$, where $d$ is the spatial dimension. Both full mass matrix and mass lumping are considered. The bound reveals that the stability condition is affected by two factors. The first one depends on the number of mesh elements and corresponds to the classic bound for the Laplace operator on a uniform mesh. The other factor reflects the effects of the interplay of the mesh geometry and the diffusion matrix. It is shown that it is not the mesh geometry itself but the mesh geometry in relation to the diffusion matrix that is crucial to the stability of explicit methods. When the mesh is uniform in the metric specified by the inverse of the diffusion matrix, the stability condition is comparable to the situation with the Laplace operator on a uniform mesh. Numerical results are presented to verify the theoretical findings.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A comparative numerical study of meshing functionals for variational mesh adaptation

We present a comparative numerical study for three functionals used for variational mesh adaptation. One of them is a generalisation of Winslow's variable diffusion functional while the others are based on equidistribution and alignment. These functionals are known to have nice theoretical properties and work well for most mesh adaptation problems either as a stand-alone variational method or combined within the moving mesh framework. Their performance is investigated numerically in terms of equidistribution and alignment mesh quality measures. Numerical results in 2D and 3D are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A geometric discretization and a simple implementation for variational mesh generation and adaptation

We present a simple direct discretization for functionals used in the variational mesh generation and adaptation. Meshing functionals are discretized on simplicial meshes and the Jacobian matrix of the continuous coordinate transformation is approximated by the Jacobian matrices of affine mappings between elements. The advantage of this direct geometric discretization is that it preserves the basic geometric structure of the continuous functional, which is useful in preventing strong decoupling or loss of integral constraints satisfied by the functional. Moreover, the discretized functional is a function of the coordinates of mesh vertices and its derivatives have a simple analytical form, which allows a simple implementation of variational mesh generation and adaptation on computer. Since the variational mesh adaptation is the base for a number of adaptive moving mesh and mesh smoothing methods, the result in this work can be used to develop simple implementations of those methods. Numerical examples are given.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

How a nonconvergent recovered Hessian works in mesh adaptation

Hessian recovery has been commonly used in mesh adaptation for obtaining the required magnitude and direction information of the solution error. Unfortunately, a recovered Hessian from a linear finite element approximation is nonconvergent in general as the mesh is refined. It has been observed numerically that adaptive meshes based on such a nonconvergent recovered Hessian can nevertheless lead to an optimal error in the finite element approximation. This also explains why Hessian recovery is still widely used despite its nonconvergence. In this paper we develop an error bound for the linear finite element solution of a general boundary value problem under a mild assumption on the closeness of the recovered Hessian to the exact one. Numerical results show that this closeness assumption is satisfied by the recovered Hessian obtained with commonly used Hessian recovery methods. Moreover, it is shown that the finite element error changes gradually with the closeness of the recovered Hessian. This provides an explanation on how a nonconvergent recovered Hessian works in mesh adaptation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Adaptive finite elements with anisotropic meshes

The paper presents a numerical study for the finite element method with anisotropic meshes. We compare the accuracy of the numerical solutions on quasi-uniform, isotropic, and anisotropic meshes for a test problem which combines several difficulties of a corner singularity, a peak, a boundary layer, and a wavefront. Numerical experiment clearly shows the advantage of anisotropic mesh adaptation. The conditioning of the resulting linear equation system is addressed as well. In particular, it is shown that the conditioning with adaptive anisotropic meshes is not as bad as generally assumed.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Conditioning of Finite Element Equations with Arbitrary Anisotropic Meshes

Bounds are developed for the condition number of the linear finite element equations of an anisotropic diffusion problem with arbitrary meshes. They depend on three factors. The first, factor proportional to a power of the number of mesh elements, represents the condition number of the linear finite element equations for the Laplacian operator on a uniform mesh. The other two factors arise from the mesh nonuniformity viewed in the Euclidean metric and in the metric defined by the diffusion matrix. The new bounds reveal that the conditioning of the finite element equations with adaptive anisotropic meshes is much better than what is commonly feared. Diagonal scaling for the linear system and its effects on the conditioning are also studied. It is shown that the Jacobi preconditioning, which is an optimal diagonal scaling for a symmetric positive definite sparse matrix, can eliminate the effects of mesh nonuniformity viewed in the Euclidean metric and reduce those effects of the mesh viewed in the metric defined by the diffusion matrix. Tight bounds on the extreme eigenvalues of the stiffness and mass matrices are obtained. Numerical examples are given.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Anisotropic Mesh Adaptation for Variational Problems Using Error Estimation Based on Hierarchical Bases

Anisotropic mesh adaptation has been successfully applied to the numerical solution of partial differential equations but little considered for variational problems. In this paper, we investigate the use of a global hierarchical basis error estimator for the development of an anisotropic metric tensor needed for the adaptive finite element solution of variational problems. The new metric tensor is completely a~posteriori and based on residual, edge jumps and the hierarchical basis error estimator. Numerical results show that it performs comparable with existing metric tensors based on Hessian recovery. A few sweeps of the symmetric Gauß-Seidel iteration for solving the global error problem prove sufficient to provide directional information necessary for successful mesh adaptation. .

0 citations0 reviewsNo code linkedOpen access