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

2D Riemann problem

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 29, 2005, 01:50
Default 2D Riemann problem
  #1
ganesh
Guest
 
Posts: n/a
Dear Friends,

Can anyone suggest the boundary conditions for 2D riemann problems? For a specific configuration wherein 4 shocks appear and interact, I get the shock curving at the boundary, while I have seen that in many simulations by others the shock is normal at the boundary.

Regards,

Ganesh
  Reply With Quote

Old   September 29, 2005, 08:53
Default Re: 2D Riemann problem
  #2
Scott Shaw
Guest
 
Posts: n/a
By boundary do you mean a wall ? If so you should treat it in the same way as you treat sub-sonic inflow or sub-sonic outflow conditions. This is because the component of velocity normal to the wall is zero (no penetration condition). I prefer to treat it as outflow as I then only have to specify one condition (U=0).

Scott
  Reply With Quote

Old   September 29, 2005, 23:14
Default Re: 2D Riemann problem
  #3
Praveen. C
Guest
 
Posts: n/a
The boundary in this case is an artificial boundary meant to truncate the computational domain. There is no rigorous way of deriving such boundary conditions. What most people do is to extrapolate from inside the computational domain. You may have to limit the extrapolated values if there are shocks nearby.
  Reply With Quote

Old   September 30, 2005, 13:53
Default Re: 2D Riemann problem
  #4
ganesh
Guest
 
Posts: n/a
Dear Friends,

I am in fact using extrapolation from the interior with limiting for shocks. The problem does not arise if the initial conditions are such that either one of the velocity components u or v is zero or both are zero in the four quadrants. ( Actually in this case, the normal velocity on the boundary actually vanishes and this makes things easy). However, when one of the quadrants does have a non-zero u and v, then I experience the problem of shock being not normal to the boundary. I have seen people using theextrapolation bc, but they seem to get normality at the boundary.

Any suggestions ?

Regards,

Ganesh

  Reply With Quote

Old   October 1, 2005, 02:38
Default Re: 2D Riemann problem
  #5
Praveen. C
Guest
 
Posts: n/a
You could try a Neumann-type boundary condition,

du/dn = 0

for all the flow quantities. Or you could do a constrained least squares interpolation, ie., fit a polynomial such that the Neumann condition is satisfied either exactly or in a least-squares sense (penalty function approach).
  Reply With Quote

Old   October 3, 2005, 02:33
Default Re: 2D Riemann problem
  #6
ganesh
Guest
 
Posts: n/a
Dear Praveen,

Thankyou for the suggestion on VN boundary condition, it does work. However, what I had to do was to revert to first order accuracy(constant poly. recons.) ie U_boundary = U_(cell centoid of the cell sharing the boundary face). I was wondering if I could do the same with second -order accuacy. Least squares is a possibility, but will it wok effectively ?

Regards,

Ganesh
  Reply With Quote

Old   October 3, 2005, 04:40
Default Re: 2D Riemann problem
  #7
Praveen. C
Guest
 
Posts: n/a
You can try to interpolate to the cell face by satisfying Neumann bc and then compute the flux using the interpolated values. Of course this will be extra work since you have to define a stencil for the cell face. You may have to limit the interpolated values especially if there are shocks. I have used a min-max type of limiting in such situations, ie., if the interpolate exceeds the interval {min, max} then chop it to the end-point of the interval. Its not a nice way to solve but it works.

PS: If you are in IISc, Bangalore, you can talk to Keshav Malagi who has recently solved 2D Riemann problem.
  Reply With Quote

Old   October 3, 2005, 05:55
Default Re: 2D Riemann problem
  #8
ganesh
Guest
 
Posts: n/a
Dear Praveen,

Thanks for the suggestions. I am in IISc and I shall talk to Keshav regarding the 2D Reimann problem.

Regards,

Ganesh

  Reply With Quote

Old   October 3, 2005, 07:58
Default Re: Keshav Malagi
  #9
rmi
Guest
 
Posts: n/a
Praveen,

Am interested Keshav Malagi work. can you introduce his URL or his department......i really need to get in touch with him.

Thanks.

  Reply With Quote

Old   October 3, 2005, 08:09
Default Re: 2D Riemann problem
  #10
A.S.
Guest
 
Posts: n/a
Ganesh,

I came across a report for boundary conditions. You can email me at shukla_apoorv1975@yahoo.co.uk. I have to did it up and will forward it to you. Hope it helps.

Regards

A.S.
  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
Problem Importing Geometry ProE to CFX fatb0y CFX 3 January 14, 2012 20:42
Can I solve this problem by Fluent? Kai_kc FLUENT 1 October 27, 2010 05:29
Velocity profiles problem behind the elbow (3D problem) kabat73 FLUENT 8 May 9, 2010 04:26
Reimann problem Shuo Main CFD Forum 5 May 28, 2007 15:44
Is this problem well posed? Thomas P. Abraham Main CFD Forum 5 September 8, 1999 14:52


All times are GMT -4. The time now is 23:18.