Abstract

The aim of this work is the study of fluid dynamics models using the CFD software
OpenFOAM, an open source software allowing meshing, manipulation, simulation
and post-processing of many problems involving fluid mechanics.
The work consists of a study with OpenFOAM of a real engineering problem, namely
to analyze hemodynamics in the thoracic aorta in collaboration with CIMNE (Centre
Internacional de Metodes Numerics a l'Enginyeria) and LABSON-UPC (Laboratorio de
Sistemas Oleohidricos y Neumcos). Speci cally, the study aims to compute the
shear stress that blood causes to aorta walls.

Abstract

The paper presents a new computational framework for the numerical simulation of fast large strain solid dynamics, with particular emphasis on the treatment of near incompressibility. A complete set of first order hyperbolic conservation equations expressed in terms of the linear momentum and the minors of the deformation (namely the deformation gradient, its co-factor and its Jacobian), in conjunction with a polyconvex nearly incompressible constitutive law, is presented. Taking advantage of this elegant formalism, alternative implementations in terms of entropy-conjugate variables are also possible, through suitable symmetrisation of the original system of conservation variables. From the spatial discretisation standpoint, modern Computational Fluid Dynamics code “OpenFOAM” [http://www.openfoam.com/] is here adapted to the field of solid mechanics, with the aim to bridge the gap between computational fluid and solid dynamics. A cell centred finite volume algorithm is employed and suitably adapted. Naturally, discontinuity of the conservation variables across control volume interfaces leads to a Riemann problem, whose resolution requires special attention when attempting to model materials with predominant nearly incompressible behaviour (). For this reason, an acoustic Riemann solver combined with a preconditioning procedure is introduced. In addition, a global a posteriori angular momentum projection procedure proposed in Haider et al. (2017) is also presented and adapted to a Total Lagrangian version of the nodal scheme of Kluth and Després (2010) used in this paper for comparison purposes. Finally, a series of challenging numerical examples is examined in order to assess the robustness and applicability of the proposed methodology with an eye on large scale simulation in future works.

Abstract

The computational cost of parametric studies currently represents the major limitation to the application of simulation-based engineering techniques in a daily industrial environment. This work presents the first nonintrusive implementation of the proper generalised decomposition (PGD) in OpenFOAM, for the approximation of parametrised laminar incompressible Navier–Stokes equations. The key feature of this approach is the seamless integration of a reduced order model (ROM) in the framework of an industrially validated computational fluid dynamics software. This is of special importance in an industrial environment because in the online phase of the PGD ROM the description of the flow for a specific set of parameters is obtained simply via interpolation of the generalised solution, without the need of any extra solution step. On the one hand, the spatial problems arising from the PGD separation of the unknowns are treated using the classical solution strategies of OpenFOAM, namely the semi-implicit method for pressure linked equations (SIMPLE) algorithm. On the other hand, the parametric iteration is solved via a collocation approach. The resulting ROM is applied to several benchmark tests of laminar incompressible Navier–Stokes flows, in two and three dimensions, with different parameters affecting the flow features. Eventually, the capability of the proposed strategy to treat industrial problems is verified by applying the methodology to a parametrised flow control in a realistic geometry of interest for the automotive industry.

Abstract

This work presents a novel proposal of a second-order accurate (in time and space) particle-based method for solving transport equations including incompressible flows problems within a mixed Lagrangian–Eulerian formulation. This methodology consists of a symmetrical operator splitting, the use of high-order operators to transfer data between the particles and the background mesh, and an improved version of the eXplicit Integration Following the Streamlines (X–IVS) method. In the case of incompressible flows, a large time-step iterative solver is employed where the momentum equation is split to improve the numerical approximation of the convective term.

New interpolation and projection operators are evaluated and quadratically accurate solutions of scalar transport tests are presented. Then, incompressible flow problems are solved where the rate of convergence of the method is assessed using both structured and unstructured background grids. The method is implemented in the open source platform OpenFOAM®allowing employing arbitrary meshes and obtaining reliable computing time comparisons with standardized solvers. The results obtained reveal that the current method is able to obtain a lower level of error than a fast Eulerian alternative, without increasing the total computing time.

Abstract

Abstract

This work presents a numerical study on the effect of flame-wall interaction (FWI) from the viewpoint of flame dynamics. For that purpose, direct numerical simulations (DNS) employing detailed calculations of reaction rates and transport coefficients have been applied to a 2D premixed methane/air flame under atmospheric condition. Free flame (FF) and side-wall quenching (SWQ) configurations are realized by defining one lateral boundary as either a symmetry plane for the FF or a cold wall with fixed temperature at 20 oC for the SWQ case. Different components of flame stretch and Markstein number regarding tangential, normal (due to curvature) and total stretch, Kas , Kac and Katot = Kas + Kac, as well as their correlations with respect to the local flame consumption speed SL have been evaluated. It has been shown that the FWI zone is dominated by negative flame stretch. In addition, SL decreases with decreasing normal stretch due to curvature Kac while approaching the cold wall. However, SL increases with decreasing Kac while approaching the symmetry boundary for the free flame case, leading to an inversion of the Markstein number Matot based on Katot from positive in the free flame case to negative in the SWQ case. The quenching distance evaluated based on wall-normal profiles of SL has been found to be approximately equal to the unstretched laminar flame thickness, which compares quantitatively well with measured data from literature. The flame speed has been confirmed to scale quasi-linearly with the stretch in the FWI zone. The results reveal a distinct correlation during transition between FWI and FF regarding flame dynamics, which brings a new perspective for modeling FWI phenomena by means of flame stretch and Markstein number. To do this, the quenching effect of the wall may be reproduced by a reversed sign of the Markstein number from positive to negative in the FWI zone and by applying the general linear Markstein correlation (SL/SL,0 = 1− Ma ·Ka), leading to a decrease of the flame speed or the reaction rate in the near-wall region.

Abstract

Abstract

Abstract

Abstract

In the present study we investigate an incompressible turbulent channel flow with heat transfer at Re_τ = 180 with a deterministic surface topography consisting of truncated cones. Two solvers for each of the two boundary handling strategies are considered. With Nek5000 and OpenFOAM the influence of the roughness elements is directly accounted for by an unstructured body fitted mesh, whereas Xcompact3d and SIMSON utilize the immersed boundary method (IBM) to deal with the 3D geometry. The main focus of this work is on an evaluation of the usability of the IBM and a comparison of the parallel performance of the different solvers. Since usability is an ambiguous definition, various quantities are compared: global statistics like Nusselt number and friction coefficient, one-dimensional wall-normal profiles for first and second order statistics, as well as three-dimensional averages over roughness sections. In addition, the computational effort for each method is documented.