Publications
386 results found
Slaughter J, Moxey D, Sherwin S, 2023, Large Eddy Simulation of an Inverted Multi-element Wing in Ground Effect, Flow, Turbulence and Combustion, Vol: 110, Pages: 917-944, ISSN: 1386-6184
Due to the proprietary nature of modern motorsport and Formula 1, current scientific literature lacks relevant studies and benchmarks that can be used to understand flow physics in this area, as well as to test and validate new simulation methodology. With the release of a new, open-source geometry (the Imperial Front Wing), we present a computational study of a multi-element aerofoil at a ride height of 0.36h/c and a Reynolds number of 2. 2× 105. A 0.16c slice of the Imperial Front Wing has been examined using high-order spectral/hp element methods. Time averaged force data is presented, finding lift and drag coefficients of - 8.33 and 0.17 respectively. Unsteady analysis of the force and surface pressure data has allowed salient feature identification with respect to the transition mechanisms of each element. The mainplane and flap laminar separation are studied and the cross-spectral phase is presented for the lower frequency modes. At St=40 an in-phase relationship is identified between mainplane and flap laminar separation bubbles, whilst at St=60 a distinct out-of-phase relationship is observed. Wake results, including wake-momentum deficit and turbulent kinetic energy are presented, which show wake meandering and subsequent breakdown due to a Kelvin–Helmholtz instability. These results, in particular the transition mechanisms, will allow for the construction of a dataset to validate novel methods in this area.
Lyu G, Chen C, Du X, et al., 2023, Stable, entropy-pressure compatible subsonic Riemann boundary condition for embedded DG compressible flow simulations, Journal of Computational Physics, Vol: 476, Pages: 1-22, ISSN: 0021-9991
One approach to reducing the computational cost of simulating transitional compressible boundary layer flow is to adopt a near body reduced domain with boundary conditions enforced to be compatible with a computationally cheaper three-dimensional RANS simulation. In such an approach it is desirable to enforce a consistent pressure distribution which is not typically the case when using the standard Riemann inflow boundary conditions. We therefore revisit the Riemann problem adopted in many DG based high fidelity formulations.Through analysis of the one-dimensional linearized Euler equations it is demonstrated that maintaining entropy compatibility with the RANS simulation is important for a stable solution. It is also necessary to maintain the invariant for the Riemann outflow boundary condition in a subsonic flow leaving one condition that can be imposed at the inflow boundary. Therefore the entropy-pressure enforcement is the only stable boundary condition to enforce a known pressure distribution. We further demonstrate that all the entropy compatible inflow Riemann boundary conditions are stable providing the invariant compatible Riemann outflow boundary condition is also respected.Although the entropy-pressure compatible Riemann inflow boundary condition is stable from the one-dimensional analysis, two-dimensional tests highlight divergence in the inviscid problem and neutrally stable wiggles in the velocity fields in viscous simulations around the stagnation point. A two-dimensional analysis about a non-uniform baseflow assumption provides insight into this stability issue (ill-posedness) and motivates the use of a mix of inflow boundary conditions in this region of the flow.As a validation we apply the proposed boundary conditions to a reduced domain of a wing section normal to the leading-edge of the CRM-NLF model taken out of a full three-dimensional RANS simulation at Mach 0.86 and a Reynolds number of 8.5 million. The results show that the entropy-pressure
Lindblad D, Sherwin SJ, Cantwell C, et al., 2023, Large Eddy simulations of isolated and installed jet noise using the high-order discontinuous Galerkin method, AIAA SCITECH 2023 Forum, Publisher: American Institute of Aeronautics and Astronautics, Pages: 1-21
A recently developed computational framework for jet noise is used to compute the noise generated by an isolated and installed jet. The framework consists of two parts. In the first part, the spectral/hp element framework Nektar++ is used to compute the near-field flow. Nektar++ solves the unfiltered Navier-Stokes equations on unstructured grids using the high-order discontinuous Galerkin method. The discrete equations are integrated in time using an implicit scheme based on the matrix-free Newton-GMRES method. In the second part, the Antares library is used to compute the far-field noise. Antares solves the Ffowcs Williams - Hawkings equation for a permeable integration surface in the time domain using a source-time dominant algorithm. The simulations are validated against experimental data obtained in the Doak Laboratory Flight Jet Rig, located at the University of Southampton. For the isolated jet, good agreement is achieved, both in terms of the flow statistics and the far-field noise. The discrepancies observed for the isolated jet are believed to be caused by an under-resolved boundary layer in the simulations. For the installed jet, the flow statistics are also well predicted. In the far-field, very good agreement is achieved for downstream observers. For upstream observers, some discrepancies are observed for very high and very low frequencies.
Arshad M, Cheng S, van Reeuwijk M, et al., 2023, Modification of the swirling well cell culture model to alter shear stress metrics, BIOTECHNOLOGY AND BIOENGINEERING, ISSN: 0006-3592
Yan ZG, Pan Y, Peiró J, et al., 2023, Eigenspectral Analysis of Preconditioners in an Adaptive Compressible Flow Solver, Pages: 521-532, ISSN: 1439-7358
An implicit-in-time discontinuous Galerkin (DG) solver has been developed for compressible flows, which adopts an adaptive time stepping strategy balancing different discretization errors. To study the effects of preconditioning for further speed-up, preconditioners based on Jacobi iteration, Gauss-Seidel iteration, incomplete LU factorization, and p-multigrid are studied using eigenspectral analysis and numerical experiments. The condition number and eigenvalue distributions are studied based on a lid-driven flow. The preliminary results show that the eigenspectral analysis can correctly reflect preconditioning effects and the p-multigrid method is highly dependent on the smoother adopted.
Isler J, Vivarelli G, Montomoli F, et al., 2023, High Fidelity Compressible and Incompressible Flow Simulations of an Engine Intake Pressure Distribution
In this work, we investigated the impact of compressibility of an engine intake pressure distribution at high Reynolds numbers in absence of transonic effects and under strong adverse pressure gradients on a flate plate by means of computational simulations. In order to assess the compressibility effects in the vicinity and over the flat plate surface, compressible and incompressible high-order Implicit Large Eddy Simulations (iLES) of the Navier-Stokes (NS) equations were performed using the Nektar++ framework. With this methodology, we were able to assert what are the physical mechanisms behind the turbulent transition processes. In addition, and crucially, what role the compressibility is playing in this flow configuration, in order to promote the design of improved geometries using high-order methods. Therefore, a consistent investigation of the density variation effects and boundary layer behaviours of the compressible and incompressible solutions provided the necessary understanding of the viscous driven separation that governs the flow in an engine intake at moderate fan speed.
Vivarelli G, Isler JA, Williams TS, et al., 2023, ON THE EFFECT OF CURVATURE AND COMPRESSIBILITY IN LAMINAR BOUNDARY LAYERS OVER FAN BLADES
During the early phases of the turbomachinery design process, it is often the case that various simplifications and assumptions are made to understand boundary layer behaviour subject to pressure gradients. For example, simplified boundary-layer models can be employed to appreciate losses incurred in flows over gas turbine components. These neglect the effect of surface curvature and/or density, while retaining the surface pressure distribution. This paper's objective is to assess the implication of those simplifications by studying a realistic fan loading and the status of the boundary layer just upstream of a shock. A set of 2D aerofoils, representing the pressure distribution found at 70% and 90% span of a modern low-speed fan, are analysed at cruise conditions. In order to investigate compressibility, the two profiles were redesigned to achieve exactly the same pressure loading at incompressible conditions. On the other hand, curvature effects were studied by means of a comparison of incompressible flow over the symmetric aerofoil and a convergent-convergent nozzle, once again with the same loading. The resulting laminar boundary layer quantities were compared explaining the necessary scaling required along with the reasons behind the discrepancies, such as the effect of density and curvature variations. All test-cases were analysed deploying Nektar++. This is a high-order spectral h/p element solver having both compressible and incompressible formulations.
Marioni YF, Adami P, Montomoli F, et al., 2023, MACHINE-LEARNT TURBULENCE CLOSURES FOR AXIAL COMPRESSOR CASCADE WITH CORNER SEPARATION
Corner separation in axial compressors is a complex phenomenon, which is hardly well captured by traditional steady RANS calculations, primarily because of the deficiencies of the Reynolds stress tensor formulations. In this work a Machine Learning (ML) framework is applied to Large-Eddy Simulation (LES) data to develop non-linear turbulence stress closures. Two linear compressor cascade LES calculations are run with the Rolls-Royce solver HYDRA: one at nominal incidence condition, for which no corner separation is observed, and one at high incidence, where more complex secondary flow structures and recirculation appear. The two cases are validated against experiments performed in the linear cascade facility at LMFA, Lyon. A transitional variant of Wilcox's k−ω SST is used as baseline turbulence model and traditional closures are found to perform well at nominal conditions, but poorly at higher incidence, as they strongly over-predict losses and secondary flows. The coefficients of an Explicit Algebraic Reynolds Stress Model (EARSM) are trained using Artificial Neural Networks (ANN) on the high incidence dataset. When tested in HYDRA, improvements are observed in the prediction of end-wall losses, as well as vorticity and Reynolds stress contours downstream of the separation region. The model also performs well at nominal incidence. The importance of a near-wall coefficient damping is discussed. Finally, end-wall loss polars are computed and compared for different non-linear constitutive relationships.
Lahooti M, Bao Y, Scott D, et al., 2023, LES/DNS fluid-structure interaction simulation of non-linear slender structures in Nektar plus plus framework, COMPUTER PHYSICS COMMUNICATIONS, Vol: 282, ISSN: 0010-4655
Moura RC, Fernandes LD, Silva AFC, et al., 2022, Spectral/hp element methods' linear mechanism of (apparent) energy transfer in Fourier space: Insights into dispersion analysis for implicit LES, JOURNAL OF COMPUTATIONAL PHYSICS, Vol: 471, ISSN: 0021-9991
- Author Web Link
- Cite
- Citations: 1
Chun S, Marcon J, Peiro J, et al., 2022, Reducing errors caused by geometrical inaccuracy to solve partial differential equations with moving frames on curvilinear domain, COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, Vol: 398, ISSN: 0045-7825
- Author Web Link
- Cite
- Citations: 1
Lyu G, Chen C, Du X, et al., 2022, Open-source Framework for Transonic Boundary Layer Natural Transition Analysis over Complex Geometries in Nektar++, AIAA Aviation 2022 Forum
Lindblad D, Sherwin S, Cantwell C, et al., 2022, Aeroacoustic analysis of a subsonic jet using the discontinuous Galerkin method, 28th AIAA/CEAS Aeroacoustics 2022 Conference, Publisher: American Institute of Aeronautics and Astronautics, Pages: 1-21
In this work, the open-source spectral/hp element framework Nektar++ is coupled with the Antares library to predict noise from a subsonic jet. Nektar++ uses the high-order discontinuous Galerkin method to solve the compressible Navier-Stokes equations on unstructured grids. Unresolved turbulent scales are modeled using an implicit Large Eddy Simulation approach. In this approach, the favourable dissipation properties of the discontinuous Galerkin method are used to remove the highest resolved wavenumbers from the solution. For time-integration, an implicit, matrix-free, Newton-Krylov method is used. To compute the far-field noise, Antares solves the Ffowcs Williams - Hawkings equation for a permeable integration surface in the time-domain using a source-time dominant algorithm. The simulation results are validated against experimental data obtained in the Doak Laboratory Flight Jet Rig, located at the University of Southampton.
Son O, Gao A, Gursul I, et al., 2022, Leading-edge vortex dynamics on plunging airfoils and wings, Journal of Fluid Mechanics, Vol: 940, Pages: 1-30, ISSN: 0022-1120
The vortex dynamics of leading-edge vortices on plunging high-aspect-ratio (AR = 10) wings and airfoils were investigated by means of volumetric velocity measurements, numerical simulations, and stability analysis in order to understand the deformation of the leading-edge vortex filament and spanwise instabilities. The vortex filaments on both the wing and airfoilexhibit spanwise waves, but with different origins. The presence of a wing tip causes the leg of the vortex to remain attached to the wing upper surface, while the initial deformation of the filament near the wing-tip resembles a helical vortex. The essential features can be modelled as the deformation of initially L-shaped semi-infinite vortex column. In contrast, the instabilityof the vortices is well captured by the instability of counter-rotating vortex pairs, which are formed either by the trailing-edge vortices or the secondary vortices rolled-up from the wing surface. The wavelengths observed in the experiments and simulations are in agreement with the stability analysis of counter-rotating vortex pairs of unequal strength.
Pan Y, Yan Z-G, Peiro J, et al., 2022, Development of a balanced adaptive time-stepping strategy based on an implicit JFNK-DG compressible flow solver, Communications on Applied Mathematics and Computation, Vol: 4, Pages: 728-757, ISSN: 2661-8893
A balanced adaptive time-stepping strategy is implemented in an implicit discontinuous Galerkin solver to guarantee the temporal accuracy of unsteady simulations. A proper relation between the spatial, temporal and iterative errors generated within one time step is constructed. With an estimate of temporal and spatial error using an embedded Runge-Kutta scheme and a higher order spatial discretization, an adaptive time-stepping strategy is proposed based on the idea that the time step should be the maximum without obviously influencing the total error of the discretization. The designed adaptive time-stepping strategy is then tested in various types of problems including isentropic vortex convection, steady-state flow past a flat plate, Taylor-Green vortex and turbulent flow over a circular cylinder at Re=3900. The results indicate that the adaptive time-stepping strategy can maintain that the discretization error is dominated by the spatial error and relatively high efficiency is obtained for unsteady and steady, well-resolved and under-resolved simulations.
Xu H, Tu G, Sherwin SJ, 2022, Theoretical advances and applications of high-fidelity computation and modelling in fluid dynamics, COMPUTERS & FLUIDS, Vol: 241, ISSN: 0045-7930
Hambli W, Slaughter J, Buscariolo FF, et al., 2022, Extension of Spectral/hp Element Methods towards Robust Large-Eddy Simulation of Industrial Automotive Geometries, FLUIDS, Vol: 7
Cassinelli A, Mateo Gabín A, Montomoli F, et al., 2022, Reynolds sensitivity of the wake passing effect on a LPT cascade using spectral/hp element methods, International Journal of Turbomachinery, Propulsion and Power, Vol: 7, Pages: 8-8, ISSN: 2504-186X
Reynolds-Averaged Navier–Stokes (RANS) methods continue to be the backbone of CFD-based design; however, the recent development of high-order unstructured solvers and meshing algorithms, combined with the lowering cost of HPC infrastructures, has the potential to allow for the introduction of high-fidelity simulations in the design loop, taking the role of a virtual wind tunnel. Extensive validation and verification is required over a broad design space. This is challenging for a number of reasons, including the range of operating conditions, the complexity of industrial geometries and their relative motion. A representative industrial low pressure turbine (LPT) cascade subject to wake passing interactions is analysed, adopting the incompressible Navier–Stokes solver implemented in the spectral/hp element framework Nektar++. The bar passing effect is modelled by leveraging a spectral-element/Fourier Smoothed Profile Method. The Reynolds sensitivity is analysed, focusing in detail on the dynamics of the separation bubble on the suction surface as well as the mean flow properties, wake profiles and loss estimations. The main findings are compared with experimental data, showing agreement in the prediction of wake traverses and losses across the entire range of flow regimes, the latter within 5% of the experimental measurements.
Liu B, Cantwell CD, Moxey D, et al., 2022, VECTORISED SPECTRAL/HP ELEMENT MATRIX-FREE OPERATOR FOR ANISOTROPIC HEAT TRANSPORT IN TOKAMAK EDGE PLASMA
A highly efficient matrix-free Helmholtz operator with single-instruction multipledata (SIMD) vectorisation is implemented in Nektar++ [1] and applied to the simulation of anisotropic heat transport in tokamak edge plasma. A tokamak is currently the leading candidate for a practical fusion reactor using the magnetic confinement approach to produce electricity through controlled thermonuclear fusion. Predicting the transport of heat in magnetized plasma is important to designing a safe tokamak design. Due to the ionized nature of plasma, the heat conduction of the magnetized plasma is highly anisotropic along the magnetic field lines. In this study, a variational form is proposed to simulate the anisotropic heat transport in magnetized plasma and the details of its mathematical derivation and implementation are presented. To accurately approximate the thermal load of plasma deposition on the wall of tokamak chamber, highly scalable and efficient algorithms are crucial. To achieve this, a matrix-free Helmholtz operator is implemented in the Nektar++ framework, utilising sum-factorisation to reduce the operation count and increase arithmetic intensity, and leveraging SIMD vectorisation to accelerate the computation on modern hardware. The performance of the implementation is assessed by measuring throughput and speed-up of the operators using deformed and regular quadrilateral and triangular elements.
Marioni YF, Cassinelli A, Adami P, et al., 2022, DEVELOPMENT OF MACHINE-LEARNT TURBULENCE CLOSURES FOR WAKE MIXING PREDICTIONS IN LOW-PRESSURE TURBINES
In this work, a DNS – Machine Learning (ML) framework is developed for low-pressure turbine (LPT) profiles to inform turbulence closures in Reynolds-Averaged Navier–Stokes (RANS) calculations. This is done by training the coefficients of Explicit Algebraic Reynolds Stress Models (EARSM) with shallow artificial neural networks (ANN) as a function of input flow features. DNS data are generated with the incompressible Navier–Stokes solver in Nektar++ and validated against experiments. All calculations include moving bars upstream of the profile to capture the effect of incoming wakes. The resulting formulations are then implemented in the Rolls-Royce solver HYDRA and tested a posteriori. The aim is to improve mixing predictions in LPT wakes, compared to the baseline model, Wilcox’s k− ω SST, in terms of velocity profiles, turbulent kinetic energy (TKE) production and mixing losses. LPT calculations are run at Reynolds numbers spanning from ≈ 80k to ≈ 300k, to cover the range of aircraft engine applications. Models for the low and high Reynolds datasets are trained separately and a method is developed to merge the two together. The resulting model is tested on an intermediate Reynolds case. This process is followed for two computational domains: one starting downstream of the profile trailing edge and one including the last portion of the profile. Finally, the developed closures are tested on the entire profile, to confirm the validity of the improvements when the additional effect of transition is included in the simulation. This work explains the methodology used to develop ML-driven closures and shows how it is possible to combine models trained on different datasets.
Vivarelli G, Isler JA, Montomoli F, et al., 2022, HIGH-ORDER SPECTRAL/HP COMPRESSIBLE AND INCOMPRESSIBLE COMPARISON OF TRANSITIONAL BOUNDARY-LAYERS SUBJECT TO A REALISTIC PRESSURE GRADIENT AND HIGH REYNOLDS NUMBER
Within the literature, there are limited high-order results concerning large Reynolds number flows under the influence of strong adverse pressure gradients, mainly due to the computational expense involved. The main advantage in employing high-order methodologies over standard second-order finite-volume solvers, relates to their ability to increase accuracy with a significantly lower number of degrees of freedom. In theory, this would permit Direct Numerical Simulation sort of analysis. Yet, there is still a significant computational cost involved. For this reason, an efficient approach to analyse such flows by means of a Nektar++ high-order Implicit Large Eddy Simulation is proposed. The flow conditions considered in this case cause a separation bubble to form with consequent turbulent transition. In particular, Tollmien-Schlichting instabilities trigger Kelvin-Helmholtz behaviour, which in turn cause the turbulent transition. The bulk of the study will be carried out with the incompressible flow solver, as it is assumed that compressibility effects are negligible within the boundary layer. An initial 2D analysis will be conducted to determine the necessary spatial resolution and whether it is possible to consider a subset of the overall simulation domain to reduce the computational expense. Once this will have been established, the 3D results will be achieved by Fourier expansion in the cross-flow direction. These results will prove the cost-effectiveness of the methodology, that could be used within an industrial setting with a limited turn-around time. Additionally, a comparison between the results achieved by means of the Nektar++ compressible flow solver in 2D and 3D will be provided, to assess any differences that may be present.
Gao AK, Sherwin SJ, Cantwell CD, 2022, THREE-DIMENSIONAL TRANSITION OF A LOW REYNOLDS NUMBER FLOW PAST A PLUNGING NACA 0012 AIRFOIL AT POST-STALL ANGLE OF ATTACK
The two-dimensional to three-dimensional transition of a flow past a plunging NACA 0012 airfoil at a Reynolds number of Re = 400, based on the chord length c, and an angle of attack of 15 degrees was investigated using global linear stability analysis and spanwise-homogeneous direct numerical simulation (DNS). The peak-to-peak plunging amplitude was fixed at A/c = 0.5 and the Strouhal number was varied from Stc = 0.10 to Stc = 1.00. This parameter regime encompasses flow phenomena of leading-edge vortex (LEV) dominated flow (0.10 ≤ Stc ≤ 0.19), almost vanishing LEV-trailing-edge vortex (TEV) interaction (0.22 ≤ Stc < 0.5), strong previous cycle LEV-TEV interaction (0.49 ≤ Stc ≤ 0.95) and aperiodic flow (Stc ≥ 0.99). For the periodic baseflow, Floquet stability analysis was conducted. Below a Strouhal number of 0.5, the Floquet multiplier is smaller than the static airfoil which indicates the plunging motion stabilises the two-dimensional baseflow. For higher frequencies, a period-doubling mode appears, which has a peak Floquet multiplier around a spanwise wavelength of 2c. This unstable mode also dominates in three-dimensional direct numerical simulations (DNS). Finally, a short-wave mode becomes unstable at Stc = 0.95, which generates more small-scale vorticies in the DNS result.
Lahooti M, Vivarelli G, Montomoli F, et al., 2022, Under-resolved direct numerical simulation of naca0012 at stall
In this work a high-order spectral-h/p element solver is employed to efficiently but accurately resolve the flow field around the NACA0012 aerofoil. In particular, the conditions considered are a Reynolds number of 150000 and three angles of attack, namely 9◦, 10◦ and 12◦. This particular study aims at providing the necessary preliminary insight into the flow dynamics of the turbulent transition at near-post stall with very well resolved Large Eddy Simulation range if not Direct Numerical Simulation. Therefore, the first step consists of determining the mesh resolution required and the spanwise length. Our results repeatedly demonstrate the possible existence of a large-scale low-frequency 2D dominant flow structure over the span where. It was found that employing a single chord length in z is not sufficient to capture it. Concerning the flow behaviour, it is observed that a laminar separation bubble forms at the aerofoil leading edge. This tends to move upstream, shorten in length and increase in height as the angle of attack is increased. Mild three-dimensional behaviour is seen right from the beginning of the aerofoil suction surface with turbulent transition occurring just after the reattachment point. In particular, this is seen to happen sooner with higher aerofoil inclination. Finally, our results indicate an interaction of large-scale structures with the boundary layer.
Reavette RM, Sherwin SJ, Tang M-X, et al., 2021, Wave intensity analysis combined with machine learning can detect impaired stroke volume in simulations of heart failure, Frontiers in Bioengineering and Biotechnology, Vol: 9, Pages: 1-13, ISSN: 2296-4185
Heart failure is treatable, but in the United Kingdom, the 1-, 5- and 10-year mortality rates are 24.1, 54.5 and 75.5%, respectively. The poor prognosis reflects, in part, the lack of specific, simple and affordable diagnostic techniques; the disease is often advanced by the time a diagnosis is made. Previous studies have demonstrated that certain metrics derived from pressure-velocity-based wave intensity analysis are significantly altered in the presence of impaired heart performance when averaged over groups, but to date, no study has examined the diagnostic potential of wave intensity on an individual basis, and, additionally, the pressure waveform can only be obtained accurately using invasive methods, which has inhibited clinical adoption. Here, we investigate whether a new form of wave intensity based on noninvasive measurements of arterial diameter and velocity can detect impaired heart performance in an individual. To do so, we have generated a virtual population of two-thousand elderly subjects, modelling half as healthy controls and half with an impaired stroke volume. All metrics derived from the diameter-velocity-based wave intensity waveforms in the carotid, brachial and radial arteries showed significant crossover between groups-no one metric in any artery could reliably indicate whether a subject's stroke volume was normal or impaired. However, after applying machine learning to the metrics, we found that a support vector classifier could simultaneously achieve up to 99% recall and 95% precision. We conclude that noninvasive wave intensity analysis has significant potential to improve heart failure screening and diagnosis.
Basso R, Hwang Y, Assi G, et al., 2021, Instabilities and sensitivities in a flow over a rotationally flexible cylinder with a rigid splitter plate, Journal of Fluid Mechanics, Vol: 928, Pages: 1-32, ISSN: 0022-1120
This paper investigates the origin of flow-induced instabilities and their sensitivities ina flow over a rotationally flexible circular cylinder with a rigid splitter plate. A linearstability and sensitivity problem is formulated in the Eulerian frame by considering thegeometric nonlinearity arising from the rotational motion of the cylinder which is notpresent in the stationary or purely translating stability methodology. This nonlinearityneeds careful and consistent treatment in the linearised problem particularly whenconsidering the Eulerian frame or reference adopted in this study and not so widelyconsidered. Two types of instabilities arising from the fluid-structure interaction arefound. The first type of the instabilities is the stationary symmetry-breaking mode, whichwas well reported in previous studies. This instability exhibits a strong correlation withthe length of the recirculation zone. A detailed analysis of the instability mode andits sensitivity reveals the importance of the flow near the tip region of the plate for thegeneration and control of this instability mode. The second type is an oscillatory torsionalflapping mode, which has not been well reported. This instability typically emerges whenthe length of the splitter plate is sufficiently long. Unlike the symmetry breaking mode,it is not so closely correlated with the length of the recirculation zone. The sensitivityanalysis however also reveals the crucial role played by the flow near the tip region inthis instability. Finally, it is found that many physical features of this instability arereminiscent of those of the flapping (or flutter instability) observed in a flow over aflexible plate or a flag, suggesting that these instabilities share the same physical origin.
Mengaldo G, Moxey D, Turner M, et al., 2021, Industry-Relevant Implicit Large-Eddy Simulation of a High-Performance Road Car via Spectral/hp Element Methods, SIAM REVIEW, Vol: 63, Pages: 723-755, ISSN: 0036-1445
- Author Web Link
- Cite
- Citations: 7
Moura R, Cassinelli A, da Silva AFC, et al., 2021, Gradient jump penalty stabilisation of spectral/hp element discretisation for under-resolved turbulence simulations, Computer Methods in Applied Mechanics and Engineering, Vol: 388, Pages: 1-29, ISSN: 0045-7825
One of the strengths of the discontinuous Galerkin (DG) method has been its balance between accuracy and robustness, which stems from DG’s intrinsic (upwind) dissipation being biased towards high frequencies/wave numbers. This is particularly useful in high Reynolds-number flow simulations wherelimitations on mesh resolution typically lead to potentially unstable under-resolved scales. In continuous Galerkin (CG) discretisations, similar properties are achievable through the addition of artificial difusion, such as spectral vanishing viscosity (SVV). The latter, although recognised as very useful in CG-based high-fidelity turbulence simulations, has been observed to be sub-optimal when compared toDG at intermediate polynomials orders (P⇡≈3). In this paper we explore an alternative stabilisation approach by the introduction of a continuous interior penalty on the gradient discontinuity at elemental boundaries, which we refer to as a gradient jump penalisation (GJP). Analogous to DG methods, this introduces a penalisation at the elemental interfaces as opposed to the interior element stabilisation of SVV. Detailed eigen analysis of the GJP approach shows its potential as equivalent (sometimes superior) to DG dissipation and hence superior to previous SVV approaches. Through eigenanalysis, a judicious choice of GJP’sP-dependent scaling parameter is made and found to be consistent with previous a-priori error analysis. The favourable properties of the GJP stabilisation approach are also supported by turbulent flow simulations of the incompressible Navier-Stokes equation, as we achieve high-quality flow solutions atP= 3 using GJP, whereas SVV performs marginally worse atP= 5 with twice as many degrees of freedom in total.
Arshad M, Rowland EM, Riemer K, et al., 2021, Improvement and validation of a computational model of flow in the swirling well cell culture model, Biotechnology and Bioengineering, ISSN: 0006-3592
Effects of fluid dynamics on cells are often studied by growing the cells on the base of cylindrical wells or dishes that are swirled on the horizontal platform of an orbital shaker. The swirling culture medium applies a shear stress to the cells that varies in magnitude and directionality from the centre to the edge of the vessel. Computational fluid dynamics methods are used to simulate the flow and hence calculate shear stresses at the base of the well. The shear characteristics at each radial location are then compared with cell behaviour at the same position. Previous simulations have generally ignored effects of surface tension and wetting, and results have only occasionally been experimentally validated. We investigated whether such idealized simulations are sufficiently accurate, examining a commonly-used swirling well configuration. The breaking wave predicted by earlier simulations was not seen, and the edge-to-centre difference in shear magnitude (but not directionality) almost disappeared, when surface tension and wetting were included. Optical measurements of fluid height and velocity agreed well only with the computational model that incorporated surface tension and wetting. These results demonstrate the importance of including accurate fluid properties in computational models of the swirling well method.
Cooke EE, Mughal MS, Sherwin S, et al., 2021, Destabilisation of stationary and travelling crossflow disturbances due to forward and backward facing steps over a swept wing, IUTAM Laminar-Turbulent Transition, Vol: 38, Pages: 713-723
The destabilisation effects of forward and backward facing steps on cross-flow (CF) disturbances on an infinite swept wing is investigated. Stationary and travelling CF-wave instability modulations, as they convect over the abrupt surface features, are investigated computationally with step heights ranging from 18% to 53% of the boundary layer thickness at chordwise locations of 10% and 20%. An embedded mesh approach is used to compute boundary layer base flow profiles over the swept wing with the high order spectral / hp element solver, Nektar++. Linear Stability Theory (LST), Parabolised Stability Equations (PSE) and Linearised Harmonic Navier-Stokes (LHNS) models are used to investigate the development of the convecting CF disturbances. LST is used to understand the instability parameter space and map out neutral curves. PSE equations fail to correctly capture the effects of the steps due to the strong short scale variations introduced whereas, the LHNS provide a rapid and more physics correct technique to ascertain flow destabilisation effects.
This data is extracted from the Web of Science and reproduced under a licence from Thomson Reuters. You may not copy or re-distribute this data in whole or in part without the written consent of the Science business of Thomson Reuters.