# Boundary condition ambiguous at edge vertices

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

 January 11, 2001, 13:50 Boundary condition ambiguous at edge vertices #1 Matthias Stripf Guest   Posts: n/a Hi, I am working on a 3D Control Volume based FEM code that uses an unstructured tetrahedral mesh. Unknowns are located at the mesh vertices. Now I came upon the following problem, setting the values at a boundary: Consider a simple channel with 1 inlet, 4 walls and 1 outlet. At the inlet, v=(1,0,0). Now the vertices at the edges that form the inlet belong to two different boundary faces. Each of them to a wall _and_ to an inlet. Which is now the correct velocity value to set at these vertices? The ones from the inlet (v=(1,0,0)) or the ones from the wall (v=(0,0,0))? Regards, Matthias

 January 11, 2001, 14:32 Re: Boundary condition ambiguous at edge vertices #2 Alton J. Reich Guest   Posts: n/a Matthias, The correct BC for any node that is on a wall is a wall (non-slip) BC. The way many codes handle BC assignement is to set the wall BCs last. That way any node on a wall (even if it is at a junction of a wall and something else) gets a wall BC. Regards, Alton

 January 12, 2001, 00:13 Re: Boundary condition ambiguous at edge vertices #3 Dan Williams Guest   Posts: n/a You will run into all kinds of difficulties like this trying to set boundary conditions directly on a vertex. Set boundary values on boundary faces and when you discretise the control volume boundary flux use the set value. On the inlet it will be the specified velocity, and on the wall the velocity is zero. This is unambiguous. Dan.

 January 12, 2001, 04:13 Re: Boundary condition ambiguous at edge vertices #4 Matthias Stripf Guest   Posts: n/a Yes, I agree that I have to set the boundary values on the boundary faces and I also use this value to compute the boundary flux. But what do I do with the other control volume faces? I use Gauss quadrature to do the integration on the other control volume faces and I need the values at the CV vertices to do so. So what do I do here? Matthias.

 January 12, 2001, 12:50 Re: Boundary condition ambiguous at edge vertices #5 Ghanshyam Singh Guest   Posts: n/a At the inlet edges, the right b.c. will be v=(1,0,0). Where as at the exit edges, it will be v=(0,0,0). At the inlet if you are interested in uniform velocity profile, you have to specify v(1,0,0) otherwise there will be a deficit in the mass. But if you are specifying a profile, e.g. a parabolic, v=(0,0,0) will have to be specified. Hope, it may help. GS

 January 13, 2001, 00:52 Re: Boundary condition ambiguous at edge vertices #6 Dan Williams Guest   Posts: n/a A boundary vertex value is "solved for", just like any other interior vertex value. You only use the boundary condition to close the flux through the boundary face. If you are running implicit you end up with coefficients (from your gauss quadrature) which multiply that vertex value at the next (implicit) time level. If you are running explicit you just multiply out the coefficients times the previous (solved for) value. The "solved for" values simply satisfy conservation for the control volume on the edge between the wall and the inlet. On the first timestep the vertex value will just be your initial guess. Dan.

 January 13, 2001, 03:07 Re: Boundary condition ambiguous at edge vertices #7 Matthias Stripf Guest   Posts: n/a Thanks! That's a good hint. I will give it a try! Matthias

 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 yating9901 FLUENT 3 June 28, 2010 12:26 HervĂ© CFX 2 June 16, 2008 06:40 Mark CFX 6 November 15, 2004 16:55 Tudor Miron CFX 15 April 2, 2004 06:18 boing Main CFD Forum 1 January 6, 2002 17:53

All times are GMT -4. The time now is 19:43.