Search

This paper uses a divergence-conforming B-spline fluid discretization to address the long-standing issue of poor mass conservation in immersed methods for computational fluid–structure interaction (FSI) that represent the influence of the structure as a forcing term in the fluid subproblem. We focus, in particular, on the immersogeometric method developed in our earlier work, analyze its convergence for linear model problems, then apply it to FSI analysis of heart valves, using divergence-conforming B-splines to discretize the fluid subproblem. Poor mass conservation can manifest as effective leakage of fluid through thin solid barriers. This leakage disrupts the qualitative behavior of FSI systems such as heart valves, which exist specifically to block flow. Divergence-conforming discretizations can enforce mass conservation exactly, avoiding this problem. To demonstrate the practical utility of immersogeometric FSI analysis with divergence-conforming B-splines, we use the methods described in this paper to construct and evaluate a computational model of an in vitro experiment that pumps water through an artificial valve. [ABSTRACT FROM AUTHOR]

We present formulations for compressible and incompressible hyperelastic thin shells which can use general 3D constitutive models. The necessary plane stress condition is enforced analytically for incompressible materials and iteratively for compressible materials. The thickness stretch is statically condensed and the shell kinematics are completely described by the first and second fundamental forms of the midsurface. We use C 1 -continuous isogeometric discretizations to build the numerical models. Numerical tests, including structural dynamics simulations of a bioprosthetic heart valve, show the good performance and applicability of the presented methods. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

There is growing interest in high fidelity simulations of buoyancy-driven flows in indoor environments. This is driven to a large extent by the push to utilize naturally occurring flows to reduce energy consumption in the built environment. Naturally occurring (buoyancy-driven) flows exhibit spatiotemporal thermal fluctuations, and accurately capturing these thermal variations is essential to calculate energy usage. In this work, we deploy a residual-based variational multiscale (VMS) finite element large-eddy simulation model for accurately modeling buoyancy-driven flows in enclosed environments. We use the canonical example of the Rayleigh–Bénard convection problem in 2D and 3D to verify and validate the VMS model. We show good comparison with benchmark numerical and experimental results across seven orders of magnitude variation in Rayleigh numbers (R a ∼ 1 0 3 to 1 0 10 ), covering laminar, transition, and turbulent regimes without any extra treatments. We additionally employ weak enforcement of Dirichlet boundary conditions for both velocity and temperature, and show that comparable results can be produced with much coarser meshes. This confirms that the VMS framework with weak imposition of Dirichlet boundary conditions is a computationally efficient approach to model buoyancy-driven flow physics in complex indoor environments. • VMS for buoyancy driven flow shows good results for a Rayleigh–Bénard problem. • VMS is able to cover laminar, transition and turbulence without any extra treatments. • Weak imposition of boundary conditions is able to drastically decrease mesh density. [ABSTRACT FROM AUTHOR]

We explore the use of the non-symmetric Nitsche method for the weak imposition of boundary and coupling conditions along interfaces that intersect through a finite element mesh. In contrast to symmetric Nitsche methods, it does not require stabilization and therefore does not depend on the appropriate estimation of stabilization parameters. We first review the available mathematical background, recollecting relevant aspects of the method from a numerical analysis viewpoint. We then compare accuracy and convergence of symmetric and non-symmetric Nitsche methods for a Laplace problem, a Kirchhoff plate, and in 3D elasticity. Our numerical experiments confirm that the non-symmetric method leads to reduced accuracy in the L 2 error, but exhibits superior accuracy and robustness for derivative quantities such as diffusive flux, bending moments or stress. Based on our numerical evidence, the non-symmetric Nitsche method is a viable alternative for problems with diffusion-type operators, in particular when the accuracy of derivative quantities is of primary interest. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Abstract Isogeometric analysis (IGA) has been a particularly impactful development in the realm of Kirchhoff–Love thin-shell analysis because the high-order basis functions employed naturally satisfy the requirement of C 1 continuity. Still, engineering models of appreciable complexity, such as wind turbine blades, are typically modeled using multiple surface patches and, often, neither rotational continuity nor conforming discretization can be practically obtained at patch interfaces. A penalty approach for coupling adjacent patches is therefore presented. The proposed method imposes both displacement and rotational continuity and is applicable to either smooth or non-smooth interfaces and either matching or non-matching discretization. The penalty formulations require only a single, dimensionless penalty coefficient for both displacement and rotation coupling terms, alleviating the problem-dependent nature of the penalty parameters. Using this coupling methodology, numerous benchmark problems encapsulating a variety of analysis types, geometrical and material properties, and matching and non-matching interfaces are addressed. The coupling methodology produces consistently accurate results throughout all tests. Furthermore, the suggested penalty coefficient of α = 1 0 3 is shown to be effective for the wide range of problem configurations addressed. Finally, a realistic wind turbine blade model, consisting of 27 patches and 51 coupling interfaces and having a chordwise- and spanwise-variant composite material definition, is subjected to buckling, vibration, and nonlinear deformation analyses using the proposed approach. [ABSTRACT FROM AUTHOR]

This work formulates frictionless contact between solid bodies in terms of a repulsive potential energy term and illustrates how numerical integration of the resulting forces is computationally similar to the “pinball algorithm” proposed and studied by Belytschko and collaborators in the 1990s. We thereby arrive at a numerical approach that has both the theoretical advantages of a potential-based formulation and the algorithmic simplicity, computational efficiency, and geometrical versatility of pinball contact. The singular nature of the contact potential requires a specialized nonlinear solver and an adaptive time stepping scheme to ensure reliable convergence of implicit dynamic calculations. We illustrate the effectiveness of this numerical method by simulating several benchmark problems and the structural mechanics of the right atrioventricular (tricuspid) heart valve. Atrioventricular valve closure involves contact between every combination of shell surfaces, edges of shells, and cables, but our formulation handles all contact scenarios in a unified manner. We take advantage of this versatility to demonstrate the effects of chordal rupture on tricuspid valve coaptation behavior. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Isogeometric analysis (IGA) fundamentally seeks to bridge the gap between engineering design and high-fidelity computational analysis by using spline functions as finite element bases. However, additional computational design paradigms must be taken into consideration to ensure that designers can take full advantage of IGA, especially within the context of design optimization. In this work, we propose a novel approach that employs IGA methodologies while still rigorously abiding by the paradigms of advanced design parameterization, analysis model validity, and interactivity. The entire design lifecycle utilizes a consistent geometry description and is contained within a single platform. Because of this unified workflow, iterative design optimization can be naturally integrated. The proposed methodology is demonstrated through an IGA-based parametric design optimization framework implemented using the Grasshopper algorithmic modeling interface for Rhinoceros 3D. The framework is capable of performing IGA-based design optimization of realistic engineering structures that are practically constructed through the use of complex geometric operations. We demonstrate the framework’s effectiveness on both an internally pressurized tube and a wind turbine blade, highlighting its applicability across a spectrum of design complexity. In addition to inherently featuring the advantageous characteristics of IGA, the seamless nature of the workflow instantiated in this framework diminishes the obstacles traditionally encountered when performing finite-element-analysis-based design optimization. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The computational modeling of thin-walled structures based on isogeometric analysis (IGA), non-uniform rational B-splines (NURBS), and Kirchhoff–Love (KL) shell formulations has attracted significant research attention in recent years. While these methods offer numerous benefits over the traditional finite element approach, including exact representation of the geometry, naturally satisfied high-order continuity within each NURBS patch, and computationally efficient rotation-free formulations, they also present a number of challenges in modeling real-world engineering structures of considerable complexity. Specifically, these NURBS-based engineering models are usually comprised of numerous patches, with discontinuous derivatives, non-conforming discretizations, and non-watertight connections at their interfaces. Moreover, the analysis of such structures often requires the full stress and strain tensors (i.e., including the transverse normal and shear components) for subsequent failure analysis and remaining life prediction. Despite the efficiency provided by the KL shell, the formulation cannot accurately predict the response in the transverse directions due to its kinematic assumptions. In this work, a penalty-based formulation for the blended coupling of KL and continuum shells is presented. The proposed approach embraces both the computational efficiency of KL shells and the availability of the full-scale stress/strain tensors of continuum shells where needed by modeling critical structural components using continuum shells and other components using KL shells. The proposed method enforces the displacement and rotational continuities in a variational manner and is applicable to non-conforming and non-smooth interfaces. The efficacy of the developed method is demonstrated through a number of benchmark studies with a variety of analysis configurations, including linear and nonlinear analyses, matching and non-matching discretizations, and isotropic and composite materials. Finally, an aircraft horizontal stabilizer is considered to demonstrate the applicability of the proposed blended shells to real-world engineering structures of significant complexity. [ABSTRACT FROM AUTHOR]

Approximately 1.6 million patients in the United States are affected by tricuspid valve regurgitation, which occurs when the tricuspid valve does not close properly to prevent backward blood flow into the right atrium. Despite its critical role in proper cardiac function, the tricuspid valve has received limited research attention compared to the mitral and aortic valves on the left side of the heart. As a result, proper valvular function and the pathologies that may cause dysfunction remain poorly understood. To promote further investigations of the biomechanical behavior and response of the tricuspid valve, this work establishes a parameter-based approach that provides a template for tricuspid valve modeling and simulation. The proposed tricuspid valve parameterization presents a comprehensive description of the leaflets and the complex chordae tendineae for capturing the typical three-leaflet structural deformation observed from medical data. This simulation framework develops a practical procedure for modeling tricuspid valves and offers a robust, flexible approach to analyze the performance and effectiveness of various valve configurations using isogeometric analysis. The proposed methods also establish a baseline to examine the tricuspid valve's structural deformation, perform future investigations of native valve configurations under healthy and disease conditions, and optimize prosthetic valve designs. [ABSTRACT FROM AUTHOR]

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- α 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. • An immersogeometric formulation for free-surface flow is proposed. • Level set method is adopted for fluid–fluid interface. • FCM with adaptive quadrature is adopted for fluid–structure interface. • An efficient computational algorithm is outlined. • Simulations of marine engineering problems are presented. [ABSTRACT FROM AUTHOR]

The transcatheter aortic valve replacement (TAVR) has emerged as a minimally invasive alternative to surgical treatments of valvular heart disease. TAVR offers many advantages, however, the safe anchoring of the transcatheter heart valve (THV) in the patient's anatomy is key to a successful procedure. In this paper, we develop and apply a novel immersogeometric fluid–structure interaction (FSI) framework for the modeling and simulation of the TAVR procedure to study the anchoring ability of the THV. To account for physiological realism, methods are proposed to model and couple the main components of the system, including the arterial wall, blood flow, valve leaflets, skirt, and frame. The THV is first crimped and deployed into an idealized ascending aorta. During the FSI simulation, the radial outward force and friction force between the aortic wall and the THV frame are examined over the entire cardiac cycle. The ratio between these two forces is computed and compared with the experimentally estimated coefficient of friction to study the likelihood of valve migration. • An immersogeometric analysis framework for transcatheter heart valves is proposed. • Components of the transcatheter heart valves are modeled and coupled. • Validation of the proposed structural and fluid mechanics methods is presented. • The framework is applied to the simulation of transcatheter aortic valve replacement. • Static friction is considered to study the risk of migration. [ABSTRACT FROM AUTHOR]