[alicezanella]

Hi everyone,

I'm having some troubles while deforming my mesh during an FSI simulation. The test case is a 2D naca0012 pitching and plunging. I'm using k-omega sst with wall functions and my mesh is unstructured + BL. The problem is that whatever displacementLaplacian method I use, I always end up with a distorsion near the trailing edge and on the pressure side. The crash happens within the first 10-20 timesteps, depending on which diffusivity I choose.

This is my dynamicMeshDict:

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      dynamicMeshDict;
}

dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs ("libfvMotionSolvers.so");

solver      displacementLaplacian;
// OpenFOAM9 or newer: rename "solver" to "motionSolver"

displacementLaplacianCoeffs {
//diffusivity quadratic inverseDistance 1 (airfoil);
// diffusivity  inverseDistance (airfoil);
diffusivity quadratic inverseDistance 1 (airfoil);
// diffusivity     exponential 0.25 inverseDistance (airfoil);
}

I'm also attaching two images: undeformed mesh and last timestep deformation before crashing.

Do you have any ideas how to solve this problem?
Thank you!