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.
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.
•Div-conforming B-splines improve immersed fluid–structure interaction (FSI) analysis.•Strong mass conservation prevents spurious non-physical leakage through barriers.•Semi-implicit time integration is shown to converge a priori for a model problem.•Div-conforming immersogeometric FSI analysis is practical for heart valve analysis.•FSI simulations reproduce qualitative features of in vitro experiments.
We introduce provably unconditionally stable mixed variational methods for phase-field models. Our formulation is based on a mixed finite element method for space discretization and a new ...second-order accurate time integration algorithm. The fully-discrete formulation inherits the main characteristics of conserved phase dynamics, namely, mass conservation and nonlinear stability with respect to the free energy. We illustrate the theory with the Cahn–Hilliard equation, but our method may be applied to other phase-field models. We also propose an adaptive time-stepping version of the new time integration method. We present some numerical examples that show the accuracy, stability and robustness of the new method.
We present an early version of a Susceptible–Exposed–Infected–Recovered–Deceased (SEIRD) mathematical model based on partial differential equations coupled with a heterogeneous diffusion model. The ...model describes the spatio-temporal spread of the COVID-19 pandemic, and aims to capture dynamics also based on human habits and geographical features. To test the model, we compare the outputs generated by a finite-element solver with measured data over the Italian region of Lombardy, which has been heavily impacted by this crisis between February and April 2020. Our results show a strong qualitative agreement between the simulated forecast of the spatio-temporal COVID-19 spread in Lombardy and epidemiological data collected at the municipality level. Additional simulations exploring alternative scenarios for the relaxation of lockdown restrictions suggest that reopening strategies should account for local population densities and the specific dynamics of the contagion. Thus, we argue that data-driven simulations of our model could ultimately inform health authorities to design effective pandemic-arresting measures and anticipate the geographical allocation of crucial medical resources.
Following a series of recent innovations, isogeometric shell analysis based on trimmed CAD surfaces is currently being developed into an accurate, efficient and mature design-through-analysis ...methodology. This work contributes to this emerging technology with respect to the following aspects. On the analysis side, we present a robust variationally consistent Nitsche-type formulation for thin shells at large deformations that weakly enforces coupling constraints at trimming curves. On the geometry side, we present a set of algorithms that enable automatic interaction of trimmed shell analysis with CAD data structures based on the STEP exchange format. We integrate these methodologies in a comprehensive framework for isogeometric trimmed shell analysis. We demonstrate that our framework is able to seamlessly perform large-deformation stress analysis of an industry-scale 76-patch surface model of a Dodge RAM hood, while delivering comparable accuracy with respect to Simulia’s commercial software package Abaqus.
•We present a variationally consistent formulation for coupling isogeometric shells at trimming curves.•We also present algorithms that enable automatic interaction with CAD data structures based on the STEP exchange format.•We integrate these methodologies in a comprehensive framework for isogeometric trimmed shell analysis.•It enables seamless and accurate large-deformation stress analysis, illustrated with a 76-patch model of a Dodge RAM hood.
We propose a framework that combines variational immersed-boundary and arbitrary Lagrangian–Eulerian methods for fluid–structure interaction (FSI) simulation of a bioprosthetic heart valve implanted ...in an artery that is allowed to deform in the model. We find that the variational immersed-boundary method for FSI remains robust and effective for heart valve analysis when the background fluid mesh undergoes deformations corresponding to the expansion and contraction of the elastic artery. Furthermore, the computations presented in this work show that the arterial wall deformation contributes significantly to the realism of the simulation results, leading to flow rates and valve motions that more closely resemble those observed in practice.
The Cahn–Hilliard equation involves fourth-order spatial derivatives. Finite element solutions are not common because primal variational formulations of fourth-order operators are only well defined ...and integrable if the finite element basis functions are piecewise smooth and globally
C
1
-continuous. There are a very limited number of two-dimensional finite elements possessing
C
1
-continuity applicable to complex geometries, but none in three-dimensions. We propose isogeometric analysis as a technology that possesses a unique combination of attributes for complex problems involving higher-order differential operators, namely, higher-order accuracy, robustness, two- and three-dimensional geometric flexibility, compact support, and, most importantly, the possibility of
C
1
and higher-order continuity. A NURBS-based variational formulation for the Cahn–Hilliard equation was tested on two- and three-dimensional problems. We present steady state solutions in two-dimensions and, for the first time, in three-dimensions. To achieve these results an adaptive time-stepping method is introduced. We also present a technique for desensitizing calculations to dependence on mesh refinement. This enables the calculation of topologically correct solutions on coarse meshes, opening the way to practical engineering applications of phase-field methodology.
We present a framework for geometric design and isogeometric analysis on unstructured quadrilateral meshes. Acknowledging the differing requirements posed by design (e.g., the convenience of an ...intuitive control net) and analysis (e.g., good approximation behavior), we propose the construction of a separate, smooth spline space for each while ensuring isogeometric compatibility – requiring the geometric models to be members of the analysis-suitable spaces. The methodology is simple and is presented for bi-cubic splines; extensions to higher degrees are possible, and are briefly discussed. The presentation has been structured to show compatibility with T-splines – a state-of-the-art CAD technology – but the approach should extend to other locally refinable spline technologies (based on local tensor-product structures). An instantiation of the framework is presented, and several numerical tests focused on geometric design and isogeometric analysis demonstrate the versatility of the developed framework, and show significantly higher convergence rates than attained previously in the considered setting.
•We present a framework for building smooth splines on meshes with extraordinary points.•The spline spaces possess several desirable properties for both CAD and IGA.•Vertex-based, smooth, linearly independent splines are used for modeling geometries.•Compatible design and analysis spaces imply exact satisfaction of all patch tests.•Optimal or almost-optimal convergence rates are achieved in typical analysis situations.