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

Axisymmetric Case 3rd Velocity Component

Register Blogs Community New Posts Updated Threads Search

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   June 4, 2014, 12:17
Default Axisymmetric Case 3rd Velocity Component
  #1
New Member
 
Join Date: Dec 2012
Posts: 12
Rep Power: 13
Mentalo is on a distinguished road
Hi Folks,

I'm solving an axisymmetric Pipe-Flow with wedge boundary conditions. On watching the residuals, I was wondering why why all 3 velocity components are solved (see Picture below). Due to the U_y velocity residual my simulation never converges.

The wedge angle is 2deg on both sides, so less than 5deg. My blockMesh looks like this:
Code:
convertToMeters 0.001;

vertices
(
    (0 0 -315)
    (5.25 -0.18333404 -315)
    (5.25 0.18333404 -315)
    (0 0 -315)
    (0 0 0)
    (5.25 -0.18333404 0)
    (5.25 0.18333404 0)
    (0 0 0)

    (0 0 1000)
    (5.25 -0.18333404 1000)
    (5.25 0.18333404 1000)
    (0 0 1000)

    (12.5 -0.43650962 0)
    (12.5 0.43650962 0)
    (12.5 -0.43650962 1000)
    (12.5 0.43650962 1000)

    (100 -3.4920769 0)
    (100 3.4920769  0)
    (100 -3.4920769 1000)
    (100 3.4920769  1000)

);

blocks
(
    hex (0 1 2 0 4 5 6 4) (3 1 158) simpleGrading (0.6 1 0.8)
    hex (4 5 6 4 8 9 10 8) (3 1 500) simpleGrading (0.6 1 1.5)
    hex (5 12 13 6 9 14 15 10) (3 1 500) simpleGrading (1.2 1 1.5)
    hex (12 16 17 13 14 18 19 15) (40 1 500) simpleGrading (1.5 1 1.5)

);

edges
(
);

boundary
(
    symmetry
    {
      type symmetryPlane;
      faces
    (

     (0 4 4 0)
     (4 8 8 4)
     
     );
    }

    inlet
    {
        type inlet;
        faces
        (           

     (2 0 0 1)
        
        );
    }

    pilot
    {
        type inlet;
        faces
        (           

     (13 6 5 12)
        
        );
    }


    coflow
    {
        type inlet;
        faces
        (           

     (17 13 12 16)
        
        );
    }

   surrounding
    {
      type patch;
      faces
    (

     (18 19 17 16)

     );

    }

    wall
    {
        type wall;
        faces
        (           

     (2 1 5 6)
        
        );
    }

    outlet
    {
        type outlet;
        faces
        (           

     (8 8 10 9)
     (10 15 14 9)
     (15 19 18 14)
        
        );
    }

    rotper1
    {
      type wedge;
      faces
    (

     (1 5 4 0)
     (5 9 8 4)
     (12 14 9 5)
     (16 18 14 12)

     );
    }
     
    rotper2
    {
      type wedge;
      faces
    (

     (6 2 0 4)
     (10 6 4 8)
     (15 13 6 10)
     (19 17 13 15)

     );
    }



 
);

mergePatchPairs
(
);


Regards
Philipp
Mentalo is offline   Reply With Quote

 


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
Swirl component in inflow BC for 2d axisymmetric flows vishalsacharya OpenFOAM Pre-Processing 2 November 2, 2018 05:33
Howe to get swirl velocity component in 2D axisymtric swirl modeling Perra Fluent UDF and Scheme Programming 2 March 6, 2018 02:14
access to velocity component in cylindrical coordinates swati_mohanty FLUENT 0 August 18, 2010 10:16
Free surface boudary conditions with SOLA-VOF Fan Main CFD Forum 10 September 9, 2006 12:24
Body force - Does it work? Jan Rusås CFX 5 August 27, 2002 09:50


All times are GMT -4. The time now is 06:08.