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

wallShearStress Monitoring

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

Like Tree1Likes
  • 1 Post By Mojtaba.a

Reply
 
LinkBack Thread Tools Display Modes
Old   June 5, 2013, 14:08
Default wallShearStress Monitoring
  #1
Senior Member
 
Mojtaba.a's Avatar
 
Mojtaba Amiraslanpour
Join Date: Jun 2011
Location: Zanjan, Iran
Posts: 233
Rep Power: 7
Mojtaba.a is on a distinguished road
Send a message via Yahoo to Mojtaba.a
Hi Guys,
I am simulating a flow over an external body in which calculating wallShearStress is my major effort. I want to be able to monitor this field, therefore I can stop simulations whenever the value of this field gets steady.
Right now i am monitoring U field without any problems using probes function object.
As you know wallShearStress is defined on surfaces, I don't know how to monitor specified points on particular surfaces.
Is there any function object similar to probes that can handle surfaces?

Thank you very much,
Mojtaba
immortality likes this.
__________________
Complex Heat & Flow Simulation Research Group
If you can't explain it simply, you don't understand it well enough. "Richard Feynman"
Mojtaba.a is offline   Reply With Quote

Old   June 8, 2013, 16:27
Default
  #2
Senior Member
 
Mojtaba.a's Avatar
 
Mojtaba Amiraslanpour
Join Date: Jun 2011
Location: Zanjan, Iran
Posts: 233
Rep Power: 7
Mojtaba.a is on a distinguished road
Send a message via Yahoo to Mojtaba.a
Any suggestions?
__________________
Complex Heat & Flow Simulation Research Group
If you can't explain it simply, you don't understand it well enough. "Richard Feynman"
Mojtaba.a is offline   Reply With Quote

Old   November 26, 2013, 12:43
Default
  #3
Senior Member
 
Artur's Avatar
 
Artur
Join Date: May 2013
Location: Southampton, UK
Posts: 226
Rep Power: 6
Artur is on a distinguished road
Hi,

Did you find an answer to this? I wrote a utility which writes out the wallShearStress at runtime using function objects and then tried using the patchCloud runtimeObject to extract the values over a set of points (see below) but it didn't work, I end up with no output files. However, if I tell it to probe pressure it works just fine...

If you haven't worked this out yet maybe together we can find a solution..

Best wishes,

Artur

Code:
samplingLine
{
    type            sets;
    functionObjectLibs ("libsampling.so");

    outputControl   outputTime;
    writeInterval   1600;

    setFormat   vtk;
    fields          ( p wallShearStress );

    interpolationScheme cellPoint;

    sets
    (
        plateProbePoints
        {
        type        patchCloud;
        axis        xyz;
        points
        (
            (0 0.01 0.05)
            (0.2 0.01 0.05)
            (0.4 0.01 0.05)
            (0.6 0.01 0.05)
            (0.8 0.01 0.05)
            (1 0.01 0.05)
        );
        maxDistance 0.1;    // maximum distance to search
        patches     (plate);
        }

    );
}
Artur is offline   Reply With Quote

Old   November 27, 2013, 05:14
Default
  #4
Senior Member
 
Mojtaba.a's Avatar
 
Mojtaba Amiraslanpour
Join Date: Jun 2011
Location: Zanjan, Iran
Posts: 233
Rep Power: 7
Mojtaba.a is on a distinguished road
Send a message via Yahoo to Mojtaba.a
Quote:
Originally Posted by Artur View Post
Hi,

Did you find an answer to this? I wrote a utility which writes out the wallShearStress at runtime using function objects and then tried using the patchCloud runtimeObject to extract the values over a set of points (see below) but it didn't work, I end up with no output files. However, if I tell it to probe pressure it works just fine...

If you haven't worked this out yet maybe together we can find a solution..

Best wishes,

Artur

Code:
samplingLine
{
    type            sets;
    functionObjectLibs ("libsampling.so");

    outputControl   outputTime;
    writeInterval   1600;

    setFormat   vtk;
    fields          ( p wallShearStress );

    interpolationScheme cellPoint;

    sets
    (
        plateProbePoints
        {
        type        patchCloud;
        axis        xyz;
        points
        (
            (0 0.01 0.05)
            (0.2 0.01 0.05)
            (0.4 0.01 0.05)
            (0.6 0.01 0.05)
            (0.8 0.01 0.05)
            (1 0.01 0.05)
        );
        maxDistance 0.1;    // maximum distance to search
        patches     (plate);
        }

    );
}
Hi,
Unfortunately I haven't found a solution yet. Well this is interesting. How have you been able to monitor a surface? You have monitored surface pressure successfully using this utility?
Can you upload your utility?
__________________
Complex Heat & Flow Simulation Research Group
If you can't explain it simply, you don't understand it well enough. "Richard Feynman"
Mojtaba.a is offline   Reply With Quote

Old   November 27, 2013, 06:08
Default
  #5
Senior Member
 
Artur's Avatar
 
Artur
Join Date: May 2013
Location: Southampton, UK
Posts: 226
Rep Power: 6
Artur is on a distinguished road
Quote:
Originally Posted by Mojtaba.a View Post
Hi,
Unfortunately I haven't found a solution yet. Well this is interesting. How have you been able to monitor a surface? You have monitored surface pressure successfully using this utility?
Can you upload your utility?
Hi,

The provided snippet, when put in the functions( ); in the controlDict will extract the pressure over a surface. The algorithm will snap the points from the specified locations to the nearest point on the surface and sample there (governed by the max distance). Say if my surface was at y=0 (which it is) it would sample at the x locations given and z=0.05. It won't do the wall shear stress though. Hence I modified the utility that calculates wallShearStess at runtime to dump the values to a file. I only coded it last night and can't spend too much time on it since it's just a side project and so it's not very well written and will dump the values for all cell centres of a specific patch (I need it for a 2D case and so that's good enough for me). If you need it only for monitoring convergence you can just pick a few points from there and look at their time trace and it should work just as well. See the dropbox link for a test case and the utility (I use OF 2.2.2) is attached..

https://www.dropbox.com/s/zuq0k17vhdik6bm/flatPlate.tgz

All the best,

Artur
Attached Files
File Type: gz wallBLRuntime.gz (5.1 KB, 5 views)
Artur 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
Monitoring Static Pressure on Rotor newbie384 CFX 0 March 30, 2013 18:20
wallShearStress with compressible flow Ruehri OpenFOAM Post-Processing 0 June 11, 2010 04:14
saving monitoring data mange FLUENT 2 March 3, 2009 06:54
monitoring cell Jane CD-adapco 2 March 4, 2004 22:01
monitoring variables Veera CD-adapco 1 February 12, 2001 18:19


All times are GMT -4. The time now is 21:10.