Metric field construction for anisotropic mesh adaptation with application to blood flow simulations.

The goal of this paper is to generate an anisotropic metric field suitable for cardiovascular geometries before a fluid simulation. Starting from a curvature map, an initial surface metric field is computed. This metric is used for anisotropic surface mesh adaptation and consecutively extended inside the volume in a frontal manner. The algorithm is based on the method proposed by Alauzet but replaces the metric intersection steps by an original metric 'blending'. This allows for a graded anisotropic volume mesh with a refinement layer close to the walls. The benefits of the resulting mesh are multiple: a reduced number of degrees of freedom, a priori refinement in areas with strong gradients of velocity and automatically increased resolution in regions with high surface curvature. The primal application of this method is in the domain of cardiovascular flows. Flow fields and derived quantities (wall shear stress) through a model bypass around a stenosed artery obtained on an adapted and standard isotropic mesh are compared. In addition, the mesh generation method is tested on a more complex patient-specific geometry. Values of computed wall shear stress are shown to be close to values obtained on anisotropic Hessian-adapted mesh, demonstrating the computational efficiency of the approach in comparison with adaptation based on error indicators derived from flow field.

Finite element analysis of thermal laser skin stimulation for a finer characterization of the nociceptive system.

Thermal laser stimulation of the skin is an efficient exploratory tool to characterize the nociceptive system. In the present study, finite element simulations are done to calculate the intra-cutaneous spatio-temporal temperature profiles following the delivery of such laser stimuli. The proposed computer-aided modeling considers a number of important parameters that have been disregarded in previous approaches: (i) variability of water content across the skin in both hairy and glabrous skin, (ii) temperature dependency of optical and thermal skin parameters, (iii) laser wavelength and corresponding absorption coefficient, (iv) beam shape (Gaussian vs. flat top) and (v) power emission (closed vs. open loop). Numerical simulations allow determining at each instant of time the volume and area of skin tissue whose temperature exceeds a given nociceptor activation threshold. This knowledge allows a finer characterization of the subpopulations of primary afferents that encode and convey nociceptive signals to the central nervous system. As an example, the approach is used to obtain an estimate of intraepidermal nerve fiber density in both physiological and pathological conditions. Moreover, a better knowledge of the heat distribution also reduces the risk of injury to the skin. Finally, in order to make the finite element simulations accessible to investigators with no prior background in numerical analysis, a specific open-source user-interface has been developed with the ONELAB software.

Cardiovascular and lung mesh generation based on centerlines.

We present a fully automatic procedure for the mesh generation of tubular geometries such as blood vessels or airways. The procedure is implemented in the open-source Gmsh software and relies on a centerline description of the input geometry. The presented method can generate different type of meshes: isotropic tetrahedral meshes, anisotropic tetrahedral meshes, and mixed hexahedral/tetrahedral meshes. Additionally, a multiple layered arterial wall can be generated with a variable thickness. All the generated meshes rely on a mesh size field and a mesh metric that is based on centerline descriptions, namely the distance to the centerlines and a local reference system based on the tangent and the normal directions to the centerlines. Different examples show that the proposed method is very efficient and robust and leads to high quality computational meshes.