[be_inspired]

Hi all,

Buid: 1.6-ext-959ec266ba5c
Solver: moveDynamicMesh

I have tried several motion solvers and always always I obtain a bad mesh from first iteration so something must be wrong from the beginning and I can not find the problem.

dynamicMeshDict

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.2.2                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      dynamicMeshDict;
}

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

dynamicFvMesh dynamicMotionSolverFvMesh;

motionSolverLibs ("libfvMotionSolvers.so");
motionSolverLibs ("libcellDecompositionMotionSolver.so");

//solver displacementSBRStress;   //displacementLaplacian;
//solver velocityComponentLaplacian z;

//solver laplaceCellDecomposition;

solver velocityLaplacian;

diffusivity  quadratic inversePointDistance 1(bottom);

//diffusivity banana;

frozenDiffusion off;

distancePatches
( 
bottom
);

/*displacementSBRStressCoeffs
{
    // diffusivity  uniform;
    // diffusivity  directional (1 200 0);
    // diffusivity  motionDirectional (1 1000 0);
    // diffusivity  file motionDiffusivity;
    diffusivity  quadratic inverseDistance 1(bottom);
}

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

pointMotionU includes the vertical distance ( 0 0 distance) between vertex of faces of the bottom patch to the real stl geometry. It has been calculated previously.

/*--------------------------------*- 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       pointVectorField;
    location    "1";
    object      Udist;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 0 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    side3
    {
        type            slip;
    }
    side1
    {
        type            slip;
    }
    side2
    {
        type            slip;
    }
    side4
    {
        type            slip;
    }
    bottom
    {
        type            fixedValue;
        value           nonuniform ( 
(0 0 375.797)
(0 0 376.355)
(0 0 376.763)
(0 0 377.228)
(0 0 377.773)
...

Surface Mesh is ok but I obtain always negative cells and so because first cells are not moved.

When using laplaceCellDecomposition, I create a motionU field where for the patch is needed to impose the vertex and face centre displacement field. I have calculated the distance between vertex of faces of the patch to the stl and distance from face centre of the patch to the stl.

/*--------------------------------*- 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       tetPointVectorField;
    location    "1";
    object      Udist;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 0 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    side3
    {
        type            slip;
    }
    side1
    {
        type            slip;
    }
    side2
    {
        type            slip;
    }
    side4
    {
        type            slip;
    }
    bottom
    {
        type            fixedValue;
        value           nonuniform ( 
(0 0 375.797)
(0 0 376.355)
(0 0 376.763)
(0 0 377.228)
....

Here is some part of the output file of the simulation:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: velocityLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Time = 1
DICPCG:  Solving for cellMotionUx, Initial residual = 1, Final residual = 5.88085e-05, No Iterations 1
DICPCG:  Solving for cellMotionUy, Initial residual = 1, Final residual = 5.37124e-05, No Iterations 1
DICPCG:  Solving for cellMotionUz, Initial residual = 1, Final residual = 9.71784e-05, No Iterations 1
volume continuity errors : volume = 1.22183e+12, old volume = 1.3896e+12, max error = 9.31323e-09, sum local = 6.84765e-15, global = 2.83277e-17
    Point usage OK.
    Upper triangular ordering OK.
    Topological cell zip-up check OK.
    Face vertices OK.
    Face-face connectivity OK.
    Mesh topology OK.
    Boundary openness (1.69157e-19 5.56167e-18 4.51055e-13) Threshold = 1e-06 OK.
 ***High aspect ratio cells found, Max aspect ratio: 1.19514e+14, number of cells 508118 Threshold = 1000
    Minumum face area = 0.000407749. Maximum face area = 223271.  Face area magnitudes OK.
 ***Zero or negative cell volume detected.  Minimum negative volume: -674665, Number of negative volume cells: 508072
    Mesh non-orthogonality Max: 179.999 average: 77.1779 Threshold = 70
   *Number of severely non-orthogonal faces: 78863.
 ***Number of non-orthogonality errors: 1541036.
 ***Error in face pyramids: 3065841 faces are incorrectly oriented.
 ***Max skewness = 125155, 9521 highly skew faces detected Threshold = 4
    Failed 5 mesh geometry checks.
    Failed 1 mesh checks.
ExecutionTime = 72.91 s  ClockTime = 77 s

Using solver displacementLaplacian, generating the pointDisplacement file:

Create time

Create dynamic mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
Selecting motion solver: displacementLaplacian
Selecting motion diffusion: quadratic
Selecting motion diffusion: inverseDistance
Time = 1
GAMG:  Solving for cellDisplacementx, Initial residual = 1, Final residual = 1.19725e-09, No Iterations 5
GAMG:  Solving for cellDisplacementy, Initial residual = 1, Final residual = 8.85574e-10, No Iterations 5
GAMG:  Solving for cellDisplacementz, Initial residual = 1, Final residual = 1.04615e-09, No Iterations 3
volume continuity errors : volume = 1.22183e+12, old volume = 1.3896e+12, max error = 1.45519e-09, sum local = 6.65964e-15, global = 3.49921e-17
    Point usage OK.
    Upper triangular ordering OK.
    Topological cell zip-up check OK.
    Face vertices OK.
    Face-face connectivity OK.
    Mesh topology OK.
    Boundary openness (-4.54089e-17 -2.27977e-18 4.52428e-13) Threshold = 1e-06 OK.
 ***High aspect ratio cells found, Max aspect ratio: 1.01474e+14, number of cells 15544 Threshold = 1000
    Minumum face area = 0.000657964. Maximum face area = 176660.  Face area magnitudes OK.
 ***Zero or negative cell volume detected.  Minimum negative volume: -329035, Number of negative volume cells: 15536
    Mesh non-orthogonality Max: 179.952 average: 21.6289 Threshold = 70
   *Number of severely non-orthogonal faces: 36111.
 ***Number of non-orthogonality errors: 46875.
 ***Error in face pyramids: 101385 faces are incorrectly oriented.
 ***Max skewness = 21738.6, 5074 highly skew faces detected Threshold = 4
    Failed 5 mesh geometry checks.
    Failed 1 mesh checks.

I have also followed the tutorial of SnakeRiver Canyon and it also give me a bad results. Maybe diffusivity field is not calculated correctly?

[pardoa]

Hello be_inspired,

Did you finally find a solution for your problem? I'm currently facing the same sort of issues.

Thanks!
Álvaro

[ffarella]

Same problem here? Did anyone solve this problem?