Direct immersogeometric fluid flow analysis using B-rep CAD models

We present a new method for immersogeometric fluid flow analysis that directly uses the CAD boundary representation (B-rep) of a complex object and immerses it into a locally refined, non-boundary-fitted discretization of the fluid domain. The motivating applications include analyzing the flow over complex geometries, such as moving vehicles, where the detailed geometric features usually require time-consuming, labor-intensive geometry cleanup or mesh manipulation for generating the surrounding boundary-fitted fluid mesh. The proposed method avoids the challenges associated with such procedures. A new method to perform point membership classification of the background mesh quadrature points is also proposed. To faithfully capture the geometry in intersected elements, we implement an adaptive quadrature rule based on the recursive splitting of elements. Dirichlet boundary conditions in intersected elements are enforced weakly in the sense of Nitsche\u27s method. To assess the accuracy of the proposed method, we perform computations of the benchmark problem of flow over a sphere represented using B-rep. Quantities of interest such as drag coefficient are in good agreement with reference values reported in the literature. The results show that the density and distribution of the surface quadrature points are crucial for the weak enforcement of Dirichlet boundary conditions and for obtaining accurate flow solutions. Also, with sufficient levels of surface quadrature element refinement, the quadrature error near the trim curves becomes insignificant. Finally, we demonstrate the effectiveness of our immersogeometric method for high-fidelity industrial scale simulations by performing an aerodynamic analysis of an agricultural tractor directly represented using B-rep

A geometric framework for immersogeometric analysis

The purpose of this dissertation is to develop a geometric framework for immersogeometric analysis that directly uses the boundary representations (B-reps) of a complex computer-aided design (CAD) model and immerses it into a locally refined, non-boundary-fitted discretization of the fluid domain. Using the non-boundary-fitted mesh which does not need to conform to the shape of the object can alleviate the challenge of mesh generation for complex geometries. This also reduces the labor-intensive and time-consuming work of geometry cleanup for the purpose of obtaining watertight CAD models in order to perform boundary-fitted mesh generation. The Dirichlet boundary conditions in the fluid domain are enforced weakly over the immersed object surface in the intersected elements. The surface quadrature points for the immersed object are generated on the parametric and analytic surfaces of the B-rep models. In the case of trimmed surfaces, adaptive quadrature rule is considered to improve the accuracy of the surface integral. For the non-boundary-fitted mesh, a sub-cell-based adaptive quadrature rule based on the recursive splitting of quadrature elements is used to faithfully capture the geometry in intersected elements. The point membership classification for identifying quadrature points in the fluid domain is based on a voxel-based approach implemented on GPUs. A variety of computational fluid dynamics (CFD) simulations are performed using the proposed method to assess its accuracy and efficiency. Finally, a fluid--structure interaction (FSI) simulation of a deforming left ventricle coupled with the heart valves shows the potential advantages of the developed geometric framework for the immersogeomtric analysis with complex moving domains

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

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

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

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