CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Triangulation condition for timeVaryingMappedFixedValue

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 27, 2009, 10:57
Default Triangulation condition for timeVaryingMappedFixedValue
  #1
Senior Member
 
sega's Avatar
 
Sebastian Gatzka
Join Date: Mar 2009
Location: Frankfurt, Germany
Posts: 729
Rep Power: 20
sega is on a distinguished road
Hello World.

I'm currently trying to understand the check for a correct local coordinate system with the use of timeVaryingMappedFixedValue boundary condition.

For simplicity I will post the relevant part of the code.
The complete code is here: http://foam.sourceforge.net/doc/Doxy...8C_source.html

Code:
    const point& p0 = samplePoints[0];

    // Find point separate from p0
    vector e1;
    label index1 = -1;

    for (label i = 1; i < samplePoints.size(); i++)
    {
        e1 = samplePoints[i] - p0;

        scalar magE1 = mag(e1);

        if (magE1 > SMALL)
        {
            e1 /= magE1;
            index1 = i;
            break;
        }
    }

    // Find point that makes angle with n1
    label index2 = -1;
    vector e2;
    vector n;

    for (label i = index1+1; i < samplePoints.size(); i++)
    {
        e2 = samplePoints[i] - p0;

        scalar magE2 = mag(e2);

        if (magE2 > SMALL)
        {
            e2 /= magE2;

            n = e1^e2;

            scalar magN = mag(n);

            if (magN > SMALL)
            {
                index2 = i;
                n /= magN;
                break;
            }
        }
    }
Now my questions are:
  • the condition for a "pass through" seems to be
    magE1 > SMALL and magN > SMALL ?
    Now I haven't found any definition of what SMALL may be?
    What is it? Where is it defined?
  • What kind of operation is e1 /= magE1 ?
  • What is the break command doing?
  • What is n = e1^e2 doing? This is not a power-operation, right?!
__________________
Schrödingers wife: "What did you do to the cat? It's half dead!"
sega is offline   Reply With Quote

Reply


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
CFX Solver : Sudden crash Hervé CFX 2 June 16, 2008 07:40
Airfoil boundary condition Frank Main CFD Forum 1 April 21, 2008 19:36
Slip boundary condition what is inside normunds OpenFOAM Running, Solving & CFD 2 June 4, 2007 07:45
Pressure Boundary Condition Matt Umbel Main CFD Forum 0 January 11, 2002 11:06
The Boundary Condition about the Flat Plate boing Main CFD Forum 1 January 6, 2002 17:53


All times are GMT -4. The time now is 14:56.