Researcher profile

Martin Ladecký

Martin Ladecký contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 19 - UnverifiedVerification L1Unclaimed author
5works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

5 published item(s)

preprint2022arXiv

Non-convex, ringing-free, FFT-accelerated solver using an incremental approximate energy functional

Fourier-accelerated micromechanical homogenization has been developed and applied to a variety of problems, despite being prone to ringing artifacts. In addition, the majority of Fourier-accelerated solvers applied to FFT-accelerated schemes only apply to convex problems. We here introduce a that allows to employ modern efficient and non-convex iterative solvers, such as trust-region solvers or LBFGS in a FFT-accelerated scheme. These solvers need the explicit energy functional of the system in their standard form. We develop a modified trust region solver, capable of handling non-convex micromechanical homogenization problems such as continuum damage employing the approximate incremental energy functional. We use the developed solver as the solver of a ringing-free FFT-accelerated solution scheme, namely the projection based scheme with finite element discretization.

preprint2022arXiv

Optimal FFT-accelerated Finite Element Solver for Homogenization

We propose a matrix-free finite element (FE) homogenization scheme that is considerably more efficient than generic FE implementations. The efficiency of our scheme follows from a preconditioned well-scaled reformulation allowing for the use of the conjugate gradient or similar iterative solvers. The geometrically-optimal preconditioner -- a discretized Green's function of a periodic homogeneous reference problem -- has a block-diagonal structure in the Fourier space which permits its efficient inversion using the fast Fourier transform (FFT) techniques for generic regular meshes. This implies that the scheme scales as $\mathcal{O}(n \log(n))$ like FFT, rendering it equivalent to spectral solvers in terms of computational efficiency. However, in contrast to classical spectral solvers, the proposed scheme works with FE shape functions with local supports and is free of the Fourier ringing phenomenon. We showcase that the scheme achieves the number of iterations that are almost independent of spatial discretisation and scales mildly with the phase contrast. Additionally, we discuss the equivalence between our displacement-based scheme and the recently proposed strain-based homogenization technique with finite-element projection.

preprint2021arXiv

Elimination of ringing artifacts by finite-element projection in FFT-based homogenization

Micromechanical homogenization is often carried out with Fourier-accelerated methods that are prone to ringing artifacts. We here generalize the compatibility projection introduced by Vondřejc, Zeman & Marek [Comput. Math. Appl. 68, 156 (2014)] beyond the Fourier basis. In particular, we formulate the compatibility projection for linear finite elements while maintaining Fourier-acceleration and the fast convergence properties of the original method. We demonstrate that this eliminates ringing artifacts and yields an efficient computational homogenization scheme that is equivalent to canonical finite-element formulations on fully structured grids.

preprint2020arXiv

FFT-based homogenisation accelerated by low-rank tensor approximations

Fast Fourier transform (FFT) based methods have turned out to be an effective computational approach for numerical homogenisation. In particular, Fourier-Galerkin methods are computational methods for partial differential equations that are discretised with trigonometric polynomials. Their computational effectiveness benefits from efficient FFT based algorithms as well as a favourable condition number. Here these kind of methods are accelerated by low-rank tensor approximation techniques for a solution field using canonical polyadic, Tucker, and tensor train formats. This reduced order model also allows to efficiently compute suboptimal global basis functions without solving the full problem. It significantly reduces computational and memory requirements for problems with a material coefficient field that admits a moderate rank approximation. The advantages of this approach against those using full material tensors are demonstrated using numerical examples for the model homogenisation problem that consists of a scalar linear elliptic variational problem defined in two and three dimensional settings with continuous and discontinuous heterogeneous material coefficients. This approach opens up the potential of an efficient reduced order modelling of large scale engineering problems with heterogeneous material.

preprint2020arXiv

Guaranteed two-sided bounds on all eigenvalues of preconditioned diffusion and elasticity problems solved by the finite element method

A method of estimating all eigenvalues of a preconditioned discretized scalar diffusion operator with Dirichlet boundary conditions has been recently introduced in T. Gergelits, K.A. Mardal, B.F. Nielsen, Z. Strakoš: Laplacian preconditioning of elliptic PDEs: Localization of the eigenvalues of the discretized operator, SIAM Journal on Numerical Analysis 57(3) (2019), 1369-1394. Motivated by this paper, we offer a slightly different approach that extends the previous results in some directions. Namely, we provide bounds on all increasingly ordered eigenvalues of a general diffusion or elasticity operator with tensor data, discretized with the conforming finite element method, preconditioned by the inverse of a matrix of the same operator with different data. Our results hold for mixed Dirichlet and Robin or periodic boundary conditions applied to the original and preconditioning problems. The bounds are two-sided, guaranteed, easily accessible, and depend solely on the material data.