Yuri Bazilevs

Computational vascular fluid–structure interaction: methodology and application to cerebral aneurysms

A computational vascular fluid–structure interaction framework for the simulation of patient-specific cerebral aneurysm configurations is presented. A new approach for the computation of the blood vessel tissue prestress is also described. Simulations of four patient-specific models are carried out, and quantities of hemodynamic interest such as wall shear stress and wall tension are studied to examine the relevance of fluid–structure interaction modeling when compared to the rigid arterial wall assumption. We demonstrate that flexible wall modeling plays an important role in accurate prediction of patient-specific hemodynamics. Discussion of the clinical relevance of our methods and results is provided.

Determination of Wall Tension in Cerebral Artery Aneurysms by Numerical Simulation

Background and Purpose— Cerebral artery aneurysms rupture when wall tension exceeds the strength of the wall tissue. At present, risk-assessment of unruptured aneurysms does not include evaluation of the lesions shape, yet clinical experience suggests that this is of importance. We aimed to develop a computational model for simulation of fluid-structure interaction in cerebral aneurysms based on patient specific lesion geometry, with special emphasis on wall tension. Methods— An advanced isogeometric fluid-structure analysis model incorporating flexible aneurysm wall based on patient specific computed tomography angiogram images was developed. Variables used in the simulation model were retrieved from a literature review. Results— The simulation results exposed areas of high wall tension and wall displacement located where aneurysms usually rupture. Conclusion— We suggest that analyzing wall tension and wall displacement in cerebral aneurysms by numeric simulation could be developed into a novel method for individualized prediction of rupture risk.

ALE-VMS and ST-VMS methods for computer modeling of wind-turbine rotor aerodynamics and fluid-structure interaction

We provide an overview of the Arbitrary Lagrangian–Eulerian Variational Multiscale (ALE-VMS) and Space–Time Variational Multiscale (ST-VMS) methods we have developed for computer modeling of wind-turbine rotor aerodynamics and fluid–structure interaction (FSI). The related techniques described include weak enforcement of the essential boundary conditions, Kirchhoff–Love shell modeling of the rotor-blade structure, NURBS-based isogeometric analysis, and full FSI coupling. We present results from application of these methods to computer modeling of NREL 5MW and NREL Phase VI wind-turbine rotors at full scale, including comparison with experimental data.

Isogeometric analysis of free-surface flow

Abstract This paper presents the first application of isogeometric analysis, a new computational technology built on higher-order and higher-continuity basis functions employed in Computer-Aided Design and computer graphics, to the computation of free-surface phenomena described using the level set approach. The method is based on the variational framework that is suitable for discretization by standard finite elements as well as the basis functions employed in isogeometric analysis. The underlying numerical formulation globally conserves mass and preserves a sharp air–water interface for the entire length of the simulation. The numerical tests indicate that the proposed methodology gives an accurate description of the free-surface behavior in both quasi-steady and dynamic regimes. Furthermore, very good per-degree-of-freedom accuracy is obtained when higher-order and higher-continuity isogeometric discretizations are employed in free-surface computations.

CHALLENGES AND DIRECTIONS IN COMPUTATIONAL FLUID–STRUCTURE INTERACTION

In this lead paper of the special issue, we provide some comments on challenges and directions in computational fluid–structure interaction (FSI). We briefly discuss the significance of computational FSI methods, their components, moving-mesh and nonmoving-mesh methods, mesh moving and remeshing concepts, and FSI coupling techniques.

Free-Surface Flow and Fluid-Object Interaction Modeling With Emphasis on Ship Hydrodynamics

Abstract : This paper presents our approach for the computation of free-surface/rigid-body interaction phenomena with emphasis on ship hydrodynamics. We adopt the level set approach to capture the free-surface. The rigid body is described using six-degree-of-freedom equations of motion. An interface-tracking method is used to handle the interface between the moving rigid body and the fluid domain. An Arbitrary Lagrangian Eulerian version of the residual-based variational multiscale formulation for the Navier Stokes and level set equations is employed in order to accommodate the fluid domain motion. The free-surface/rigid body problem is formulated and solved in a fully coupled fashion. The numerical results illustrate the accuracy and robustness of the proposed approach.

STRUCTURAL MECHANICS MODELING AND FSI SIMULATION OF WIND TURBINES

A fluid–structure interaction (FSI) validation study of the Micon 65/13M wind turbine with Sandia CX-100 composite blades is presented. A rotation-free isogeometric shell formulation is used to model the blade structure, while the aerodynamics formulation makes use of the FEM-based ALE-VMS method. The structural mechanics formulation is validated by means of eigenfrequency analysis of the CX-100 blade. For the coupling between the fluid and structural mechanics domains, a nonmatching discretization approach is adopted. The simulations are done at realistic wind conditions and rotor speeds. The rotor-tower interaction that influences the aerodynamic torque is captured. The computed aerodynamic torque generated by the Micon 65/13M wind turbine compares well with that obtained from on-land experimental tests.

A conservative level set method suitable for variable-order approximations and unstructured meshes

This paper presents a formulation for free-surface computations capable of handling complex phenomena, such as wave breaking, without excessive mass loss or smearing of the interface. The formulation is suitable for discretizations using finite elements of any topology and order, or other approaches such as isogeometric and finite volume methods. Furthermore, the approach builds on standard level set tools and can therefore be used to augment existing implementations of level set methods with discrete conservation properties. Implementations of the method are tested on several difficult two- and three-dimensional problems, including two incompressible air/water flow problems with available experimental results. Linear and quadratic approximations on unstructured tetrahedral and trilinear approximations on hexahedral meshes were tested. Global conservation and agreement with experiments as well as computations by other researchers are obtained.

Fluid–Structure Interaction Modeling of Vertical-Axis Wind Turbines

Full-scale, 3D, time-dependent aerodynamics and fluid‐structure interaction (FSI) simulations of a Darrieus-type vertical-axis wind turbine (VAWT) are presented. A structural model of the Windspire VAWT (Windspire energy, http://www.windspireenergy.com/ )i s developed, which makes use of the recently proposed rotation-free Kirchhoff‐Love shell and beam/cable formulations. A moving-domain finite-element-based ALE-VMS (arbitrary Lagrangian‐Eulerian-variational-multiscale) formulation is employed for the aerodynamics in combination with the sliding-interface formulation to handle the VAWT mechanical components in relative motion. The sliding-interface formulation is augmented to handle nonstationary cylindrical sliding interfaces, which are needed for the FSI modeling of VAWTs. The computational results presented show good agreement with the field-test data. Additionally, several scenarios are considered to investigate the transient VAWT response and the issues related to self-starting. [DOI: 10.1115/1.4027466]

Aerodynamic Simulation of Vertical-Axis Wind Turbines

Full-scale, 3D, time-dependent aerodynamics modeling and simulation of a Darrieus-type vertical-axis wind turbine (VAWT) is presented. The simulations are performed using a moving-domain finite-element-based ALE-VMS technique augmented with a sliding-interface formulation to handle the rotor-stator interactions present. We simulate a single VAWT using a sequence of meshes with increased resolution to assess the computational requirements for this class of problems. The computational results are in good agreement with experimental data. We also perform a computation of two side-by-side counterrotating VAWTs to illustrate how the ALE-VMS technique may be used for the simulation of multiple turbines placed in arrays.