Researcher profile

Jean-François Remacle

Jean-François Remacle contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
5topics
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

7 published item(s)

preprint2021arXiv

Multidirectionnal sweeping preconditioners with non-overlapping checkerboard domain decomposition for Helmholtz problems

This paper explores a family of generalized sweeping preconditionners for Helmholtz problems with non-overlapping checkerboard partition of the computational domain. The domain decomposition procedure relies on high-order transmission conditions and cross-point treatments, which cannot scale without an efficient preconditioning technique when the number of subdomains increases. With the proposed approach, existing sweeping preconditioners, such as the symmetric Gauss-Seidel and parallel double sweep preconditioners, can be applied to checkerboard partitions with different sweeping directions (e.g. horizontal and diagonal). Several directions can be combined thanks to the flexible version of GMRES, allowing for the rapid transfer of information in the different zones of the computational domain, then accelerating the convergence of the final iterative solution procedure. Several two-dimensional finite element results are proposed to study and to compare the sweeping preconditioners, and to illustrate the performance on cases of increasing complexity.

preprint2021arXiv

Quasi-structured quadrilateral meshing in Gmsh -- a robust pipeline for complex CAD models

We propose an end-to-end pipeline to robustly generate high-quality quadrilateral meshes for complex CAD models. An initial quad-dominant mesh is generated with frontal point insertion guided by a locally integrable cross field and a scalar size map adapted to the small CAD features. After triangle combination and midpoint-subdivision into an all-quadrilateral mesh, the topology of the mesh is modified to reduce the number of irregular vertices. The idea is to preserve the irregular vertices matching cross-field singularities and to eliminate the others. The topological modifications are either local and based on disk quadrangulations, or more global with the remeshing of patches of quads according to predefined patterns. Validity of the quad mesh is guaranteed by monitoring element quality during all operations and reverting the changes when necessary. Advantages of our approach include robustness, strict respect of the CAD features and support for user-prescribed size constraints. The quad mesher, which is available in Gmsh, is validated and illustrated on two datasets of CAD models.

preprint2020arXiv

Automatic feature-preserving size field for 3D mesh generation

This paper presents a methodology aiming at easing considerably the generation of high-quality meshes for complex 3D domains. We show that the whole mesh generation process can be controlled with only five parameters to generate in one stroke quality meshes for arbitrary geometries. The main idea is to build a meshsize field $h(x)$ taking local features of the geometry, such as curvatures, into account. Meshsize information is then propagated from the surfaces into the volume, ensuring that the magnitude of $\vert \nabla h \vert$ is always controlled so as to obtain a smoothly graded mesh. As the meshsize field is stored in an independent octree data structure, the function h can be computed separately, and then plugged in into any mesh generator able to respect a prescribed meshsize field. The whole procedure is automatic, in the sense that minimal interaction with the user is required. Applications examples based on models taken from the very large ABC dataset, are then presented, all treated with the same generic set of parameter values, to demonstrate the efficiency and the universality of the technique.

preprint2020arXiv

Automatic surface mesh generation for discrete models: A complete and automatic pipeline based on reparameterization

Triangulations are an ubiquitous input for the finite element community. However, most raw triangulations obtained by imaging techniques are unsuitable as-is for finite element analysis. In this paper, we give a robust pipeline for handling those triangulations, based on the computation of a one-to-one parametrization for automatically selected patches of input triangles, which makes each patch amenable to remeshing by standard finite element meshing algorithms. Using only geometrical arguments, we prove that a discrete parametrization of a patch is one-to-one if (and only if) its image in the parameter space is such that all parametric triangles have a positive area. We then derive a non-standard linear discretization scheme based on mean value coordinates to compute such one-to-one parametrizations, and show that the scheme does not discretize a Laplacian on a structured mesh. The proposed pipeline is implemented in the open source mesh generator Gmsh, where the creation of suitable patches is based on triangulation topology and parametrization quality, combined with feature edge detection. Several examples illustrate the robustness of the resulting implementation.

preprint2020arXiv

Computing cross fields -- A PDE approach based on the Ginzburg-Landau theory

This paper proposes a method to compute crossfields based on the Ginzburg-Landau theory. The Ginzburg-Landau functional has two terms: the Dirichlet energy of the distribution and a term penalizing the mismatch between the fixed and actual norm of the distribution. Directional fields on surfaces are known to have a number of critical points, which are properly identified with the Ginzburg-Landau approach: the asymptotic behavior of Ginzburg-Landau problem provides well-distributed critical points over the 2-manifold, whose indices are as low as possible. The central idea in this paper is to exploit this theoretical background for crossfield computation on arbitrary surfaces. Such crossfields are instrumental in the generation of meshes with quadrangular elements. The relation between the topological properties of quadrangular meshes and crossfields are hence first recalled. It is then shown that a crossfield on a surface can be represented by a complex function of unit norm with a number of critical points, i.e., a nearly everywhere smooth function taking its values in the unit circle of the complex plane. As maximal smoothness of the crossfield is equivalent with minimal energy, the crossfield problem is equivalent to an optimization problem based on Ginzburg-Landau functional. A discretization scheme with Crouzeix-Raviart elements is applied and the correctness of the resulting finite element formulation is validated on the unit disk by comparison with an analytical solution. The method is also applied to the 2-sphere where, surprisingly but rightly, the computed critical points are not located at the vertices of a cube, but at those of an anticube.

preprint2020arXiv

Quality tetrahedral mesh generation with HXT

We proposed, in a recent paper (10.1002/nme.5987), a fast 3D parallel Delaunay kernel for tetrahedral mesh generation. This kernel was however incomplete in the sense that it lacked the necessary mesh improvement tools. The present paper builds on that previous work and proposes a fast parallel mesh improvement stage that delivers high-quality tetrahedral meshes compared to alternative open-source mesh generators. Our mesh improvement toolkit includes edge removal and improved Laplacian smoothing as well as a brand new operator called the Growing SPR Cavity, which can be regarded as the mother of all flips. The paper describes the workflow of the new mesh improvement schedule, as well as the details of the implementation. The result of this research is a series of open-source scalable software components, called HXT, whose overall efficiency is demonstrated on practical examples by means of a detailed comparative benchmark with two open-source mesh generators: Gmsh and TetGen.

preprint2018arXiv

Representing three-dimensional cross fields using 4th order tensors

This paper presents a new way of describing cross fields based on fourth order tensors. We prove that the new formulation is forming a linear space in $\mathbb{R}^9$. The algebraic structure of the tensors and their projections on $\mbox{SO}(3)$ are presented. The relationship of the new formulation with spherical harmonics is exposed. This paper is quite theoretical. Due to pages limitation, few practical aspects related to the computations of cross fields are exposed. Nevetheless, a global smoothing algorithm is briefly presented and computation of cross fields are finally depicted.