Weakly Enforced Boundary Conditions for the NURBS-Based Finite Cell Method

In this paper, we present a variationally consistent formulation for the weak enforcement of essential boundary conditions as an extension to the finite cell method, a fictitious domain method of higher order. The absence of boundary fitted elements in fictitious domain or immersed boundary methods significantly restricts a strong enforcement of essential boundary conditions to models where the boundary of the solution domain coincides with the embedding analysis domain. Penalty methods and Lagrange multiplier methods are adequate means to overcome this limitation but often suffer from various drawbacks with severe consequences for a stable and accurate solution of the governing system of equations. In this contribution, we follow the idea of NITSCHE [29] who developed a stable scheme for the solution of the Laplace problem taking weak boundary conditions into account. An extension to problems from linear elasticity shows an appropriate behavior with regard to numerical stability, accuracy and an adequate convergence behavior. NURBS are chosen as a high-order approximation basis to benefit from their smoothness and flexibility in the process of uniform model refinement

A generalized finite element formulation for arbitrary basis functions : from isogeometric analysis to XFEM

Many of the formulations of cm-rent research interest, including iosogeometric methods and the extended finite element method, use nontraditional basis functions. Some, such as subdivision surfaces, may not have convenient analytical representations. The concept of an element, if appropriate at all, no longer coincides with the traditional definition. Developing a new software for each new class of basis functions is a large research burden, especially, if the problems involve large deformations, non-linear materials, and contact. The objective of this paper is to present a method that separates as much as possible the generation and evaluation of the basis functions from the analysis, resulting in a formulation that can be implemented within the traditional structure of a finite clement program but that permits the use of arbitrary sets of basis functions that are defined only through the input file. Elements ranging from a traditional linear four-node tetrahedron through a higher-order element combining XFEM and isogeometric analysis may be specified entirely through an input file without any additional programming. Examples of this framework to applications with Lagrange elements, isogeometric elements, and XFEM basis functions for fracture are presented

X-FEM in isogeometric analysis for linear fracture mechanics

The extended finite element method (X-FEM) has proven to be an accurate, robust method for solving problems in fracture mechanics. X-FEM has typically been used with elements using linear basis functions, although some work has been performed using quadratics. In the current work, the X-FEM formulation is incorporated into isogeometric analysis to obtain solutions with higher order convergence rates for problems in linear fracture mechanics. In comparison with X-FEM with conventional finite elements of equal degree, the NURBS-based isogeometric analysis gives equal asymptotic convergence rates and equal accuracy with fewer degrees of freedom (DOF). Results for linear through quartic NURBS basis functions are presented for a multiplicity of one or a multiplicity equal the degree

Fluid-structure interaction in blood flow capturing non-zero longitudinal structure displacement

We present a new model and a novel loosely coupled partitioned numerical scheme modeling fluid-structure interaction (FSI) in blood flow allowing non-zero longitudinal displacement. Arterial walls are modeled by a {linearly viscoelastic, cylindrical Koiter shell model capturing both radial and longitudinal displacement}. Fluid flow is modeled by the Navier-Stokes equations for an incompressible, viscous fluid. The two are fully coupled via kinematic and dynamic coupling conditions. Our numerical scheme is based on a new modified Lie operator splitting that decouples the fluid and structure sub-problems in a way that leads to a loosely coupled scheme which is {unconditionally} stable. This was achieved by a clever use of the kinematic coupling condition at the fluid and structure sub-problems, leading to an implicit coupling between the fluid and structure velocities. The proposed scheme is a modification of the recently introduced "kinematically coupled scheme" for which the newly proposed modified Lie splitting significantly increases the accuracy. The performance and accuracy of the scheme were studied on a couple of instructive examples including a comparison with a monolithic scheme. It was shown that the accuracy of our scheme was comparable to that of the monolithic scheme, while our scheme retains all the main advantages of partitioned schemes, such as modularity, simple implementation, and low computational costs

The Influence of Quadrature Errors on Isogeometric Mortar Methods

Mortar methods have recently been shown to be well suited for isogeometric analysis. We review the recent mathematical analysis and then investigate the variational crime introduced by quadrature formulas for the coupling integrals. Motivated by finite element observations, we consider a quadrature rule purely based on the slave mesh as well as a method using quadrature rules based on the slave mesh and on the master mesh, resulting in a non-symmetric saddle point problem. While in the first case reduced convergence rates can be observed, in the second case the influence of the variational crime is less significant

Proper Orthogonal Decomposition Closure Models For Turbulent Flows: A Numerical Comparison

This paper puts forth two new closure models for the proper orthogonal decomposition reduced-order modeling of structurally dominated turbulent flows: the dynamic subgrid-scale model and the variational multiscale model. These models, which are considered state-of-the-art in large eddy simulation, together with the mixing length and the Smagorinsky closure models, are tested in the numerical simulation of a 3D turbulent flow around a circular cylinder at Re = 1,000. Two criteria are used in judging the performance of the proper orthogonal decomposition reduced-order models: the kinetic energy spectrum and the time evolution of the POD coefficients. All the numerical results are benchmarked against a direct numerical simulation. Based on these numerical results, we conclude that the dynamic subgrid-scale and the variational multiscale models perform best.Comment: 28 pages, 6 figure

A full Eulerian finite difference approach for solving fluid-structure coupling problems

A new simulation method for solving fluid-structure coupling problems has been developed. All the basic equations are numerically solved on a fixed Cartesian grid using a finite difference scheme. A volume-of-fluid formulation (Hirt and Nichols (1981, J. Comput. Phys., 39, 201)), which has been widely used for multiphase flow simulations, is applied to describing the multi-component geometry. The temporal change in the solid deformation is described in the Eulerian frame by updating a left Cauchy-Green deformation tensor, which is used to express constitutive equations for nonlinear Mooney-Rivlin materials. In this paper, various verifications and validations of the present full Eulerian method, which solves the fluid and solid motions on a fixed grid, are demonstrated, and the numerical accuracy involved in the fluid-structure coupling problems is examined.Comment: 38 pages, 27 figures, accepted for publication in J. Comput. Phy

Variational multiscale approximation of the one-dimensional forced Burgers equation: the role of orthogonal subgrid scales in turbulence modeling

This is the accepted version of the following article: [Bayona C, Baiges J, Codina R. Variational multi-scale approximation of the one-dimensional forced Burgers equation: the role of orthogonal subgrid scales in turbulence modeling. Int J Numer Meth Fluids. 2018;86:313–328. https://doi.org/10.1002/fld.4420], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/fld.4420/fullA numerical approximation for the one-dimensional Burgers equation is proposed by means of the orthogonal subgrid scales–variational multiscale (OSGS-VMS) method. We evaluate the role of the variational subscales in describing the Burgers “turbulence” phenomena. Particularly, we seek to clarify the interaction between the subscales and the resolved scales when the former are defined to be orthogonal to the finite-dimensional space. Direct numerical simulation is used to evaluate the resulting OSGS-VMS energy spectra. The comparison against a large eddy simulation model is presented for numerical discretizations in which the grid is not capable of solving the small scales. An accurate approximation to the phenomena of turbulence is obtained with the addition of the purely dissipative numerical terms given by the OSGS-VMS method without any modification of the continuous problem.Peer ReviewedPostprint (author's final draft

The role of the Bézier extraction operator for T-splines of arbitrary degree: linear dependencies, partition of unity property, nesting behaviour and local refinement

We determine linear dependencies and the partition of unity property of T-spline meshes of arbitrary degree using the Bézier extraction operator. Local refinement strategies for standard, semi-standard and nonstandard T-splines – also by making use of the Bézier extraction operator – are presented for meshes of even and odd polynomial degrees. A technique is presented to determine the nesting between two T-spline meshes, again exploiting the Bézier extraction operator. Finally, the hierarchical refinement of standard, semi-standard and non-standard T-spline meshes is discussed. This technique utilises the reconstruction operator, which is the inverse of the Bézier extraction operator

Geometry and Mechanics

IASS-IACM 2008 Session: Geometry and Mechanics -- Session Organizers: Kai-Uwe BLETZINGER (TU Munich), Fehmi CIRAK (University of Cambridge) -- Keynote Lecture: "Modeling and computation of patient-specific vascular fluid-structure interaction using Isogeometric Analysis" by Yuri BAZILEVS , Victor M. CALO, Thomas J. R. HUGHES (University of Texas at Austin), Yongie ZHANG (Carnegie Mellon University) -- Keynote Lecture: "Optimal shapes of mechanically motivated surfaces" by Kai-Uwe BLETZINGER , Matthias FIRL, Johannes LINHARD, Roland WUCHNER (TU Munich) -- "Subdivision shells for nonsmooth and branching geometries" by Quan LONG, Fehmi CIRAK (University of Cambridge) -- "Water landing analyses with explicit finite element method" by John T. WANG (NASA Langley Research Center) -- "On a geometrically exact contact description for shells: From linear approximations for shells to high-order FEM" by Alexander KONYUKHOV, Karl SCHWEIZERHOF (University of Karlsruhe