|
[Sponsors] |
August 17, 2015, 16:24 |
Experience or Caveats using CE/SE Method?
|
#1 |
Senior Member
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 25 |
Most of my experience has been with pressure-based CFD methods. I have recently been investigating density-based...especially explicit methods in earnest, and I bumped into the CE/SE method recently, and it seems impressive--dissipation free even on unstructured meshes, time-accurate local timestepping, no need for directional splitting or Reimannian solvers, etc. It makes me question why it is not more popular.
I've found other short CE/SE discussions on this board in the past, but it seems like it is more of a curiosity to those who have heard of it. Does anyone here have direct experience with the method? Are there any gotchas here that I am missing? Or it is that RK-TVD-Reimann solvers are just too entrenched and this can't gain traction? Thanks |
|
August 19, 2015, 01:02 |
|
#2 |
Super Moderator
|
I have also wondered why there is not much discussion of this. Whenever I read papers on high order methods, people only discuss DG, spectral volume, spectral difference, flux reconstruction, weno, etc. Almost nobody mentions CE/SE method. I did not see any specific advantages of CE/SE method that would make it better than these other methods other than not requiring riemann solvers. When I look at the literature on high order methods in aerospace CFD, it seems to me that DG schemes are the dominant ones (for high Re problems). See e.g., the HIOCFD workshops.
Here is a talk on CE/SE http://nia-mediasite.nianet.org/NIAM...2-0d8a1a25ea58 http://ossanworld.com/hiroakinishika...Chang_CESE.pdf |
|
August 19, 2015, 01:21 |
|
#3 | |
Senior Member
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 25 |
Quote:
I see advantages over DG, mainly in the much less restrictive CFL--for DG, the CFL limit gets increasingly worse as order increases. WENO is an implementation and computational nightmare on unstructured meshes. Even TVD limiting on unstructured meshes...if done in a vector context per Berger's paper ...is a very difficult implementation/computational task (solving a nontrivial linear programming problem per cell per field). For vector implementations (AVX2/AVX512/OpenCL), avoiding the limiting and Reimannian solvers could be a big performance win just because the code is not loaded with conditionals and iterative stages. I am specifically looking at time-accurate methods and the prospect of local timestepping is very attractive. We likely won't have to deal with boundary layer type cell-size variations, but I would like to avoid having my global timestep limited by the smallest cell. I honestly wonder if "element" in the name is throwing people off. This is a Finite Volume-type method, even moreso than DG. |
||
August 19, 2015, 01:32 |
|
#4 |
Super Moderator
|
There are explicit space-time DG schemes available, where you can do time accurate local time stepping, see papers of Gassner, Munz, Dumbser, etc. You can go to very high orders of accuracy.
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
addCmptAvBoundaryDiag method | santiagomarquezd | OpenFOAM Programming & Development | 3 | December 26, 2020 08:29 |
Lattice Boltzmann method vs Finite Element Method and Finite Volume Method | solemnpriest | Main CFD Forum | 3 | August 12, 2013 11:00 |
Implicit iterative Gauss-Seidel method | dearboy | Main CFD Forum | 5 | November 29, 2010 09:18 |
Multigrid method | Suresh kumar | Main CFD Forum | 1 | April 24, 2006 20:07 |
Has the pressure correction method problems? | X. Ye | Main CFD Forum | 17 | September 1, 1999 16:14 |