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

smkhan · July 10, 2010, 05:29

Hi!
i am using own code. actually i m working on 1_D Riemann Problems, 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 and t=0.012.

DoHander · July 10, 2010, 09:24

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 10, 2010, 05:29	shock detection in the Riemann problems(specially in the stationary contact problems)	#1
smkhan New Member S M Khan Join Date: Sep 2009 Posts: 17 Rep Power: 17	Hi! i am using own code. actually i m working on 1_D Riemann Problems, 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 and t=0.012.

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
shock detection in the Riemann problems(specially in the stationary contact problems)	smkhan	Main CFD Forum	4	July 10, 2010 09:25
problems with UDF to set contact angle	poiuy219	Main CFD Forum	0	April 30, 2009 11:43
How to set stationary shock?	Paul	Main CFD Forum	1	October 14, 2001 11:31

July 10, 2010, 09:24		#2
DoHander Senior Member Join Date: Nov 2009 Posts: 411 Rep Power: 19	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