Immersogeometric analysis of compressible flows

This dissertation presents the development of a novel immersogeometric method for the simulation of turbulent compressible flows around complex geometries. The immersogeometric analysis is first extended into the version of tetrahedral finite cell method, in order to handle complex geometries flexibly and accurately. The developed method immerses complex objects into non-boundary-fitted meshes of tetrahedral finite elements which can be easily refined in interesting regions. Adaptively-refined quadrature rules faithfully capture the flow do- main geometry in the discrete problem without modifying the non-boundary-fitted finite element mesh. Particular emphasis is placed on studying the importance of the geometry resolution in in- tersected elements. Aligning with the immersogeometric concept, the results show that the faithful representation of the geometry in intersected elements is critical for accurate flow analysis. To simulate the compressible flows in an accurate and and robust way, a novel stabilized finite element formulation is developed. New weak imposition of essential boundary conditions and sliding-interface formulations are also proposed in the context of moving-domain compressible flows. The new formulation is successfully tested on a set of examples spanning a wide range of Reynolds and Mach numbers showing its superior robustness. Experimental validation of the new formulation is also carried out with good success. The developments of tetrahedral finite cell method and the stabilized finite element formulations are combined to further develop the immersogeometric method for compressible flows. Non-symmetric Nitsche method is used in the weak-boundary-condition operator, to offer good performance in the context of non-boundary-fitted discretization. The developed immersogeomet- ric method is tested against several benchmark problems, to prove its comparable accuracy to its boundary-fitted counterpart. Finally, the aerodynamic analysis of a UH-60 helicopter is carried outusing the developed method, to illustrate its potential to support design of real engineering systems through high-fidelity aerodynamic analysis

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. In this paper, we develop a novel IMGA approach for simulating incompressible and compressible flows around complex geometries represented by point clouds. 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.Comment: 30 pages + references; Accepted in Computer Methods in Applied Mechanics and Engineerin

An immersogeometric formulation for free-surface flows with application to marine engineering problems

An immersogeometric formulation is proposed to simulate free-surface flows around structures with complex geometry. The fluid–fluid interface (air–water interface) is handled by the level set method, while the fluid–structure interface is handled through an immersogeometric approach by immersing structures into non-boundary-fitted meshes and enforcing Dirichlet boundary conditions weakly. Residual-based variational multiscale method (RBVMS) is employed to stabilize the coupled Navier–Stokes equations of incompressible flows and level set convection equation. Other level set techniques, including re-distancing and mass balancing, are also incorporated into the immersed formulation. Adaptive quadrature rule is used to better capture the geometry of the immersed structure boundary by accurately integrating the intersected background elements. Generalized-α role= presentation style= box-sizing: border-box; margin: 0px; padding: 0px; display: inline-block; line-height: normal; font-size: 16.2px; word-spacing: normal; overflow-wrap: normal; white-space: nowrap; float: none; direction: ltr; max-width: none; max-height: none; min-width: 0px; min-height: 0px; border: 0px; position: relative; \u3eα method is adopted for time integration, which results in a two-stage predictor multi-corrector algorithm. GMRES solver preconditioned with block Jacobian matrices of individual fluid and level set subproblems is used for solving the coupled linear systems arising from the multi-corrector stage. The capability and accuracy of the proposed method are assessed by simulating three challenging marine engineering problems, which are a solitary wave impacting a stationary platform, dam break with an obstacle, and planing of a DTMB 5415 ship model. A refinement study is performed. The predictions of key quantities of interest by the proposed formulation are in good agreement with experimental results and boundary-fitted simulation results from others. The proposed formulation has great potential for wide applications in marine engineering problems

Heart valve isogeometric sequentially-coupled FSI analysis with the space–time 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

A note on the penalty parameter in Nitsche's method for unfitted boundary value problems

Nitsche's method is a popular approach to implement Dirichlet-type boundary conditions in situations where a strong imposition is either inconvenient or simply not feasible. The method is widely applied in the context of unfitted finite element methods. From the classical (symmetric) Nitsche's method it is well-known that the stabilization parameter in the method has to be chosen sufficiently large to obtain unique solvability of discrete systems. In this short note we discuss an often used strategy to set the stabilization parameter and describe a possible problem that can arise from this. We show that in specific situations error bounds can deteriorate and give examples of computations where Nitsche's method yields large and even diverging discretization errors

Industrial scale large eddy simulations (LES) 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 $\mathcal{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-10^6$ encompassing the drag crisis regime; simulating flow features across a semi-truck for investigating the effect of platooning on efficiency.Comment: Accepted for publication at Computer and Mathematics with Application

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