In recent years, isogeometric analysis (IGA) has received great attention in many fields of computational mechanics research. Especially for computational contact mechanics, an exact and smooth ...surface representation is highly desirable. As a consequence, many well-known finite element methods and algorithms for contact mechanics have been transferred to IGA. In the present contribution, the so-called dual mortar method is investigated for both contact mechanics and classical domain decomposition using NURBS basis functions. In contrast to standard mortar methods, the use of dual basis functions for the Lagrange multiplier based on the mathematical concept of biorthogonality enables an easy elimination of the additional Lagrange multiplier degrees of freedom from the global system. This condensed system is smaller in size, and no longer of saddle point type but positive definite. A very simple and commonly used element-wise construction of the dual basis functions is directly transferred to the IGA case. The resulting Lagrange multiplier interpolation satisfies discrete inf–sup stability and biorthogonality, however, the reproduction order is limited to one. In the domain decomposition case, this results in a limitation of the spatial convergence order to O(h32) in the energy norm, whereas for unilateral contact, due to the lower regularity of the solution, optimal convergence rates are still met. Numerical examples are presented that illustrate these theoretical considerations on convergence rates and compare the newly developed isogeometric dual mortar contact formulation with its standard mortar counterpart as well as classical finite elements based on first and second order Lagrange polynomials.
•A dual mortar method for NURBS-based isogeometric analysis is developed.•Spatial convergence orders are analyzed for mesh tying and contact mechanics.•A lack of reproduction properties limits the convergence in mesh tying applications.•Optimal convergence results are achieved for contact applications.•The higher smoothness of NURBS delivers smoother results for contact forces.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
For the simulation of fractured porous media, a common approach is the use of co-dimension one models for the fracture description. In order to simulate correctly the behavior at fracture crossings, ...standard models are not sufficient because they either cannot capture all important flow processes or are computationally inefficient. We propose a new concept to simulate co-dimension one fracture crossings and show its necessity and accuracy by means of an example and a comparison to a literature benchmark. From the application point of view, often the pressure is known only at a limited number of discrete points and an interpolation is used to define the boundary condition at the remaining parts of the boundary. The quality of the interpolation, especially in fracture models, influences the global solution significantly. We propose a new method to interpolate boundary conditions on boundaries that are intersected by fractures and show the advantages over standard interpolation methods.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
The mortar finite element method allows the coupling of different discretization schemes and triangulations across subregion boundaries. In the original mortar approach the matching at the interface ...is realized by enforcing an orthogonality relation between the jump and a modified trace space which serves as a space of Lagrange multipliers. In this paper, this Lagrange multiplier space is replaced by a dual space without losing the optimality of the method. The advantage of this new approach is that the matching condition is much easier to realize. In particular, all the basis functions of the new method are supported in a few elements. The mortar map can be represented by a diagonal matrix; in the standard mortar method a linear system of equations must be solved. The problem is considered in a positive definite nonconforming variational as well as an equivalent saddle-point formulation.
Full text
Available for:
BFBNIB, CEKLJ, DOBA, INZLJ, IZUM, KILJ, NMLJ, NUK, PILJ, PNG, SAZU, UILJ, UKNU, UL, UM, UPUK, ZRSKP
In this work, we introduce an algorithmic approach to generate microvascular networks starting from larger vessels that can be reconstructed without noticeable segmentation errors. Contrary to larger ...vessels, the reconstruction of fine‐scale components of microvascular networks shows significant segmentation errors, and an accurate mapping is time and cost intense. Thus there is a need for fast and reliable reconstruction algorithms yielding surrogate networks having similar stochastic properties as the original ones. The microvascular networks are constructed in a marching way by adding vessels to the outlets of the vascular tree from the previous step. To optimise the structure of the vascular trees, we use Murray's law to determine the radii of the vessels and bifurcation angles. In each step, we compute the local gradient of the partial pressure of oxygen and adapt the orientation of the new vessels to this gradient. At the same time, we use the partial pressure of oxygen to check whether the considered tissue block is supplied sufficiently with oxygen. Computing the partial pressure of oxygen, we use a 3D‐1D coupled model for blood flow and oxygen transport. To decrease the complexity of a fully coupled 3D model, we reduce the blood vessel network to a 1D graph structure and use a bi‐directional coupling with the tissue which is described by a 3D homogeneous porous medium. The resulting surrogate networks are analysed with respect to morphological and physiological aspects.
In this article, we present an algorithm to generate microvascular networks out of well‐segmented vessels. Our algorithm is based on 3D‐1D coupled models for simulating flow and oxygen transport as well as Murray's law to determine the radii of new vessels and bifurcations. The resulting networks are in good agreement with a segmented network related to morphological features.
Full text
Available for:
FZAB, GIS, IJS, IZUM, KILJ, NLZOH, NUK, OILJ, PILJ, PNG, SAZU, SBCE, SBMB, UL, UM, UPUK
We present novel H(div) and H1 liftings of given piecewise polynomials over a hierarchy of simplicial meshes, based on a global solve on the coarsest mesh and on local solves on patches of mesh ...elements around vertices on subsequent mesh levels. This in particular allows to lift a given algebraic residual. In connection with approaches lifting the total residual, we show how to obtain guaranteed, fully computable, and constant-free upper and lower a posteriori bounds on the algebraic, total, and discretization errors; here we consider the model Poisson equation discretized by the conforming finite element method of arbitrary order and including an arbitrary iterative solver. We next formulate safe stopping criteria ensuring that the algebraic error does not dominate the total error. We also prove efficiency, i.e., equivalence of our upper total and algebraic estimates with the total and algebraic errors, respectively, up to a generic constant; this constant is polynomial-degree-independent for the total error. Numerical experiments illustrate sharp control of all error components and accurate prediction of their spatial distribution in several test problems, including cases where some classical estimators fail. The H(div)-liftings at the same time allow to recover mass balance for any problem, any numerical discretization, and any situation such as inexact solution of (nonlinear) algebraic systems or algorithm failure, which we believe is of independent interest. We demonstrate this mass balance recovery in a simulation of immiscible incompressible two-phase flow in porous media.
•Novel H(div) and H1 liftings of given piecewise polynomials.•Designed over a hierarchy of general simplicial meshes.•Upper and lower a posteriori bounds on algebraic, total, and discretization errors.•Guaranteed, fully computable, constant-free, and efficient estimators.•Safe stopping criteria: the algebraic error does not dominate the total error.•Mass balance recovery for any problem, numerical discretization, and situation.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
This paper presents an algorithm for solving quasi-static, non-linear elasticity contact problems with friction in the context of rough surfaces. Here, we want to model the transition from sticking ...to slipping also called micro slip in a physically correct way in order to reproduce measured frictional damping. The popular dual Mortar method is used to enforce the contact constraints in a variationally consistent way without increasing the algebraic system size. The algorithm is deduced from a perturbed Lagrange formulation and combined with a serial–parallel Iwan model. This leads to a regularized saddle point problem, for which a non-linear complementary function and thus a semi-smooth Newton method can be derived. Numerical examples demonstrate the applicability to industrial problems and show good agreement to experimentally obtained results.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UL, UM, UPCLJ, UPUK
We consider a Cauchy problem for the inhomogeneous differential equation given in terms of an unbounded linear operator A and the Caputo fractional derivative of order α∈(0,2) in time. The previously ...known representation of the mild solution to such a problem does not have a conventional variation-of-constants like form, with the propagator derived from the associated homogeneous problem. Instead, it relies on the existence of two propagators with different analytical properties. This fact limits theoretical and especially numerical applicability of the existing solution representation. Here, we propose an alternative representation of the mild solution to the given problem that consolidates the solution formulas for sub-parabolic, parabolic and sub-hyperbolic equations with a positive sectorial operator A and non-zero initial data. The new representation is solely based on the propagator of the homogeneous problem and, therefore, can be considered as a more natural fractional extension of the solution to the classical parabolic Cauchy problem. By exploiting a trade-off between the regularity assumptions on the initial data in terms of the fractional powers of A and the regularity assumptions on the right-hand side in time, we show that the proposed solution formula is strongly convergent for t≥0 under considerably weaker assumptions compared to the standard results from the literature. Crucially, the achieved relaxation of space regularity assumptions ensures that the new solution representation is practically feasible for any α∈(0,2) and is amenable to the numerical evaluation using uniformly accurate quadrature-based algorithms.
We present an exponentially convergent numerical method to approximate the solution of the Cauchy problem for the inhomogeneous fractional differential equation with an unbounded operator coefficient ...and Caputo fractional derivative in time. The numerical method is based on the newly obtained solution formula that consolidates the mild solution representations of sub-parabolic, parabolic and sub-hyperbolic equations with sectorial operator coefficient A and non-zero initial data. The involved integral operators are approximated using the sinc-quadrature formulas that are tailored to the spectral parameters of A, fractional order α and the smoothness of the first initial condition, as well as to the properties of the equation’s right-hand side f(t). The resulting method possesses exponential convergence for positive sectorial A, any finite t, including t=0 and the whole range α∈(0,2). It is suitable for a practically important case, when no knowledge of f(t) is available outside the considered interval t∈0,T. The algorithm of the method is capable of multi-level parallelism. We provide numerical examples that confirm the theoretical error estimates.
Time-fractional partial differential equations are nonlocal-in-time and show an innate memory effect. Previously, examples like the time-fractional Cahn-Hilliard and Fokker-Planck equations have been ...studied. In this work, we propose a general framework of time-fractional gradient flows and we provide a rigorous analysis of well-posedness using the Faedo-Galerkin approach. Furthermore, we investigate the monotonicity of the energy functional of time-fractional gradient flows. Interestingly, it is still an open problem whether the energy is dissipating in time. This property is essential for integer-order gradient flows and many numerical schemes exploit this steepest descent characterization. We propose an augmented energy functional, which includes the history of the solution. Based on this new energy, we prove the equivalence of a time-fractional gradient flow to an integer-order one. This correlation guarantees the dissipating character of the augmented energy. The state function of the integer-order gradient flow acts on an extended domain similar to the Caffarelli-Silvestre extension for the fractional Laplacian. Additionally, we present a numerical scheme for solving time-fractional gradient flows, which is based on kernel compressing methods and reduces the problem to a system of ordinary differential equations. We illustrate the behavior of the original and augmented energy in the case of the Ginzburg-Landau energy.