For industrial LES, in addition to accuracy and efficiency, there are several other important factors to consider:
- Ability to handle complex geometries, and ease of mesh generation
- Robustness for a wide variety of flow problems
- Scalability on supercomputers
Our own research effort has led to the development of a high-order solver based on the FR/CPR method called hpMusic. We recently performed a benchmark LES comparison between hpMusic and a leading commercial solver, on the same family of hybrid meshes at a transonic condition with a Reynolds number more than 1M. The 3rd order hpMusic simulation has 9.6M degrees of freedom (DOFs), and costs about 1/3 the CPU time of the 2nd order simulation, which has 28.7M DOFs, using the commercial solver. Furthermore, the 3rd order simulation is much more accurate as shown in Figure 1. It is estimated that hpMusic would be an order magnitude faster to achieve a similar accuracy. This study will be presented at AIAA's SciTech 2018 conference next week.
(a) hpMusic 3rd Order, 9.6M DOFs
(b) Commercial Solver, 2nd Order, 28.7M DOFs
Figure 1. Comparison of Q-criterion and Schlieren
I certainly believe high-order solvers are ready for industrial LES. In fact, the commercial version of our high-order solver, hoMusic (pronounced hi-o-music), is announced by hoCFD LLC (disclaimer: I am the company founder). Give it a try for your problems, and you may be surprised. Academic and trial uses are completely free. Just visit hocfd.com to download the solver. A GUI has been developed to simplify problem setup. Your thoughts and comments are highly welcome.
Happy 2018!
Strange that in an energy preserving 2nd order method, with a much more refined grid, the compressive shocks are absent...
ReplyDeleteIt is very difficult to know what is being used exactly with a commercial code. But I believe the SGS model and the numerical method play a critical part on the quality of numerical results.
ReplyDelete2nd order upwind is too dissipative for LES, the code should have used a hybrid (upwind/central) scheme. I find also strange that those compression waves are not captured, pointing again to probable setup issue. This problem is neither a complex geometry nor very challenging And I expect both codes to a very reasonable job for this LES simulation . It is a very used benchmark problem, which can be run very efficiently even with legacy FV codes.
ReplyDeleteWhat I meant by “.. not very challenging” is that it does not involve a very complex procedure to “setup” the correct unsteady upstream turbulent field, but it is as challenging as one may expect in terms of domain size, mesh and solver setup for a “proper” LES...
ReplyDeleteAdding incoming turbulence will add uncertainty in this problem. We may consider that in the future.
DeleteDid you make a comparison with OpenFOAM? It has high order convective schemes.
ReplyDeleteWe did not. But if somebody is interested, please let me know.
DeleteThis comment has been removed by the author.
ReplyDeleteA hybrid upwind/central scheme was indeed used for the commercial silver. We did see some spurious oscillations, but not able to remove them. Would be interested in seeing other codes’ performance for this case. If you can suggest another case, we would be happy to take a look.
ReplyDeleteI know this is an oldish thread now but I'm really interested in the status of higher order methods. Let's keep the discussion going!
ReplyDelete