shock detection in the Riemann problems(specially in the stationary contact problems)

smkhan · July 9, 2010, 12:22

Hello!
i am not able to detect the shock correctly in the Riemann problems and getting large oscillations due to this reason.
thanks

smkhan · July 9, 2010, 12:25

Hello!
i am not able to detect the shock correctly in the Riemann problems and getting large oscillations due to this reason. pliz let me know what conditions should be put to eliminate the oscillations.
thanks

DoHander · July 9, 2010, 16:27

Are you using a commercial code or your own code ?
Also can you describe the initial settings for your problem:
1d/2d
T,p,v ... at t=0.

Do

smkhan · July 10, 2010, 04:28

Hi DoHander!
i am using own code. actually i m working on Riemann Problem, conditions are
domain length=1.0, interface x=0.8,RHO_L=1.0,U_L=-19.59745,P_L=1000and
RHO_R=1.0,U_R=-19.59745,P_R=0.01. L stands for left and R stands for right.

DoHander · July 10, 2010, 09:25

First you should have an analytical solution for your problem if you don't have one (you can compare this with the numerical one):

J.D. Anderson - Modern compressible flow.

Second - if your numerical solutions oscillates you probably have an implementation problem: maybe with your boundary conditions or maybe with the order of the numerical scheme you use. For a second order scheme (and high orders) you must use limiters functions in order to remove the oscillations (or you can try with a very fine mesh).

Do

July 9, 2010, 12:22	shock detection in the Riemann problems	#1
smkhan New Member S M Khan Join Date: Sep 2009 Posts: 17 Rep Power: 17	Hello! i am not able to detect the shock correctly in the Riemann problems and getting large oscillations due to this reason. thanks

July 9, 2010, 12:25	shock detection in the Riemann problems(specially in the stationary contact problems)	#2
smkhan New Member S M Khan Join Date: Sep 2009 Posts: 17 Rep Power: 17	Hello! i am not able to detect the shock correctly in the Riemann problems and getting large oscillations due to this reason. pliz let me know what conditions should be put to eliminate the oscillations. thanks

July 10, 2010, 04:28	shock detection	#4
smkhan New Member S M Khan Join Date: Sep 2009 Posts: 17 Rep Power: 17	Hi DoHander! i am using own code. actually i m working on Riemann Problem, conditions are domain length=1.0, interface x=0.8,RHO_L=1.0,U_L=-19.59745,P_L=1000and RHO_R=1.0,U_R=-19.59745,P_R=0.01. L stands for left and R stands for right.

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
HLL Riemann Shock Tube Matlab Problem	Luke F	Main CFD Forum	2	May 20, 2016 03:10
Slug Flow, interFoam, problems with Contact Angle	PrzemekPL	OpenFOAM Running, Solving & CFD	13	February 18, 2014 23:10
contact problems for high accuracy Euler algorithm	igor	Main CFD Forum	0	May 28, 2003 08:24
How to set stationary shock?	Paul	Main CFD Forum	1	October 14, 2001 11:31

July 9, 2010, 16:27		#3
DoHander Senior Member Join Date: Nov 2009 Posts: 411 Rep Power: 20	Are you using a commercial code or your own code ? Also can you describe the initial settings for your problem: 1d/2d T,p,v ... at t=0. Do

July 10, 2010, 09:25		#5
DoHander Senior Member Join Date: Nov 2009 Posts: 411 Rep Power: 20	First you should have an analytical solution for your problem if you don't have one (you can compare this with the numerical one): J.D. Anderson - Modern compressible flow. Second - if your numerical solutions oscillates you probably have an implementation problem: maybe with your boundary conditions or maybe with the order of the numerical scheme you use. For a second order scheme (and high orders) you must use limiters functions in order to remove the oscillations (or you can try with a very fine mesh). Do