Finite volume inflow/outflow BCs

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

 June 5, 2008, 06:26 Finite volume inflow/outflow BCs #1 Niels Guest   Posts: n/a Hello, I am trying to figure out what boundary conditions to apply at inflow and outflow boundaries for my 3D finite volume solver. I came across this old paper: Ph. Rostand, B. Stoufflet, "Finite Volume Galerkin Methods for Viscous Gas Dynamics", INRIA, 1988 http://hal.inria.fr/inria-00075691/en/ Section II.3 mentions a mechanism whereby such BC's can be imposed by means of a split of the inviscid flux Jacobian A = dF/dU according to the sign of the eigenvalues: A = A+ + A- where A+ = R(D+)L, A- = R(D-)L and D+, D- contain only the non-negative and non-positive eigenvalues respectively. I realize that for a homogenous function of degree 1, as is the case for the Euler flux function, the fluxes can now be calculated as F(U) = (A+)U + (A-)U However, the paper says that the boundary fluxes at an inflow or outflow can be calculated as: F(U) = (A-)U_external + (A+)U_internal U_external is the state exterior to the computational domain, U_internal the state just inside the boundary and where the flux Jacobians are evaluated at U_internal. I don't understand. Intuitively, I would have it this way: F(U) = (A-)U_internal + (A+)U_external because A+ and A- correspond to characteristics entering and leaving the computational domain. For example, at an inflow I would like extererior information to propagate into the domain. I tried to find the background paper by Steger and Warming, but failed. So if anyone can explain why I am wrong, that would be great. Thanks, Martin

 June 5, 2008, 09:20 Re: Finite volume inflow/outflow BCs #2 Vinayender Guest   Posts: n/a Hi Niels, instead of "F(U) = (A-)U_external + (A+)U_internal " consider "F(U) = (A-)U_right_cell + (A+)U_left_cell " for time being, and convensions are face normal is defined such that it propogates in to the right_cell. Now see, INCELL (A+)-> |(face)| <-(A-) OUTCELL at a given face A- will always points towards face from outcell and A+ will always points towards face from incell Finally , at inflow, INCELL will be the exterior domain and hence you use A+, you are right here but you confused because in that formula (F(U) = (A-)U_exterior + (A+)U_interior) external cell means external to the face with the normal direction I mentioned and hence external cell what he is refering is OUTCELL which is a interior zone not the exterior zone with respect to the computational domain.......... You have just confused with the word external which does not mean exterior domain to the computatinal domain but cell which is right side to the face (with face normal that i defined) and all your "physics of upwinding" understanding is right. If you dont understand what I said please draw it on a paper while reading, and even then you dont understand it means that I dint expressed thing properly and I am sorry for that. Thanks, Vinayender.

 June 6, 2008, 01:34 Re: Finite volume inflow/outflow BCs #3 Vinayender Guest   Posts: n/a Hi, Sorru in above discription Incell corresponds to left cell and Outcell correspoond to right cell for a particular face. Thanks, Vinayender.

 June 6, 2008, 03:50 Re: Finite volume inflow/outflow BCs #4 Martin Guest   Posts: n/a Hi Vinayender, Thank you very much for your answer. It totally cleared things up! It seems that my text books, although they never explicitly mention it, consistently adhere to the convention of traversing boundaries so that the interior of the volume always lies on the left hand side of the boundary. Just like you say. Regards, Martin

 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 Ganesh FLUENT 13 January 22, 2014 05:11 gaottino OpenFOAM Native Meshers: blockMesh 7 July 19, 2010 14:11 paean OpenFOAM Running, Solving & CFD 0 November 14, 2008 22:14 SSL FLUENT 2 January 26, 2008 12:55 Rasmus Gjesing (Gjesing) OpenFOAM Native Meshers: blockMesh 10 April 2, 2007 14:00

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