Immersogeometric Analysis Research Articles

Interpolation-based immersogeometric analysis methods for multi-material and multi-physics problems

AbstractImmersed boundary methods are high-order accurate computational tools used to model geometrically complex problems in computational mechanics. While traditional finite element methods require the construction of high-quality boundary-fitted meshes, immersed boundary methods instead embed the computational domain in a structured background grid. Interpolation-based immersed boundary methods augment existing finite element software to non-invasively implement immersed boundary capabilities through extraction. Extraction interpolates the structured background basis as a linear combination of Lagrange polynomials defined on a foreground mesh, creating an interpolated basis that can be easily integrated by existing methods. This work extends the interpolation-based immersed isogeometric method to multi-material and multi-physics problems. Beginning from level-set descriptions of domain geometries, Heaviside enrichment is implemented to accommodate discontinuities in state variable fields across material interfaces. Adaptive refinement with truncated hierarchically refined B-splines (THB-splines) is used to both improve interface geometry representations and to resolve large solution gradients near interfaces. Multi-physics problems typically involve coupled fields where each field has unique discretization requirements. This work presents a novel discretization method for coupled problems through the application of extraction, using a single foreground mesh for all fields. Numerical examples illustrate optimal convergence rates for this method in both 2D and 3D, for partial differential equations representing heat conduction, linear elasticity, and a coupled thermo-mechanical problem. The utility of this method is demonstrated through image-based analysis of a composite sample, where in addition to circumventing typical meshing difficulties, this method reduces the required degrees of freedom when compared to classical boundary-fitted finite element methods.

Photogrammetry-based computational fluid dynamics

Computational fluid dynamics (CFD) is the cornerstone of the design and analysis process in many engineering applications. Not only has it been applied in the design phase, but it has also been employed for analyzing the fluid flow phenomena during the operation phase for many in-use structures, such as vehicles, buildings, and landscapes. However, creating a 3D mesh-based model of in-use structures that can be used by conventional boundary-fitted CFD methods is labor-intensive, time-consuming, and sometimes impossible. Due to the challenges introduced by geometry complexity and lack of design information, it is often difficult to perform an accurate and efficient CFD analysis of these objects. This paper aims to overcome such challenges by proposing a novel photogrammetry-based CFD framework for simulating in-use structures whose design models and analysis meshes are hard to obtain. The proposed framework integrates machine learning-based 3D point cloud reconstruction of structures from 2D images obtained from portable devices (e.g., cell phones and drones) and an immersogeometric approach that can carry out flow analysis directly on reconstructed point clouds. We first present the technical details of point cloud reconstruction techniques and the immersogeometric analysis method. We then simulate the flow past a standard 12 oz soda can reconstructed using photogrammetry and compare the results with reference solutions to assess the accuracy of the approach. Finally, the proposed photogrammetry-based CFD is applied to simulations of a bell tower and the Kavita and Lalit Bahl Smart Bridge on the University of Illinois Urbana-Champaign (UIUC) campus to demonstrate the robustness of the framework and its applicability to real-world in-use civil structures.

Direct immersogeometric fluid flow and heat transfer analysis of objects represented by point clouds

Immersogeometric analysis (IMGA) is a geometrically flexible method that enables one to perform multiphysics analysis directly using complex computer-aided design (CAD) models. While the IMGA approach is well-studied and has a remarkable advantage over traditional CFD, IMGA still requires a well-defined B-rep model to represent the geometry. Obtaining such a model can sometimes be equally as challenging as creating a body-fitted mesh. To address this issue, we develop a novel IMGA approach for the simulation of incompressible and compressible flows around complex geometries represented by point clouds in this work. The point cloud representation of geometries is a direct method for digitally acquiring geometric information using LiDAR scanners, optical scanners, or other passive methods such as multi-view stereo images. The point cloud object’s geometry is represented using a set of unstructured points in the Euclidean space with (possible) orientation information in the form of surface normals. Due to the absence of topological information in the point cloud model, there are no guarantees for the geometric representation to be watertight or 2-manifold or to have consistent normals. To perform IMGA directly using point cloud geometries, we first develop a method for estimating the inside–outside information and the surface normals directly from the point cloud. We also propose a method to compute the Jacobian determinant for the surface integration (over the point cloud) necessary for the weak enforcement of Dirichlet boundary conditions. We validate these geometric estimation methods by comparing the geometric quantities computed from the point cloud with those obtained from analytical geometry and tessellated CAD models. In this work, we also develop thermal IMGA to simulate heat transfer in the presence of flow over complex geometries. The proposed framework is tested for a wide range of Reynolds and Mach numbers on benchmark problems of geometries represented by point clouds, showing the robustness and accuracy of the method. Finally, we demonstrate the applicability of our approach by performing IMGA on large industrial-scale construction machinery represented using a point cloud of more than 12 million points.

Computational Cardiovascular Medicine With Isogeometric Analysis

Isogeometric analysis (IGA) brought superior accuracy to computations in both fluid and solid mechanics. The increased accuracy has been in representing both the problem geometry and the variables computed. Beyond using IGA basis functions in space, with IGA basis functions in time in a space–time (ST) context, we can have increased accuracy also in representing the motion of solid surfaces. Around the core methods such as the residual-based variational multiscale (VMS), ST-VMS and arbitrary Lagrangian–Eulerian VMS methods, with complex-geometry IGA mesh generation methods and immersogeometric analysis, and with special methods targeting specific classes of computations, the IGA has been very effective in computational cardiovascular medicine. We provide an overview of these IGA-based computational cardiovascular-medicine methods and present examples of the computations performed.This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium provided the original work is properly cited.

High-Accuracy Mesh-Free Quadrature for Trimmed Parametric Surfaces and Volumes

This work presents a high-accuracy, mesh-free, generalized Stokes theorem-based numerical quadrature scheme for integrating functions over trimmed parametric surfaces and volumes. The algorithm relies on two fundamental steps: (1) We iteratively reduce the dimensionality of integration using the generalized Stokes theorem to line integrals over trimming curves, and (2) we employ numerical antidifferentiation in the generalized Stokes theorem using high-order quadrature rules. The scheme achieves exponential convergence up to trimming curve approximation error and has applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order curvilinear meshes, and initialization of multi-material simulations. We compare the quadrature scheme to commonly-used quadrature schemes in the literature and show that our scheme is much more efficient in terms of number of quadrature points used. We provide an open-source implementation of the scheme in MATLAB as part of QuaHOG, a software package for Quadrature of High-Order Geometries.

Industrial scale Large Eddy Simulations with adaptive octree meshes using immersogeometric analysis

We present a variant of the immersed boundary method integrated with octree meshes for highly efficient and accurate Large Eddy Simulations (LES) of flows around complex geometries. We demonstrate the scalability of the proposed method up to O(32K) processors. This is achieved by (a) rapid in-out tests; (b) adaptive quadrature for an accurate evaluation of forces; (c) tensorized evaluation during matrix assembly. We showcase this method on two non-trivial applications: accurately computing the drag coefficient of a sphere across Reynolds numbers 1−106 encompassing the drag crisis regime; simulating flow features across a semi-truck for investigating the effect of platooning on efficiency.

An octree-based immersogeometric approach for modeling inertial migration of particles in channels

In this paper, we develop a scalable, adaptively refined, octree-based finite element approach with immersogeometric analysis to track inertial migration of particles in microchannels. Fluid physics is modeled using a residual-based variational multiscale method, and the particle movement is modeled as rigid body motion. A parallel, hierarchically refined octree mesh is employed as the background mesh, on which a variationally consistent immersogeometric formulation is adopted for tracking the particle motion in the fluid. Adaptations of immersogeometric analysis on an octree background mesh are developed to enable efficient searching of background element for a given surface quadrature point, as well as a distribution of surface quadrature points over processors to reduce memory overhead and better parallelize the surface assembly. An octree-based adaptive mesh refinement algorithm adapted to in-out test in the immersogeometric approach is also developed. The validation of our octree-based immersogeometric approach is carried out using a benchmark case of a sphere settling in quiescent fluid, with good agreement presented. In addition, good strong (and weak) scalability on supercomputing resources for this benchmark case up to 16,384 processes is demonstrated. The proposed method is further deployed for exploring particle migration in straight and converging-diverging channels. This example illustrates the potential of the octree-based immersogeometric approach for efficiently tracking particle motion in complex channel flows – a problem with a diverse array of applications.

Spectral Mesh-Free Quadrature for Planar Regions Bounded by Rational Parametric Curves

This work presents spectral, mesh-free, Green’s theorem-based numerical quadrature schemes for integrating functions over planar regions bounded by rational parametric curves. Our algorithm proceeds in two steps: (1) We first find intermediate quadrature rules for line integrals along the region’s boundary curves corresponding to Green’s theorem. (2) We then use a high-order quadrature rule to compute the numerical antiderivative of the integrand along a coordinate axis, which is used to evaluate the Green’s theorem line integral. We present two methods to compute the intermediate quadrature rule. The first is spectrally accurate (it converges faster than any algebraic order with respect to number of quadrature points) and is relatively easy to implement, but has no guarantee of polynomial exactness. The second guarantees exactness for polynomial integrands up to a pre-specified degree k with an a priori-known number of quadrature points and retains the convergence properties of the first, but is slightly more complicated. The quadrature schemes have applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order meshes, and initialization of multi-material simulations with rational geometry. We compare the quadrature schemes produced using our method to other methods in the literature and show that they are much more efficient both in terms of number of quadrature points and computational time. We provide an open-source implementation of the algorithm in MATLAB.

Heart valve isogeometric sequentially-coupled FSI analysis with the space\u2013time topology change method

Heart valve fluid–structure interaction (FSI) analysis is one of the computationally challenging cases in cardiovascular fluid mechanics. The challenges include unsteady flow through a complex geometry, solid surfaces with large motion, and contact between the valve leaflets. We introduce here an isogeometric sequentially-coupled FSI (SCFSI) method that can address the challenges with an outcome of high-fidelity flow solutions. The SCFSI analysis enables dealing with the fluid and structure parts individually at different steps of the solutions sequence, and also enables using different methods or different mesh resolution levels at different steps. In the isogeometric SCFSI analysis here, the first step is a previously computed (fully) coupled Immersogeometric Analysis FSI of the heart valve with a reasonable flow solution. With the valve leaflet and arterial surface motion coming from that, we perform a new, higher-fidelity fluid mechanics computation with the space–time topology change method and isogeometric discretization. Both the immersogeometric and space–time methods are variational multiscale methods. The computation presented for a bioprosthetic heart valve demonstrates the power of the method introduced.

Computational Cardiovascular Flow Analysis with the Variational Multiscale Methods

Computational cardiovascular flow analysis can provide valuable information to medical doctors in a wide range of patientspecific cases, including cerebral aneurysms, aortas and heart valves. The computational challenges faced in this class of flow analyses also have a wide range. They include unsteady flows, complex cardiovascular geometries, moving boundaries and interfaces, such as the motion of the heart valve leaflets, contact between moving solid surfaces, such as the contact between the leaflets, and the fluid–structure interaction between the blood and the cardiovascular structure. Many of these challenges have been or are being addressed by the Space–Time Variational Multiscale (ST-VMS) method, Arbitrary Lagrangian–Eulerian VMS (ALE-VMS) method, and the VMS-based Immersogeometric Analysis (IMGA-VMS), which serve as the core computational methods, and the special methods used in combination with them. We provide an overview of the core and special methods and present examples of challenging computations carried out with these methods, including aorta and heart valve flow analyses. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium provided the original work is properly cited.

Immersogeometric analysis of moving objects in incompressible flows

We deploy the immersogeometric approach for tracking moving objects. The method immerses objects into non-boundary-fitted meshes and weakly enforces Dirichlet boundary conditions on the object boundaries. The object motion is driven by the integrated surface force and external body forces. A residual-based variational multiscale method is employed to stabilize the finite element formulation for incompressible flows. Adaptively refined quadrature rules are used to better capture the geometry of the immersed boundaries by accurately integrating the intersected background elements. Treatment for the freshly-cleared nodes (i.e. background mesh nodes that are inside the object at one time step, but are in the fluid domain at the next time step) is considered. We assess the accuracy of the method by analyzing object motion in different flow structures including objects freely dropping in viscous fluids and particle focusing in unobstructed and obstructed micro-channels. We show that key quantities of interest are in very good agreements with analytical, numerical and experimental solutions. We also show a much better computational efficiency of this framework than current commercial codes using adaptive boundary-fitted approaches. We anticipate deploying this framework for applications of particle inertial migration in microfluidic channels.

Immersogeometric analysis of compressible flows with application to aerodynamic simulation of rotorcraft

We present an immersogeometric analysis (IMGA) approach for the simulation of compressible flows around complex geometries. In this method, compressible flow simulations are performed directly on various boundary representations (B-reps) of mechanical designs, circumventing the labor-intensive and time-consuming cleanup of complex geometric models. A new formulation for the weak imposition of essential boundary conditions in the context of non-body-fitted meshes is proposed. The formulation employs the non-symmetric Nitsche method, which yields good performance especially when the penalty parameters are difficult to estimate. We test the proposed immersogeometric formulation on benchmark problems for a wide range of Reynolds and Mach numbers, showing its robustness and accuracy. Finally, the methodology is applied to the simulation of a UH-60 Black Hawk helicopter in forward flight, illustrating the ability of the proposed approach to support the design of real-world engineering systems through high-fidelity aerodynamic analysis.

Skeleton-stabilized immersogeometric analysis for incompressible viscous flow problems

A Skeleton-stabilized ImmersoGeometric Analysis technique is proposed for incompressible viscous flow problems with moderate Reynolds numbers. The proposed formulation fits within the framework of the finite cell method, where essential boundary conditions are imposed weakly using a Nitsche-type method. The key idea of the proposed formulation is to stabilize the jumps of high-order derivatives of variables over the skeleton of the background mesh. The formulation allows the use of identical finite-dimensional spaces for the approximation of the pressure and velocity fields in immersed domains. The stability issues observed for inf–sup stable discretizations of immersed incompressible flow problems are avoided with this formulation. For B-spline basis functions of degree k with highest regularity, only the derivative of order k has to be controlled, which requires specification of only a single stabilization parameter for the pressure field. The Stokes and Navier–Stokes equations are studied numerically in two and three dimensions using various immersed test cases. Oscillation-free solutions and high-order optimal convergence rates can be obtained. The formulation is shown to be stable even in limit cases where almost every elements of the physical domain is cut, and hence it does not require the existence of interior cells. In terms of the sparsity pattern, the algebraic system has a considerably smaller stencil than counterpart approaches based on Lagrange basis functions. This important property makes the proposed skeleton-stabilized technique computationally practical. To demonstrate the stability and robustness of the method, we perform a simulation of fluid flow through a porous medium, of which the geometry is directly extracted from 3D μCT scan data.

Rapid B-rep model preprocessing for immersogeometric analysis using analytic surfaces

Computational fluid dynamics (CFD) simulations of flow over complex objects have been performed traditionally using fluid-domain meshes that conform to the shape of the object. However, creating shape conforming meshes for complicated geometries such as automobiles require extensive geometry preprocessing. This process is usually tedious and requires modifying the geometry, including specialized operations such as defeaturing and filling of small gaps. Hsu et al. (2016) developed a novel immersogeometric fluid-flow method that does not require the generation of a boundary-fitted mesh for the fluid domain. However, their method used the NURBS parameterization of the surfaces for generating the surface quadrature points to enforce the boundary conditions, which required the B-rep model to be converted completely to NURBS before analysis can be performed. This conversion usually leads to poorly parameterized NURBS surfaces and can lead to poorly trimmed or missing surface features. In addition, converting simple geometries such as cylinders to NURBS imposes a performance penalty since these geometries have to be dealt with as rational splines. As a result, the geometry has to be inspected again after conversion to ensure analysis compatibility and can increase the computational cost. In this work, we have extended the immersogeometric method to generate surface quadrature points directly using analytic surfaces. We have developed quadrature rules for all four kinds of analytic surfaces: planes, cones, spheres, and tori. We have also developed methods for performing adaptive quadrature on trimmed analytic surfaces. Since analytic surfaces have frequently been used for constructing solid models, this method is also faster to generate quadrature points on real-world geometries than using only NURBS surfaces. To assess the accuracy of the proposed method, we perform simulations of a benchmark problem of flow over a torpedo shape made of analytic surfaces and compare those to immersogeometric simulations of the same model with NURBS surfaces. We also compare the results of our immersogeometric method with those obtained using boundary-fitted CFD of a tessellated torpedo shape, and quantities of interest such as drag coefficient are in good agreement. Finally, we demonstrate the effectiveness of our immersogeometric method for high-fidelity industrial scale simulations by performing an aerodynamic analysis of a truck that has a large percentage of analytic surfaces. Using analytic surfaces over NURBS avoids unnecessary surface type conversion and significantly reduces model-preprocessing time, while providing the same accuracy for the aerodynamic quantities of interest.

Optimizing fluid–structure interaction systems with immersogeometric analysis and surrogate modeling: Application to a hydraulic arresting gear

This work describes a fluid–structure interaction (FSI) design optimization framework and applies it to improving the structural performance of a water brake used to stop aircraft landing on short runways. Inside the water brake, a dissipative torque is exerted on a rotor through interactions between rotor blades and a surrounding fluid. We seek to optimize blade shape over a parameterized design space, to prevent potentially-damaging stress concentrations without compromising performance. To avoid excessive numbers of costly simulations while exploring the design space, we use a surrogate management framework that combines derivative-free pattern search optimization with automated construction of a low-fidelity surrogate model, requiring only a handful of high-fidelity FSI simulations. We avoid the difficult problem of generating fluid and structure meshes at new points in the design space by using immersogeometric FSI analysis. The structure is analyzed isogeometrically: its design geometry also serves as a computational mesh. This geometry is then immersed in an unfitted fluid mesh that does not depend on the structure’s design parameters. We use this framework to make significant improvements to a baseline design found in the literature. Specifically, there is a 35% reduction of von Mises stress variance and a 25% reduction of maximum of stress, while the resisting torque and mass of the optimized blades remain uncompromised.

The tetrahedral finite cell method: Higher‐order immersogeometric analysis on adaptive non‐boundary‐fitted meshes

SummaryThe finite cell method (FCM) is an immersed domain finite element method that combines higher‐order non‐boundary‐fitted meshes, weak enforcement of Dirichlet boundary conditions, and adaptive quadrature based on recursive subdivision. Because of its ability to improve the geometric resolution of intersected elements, it can be characterized as an immersogeometric method. In this paper, we extend the FCM, so far only used with Cartesian hexahedral elements, to higher‐order non‐boundary‐fitted tetrahedral meshes, based on a reformulation of the octree‐based subdivision algorithm for tetrahedral elements. We show that the resulting TetFCM scheme is fully accurate in an immersogeometric sense, that is, the solution fields achieve optimal and exponential rates of convergence for h‐refinement and p‐refinement, if the immersed geometry is resolved with sufficient accuracy. TetFCM can leverage the natural ability of tetrahedral elements for local mesh refinement in three dimensions. Its suitability for problems with sharp gradients and highly localized features is illustrated by the immersogeometric phase‐field fracture analysis of a human femur bone. Copyright © 2016 John Wiley & Sons, Ltd.

The tetrahedral finite cell method for fluids: Immersogeometric analysis of turbulent flow around complex geometries

Abstract We present a tetrahedral finite cell method for the simulation of incompressible flow around geometrically complex objects. The method immerses such objects into non-boundary-fitted meshes of tetrahedral finite elements and weakly enforces Dirichlet boundary conditions on the objects’ surfaces. Adaptively-refined quadrature rules faithfully capture the flow domain geometry in the discrete problem without modifying the non-boundary-fitted finite element mesh. A variational multiscale formulation provides accuracy and robustness in both laminar and turbulent flow conditions. We assess the accuracy of the method by analyzing the flow around an immersed sphere for a wide range of Reynolds numbers. We show that quantities of interest such as the drag coefficient, Strouhal number and pressure distribution over the sphere are in very good agreement with reference values obtained from standard boundary-fitted approaches. We place particular emphasis on studying the importance of the geometry resolution in intersected elements. Aligning with the immersogeometric concept, our results show that the faithful representation of the geometry in intersected elements is critical for accurate flow analysis. We demonstrate the potential of our proposed method for high-fidelity industrial scale simulations by performing an aerodynamic analysis of an agricultural tractor.

An immersogeometric variational framework for fluid–structure interaction: Application to bioprosthetic heart valves

In this paper, we develop a geometrically flexible technique for computational fluid–structure interaction (FSI). The motivating application is the simulation of tri-leaflet bioprosthetic heart valve function over the complete cardiac cycle. Due to the complex motion of the heart valve leaflets, the fluid domain undergoes large deformations, including changes of topology. The proposed method directly analyzes a spline-based surface representation of the structure by immersing it into a non-boundary-fitted discretization of the surrounding fluid domain. This places our method within an emerging class of computational techniques that aim to capture geometry on non-boundary-fitted analysis meshes. We introduce the term “immersogeometric analysis” to identify this paradigm.The framework starts with an augmented Lagrangian formulation for FSI that enforces kinematic constraints with a combination of Lagrange multipliers and penalty forces. For immersed volumetric objects, we formally eliminate the multiplier field by substituting a fluid–structure interface traction, arriving at Nitsche’s method for enforcing Dirichlet boundary conditions on object surfaces. For immersed thin shell structures modeled geometrically as surfaces, the tractions from opposite sides cancel due to the continuity of the background fluid solution space, leaving a penalty method. Application to a bioprosthetic heart valve, where there is a large pressure jump across the leaflets, reveals shortcomings of the penalty approach. To counteract steep pressure gradients through the structure without the conditioning problems that accompany strong penalty forces, we resurrect the Lagrange multiplier field. Further, since the fluid discretization is not tailored to the structure geometry, there is a significant error in the approximation of pressure discontinuities across the shell. This error becomes especially troublesome in residual-based stabilized methods for incompressible flow, leading to problematic compressibility at practical levels of refinement. We modify existing stabilized methods to improve performance.To evaluate the accuracy of the proposed methods, we test them on benchmark problems and compare the results with those of established boundary-fitted techniques. Finally, we simulate the coupling of the bioprosthetic heart valve and the surrounding blood flow under physiological conditions, demonstrating the effectiveness of the proposed techniques in practical computations.