key: cord-0045933-n2jz4yzv authors: Geronzi, Leonardo; Gasparotti, Emanuele; Capellini, Katia; Cella, Ubaldo; Groth, Corrado; Porziani, Stefano; Chiappa, Andrea; Celi, Simona; Biancolini, Marco Evangelos title: Advanced Radial Basis Functions Mesh Morphing for High Fidelity Fluid-Structure Interaction with Known Movement of the Walls: Simulation of an Aortic Valve date: 2020-05-25 journal: Computational Science - ICCS 2020 DOI: 10.1007/978-3-030-50433-5_22 sha: 1bebb0d96018c1e1e273dcd2e52a3c08f773b820 doc_id: 45933 cord_uid: n2jz4yzv High fidelity Fluid-Structure Interaction (FSI) can be tackled by means of non-linear Finite Element Models (FEM) suitable to capture large deflections of structural parts interacting with fluids and by means of detailed Computational Fluid Dynamics (CFD). High fidelity is gained thanks to the spatial resolution of the computational grids and a key enabler to have a proper exchange of information between the structural solver and the fluid one is the management of the interfaces. A class of applications consists in problems where the complex movement of the walls is known in advance or can be computed by FEM and has to be transferred to the CFD solver. The aforementioned approach, known also as one-way FSI, requires effective methods for the time marching adaption of the computation grid of the CFD model. A versatile and well established approach consists in a continuum update of the mesh that is regenerated so to fit the evolution of the moving walls. In this study, an innovative method based on Radial Basis Functions (RBF) mesh morphing is proposed, allowing to keep the same mesh topology suitable for a continuum update of the shape. A set of key configurations are exactly guaranteed whilst time interpolation is adopted between frames. The new framework is detailed and then demonstrated, adopting as a reference the established approach based on remeshing, for the study of a Polymeric-Prosthetic Heart Valve (P-PHV). In all the engineering areas of development, multiphysics analysis appears highly difficult to be carried out because of the interactions between more than one physics involved. During the achievement and the analysis of coupled systems, the way in which the shape of the object influences its performances is required to be carefully taken into account. In the biomedical engineering field for example, the design and the evaluation of the behaviour of prosthetic valves, stent-grafts and ventricular assist devices are related to both the structural and fluid mechanics physics [1] [2] [3] . Generally speaking, numerical meshes need to be created for the specific kind of analysis and, in a multi-physics context, two or more meshes have to be realised. In this environment, each geometric change is applied to all the numerical models involved in the analysis: such update has to be performed in a rapid way and as easy as possible. This task, usually carried out through remeshing methods, may be also obtained faster with the use of mesh morphing techniques. This approach allows the changing of the shape of a meshed surface so that the topology is preserved while nodal positions are updated [4, 5] : modifications are applied on a baseline grid by moving the surface nodes and propagating displacements inside the surrounding volume mesh nodes. Concerning biomedical applications, morphing methods have been applied in both bone and cardiovascular fields. Recently, in [6] , an interactive sculpting and RBF mesh morphing approach has been proposed to address geometry modifications using a force-feedback device, while in [7] RBF have been applied to improve cranioplasty applications. In the cardiovascular field, morphing approaches were employed in [8] for the registration procedures of the cardiac muscle and to model an aorta aneurysm carrying a one-way FSI [9] [10] [11] . In literature RBF have been extensively employed to tackle FSI problems [12] , using the modal method for both static [13, 14] and transient simulations [15] [16] [17] , as well as the two-way approach [18, 19] with RBF-based mapping methods [20] . In this work the state of the art regarding biomedical one-way FSI applications, as shown in [11] , is furtherly improved for a transient simulation by taking into account the non-linear deformations of the wetted surfaces during motion. An ad-hoc workflow was developed to transfer and to update, using an RBF-based morphing technique, the CFD mesh, incrementally adjusting the geometry according to the non-linear evolution predicted by the FEM solver. To demonstrate the effectiveness of the proposed approach, it was applied to a tailor-made Polymeric-Prosthetic Heart Valve; these devices [21] proved to significantly reduce blood coagulation problems, maintaining excellent properties in term of strength, efficient function and long-term durability [22, 23] . This work is arranged as follows: at first, an introduction on FSI coupling is given, comparing remeshing and morphing workflows. The proposed procedure is implemented in the following paragraph, in which the incremental approach is shown on the P-PHV case. Results are finally discussed and compared to those obtained by remeshing. One of the most delicate processes in a Fluid-Structure Interaction analysis concerns the management of the fluid-solid interfaces; at the boundary surfaces between these two domains, solution data is shared between the fluid solver and the structural one. A useful and employed kind of FSI analysis concerns all the simulations in which a moving body interacts with an internal or external surrounding fluid in a unidirectional way. This one-way coupling does not guarantee energy conservation at the fluid-solid interface but holds the benefit of lower computational time in comparison to the bidirectional one, in which a continuous data exchange between Computational Structural Mechanics (CSM) and CFD solvers is required. In this work two Arbitrary Lagrangian Eulerian (ALE) methods [24] [25] [26] for moving meshes are employed: the first based on remeshing algorithms, called in this paper "standard" approach and the novel one using RBF mesh morphing procedures. In this first kind of analysis, as shown in Fig. 1 , a specific component is responsible to transfer the deformation of the CSM mesh to the grid of the CFD solver. In fact, one of the most applied strategies to manage body meshes in FSI applications is simply to move the grid for as long as possible and, when the quality of the mesh becomes critical, to enable the remeshing tools updating the lowquality cells generated by high displacements in the fluid domain [27] . Larger are the displacements of the grid, wider are the distortions of the cells; remeshing method agglomerates cells that violate the initially defined Skewness [28] or size criteria and remeshes them. If the new cells or faces satisfy the Skewness criterion, the mesh is locally updated with the new cells, interpolating the solution from the old cells. Otherwise, the new cells are discarded and another remeshing step is required. Obviously, when small displacements are involved and just few remeshing steps are necessary, this approach may be considered efficient. However, for larger displacements, the number of remeshing steps increases to avoid the presence of negative and invalid elements; as a result, the simulation is slowed down because a new mesh with different numbers and positions of nodes and cells has to be generated more frequently. To overcome the problems related to remeshing, a procedure based on mesh morphing to adapt the shape according to a target one [29] is here described. Among the morphing methods available in literature, RBF are well known for their interpolation quality also on very large meshes [4, 30] . RBF allow to interpolate everywhere in the space a scalar functions known at discrete points, called Source points (Sp). By interpolating three scalar values it is possible, solving a linear system of order equal to the number of Sp employed [4] , to describe a displacement of the Sp in the three directions in space. The interpolation function is defined as follows: where x is a generic position in the space, x si the Sp position, s(·) the scalar function which represents a transformation R n → R, ϕ (·) the radial function of order m, γ i the weight and h(x) a polynomial term with degree m − 1 added to improve the fit assuring uniqueness of the problem and polynomial precision. The unknowns of the system, namely the polynomial coefficients and the weights γ i of the radial functions, are retrieved by imposing the passage of the function on the given values and an orthogonality condition on the polinomials. If the RBF is conditionally positive definite, it can be demonstrated that a unique interpolant exists and in 3D, if the order is equal or less than 2, a linear polynomial in the form h(x) = β 1 + β 2 x + β 3 y + β 4 z can be used. The linear problem can be also written in matrix form: in which M is the interpolation matrix containing all the distances between RBF centres M ij = ϕ ( x i − x j ), P the matrix containing the polynomial terms that has for each row j the form P j = [1 x 1j x 2j ... x nj ] and g the known values at Sp. The new nodal positions, if interpolating the displacements, can be retrieved for each node as: Being the problem solved pointwise, the approach is meshless and able to manage every kind element (tetrahedral, hexahedral, polyhedral and others), both for surface and volume mesh smoothing ensuring the preservation of their topology. In general, a morphing operation can introduce a reduction of the mesh quality but a good morpher has to minimize this effect and to maximize the possible shape modifications. If mesh quality is well preserved, morphing has a clear benefit over remeshing because it avoids introducing noise [4] . This procedure To apply the shape deformation without remeshing, the morpher RBF Morph is in this study employed [31] . To delimit the morphing action, an encapsulation technique has to be implemented, defining in this way sub-domains or parts of the fluid domain within which the morpher action is applied [32] . According to [33], S p effectively used are a random subgroup which includes a percentage close to 2% of the total extracted surface nodes of the structural model. A Radial Basis Function is selected and, when the settings of the procedure are completed, morphing solutions are calculated, each per morphing transitions, from the starting configuration to the final position of the model. Every transition, a new deformed mesh is generated and the solution files are stored. CFD Set-Up. Once all the morphed meshes are obtained, the CFD analysis is carried out every time step on a new deformed mesh, by using a specific script in Scheme Programming Language. This script is recalled before the starting of the CFD simulation and it allows the automatization of the MultiSol tool of RBF Morph software in order to guarantee the synchronous implementation of the deformed shape in fluid environment. The nodal displacement inside two subsequent morphed positions is simply handled inside the Scheme script by means of a parameter called amplification A(t) which varies over time from 0 to A 0 (usually equal to +1 or −1). Because of the samplings available at different times, an accurate and effective solution can be obviously achieved, handling all the morphing deformations separately and linearly modulating the amplification within each of the sampled grids: calling i the index that concerns all the mesh conformations P 0 ...P n , the amplification A i (t) inside each time step can be evaluated as: In this way, the fluid grid is able to follow the displacements of the CSM mesh and solution can be quickly calculated every time on a new deformed grid. In this paper, both the presented FSI solution approaches are applied to study the opening phase of a P-PHV. This kind of analysis is conducted on a prosthetic valve to study its haemodynamic behaviour during the opening phase in which the leaflets or cusps [34, 35] are pushed open to allow the ejection of the blood flow and to investigate the effects of the solid domain (P-PHV) towards the fluid domain (blood flow). The design of a prosthetic aortic valve is realised using the software Space-Claim, following the geometry proposed in [36] , and then imported inside ANSYS Workbench (v193). The workflow of the first approach based on remeshing and implemented inside ANSYS Workbench is proposed in Fig. 4 . In this case, the displacements of a Concerning the loads, a physiological time-varying pressure is uniformly applied to the ventricular portion of the valve to simulate the opening pressure. Fluent was used as CFD solver, activating the Dynamic Mesh to deform and regenerate the computational grid [38]: the Spring-Laplace based Smoothing Method and the Remeshing Method, two tools to solve the analysis of deformed domains due to boundary movement over time [39] , are applied. Inside the System Coupling component, all the surfaces of the valve in contact with the blood are set as fluid-solid interfaces and a time step of 1e-5 s is selected. The novel workflow, based on a RBF mesh morphing, also in this case developed into the ANSYS Workbench environment is shown in Fig. 5 which depicts the overall incremental procedure; in particular, the System Coupling component, fundamental in the previously shown approach, here is not employed. Results obtained with the use of RBF Morph software are compared to those achieved using Fluent remeshing tools. They are reported in terms of pressure map, velocity streamlines and computational time. Structural Simulations. The same CSM analysis implemented in the remeshing approach is conducted. A specific Application Customization Toolkit (ACT) Extension, i.e. a method to achieve custom applications inside ANSYS Mechanical, is developed to collect Sp and Tp according to a sampling time chosen by the user: in this case, ten different surface nodal positions (P 0 ...P 9 ) are collected. The RBF Morph Fluent (v193) Add-On is adopted: the RBF are used to step by step project the S p to the T p and interpolate the displacements of the volume nodes inside the blood domain over time. RBF Morph tools may impose a blood domain adaptation in the CFD analysis following the deformed shapes of the valvular domain, extracted from the previous CSM simulation. In the case of the valve, to control the morphing action in the space, no encapsulation volumes [32] are generated but a fixed Surf Set (with null displacement during the movement of the leaflets) on the wall of the Sinotubular Junction (STJ) [40] is preferred (Fig. 6) . In this way, since the distortion of the mesh is extremely high, such deformations can be distributed also in the other zones of the fluid domain and not only in proximity to the valve. To prevent mesh quality reduction, in this work the bi-harmonic kernel ϕ(r) = r is used: it is the RBF that guarantees to keep the mesh quality degradation at minimum [41] . The nine solutions are calculated using the Solve Panel of the software, each per morphing action. Two different strategies are evaluated: -Forwards: from P 0 (initial position) directly to P 9 (totally open position). -Backwards/Forwards (Fig. 7) : extracting a new already deformed valve in position P 2 of which a new mesh is achieved, from P 2 going back to P 0 (backwards -step 1) and then from P 0 moving on to P 9 (forwards -step 2). The first strategy produced an excessive mesh distortion with negative cells volume in particular at the centre of the leaflets. No mesh problems are detected with the adopted second approach. Fluent Setting. The same time step of the standard approach (1e-5 s) is selected. The flow is initialized with the valve in position P 0 and the Scheme file with the amplification of each of the nine different transitions is recalled inside Fluent. Thanks to this workflow, the FSI simulation can be conducted only by means of this CFD analysis in which the mesh is updated every time step. The first important parameter to check in the implementation of both the approaches is the Skewness: no negative cells are detected and maximum Sk values are 0.858 and 0.961 respectively for the procedure implemented through remeshing and morphing. No convergence issues due to mesh degradation are recorded. In this work, only pressure and velocity results of the 3D-model at a specific and significant time step are shown using the CFD-Post software. Pressure Map -A 3D-pressure rendering at the time in which the valve is for the first time open (7 ms) is represented in Fig. 8 ; as it can be observed, no significant differences in terms of pressure values are recorded between the two approaches: the mean pressure are 69.64 mmHg and 70.31 mmHg respectively for the simulation based on remeshing and the analysis using mesh morphing (difference of less than 1%). Velocity Streamlines -Streamlines of both the solution methods are reported in Fig. 9 . It is possible to observe how in this 3D-representation the local maximum values are placed in the same zones. Maximum local fluid velocities are reached right at t = 7 ms because the fast movement of the leaflets pushes away the blood: maximum value for the FSI implemented with the novel approach is 3.011 m/s, about 5.58% less than the value detected in the analysis based on remeshing, that is 3.189 m/s. Computational Time -The time required to calculate the solution is the most interesting difference between these two solution methods: with the same number of elements and to parity of time step (1e-5 s), the standard workflow requires 6283 min to solve the problem while the CFD analysis with the morphing application only 396 min, approximately 16 times faster. Furthermore, if time step is increased, Fluent solver using remeshing starts having troubles and the obvious consequence is the no-convergence of the model. This phenomenon does not happen with morphing by which time step can be highly increased without problems and solution time can be reduced up to 60 times compared to that of remeshing. In this work, a fast high fidelity workflow for a multi-physic Fluid-Structure Interaction analysis exploiting ANSYS Workbench and RBF Morph is presented. In particular, this powerful system has proved to be useful and accurate in all those applications in which the movement of the fluid domain is strongly controlled by the deformation of the solid domain. The possibility to replicate the non-linear analysis in which the motion of a solid model is transferred to a fluid domain carefully handling the interfaces is demonstrated. In a specific application case, the attempt to transfer aortic valve opening kinematics to the CFD analysis through a mesh morphing technique results to be successful and the effectiveness of the FSI simulation implemented through mesh morphing to reproduce the fluid dynamic field of the aortic valve is here shown and validated. Morphing has turned out to be very consistent in comparison to remeshing. Further research is needed to asses in more detail the differences between these two distinct approaches. Testing the novel workflow on a FSI analysis about the study of the opening of a polymeric aortic valve, it has allowed to reduce the simulation time up to 16 times in comparison to that required by remeshing methods using the same time step. Until now, numerical FSI simulations have been used in a restricted manner precisely because of the long time required to solve them, as in the case of heart valve computational analyses. Dealing with clinical trials, this approach could ensure a considerable time saving with a decrease of the final mesh quality which results acceptable for the numerical solver here adopted. This new workflow has definitely broken down this limit, granting in this case the possibility to test several models of P-PHV in a lot of flow conditions. In addition, a parametric design of P-PHV combined with this method could finally make a strong contribution to the patient-specific aortic valve replacement. Numerical methods for fluid-structure interaction models of aortic valves A literature review of the numerical analysis of abdominal aortic aneurysms treated with endovascular stent grafts Numerical model of flow in a sac-type ventricular assist device Fast Radial Basis Functions for Engineering Applications A comparison of mesh morphing methods for 3D shape optimization Virtual human bone modelling by interactive sculpting, mesh morphing and force-feedback Improvement of radial basis function interpolation performance on cranial implant design One mesh to rule them all: registration-based personalized cardiac flow simulations An image-based and RBF mesh morphing CFD simulation for parametric aTAA hemodynamics Computational fluid dynamic study for aTAA hemodynamics: an integrated image-based and radial basis functions mesh morphing approach A coupled CFD and RBF mesh morphing technique as surrogate for one-way FSI study Fast high fidelity CFD/CSM fluid structure interaction using RBF mesh morphing and modal superposition method Static aeroelastic analysis of an aircraft wind-tunnel model by means of modal RBF mesh updating Validation of high fidelity computational methods for aeronautical FSI analyses Higher-order time integration through smooth mesh deformation for 3D fluid-structure interaction simulations Fluid structure interaction analysis: vortex shedding induced vibrations Unsteady FSI analysis of a square array of tubes in water crossflow Aeroelastic analysis of aircraft wind-tunnel model coupling structural and fluid dynamic codes Fluid-structure coupled analysis of a transport aircraft and flight-test validation A balanced load mapping method based on radial basis functions and fuzzy sets Comparative fluid-structure interaction analysis of polymeric transcatheter and surgical aortic valves' hemodynamics and structural mechanics Polymeric heart valves for surgical implantation, catheter-based technologies and heart assist devices Polymeric heart valves: new materials, emerging hopes An arbitrary Lagrangian-Eulerian computing method for all flow speeds Analysis of fluid-structure interaction problems with structural buckling and large domain changes by ALE finite element method A new ALE adaptive unstructured methodology for the simulation of moving bodies A study of recirculating flow computation using body-fitted coordinates: consistency aspects and mesh skewness Mesh morphing techniques in CFD Advantages of radial basis function networks for dynamic system design Mesh morphing and smoothing by means of radial basis functions (RBF): a practical example using fluent and RBF morph RBF Morph for FLUENT: User's Guide. Release V1 Clinical anatomy of the aortic root Handbook of Cardiac Anatomy, Physiology, and Devices Design and manufacture of a polyvinyl alcohol (PVA) cryogel tri-leaflet heart valve prosthesis The anatomy of the sinus of Valsalva Numerical simulation of 3D unsteady flow in centrifugal pump by dynamic mesh technique Sinotubular junction size affects aortic root geometry and aortic valve function in the aortic valve reimplantation procedure: an in vitro study using the Valsalva graft Mesh deformation using the biharmonic operator Acknowledgment. The research has received funding from the European Union's Horizon 2020 research and innovation programme under the Marie Sk lodowska-Curie grant agreement No. 859836 and has been partially supported by RBF Morph R .