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

[swak4Foam] Dynamically morphing airfoil using groovyBC

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 20, 2022, 05:35
Default Dynamically morphing airfoil using groovyBC
  #1
Member
 
Join Date: Apr 2012
Location: Trivandrum
Posts: 34
Rep Power: 12
toolpost is on a distinguished road
Hello

I am trying to simulate 2D flow around a dynamically morphing airfoil using OpenFOAM-v2106. The airfoil surface near the trailing edge is supposed to morphed downwards with time. I believe this can be achieved using groovyBC. I set everything up and the mesh seems to move, but I am facing a weird issue.

The original C-type mesh is created in ICEM CFD, and exported to OpenFOAM. I want to keep the first cell height near the wall as constant throughout the simulation in order to keep a constant y+. But as you can see from the screenshots, in the very first iteration something weird happens and the mesh moves in an abnormal way. In the first time step, the airfoil wall didn't move, but the mesh around it moved. This results in a deformed mesh where the wall distance become reduced on the top side but keeps as expected in the bottom side. Later from the next time-step on-wards, the mesh deforms in the same manner I expected, and from thereon the wall distance no longer changes.

I cannot post the case file here, because the mesh size is large. A minimalistic cellDisplacementDict to reproduce this is given below.

Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      cellDisplacement;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 0 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    TOP
    {
           type          groovyBC;
         value        uniform (0 0 0);
          valueExpression    "vector(0,df,0)"; 
    
           variables
           (
        "df=-0.3*time();" 
           );
    }

    BOTTOM
    {
           type          groovyBC;
         value        uniform (0 0 0);
          valueExpression    "vector(0,df,0)"; 
    
           variables
           (
        "df=-0.3*time();" 
           );
    }
    
    FIXEDWALLS
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    INLET
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }

    OUTLET
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }

    frontAndBack
    {
        type            empty;
    }
}
This is dynamicMeshDict
Code:
FoamFile
{
    format      ascii;
    class       dictionary;
    location    "system";
    object      dynamicMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

 dynamicFvMesh   dynamicMotionSolverFvMesh;

 motionSolverLibs ("libfvMotionSolvers.so");

 motionSolver    displacementLaplacian;

 diffusivity  directional (20 0.1 0);
What am I doing wrong? Any suggestions would be greatly appreciated.
Attached Images
File Type: jpg timeZero.jpg (56.8 KB, 15 views)
File Type: jpg timedelt.jpg (53.5 KB, 16 views)
toolpost is offline   Reply With Quote

Old   May 20, 2022, 12:47
Default
  #2
Member
 
Join Date: Apr 2012
Location: Trivandrum
Posts: 34
Rep Power: 12
toolpost is on a distinguished road
Deleted.....

Last edited by toolpost; June 17, 2022 at 00:28.
toolpost is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with airfoil shape optimization robyTKD SU2 Shape Design 7 March 7, 2022 17:18
SU2 AOA optimization 454514566@qq.com SU2 9 March 7, 2022 17:17
Ffd_control_point_2d feiyi SU2 4 September 30, 2019 13:42
High drag for airfoil compared to XFOIL and wind tunnel data Ry10 SU2 15 October 30, 2016 18:27
2D FFD Optimization RLangtry SU2 2 August 5, 2014 10:48


All times are GMT -4. The time now is 07:12.