Search

For wind turbines operating in cold weather conditions, ice accretion is an established issue that remains an obstacle in effective turbine operation. While the aerodynamic performance of wind turbine blades with ice accretion has received considerable research attention, few studies have investigated the structural impact of blade ice accretion. This work proposes an adaptable projection-based method to superimpose complex ice configurations onto a baseline structure. The proposed approach provides an efficient methodology to include ice accretion in the high-fidelity isogeometric shell analysis of a realistic wind turbine blade. Linear vibration and nonlinear deflection analyses of the blade are performed for various ice configurations to demonstrate the impact of different ice accretion distributions on structural performance. These analyses indicate decreases in the blade natural frequencies and deflection under icing conditions. Such ice-induced changes clearly reveal the need for structural design consideration for turbines operating under icing conditions.

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.

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.

Many of the challenges encountered in computational analysis of wind turbines and turbomachinery are being addressed by the Arbitrary Lagrangian–Eulerian (ALE) and Space–Time (ST) Variational Multiscale (VMS) methods and isogeometric discretization. The computational challenges include turbulent rotational flows, complex geometries, moving boundaries and interfaces, such as the rotor motion, and the fluid–structure interaction (FSI), such as the FSI between the wind turbine blade and the air. The core computational methods are the ALE-VMS and ST-VMS methods. These are supplemented with special methods like the Slip Interface (SI) method and ST Isogeometric Analysis with NURBS basis functions in time. We describe the core and special methods and present, as examples of challenging computations performed, computational analysis of horizontal- and vertical-axis wind turbines and flow-driven string dynamics in pumps.

This paper considers an anisotropic hyperelastic soft tissue model, originally proposed for native valve tissue and referred to herein as the Lee–Sacks model, in an isogeometric thin shell analysis framework that can be readily combined with immersogeometric fluid–structure interaction (FSI) analysis for high-fidelity simulations of bioprosthetic heart valves (BHVs) interacting with blood flow. We find that the Lee–Sacks model is well-suited to reproduce the anisotropic stress–strain behavior of the cross-linked bovine pericardial tissues that are commonly used in BHVs. An automated procedure for parameter selection leads to an instance of the Lee–Sacks model that matches biaxial stress–strain data from the literature more closely, over a wider range of strains, than other soft tissue models. The relative simplicity of the Lee–Sacks model is attractive for computationally-demanding applications such as FSI analysis and we use the model to demonstrate how the presence and direction of material anisotropy affect the FSI dynamics of BHV leaflets. © 2018. This is the authors’ accepted and refereed manuscript to the article. Locked until 12.4.2019 due to copyright restrictions. This manuscript version is made available under the CC-BY-NC-ND 4.0 license http://creativecommons.org/licenses/by-nc-nd/4.0/

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.

This article provides Abaqus input files and user subroutines for performing finite element simulations of the tricuspid heart valve with an idealized geometry. Additional post-processing steps to obtain a ParaView visualization file (*.vtk) of the deformed geometry are also provided to allow the readers to use the included ParaView state file (*.pvsm) for customizable visualization and evaluation of the simulation results. We expect this first-of-its-kind in-silico benchmark dataset will facilitate user-friendly simulations considering material nonlinearity, leaflet-to-leaflet contact, and large deformations. Additionally, the information included herein can be used to rapidly evaluate other novel in-silico approaches developed for simulating cardiac valve function. The benchmark can be expanded to consider more complex features of the tricuspid valve function, such as the dynamic annulus motion or the time-varying transvalvular pressure. Interested readers are referred to the companion article (Johnson et al., 2021) for an example application of this in-silico tool for isogeometric analysis of tricuspid valves.

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.

This paper discusses a method of stabilizing Lagrange multiplier fields used to couple thin immersed shell structures and surrounding fluids. The method retains essential conservation properties by stabilizing only the portion of the constraint orthogonal to a coarse multiplier space. This stabilization can easily be applied within iterative methods or semi-implicit time integrators that avoid directly solving a saddle point problem for the Lagrange multiplier field. Heart valve simulations demonstrate applicability of the proposed method to 3D unsteady simulations. An appendix sketches the relation between the proposed method and a high-order-accurate approach for simpler model problems.

In this paper, we present an isogeometric analysis framework for design space exploration. While the methodology is presented in the setting of structural mechanics, it is applicable to any system of parametric partial differential equations. The design space exploration framework elucidates design parameter sensitivities used to inform initial and early-stage design. Moreover, this framework enables the visualization of a full system response, including the displacement and stress fields throughout the domain, by providing an approximation to the system solution vector. This is accomplished through a collocation-like approach where various geometries throughout the design space under consideration are sampled. The sampling scheme follows a quadrature rule while the physical solutions to these sampled geometries are obtained through an isogeometric method. A surrogate model to the design space solution manifold is constructed through either an interpolating polynomial or pseudospectral expansion. Examples of this framework are presented with applications to the Scordelis–Lo roof, a Flat L-Bracket, and an NREL 5 MW wind turbine blade.

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.

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.

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.

Fluid–structure interaction (FSI) analysis of a full-scale hydraulic arresting gear used to retard the forward motion of an aircraft landing on an aircraft-carrier deck is performed. The simulations make use of the recently developed core and special-purpose FSI techniques for other problem classes, specialized to the present application. A recently proposed interactive geometry modeling and parametric design platform for isogeometric analysis (IGA) is directly employed to create the arresting gear model, and illustrates a natural application of IGA to this problem class. The fluid mechanics and FSI simulation results are reported in terms of the arresting-gear rotor loads and blade structural deformation and vibration. Excellent agreement is achieved with the experimental results for the arresting gear design simulated in this work.

Immersogeometric analysis that directly uses the B-rep CAD model is proposed.A GPU-accelerated point membership classification is performed.Distribution of the surface quadrature points is crucial for accuracy.The quadrature error near the trim curves is relatively insignificant.The methodology is found effective on a 3D benchmark and an industrial problem. 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's 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.

The challenges encountered in computational analysis of wind turbines and turbomachinery include turbulent rotational flows, complex geometries, moving boundaries and interfaces, such as the rotor motion, and the fluid-structure interaction (FSI), such as the FSI between the wind turbine blade and the air. The Arbitrary Lagrangian-Eulerian (ALE) and Space-Time (ST) Variational Multiscale (VMS) methods and isogeometric discretization have been effective in addressing these challenges. The ALE-VMS and ST-VMS serve as core computational methods. They are supplemented with special methods like the Slip Interface (SI) method and ST Isogeometric Analysis with NURBS basis functions in time. We describe the core and special methods and present, as examples of challenging computations performed, computational analysis of horizontal and vertical-axis wind turbines and flow-driven 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.

In the mid-2010s, we began applying a combination of isogeometric analysis and immersed boundary methods to the problem of bioprosthetic heart valve (BHV) fluid–structure interaction (FSI). This chapter reviews how our research on BHV FSI (1) crystallized the emerging concept of immersogeometric analysis, (2) introduced a new semi-implicit numerical method for weakly enforcing constraints in time dependent problems, which we refer to as the dynamic augmented Lagrangian approach, and (3) clarified the important role of mass conservation in immersed FSI analysis. We illustrate these contributions with selected numerical results and discuss future improvements to, and applications of, the computational FSI techniques we have developed.

This paper builds on a recently developed computationally tractable material model merged with an immersogeometric fluid–structure interaction methodology for bioprosthetic heart valve modeling and simulation. Our main objective is to enable improved application of the use of exogenous crosslinked tissues in prosthesis design through computational methods by utilizing physically realistic constitutive models. To enhance constitutive modeling, valve leaflets are modeled with a computationally efficient phenomenological constitutive relation stemmed from a full structural model to explore the influence of incorporating a high-fidelity material model for the leaflets. We call this phenomenological version as the effective model. This effective model constitutive form is incorporated in the context of the isogeometric analysis to develop an efficient fluid–structure interaction method for thin shell structure of the leaflet tissues. The implementation is supported by representative simulations showing the applicability and usefulness of our effective material model in heart valve simulation framework.

In this paper an interactive parametric design-through-analysis platform is proposed to help design engineers and analysts make more effective use of Isogeometric Analysis (IGA) to improve their product design and performance. We develop several Rhinoceros (Rhino) plug-ins to take input design parameters through a user-friendly interface, generate appropriate surface and/or volumetric models, perform mechanical analysis, and visualize the solution fields, all within the same Computer-Aided Design (CAD) program. As part of this effort we propose and implement graphical generative algorithms for IGA model creation and visualization based on Grasshopper, a visual programming interface to Rhino. The developed platform is demonstrated on two structural mechanics examples-an actual wind turbine blade and a model of an integrally bladed rotor (IBR). In the latter example we demonstrate how the Rhino functionality may be utilized to create conforming volumetric models for IGA.

The purpose of this study is to enhance the stability properties of our recently-developed numerical method [D. Kamensky, M.-C. Hsu, D. Schillinger, J.A. Evans, A. Aggarwal, Y. Bazilevs, M.S. Sacks, T.J.R. Hughes, “An immersogeometric variational framework for fluid-structure interaction: Application to bioprosthetic heart valves”, Comput. Methods Appl. Mech. Engrg., 284 (2015) 1005–1053] for immersing spline-based representations of shell structures into unsteady viscous incompressible flows. In the cited work, we formulated the fluid-structure interaction (FSI) problem using an augmented Lagrangian to enforce kinematic constraints. We discretized this Lagrangian as a set of collocated constraints, at quadrature points of the surface integration rule for the immersed interface. Because the density of quadrature points is not controlled relative to the fluid discretization, the resulting semi-discrete problem may be over-constrained. Semi-implicit time integration circumvents this difficulty in the fully-discrete scheme. If this time-stepping algorithm is applied to fluid-structure systems that approach steady solutions, though, we find that spatially-oscillating modes of the Lagrange multiplier field can grow over time. In the present work, we stabilize the semi-implicit integration scheme to prevent potential divergence of the multiplier field as time goes to infinity. This stabilized time integration may also be applied in pseudo-time within each time step, giving rise to a fully implicit solution method. We discuss the theoretical implications of this stabilization scheme for several simplified model problems, then demonstrate its practical efficacy through numerical examples.

This paper builds on a recently developed immersogeometric fluid---structure interaction (FSI) methodology for bioprosthetic heart valve (BHV) modeling and simulation. It enhances the proposed framework in the areas of geometry design and constitutive modeling. With these enhancements, BHV FSI simulations may be performed with greater levels of automation, robustness and physical realism. In addition, the paper presents a comparison between FSI analysis and standalone structural dynamics simulation driven by prescribed transvalvular pressure, the latter being a more common modeling choice for this class of problems. The FSI computation achieved better physiological realism in predicting the valve leaflet deformation than its standalone structural dynamics counterpart.

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.

Fluid–structure interaction (FSI) analysis of a full-scale hydraulic arresting gear used to retard the forward motion of an aircraft landing on an aircraft-carrier deck is performed. The simulations make use of the recently developed core and special-purpose FSI techniques for other problem classes, specialized to the present application. A recently proposed interactive geometry modeling and parametric design platform for isogeometric analysis (IGA) is directly employed to create the arresting gear model, and illustrates a natural application of IGA to this problem class. The fluid mechanics and FSI simulation results are reported in terms of the arresting gear rotor loads and blade structural deformation and vibration. Excellent agreement is achieved with the experimental results for the arresting gear design simulated in this work.

Numerous studies have suggested that medical image derived computational mechanics models could be developed to reduce mortality and morbidity due to cardiovascular diseases by allowing for patient-specific surgical planning and customized medical device design. In this work, we present a novel framework for designing prosthetic heart valves using a parametric design platform and immersogeometric fluid-structure interaction (FSI) analysis. We parameterize the leaflet geometry using several key design parameters. This allows for generating various perturbations of the leaflet design for the patient-specific aortic root reconstructed from the medical image data. Each design is analyzed using our hybrid arbitrary Lagrangian-Eulerian/immersogeometric FSI methodology, which allows us to efficiently simulate the coupling of the deforming aortic root, the parametrically designed prosthetic valves, and the surrounding blood flow under physiological conditions. A parametric study is performed to investigate the influence of the geometry on heart valve performance, indicated by the effective orifice area and the coaptation area. Finally, the FSI simulation result of a design that balances effective orifice area and coaptation area reasonably well is compared with patient-specific phase contrast magnetic resonance imaging data to demonstrate the qualitative similarity of the flow patterns in the ascending aorta.