COMPLAS 2021 is the 16th conference of the COMPLAS Series.
The COMPLAS conferences started in 1987 and since then have become established events in the field of computational plasticity and related topics. The first fifteen conferences in the COMPLAS series were all held in the city of Barcelona (Spain) and were very successful from the scientific, engineering and social points of view. We intend to make the 16th edition of the conferenceanother successful edition of the COMPLAS meetings.
The objectives of COMPLAS 2021 are to address both the theoretical bases for the solution of nonlinear solid mechanics problems, involving plasticity and other material nonlinearities, and the numerical algorithms necessary for efficient and robust computer implementation. COMPLAS 2021 aims to act as a forum for practitioners in the nonlinear structural mechanics field to discuss recent advances and identify future research directions.
Scope
COMPLAS 2021 is the 16th conference of the COMPLAS Series.
In this contribution, we shall be concerned with the question of wellposedness of thixoviscoplastic flow problems in context of FEM approximations.We restrict our analysis to a quasiNewtonian modeling approach with the aim to set foundations for an efficient monolithic Newtonmultigrid solver. We present the wellposedness of viscoplastic subproblems and structure subproblems in parallel/independent fashion showing the possibility for a combined treatment. Then, we use the fixed point theorem for the coupled problem. For the numerical solutions, we choose 4:1 contraction configuration and use monolithic Newton-multigrid solver. We analyse the effect of taking into consideration thixotropic phenomena in viscoplastic material and opening up for more different coupling by inclusions of shear thickening and shear thinning behaviors for plastic viscosity and/or elastic behavior below the critical yield stress limit in more a general thixotropic models.
Abstract In this contribution, we shall be concerned with the question of wellposedness of thixoviscoplastic flow problems in context of FEM approximations.We restrict our analysis [...]
We formulate a material model for micro-magneto-mechanics based on the generalized standard material approach. Our model includes exchange, elastic, anisotropy, demagnetizing and Zeeman energy. Furthermore we account for dissipative micro-magnetic behavior by means of a dissipation potential. For the constrained optimization problem w.r.t. magnetization we rely on the exponential map algorithm. We demonstrate our ideas with numerical examples. In particular we apply our model to a thin film composite. With this composite we represent the magneto-mechanical part of a magneto-electric composite sensor (resp. small sensor segment). Our numerical experiments focus on FeCoSiB as the magnetostrictive material. We discuss the coupling effects for the considered thin film composite in detail.
Abstract We formulate a material model for micro-magneto-mechanics based on the generalized standard material approach. Our model includes exchange, elastic, anisotropy, demagnetizing [...]
Sensitivity analysis is considered a fundamental tool in aerospace engineering, allowing to evaluate the impact of parameters variations, and optimize the aerodynamic and structural design. There has been much effort in the development of both theoretical and numerical frameworks for sensitivity calculations through adjoint solutions. Regarding the implementation of these analyses into agile and versatile numerical tools, the use of scalable and transferable libraries has become of particular importance. In this work, FEniCS is used to calculate the sensitivity of aerodynamic observables in different flow condition. A comparison with different theoretical and benchmark cases is used to validate the methodology before applying it to further and more complex configurations, expanding the scope of the analyses to unsteady flows.
Abstract Sensitivity analysis is considered a fundamental tool in aerospace engineering, allowing to evaluate the impact of parameters variations, and optimize the aerodynamic and [...]
We derive an efficient numerical scheme for transient wave propagation in configurations where a fluid domain and a solid domain are separated by a thin coating material. These type of configurations are numerically challenging for multiple factors: managing fluid solid coupling, enabling non-conform space discretizations, and rendering robust time discrete algorithm w.r.t the thin layer thickness. By combining the mortar element method with effective transmission conditions we are able to address these challenges. We illustrate our approach by proposing relevant 2D numerical illustrations inspired from simple ultrasonic testing experiments.
Abstract We derive an efficient numerical scheme for transient wave propagation in configurations where a fluid domain and a solid domain are separated by a thin coating material. [...]
Simulating multi-scale phenomena such as turbulent fluid flows is typically computationally very expensive. Filtering the smaller scales allows for using coarse discretizations, however, this requires closure models to account for the effects of the unresolved on the resolved scales. The common approach is to filter the continuous equations, but this gives rise to several commutator errors due to nonlinear terms, non-uniform filters, or boundary conditions. We propose a new approach to filtering, where the equations are discretized first and then filtered. For a non-uniform filter applied to the linear convection equation, we show that the discretely filtered convection operator can be inferred using three methods: intrusive (`explicit reconstruction') or non-intrusive operator inference, either via `derivative fitting' or `trajectory fitting' (embedded learning). We show that explicit reconstruction and derivative fitting identify a similar operator and produce small errors, but that trajectory fitting requires significant effort to train to achieve similar performance. However, the explicit reconstruction approach is more prone to instabilities.
Abstract Simulating multi-scale phenomena such as turbulent fluid flows is typically computationally very expensive. Filtering the smaller scales allows for using coarse discretizations, [...]
In this paper a robust topology optimization algorithm for linear elastic structures in unilateral contact is presented. The deformation of the linear elastic structure is constrained by support structures that are modeled with the help of Signorini's contact conditions. The contact conditions in turn are enforced with the augmented Lagrangian approach. Doing so, the robust optimization considers uncertainties at the support such as manufacturing tolerances and its local friction behavior. Due to high numerical costs in robust optimization, the firstorder second-moment approach is used to approximate the mean and variance of the objective. This approximation results in minimal additional costs to approximate the mean, the variance and their gradients. Consequentially, a gradient-based optimization algorithm can be used to minimize a weighted sum of both. The results show that the presented approach indeed improves the robustness with respect to uncertain contact conditions compared to a deterministic optimization.
Abstract In this paper a robust topology optimization algorithm for linear elastic structures in unilateral contact is presented. The deformation of the linear elastic structure is [...]
This paper reports on numerical experiments on arterial bypass-graft anastomoses. Bypass-grafts are oftentimes used in surgical procedures to divert blood around narrowed or occluded parts of an artery. The diverted blood flow is crucial to the success of the operation as it may lead to undesirable peculiarities that can result to a renewed occlusion in the distal connection of the graft. However, an a priori prediction of detrimental hemodynamic aspects due to undesirable flow properties is difficult to perform in vitro or in vivo conditions. To this end, this work targets to enhance our understanding of harming mechanisms through in silico experiments using computational fluid dynamics (CFD) and fluid-structure interaction (FSI) simulations. The latter are realized through a partitioned coupled approach which is verified for a 2D benchmark case against literature-reported results. Finally, we present numerical results on grafts with different cuff sizes. Wall shear stress (WSS), oscillatory shear index (OSI) and hemolysis are monitored and compared in the context of either rigid or elastic walls and cuff sizes. Special interest is given to the prediction of hemolysis induction which is often not considered in such studies. We show that wall elasticity is the key parameter in terms of WSS prediction while cuff size mainly affects the estimation of OSI.
Abstract This paper reports on numerical experiments on arterial bypass-graft anastomoses. Bypass-grafts are oftentimes used in surgical procedures to divert blood around narrowed [...]
In computational fluid dynamics (CFD), unsteady computations are cost-intensive. The Harmonic Balance (HB) method [7] represents a cost-efficient alternative. Here, for timeperiodic flows, the governing equations are recasted in the Fourier domain. In the low Mach regime, the compressible governing equations are stiff. Therefore, densitybased solvers converge slowly. Low Mach preconditioning equalizes the eigenvalues of the system of equations, to improve the condition number and remove the stiffness of the system [14]. In this paper, low Mach preconditioning is applied to the HB method, with emphasis on the non-reflecting boundary conditions (NRBCs). These boundary conditions have a crucial impact on the flow inside the truncated computational domains used in CFD. Improper boundary conditions reflect waves exiting the computational domain and deteriorate the quality of the solution. However, NRBCs [6] avoid spurious reflections. We explain that to precondition the NRBC its formulation in terms of characteristics has to be adapted. An academic wave propagation test case is computed for different wave configurations to validate the preconditioned boundary conditions. The use of non-preconditioned NRBCs in a preconditioned computation leads to instabilities and reflections at the boundaries of the domain. A consistent setup with preconditioned NRBCs improves the stability and leads to good non-reflecting properties for all presented wave configurations.
Abstract In computational fluid dynamics (CFD), unsteady computations are cost-intensive. The Harmonic Balance (HB) method [7] represents a cost-efficient alternative. Here, for timeperiodic [...]
N. Pynaert, J. Wauters, G. Crevecoeur, J. Degroote
eccomas2022.
Abstract
Airborne wind energy (AWE) is an emerging technology for the conversion of wind energy into electricity by flying crosswind patterns with a tethered aircraft connected to a generator either on board or on the ground. Having a proper understanding of the unsteady interaction of the air with the flexible and dynamic system during operation is key to developing viable AWE systems. The research goal is to simulate the time-varying fluid-structure interaction (FSI) of an AWE system in a crosswind flight maneuver using high fidelity simulation tools. In this work a framework is presented that serves as a proof of concept to perform high fidelity simulations of airborne wind energy systems. This is done using a partitioned and explicit approach in the open-source coupling tool CoCoNuT. An existing finite element method (FEM) model of the wing structure is coupled with a newly developed computational fluid dynamics (CFD) model of the wing aerodynamics including rigid body motion. It has been found that the mesh deformation is quite sensitive to dynamic mesh parameters. On the other hand, the overset/Chimera technique has been proven to be a robust approach to simulate the motion of an AWE system in CFD simulations.
Abstract Airborne wind energy (AWE) is an emerging technology for the conversion of wind energy into electricity by flying crosswind patterns with a tethered aircraft connected to [...]
F. Kristoffersen, M. Larsson, S. Johnsen, W. Schröder, B. Müller
eccomas2022.
Abstract
A 3D fluid-structure interaction (FSI) code is under development. The fluid domain (Navier-Stokes) solver will employ a sharp interface ghost node immersed boundary method (IBM) to apply the boundary conditions at fluid-solid interfaces. The Navier-Stokes (N-S) solver has been verified using a classic Poiseuille channel flow. The current version of the immersed boundary method is being tested by solving a heat conduction problem. The order of accuracy of the IBM was shown to be just above second order.
Abstract A 3D fluid-structure interaction (FSI) code is under development. The fluid domain (Navier-Stokes) solver will employ a sharp interface ghost node immersed boundary method [...]