|
[Sponsors] |
Boundary condition ambiguous at edge vertices |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 11, 2001, 12:50 |
Boundary condition ambiguous at edge vertices
|
#1 |
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, 13:32 |
Re: Boundary condition ambiguous at edge vertices
|
#2 |
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 11, 2001, 23:13 |
Re: Boundary condition ambiguous at edge vertices
|
#3 |
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, 03:13 |
Re: Boundary condition ambiguous at edge vertices
|
#4 |
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, 11:50 |
Re: Boundary condition ambiguous at edge vertices
|
#5 |
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 12, 2001, 23:52 |
Re: Boundary condition ambiguous at edge vertices
|
#6 |
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, 02:07 |
Re: Boundary condition ambiguous at edge vertices
|
#7 |
Guest
Posts: n/a
|
Thanks! That's a good hint. I will give it a try!
Matthias |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
How exactly the "pressure outlet" bdry condition compute properties on the boundary? | yating9901 | FLUENT | 3 | June 28, 2010 12:26 |
CFX Solver : Sudden crash | Hervé | CFX | 2 | June 16, 2008 06:40 |
Convective Heat Transfer - Heat Exchanger | Mark | CFX | 6 | November 15, 2004 15:55 |
New topic on same subject - Flow around race car | Tudor Miron | CFX | 15 | April 2, 2004 06:18 |
The Boundary Condition about the Flat Plate | boing | Main CFD Forum | 1 | January 6, 2002 16:53 |