CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

Solution of SOD Shock Tube Problem by using HARTEN Modified Flux TVD Scheme

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   December 13, 2012, 01:40
Default Solution of SOD Shock Tube Problem by using HARTEN Modified Flux TVD Scheme
  #1
New Member
 
Join Date: Nov 2012
Posts: 5
Rep Power: 4
xue sheng is on a distinguished road
I am solving 1D shock tube SOD problem by using Harten modified flux TVD scheme. My code is stable till CFL less than 0.6, after that it is unstable.
It should be stable till CFL less than 1.

-1=<x<0 P=1, RHO=1, V=0
0=<x<1 P=0.1, RHO=0.125, V=0

I have following question,

1. How should calculate Time Step Size (DT)? Is this formula correct?
DT = (CFL x DX)/Cmax

2. How should calculate DX? Is this formula correct?
DX=L/(200) where L = 2 (201=Grid points)

3. How should calculate Cmax? Is this formula correct?
Cmax=0
Do i=1 to 200
C=DSQRT(GAMA*Pressure(i)/Density(i))
Cmax=DMAX1(C,Cmax)
Continue

Anyone have solved this problem by using Harten modified flux TVD scheme?
xue sheng is offline   Reply With Quote

Old   December 13, 2012, 04:20
Default
  #2
Super Moderator
 
praveen's Avatar
 
Praveen. C
Join Date: Mar 2009
Location: Bangalore
Posts: 244
Blog Entries: 6
Rep Power: 9
praveen is on a distinguished road
You must also include fluid velocity. It might increase the wave speed

C = ABS(V(i)) + DSQRT(GAMA*Pressure(i)/Density(i))
praveen is offline   Reply With Quote

Old   December 13, 2012, 05:23
Default
  #3
New Member
 
Join Date: Nov 2012
Posts: 5
Rep Power: 4
xue sheng is on a distinguished road
Thanks for correction.
But it would further decrease CFL criterion.
HARTEN TVD scheme should be stable below CFL <1.
But in my case it is stable around CFL<0.5.

May you please point out the reason?
xue sheng is offline   Reply With Quote

Old   December 13, 2012, 05:31
Default
  #4
Super Moderator
 
praveen's Avatar
 
Praveen. C
Join Date: Mar 2009
Location: Bangalore
Posts: 244
Blog Entries: 6
Rep Power: 9
praveen is on a distinguished road
As I wrote your CFL condition is wrong. You must use largest wave speed to compute time step which is what I wrote in my last post.
praveen is offline   Reply With Quote

Old   December 13, 2012, 06:11
Default
  #5
New Member
 
Join Date: Nov 2012
Posts: 5
Rep Power: 4
xue sheng is on a distinguished road
Thanks.
After your recommended correction Its working well.
Now my Harten modified flux TVD scheme, ROE scheme, and LAX_WENDROFF scheme Code working well till CFL <1.

But my HARTEN-YEE TVD scheme code is stable till CFL<0.6. After that it become unstable.

May you please inform me that HARTEN modified flux TVD scheme and HARTEN-YEE TVD scheme has same CFL restriction or different?

I think both have same CFL restriction. But in my code HARTEN-YEE become unstable at CFL>0.6.
May you please point out the reason?
xue sheng is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
analytical solution of 1D shock tube problem ma Main CFD Forum 10 May 11, 2015 05:20
Sod shock tube problem with real gas effects ganesh Main CFD Forum 0 October 17, 2006 01:22
SOD shock tube problem Luuk Main CFD Forum 3 September 6, 2005 05:08
help: analysis solution of gas-water sod problem gooo Main CFD Forum 0 April 19, 2005 01:05
Shock tube problem ger Main CFD Forum 1 September 24, 2003 07:41


All times are GMT -4. The time now is 22:36.