# Riemann Solvers in multidimension

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

 March 26, 2007, 12:08 Riemann Solvers in multidimension #1 DSS Guest   Posts: n/a Hi, Given two state vectors at left and right side of the face, can anyone explain me how to proceed with the Riemann Solvers? In case of cartesian mesh, I guess it is straight forward to treat it dimension by dimension. But, when the face normal is not aligned with the axis, should we do a coordinate transform for the state vectors? How do we calculate the left and right wave speeds? I am a bit confused becasue, the HLLC solver works for cartesian coordinates, but, fails with unstructed grid. If I use a simple averaging, it works fine. So, I am sure I have some error in implementation of the Riemann solver. Thanks in advance for your help.

 March 26, 2007, 13:04 Re: Riemann Solvers in multidimension #2 ag Guest   Posts: n/a Any Riemann solver can be used for unstructured grids. The Riemann problem is considered as taking place in the direction normal to the cell face. Use the cell quantities (with higher-order interpolation if desired) to construct the left and right states. The flux passing through the face is then determined. Just construct a local coordinate system at the cell face.

 March 26, 2007, 23:19 Re: Riemann Solvers in multidimension #3 DSS Guest   Posts: n/a Hi ag, Thanks for your response. I figured out some thing... In the HLL flux calculation, when I compute the fluxes (left and right), I calculated them in X and Y directions. But, I used the wave amplitudes to be in the normal direction. This was unstable. However, when I used the FL and FR to be along the normal direction, the problem was solved. I dont know if it was a bug in my code, which I rectified without my knowledge, or the treatment of fluxes made a difference. Anyway, for now, I guess it works fine.

 April 20, 2007, 15:50 Re: Riemann Solvers in multidimension #4 e f toro Guest   Posts: n/a I am surprised about HLLC not working for non Cartesian meshes. For a rotational invariant system, you first rotate the data in the direction normal to the face, calculate the HLLC, or any other, and then rotate back. See chapt 16 of Toro E F. Riemann solbvers and numerical methods for fluid dynamics. Springer 1999. EF Toro

 April 26, 2007, 04:05 Re: Riemann Solvers in multidimension #5 DSS Guest   Posts: n/a Hi, Thanks for your input. I corrected my code, and it works!! The book also gives multiple options for wave speed estimates, and I find that for a solution containing shock, these wave speed estimates play a significant role. Can you suggest me a good wave speed estimate? Regards, DSS

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post megacrout OpenFOAM 5 July 8, 2011 10:54 pbrane11 Fluent UDF and Scheme Programming 2 January 27, 2010 17:27 shuo Main CFD Forum 1 November 27, 2009 11:45 roadracer Main CFD Forum 1 September 6, 2008 05:56 sanjeev Main CFD Forum 3 July 22, 2000 16:59

All times are GMT -4. The time now is 20:00.