[VLima]

Thank you!

Quick question:

In the breakwater case, is it ok that the file GaugesP_p.xy is not present in the folder postProcessing/sets/0/ ? A warning is issued concerning the file.

[Phicau]

Yes, it is totally normal, p does not exist at t = 0 s, it is created on the first step.

[VLima]

Dear Pablo,

For some reason when running runCaseKEpsilon (or SST) I am now retrieving the following errors :

Processing and plotting variables
File "postSensVOF.py", line 42
print 'Sensor ' + '%i' % int(i+1) + ' of ' + '%i' % nSens + '.'
^
SyntaxError: invalid syntax
File "postSensPres.py", line 40
print 'Sensor ' + '%i' % int(i+1) + ' of ' + '%i' % nSens + '.'
^
SyntaxError: invalid syntax
Traceback (most recent call last):
File "plotSensVOF.py", line 8, in <module>
a = os.listdir(readPath)
FileNotFoundError: [Errno 2] No such file or directory: '/home/user/OpenFOAM/user-v1606+/run/tutorials/olaFlow/tutorials/breakwater/gaugesVOF'
Traceback (most recent call last):
File "plotSensPres.py", line 8, in <module>
a = os.listdir(readPath)
FileNotFoundError: [Errno 2] No such file or directory: '/home/user/OpenFOAM/user-v1606+/run/tutorials/olaFlow/tutorials/breakwater/sensorsPres


I did not perform any changes in any files.

Any idea why this error now and how to solve this?

Thank you very much in advance.

[Phicau]

Hi Vania,

I think you are using Python 3. The postprocessing files are prepared to be run with Python 2.

Please run them manually with Python 2 or change the python executable in lines 67-70 of the runCase... script.

Best,

Pablo

[VLima]

Dear Pablo,

I adapted the wavemaker input, as in wavemakerFlume with Piston, for the baseWaveFlume case.

I sampled several gauges along the flume, unfortunately there is a problem in the sampling and I get an error.

The log for the sampling :


Create time

Create mesh for time = 0

Reading set description:
GaugeVOF01
GaugeVOF02
GaugeVOF03
GaugeVOF04
GaugeVOF05

Time = 0

Time = 0.05


The terminal error:

...
Simulation complete.
Sampling free surface...
#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
#1 Foam::sigSegv::sigHandler(int) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so"
#2 ? in "/lib64/libc.so.6"
#3 Foam::HashTable<Foam::regIOobject*, Foam::word, Foam::string::hash>::find(Foam::word const&) const in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so"
#4 Foam::MeshObject<Foam:ointMesh, Foam::UpdateableMeshObject, Foam:ointConstraints>::New(Foam:ointMesh const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so"
#5 void Foam::volPointInterpolation::interpolate<double>(F oam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam:ointPatchField, Foam:ointMesh>&) const in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so"
#6 Foam::tmp<Foam::GeometricField<double, Foam:ointPatchField, Foam:ointMesh> > Foam::volPointInterpolation::interpolate<double>(F oam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&, bool) const in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so"
#7 Foam::interpolationCellPoint<double>::interpolatio nCellPoint(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so"
#8 Foam::interpolation<double>::adddictionaryConstruc torToTable<Foam::interpolationCellPoint<double> >::New(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so"
#9 Foam::interpolation<double>::New(Foam::word const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so"
#10 Foam::sampledSets::volFieldSampler<double>::volFie ldSampler(Foam::word const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::PtrList<Foam::sampledSet> const&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libsampling.so"
#11 void Foam::sampledSets::sampleAndWrite<double>(Foam::sa mpledSets::fieldGroup<double>&) in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libsampling.so"
#12 Foam::sampledSets::write() in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/lib/libsampling.so"
#13 ? in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/sample"
#14 __libc_start_main in "/lib64/libc.so.6"
#15 ? in "/opt/OpenFOAM/OpenFOAM-v1606+/platforms/linux64GccDPInt32Opt/bin/sample"
runCase: line 63: 784 Segmentation fault (core dumped) sample -dict system/sampleDictVOF > sampleVOF.log


Any idea what might be wrong in my settings?

Thank you very much in advance for your help.

[Phicau]

Hi Vania,

did I understand correctly, your mesh is moving, right?

If I am not wrong the sampling utility fails because your mesh is moving and the utility does not update it. Therefore, the initial points in which you were sampling are possibly in new cells right now and this triggers an error.

You might solve this by sampling at runtime with a function object or by creating your own script in which you need to do:

sample -dict system/sampleDictVOF -time $timeFolderName

for all the times (timeFolderName) of your case.

Best,

Pablo

[Phicau]

Dear all,

I have just included olaFlow in the Zenodo database to obtain a DOI. This identification system works in a similar way as the DOI of papers.

The DOI identifier for olaFlow is: https://doi.org/10.5281/zenodo.1297012

Important notice

Although no significant functionality changes have been introduced in the code during the preparation of the code for the DOI minting, the git history has been rewritten. This means that your local copy of olaFlow will not be able to obtain the latest updates in the usual manner.

Users are advised to download a fresh copy of the repository by running the following instruction in the terminal:

git clone git://github.com/phicau/olaFlow.git

Alternatively, advanced users also correct this by running:

git pull
git reset --hard origin/master
git pull origin master

After this, future updates will be available again.

Users who have made custom changes to the code and wish to continue getting updates may integrate them in the usual simple (e.g. copying and pasting their files) and advanced (e.g. git rebase) ways.

This is the official release announcement: https://sites.google.com/view/olaflo...og/olaflow-doi

[VLima]

Dear Pablo,

I am trying to visualize the vertex numbers/indices of the wavemakerTank case (tried with Flume case too).

Doing paraFoam -block I get :

bash-4.1$ paraFoam -block
Created temporary 'wavemakerTank_teste.blockMesh'
fileName::stripInvalid() called for invalid fileName /home/user/Desktop/LinktoOpenFOAM/user-v1606+/run/tutorials/olaFlow/tutorials/wavemakerTank_teste/wavemakerTank_teste.blockMesh
For debug level (= 2) > 1 this is considered fatal


It is supposed to be possible to visualize with this option.

Any idea what might be causing this?

I would appreciate help or any other feasible suggestions to do this.

Many thanks in advance!

[Phicau]

Hi Vania,

no idea, I have never used that option. However, this is a problem with paraFoam rather than with olaFlow, so I suggest that you open a thread here: https://www.cfd-online.com/Forums/openfoam-paraview/

Just a suggestion, try to move your case to the /home/user/ folder and try again.

Best,

Pablo

[VLima]

Well that surprisingly worked!

Thank you!

[dietlind]

Dear Pablo,


I am new to olaFlow and OpenFoam and just started running a few of the tutorials.
I have some issues when it comes to the wavemakerFlume tutorial (Piston wavemaker). I did not change anything in the files provided. I only added the sampling method from the breakwater tutorial to export the surface elevation.

I am able to run the simulation and there is no error occuring.


When I look at the results in Paraview and also in the surface elevation file, there is no wave generated.
The wavemakermovement.txt file is created and shows a changing in paddle position but still there is no change in surface elevation.


Do you have any idea what the problem might be?




Furthermore, when I want to adapt the data and use experimental data (wavemaker position), can I just change the data in the wavemakermovement.txt file and insert my data andthen run the case again or is there anything else that needs to be adjusted?



Kind regards,
Dietlind

[Phicau]

Hi Dietlind,

which OpenFOAM version are you using? Are you sure that you are running the case with the olaDyMFlow solver?

If the case without changes works then you know that there is something wrong with the way you introduce the sampling procedure. Otherwise please come back to me with more details, there might be a bug somewhere.

Yes, you just need to edit the wavemakermovement.txt file and depending on the duration of your own signal, controlDict to set the endTime.

Best,

Pablo

[dietlind]

Hi Pablo,


thanks for the quick reply. I am using Open Foam Version 2.4. I used the olaDyMFlow command to run the case.

Actually there was one minor change that I did to the transportProperties file because I always got an error saying "wrong token, expected word, found '['" in the line where sigma is specified. So I changed this line to


sigma sigma [1 0 -2 0 0 0 0] 0.07;


Then it worked. It was the same for the other tutorials I conducted. I run the baseWaveFlume and breakwater tutorial and there I also had to make this change in sigma but these two tutorials worked.

Furthermore, for the sampling I implemented the breakwater sampling to the baseWaveFlume tutorial to check if I had done it correctly and it worked.



In paraview also waves appeared for the baseWaveFlume and breakwater tutorial but not for the piston wave maker. The sampling was just to check if it was more an illustrative problem with paraview but that seems not to be the case.



I am not sure which information you need exactly to find out what the solution to that problem might be.


I attached the wavemakermovement.txt file which was created during the run aswell as the sampleDictVOF file I used. As an example of what I got from the sampling I added the file for gauge04.

I hope you can help me to solve this.
Thanks in advance!


wavemakerMovement.txt

sampleDictVOF.txt

GaugeVOF04.txt

Originally Posted by Phicau

Hi Dietlind,

which OpenFOAM version are you using? Are you sure that you are running the case with the olaDyMFlow solver?

If the case without changes works then you know that there is something wrong with the way you introduce the sampling procedure. Otherwise please come back to me with more details, there might be a bug somewhere.

Yes, you just need to edit the wavemakermovement.txt file and depending on the duration of your own signal, controlDict to set the endTime.


Best,

Pablo

[Phicau]

Hi Dietlind,

the olaFlow tutorials for version 2.4 and any older versions are in another repository, try those:

https://github.com/phicau/olaFlow_oldVersionTutorials

Pablo

[arashghgood]

Dear Pablo

is it possible to plot the following in paraveiw after running olaflow examples?

the position of the crest of the wave
the position of the crest of the interacted waves
the changes of the width and so on

[Phicau]

Hi Arash,

although I have never done that there is probably a way to do it, preview can do pretty complex analysis with the elevation, threshold and python programming filters.

I would suggest, however, that you export free surface elevation and perform such analysis in python/matlab, because it might be easier and more universal to apply to other of your cases.

Best,

Pablo

[arashghgood]

Dear Pablo
How can I make two solitary waves with the same direction but with different properties such as amplitude and velocity so that they can interact with each other?

[arashghgood]

Dear Pablo
is it possible to write my own solver based on my 1D or 2D equations in OpenFOAM or olaFlow to study a given problem?
In fact, I have an equation which is KdV type in 1D. I would like to solve and study my problem based on this equation.

[Phicau]

Hi Arash,

there is a way to simulate this using olaFlow, solving the RANS equations. The process is as follows:
- Simulate the first solitary wave.
- See the log and stop the simulation at time equal to the reported equivalent wave period of the solitary wave.
- Duplicate the original simulation, remove all the time steps except the last one and rename it to be time 0.
- Change the solitary wave parameters and set rereadAlpha and rereadU to true in waveDict


Regarding solving KdV equation, yes, it can be implemented in a custom solver. You will also need to modify the olaFlow wave generation BCs to work with the new solver.

Good luck

Pablo

[arashghgood]

Dear Pablo
Regarding solving KdV equation, which files and which parts of the files should be modified or changed?