|
[Sponsors] |
June 28, 2013, 07:16 |
WENO accuracy study, a 2D test case?
|
#1 |
New Member
Faraz
Join Date: Jan 2011
Location: Netherlands
Posts: 10
Rep Power: 15 |
Hi,
Is there a standard 2D test case to perform accuracy study of WENO without using periodic boundary conditions? I know that people do tests with periodic BC, but I would like to understand if the boundary effects could decrease convergence rate significantly. Say, for example in a scalar convection case? Is there a paper or a reference to check with that? Thanks |
|
October 16, 2017, 16:45 |
|
#2 |
New Member
Aranya Dan
Join Date: Oct 2017
Posts: 19
Rep Power: 8 |
Try with 2D Reimann problems.
A large list can be found in the paper "Solution of 2d Reimann problems for Gas Dynamics without Reimann Problem Solvers" -A. Kurganov, E. Tadmor |
|
October 17, 2017, 03:01 |
|
#3 |
Senior Member
Join Date: Oct 2011
Posts: 239
Rep Power: 16 |
Hello,
The advantage of periodic bc is that you can simulate a long time/distance advection case and then more easily show the less diffusive properties of high order schemes. Without periodic bc, you would need a much wider grid and computational time would be much higher. If you do not want to work with such a bc, you can still perform accuracy test with advection equation to check the formal order of accuracy of your scheme. But you have to be very careful on the treatment of the limit as indeed it will spoil the observed order of accuracy for finer grids. What I am used to do is to prescribe the exact Riemann state (from the exact solution) at the face gauss nodes and let the Riemann solver do the job. François |
|
October 24, 2017, 15:32 |
More test cases for viscous flows
|
#4 |
New Member
Aranya Dan
Join Date: Oct 2017
Posts: 19
Rep Power: 8 |
Hi,
I have a 2d WENO code for compressible viscous flows and looking for test cases to run it on. I tried running on the lid driven cavity problem but couldnt get it to converge. Can you suggest some other test cases which are simple to set up like the Reimann problems but for viscous flows? |
|
October 24, 2017, 16:10 |
|
#5 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71 |
Quote:
Start first with a classical test on the Euler equations, you can also have a look to the book of LeVeque. |
||
October 24, 2017, 16:28 |
|
#6 | |
New Member
Aranya Dan
Join Date: Oct 2017
Posts: 19
Rep Power: 8 |
Quote:
I have ran tests with inviscid flow on 2d reimann problems and have got acceptable results. I just wanted to check my results after adding the viscous terms, thereby making it a full Navier Stokes solution. So yes, i need some simple test cases for the NS equations. 0012 is a widely used test case which i dont want to attempt as it would take a lot of time to just set up the grid and transformations. I will have a look at the book you mentioned but please do tell me if you know of any easy test cases for NS solver. Thanks. |
||
October 24, 2017, 16:39 |
|
#7 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71 |
Quote:
That strongly depends if you are interested in testing the capability of your code in solving flows with relevant compressibility effects or not. If your code can solve low Mach flows, you can use standard test-cases for incompressible flows (Taylor solution, backward facing step flow, etc). |
||
October 24, 2017, 16:50 |
|
#8 | |
New Member
Aranya Dan
Join Date: Oct 2017
Posts: 19
Rep Power: 8 |
Quote:
Also, I am using rho_inf*u_inf^2 to non dimensionalise the pressure, so Mach number doesnt explicitly show up in the equations I am solving. Would that cause any problems or should I go for the conventional 1/(gamma*M_inf^2) which is typically followed in high-speed flow literature? My final aim is to study shock wave boundary layer interactions so I would be running the final code at high Mach numbers. |
||
October 24, 2017, 17:07 |
|
#9 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71 |
Mach number appears intrinsically as a response between the density and pressure variation. For a low Mach number flow, some codes use pre-conditioning techniques. But I see some codes working without that.
Try first a steady laminar flow in a channel, for example as described here http://acoustics.ae.illinois.edu/pdf...-lele-1992.pdf |
|
October 24, 2017, 17:15 |
|
#10 |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71 |
try also searching here https://www.cfd-online.com/Wiki/Vali...2-D_test_cases
|
|
October 25, 2017, 17:12 |
|
#11 | |
New Member
Aranya Dan
Join Date: Oct 2017
Posts: 19
Rep Power: 8 |
Quote:
How good would be a Taylor Green Vortex (2d) as a test case? With periodic boundary conditions i am getting a drift of the vortex and it doesnt match at all with the theoretical (exponential decay) solution. So i am confused as some test case (2d Reimann inviscid gives good results while others (Taylor vortex) isnt even near to what i expect. |
||
October 25, 2017, 17:20 |
|
#12 | |
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,768
Rep Power: 71 |
Quote:
|
||
October 26, 2017, 10:45 |
Weno
|
#13 |
Senior Member
Selig
Join Date: Jul 2016
Posts: 213
Rep Power: 10 |
I don't mean to be off topic, but when using WENO for incompressible NS. If I want to discretize the nonlinear advection term, do I need to use a hyperbolic scheme on top of WENO or can I just perform WENO on the nonlinear terms? In compressible cases you see Local Lax Friedrichs + WENO (for example.)
|
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Pressure instabilities with interDyMFoam for the floatingObject case | nbadano | OpenFOAM Running, Solving & CFD | 15 | October 15, 2021 06:35 |
Help Required With Simple Test Case | steph79 | OpenFOAM Pre-Processing | 4 | August 3, 2010 07:45 |
Exact unsteady Stokes test case | neazen | Main CFD Forum | 5 | April 13, 2009 14:15 |
Porous Media test case | Alex | FLUENT | 0 | April 9, 2006 08:23 |
c1 body test case | Eric Lenormand | Main CFD Forum | 0 | March 2, 2000 06:54 |