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

dynamic mesh motion (angularOscillatingDisplacement) in OF1.6

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

Reply
 
LinkBack Thread Tools Display Modes
Old   May 12, 2010, 11:29
Default dynamic mesh motion (angularOscillatingDisplacement) in OF1.6
  #1
New Member
 
Tamás Baksa
Join Date: Dec 2009
Location: Budapest, Hungary
Posts: 2
Rep Power: 0
tamas.baksa.hu is on a distinguished road
Dear Foamers,

I'd like to simulate the flow around an oscillating cross-section. The shape is very simple it is almost like a rectangle. I had to define very fine boundary layer resolution in order to stay below y+<1. The chord length of my shape is 1 [m] and the first cell next to the wall is 0.04e-3 [m].

My constant/dynamicMeshDict looks like:
Code:
dynamicFvMesh     dynamicMotionSolverFvMesh;
solver                  displacementLaplacian;
motionSolverLibs    ("libfvMotionSolvers.so");
diffusivity             quadratic inverseDistance 1(fixedWalls);
// **************************************** //
I have defined 0/pointDisplacement like the following:
Code:
...
boundaryField
{
    [inlet, outlet, frontAndBack description]
    fixedWalls
    {
        type            angularOscillatingDisplacement;
        //  angle = angle0_ + amplitude_*sin(omega_*t.value());
        axis            (0 0 1);                //  direction
        origin          (0 0 0);
        angle0         0;
        amplitude     0.17453;              //  10 deg in rad
        omega         6.2832;               //  1Hz * (2pi)
        value           uniform (0 0 0);    //  ???
    }
}
// ********************************************* //
Using this configuration and moveDynamicMesh utililty severely non-orthogonal faces appear at angle = 0.59° and runs until the following error at angle = 1.83°:
Code:
Mesh non-orthogonality Max: 179.989 average: 7.45988
   *Number of severely non-orthogonal faces: 318.
 ***Number of non-orthogonality errors: 50.
 ***Error in face pyramids: 204 faces are incorrectly oriented.
In the next timestep OF exits and the mesh looks like in the following figure. The boundary layer mesh should remain parallel to the wall during the whole cycle.
http://img526.imageshack.us/img526/7...edmesha183.png

Do you have any idea how to avoid this?
1) Maybe using different diffusivity fields? According to Jasak and Tukovic paper [1] I should try diffusivity proportional to the distortion energy. But none of the available diffusivity keywords sound like dist. energy in OF1.6:
Code:
10
(
uniform
quadratic
file
motionDirectional
inverseDistance
inversePointDistance
inverseFaceDistance
directional
inverseVolume
exponential
)
2) Shall I define different sets for near field (BL mesh) and far field then rotate the near field set rigidly? I have found a similar case [2] but it was used in OF1.3. There is no subsetMotionSolverFvMesh in dynamicFvMesh keyword list in OF1.6:
Code:
5
(
dynamicMotionSolverFvMesh     // I'm using this one
solidBodyMotionFvMesh
dynamicRefineFvMesh
dynamicInkJetFvMesh
staticFvMesh                          //  mesh has no motion
)
Any idea is appreciated! Thanks in advance.

ps: my time step is sufficiently small, I think: 2e-4


[1] http://powerlab.fsb.hr/ped/kturbo/Op...tionFAMENA.pdf
[2] http://www.aero.lr.tudelft.nl/~frank/index.php?id=research/cfd/problems/2D_movingmesh_wing
tamas.baksa.hu is offline   Reply With Quote

Old   May 13, 2010, 05:43
Default
  #2
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 197
Blog Entries: 1
Rep Power: 9
egp is on a distinguished road
Why not look at 1.5-dev? subsetMotionSolverFvMesh is there for you to use, or at least learn what to do and then re-implement in 1.6.
egp is offline   Reply With Quote

Old   May 13, 2010, 10:53
Default
  #3
Senior Member
 
Kevin Smith
Join Date: Mar 2009
Posts: 103
Rep Power: 8
kev4573 is on a distinguished road
If you can conform your mesh around the body to a circle, you could use GGI to oscillate the body, which would avoid mesh deformation. Take a look at the mixerGGI example and code, it is fairly easy to change the driving function to a cosine signal.

If you don't want to change your mesh or simply cannot conform the mesh surrounding the body to a circle, I would look into the subset motion as Eric mentioned or RBF mesh motion.
kev4573 is offline   Reply With Quote

Old   May 13, 2010, 11:16
Default
  #4
egp
Senior Member
 
egp's Avatar
 
Eric Paterson
Join Date: Mar 2009
Location: Blacksburg, VA
Posts: 197
Blog Entries: 1
Rep Power: 9
egp is on a distinguished road
But Kevin, GGI is only in 1.5-dev, and the OP is using 1.6....
egp is offline   Reply With Quote

Old   May 13, 2010, 12:22
Default
  #5
Senior Member
 
Kevin Smith
Join Date: Mar 2009
Posts: 103
Rep Power: 8
kev4573 is on a distinguished road
Ah, you are correct.
kev4573 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
dynamic mesh on a hexa grid Manoj Kumar FLUENT 4 August 8, 2010 08:28
Turbulent problems during Dynamic Mesh Motion arturojortega FLUENT 2 June 2, 2009 10:33
Discrete Values for Mesh Motion UnderwaterAlex CFX 2 February 21, 2008 19:21
dynamic mesh on a hexa grid Manoj Kumar FLUENT 0 August 21, 2007 07:41
Dynamic mesh + grid adapt = Crash! (Files included BillH FLUENT 4 July 24, 2007 15:31


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