We extend the recently introduced divergence-conforming immersed boundary (DCIB) method [1] to fluid-structure interaction (FSI) problems involving closed co-dimension one solids. We focus on capsules and vesicles, whose discretization is particularly challenging due to the higher-order derivatives that appear in their formulations. In two-dimensional settings, we employ cubic B-splines with periodic knot vectors to obtain discretizations of closed curves with C² inter-element continuity. In three-dimensional settings, we use analysis-suitable bi-cubic T-splines to obtain discretizations of closed surfaces with at least C¹ inter-element continuity. Large spurious changes of the fluid volume inside closed co-dimension one solids are a well-known issue for IB methods. The DCIB method results in volume changes orders of magnitude lower than conventional IB methods. This is a byproduct of discretizing the velocity-pressure pair with divergence-conforming B-splines, which lead to negligible incompressibility errors at the Eulerian level. The higher inter-element continuity of divergence-conforming B-splines is also crucial to avoid the quadrature/interpolation errors of IB methods becoming the dominant discretization error. Benchmark and application problems of vesicle and capsule dynamics are solved, including mesh-independence studies and comparisons with other numerical methods.

@en

Spaces of discrete differential forms can be applied to numerically solve the partial differential equations that govern phenomena such as electromagnetics and fluid mechanics. Robustness of the resulting numerical methods is complemented by pointwise satisfaction of conservation laws (e.g., mass conservation) in the discrete setting. Here we present the construction of isogeometric discrete differential forms, i.e., differential form spaces built using smooth splines. We first present an algorithm for computing Bézier extraction operators for univariate spline differential forms that allow local degree elevation. Then, using tensor-products of the univariate splines, a complex of discrete differential forms is built on meshes that contain polar singularities, i.e., edges that are singularly mapped onto points. We prove that the spline complexes share the same cohomological structure as the de Rham complex. Several examples are presented to demonstrate the applicability of the proposed methodology. In particular, the splines spaces derived are used to simulate generalized Stokes flow on arbitrarily curved smooth surfaces and to numerically demonstrate (a) optimal approximation and inf–sup stability of the spline spaces; (b) pointwise incompressible flows; and (c) flows on deforming surfaces.@en

Polynomial splines are ubiquitous in the fields of computer-aided geometric design and computational analysis. Splines on T-meshes, especially, have the potential to be incredibly versatile since local mesh adaptivity enables efficient modeling and approximation of local features. Meaningful use of such splines for modeling and approximation requires the construction of a suitable spanning set of linearly independent splines, and a theoretical understanding of the spline space dimension can be a useful tool when assessing possible approaches for building such splines. Here, we provide such a tool. Focusing on T-meshes, we study the dimension of the space of bivariate polynomial splines, and we discuss the general setting where local mesh adaptivity is combined with local polynomial degree adaptivity. The latter allows for the flexibility of choosing non-uniform bi-degrees for the splines, i.e., different bi-degrees on different faces of the T-mesh. In particular, approaching the problem using tools from homological algebra, we generalize the framework and the discourse presented by Mourrain (Math. Comput. 83(286):847–871, 2014) for uniform bi-degree splines. We derive combinatorial lower and upper bounds on the spline space dimension and subsequently outline sufficient conditions for the bounds to coincide.@en

With the aim of a seamless integration with Computer-Aided Design, Isogeometric Analysis has been proposed by Hughes et al. (2005) as a numerical technique for the solution of partial differential equations. Indeed, isogeometric analysis is based on splines, the same functions that are adopted for geometry parametrizations in CAD. Smooth splines yield two important benefits when compared to C⁰ piecewise polynomial approximations: superior accuracy and stability, and the possibility to directly discretize high-order differential equations, such as the ones arising in thin-shell theory, in fracture models, in phase-field based multiphase flows, and in geometric flows on surfaces. In this chapter we review three different methods to construct C¹ isogeometric spaces on multi-patch domains or unstructured quadrilateral meshes. The first, proposed in Collin et al. (2016), is based on the concept of geometric continuity, well-known in geometric design. The second, from Toshniwal et al. (2017c), uses a specific singular construction (the D-patch construction) at extraordinary points. The third is a polar construction from Toshniwal et al. (2017a). Such constructions possess properties that make them suitable for performing both computational analysis and geometric modeling.

@en

In this paper, we present an efficient and robust approach to compute a normalized B-spline-like basis for spline spaces with pieces drawn from extended Tchebycheff spaces. The extended Tchebycheff spaces and their dimensions are allowed to change from interval to interval. The approach works by constructing a matrix that maps a generalized Bernstein-like basis to the B-spline-like basis of interest. The B-spline-like basis shares many characterizing properties with classical univariate B-splines and may easily be incorporated in existing spline codes. This may contribute to the full exploitation of Tchebycheffian splines in applications, freeing them from the restricted role of an elegant theoretical extension of polynomial splines. Numerical examples are provided that illustrate the procedure described.

@en

Analysis-suitable T-splines (ASTS) including both extraordinary points and T-junctions are used to solve Kirchhoff–Love shell problems. Extraordinary points are required to represent surfaces with arbitrary topological genus. T-junctions enable local refinement of regions where increased resolution is needed. The benefits of using ASTS to define shell geometries are at least two-fold: (1) The manual and time-consuming task of building a new mesh from scratch using the CAD geometry as an input is avoided and (2) C¹ or higher inter-element continuity enables the discretization of shell formulations in primal form defined by fourth-order partial differential equations. A complete and state-of-the-art description of the development of ASTS, including extraordinary points and T-junctions, is presented. In particular, we improve the construction of C¹-continuous non-negative spline basis functions near extraordinary points to obtain optimal convergence rates with respect to the square root of the number of degrees of freedom when solving linear elliptic problems. The applicability of the proposed technology to shell analysis is exemplified by performing geometrically nonlinear Kirchhoff–Love shell simulations of a pinched hemisphere, an oil sump of a car, a pipe junction, and a B-pillar of a car with 15 holes. Building ASTS for these examples involves using T-junctions and extraordinary points with valences 3, 5, and 6, which often suffice for the design of free-form surfaces. Our analysis results are compared with data from the literature using either a seven-parameter shell formulation or Kirchhoff–Love shells. We have also imported both finite element meshes and ASTS meshes into the commercial software LS-DYNA, used Reissner–Mindlin shells, and compared the result with our Kirchhoff–Love shell results. Excellent agreement is found in all cases. The complexity of the shell geometries considered in this paper shows that ASTS are applicable to real-world industrial problems.

@en