Documents published in Scipedia

• An Immersed Boundary Method for the CFD Solver Airbus-CODA

  V. Llorente, D. Lodares, E. Ferrer, E. Valero

  eccomas2022.

  
  

  Abstract

  This work implements and analyses an Immersed Boundary Method based on Volume Pezalization for the flow simulator Airbus-CODA (CFD for ONERA, DLR, and AIRBUS). The Immersed Boundary Volume Penalization has unique advantages, e.g. easy to implement, straightforward formulation for moving geometries, and numerical errors can be controlled apriori [1, 2], showing the potential for aeronautical applications. Numerical experiments will assess the accuracy of the Immersed Boundary Volume Penalization in CODA.

  Abstract
  This work implements and analyses an Immersed Boundary Method based on Volume Pezalization for the flow simulator Airbus-CODA (CFD for ONERA, DLR, and AIRBUS). The Immersed [...]

  • 33
  •  read

• Momentum-conserving ROMs for the incompressible Navier-Stokes equations

  H. Rosenberger, B. Sanderse

  eccomas2022.

  
  

  Abstract

  Projection-based model order reduction of an ordinary differential equation (ODE) results in a projected ODE. Based on this ODE, an existing reduced-order model (ROM) for finite volume discretizations satisfies the underlying conservation law over arbitrarily chosen subdomains. However, this ROM does not satisfy the projected ODE exactly but introduces an additional perturbation term. In this work, we propose a novel ROM with the same subdomain conservation properties which also satisfies the perturbed ODE exactly. We apply this ROM to the incompressible Navier-Stokes equations and show with regard to the mass equation how the novel ROM can be constructed to satisfy algebraic constraints. Furthermore, we show that the resulting mass-conserving ROM allows us to derive kinetic energy conservation and consequently nonlinear stability, which was not possible for the existing ROM due to the presence of the perturbation term.

  Abstract
  Projection-based model order reduction of an ordinary differential equation (ODE) results in a projected ODE. Based on this ODE, an existing reduced-order model (ROM) for [...]

  • 4
  •  read

• Numerical simulation of the micro-extrusion process of printable biomaterials

  A. Amani, D. Kizildag, J. Castro, L. del Mazo, M. Pegueroles, M. Ginebra

  eccomas2022.

  
  

  Abstract

  This work aims to gain a better understanding of how the rheological properties of printable materials affect their processability, as well as the quality of the final product, which at the end can lead to reducing time and costs of the process and increase product development. As the first step, the proper rheological non-Newtonian models are extracted from experimental studies. Later, three-dimensional numerical simulation of extrusion process is performed in the context of Direct Numerical Simulation (DNS) of governing equations, where the whole physics of fluid motion is taken into account. A finite-volume fractional step approach is used to solve the Navier-Stocks equations on collocated arbitrary meshes. Geometrical volume-of-fluid (GVOF) interface capturing approach is used to resolve the topological changes of the moving interface. The governing equations are solved using High-Performance Computing (HPC) parallel approaches. Besides the contribution of this work to the advancement of numerical techniques applied to multiphase complex flows, obtained results will shed light on the nature of non-Newtonian extrusion process with vast applications in the 3D printer industrial sectors.

  Abstract
  This work aims to gain a better understanding of how the rheological properties of printable materials affect their processability, as well as the quality of the final product, [...]

  • 1
  •  read

• The X-Mesh method applied to Multiphase Flows

  Q. Antoine, J. Remacle, J. Lambrechts, N. Moes

  eccomas2022.

  
  

  Abstract

  The accurate modeling of moving boundaries and interfaces is a difficulty present in many situations in computational mechanics. In this paper we use a new approach, X-Mesh, to simulate with the finite element method the interaction between two immiscible fluids while keeping an accurate description of the interface without mesh regeneration. The method is validated with complex problems such as Rayleigh-Taylor instabilities, sloshing and dambreak. The quality of the results and the efficiency of the method show the potential of this approach to simulate such physical phenomena.

  Abstract
  The accurate modeling of moving boundaries and interfaces is a difficulty present in many situations in computational mechanics. In this paper we use a new approach, X-Mesh, [...]

  • 5
  •  read

• DNS and LES on unstructured grids: playing with matrices to preserve symmetries using a minimal set of algebraic kernels

  F. Trias, X. Álvarez-Farré, À. Alsalti-Baldellou, A. Gorobets, A. Oliva

  eccomas2022.

  
  

  Abstract

  The essence of turbulence are the smallest scales of motion. They result from a subtle balance between two differential operators differing in symmetry: the convective operator is skew-symmetric, whereas the diffusive is symmetric and negative-definite. On the other hand, accuracy and stability need to be reconciled for numerical simulations of turbulent flows in complex configurations. With this in mind, a fully-conservative discretization method for collocated unstructured grids was proposed [Trias et al., J.Comp.Phys. 258, 246-267, 2014]: it preserves the symmetries of the differential operators and it has shown to be a very suitable approach for DNS and LES. On the other hand, an efficient cross-platform portability is nowadays one of the greatest challenges for CFD codes. In this regard, our leitmotiv reads: relying on a minimal set of (algebraic) kernels is crucial for code portability and maintenance! In this context, this work focuses on the computation of eigenbounds for the above-mentioned convection and diffusion matrices which are needed to determine the time-step `a la CFL. A new inexpensive method that allows this, without explicitly constructing these time-dependent matrices is proposed and tested. It only requires a sparse-matrix vector product where only the vector changes on time. Hence, apart from being significantly more efficient than the standard CFL condition, cross-platform portability is straightforward.

  Abstract
  The essence of turbulence are the smallest scales of motion. They result from a subtle balance between two differential operators differing in symmetry: the convective operator [...]

  • 24
  •  read

• The interaction between an isolated roughness element and free-stream turbulence

  M. Gholamisheeri, K. Durovic, S. Mamidala, J. Fransson, A. Hanifi, D. Henningson

  eccomas2022.

  
  

  Abstract

  Control and delay of the laminar-turbulent transition is a key parameter in reducing skin friction and drag. The flow characteristics, surface roughness, and environmental noise can affect the onset of transition. The present work investigates, numerically and experimentally, the interaction of the free-stream turbulence (FST) and an isolated cylindrical roughness element, and the resulting impact on the transition onset in a flat-plate boundary layer. High-fidelity direct numerical simulations (DNS) are performed for a roughness element immersed in the boundary layer over a flat plate with an asymmetrical leading edge, with and without FST. The numerical results are compared to hot-wire anemometry measurements performed in the Minimum Turbulence Level wind tunnel at KTH. The initial numerical and experimental results show that in the absence of FST, for the chosen flow parameters, highand low-speed streaks are generated downstream of the roughness element while the flow remains laminar and globally stable. When FST is added, the spanwise spacing of the streaky structures changes and the transition location of the boundary layer moves upstream. It was found that the aspect ratio of the streaky structures does not vary significantly.

  Abstract
  Control and delay of the laminar-turbulent transition is a key parameter in reducing skin friction and drag. The flow characteristics, surface roughness, and environmental [...]

  • 4
  •  read

• Non-newtonian viscous elongation and shear fluid model based on optimal triple tensor decomposition

  M. Rütten, T. Schomberg

  eccomas2022.

  
  

  Abstract

  The modeling of the distinct non-Newtonian fluid properties is an essential prerequisite for the computational simulation of associated flow fields. In particular, some non-Newtonian fluids reveal strong diverse viscosity response behaviours to pure elongational and simple shearing flows. Therefore, it is necessary to be able to distinguish between these flow types even in complex flow configuration. Unfortunately flow types are naturally mixed and this distinction becomes quite difficult. Only in a Lagrangian framework the tracking of the Lagrangian fluid element deformation allows an accurate strain related deformation type assignment. However, most CFD approaches prefer the Eulerian framework accepting the loss of the natural flow path alignment of the moving fluid particles. Consequently shear and elongation rates are barely separable without particular assignment methods. In this work a tensor decomposition method from vortex dynamics is discussed which allows to distinguish between these flow types. In vortex dynamics the problem occurred to separate shearing from purely rotational flows because different hydro- and aerodynamic flow phenomena are caused by shear and vortex related flow types. Thereto, various methods were proposed, among others the optimal triple tensor decomposition method which is able to separate vortical from shearing flows but also, after some modification, elongational from shearing flows. This tensor decomposition is now used to calculate elongation and shearing rates as input variables into non-Newtonian fluid models for the calculation of the local elongational and shear viscosity. The application case is a cross slot channel flow often used as reference. In this numerical simulation study the impact of the elongation rate modeling on the contraction flow topology is shown and discussed. It is shown that the modeling of different viscous elongational and shear-thinning affects the resulting flow significantly  

  Abstract
  The modeling of the distinct non-Newtonian fluid properties is an essential prerequisite for the computational simulation of associated flow fields. In particular, some non-Newtonian [...]

  • 1
  •  read

• Neural network prediction of the flow field in a periodic domain with hyper-neural network parametrization

  O. Bublík, V. Heidler, A. Pecka, J. Vimmr

  eccomas2022.

  
  

  Abstract

  This paper is concerned with fast flow field prediction in a blade cascade for variable blade shapes as well as variable Reynolds number using the machine-learning architecture called convolutional neural network. To generate flow field for a specific Reynolds number, an encoder-decoder convolutional neural network, also called U-Net, is used. The values 500, 1000 and 1500 of the Reynolds number are chosen as the training set. Three U-Nets were trained on CFD results for 100 blade profiles, each U-Net for a different Reynolds number. In order to get a prediction for variable Reynolds number, a so-called hypernetwork in employed. The hypernetwork essentially interpolates between the two trained U-Nets. The architecture of the hypernetwork is fully-connected feedforward neural network with one input neuron corresponding to the Reynolds number, one hidden layer and the output layer corresponds to the weights for the interpolated U-Net. The concept of the hypernetwork-based parametrization is tested on a problem of compressible fluid flow through a blade cascade with three unseen blade profiles and unseen Reynolds number.

  Abstract
  This paper is concerned with fast flow field prediction in a blade cascade for variable blade shapes as well as variable Reynolds number using the machine-learning architecture [...]

  • 1
  •  read

• Symmetry-preserving discretisation methods for magnetohydrodynamics

  J. Hopman, X. Trias, J. Rigola

  eccomas2022.

  
  

  Abstract

  In this work, the symmetry-preserving method [1, 2, 3] is extended to include magnetohydrodynamic effects, using the collocated grid arrangement of Ni et al. [4, 5]. The electromagnetic part is solved explicitly using the induction-less approximation and an electric potential Poisson equation. The proposed solver is implemented in OpenFOAM and tested for accuracy and stability, and compared to the method of Ni et al. [4, 5]. A new benchmark case using a Taylor-Green vortex in a transverse magnetic field is used, for which kinetic energy budget terms are compared to the analytical solutions. Finally, Hunt's case is used to compare flow profiles to the analytical solutions. Influence of the spatial discretisation on accuracy and stability is also examined by solving both cases on meshes with variable degrees of distortion. The symmetry-preserving method showed accuracy on Cartesian meshes and stability even on extremely distorted meshes, whereas the method of Ni et al. [4, 5] showed less accurate conservation of current density and was not able to produce stable solutions on the extremely distorted meshes.

  Abstract
  In this work, the symmetry-preserving method [1, 2, 3] is extended to include magnetohydrodynamic effects, using the collocated grid arrangement of Ni et al. [4, 5]. The electromagnetic [...]

  • 2
  •  read

• Fluidization by gas pore pressure of dense granular flows: numerical simulations versus experiments

  A. Aravena, L. Chupin, T. Dubois, O. Roche

  eccomas2022.

  
  

  Abstract

  A model accounting for fluidisation by pore gas pressure in dense granular flows is presented. A viscoplastic rheology, based on the Drucker-Prager criterium, is used to describe the granular medium which is a mixture of air and glass beads. The pore gas pressure, which satisfies an advection-diffusion equation, reduces the friction between the particles and thus the value of the apparent viscosity. As a consequence, dense fluidised granular flows can travel longer distances. In laboratory experiments, the run-out distance reached by dense granular columns when collapsing is almost doubled when fluidisation is applied. This fundamental result, in the context of pyroclastic density currents, is reproduced by numerical simulations performed with the fluidised model.

  Abstract
  A model accounting for fluidisation by pore gas pressure in dense granular flows is presented. A viscoplastic rheology, based on the Drucker-Prager criterium, is used to describe [...]

  • 14
  •  read