Moving Domain Computational Fluid Dynamics to Interface with an Embryonic Model of Cardiac Morphogenesis

Peristaltic contraction of the embryonic heart tube produces time- and spatial-varying wall shear stress (WSS) and pressure gradients (∇P) across the atrioventricular (AV) canal. Zebrafish (Danio rerio) are a genetically tractable system to investigate cardiac morphogenesis. The use of Tg(fli1a:EGFP)y1 transgenic embryos allowed for delineation and two-dimensional reconstruction of the endocardium. This time-varying wall motion was then prescribed in a two-dimensional moving domain computational fluid dynamics (CFD) model, providing new insights into spatial and temporal variations in WSS and ∇P during cardiac development. The CFD simulations were validated with particle image velocimetry (PIV) across the atrioventricular (AV) canal, revealing an increase in both velocities and heart rates, but a decrease in the duration of atrial systole from early to later stages. At 20-30 hours post fertilization (hpf), simulation results revealed bidirectional WSS across the AV canal in the heart tube in response to peristaltic motion of the wall. At 40-50 hpf, the tube structure undergoes cardiac looping, accompanied by a nearly 3-fold increase in WSS magnitude. At 110-120 hpf, distinct AV valve, atrium, ventricle, and bulbus arteriosus form, accompanied by incremental increases in both WSS magnitude and ∇P, but a decrease in bi-directional flow. Laminar flow develops across the AV canal at 20-30 hpf, and persists at 110-120 hpf. Reynolds numbers at the AV canal increase from 0.07±0.03 at 20-30 hpf to 0.23±0.07 at 110-120 hpf (p< 0.05, n=6), whereas Womersley numbers remain relatively unchanged from 0.11 to 0.13. Our moving domain simulations highlights hemodynamic changes in relation to cardiac morphogenesis; thereby, providing a 2-D quantitative approach to complement imaging analysis. © 2013 Lee et al

Comparative patient-specific fsi modeling of cerebral aneurysms

We consider a total of ten cases, at three different locations, half of which ruptured sometime after the images were taken. We use the stabilized space–time FSI technique developed by the Team for Advanced Flow Simulation and Modeling, together with a number of special techniques targeting arterial FSI modeling. We compare the ten cases based on the wall shear stress, oscillatory shear index, and the arterial-wall stress. We also investigate how simpler approaches to computer modeling of cerebral aneurysms perform compared to FSI modeling

Computational analysis of performance deterioration of a wind turbine blade strip subjected to environmental erosion

Wind-turbine blade rain and sand erosion, over long periods of time, can degrade the aerodynamic performance and therefore the power production. Computational analysis of the erosion can help engineers have a better understanding of the maintenance and protection requirements. We present an integrated method for this class of computational analysis. The main components of the method are the streamline-upwind/Petrov–Galerkin (SUPG) and pressure-stabilizing/Petrov–Galerkin (PSPG) stabilizations, a finite element particle-cloud tracking method, an erosion model based on two time scales, and the solid-extension mesh moving technique (SEMMT). The turbulent-flow nature of the analysis is handled with a Reynolds-averaged Navier–Stokes model and SUPG/PSPG stabilization, the particle-cloud trajectories are calculated based on the computed flow field and closure models defined for the turbulent dispersion of particles, and one-way dependence is assumed between the flow and particle dynamics. Because the geometry update due to the erosion has a very long time scale compared to the fluid–particle dynamics, the update takes place in a sequence of “evolution steps” representing the impact of the erosion. A scale-up factor, calculated in different ways depending on the update threshold criterion, relates the erosions and particle counts in the evolution steps to those in the fluid–particle simulation. As the blade geometry evolves, the mesh is updated with the SEMMT. We present computational analysis of rain and sand erosion for a wind-turbine blade strip, including a case with actual rainfall data and experimental aerodynamic data for eroded airfoil geometries

Parallel finite element simulation of large ram-air parachutes

In the near future, large ram-air parachutes are expected to provide the capability of delivering 21 ton payloads from altitudes as high as 25,000 ft. In development and test and evaluation of these parachutes the size of the parachute needed and the deployment stages involved make high-performance computing (HPC) simulations a desirable alternative to costly airdrop tests. Although computational simulations based on realistic, 3D, time-dependent models will continue to be a major computational challenge, advanced finite element simulation techniques recently developed for this purpose and the execution of these techniques on HPC platforms are significant steps in the direction to meet this challenge. In this paper, two approaches for analysis of the inflation and gliding of ram-air parachutes are presented. In one of the approaches the point mass flight mechanics equations are solved with the time-varying drag and lift areas obtained from empirical data. This approach is limited to parachutes with similar configurations to those for which data are available. The other approach is 3D finite element computations based on the Navier-Stokes equations governing the airflow around the parachute canopy and Newton's law of motion governing the 3D dynamics of the canopy, with the forces acting on the canopy calculated from the simulated flow field. At the earlier stages of canopy inflation the parachute is modelled as an expanding box, whereas at the later stages, as it expands, the box transforms to a parafoil and glides. These finite element computations are carried out on the massively parallel supercomputers CRAY T3D and Thinking Machines CM-5, typically with millions of coupled, non-linear finite element equations solved simultaneously at every time step or pseudo-time step of the simulation

Physics based GMRES preconditioner for compressible and incompressible Navier-Stokes equations

This paper presents the implementation of a local physics preconditioning mass matrix [8] for an unified approach of 3D compressible and incompressible Navier-Stokes equations using an SUPG finite element formulation and GMRES implicit solver. During the last years a lot of effort has been dedicated to finding a&nbsp;unified approach&nbsp;for compressible and incompressible flow in order to treat fluid dynamic problems with a very wide range of Mach and Reynolds numbers [10,26,37]. On the other hand, SUPG finite element formulation and GMRES implicit solver is one of the most robust combinations to solve state of the art CFD problems [1,6,9,22,29,30,31]. The selection of a good&nbsp;preconditioner&nbsp;and its performance on parallel architecture is another open problem in CFD community. The&nbsp;local&nbsp;feature of the preconditioner presented here means that no communication among processors is needed when working on parallel architectures. Due to these facts we consider that this research can make some contributions towards the development of a unified fluid dynamic model with high rates of convergence for any combination of Mach and Reynolds numbers, being very suitable for massively parallel computations. Finally, it is important to remark that while this kind of preconditioning produces stabilized results in nearly incompressible regimes the standard version exhibits some numerical drawbacks that lead to solutions without physical meaning

Massively parallel finite element simulation of compressible and incompressible flows

We present a review of where our research group stands in parallel finite element simulation of flow problems on the Connection Machines, an effort that started for our group in the fourth quarter of 1991. This review includes an overview of our work on computation of flow problems involving moving boundaries and interfaces, such as free surfaces, two-liquid interfaces, and fluid-structure and fluid-particle interactions. With numerous examples, we demonstrate that, with these new computational capabilities, today we are at a point where we routinely solve practical flow problems, including those in 3D and those involving moving boundaries and interfaces. We solve these problems with unstructured grids and implicit methods, with some of the problem sizes exceeding 5 000 000 equations, and with computational speeds up to two orders of magnitude higher than what was previously available to us on the traditional vector supercomputers

Dynamic and fluid–structure interaction simulations of bioprosthetic heart valves using parametric design with T-splines and Fung-type material models

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 Wind-Turbine Analysis with the ALE-VMS and ST-VMS Methods

We provide an overview of the aerodynamic and FSI analysis of wind turbines the first three authors’ teams carried out in recent years with the ALE-VMS and ST-VMS methods. The ALE-VMS method is the variational multiscale version of the Arbitrary Lagrangian–Eulerian (ALE) method. The VMS components are from the residual-based VMS (RBVMS) method. The ST-VMS method is the VMS version of the Deforming-Spatial-Domain/Stabilized Space–Time (DSD/SST) method. The techniques complementing these core methods include weak enforcement of the essential boundary conditions, NURBS-based isogeometric analysis, using NURBS basis functions in temporal representation of the rotor motion, mesh motion and also in remeshing, rotation representation with constant angular velocity, Kirchhoff–Love shell modeling of the rotor-blade structure, and full FSI coupling. The analysis cases include the aerodynamics of wind-turbine rotor and tower and the FSI that accounts for the deformation of the rotor blades. The specific wind turbines considered are NREL 5MW, NREL Phase VI and Micon 65/13M, all at full scale, and our analysis for NREL Phase VI and Micon 65/13M includes comparison with the experimental data