title={Spectral/hp large eddy simulation of vortex-dominated automotive flows around bluff bodies with diffuser and complex front wing geometries},
author={Fabian Buscariolo, Filipe},
year={2019},
school={Imperial College London},
groups={thesis}
abstract={Inthisresearchproject,itisdemonstratedtheuseofspectral/hpelementmethodforsimulationsoffully3Dcomplexgeometries.SuchsolutionsathighReynoldsnumbersandwithhigherorderpolynomialswerepreviouslyintractableduetonumericalstabilityissuesaffectingtheconvergenceofthescheme.Forthisapproach,wehaveemployedthelatestdevelopmentofcontinuousGalerkinspectralvanishingviscosity(CG-SVV)with a discontinuous Galerkin (DG) mimicking kernel. Together with dealiasing techniques, the numerical stability and convergence characteristics of the spectral/hp element method have been greatly improved. These advances in numerical methods are also supported by novel meshing strategies, taking advantage of the additional flexibility in changing the uniform polynomial orders of the mesh and the solution. As a result, efficient simulations can be formulated with consistent and highly accurate solutions obtainable. Specific for this work, the focus is on complex geometries often found in automotive engineering. To reduce the computational demands, this research explored the use of symmetry boundary conditions for large eddy simulations (LES) using a half model. It is found that if only the average flow properties near the body are of interest, such an approach can provide more than 50\% reduction in simulation time while maintaining the solution quality. In terms of improving the solution resolution, as one might expect from a p-type method we have observed that increasing the polynomial order can be a more effective approach in comparison to conventional mesh refinement. In the three test cases, we have successfully exploited the use of polynomial accuracy of 4th, 5th and 6th order. This is the first comprehensive study using polynomials of such high orders, and the corresponding solutions are obtained for fully 3D geometries using spectral/hp element method. Three test cases have been considered, the first being the simulations of the original Ahmed Body serves as a validation study for 3D simulations of the spectral/hp element method. The Ahmed Body is one of the most widely studied bluff bodies used for automotive conceptual studies and computational fluid dynamics (CFD) software validation. For this validation study, the differences in results obtained using various polynomial orders for the mesh as well as for the solution interpolation have been examined in detail. With the proposed approach, we were able to obtain fairly good correlations with the aerodynamic quantities for polynomial orders of 5 and above. Regarding the flow features around the body, solutions from the 6th order polynomial showed clear advantage in the slant vortex intensity. With the computational facilities, further increase of solution polynomial order is not feasible; however, the required solution resolution can also be obtained via the use of local mesh refinement. We determine that this level of solution accuracy, after comparing with various studies in the literature, cannot be obtainable using steady-state simulations such as the very popular Reynolds averaged Navier-Stokes (RANS) method. Based on the validation result, the second test case involved the simulations of Ahmed Body geometries with a simplified diffuser using the proposed method. This case serves as an independent study examining the suitability of the method for design analysis. Using the same 6th order polynomial and Refined mesh, the solution successfully identified the flow features consistent with to past literature on underbody diffusers. Additionally, we have found that the geometry of the reference body imposes a quite significant influence on the performance of the diffuser, as well as identified some strong interplay between the lower-side vortex and the diffuser flow. The toolchain has clearly demonstrated its capability in assisting integrated design analysis for a simplified road vehicle equipped with a diffuser. In the final test case, a new benchmark study case for aerodynamic design of high-performance vehicles and racing cars, the Imperial front wing is proposed. This study consists of a multi-element front-wing based on a Formula One front wing design. It generates complex flow features including ground effects, and multiple vortex system development and interaction. We used this test case as a challenging examination of our proposed method and simulation strategy using the spectral/hp element method. The simulations were also supported by an independent experimental study and results obtained for comparison achieved a high level of agreement. Using a polynomial order of 4th and above have successfully correlated the flow velocity fields at various planes downstream, while increasing the polynomial order to 5th will further result in a good matching of flow visualisation details. From all three test cases, the spectral/hp element method when applied to suitable meshes at reasonably polynomial orders has been able to accurately and consistently yield reliable solutions in good agreement with experiment. The benefits of using high order polynomials for mesh generation of complex geometries, and for solution interpolation of higher accuracy have enabled the use of much coarser meshes than would typically be applied in commercial CFD codes. The progress made in this research is a solid step forward for the adaptation of the spectral/hp element for industrial level applications.}
author={Dimitrios S. Lampropoulos, Ioannis D. Boutopoulos, George C. Bourantas, Karol Miller, Petros E. Zampakis and Vassilios C. Loukopoulos},
author={Dimitrios S Lampropoulos and Ioannis D. Boutopoulos and George C. Bourantas and Karol Miller and Petros E. Zampakis and Vassilios C. Loukopoulos},
title={Hemodynamics of anterior circulation intracranial aneurysms with daughter blebs: investigating the multidirectionality of blood flow fields},
journal={Computer Methods in Biomechanics and Biomedical Engineering},
author = {R.C. Moura and L.D. Fernandes, and A.F. C. da Silva and S. J. Sherwin},
doi = {10.1016/j.cma.2024.117025},
groups = {core},
abstract = {We present a new linear eigensolution analysis technique that provides superior estimates of dissipation distribution in wavenumber space for the continuous Galerkin (CG) method. The technique builds upon traditional dispersion–diffusion analyses that have been applied to spectral/hp element methods, but in particular is an improvement upon the non-modal eigenanalysis approach proposed by Fernandez et al. (2019). The present technique takes into account the indirect effects that dispersion may have on dissipation, as recently discussed by Moura et al. (2022), in order to better represent dissipation itself. Also, a concept used by the dynamic mode decomposition (DMD) community is invoked to weight the relative contribution of the multiple diffusion curves that stem from temporal eigenanalysis. This allows for obtaining a single dissipation profile in wavenumber space, so that the proposed technique is named joint-mode analysis. Although the non-modal approach also provides a single diffusion curve, the joint-mode dissipation curve is shown to correlate significantly better with the energy spectrum of Burgers’ turbulence at large and intermediate scales, which is particularly relevant for implicit large-eddy simulation (LES). The proposed technique is readily extensible to other spectral/hp element methods.}
