|March 29, 2016, 05:29||
Taylor-Green vortex sheet
Join Date: Jan 2011
Location: ETSIN, Madrid
Posts: 16Rep Power: 7
I have been simulating the Taylor-Green vortex sheet with OpenFOAM. My main target is to establish how the error decreases with the time-step Dt, for fixed Courant number. I am guessing it should be proporcional to Dt^2, as shown in this Example.
However, I am getting this puzzling result
It seems to show a very rapid decay, as Dt^4 for longer times, but then the error stalls completely. I thought it may have to do with the precision used to write the results (I have had those issues in the past), but increasing the number of decimals does not change much.
I would really appreciate if anyone has an explanation for this. It could be that the method does not improve below some spacial resolution? I am really interested in this because I am implementing a new CFD particle method, and it would be great to compare it against OpenFOAM.
I am uploading a run directory for the second smallest time-step. The rest can be easily re-created by doubling the time-step in controlDict and the cell size in blockMeshDict (or the opposite). Notice I am running at Re=200, and Co=0.5, at variance with the example quoted above. The python utility I have written to obtain the error is also included. The final error lacks a sqrt, which is included in the figure above.
|Thread||Thread Starter||Forum||Replies||Last Post|
|Presenation of vortex strength and velocity concept?||fruitkiwi||Main CFD Forum||0||September 26, 2012 22:08|
|vortex cause pressure gradient or pressure gradient induce vortex?||fruitkiwi||Main CFD Forum||4||June 12, 2012 01:12|
|creating vortex core line||aay023||Main CFD Forum||0||September 15, 2010 21:49|
|turbulent scales of forced vortex and free vortex||lcw||Main CFD Forum||3||September 1, 2005 13:40|
|help in galvanising sheet problem||MANOJ||FLUENT||8||August 22, 2005 02:29|