CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

probing vorticity at a point

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

Like Tree1Likes
  • 1 Post By gschaider

Reply
 
LinkBack Thread Tools Display Modes
Old   May 29, 2012, 10:57
Default probing vorticity at a point
  #1
New Member
 
Join Date: Mar 2009
Location: Pune, India
Posts: 27
Rep Power: 8
suryawanshi_nitin is on a distinguished road
Hi
Im doing 2D LES simulatin for flow over flat plate using pisoFoam algorithm. I want to monitor vorticity at a particular point, like we use probe function in openfoam for velocity/ pressure at point.
Even i tried with available post utilities feature vorticity in openfoam but it is giving me only vorticity values at point but without giving the corresponding x,y,z coordinates of that point... plz help me on this
I want to monitor vorticity within boundary layer about 40 points.

thanks in advance...

Regards
Nitin Suryawanshi
suryawanshi_nitin is offline   Reply With Quote

Old   May 29, 2012, 12:12
Default
  #2
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by suryawanshi_nitin View Post
Hi
Im doing 2D LES simulatin for flow over flat plate using pisoFoam algorithm. I want to monitor vorticity at a particular point, like we use probe function in openfoam for velocity/ pressure at point.
Even i tried with available post utilities feature vorticity in openfoam but it is giving me only vorticity values at point but without giving the corresponding x,y,z coordinates of that point... plz help me on this
I want to monitor vorticity within boundary layer about 40 points.
"only vorticity values at point": you mean you get a full field-file but you want only values at probe-locations? Then you might want to run the probeLocations-utility afterwards (it needs a dict in which you can specify the name of the field and the locations of the probes).

If you want to monitor the vorticity during the run you might want to try swak4Foam with this in controlDict:
Code:
functions {
    vortField
    {
        type expressionField;
        outputControl timeStep;
        outputInterval 1;
        fieldName vorticity;
        expression "curl(U)";
        autowrite true;        
    } 
    probeVort
    {
        type probes;
        fields ( 
           vorticity
        );
....
Maybe there is a functionObject in OF that has the functionality of the vorticity-utility then you can use that instead of the expressionField (and avoid using swak4Foam at all)
wyldckat likes this.
gschaider is offline   Reply With Quote

Old   June 16, 2012, 01:19
Default
  #3
New Member
 
Join Date: Mar 2009
Location: Pune, India
Posts: 27
Rep Power: 8
suryawanshi_nitin is on a distinguished road
Thanks for your reply sir and very very sorry for delayed reply as i was busy with other things.
I have tried with probeLocations sampling utility. it has worked well..
but for that first we should create vorticity field for each time step by typing vorticity in termianl within working directory in which time step wise results are there

Thanks again...
suryawanshi_nitin is offline   Reply With Quote

Old   May 13, 2014, 12:37
Default
  #4
Member
 
Join Date: Jun 2011
Posts: 65
Rep Power: 6
maalan is on a distinguished road
Quote:
If you want to monitor the vorticity during the run you might want to try swak4Foam with this in controlDict:
Code:
functions {
    vortField
    {
        type expressionField;
        outputControl timeStep;
        outputInterval 1;
        fieldName vorticity;
        expression "curl(U)";
        autowrite true;        
    } 
    probeVort
    {
        type probes;
        fields ( 
           vorticity
        );
....
Maybe there is a functionObject in OF that has the functionality of the vorticity-utility then you can use that instead of the expressionField (and avoid using swak4Foam at all)
I do exactly the same in order to compute the drag or lift coefficients for the flow past a square cylinder 2D via the vorticity formulation, i.e.:

Code:
    IpressureLift
    {
        type swakExpression;
        valueType patch;
        patchName SQUARE;
        accumulations (
            sum
        );
        expression "p*normal()&vector(0,1,0)*area()";
        verbose true;
    }

    IviscousLift
    {
        type swakExpression;
        valueType patch;
        patchName SQUARE;
        accumulations (
            sum
        );
        expression "-(2/200)*(vorticity^normal())&vector(0,1,0)*area()";
        verbose true;
    }
When I compared this results versus the forces obtained via the function forces in OF, I find differences in the viscous term... The pressure term is exactly the same.

I guess that OF calculates the viscosity term by using the viscous stress tensor, isn't it? Anyway, how could I improve my results bernhard?
maalan is offline   Reply With Quote

Old   May 13, 2014, 19:56
Default
  #5
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by maalan View Post
I do exactly the same in order to compute the drag or lift coefficients for the flow past a square cylinder 2D via the vorticity formulation, i.e.:

Code:
    IpressureLift
    {
        type swakExpression;
        valueType patch;
        patchName SQUARE;
        accumulations (
            sum
        );
        expression "p*normal()&vector(0,1,0)*area()";
        verbose true;
    }

    IviscousLift
    {
        type swakExpression;
        valueType patch;
        patchName SQUARE;
        accumulations (
            sum
        );
        expression "-(2/200)*(vorticity^normal())&vector(0,1,0)*area()";
        verbose true;
    }
When I compared this results versus the forces obtained via the function forces in OF, I find differences in the viscous term... The pressure term is exactly the same.

I guess that OF calculates the viscosity term by using the viscous stress tensor, isn't it? Anyway, how could I improve my results bernhard?
Yep. And part of the problem might be that expressionField does its calculations on the cells and the boundary values are basically zeroGradient. The stuff that the turbulence-model calculates (R etc) are available through function plugins (only to the internalField-parser but if I remember correctly with the right BCs). To use function plugins add the right library (for instance "libswakThermoTurbFunctionPlugin.so") to the libs-entry and the parser will know the new functions (also report a list of functions the first time it is used)
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   May 14, 2014, 01:58
Default
  #6
Member
 
Join Date: Jun 2011
Posts: 65
Rep Power: 6
maalan is on a distinguished road
Quote:
Yep. And part of the problem might be that expressionField does its calculations on the cells and the boundary values are basically zeroGradient. The stuff that the turbulence-model calculates (R etc) are available through function plugins (only to the internalField-parser but if I remember correctly with the right BCs). To use function plugins add the right library (for instance "libswakThermoTurbFunctionPlugin.so") to the libs-entry and the parser will know the new functions (also report a list of functions the first time it is used)
Thank you for your reply, bernhard, but I forgot to say that my case is laminar (Re_L=200). How do you think I could manage the vorticity issue on the square surface?? (I assume that the critical point is on the corners).

Thanks again!
maalan is offline   Reply With Quote

Old   May 14, 2014, 09:35
Default
  #7
Assistant Moderator
 
Bernhard Gschaider
Join Date: Mar 2009
Posts: 3,915
Rep Power: 40
gschaider will become famous soon enoughgschaider will become famous soon enough
Quote:
Originally Posted by maalan View Post
Thank you for your reply, bernhard, but I forgot to say that my case is laminar (Re_L=200). How do you think I could manage the vorticity issue on the square surface?? (I assume that the critical point is on the corners).

Thanks again!
Sorry. Haven't time to think it properly through, but my gut feeling says that it has something to do (or can be replaced) by snGrad(U) (that'd be the gradient of U in the direction of normal())
__________________
Note: I don't use "Friend"-feature on this forum out of principle. Ah. And by the way: I'm not on Facebook either. So don't be offended if I don't accept your invitation/friend request
gschaider is offline   Reply With Quote

Old   May 15, 2014, 05:46
Default
  #8
Member
 
Join Date: Jun 2011
Posts: 65
Rep Power: 6
maalan is on a distinguished road
Quote:
Sorry. Haven't time to think it properly through, but my gut feeling says that it has something to do (or can be replaced) by snGrad(U) (that'd be the gradient of U in the direction of normal())
Sorry but I am little confused... I'm trying to calculate the viscous term via the vorticity integral as shown:

Code:
    IviscousLift
    {
        type swakExpression;
        valueType patch;
        patchName SQUARE;
        accumulations (
            sum
        );
        expression "-(2/200)*(vorticity^normal())&vector(0,1,0)*area()";
        verbose true;
    }
Where the vorticity field was previously calculated in the next manner:

Code:
    
    vorticity
    {
        type expressionField;
        outputControl timeStep;
        fieldName vorticity;
        expression "curl(U)";
        verbose true;
        autowrite true;
    }
So, I don't understand why it has to do with snGrad(U).
Would you mind to explain a little more, bernhard?

Thanks!
maalan is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
problem in toposet Ahmed Khattab OpenFOAM Pre-Processing 10 March 24, 2013 16:18
error EOF in blockMesh Ahmed Khattab OpenFOAM Meshing & Mesh Conversion 7 May 17, 2012 00:37
block-structured mesh for t-junction Robert@cfd ANSYS Meshing & Geometry 20 November 11, 2011 05:59
Gmsh and samplesurface touf Open Source Meshers: Gmsh, Netgen, CGNS, ... 2 December 10, 2007 03:27
CFX4.3 -build analysis form Chie Min CFX 5 July 12, 2001 23:19


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