297 research outputs found
Error-estimate-based adaptive integration for immersed isogeometric analysis
The Finite Cell Method (FCM) together with Isogeometric analysis (IGA) has been applied successfully in various problems in solid mechanics, in image-based analysis, fluid–structure interaction and in many other applications. A challenging aspect of the isogeometric finite cell method is the integration of cut cells. In particular in three-dimensional simulations the computational effort associated with integration can be the critical component of a simulation. A myriad of integration strategies has been proposed over the past years to ameliorate the difficulties associated with integration, but a general optimal integration framework that suits a broad class of engineering problems is not yet available. In this contribution we provide a thorough investigation of the accuracy and computational effort of the octree integration scheme. We quantify the contribution of the integration error using the theoretical basis provided by Strang's first lemma. Based on this study we propose an error-estimate-based adaptive integration procedure for immersed isogeometric analysis. Additionally, we present a detailed numerical investigation of the proposed optimal integration algorithm and its application to immersed isogeometric analysis using two- and three-dimensional linear elasticity problems
Error-estimate-based Adaptive Integration For Immersed Isogeometric Analysis
The Finite Cell Method (FCM) together with Isogeometric analysis (IGA) has
been applied successfully in various problems in solid mechanics, in
image-based analysis, fluid-structure interaction and in many other
applications. A challenging aspect of the isogeometric finite cell method is
the integration of cut cells. In particular in three-dimensional simulations
the computational effort associated with integration can be the critical
component of a simulation. A myriad of integration strategies has been proposed
over the past years to ameliorate the difficulties associated with integration,
but a general optimal integration framework that suits a broad class of
engineering problems is not yet available. In this contribution we provide a
thorough investigation of the accuracy and computational effort of the octree
integration scheme. We quantify the contribution of the integration error using
the theoretical basis provided by Strang's first lemma. Based on this study we
propose an error-estimate-based adaptive integration procedure for immersed
isogeometric analysis. Additionally, we present a detailed numerical
investigation of the proposed optimal integration algorithm and its application
to immersed isogeometric analysis using two- and three-dimensional linear
elasticity problems.Comment: To CAMW
Isogeometric Analysis on V-reps: first results
Inspired by the introduction of Volumetric Modeling via volumetric
representations (V-reps) by Massarwi and Elber in 2016, in this paper we
present a novel approach for the construction of isogeometric numerical methods
for elliptic PDEs on trimmed geometries, seen as a special class of more
general V-reps. We develop tools for approximation and local re-parametrization
of trimmed elements for three dimensional problems, and we provide a
theoretical framework that fully justify our algorithmic choices. We validate
our approach both on two and three dimensional problems, for diffusion and
linear elasticity.Comment: 36 pages, 44 figures. Reviewed versio
Scan-based immersed isogeometric analysis
Scan-based simulations contain innate topologically complex three-dimensional geometries, represented by large data sets in formats which are not directly suitable for analysis. Consequently, performing high-fidelity scan-based simulations at practical computational costs is still very challenging. The main objective of this dissertation is to develop an efficient and robust scan-based simulation strategy by acquiring a profound understanding of three prominent challenges in scan-based IGA, viz.: i) balancing the accuracy and computational effort associated with numerical integration; ii) the preservation of topology in the spline-based segmentation procedure; and iii) the control of accuracy using error estimation and adaptivity techniques.
In three-dimensional immersed isogeometric simulations, the computational effort associated with integration can be the critical component. A myriad of integration strategies has been proposed over the past years to ameliorate the difficulties associated with integration, but a general optimal integration framework that suits a broad class of engineering problems is not yet available. In this dissertation we provide a thorough investigation of the accuracy and computational effort of the octree integration technique. We quantify the contribution of the integration error using the theoretical basis provided by Strang’s first lemma. Based on this study we propose an error-estimate-based adaptive integration procedure for immersed IGA.
To exploit the advantageous properties of IGA in a scan-based setting, it is important to extract a smooth geometry. This can be established by convoluting the voxel data using B-splines, but this can induce problematic topological changes when features with a size similar to that of the voxels are encountered. This dissertation presents a topology-preserving segmentation procedure using truncated hierarchical (TH)B-splines. A moving-window-based topological anomaly detection algorithm is proposed to identify regions in which (TH)B-spline refinements must be performed. The criterion to identify topological anomalies is based on the Euler characteristic, giving it the capability to distinguish between topological and shape changes. A Fourier analysis is presented to explain the effectiveness of the developed procedure.
An additional computational challenge in the context of immersed IGA is the construction of optimal approximations using locally refined splines. For scan-based volumetric domains, hierarchical splines are particularly suitable, as they optimally leverage the advantages offered by the availability of a geometrically simple background mesh. Although truncated hierarchical B-splines have been successfully applied in the context of IGA, their application in the immersed setting is largely unexplored. In this dissertation we propose a computational strategy for the application of error estimation-based mesh adaptivity for stabilized immersed IGA.
The conducted analyses and developed computational techniques for scan-based immersed IGA are interrelated, and together constitute a significant improvement in the efficiency and robustness of the analysis paradigm. In combination with other state-of-the-art developments regarding immersed FEM/IGA (\emph{e.g.}, iterative solution techniques, parallel computing), the research in this thesis opens the doors to scan-based simulations with more sophisticated physical behavior, geometries of increased complexity, and larger scan-data sizes.Scan-based simulations contain innate topologically complex three-dimensional geometries, represented by large data sets in formats which are not directly suitable for analysis. Consequently, performing high-fidelity scan-based simulations at practical computational costs is still very challenging. The main objective of this dissertation is to develop an efficient and robust scan-based simulation strategy by acquiring a profound understanding of three prominent challenges in scan-based IGA, viz.: i) balancing the accuracy and computational effort associated with numerical integration; ii) the preservation of topology in the spline-based segmentation procedure; and iii) the control of accuracy using error estimation and adaptivity techniques.
In three-dimensional immersed isogeometric simulations, the computational effort associated with integration can be the critical component. A myriad of integration strategies has been proposed over the past years to ameliorate the difficulties associated with integration, but a general optimal integration framework that suits a broad class of engineering problems is not yet available. In this dissertation we provide a thorough investigation of the accuracy and computational effort of the octree integration technique. We quantify the contribution of the integration error using the theoretical basis provided by Strang’s first lemma. Based on this study we propose an error-estimate-based adaptive integration procedure for immersed IGA.
To exploit the advantageous properties of IGA in a scan-based setting, it is important to extract a smooth geometry. This can be established by convoluting the voxel data using B-splines, but this can induce problematic topological changes when features with a size similar to that of the voxels are encountered. This dissertation presents a topology-preserving segmentation procedure using truncated hierarchical (TH)B-splines. A moving-window-based topological anomaly detection algorithm is proposed to identify regions in which (TH)B-spline refinements must be performed. The criterion to identify topological anomalies is based on the Euler characteristic, giving it the capability to distinguish between topological and shape changes. A Fourier analysis is presented to explain the effectiveness of the developed procedure.
An additional computational challenge in the context of immersed IGA is the construction of optimal approximations using locally refined splines. For scan-based volumetric domains, hierarchical splines are particularly suitable, as they optimally leverage the advantages offered by the availability of a geometrically simple background mesh. Although truncated hierarchical B-splines have been successfully applied in the context of IGA, their application in the immersed setting is largely unexplored. In this dissertation we propose a computational strategy for the application of error estimation-based mesh adaptivity for stabilized immersed IGA.
The conducted analyses and developed computational techniques for scan-based immersed IGA are interrelated, and together constitute a significant improvement in the efficiency and robustness of the analysis paradigm. In combination with other state-of-the-art developments regarding immersed FEM/IGA (\emph{e.g.}, iterative solution techniques, parallel computing), the research in this thesis opens the doors to scan-based simulations with more sophisticated physical behavior, geometries of increased complexity, and larger scan-data sizes
Isogeometric FEM-BEM coupled structural-acoustic analysis of shells using subdivision surfaces
We introduce a coupled finite and boundary element formulation for acoustic
scattering analysis over thin shell structures. A triangular Loop subdivision
surface discretisation is used for both geometry and analysis fields. The
Kirchhoff-Love shell equation is discretised with the finite element method and
the Helmholtz equation for the acoustic field with the boundary element method.
The use of the boundary element formulation allows the elegant handling of
infinite domains and precludes the need for volumetric meshing. In the present
work the subdivision control meshes for the shell displacements and the
acoustic pressures have the same resolution. The corresponding smooth
subdivision basis functions have the continuity property required for the
Kirchhoff-Love formulation and are highly efficient for the acoustic field
computations. We validate the proposed isogeometric formulation through a
closed-form solution of acoustic scattering over a thin shell sphere.
Furthermore, we demonstrate the ability of the proposed approach to handle
complex geometries with arbitrary topology that provides an integrated
isogeometric design and analysis workflow for coupled structural-acoustic
analysis of shells
Suitably graded THB-spline refinement and coarsening: Towards an adaptive isogeometric analysis of additive manufacturing processes
In the present work we introduce a complete set of algorithms to efficiently
perform adaptive refinement and coarsening by exploiting truncated hierarchical
B-splines (THB-splines) defined on suitably graded isogeometric meshes, that
are called admissible mesh configurations. We apply the proposed algorithms to
two-dimensional linear heat transfer problems with localized moving heat
source, as simplified models for additive manufacturing applications. We first
verify the accuracy of the admissible adaptive scheme with respect to an
overkilled solution, for then comparing our results with similar schemes which
consider different refinement and coarsening algorithms, with or without taking
into account grading parameters. This study shows that the THB-spline
admissible solution delivers an optimal discretization for what concerns not
only the accuracy of the approximation, but also the (reduced) number of
degrees of freedom per time step. In the last example we investigate the
capability of the algorithms to approximate the thermal history of the problem
for a more complicated source path. The comparison with uniform and
non-admissible hierarchical meshes demonstrates that also in this case our
adaptive scheme returns the desired accuracy while strongly improving the
computational efficiency.Comment: 20 pages, 12 figure
A note on the penalty parameter in Nitsche's method for unfitted boundary value problems
Nitsche's method is a popular approach to implement Dirichlet-type boundary
conditions in situations where a strong imposition is either inconvenient or
simply not feasible. The method is widely applied in the context of unfitted
finite element methods. From the classical (symmetric) Nitsche's method it is
well-known that the stabilization parameter in the method has to be chosen
sufficiently large to obtain unique solvability of discrete systems. In this
short note we discuss an often used strategy to set the stabilization parameter
and describe a possible problem that can arise from this. We show that in
specific situations error bounds can deteriorate and give examples of
computations where Nitsche's method yields large and even diverging
discretization errors
Residual-based error estimation and adaptivity for stabilized immersed isogeometric analysis using truncated hierarchical B-splines
We propose an adaptive mesh refinement strategy for immersed isogeometric
analysis, with application to steady heat conduction and viscous flow problems.
The proposed strategy is based on residual-based error estimation, which has
been tailored to the immersed setting by the incorporation of appropriately
scaled stabilization and boundary terms. Element-wise error indicators are
elaborated for the Laplace and Stokes problems, and a THB-spline-based local
mesh refinement strategy is proposed. The error estimation .and adaptivity
procedure is applied to a series of benchmark problems, demonstrating the
suitability of the technique for a range of smooth and non-smooth problems. The
adaptivity strategy is also integrated in a scan-based analysis workflow,
capable of generating reliable, error-controlled, results from scan data,
without the need for extensive user interactions or interventions.Comment: Submitted to Journal of Mechanic
Condition number analysis and preconditioning of the finite cell method
The (Isogeometric) Finite Cell Method - in which a domain is immersed in a
structured background mesh - suffers from conditioning problems when cells with
small volume fractions occur. In this contribution, we establish a rigorous
scaling relation between the condition number of (I)FCM system matrices and the
smallest cell volume fraction. Ill-conditioning stems either from basis
functions being small on cells with small volume fractions, or from basis
functions being nearly linearly dependent on such cells. Based on these two
sources of ill-conditioning, an algebraic preconditioning technique is
developed, which is referred to as Symmetric Incomplete Permuted Inverse
Cholesky (SIPIC). A detailed numerical investigation of the effectivity of the
SIPIC preconditioner in improving (I)FCM condition numbers and in improving the
convergence speed and accuracy of iterative solvers is presented for the
Poisson problem and for two- and three-dimensional problems in linear
elasticity, in which Nitche's method is applied in either the normal or
tangential direction. The accuracy of the preconditioned iterative solver
enables mesh convergence studies of the finite cell method
- …