CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Community Contributions (https://www.cfd-online.com/Forums/openfoam-community-contributions/)
-   -   [waves2Foam] Problems in Operating waveDyMFoam in openfoam 2.3.1 (https://www.cfd-online.com/Forums/openfoam-community-contributions/207966-problems-operating-wavedymfoam-openfoam-2-3-1-a.html)

zhxter May 18, 2015 02:26

Problems in Operating waveDyMFoam in openfoam 2.3.1
 
Dear all,
I am now calculating the dynamic response of a rectangular box in waves. But I get some problems when using the waveDyMFoam application.

Code:

zhxt@ubuntu:~/waves2Foam/tutorials/waveFoam/2Dgreenwatermove$ waveDyMFoam
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 2.3.1-bcfaaa7b8660
Exec  : waveDyMFoam
Date  : May 17 2015
Time  : 21:51:35
Host  : "ubuntu"
PID    : 49867
Case  : /home/zhxt/waves2Foam/tutorials/waveFoam/2Dgreenwatermove
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

Selecting dynamicFvMesh dynamicMotionSolverFvMesh
// using new solver syntax:
U
{
    solver          PBiCG;
    preconditioner  DILU;
    tolerance      1e-09;
    relTol          0;
}

// using new solver syntax:
UFinal
{
    solver          PBiCG;
    preconditioner  DILU;
    tolerance      1e-09;
    relTol          0;
}

// using new solver syntax:
gamma
{
    solver          PBiCG;
    preconditioner  DILU;
    tolerance      1e-07;
    relTol          0;
}

Selecting motion solver: sixDoFRigidBodyMotion
Translational constraint tensor (0 0 0 0 1 0 0 0 0)
Rotational constraint tensor (0 0 0 0 0 0 0 0 1)

Reading g

Reading waveProperties

PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Reading transportProperties

Selecting incompressible transport model Newtonian
Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Calculating field g.h

No finite volume options present

GAMG:  Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
GAMG:  Solving for pcorr, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 0, global = 0, cumulative = 0
Reading/calculating face velocity Uf

Courant Number mean: 0 max: 0

Starting time loop

Interface Courant Number mean: 0 max: 0
Courant Number mean: 0 max: 0
deltaT = 0.00119048
Time = 0.00119048

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::sixDoFRigidBodyMotion::updatePosition(double, double) at ??:?
#4  Foam::sixDoFRigidBodyMotionSolver::solve() at ??:?
#5  Foam::motionSolver::newPoints() at ??:?
#6  Foam::dynamicMotionSolverFvMesh::update() at ??:?
#7 
 at ??:?
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 
 at ??:?
Floating point exception (core dumped)\

The setup is listed as follows:
1, alpha.water:
Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  1.5-dev                              |
|  \\  /    A nd          | Revision: 1740                                  |
|    \\/    M anipulation  | Web:      http://www.OpenFOAM.org              |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    location    "0";
    object      alpha.org;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 0 0 0 0];

internalField  uniform 0;

boundaryField
{
    inlet
    {
        type            waveAlpha;
        refValue        uniform 0;
        refGrad        uniform 0;
        valueFraction  uniform 1;
        value          uniform 0;
    }
    bottom
    {
        type            zeroGradient;
    }
    outlet
    {
        type            zeroGradient;
    }
    atmosphere
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value          uniform 0;
    }
    vessel
    {
        type            zeroGradient;
    }
    frontBack
    {
        type            empty;
    }
}

2. pointDisplacement
Code:

FoamFile
{
    version    2.0;
    format      ascii;
    class      pointVectorField;
    location    "0";
    object      pointDisplacement;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 0 0 0 0 0];

internalField  uniform (0 0 0);

boundaryField
{
    //- Set patchGroups for constraint patches
    #include "${WM_PROJECT_DIR}/etc/caseDicts/setConstraintTypes"

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

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

    atmosphere
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }
   
    bottom
    {
        type            fixedValue;
        value          uniform (0 0 0);
    }

    vessel
    {
        type            calculated;
    }
}


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

3. P_rgh
Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  1.5-dev                              |
|  \\  /    A nd          | Revision: 1740                                  |
|    \\/    M anipulation  | Web:      http://www.OpenFOAM.org              |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    location    "0";
    object      pd;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [1 -1 -2 0 0 0 0];

internalField  uniform 0;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }
    bottom
    {
        type            zeroGradient;
    }
    outlet
    {
        type            zeroGradient;
    }
    atmosphere
    {
        type            totalPressure;
        U              U;
        phi            phi;
        rho            none;
        psi            none;
        gamma          1;
        p0              uniform 0;
        value          uniform 0;
    }
/*    vessel
    {
        type            fixedFluxPressure;
        gradient        uniform 0;
        value          uniform 0;
    }
*/
    vessel
    {
        type            zeroGradient;
    }
    frontBack
    {
        type            empty;
    }
}


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

4. U
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      volVectorField;
    location    "0";
    object      U.org;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField  uniform ( 0.0 0.0 0.0);

boundaryField
{
    inlet
    {
        type            waveVelocity;
        refValue        uniform ( 0 0 0 );
        refGradient    uniform ( 0 0 0 );
        valueFraction  uniform 1;
        value          uniform ( 0 0 0 );
    }
    bottom
    {
        type            fixedValue;
        value          uniform ( 0 0 0 );
    }
    outlet
    {
        type            fixedValue;
        value          uniform ( 0 0 0 );
    }
    atmosphere
    {
        type            pressureInletOutletVelocity;
        value          uniform ( 0 0 0 );
    }
    vessel
    {
        type            movingWallVelocity;
        value          uniform (0 0 0);
    }
    frontBack
    {
        type            empty;
    }
}


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

5. dynamicMeshDict
Code:

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

dynamicFvMesh    dynamicMotionSolverFvMesh;

motionSolverLibs    ("libsixDoFRigidBodyMotion.so");

solver              sixDoFRigidBodyMotion;

sixDoFRigidBodyMotionCoeffs
{
    patches        (vessel);
    innerDistance  0.05;
    outerDistance  0.11;

    centreOfMass    (6.29 0.9 0);
    mass            150;
    momentOfInertia (0 0 1.5);
    rhoInf          1;
    report          on;
    accelerationRelaxation 0.3;
    value          uniform (0 0 0);

    constraints
    {
        Heave
        {
            sixDoFRigidBodyMotionConstraint line;
            direction    (0 1 0);
        }
        Pitch
        {
            sixDoFRigidBodyMotionConstraint axis;
            axis          (0 0 1);
        }
    }
}

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

Could anyone help me with that ??

hchen May 18, 2015 05:04

Hi Xintao,

Just a question, waves2Foam now can be used in OpenFOAM 2.3.1? and also one more question, do you meet any problem when you use dynamicMotionSolverFvmesh in parallel?

Best,
Hao

zhxter May 18, 2015 21:57

Hi, hao

Yeah, I successfully compiled waves2foam on OF 2.3.1. And waveFoam solver works well. The only problem is about the solver waveDyMFoam, which still cannot work now.


All times are GMT -4. The time now is 16:58.