|
[Sponsors] |
January 27, 2006, 10:27 |
Jameson Dissipation.
|
#1 |
Guest
Posts: n/a
|
I have a quasi-one-dimension nozzle flow problem and I need to solve the equation using Jameson's scheme coding a cell vertex version as well as suitable second and fourth order disspipation. I can't seem to find the right equations for the dissipation anywhere and could do with some help in finding them.
Any suggestions? Any more info needed just let me know. Pete |
|
January 27, 2006, 12:20 |
Re: Jameson Dissipation.
|
#2 |
Guest
Posts: n/a
|
Try looking through the extensive list of pdf files (mostly near the bottom) at
http://aero-comlab.stanford.edu/jame...tion_list.html |
|
January 28, 2006, 08:22 |
Re: Jameson Dissipation.
|
#3 |
Guest
Posts: n/a
|
Ok I have had a look through those and I still can't seem to get it right. Maybe if I post up the problem and the code I have got so far, maybe someone could let me know how to code what I need. The area I need coded is currently filled with stars:
<https://www.bris.ac.uk/fluff/u/pc583...lk4ahwCr9fgi7/> - Problem (pdf) <https://www.bris.ac.uk/fluff/u/pc583...ClluM6LIigwi7/> - Code (cpp) <https://www.bris.ac.uk/fluff/u/pc583..._OHZJoDJZDQi7/> - Include / Library File <https://www.bris.ac.uk/fluff/u/pc583...e4UgDc9etOwi7/> - Include / Library File <https://www.bris.ac.uk/fluff/u/pc583...RuIV0agK4Igi7/> - Include / Library File <https://www.bris.ac.uk/fluff/u/pc583...0h_WzCDWaXgi7/> - Include / Library File |
|
January 28, 2006, 20:06 |
Re: Jameson Dissipation.
|
#4 |
Guest
Posts: n/a
|
The reference you are looking for is number 62 in the list for which Jeff provided the link. That paper describes the original JST scheme in all the detail that you will need for implementation. Read it carefully and all your questions should be answered.
|
|
January 29, 2006, 14:31 |
Re: Jameson Dissipation.
|
#5 |
Guest
Posts: n/a
|
I enter this thread with a question. The paper 62 is ok. But I found nowhere (after extensive searching) the contruction of the di-1/2,j therm. Is there (rho i,j -rho i-1,j) or (rho i-1,j - rho i,j) ? And the second therm in eq.(11)? I tray now to implement the Jameson's diffusion in a 3D transonic compressor code, but it does not work. I reduced the cfl from 2.8 down to 0.5. There are this constants 1/4, 1/256 free for variation. I do something wrong but I do not know why. Is somebody experienced with writing such a code?
|
|
February 1, 2006, 04:52 |
Re: Jameson Dissipation.
|
#6 |
Guest
Posts: n/a
|
That's exactly what I am looking for. I need someone to treat me like the idiot that I am and show me step by step how to implement Jameson's equations into code as everything I am doing at the moment doesn't seem to be working.
|
|
February 1, 2006, 07:31 |
Re: Jameson Dissipation.
|
#7 |
Guest
Posts: n/a
|
Try using 1/2 and 1/64 they are more stable as per my experience
|
|
February 1, 2006, 13:52 |
Re: Jameson Dissipation.
|
#8 |
Guest
Posts: n/a
|
Thanks,
I apply the diffusion term something like this [simplyfied in 1D, u1 is density, x momentum, etc., in general u1(i,j,k)] d(i)= 1scalep*(eps2p*(u1(i+1)-u1(i))- 1 eps4p*(u1(i+2)-3.*u1(i+1)+3.*u1(i)-u1(i-1)))- 1scalem*(eps2m*(u1(i)-u1(i-1))- 1 eps4m*(u1(i-2)-3.*u1(i-1)+3.*u1(i)-u1(i+1))) This gives in the result something similar like 2nd derivative minus 3rd derivative (?????). The scalep and scalem are scaling factors. I tried to set the ratio of (volume(i+1)+volume(i))/2./timestep for scalep, etc. but I do not know if it is ok. eps2p,etc. are Jamesons factors The d(i) values I calculate in all 3 dimesions and add it to one number. For each control volume I sum the fluxes in x, y, z directions and substract d. That is all. But will not work. |
|
February 3, 2006, 10:31 |
Re: Jameson Dissipation.
|
#9 |
Guest
Posts: n/a
|
Other interpretation of Jamesons' AIAA paper is
d(i)= 1 scalep*(eps2p*(u1(i+1)-u1(i))- 1 ........eps4p*(u1(i+2)-3.*u1(i+1)+3.*u1(i)-u1(i-1)))- 1 scalem*(eps2m*(u1(i)-u1(i-1))+ plus here!!!!!! 1 ........eps4m*(u1(i-2)-3.*u1(i-1)+3.*u1(i)-u1(i+1))) This gives in the result something similar like 2nd derivative minus 4th derivative (?????). And it seems to work with 1/2 and 1/64 excelent for Q1D. But in 3D case there must be some trick else. It works but not satisfactory. Who knows how are treated the cells at wall, where the points i+1, i+2 are not present. The next trick maz be the scaling scalep,...etc. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
CFX12 rif errors | romance | CFX | 4 | October 26, 2009 13:41 |
Dissipation in Harten-Yee and Davis-Yee TVD Scheme | cfdxue | Main CFD Forum | 0 | April 29, 2008 10:45 |
Dissipation rate access | K. Kevala | FLUENT | 0 | July 14, 2004 11:11 |
artifical dissipation on unstructured meshes | Frederic | Main CFD Forum | 0 | October 31, 2001 17:42 |
calculation of (turbulent) dissipation energy ? | max | Main CFD Forum | 3 | August 18, 1999 09:42 |