📣 Help Shape the Future of UKRI's Gateway to Research (GtR)

We're improving UKRI's Gateway to Research and are seeking your input! If you would be interested in being interviewed about the improvements we're making and to have your say about how we can make GtR more user-friendly, impactful, and effective for the Research and Innovation community, please email gateway@ukri.org.

Adaptively Tuned High-Order Unstructured Finite-Volume Methods for Turbulent Flows

Lead Research Organisation: CRANFIELD UNIVERSITY
Department Name: Sch of Aerospace, Transport & Manufact

Abstract

Turbulent fluid flow is characterised by the seemingly disordered and unpredictable movement of fluid particles that begin to swirl and create eddies. Several processes and products such as wind-turbines, combustion, cars, aircrafts are fated to coexist with turbulence. Therefore, understanding turbulent flows, is of paramount importance to improve the performance of many processes and products. Until the full secrets of the Navier-Stokes equations that describe fluid flows are unlocked, researchers across several disciplines will continue to rely on experiments and computational fluid dynamics (CFD) for casting new light on turbulent flows.

Most turbulent flows of practical importance to scientists and engineers are characterised by high-Reynolds numbers, with a wide range of spatial and temporal scales available. The resolution required to resolve all the scales present using direct numerical simulation (DNS) of the Navier-Stokes equations is not feasible due to the computational cost, even with the most powerful supercomputers available today. Large eddy simulation (LES) family of methods have established themselves as an elegant alternative for transient simulations, where the large scales are resolved, and the effects of the unresolved small scales is modelled.

Non-linear numerical methods are widely used for compressible LES simulations since they can exhibit two important properties. Firstly, a non-oscillatory behaviour across flow discontinuities and secondly a low-numerical dissipation in smooth regions of the flow. A considerable body of research work has been devoted in striking the perfect balance between these two properties. However, their tuning has been limited to controlling the numerical dissipation and dispersion on uniform Cartesian meshes, since there are no established techniques to expand their tuning to meshes consisting of different element types.

This research seeks to develop a new class of non-linear methods, that can offer automatic dissipation and dispersion adjustment (ADDA), through the introduction of suitable novel metrics, that are applicable to any type of unstructured mesh. These methods will be applied to a series of carefully selected turbulent flow problems including moderate and high Reynolds numbers using unstructured meshes, and using the national ARCHER2 High Performance Computing (HPC) facility. The developed methods will be available in the free open source ucns3d CFD software and will be supported by the generated datasets that will be available in a open repositories.


Finally this research will make a considerable impact on the UK fluid dynamics industry that generates £13.9 billion worth of output from over 2,200 firms and employs 45,000 people, while the total UK turnover of firms engaged in fluid dynamics exceeds £200 billion and employ over half a million people according to the latest UK Fluids Network Sept 2021 report. This research will contribute towards UK-led research that enables the development of CFD software that is routinely employed in all industry sectors and has established the UK as a world leader in CFD that in turn drives global companies to engage with UK experts for driving innovation for improving the quality of life.
 
Description We have introduced a new metric being based on kinetic Energy dissipation rate to understand fundamental turbulent flow modelling properties of high-order method. This metric is then used to drive the adaptation of numerical methods for the correct modelling of under-resolved turbulent flows. Initial results demonstrate that we can consistently obtain the correct energy transfer from large scales to small scales, supported by the Kolmogorov theory.
Exploitation Route I envisage that the outcome of this research will lead to new exciting directions of research using the developed methods and in particular the following applications:
1) High-fidelity turbulence modelling of compressible flows (supersonic, hypersonic) over complex geometries
2) Use of the metrics to train AI techniques to detect state of turbulence (resolved, unresolved, artificial due to numerical method dissipation etc) from off-line and online datasets
3) Further advance the techniques for multi physics modelling of compressible flows
Sectors Aerospace

Defence and Marine

Energy

 
Description Our findings have already being used by the scientific community using computational fluid dynamics. All the developed methods and tools so far have been implemented and made available in the established open-source high-order CFD code (UCNS3D). The methods and tools developed now form a valuable research and teaching resource for the scientific computing community, and higher fidelity of simulations can be performed at reduced computational cost at large HPC facilities that imply a reduced economic footprint.
First Year Of Impact 2022
Sector Aerospace, Defence and Marine,Digital/Communication/Information Technologies (including Software),Education
Impact Types Economic

 
Title Data supporting: "High-order hybrid DG-FV framework for compressible multi-fluid problems on unstructured meshes" 
Description This dataset contains binary output in Tecplot format for the test problems analysed in the "High-order hybrid DG-FV framework for compressible multi-fluid problems on unstructured meshes" JCP paper. Test cases included are: - Gas-water isolated material interface advection - 2D and 3D helium bubble interaction with shock wave - 2D shock driven air bubble collapse in water - 2D and 3D shock driven air bubble array collapse in water - 2D underwater explosion 
Type Of Material Database/Collection of data 
Year Produced 2024 
Provided To Others? Yes  
URL https://cord.cranfield.ac.uk/articles/dataset/Data_supporting_High-order_hybrid_DG-FV_framework_for_...
 
Title Data supporting: "High-order hybrid DG-FV framework for compressible multi-fluid problems on unstructured meshes" 
Description This dataset contains binary output in Tecplot format for the test problems analysed in the "High-order hybrid DG-FV framework for compressible multi-fluid problems on unstructured meshes" JCP paper. Test cases included are: - Gas-water isolated material interface advection - 2D and 3D helium bubble interaction with shock wave - 2D shock driven air bubble collapse in water - 2D and 3D shock driven air bubble array collapse in water - 2D underwater explosion 
Type Of Material Database/Collection of data 
Year Produced 2024 
Provided To Others? Yes  
Impact This dataset supports a recently accepted journal publication where a series of algorithms are assessed for switching between DG-FV for well established mutliphysics problems. This algorithm performed in a-priori fashion is one of the most crucial building blocks for the adaptive tuned methods used for turbulent flows. Part of this research involves testing these algorithms for compressible flow simulations. 
URL https://cord.cranfield.ac.uk/articles/dataset/Data_supporting_High-order_hybrid_DG-FV_framework_for_...
 
Title Data supporting: 'A short note on a 3D spectral analysis for turbulent flows on unstructured meshes' 
Description This compressed file includes the following folders: POST_X: directory with the FFT from all the simulations postprocessed with the fortran interpolator SPECTRA FORTRAN INTERPOLATOR: A simple fortran program to interpolate data sets from unstructured meshes to the auxiliarry structured mesh (ASM) TECPLOT_LAYOUTS_DATA: All the tecplot layouts and datasets used for the images in the article. 
Type Of Material Database/Collection of data 
Year Produced 2022 
Provided To Others? Yes  
Impact We provide a new technique for obtaining the 3D FFT for turbulent flows, something that can impact other disciplines where ustructured datasets are used to obtained 3D FFT. 
URL https://cord.cranfield.ac.uk/articles/dataset/Data_supporting_A_short_note_on_a_3D_spectral_analysis...
 
Title Data supporting: 'Hybrid discontinuous Galerkin-finite volume techniques for compressible flows on unstructured meshes' 
Description This dataset includes the following test problems binary outputs in tecplot format: Sonic Boom Schardin TGV-Subsonic TGV-Supersonic RUNS_2D 
Type Of Material Database/Collection of data 
Year Produced 2022 
Provided To Others? Yes  
Impact We provide the datasets for a series of well-established test problems where the substantially improved accuracy of the hybrid DG-FV methodology introduced is demonstrated. All the methods have been implemented and deployed within the UCNS3D open-source high-order unstructured Computational Fluid Dynamics (CFD) solver. The present coupling has the potential to improve the shortcomings of both FV-DG in a computational efficient manner. The improved accuracy and robustness provided is a characteristic of paramount importance for industrial-scale CFD applications, and favours the extension to other systems of governing equations. 
URL https://cord.cranfield.ac.uk/articles/dataset/Data_supporting_Hybrid_discontinuous_Galerkin-finite_v...
 
Title Dataset for Arbitrary-Order Unstructured Finite-Volume Methods for Implicit Large Eddy Simulation of Turbulent Flows with Adaptive Dissipation/Dispersion Adjustment (ADDA) 
Description This dataset contains the folders for each of the simulations of this article including: 4.1 (1_CALIBRATION_ER) 5.1 (2 TGV_SUBSONIC) 5.2 (3_SD7003) 5.3 (4_TGV_SUPERSONIC) 5.4 (5_CYLINDER) Each folder contains the following: 1) meshes used to run the simulations 2) UCNS3D.DAT parameter files for UCNS3D CFD SOFTWARE (freely available at https://github.com/ucns3d-team/UCNS3D) 3) output files in tecplot format/or vtu (that can be used in tecplot/paraview/visit post-processing software packages) 4) restart/checkpoint files 
Type Of Material Database/Collection of data 
Year Produced 2024 
Provided To Others? Yes  
Impact This dataset encapsulates the ADDA technique developed in this work and is associated with its calibration on the Taylor-Green vortex flow problem 4 other flow problems. 
URL https://dspace.lib.cranfield.ac.uk/handle/1826/22991
 
Title Numerical simulation of the inviscid Taylor-Green Vortex using a Discontinuous Galerkin method with adaptive filtering to stabilize the solution: data 
Description Numerical simulation of the inviscid Taylor-Green Vortex using a Discontinuous Galerkin method with adaptive filtering to stabilize the solution. 
Type Of Material Database/Collection of data 
Year Produced 2024 
Provided To Others? Yes  
Impact This dataset supports a recently accepted journal publication where the development of an adaptive dissipation algorithm (ADDA) is detailed, and deployed for the iLES of the Taylor-Green vortex flow problem, and several statistics are obtained. The results from the simulations with the ADDA algorithm obtained with the open-source freely available UCNS3D CFD code are provided in this dataset. 
URL https://cord.cranfield.ac.uk/articles/dataset/Numerical_simulation_of_the_inviscid_Taylor-Green_Vort...
 
Title Numerical simulation of the inviscid Taylor-Green Vortex using a Discontinuous Galerkin method with adaptive filtering to stabilize the solution: data 
Description Numerical simulation of the inviscid Taylor-Green Vortex using a Discontinuous Galerkin method with adaptive filtering to stabilize the solution. 
Type Of Material Database/Collection of data 
Year Produced 2024 
Provided To Others? Yes  
URL https://cord.cranfield.ac.uk/articles/dataset/Numerical_simulation_of_the_inviscid_Taylor-Green_Vort...