|
[Sponsors] |
![]() |
![]() |
#61 |
Member
Hao Chen
Join Date: Aug 2014
Posts: 66
Rep Power: 0 ![]() |
I reinstall the system to make sure foam extend install correctly with openmpi.
But then I rerun irregular45degreeTanks, then I get: Code:
[hchen-OptiPlex-9020:22570] *** An error occurred in MPI_Recv [hchen-OptiPlex-9020:22570] *** on communicator MPI_COMM_WORLD [hchen-OptiPlex-9020:22570] *** MPI_ERR_TRUNCATE: message truncated [hchen-OptiPlex-9020:22570] *** MPI_ERRORS_ARE_FATAL: your MPI job will now abort -------------------------------------------------------------------------- mpirun has exited due to process rank 2 with PID 22569 on node hchen-OptiPlex-9020 exiting improperly. There are two reasons this could occur: 1. this process did not call "init" before exiting, but others in the job did. This can cause a job to hang indefinitely while it waits for all processes to call "init". By rule, if one process calls "init", then ALL processes must call "init" prior to termination. 2. this process called "init", but exited without calling "finalize". By rule, all processes that call "init" MUST call "finalize" prior to exiting or it will be considered an "abnormal termination" This may have caused other processes in the application to be terminated by signals sent by mpirun (as reported here). -------------------------------------------------------------------------- [hchen-OptiPlex-9020:22566] 3 more processes have sent help message help-mpi-errors.txt / mpi_errors_are_fatal [hchen-OptiPlex-9020:22566] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages Simulation complete. |
|
![]() |
![]() |
![]() |
![]() |
#62 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 ![]() |
Hi Hao,
regarding your first post, I think OpenFOAM is not correctly installed, there are some libraries that are not compiled, so dummy ones get called. Please, follow the instructions provided by OpenFOAM or found in this forum to install the version you have on your system. Regarding the second post, this is a known issue, check this out: http://openfoamwiki.net/index.php/Contrib/IHFOAM#Supported_Versions The solution is easy: edit $WM_PROJECT_DIR/etc/controlDict and change commsType to nonBlocking Best, Pablo |
|
![]() |
![]() |
![]() |
![]() |
#63 | |
Member
Hao Chen
Join Date: Aug 2014
Posts: 66
Rep Power: 0 ![]() |
Hi Pablo:
Thanks a lot! I should check the wikki page first. Best regards Hao Quote:
|
||
![]() |
![]() |
![]() |
![]() |
#64 | |
New Member
Bo Terp Paulsen
Join Date: Oct 2010
Posts: 13
Rep Power: 15 ![]() |
Quote:
Thank you for clarifying this, it is appreciated. Kind Regards, Bo Terp |
||
![]() |
![]() |
![]() |
![]() |
#65 |
New Member
Remi Carmi
Join Date: Jul 2014
Posts: 15
Rep Power: 12 ![]() |
Hi All,
Is IHFOAM supposed to work with backward time integration? If yes great! if not any explanation? I am getting a linear instability when I try. Thanks |
|
![]() |
![]() |
![]() |
![]() |
#66 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 ![]() |
Hi Remi,
yes, and no. There are no mathematical restrictions for IHFOAM working with negative time steps. However, it is not possible to generate waves that propagate in a direction outwards of the domain (not in IHFOAM, not in any other models I know about). Read the irregular waves section in the manual, that is why the cosine function was implemented. Without further information, my best guess is that this is the problem you are experiencing. Best, Pablo |
|
![]() |
![]() |
![]() |
![]() |
#67 |
New Member
Remi Carmi
Join Date: Jul 2014
Posts: 15
Rep Power: 12 ![]() |
Well I am still not using the wave generator. Just the solver for a standing wave.
So i start with a cosine deformation and a uniform 2D mesh and I set atmospher bc at the top and slip everywhere else (but empty at front and back of course). I use linear, linear upwind for the div(phi,u). And I change the refinement and the time step. I am not able to make the backward to work (being stable and actually better than the Euler) but with ridiculously small time step (courant smaller than 0.1!). ... Thanks Remi |
|
![]() |
![]() |
![]() |
![]() |
#68 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 ![]() |
Hi Remi,
if no wave generation/absorption or porosity is involved you can run the same case with interFoam/interDyMFoam and the results (and limitations) will be the same. Therefore, if you want to maximize the options of someone answering your questions, you should post them in a separate thread, as a number of people working with interFoam may not be following this particular thread. Best, Pablo |
|
![]() |
![]() |
![]() |
![]() |
#69 |
New Member
Join Date: Jul 2014
Location: Aachen, Germany
Posts: 3
Rep Power: 12 ![]() |
Hello,
Thanks a lot for developping the IHFoam tool and also offering such a detailed and straightforward documentation! Anyway, I got a question before starting to get too deep into the topic: Did I get it right that I can only use the IHFoam boundaries and solve with interFoam by solely linking the libs in the controlDict? Or doesn't this make any sense if I really want to adopt wave asorption at one of the boundaries? Thanks a lot for your help! |
|
![]() |
![]() |
![]() |
![]() |
#70 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 ![]() |
Hi Reni,
yes, you are right. If you don't need the porous media flow capabilities, you can link the wave generation / absorption libraries inside controlDict and use the regular interFoam / interDyMFoam solvers. Take a look at the wiki for the instructions. Best, Pablo |
|
![]() |
![]() |
![]() |
![]() |
#71 |
New Member
Join Date: Jul 2014
Location: Aachen, Germany
Posts: 3
Rep Power: 12 ![]() |
Hi Pablo,
great. Thanks for the fast reply! Best Verena |
|
![]() |
![]() |
![]() |
![]() |
#72 |
Member
António Pires
Join Date: Oct 2014
Posts: 33
Rep Power: 11 ![]() |
Greetings,
When i try to run the tutorial basewaveflume by running ./runCase i get these lines in termnial: blockMesh meshing... Preparing 0 folder... Setting the fields... Running... #0 Foam::error:: printStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linuxGccDPOpt/lib/libOpenFOAM.so" #2 Uninterpreted: #3 Elliptic::ellipticIntegralsKE(double, double*, double*) in "/home/antonio/OpenFOAM/antonio-2.2.2/platforms/linuxGccDPOpt/lib/libIHwaveGeneration.so" #4 cnoidalFun::calculations(double, double, double, double*, double*) in "/home/antonio/OpenFOAM/antonio-2.2.2/platforms/linuxGccDPOpt/lib/libIHwaveGeneration.so" #5 Foam::IH_Waves_InletAlphaFvPatchScalarField::updat eCoeffs() in "/home/antonio/OpenFOAM/antonio-2.2.2/platforms/linuxGccDPOpt/lib/libIHwaveGeneration.so" #6 Foam::fvPatchField<double>::evaluate(Foam::UPstrea m::commsTypes) in "/home/antonio/OpenFOAM/antonio-2.2.2/platforms/linuxGccDPOpt/bin/ihFoam" #7 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::evaluate() in "/home/antonio/OpenFOAM/antonio-2.2.2/platforms/linuxGccDPOpt/bin/ihFoam" #8 in "/home/antonio/OpenFOAM/antonio-2.2.2/platforms/linuxGccDPOpt/bin/ihFoam" #9 __libc_start_main in "/lib/i386-linux-gnu/libc.so.6" #10 in "/home/antonio/OpenFOAM/antonio-2.2.2/platforms/linuxGccDPOpt/bin/ihFoam" ./runCase: line 16: 2861 Floating point exception(core dumped) ihFoam > ihFoam.log Simulation complete. And it only creates de folder "0", no other time step is created and i cannot see any results. Anyone knows how to solve this? Maybe there is some mistake with the installation process? Thanks, António |
|
![]() |
![]() |
![]() |
![]() |
#73 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 ![]() |
Hi António,
apparently something is wrong with the elliptic functions. Can you submit your IHWavesDict file and water depth to check it myself? Have you made any other changes in the case? It looks very similar to what happened here: http://www.cfd-online.com/Forums/ope...tml#post503610 , can I have full details of your installation (e.g. Linux distribution, architecture...)? Best, Pablo Last edited by Phicau; November 3, 2014 at 10:52. Reason: Add reference to previous case |
|
![]() |
![]() |
![]() |
![]() |
#74 |
Member
António Pires
Join Date: Oct 2014
Posts: 33
Rep Power: 11 ![]() |
Thank you for the reply Pablo,
I'm using UBUNTU 12.04 and OpenFoam 2.2.2 I followed the installation instructions located at http://openfoamwiki.net/index.php/Contrib/IHFOAM I didn't make the step 2.4 but i think it is not needed for this tutorial am i right? Here are the files https://www.dropbox.com/sh/1exty2zf5...2rj4h9fma?dl=0 Thanks, António |
|
![]() |
![]() |
![]() |
![]() |
#75 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 ![]() |
Hi António,
64 bits system? It is very strange, because I cannot reproduce this error. The case runs flawlessly in my two computers (Ubuntu 12.04, 64 bits, OpenFOAM 2.2.2-9240f8b967db). Let's see if we can find this error, first reported by Dmitrjs. Can you test that the cnoidal and elliptical functions work on your machine? Download: https://www.dropbox.com/s/6aggsqfoc1...og.tar.gz?dl=0 and run: Code:
./compile ./testProgram Code:
Wave theory: cnoidal H: 0.1 T: 3 h: 0.4 L: 6.0043 m parameter: 0.9771 Elliptic integrals K: 3.28778 E: 1.03199 Best, Pablo |
|
![]() |
![]() |
![]() |
![]() |
#76 |
Member
António Pires
Join Date: Oct 2014
Posts: 33
Rep Power: 11 ![]() |
Hello Pablo,
I have a 32-bit system I don't think the functions are working, I ran the command ./compile and i got this lines: antonio@antonio-P5QL-PRO:~/IHFOAM/IHFOAM_materials$ ./compile waveFun.C: In function ‘int cnoidalFun::calculations(double, double, double, double*, double*)’: waveFun.C:602:16: warning: unused variable ‘KElliptic’ [-Wunused-variable] waveFun.C:603:16: warning: unused variable ‘EElliptic’ [-Wunused-variable] And after this, when i run ./testProgram i dont get any response from terminal.. |
|
![]() |
![]() |
![]() |
![]() |
#77 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 ![]() |
Dear António, dear all,
I think the problem was related to an infinite loop due to a hardcoded tolerance. The tolerance checks now the precision of the computer, so I hope to have fixed it. António, please update the code from git, following the instructions of the wiki. Then re-compile it and try the case again. I kindly ask you to report back to check if this really works. Best, Pablo |
|
![]() |
![]() |
![]() |
![]() |
#78 |
Member
António Pires
Join Date: Oct 2014
Posts: 33
Rep Power: 11 ![]() |
Hi Pablo,
I updated the code and everything is working perfectly now! I've only done the first tutorial but it went with no problem at all, time steps we're created as they should be. Thank you very much once again! António |
|
![]() |
![]() |
![]() |
![]() |
#79 |
Senior Member
Join Date: Jul 2011
Posts: 120
Rep Power: 15 ![]() |
Hi All, I am curious and would like to contribute some of my findings after several months of using IHFoam. Correct me if I'm wrong but as of now, there is no way of extracting the time series of the water depths of surface elevations at several positions of the domain built in to IHFoam?
I am currently using this in my controlDict file to extract the surface @ runtime: Code:
functions ( freeSurface { type surfaces; functionObjectLibs ( "libsampling.so" ); outputControl timeStep; outputInterval 10; surfaceFormat vtk; fields ( alpha1 ); surfaces ( freeSurface { type isoSurfaceCell; isoField alpha.water; isoValue 0.5; interpolate false; regularise false; } ); interpolationScheme cell; } ); Another issue is that since this is done at runtime, I have to run this again if I were to tune the isovalue. Does anyone know of a more elegant way of doing this as a post processing step instead? |
|
![]() |
![]() |
![]() |
![]() |
#80 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 ![]() |
Hi Haze,
thanks for your code that extracts free surface as a 3D surface along the whole domain. It has been around other threads, but it can be handy for users, as sometimes things are difficult to look for when you don't know that they exist. There is a way to run this function object after the case has finished, to work like the sample tool. You just need to run execFlowFunctionObjects. You are right, alpha = 0.5 is just an arbitrary convention, but that cell is really half air and half water, so there is plenty of water above it yet. In our papers we found out that to estimate run-up it is reasonable to consider alpha = 0.9. Moreover, there are lots of ways to obtain free surface elevation, we cover 4 or 5 in the training courses. Did you already register in http://ihfoam.ihcantabria.com/source-download/ ? We do provide tools to obtain free surface elevation at any location in the free materials that you can download when you register. Check the breakwater case. There you can find a procedure to sample the pressure and free surface elevation, and some python tools to convert the raw data into a time series. You also have a python tool to plot it. Best, Pablo |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
Divergence detected in AMG solver: k when udf loaded | google9002 | Fluent UDF and Scheme Programming | 3 | November 7, 2019 23:34 |
udf problem | jane | Fluent UDF and Scheme Programming | 37 | February 20, 2018 04:17 |
UDF velocity profile | willroca | Fluent UDF and Scheme Programming | 2 | January 10, 2016 03:13 |
Error messages | atg | enGrid | 7 | August 30, 2013 11:16 |
Phase locked average in run time | panara | OpenFOAM | 2 | February 20, 2008 14:37 |