CFD Online Discussion Forums - RBF motion solver implementation

Hi all,

I am currently simulating a foil with sinusoidal pitching and heaving motion using dynamic mesh method. I am using Solid Body Rotation (SBR) Stress dynamic mesh motion solver, and I used the Laplacian dynamic mesh before as well. I was able to conduct the simulation and validate the results with the non-inertial reference frame method.

However, I would like to increase the mesh resolution of the current simulation, and both the SBR and Laplacian solvers are not robust enough, and mesh cells got too skewed. I would like to try the Radial Basis Function (RBF) motion solver and see whether it can cope with the large mesh cell deformation. I am looking into the tutorial and also source codes, but I am not sure quite what are the variables in the dynamicMeshDict refer to. Could someone show some insights or recommend some documentations?

I know there is the PhD thesis by Frank Bos, but it does not go into details about the implementation in OpenFOAM.

The following is the link to the tuturial case:
http://sourceforge.net/p/openfoam-ex...ovingBlockRBF/

The dynamicMeshDict is as follows:

Code:

/*--------------------------------*- C++ -*----------------------------------*\

| =========                 |                                                 |

| \\      /  F ield         | OpenFOAM Extend Project: Open Source CFD        |

|  \\    /   O peration     | Version:  1.6-ext                               |

|   \\  /    A nd           | Web:      www.extend-project.de                 |

|    \\/     M anipulation  |                                                 |

\*---------------------------------------------------------------------------*/

FoamFile

{

    version     2.0;

    format      ascii;

    class       dictionary;

    object      dynamicMeshDict;

}

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



dynamicFvMesh dynamicMotionSolverFvMesh;



solver RBFMotionSolver;



movingPatches ( block );



staticPatches ( left right top bottom );



coarseningRatio 5;



includeStaticPatches   no;

frozenInterpolation    yes;



interpolation

{

    RBF  IMQB;

    focalPoint (0 0 0);

    innerRadius 2.5;//5.0

    outerRadius 12.5;

    polynomials true;



    W2Coeffs

    {

        radius     1.0;

    }

    TPSCoeffs

    {

        radius     5.0;

    }

    GaussCoeffs

    {

        radius     0.1;

    }

    IMQBCoeffs

    {

        radius     0.001;

    }

}





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

Thank you very much!