[irengclenk]

Dear All,

I just noticed a new feature from waves2Foam, the GABC, and found that interesting. I would like to use that.

I already compiled the updated waves2Foam code, all is working perfectly with the relaxationZone boundary condition tutorial. I am using OpenFOAM-v2012. But when I run the GABC boundary condition I got this kind of error

*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2012                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : _7bdb509494-20201222 OPENFOAM=2012
Arch   : "LSB;label=32;scalar=64"
Exec   : setWaveField
Date   : Dec 17 2021
Time   : 15:11:57
Host   : coastal4
PID    : 43687
I/O    : uncollated
Case   : /home/u0142217/Documents/openfoam_secondyear/GABC_trial/regularWave
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

// using new solver syntax:
p_rgh
{
    solver          GAMG;
    tolerance       1e-10;
    relTol          0;
    smoother        DIC;
    nPreSweeps      0;
    nPostSweeps     2;
    nFinestSweeps   2;
    cacheAgglomeration true;
    nCellsInCoarsestLevel 10;
    agglomerator    faceAreaPair;
    mergeLevels     1;
}

// using new solver syntax:
p_rghFinal
{
    solver          GAMG;
    tolerance       1e-10;
    relTol          0;
    smoother        DIC;
    nPreSweeps      0;
    nPostSweeps     2;
    nFinestSweeps   2;
    cacheAgglomeration true;
    nCellsInCoarsestLevel 10;
    agglomerator    faceAreaPair;
    mergeLevels     1;
}


Reading g

Reading waveProperties

Reading waveProperties

Reading field alpha

Reading field U

Reading field p

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigSegv::sigHandler(int) at ??:?
#2  ? in /lib/x86_64-linux-gnu/libpthread.so.0
#3  Foam::HashTable<Foam::entry*, Foam::word, Foam::string::hash>::Iterator<true>::Iterator(Foam::HashTable<Foam::entry*, Foam::word, Foam::string::hash> const*, Foam::word const&) at ??:?
#4  Foam::dictionary::csearch(Foam::word const&, Foam::keyType::option) const at ??:?
#5  Foam::dictionary::subDict(Foam::word const&, Foam::keyType::option) const at ??:?
#6  Foam::celerityShapeFunctions::New(Foam::fvMesh const&, Foam::word const&, Foam::dictionary const&, Foam::Vector<double> const&, double const&) at ??:?
#7  Foam::gabcPressureRobinVFvPatchScalarField::gabcPressureRobinVFvPatchScalarField(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#8  Foam::fvPatchField<double>::adddictionaryConstructorToTable<Foam::gabcPressureRobinVFvPatchScalarField>::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#9  Foam::fvPatchField<double>::New(Foam::fvPatch const&, Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in ~/OpenFOAM/u0142217-v2012/platforms/linux64GccDPInt32Opt/bin/setWaveField
#10  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary::readField(Foam::DimensionedField<double, Foam::volMesh> const&, Foam::dictionary const&) in ~/OpenFOAM/u0142217-v2012/platforms/linux64GccDPInt32Opt/bin/setWaveField
#11  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) in ~/OpenFOAM/u0142217-v2012/platforms/linux64GccDPInt32Opt/bin/setWaveField
#12  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::readFields() in ~/OpenFOAM/u0142217-v2012/platforms/linux64GccDPInt32Opt/bin/setWaveField
#13  Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&, bool) in ~/OpenFOAM/u0142217-v2012/platforms/linux64GccDPInt32Opt/bin/setWaveField
#14  main in ~/OpenFOAM/u0142217-v2012/platforms/linux64GccDPInt32Opt/bin/setWaveField
#15  __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#16  ? in ~/OpenFOAM/u0142217-v2012/platforms/linux64GccDPInt32Opt/bin/setWaveField
Segmentation fault (core dumped)

What is possibly the cause of the error?

Thank you

Kind regards,

[ngj]

Hi Ivandito,


Have you tried to run any of the tutorials (e.g. tutorials/waveFoamGABC/regularWave)? Here, there is no problem with setWaveField, when I am running on 2106. I guess that your problem is that the wave properties are not correctly set.



I am looking forward to hear your experience, since the GABC is a very strong tool for many wave-type problems.


Kind regards


Niels

[irengclenk]

Hi Niels,

Thanks for your response. Yes I also installed it in OF2106 and it was working properly, then I checked the installation in my OF2012 again because I only had problem with OF2012.

I dont know wether this is the correct solution but it works.

So apparently, if I want to install wave2Foam in my OF2012, I need to comment the primitivePatchBase.H in all dep.C file from each folder in preProcessing, so that the preProcessing feature are all compiled. All tutorials are now working in OF2012 as well. I did comparison between analytical focused waves and what we have in the tutorial and it gives the correct result. So, I would say this solves my problem.

Now working on implementing waves-current with GABC

Thanks!

[ngj]

Hi Ivandito,


Hmmm, that definitely does not sound like a healthy approach to get it compiled. I do not have a version 2012 myself, but as long as it works for 2106, I am happy with that.


It is good to hear that you have made it work and I would love to hear some comments on the usability, accuracy, etc. To the better of my knowledge, you are the first to give the GABC a spin, where I am not sitting next to the person.


Good luck on the waves and current approach. It is not super simple, but it should be possible. Unfortunately, I have not managed to derive anything that I am happy with to this point.



Cheers


Niels

[MFS]

Hi irenglenk and Niels,


Have you managed to implement wave+current with GABC? I've been trying to implement "combinedWaves" but I didn't manage to get it running...


Kind regards,


Manu