# Waves2Foam Related Topics

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

 October 9, 2013, 12:57 #581 Senior Member   Arthur Loginow Join Date: Aug 2012 Posts: 102 Rep Power: 4 Mr.Ngj would you mind if I send you a PM?

 October 10, 2013, 08:11 #582 Member   Join Date: Dec 2009 Posts: 42 Rep Power: 7 Ni Niels, I'm trying to simulate wave breaking/vertical jetting due to focusing waves event (very high steepness). So far, I have managed to go to steepness, Akc of 0.2 and 0.3 without breaking. However, as I reach Akc 0.4, the simulation stop with error code as follows : Code: ```MULES: Solving for alpha1 Phase-1 volume fraction = 0.699903 Min(alpha1) = -3.7219e-19 Max(alpha1) = 1.36268 DILUPBiCG: Solving for Ux, Initial residual = 0.00387945, Final residual = 1.89691e-10, No Iterations 4 DILUPBiCG: Solving for Uy, Initial residual = 0.00468833, Final residual = 6.7893e-10, No Iterations 3 GAMG: Solving for p_rgh, Initial residual = 0.0247348, Final residual = 8.56019e-08, No Iterations 8 GAMG: Solving for p_rgh, Initial residual = 8.55912e-08, Final residual = 8.55912e-08, No Iterations 0 time step continuity errors : sum local = 1.15524e-09, global = -6.6881e-10, cumulative = -5.17334e-07 GAMG: Solving for p_rgh, Initial residual = 0.000950645, Final residual = 6.40902e-08, No Iterations 6 GAMG: Solving for p_rgh, Initial residual = 6.38803e-08, Final residual = 6.38803e-08, No Iterations 0 time step continuity errors : sum local = 8.65e-10, global = 3.06962e-11, cumulative = -5.17304e-07 GAMG: Solving for p_rgh, Initial residual = 4.12761e-05, Final residual = 5.01357e-08, No Iterations 3 GAMG: Solving for p_rgh, Initial residual = 5.01163e-08, Final residual = 7.33679e-09, No Iterations 2 time step continuity errors : sum local = 9.93841e-11, global = -1.08676e-11, cumulative = -5.17315e-07 ExecutionTime = 11862.1 s ClockTime = 11913 s Courant Number mean: 0.00701962 max: 0.307835 Interface Courant Number mean: 9.04814e-05 max: 0.21481 deltaT = 0.00117681 Time = 22.3717 MULES: Solving for alpha1 Phase-1 volume fraction = 0.6999 Min(alpha1) = -1.76375e-19 Max(alpha1) = 1.36617 DILUPBiCG: Solving for Ux, Initial residual = 0.00313601, Final residual = 1.90061e-10, No Iterations 4 DILUPBiCG: Solving for Uy, Initial residual = 0.00375741, Final residual = 1.75008e-10, No Iterations 4 [0] #0 Foam::error::printStack(Foam::Ostream&) in "/home/ehk112/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #1 Foam::sigFpe::sigHandler(int) in "/home/ehk112/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #2 in "/lib64/libc.so.6" [0] #3 Foam::GAMGSolver::scale(Foam::Field&, Foam::Field&, Foam::lduMatrix const&, Foam::FieldField const&, Foam::UPtrList const&, Foam::Field const&, unsigned char) const in "/home/ehk112/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #4 Foam::GAMGSolver::Vcycle(Foam::PtrList const&, Foam::Field&, Foam::Field const&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::PtrList >&, Foam::PtrList >&, unsigned char) const in "/home/ehk112/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #5 Foam::GAMGSolver::solve(Foam::Field&, Foam::Field const&, unsigned char) const in "/home/ehk112/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so" [0] #6 Foam::fvMatrix::solveSegregated(Foam::dictionary const&) in "/home/ehk112/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so" [0] #7 Foam::fvMatrix::solve(Foam::dictionary const&) in "/home/ehk112/OpenFOAM/ehk112-2.2.x/platforms/linux64GccDPOpt/bin/waveFoam" [0] #8 [0] in "/home/ehk112/OpenFOAM/ehk112-2.2.x/platforms/linux64GccDPOpt/bin/waveFoam" [0] #9 __libc_start_main in "/lib64/libc.so.6" [0] #10 [0] in "/home/ehk112/OpenFOAM/ehk112-2.2.x/platforms/linux64GccDPOpt/bin/waveFoam" [cx1-5-15-2:16953] *** Process received signal *** [cx1-5-15-2:16953] Signal: Floating point exception (8) [cx1-5-15-2:16953] Signal code: (-6) [cx1-5-15-2:16953] Failing at address: 0x7a2b300004239 [cx1-5-15-2:16953] [ 0] /lib64/libc.so.6(+0x32920) [0x2ac14c5af920] [cx1-5-15-2:16953] [ 1] /lib64/libc.so.6(gsignal+0x35) [0x2ac14c5af8a5] [cx1-5-15-2:16953] [ 2] /lib64/libc.so.6(+0x32920) [0x2ac14c5af920] [cx1-5-15-2:16953] [ 3] /home/ehk112/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0xba) [0x2ac14b66be7a] [cx1-5-15-2:16953] [ 4] /home/ehk112/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_RNS1_IS8_EESD_h+0x1186) [0x2ac14b66ecc6] [cx1-5-15-2:16953] [ 5] /home/ehk112/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x3a8) [0x2ac14b670308] [cx1-5-15-2:16953] [ 6] /home/ehk112/OpenFOAM/OpenFOAM-2.2.x/platforms/linux64GccDPOpt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x118) [0x2ac1495bf168] [cx1-5-15-2:16953] [ 7] waveFoam(_ZN4Foam8fvMatrixIdE5solveERKNS_10dictionaryE+0x11c) [0x459c9c] [cx1-5-15-2:16953] [ 8] waveFoam() [0x48694e] [cx1-5-15-2:16953] [ 9] /lib64/libc.so.6(__libc_start_main+0xfd) [0x2ac14c59bcdd] [cx1-5-15-2:16953] [10] waveFoam() [0x4334d5] [cx1-5-15-2:16953] *** End of error message ***``` I'm not sure what the problem is. Is the error due to problem with the probe? At breaking point, alpha1 would have multiple values causing the interpolation operation for the probe to crash. Or is it due to the simulation has reach minimum time step? # UPDATE : A closer inspection for the alpha1 value from the log reveals the following. 1. For wave with steepness, Akc = 0.2, the alpha1 has a limit of 1 and 0 for every time step as follows: 2. For wave with steepness, Akc = 0.3, the alpha1 does not sustain the limit of 1 and 0 between 22 to 27 s. The simulation however managed to get stable and revert alpha1 to the limit of 1 and 0. 3. For wave reaching the breaking limit Akc=0.4, alpha1 shoot up higher than 1 and the simulation blows up (blue and red line correspond to minimum and maximum value of alpha1 respectively) Any idea on tackling this problem with unboundedness of alpha1 is highly appreciated. Kind regards, katakgoreng Last edited by katakgoreng; October 13, 2013 at 13:27.

 October 10, 2013, 14:21 #583 Senior Member   Arthur Loginow Join Date: Aug 2012 Posts: 102 Rep Power: 4 Hello, I am trying to install waveFoam, so far I have some doubts: 1. Which is the best OF version to work with waveFoam? 2. I did the step 3 from the waves2Foam wiki (This is copy the interFoam source code and modify it) 3. Now I have to proceed with the step number 8 (Is this correct or should I do this step before the step number 3?) 4.At this point I am not sure what to do when it says: -Obtain the source code via SVN as described above -Execute the Allwmake script in the folder waves2Foam What really means by source code?

 October 11, 2013, 14:48 #584 Senior Member   Arthur Loginow Join Date: Aug 2012 Posts: 102 Rep Power: 4 I made the changes to both the interFoam and the interDyMFoam files and when I try to finish the installation with the wmake command I get the following error: could not open file relaxationZone.H for source file waveDyMFoam.C could not open file readWaveProperties.H for source file waveDyMFoam.C SOURCE=waveDyMFoam.C ; g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -I.. -I/opt/openfoam221/src/transportModels/twoPhaseMixture/lnInclude -I/opt/openfoam221/src/transportModels -I/opt/openfoam221/src/transportModels/incompressible/lnInclude -I/opt/openfoam221/src/transportModels/interfaceProperties/lnInclude -I/opt/openfoam221/src/turbulenceModels/incompressible/turbulenceModel -I/opt/openfoam221/src/finiteVolume/lnInclude -I/opt/openfoam221/src/dynamicMesh/lnInclude -I/opt/openfoam221/src/dynamicFvMesh/lnInclude -I./../../../../../src/lnInclude -I/opt/openfoam221/src/meshTools/lnInclude -I/opt/openfoam221/src/fvOptions/lnInclude -I/opt/openfoam221/src/sampling/lnInclude -IlnInclude -I. -I/opt/openfoam221/src/OpenFOAM/lnInclude -I/opt/openfoam221/src/OSspecific/POSIX/lnInclude -fPIC -c \$SOURCE -o Make/linux64GccDPOpt/waveDyMFoam.o waveDyMFoam.C:44:28: fatal error: relaxationZone.H: No such file or directory compilation terminated. make: *** [Make/linux64GccDPOpt/waveDyMFoam.o] Error 1 Any help?

 October 12, 2013, 04:44 #585 Senior Member   Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Deltares, Delft, The Netherlands Posts: 1,608 Rep Power: 25 Hi Arthur, You say above that you are unsure on how to use SVN and obtain the source code. I suppose that you are experiencing problems with compiling waveDymFoam simply because you have not obtained the source code from the svn-repository. This is done by executing the command stated on the wiki, which starts with something like: Code: `svn co` After this step, you have to compile the source code as also described on the wiki by executing the Allwmake script as in Code: `./Allwmake` You do not need to make your own waveFoam solver for OF221, since it already comes along with waves2Foam, but if you need to use the moving meshes in OpenFoam, you need to create your own waveDymFoam solver. In that regard please note the NB-statement, which tells you that you should not do exactly as given on the wiki, as the described is partly outdated. Use the existing solvers to see, which modifications are required. Kind regards Niels __________________ Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.

 October 12, 2013, 13:39 #586 Senior Member   Arthur Loginow Join Date: Aug 2012 Posts: 102 Rep Power: 4 Neil thanks for your reply, when I run the ./Allwmake command I get the following error: ===================================== ENVIRONMENTAL VARIABLES ===================================== WAVES_XVERSION=0 WM_PROJECT_VERSION_NUMBER=221 WAVES_SOL=/home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam/applications/solvers/solvers221 WAVES_POST=/home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam/applications/utilities/postProcessing WAVES_UTIL=/home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam/applications/utilities EXTBRANCH=0 WAVES_APPBIN=/home/arthur/OpenFOAM/arthur-2.2.1/platforms/linux64GccDPOpt/bin WAVES_TUT=/home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam/tutorials WAVES_GSL_INCLUDE=/usr/include WAVES_GSL_LIB=/usr/lib64 WAVES_PRE=/home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam/applications/utilities/preProcessing WAVES_SRC=/home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam/src WAVES_DIR=/home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam WAVES_LIBBIN=/home/arthur/OpenFOAM/arthur-2.2.1/platforms/linux64GccDPOpt/lib FATAL ERROR. The directory path /home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam does not exist. Correct the path in bin/bashrc Please note that once bin/bashrc is created, bin/bashrc.org is only an inactive file. The latter is also the only of the two files, which is updated through the SVN-repository. EXITING I tried the same creating the /home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam Directory but still getting another error, what I am doing wrong?

 October 12, 2013, 15:08 #587 Senior Member   Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Deltares, Delft, The Netherlands Posts: 1,608 Rep Power: 25 Note that the error message does not advice you to create the given directory, but to adjust bin/bashrc, such that it is pointing to the location, where you have placed waves2Foam. Kind regards Niels __________________ Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.

 October 12, 2013, 18:33 #588 Senior Member   Arthur Loginow Join Date: Aug 2012 Posts: 102 Rep Power: 4 I was able to compile the application but now when I try to run the solver (waveFoam) at the tutorial (periodicSolitary) I get the following error: arthur@ubuntu:~/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam/tutorials/waveFoam/periodicSolitary\$ waveFoam /*---------------------------------------------------------------------------*\ | ========= | | | \\ / F ield | OpenFOAM: The Open Source CFD Toolbox | | \\ / O peration | Version: 2.2.1 | | \\ / A nd | Web: www.OpenFOAM.org | | \\/ M anipulation | | \*---------------------------------------------------------------------------*/ Build : 2.2.1-57f3c3617a2d Exec : waveFoam Date : Oct 12 2013 Time : 18:31:03 Host : "ubuntu" PID : 10165 Case : /home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam/tutorials/waveFoam/periodicSolitary nProcs : 1 sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE). fileModificationChecking : Monitoring run-time modified files using timeStampMaster allowSystemOperations : Disallowing user-supplied system call operations // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 --> FOAM FATAL IO ERROR: cannot open file file: /home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam/tutorials/waveFoam/periodicSolitary/system/fvSchemes at line 0. From function regIOobject::readStream() in file db/regIOobject/regIOobjectRead.C at line 87. FOAM exiting Any suggestion?

 October 13, 2013, 06:42 #589 Senior Member   Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Deltares, Delft, The Netherlands Posts: 1,608 Rep Power: 25 Hi Arthur, Execute the Allrun script. This will take care of cross-version compatibility of the tutorial cases. Kind regards Niels __________________ Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.

 October 13, 2013, 16:55 #590 Senior Member   Arthur Loginow Join Date: Aug 2012 Posts: 102 Rep Power: 4 Neils, thanks for your reply, I did what you said, I was able to generate the mesh without any problem, however I get this error when I try to run the solver: --> FOAM FATAL IO ERROR: Cannot find patchField entry for cyclic cyclic1_half0 Is your field uptodate with split cyclics? Run foamUpgradeCyclics to convert mesh and fields to split cyclics. file: /home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam/tutorials/waveFoam/periodicSolitary/0/p_rgh.boundaryField from line 25 to line 47. From function GeometricField::GeometricBoundaryField::readField(const DimensionedField&, const dictionary&) in file /opt/openfoam221/src/OpenFOAM/lnInclude/GeometricBoundaryField.C at line 181. FOAM exiting By this I understand that I have to specify the cyclic at the cyclic boundary, is this correct? Last edited by Maralady; October 13, 2013 at 18:53.

 October 14, 2013, 00:41 #591 Senior Member   Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Deltares, Delft, The Netherlands Posts: 1,608 Rep Power: 25 Good morning, I had forgotten that the definition of the cyclic patches has changed. Please search the forum to identify how to convert the old (pre OF2.0) format to the present formulation. There are plenty of threads on that topic - or simply read through the release notes on openfoam.org. Kind regards, Niels __________________ Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.

 October 14, 2013, 23:55 #592 Senior Member   Arthur Loginow Join Date: Aug 2012 Posts: 102 Rep Power: 4 Hello Niels, thanks for your reply, I have been looking the tutorials folders and found that at the 0 folder for the waveFlume all the files are .org, now I understand from the interFoam tutorial that I only need 1 .org file and that one is the alpha1 file, so my question is: It is the same with this case? Should all the files be regular files and should I have only 1 .org file (the alpha1.org file) When I try to use setFields I get the following error: --> FOAM FATAL IO ERROR: cannot open file file: /home/arthur/OpenFOAM/arthur-2.2.1/applications/utilities/waves2Foam/tutorials/waveFoam/waveFlume/system/setFieldsDict at line 0. From function regIOobject::readStream() in file db/regIOobject/regIOobjectRead.C at line 87. FOAM exiting However I can create the mesh and even check the mesh without having any warning Any suggestion??

 October 15, 2013, 00:25 #593 Senior Member   Arthur Loginow Join Date: Aug 2012 Posts: 102 Rep Power: 4 Just find out at this forum what I was doing wrong, the general process is: blockMesh, setWaveParameters, setWaveField, waveFoam My bad, I guess I will read the whole forum and re-read the wiki again before ask any other question But as 1 last thing Mr.Niels, I want to simulate 2D waves in different geometries with NO floating bodies on the water, what tutorials would you recommend me? (So far I think that relaxationZone and waveFlume are enough to work on my project but not quite sure) thanks for your time and THANKS for doing such a nice work with this application

October 16, 2013, 00:46
#594
New Member

Peter
Join Date: Oct 2013
Posts: 1
Rep Power: 0
Quote:
 Originally Posted by ngj Hi Sri, In your specific case it should read Code: `export WAVES_DIR=\$HOME/Desktop/waves2Foam` Kind regards Niels
Hi mate,
Where do you type that export... line? in the terminal? Or as a line in the bashrc file? Or as a line in the Allwmake file?
Peter

 October 17, 2013, 15:37 #595 Senior Member   Arthur Loginow Join Date: Aug 2012 Posts: 102 Rep Power: 4 Neils I need to change the wavesProerties file, and I was wondering if I change the depth of the water at this file, do I have to modify the geometry? Because I am not sure if those 2 are truly linked, hope you understand my question greetings

 October 22, 2013, 12:19 modelling waves crashing into a fixed structure #596 New Member   Imogen Wright Join Date: Oct 2013 Posts: 2 Rep Power: 0 I'd like to model open ocean waves in shallow water crashing into and passing a fixed wall for a personal project. I come from a computer science/physics background but am very new to CFD. Could anyone point me to a tutorial that would help with this? I've been playing with the 3Dwaves tutorial but I can't seem to draw a solid wall - however I set up the mesh and initial conditions the waves seem to go straight though my wall. What am I missing?

 October 22, 2013, 13:15 #597 Senior Member   Niels Gjoel Jacobsen Join Date: Mar 2009 Location: Deltares, Delft, The Netherlands Posts: 1,608 Rep Power: 25 Good evening, I will try to catch up after some time offline: @Peter: You should modify the file waves2Foam/bin/bashrc. @Arthur: Please read the posts by Katakgoreng above, as he recently asked the exact same question. @Imogen: I suppose that you are using OpenFoam-2.2? The 3Dwaves tutorial is not fully functioning under more recent versions, because the faceSet-utility was replaced by topoSet. This means that it is not possible to create the inner wall without modifying the tutorial. At least the part with faceSet needs to be modified. See the Allrun script for details on the execution. Kind regards, Niels __________________ Please note that I do not use the Friend-feature, so do not be offended, if I do not accept a request.

 October 22, 2013, 17:51 #598 New Member   Imogen Wright Join Date: Oct 2013 Posts: 2 Rep Power: 0 Dear Niels, Thank you so much for your prompt reply, and thank you for maintaining this incredibly useful tool! I've managed to convert the faceSet piece of the 3DWaves example to a working topoSet. However, when I run createBaffles I now get the following error: Code: ```--> FOAM FATAL IO ERROR: incorrect first token, expected

 October 23, 2013, 17:52 #599 Senior Member   Arthur Loginow Join Date: Aug 2012 Posts: 102 Rep Power: 4 It is Tsoft the period? How can I define the wave Period?

October 23, 2013, 21:57
#600
Member

Gregoire Junqua
Join Date: Jun 2011
Location: China
Posts: 58
Blog Entries: 1
Rep Power: 6
Quote:
 Originally Posted by Maralady It is Tsoft the period? How can I define the wave Period?
Hi
Nope it's should be set with the wave period depending your experiment but it's not the period

From the wiki : http://openfoamwiki.net/index.php/Contrib/waves2Foam
"Tsoft should be given a value of 0.0; Tsoft is typically set to the wave period as default and is not required by most of the wave theories. "

It is the same file waveProperties
Just after you would found
// Water depth at the boundary and in the relaxation zone depth 0.400000; // Period - needed by setWaveParameters -> yields omega period 2.0; // Cyclic wave frequency omega 3.141593;

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post hjasak OpenFOAM 2 October 26, 2013 04:33 Roberthealy1 CFD-Wiki 6 August 23, 2007 17:58 Jonas Larsson Main CFD Forum 3 February 9, 2007 11:11 vivekanand CFX 0 October 27, 2004 05:17 Antonio Filippone Main CFD Forum 0 August 28, 1999 12:16

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