| 1. |
- Siemers, Alexander, 1970-
(författare)
-
Contributions to Modelling and Visualisation of Multibody Systems Simulations with Detailed Contact Analysis
- 2010
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- The steadily increasing performance of modern computer systems is having a large influence on simulation technologies. It enables increasingly detailed simulations of larger and more comprehensive simulation models. Increasingly large amounts of numerical data are produced by these simulations.This thesis presents several contributions in the field of mechanical system simulation and visualisation. The work described in the thesis is of practical relevance and results have been tested and implemented in tools that are used daily in the industry i.e., the BEAST (BEAring Simulation Tool) tool box. BEAST is a multibody system (MBS) simulation software with special focus on detailed contact calculations. Our work is primarily focusing on these types of systems.focusing on these types of systems. Research in the field of simulation modelling typically focuses on one or several specific topics around the modelling and simulation work process. The work presented here is novel in the sense that it provides a complete analysis and tool chain for the whole work process for simulation modelling and analysis of multibody systems with detailed contact models. The focus is on detecting and dealing with possible problems and bottlenecks in the work process, with respect to multibody systems with detailed contact models.The following primary research questions have been formulated:How to utilise object-oriented techniques for modelling of multibody systems with special reference tocontact modelling?How to integrate visualisation with the modelling and simulation process of multibody systems withdetailed contacts.How to reuse and combine existing simulation models to simulate large mechanical systems consistingof several sub-systems by means of co-simulation modelling?Unique in this work is the focus on detailed contact models. Most modelling approaches for multibody systems focus on modelling of bodies and boundary conditions of such bodies, e.g., springs, dampers, and possibly simple contacts. Here an object oriented modelling approach for multibody simulation and modelling is presented that, in comparison to common approaches, puts emphasis on integrated contact modelling and visualisation. The visualisation techniques are commonly used to verify the system model visually and to analyse simulation results. Data visualisation covers a broad spectrum within research and development. The focus is often on detailed solutions covering a fraction of the whole visualisation process. The novel visualisation aspect of the work presented here is that it presents techniques covering the entire visualisation process integrated with modeling and simulation. This includes a novel data structure for efficient storage and visualisation of multidimensional transient surface related data from detailed contact calculations.Different mechanical system simulation models typically focus on different parts (sub-systems) of a system. To fully understand a complete mechanical system it is often necessary to investigate several or all parts simultaneously. One solution for a more complete system analysis is to couple different simulation models into one coherent simulation. Part of this work is concerned with such co-simulation modelling. Co-simulation modelling typically focuses on data handling, connection modelling, and numerical stability. This work puts all emphasis on ease of use, i.e., making mechanical system co-simulation modelling applicable for a larger group of people. A novel meta-model based approach for mechanical system co-simulation modelling is presented. The meta-modelling process has been defined and tools and techniques been created to fully support the complete process. A component integrator and modelling environment are presented that support automated interface detection, interface alignment with automated three-dimensional coordinate translations, and three dimensional visual co-simulation modelling. The integrated simulator is based on a general framework for mechanical system co-simulations that guarantees numerical stability.
|
|
| 2. |
- Araujo-Cabarcas, Juan Carlos, 1981-
(författare)
-
Reliable hp finite element computations of scattering resonances in nano optics
- 2019
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Eigenfrequencies are commonly studied in wave propagation problems, as they are important in the analysis of closed cavities such as a microwave oven. For open systems, energy leaks into infinity and therefore scattering resonances are used instead of eigenfrequencies. An interesting application where resonances take an important place is in whispering gallery mode resonators.The objective of the thesis is the reliable and accurate approximation of scattering resonances using high order finite element methods. The discussion focuses on the electromagnetic scattering resonances in metal-dielectric nano-structures using a Drude-Lorentz model for the description of the material properties. A scattering resonance pair satisfies a reduced wave equationand an outgoing wave condition. In this thesis, the outgoing wave condition is replaced by a Dirichlet-to-Neumann map, or a Perfectly Matched Layer. For electromagnetic waves and for acoustic waves, the reduced wave equation is discretized with finite elements. As a result, the scattering resonance problem is transformed into a nonlinear eigenvalue problem.In addition to the correct approximation of the true resonances, a large number of numerical solutions that are unrelated to the physical problem are also computed in the solution process. A new method based on a volume integral equation is developed to remove these false solutions.The main results of the thesis are a novel method for removing false solutions of the physical problem, efficient solutions of non-linear eigenvalue problems, and a new a-priori based refinement strategy for high order finite element methods. The overall material in the thesis translates into a reliable and accurate method to compute scattering resonances in physics and engineering.
|
|
| 3. |
- Johansson, August, 1978-
(författare)
-
Duality-based adaptive finite element methods with application to time-dependent problems
- 2010
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- To simulate real world problems modeled by differential equations, it is often not sufficient to consider and tackle a single equation. Rather, complex phenomena are modeled by several partial dierential equations that are coupled to each other. For example, a heart beat involve electric activity, mechanics of the movement of the walls and valves, as well as blood fow - a true multiphysics problem. There may also be ordinary differential equations modeling the reactions on a cellular level, and these may act on a much finer scale in both space and time. Determining efficient and accurate simulation tools for such multiscalar multiphysics problems is a challenge. The five scientific papers constituting this thesis investigate and present solutions to issues regarding accurate and efficient simulation using adaptive finite element methods. These include handling local accuracy through submodeling, analyzing error propagation in time-dependent multiphysics problems, developing efficient algorithms for adaptivity in time and space, and deriving error analysis for coupled PDE-ODE systems. In all these examples, the error is analyzed and controlled using the framework of dual-weighted residuals, and the spatial meshes are handled using octree based data structures. However, few realistic geometries fit such grid and to address this issue a discontinuous Galerkin Nitsche method is presented and analyzed.
|
|
| 4. |
|
|
| 5. |
- Jonsson, Tobias, 1991-
(författare)
-
Cut isogeometric methods on trimmed multipatch surfaces
- 2022
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Partial differential equations (PDE) on surfaces appear in a variety of applications, such as image processing, modeling of lubrication, fluid flows, diffusion, and transport of surfactants. In some applications, surfaces are drawn and modeled by using CAD software, giving a very precise patchwise parametric description of the surface. This thesis deals with the development of methods for finding numerical solutions to PDE posed on such parametrically described multipatch surfaces. The thesis consists of an introduction and five papers.In the first paper, we develop a general framework for the Laplace-Beltrami operator on a patchwise parametric surface. Each patch map induces a Riemannian metric, which we utilize to compute quantities in the simpler reference domain. We use the cut finite element method together with Nitsche’s method to enforce continuity over the interfaces between patches.In the second paper, we extend the framework to be able to handle geometries that consist of an arrangement of surfaces, i.e., more than two per interface. By using a Kirchhoff's condition this method avoids defining any co-normal to each surface and can deal with sharp edges. This approach is shown to be equivalent to standard Nitsche interface method for flat geometries.In the third paper, we developed a cut finite element method for elliptic problems with corner singularities. The main idea is to use an appropriate radial map that grades the finite element mesh towards the corner that counter-acts the solution's singularity.In the fourth paper, we present a new robust isogeometric method for surfaces described by CAD patches with gaps or overlaps. The main approach here is to cover all interfaces with a three-dimensional mesh and then use a hybrid variable in a Nitsche-type formulation to transfer data over the gaps. Using this hybridized approach leads to a convenient and easy to implement method with no restriction on the number of coupled patches per interface.In the fifth paper, we present a routine to the multipatch isogeometric framework for dealing with singular maps. To exemplify this, we consider a specific type of singular parametrization which essentially maps a square onto a triangle. One part of the boundary of the square will be transformed into a single point and the metric tensor becomes singular as we approach this boundary. In this work we propose a regularization procedure which is based on eigenvalue decomposition of the metric tensor.
|
|
| 6. |
- Larsson, Karl, 1981-
(författare)
-
Finite element methods for threads and plates with real-time applications
- 2010
-
Licentiatavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Thin and slender structures are widely occurring both in nature and in human creations. Clever geometries of thin structures can produce strong constructions while using a minimal amount of material. Computer modeling and analysis of thin and slender structures has its own set of problems stemming from assumptions made when deriving the equations modeling their behavior from the theory of continuum mechanics. In this thesis we consider two kinds of thin elastic structures; threads and plates. Real-time simulation of threads are of interest in various types of virtual simulations such as surgery simulation for instance. In the first paper of this thesis we develop a thread model for use in interactive applications. By viewing the thread as a continuum rather than a truly one dimensional object existing in three dimensional space we derive a thread model that naturally handles both bending, torsion and inertial effects. We apply a corotational framework to simulate large deformation in real-time. On the fly adaptive resolution is used to minimize corotational artifacts. Plates are flat elastic structures only allowing deflection in the normal direction. In the second paper in this thesis we propose a family of finite elements for approximating solutions to the Kirchhoff-Love plate equation using a continuous piecewise linear deflection field. We reconstruct a discontinuous piecewise quadratic deflection field which is applied in a discontinuous Galerkin method. Given a criterion on the reconstruction operator we prove a priori estimates in energy and L2 norms. Numerical results for the method using three possible reconstructions are presented.
|
|
| 7. |
- Lindh, Magnus, 1973-
(författare)
-
Evolution of Plants : a mathematical perspective
- 2016
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- The Earth harbors around 300 000 plant species. The rich and complex environment provided by plants is considered a key factor for the extraordinary diversity of the terrestrial fauna by, for example, providing food and shelter. This thesis contributes to the understanding of these questions by investigating how the interplay of physiology, demography, and evolution gives rise to variation and diversity in fundamental plant traits. This will help us answer questions such as: How has this amazing diversity of plant species emerged? Which mechanisms maintain diversity? How are plant strategies and plant diversity influenced by variations in the environment?A plant faces multiple problems to survive and reproduce successfully. These problems can be modeled by considering traits, trade-offs and a fitness measure. For example: How to maximize growth rate, while maximizing structural stability? I will investigate four plant models in order to understand the function of plants, and mechanisms promoting diversity. Paper I: We study how annual plants with and without growth constraints should optimize their flowering time when productivity or season length changes. With a dynamic ontogenetic growth model and optimal control theory we prove that a bang-bang reproductive control is optimal under constrained growth and constant mortality rate. We find that growth constraints can flip the direction of optimal phenological response for increasing productivity. The reason is that the growth rate of vegetative mass saturates at high productivity and therefore it is better to flower earlier and take advantage of a longer reproductive period. If season length extends equally both in the beginning and the end of the season, growth constraints control the direction of the optimal response as well. Our theory can help explaining phenological patterns along productivity gradients, and can be linked to empirical observations made on a calendar scale.Paper II: We introduce a new measure of tree crown-rise efficiency based on the loss of biomass of the tree during growth. The more mass the tree looses during growth, the less crown-rise efficient it is. Top-heavy shapes loose more biomass than bottom-heavy shapes. Light-use efficiency is defined as the mean light assimilation of the leaves in the crown times the ratio of leaf mass and total mass. We then study the trade-off between light-use efficiency to crown-rise efficiency for tree crown shapes. We assume that the total tree mass is constant, and a constant vertical light gradient represent the shading from a surrounding forest. We find large differences in crown shapes at intermediate vertical light gradient, when both self-shading and mean-field shading are important, suggesting light-use vs crown-rise efficiency as a new trade-off that can explain tree diversity. Our crown-rise efficiency measure could easily be integrated into existing forest models.Paper III: We extend an evolutionary tree crown model, where trees with different heights compete for light, with drought-induced mortality rates depending on ground-water availability and the depth of an optional taproot. The model does not include competition for ground water. Our model explains how ground-water availability can shape plant communities, when taproot and non-taproot strategies can coexist, and when only one of these strategies can persist. We investigate how emerging plant diversity varies with water table depth, soil water gradient and drought-induced mortality rate. The taproot enables plants to reach deep water, thus reducing mortality, but also carries a construction cost, thus inducing a trade-off. We find that taproots maintain plant diversity under increasing drought mortality, and that taproots evolve when groundwater is accessible at low depths. There are no viable strategies at high drought mortality and deep water table. Red Queen evolutionary dynamics appear at intermediate drought mortality in mixed communities with and without taproots, as the community never reaches a final evolutionarily stable composition.Paper IV: We extend a size-structured plant model, with self-shading and two evolving traits, crown top-heaviness and crown width-to-height ratio. The model allows us to identify salient trade-offs for the crown shape. The most important trade-off for top-heaviness is light-use efficiency vs crownrise efficiency, and the most important trade-off for width-to-height ratio is self-shading vs branch costs. We find that when the two traits coevolve; the outcome is a single common evolutionarily stable strategy (ESS), far away from the highest net primary production (NPP). When only sun angle is decreasing with increasing latitude both the crown width-to-height ratio and crown top-heaviness decrease. However, when light response in addition to the sun angle decreases with increasing latitude, the crown width-to-height ratio is nearly invariant of latitude except at low site productivity when the ratio decreases with latitude. Top-heaviness is always decreasing with increasing latitude. Finally, we find that crown top-heaviness increases with the NPP or leaf-area index (LAI) at ESS, but crown width-to-height ratio is maximal at an intermediate NPP or LAI.
|
|
| 8. |
- Sticko, Simon, 1988-
(författare)
-
High Order Cut Finite Element Methods for Wave Equations
- 2018
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- This thesis considers wave propagation problems solved using finite element methods where a boundary or interface of the domain is not aligned with the computational mesh. Such methods are usually referred to as cut or immersed methods. The motivation for using immersed methods for wave propagation comes largely from scattering problems when the geometry of the domain is not known a priori. For wave propagation problems, the amount of computational work per dispersion error is generally lower when using a high order method. For this reason, this thesis aims at studying high order immersed methods.Nitsche's method is a common way to assign boundary or interface conditions in immersed finite element methods. Here, penalty terms that are consistent with the boundary/interface conditions are added to the weak form. This requires that special quadrature rules are constructed on the intersected elements, which take the location of the immersed boundary/interface into account. A common problem for all immersed methods is small cuts occurring between the elements in the mesh and the computational domain. A suggested way to remedy this is to add terms penalizing jumps in normal derivatives over the faces of the intersected elements.Paper I and Paper II consider the acoustic wave equation, using first order elements in Paper I, and using higher order elements in Paper II. High order elements are then used for the elastic wave equation in Paper III. Papers I to III all use continuous Galerkin, Nitsche's method, and jump-stabilization. Paper IV compares the errors of this type of cut finite element method with two other numerical methods. One result from Paper II is that the added jump-stabilization results in a mass matrix with a high condition number. This motivates the investigation of alternatives. Paper V considers a hybridizable discontinuous Galerkin method. This paper investigates to what extent local time stepping in combination with cell-merging can be used to overcome the problem of small cuts.
|
|
| 9. |
- Larsson, Karl, 1981-
(författare)
-
Finite Element Methods for Thin Structures with Applications in Solid Mechanics
- 2013
-
Doktorsavhandling (övrigt vetenskapligt/konstnärligt)abstract
- Thin and slender structures are widely occurring both in nature and in human creations. Clever geometries of thin structures can produce strong constructions while requiring a minimal amount of material. Computer modeling and analysis of thin and slender structures have their own set of problems, stemming from assumptions made when deriving the governing equations. This thesis deals with the derivation of numerical methods suitable for approximating solutions to problems on thin geometries. It consists of an introduction and four papers.In the first paper we introduce a thread model for use in interactive simulation. Based on a three-dimensional beam model, a corotational approach is used for interactive simulation speeds in combination with adaptive mesh resolution to maintain accuracy.In the second paper we present a family of continuous piecewise linear finite elements for thin plate problems. Patchwise reconstruction of a discontinuous piecewise quadratic deflection field allows us touse a discontinuous Galerkin method for the plate problem. Assuming a criterion on the reconstructions is fulfilled we prove a priori error estimates in energy norm and L2-norm and provide numerical results to support our findings.The third paper deals with the biharmonic equation on a surface embedded in R3. We extend theory and formalism, developed for the approximation of solutions to the Laplace-Beltrami problem on an implicitly defined surface, to also cover the biharmonic problem. A priori error estimates for a continuous/discontinuous Galerkin method is proven in energy norm and L2-norm, and we support the theoretical results by numerical convergence studies for problems on a sphere and on a torus.In the fourth paper we consider finite element modeling of curved beams in R3. We let the geometry of the beam be implicitly defined by a vector distance function. Starting from the three-dimensional equations of linear elasticity, we derive a weak formulation for a linear curved beam expressed in global coordinates. Numerical results from a finite element implementation based on these equations are compared with classical results.
|
|
