CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

interDymFOAM with dynamicRefineFvMesh for DTMB5415 ship

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 15, 2018, 11:27
Default interDymFOAM with dynamicRefineFvMesh for DTMB5415 ship
  #1
New Member
 
Luc Bordier
Join Date: Feb 2010
Posts: 11
Rep Power: 13
lbordier is on a distinguished road
Hi,

I am trying to use dynamicRefineFvMesh within a interDyMFoam computation to refine based on free surface (alpha.water field).

When my computation is in PISO mode, the refinement is done but then, pressure number of iteration increase and forces on body become instable such that finally computation crashes.

When my computation runs in PIMPLE mode, solver crashes just after the first solve of p_rgh :

Code:
PIMPLE: iteration 1
Selected 377308 cells for refinement out of 11056721.
Refined from 11056721 to 13697877 cells.
Selected 0 split points out of a possible 377308.
Execution time for mesh.update() = 18.03 s
GAMG:  Solving for pcorr, Initial residual = 1, Final residual = 9.78769e-05, No Iterations 74
time step continuity errors : sum local = 5.99048e-12, global = -1.6006e-14, cumulative = 3.108e-10
smoothSolver:  Solving for alpha.water, Initial residual = 1.8908e-08, Final residual = 8.08454e-11, No Iterations 1
Phase-1 volume fraction = 0.75008  Min(alpha.water) = -1.56267e-05  Max(alpha.water) = 1.16759
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase-1 volume fraction = 0.75008  Min(alpha.water) = -0.00340939  Max(alpha.water) = 1.16624
GAMG:  Solving for p_rgh, Initial residual = 0.0813079, Final residual = 0.000809642, No Iterations 7
[13] #0  Foam::error::printStack(Foam::Ostream&)[7] #0  Foam::error::printStack(Foam::Ostream&) addr2line failed
[13] #1  Foam::sigFpe::sigHandler(int) addr2line failed
[13] #2  ? in "/home/applis/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
[7] #1  Foam::sigFpe::sigHandler(int) addr2line failed
[13] #3   in "/home/applis/OpenFOAM/OpenFOAM-4.0/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
[7] #2  ??
[13] #4   in "/lib64/libc.so.6"
[7] #3  ?
[13] #5  ?
[13] #6  __libc_start_main? addr2line failed
[13] #7  Foam::UOPstream::write(float)
My fvSolution is as following :
Code:
// ************************************************************************* //
[bordier@visu: dtmb5415]$ cat system/fvSolution
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.0                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    "alpha.water.*"
    {
        nAlphaCorr      2;
        nAlphaSubCycles 1;
        cAlpha          1;
        icAlpha         0;

        MULESCorr       yes;
        nLimiterIter    10;
        alphaApplyPrevCorr  yes;

        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-8;
        relTol          0;
        minIter         1;
    }

  /*  p_rgh
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-7;
        relTol          0.01;
    };*/

   p_rgh
    {
        solver          GAMG;

        smoother        DIC;

        tolerance       1e-7;
        relTol          0.01;
        cacheAgglomeration false;
    };

    p_rghFinal
    {
        $p_rgh;
        relTol          0;
    }

    "pcorr.*"
    {
        $p_rghFinal;
        relTol          0.0001;
    }


    "(U|k|omega).*"
    {
        solver          smoothSolver;

        smoother        symGaussSeidel;
        nSweeps         1;

        tolerance       1e-7;
        relTol          0.1;
        minIter         1;
    };
}

PIMPLE
{
    momentumPredictor   no;

    nOuterCorrectors    25;
    nCorrectors         2; //2
    nNonOrthogonalCorrectors 0;

    maxCo               15;
    maxAlphaCo          5;

    rDeltaTSmoothingCoeff 0.05;
    rDeltaTDampingCoeff 0.5;
    nAlphaSpreadIter    0;
    nAlphaSweepIter     0;
    maxDeltaT           0.005;

    residualControl
    {
        p_rgh
        {
                relTol 0;
                tolerance 5e-5;
        }

        U
        {
                relTol 0;
                tolerance 1e-4;
        }
    }
}

relaxationFactors
{
    fields
    {
        p_rgh 0.4;
        p_rghFinal 0.4;
    }
    equations
    {
        U 0.6;
        UFinal 0.6;
    }
}

cache
{
    grad(U);
}


// ************************************************************************* //
My dynamicMeshDict is :

Code:
dynamicRefineFvMeshCoeffs
{
    // How often to refine
    refineInterval  5;

    // Field to be refinement on
    field           alpha.water;

    // Refine field inbetween lower..upper
    lowerRefineLevel 0.001;
    upperRefineLevel 0.999;

    // If value < unrefineLevel unrefine
    unrefineLevel   10;

    // Have slower than 2:1 refinement
    nBufferLayers   2;

    // Refine cells only up to maxRefinement levels
    maxRefinement   1;

    // Stop refinement if maxCells reached
    maxCells        30000000;

    // Flux field and corresponding velocity field. Fluxes on changed
    // faces get recalculated by interpolating the velocity. Use 'none'
    // on surfaceScalarFields that do not need to be reinterpolated.
    correctFluxes
    (
        (phi none)
        (nHatf none)
        (rhoPhi none)
        (alphaPhi none)
        (ghf none)
        ("(alphaPhi-flux(alpha.water))" none)
    );

    // Write the refinement level as a volScalarField
    dumpLevel       true;
}
Does anybody have an idea how to solve this issue ?
Any help is appreciated.
lbordier is offline   Reply With Quote

Reply

Tags
dynamic, interdymfoam, refine, ship

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
interDymFoam: the sinking ship smiles wouter OpenFOAM Running, Solving & CFD 10 March 15, 2021 19:24
Full-scale ship case with velocity using interDyMFoam show unrealistic trim angle!! lil_star OpenFOAM 0 December 4, 2017 14:00
Long simulation time of interFoam and interDyMFoam; ship hull simulation Mueye OpenFOAM Running, Solving & CFD 0 December 7, 2016 04:13
Problem with p_rgh with interDyMFoam in a ship roll decay test karmo OpenFOAM Running, Solving & CFD 1 August 23, 2012 18:07
OpenFOAM15dev interDyMFoam problem using dynamicRefineFvMesh eberberovic OpenFOAM Bugs 6 January 14, 2010 06:18


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