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.
Machine learning entails powerful information processing algorithms that are relevant for modelling, optimization, and control of fluids. Currently, machine-learning capabilities are advancing at an incredible rate, and fluid mechanics is beginning to tap into the full potential of these powerful methods. Many tasks in fluid mechanics, such as reduced-order modelling, shape optimization and uncertainty quantification, may be posed as optimization and regression tasks. Machine learning can dramatically improve optimization performance and reduce convergence time. In this paper, the potential of tree-based machine learning techniques for the aerodynamic prediction of pressure coefficients of an AIRBUS XRF1 aircraft wing-body configuration has been assessed. For this purpose, a dataset including computational fluid dynamics (CFD) simulations has been employed to train the different models, with and without the use of proper orthogonal decomposition (POD) and having their hyperparameters values optimized to obtain the optimal subspace. A deep comparison of decision tree regressors and random forest algorithms has been performed, showing that the random forest regressor model performs better on all configurations.
Abstract Machine learning entails powerful information processing algorithms that are relevant for modelling, optimization, and control of fluids. Currently, machine-learning capabilities [...]
Powder-based additive manufacturing technologies, specifically selective laser melting, are challenging to model due to the complex, interrelated physical phenomena that are present on multiple spatial scales, during the process. A key element of such models will be the detailed simulation of flow and heat transfer throughout the melt pool that is formed when the powder particles melt. Due to the high-temperature gradients that are generated inside the melt pool, the Marangoni force plays a key role in governing the flows inside the melt pool and deciding its shape and dimensions. On the other hand, the mass and heat transfer between the melt and the powder also has a significant role in shaping the melt pool at the edges. In this study, we modified an OpenFOAM solver (icoReactingMultiphaseInterFoam) coupled with an in-house developed DEM code known as eXtended Discrete Element Method or XDEM which models the dynamics and thermodynamics of the particles. By adding the Marangoni force to the momentum equation and also defining a laser model as a boundary condition for liquid-gas interface, the solver is capable of modeling the selective laser melting process from the moment of particle melting to the completion of the solidified track. The coupled solver was validated with an ice packed bed melting case and was used to simulate a multi-track selective laser melting process.
Abstract Powder-based additive manufacturing technologies, specifically selective laser melting, are challenging to model due to the complex, interrelated physical phenomena that are [...]
A fluid-structure interaction model is employed to numerically investigate the interaction between the pressurized thin lubricant film and the radial, plastically deformed steel wire in a dry wire drawing process. A transient simulation is presented, with the implementation of a sliding fluid-structure interaction interface. Moreover, the fluid film has been calculated by the Navier-Stokes equations and the coupling with the wire model is performed by the IQN-ILS technique. This results on the one hand in the monitoring of the stresses and displacements of the structure and on the other hand in an observation of the hydrodynamic pressure build-up and wall shear stresses in the lubricant. Additionally, the evolution of the thickness of the fluid film is presented.
Abstract A fluid-structure interaction model is employed to numerically investigate the interaction between the pressurized thin lubricant film and the radial, plastically deformed [...]
Aeroacoustics is the field of studying flow-induced sound, which results from the interaction of unsteady flow with solid structures, such as aircraft and automobiles. Different methods are available to achieve this, including theoretical, experimental, and computational methods. Due to the high costs of experiments, the concentration on computational methods has increased. Computational aeroacoustics (CAA), based on computational fluid dynamics (CFD), has received special attention from researchers because of its outstanding capability to get acceptable results with reasonable computational costs. The partially-Averaged Navier Stokes (PANS) method is a hybrid LES/RANS method based on dynamic resolution parameters. The SSVPANS method is a k---f based PANS method with an additional modeled equation for the resolved kinetic energy. This method has been implemented in FASTEST, an in-house finitevolume solver to compute the flows in complex applications. This study aims to investigate the aeroacoustic performance of the SSV-PANS method compared to a reference Large-eddy Simulation [1] regarding the computational accuracy and costs. To do this, hybrid method based on decomposing the fluid variables into incompressible hydrodynamics and compressible perturbation equations is used to study the aerodynamic noise. The aeroacoustic sources are computed from the incompressible flow field using the SSV-PANS method. In addition, the Kirchhoff wave extrapolation method is used to have an efficient evaluation of the far-field noise.
Abstract Aeroacoustics is the field of studying flow-induced sound, which results from the interaction of unsteady flow with solid structures, such as aircraft and automobiles. Different [...]
High-accuracy simulations of internal combustion engines (ICE) allow deep insight into the physical processes of the different phases of the engine cycle: gas exchange, mixture formation, compression, combustion and emission formation. The commercial solvers for ICE simulations provide a full package which covers these areas. However, the user of such software is unable to look into the source code, making it impossible to implement new models or investigate possible implementation errors in the code, and costs arise due to licensing requirements for commercial solvers. Although the open source framework OpenFOAM already includes multiple classes and two solvers dedicated to internal combustion engine simulations, there is no way to move engine valves and piston simultaneously with its standard tools. Thus, this paper presents a new engine library for ICE simulations written for OpenFOAM. The new framework is capable of simulating a complete fired engine cycle. The piston and the valves are moved simultaneously. To address large deformations in the mesh, a methodology to avoid insufficient mesh quality was developed. Ignition and combustion is modeled with standard tools from OpenFOAM. To validate the method, the simulation results for the averaged in-cylinder quantities pressure, temperature and mass are compared with experimental data.
Abstract High-accuracy simulations of internal combustion engines (ICE) allow deep insight into the physical processes of the different phases of the engine cycle: gas exchange, mixture [...]
D. Santos Serrano, F. Trias Miquel, G. Colomer Rey, C. Pérez Segarra
eccomas2022.
Abstract
Preservation of energy is fundamental in order to avoid the introduction of unphysical energy that can lead to unstable simulations. In this work, an energy-preserving unconditionally stable fractional step method on collocated grids is presented as a method which guarantees both preservation of energy and stability of our simulation. Using an algebraic (matrix-vector) representation of the classical incompressible Navier-Stokes equations mimicking the continuous properties of the differential operators, conservation of energy is formally proven. Furthermore, the appearence of unphysical velocities in highly distorted meshes is also adressed. This problem comes from the interpolation of the pressure gradient from faces to cells in the velocity correction equation, and can be corrected by using a proper interpolation.
Abstract Preservation of energy is fundamental in order to avoid the introduction of unphysical energy that can lead to unstable simulations. In this work, an energy-preserving unconditionally [...]
The formulation of multiphase flows emanates from basic conservation laws: mass, momentum and energy. While these are embedded in the celebrated Navier-Stokes equations, none of these properties do necessarily hold when constructing a computational model, unless special care is taken in discretizing the different terms of the governing equations. The conservation of both primary (mass, momentum) and secondary (energy) quantities is not only relevant to mimic the dynamics of the system, but also computationally beneficial. Conservation of such quantities produce an enhanced physical reliability, removing most of the need for stabilization artifacts. In addition, discrete conservation implies numerical stability as well, producing inherently stable problems. Focusing on the capillary force, which is one of the most distinguishable features of multiphase flows, we present here our most recent developments in the quest for conservation. Departing from an inherently mass conservative method, in this work we sketch our previous developments to obtain an energy conservation and next we present our attempt at momentum. By carefully assessing the continuum formulation, we delve into the mathematical properties responsible for the conservation of linear momentum, which we then mimic in the regularized and discrete formulations.
Abstract The formulation of multiphase flows emanates from basic conservation laws: mass, momentum and energy. While these are embedded in the celebrated Navier-Stokes equations, none [...]
Nonlinear instabilities are one of the major problems in turbulence simulations. One reason behind this problem is the accumulation of aliasing errors produced by the discrete evaluation of the convective term. This can be improved by preserving the quadratic invariants in a discrete sense. However, another source of instabilities is the error due to an incorrect evolution of thermodynamic variables, such as entropy. An appropriate discretization of the energy equation is needed to address this issue. An analysis of the preservation properties of various discretizations of the compressible Euler equations is reported, which includes some of the most common approaches used in the literature, together with some new formulations. Two main factors have been identified and studied: one is the choice of the energy equation to be directly discretized; the other is the particular splitting of the convective terms, chosen among the Kinetic Energy Preserving (KEP) forms. The energy equations analyzed in this paper are total and internal energy, entropy, and speed of sound. All the cases examined are locally conservative and KEP, since this is considered an essential condition for a robust simulation. The differences among the formulations have been theoretically investigated through the study of the discrete evolution equations induced by the chosen energy variable, showing which quantities may be preserved. Both one-dimensional and two-dimensional tests have been performed to assess the advantages and disadvantages of the various options in different cases.
Abstract Nonlinear instabilities are one of the major problems in turbulence simulations. One reason behind this problem is the accumulation of aliasing errors produced by the discrete [...]
A. Colombo, A. Ghidoni, E. Mantecca, G. Noventa, S. Rebay, D. Pasquale
eccomas2022.
Abstract
A high-order Discontinuous Galerkin (DG) solver is assessed in the computation of the flow through an Organic Rankine Cycle turbine nozzle and stage. The flow features are predicted with a RANS (Reynolds averaged NavierStoke) approach and the k-log() turbulence model in a multi reference frame, where interfaces between fixed and rotating zones are treated with a mixing plane approach, and non reflecting boundary conditions are used. Primitive variables based on pressure and temperature logarithms are adopted to ensure non-negative thermodynamic variables at a discrete level. The fluid can be modeled with the polytropic ideal gas law and the Peng-Robinson equation of state.
Abstract A high-order Discontinuous Galerkin (DG) solver is assessed in the computation of the flow through an Organic Rankine Cycle turbine nozzle and stage. The flow features are [...]
A. Colanera, A. Della Pia, M. Acquaviva, M. Chiatto, L. de Luca
eccomas2022.
Abstract
Modal analysis of three-dimensional gravitational thin liquid sheet flows, interacting with unconfined gaseous environments located on both sides of the liquid phase, is performed in the present work. Numerical data of this relevant two-phase flow configuration are obtained through the single-phase formulation and the Volume-of-Fluid (VOF) technique implemented in the flow solver Basilisk. This class of flows exhibits a variety of spatial and temporal relevant structures, both in free and forced configurations, that are investigated through the Spectral Proper Orthogonal Decomposition (SPOD). By means of such methodology, we explore the effect of two main governing parameters on the flow dynamics, namely the liquid sheet aspect ratio, AR = W/H, where H and W are the sheet inlet thickness and width, and the Weber number, We = lU2H/(2), in which U is the inlet liquid velocity, lthe liquid density, and the surface tension coefficient. Finally, for the highest aspect ratio value considered (AR = 40), we investigate the forced dynamics of the system excited by a harmonic perturbation in transverse velocity component applied at the inlet section, comparing results with ones arising from a purely two-dimensional analysis of the flow. The obtained results highlight the low rank behavior exhibited by the flow, suggesting that reduced order modeling could be particularly appealing to reduce complexity and computational effort in numerical simulation of this class of flows.
Abstract Modal analysis of three-dimensional gravitational thin liquid sheet flows, interacting with unconfined gaseous environments located on both sides of the liquid phase, is performed [...]