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

OLAFOAM solver stops working

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 5, 2017, 22:30
Question OLAFOAM solver stops working
  #1
New Member
 
Join Date: Dec 2017
Posts: 2
Rep Power: 0
Chubo is on a distinguished road
Dear Foamers,

I am new in OpenFOAM simulation. Recently, I am building a numerical wave flume by using OLAFOAM library, to simulate the motion of a wave energy converter (olaDyMFoam). The flume works well. But when I add a function to the controlDict to run-time sample alpha.water at fixed points, the solver "olaDyMFoam" only runs 9.7 sec and then stop updating the log file (log.olaDyMFoam). The odd thing is, this job will be still running until timeout without any errors.

I use OpenFOAM-3.0.1 and my job was calculated on a cluster. The controlDict setting is listed below. Many thanks for your help.

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

application     olaDyMFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         240;

deltaT          0.1;

writeControl    adjustableRunTime;

writeInterval   0.05;

purgeWrite      0;

writeFormat     ascii;

writePrecision  12;

writeCompression uncompressed;

timeFormat      general;

timePrecision   6;

runTimeModifiable yes;

adjustTimeStep  yes;

maxCo           .5;
maxAlphaCo      .5;
maxDeltaT       0.02; 

functions
{
    waveElevation
    {
        type                   sets;
        functionObjectLibs     ("libsampling.so");
        setFormat              raw;
        interpolationScheme    cellPoint;
        outputControl          outputTime; // Don't use timeStep
        outputInterval 1;
        sets
        (
            Probe_50m
            {
                type        uniform;
                axis        z;
                start       (50 0 2);
                end         (50 0 -2);
                nPoints     401;
            }
            Probe_100m
            {
                type        uniform;
                axis        z;
                start       (100 0 2);
                end         (100 0 -2);
                nPoints     401;
            }

        );
        fields
        (
           alpha.water
        );
    }
};


libs
(
    "libwaveGeneration.so"
    "libwaveAbsorption.so"
    "libsampling.so"
);
Chubo is offline   Reply With Quote

Old   December 5, 2017, 23:45
Default
  #2
Senior Member
 
Pablo Higuera
Join Date: Jan 2011
Location: Singapore
Posts: 399
Rep Power: 11
Phicau is on a distinguished road
Hi Chubo,

I have faced this problem too.

This is caused by a limitation of OpenFOAM sampling procedure. From what I remember, as the time and mesh deformation progress, the initial sampling points that were inside a given cell will belong to another one (or even to a face, which is less likely, but has also happened to me), and sampling will fail.

I suggest that you open a bug in the OpenFOAM system to see if they can do something to update the points at runtime for 'DyM' solvers.

https://bugs.openfoam.org/rules.php

Best,

Pablo
Phicau is offline   Reply With Quote

Old   December 6, 2017, 03:24
Default
  #3
New Member
 
Join Date: Dec 2017
Posts: 2
Rep Power: 0
Chubo is on a distinguished road
Hi Pablo,

Thank you for your information.

Not sure the problem will be solved if I change the interpolationScheme from 'cellPoint' to 'cellPointFace'. I will let you know the result.

Regards
Chubo
Chubo is offline   Reply With Quote

Reply

Tags
olafoam, sampling, solving

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
icoFoam Courant number growth problem tdof OpenFOAM Running, Solving & CFD 11 June 27, 2017 07:19
OF23x: const labelList& cZone() and const cellZone& cZone() not working A_Pete OpenFOAM Programming & Development 2 October 24, 2014 04:47
thobois class engineTopoChangerMesh error Peter_600 OpenFOAM 4 August 2, 2014 09:52
3d vof Smaras FLUENT 2 February 19, 2013 07:58
Getting too many iterations by velocity solving (aborting). Changing U - Solver? suitup OpenFOAM Running, Solving & CFD 0 January 20, 2010 08:45


All times are GMT -4. The time now is 10:37.