|
[Sponsors] |
March 13, 2018, 21:01 |
|
#21 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Hi Chen,
olaFlow is written and tested to work in parallel, if you experience any differences in the unmodified version, please report it, as it would most probably be a bug. However, I have never experience such. Probably there are some glitches in your new implementation that might need to be programmed to be parallel-aware. Best, Pablo |
|
April 16, 2018, 09:15 |
|
#22 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Dear all,
I have just released a library that provides modified versions of k-ε and k-ω SST turbulence models to simulate correctly multiphase systems and mitigate turbulence build-up effect. You can find all the information in the release post: https://sites.google.com/view/olaflo...ve-simulations I hope that you find it useful. Best, Pablo |
|
April 16, 2018, 13:28 |
|
#23 | |
Senior Member
Hua Zen
Join Date: Mar 2009
Posts: 138
Rep Power: 17 |
Quote:
|
||
April 16, 2018, 21:32 |
|
#24 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Hi Hua Zen,
Check this reference out: Application of a buoyancy-modified k-ω SST turbulence model to simulate wave run-up around a monopile subjected to regular waves using OpenFOAM Brecht Devolder, Pieter Rauwoens & Peter Troch Coastal Engineering (2017), vol. 125, pp. 81–94 https://doi.org/10.1016/j.coastaleng.2017.04.004 Best, Pablo |
|
April 17, 2018, 12:17 |
|
#25 |
Senior Member
Hua Zen
Join Date: Mar 2009
Posts: 138
Rep Power: 17 |
Hi Phicau, I've read the paper by Devolder et al. (2017), in which an additional buoyancy term is used to suppress the spurious turbulence generation at the interface. From your blog, two additional turbulence models are implemented(kEpsilonMultiphase, kOmegaSSTMultiphase), what is the difference between these two implementation and that of Devolder et al. (2017)?
|
|
April 17, 2018, 23:37 |
|
#26 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Hi Hua Zen,
the problem in OpenFOAM is that incompressible turbulence models are not prepared for multiphase flows. They assume a constant rho, since under the incompressibility assumption, it will go out of the differential operators. In incompressible multiphase flows, the density of each phase is constant, but the density of the system changes throughout the domain depending on the VOF function value in each cell. kEpsilonMultiphase and kOmegaSSTMultiphase implementations are exactly the same as provided by OpenFOAM, but the density (rho) has been modified so that the one provided by the multiphase solver (e.g. olaFlow) is used. Best, Pablo |
|
April 25, 2018, 05:07 |
|
#27 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Dear all,
olaFlow is now compatible with OpenFOAM-dev. I have made some changes in the compilation scripts and added the olaFlow solver. Please note that olaFlow and olaDyMFlow have been merged in the dev version, so olaFlow can handle moving meshes too. As usual, the latest version is in: https://github.com/phicau/olaFlow Best, Pablo |
|
April 29, 2018, 03:34 |
|
#28 |
New Member
aref shahryari
Join Date: Jul 2017
Posts: 10
Rep Power: 9 |
Dear Pablo,
Hi I am wondering if it is possible to change the waveTheory in the wavemakerFlume tutorial (since it doesn't have any waveDicts)? If yes, how? |
|
April 29, 2018, 22:02 |
|
#29 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Hi Aref,
the waveDict is generated by the runCase script, in which the linear wave motion of the wavemaker (piston and flap) is programmed. Once you have run the case you can edit the waveDict file as you desire, but runCase only provides the linear wave theory transfer functions. Best, Pablo |
|
April 30, 2018, 05:39 |
|
#30 |
New Member
aref shahryari
Join Date: Jul 2017
Posts: 10
Rep Power: 9 |
Hi Pablo,
you mean the waveDict file can be edited while running the case? thanks for your reply. |
|
April 30, 2018, 05:42 |
|
#31 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Hi Aref,
no, what I mean is that the waveDict file is generated by the runCase script. Once it has been created you can modify it as you wish and use it in any case. However, the file is only read by olaFlow during the initial time step, therefore, any modifications that you make while the case is running will not have any effects on the simulation that is running. Best, Pablo |
|
April 30, 2018, 06:18 |
|
#32 |
New Member
aref shahryari
Join Date: Jul 2017
Posts: 10
Rep Power: 9 |
thanks again Pablo,
it was really helpful. |
|
May 17, 2018, 06:51 |
Reference olaFlow
|
#33 |
New Member
Diogo R.C.B. Neves
Join Date: Aug 2014
Posts: 16
Rep Power: 12 |
Hi Pablo, how do I properly reference olaFlow?
I enhance your novelties regarding the turbulence models. Beste Regards, Diogo |
|
May 17, 2018, 07:38 |
|
#34 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Hi Diogo,
it might be too hidden in the website: https://sites.google.com/view/olaflo...del/references "If you want to reference the model in your publications it should be called olaFlow, and include any of the following references when citing implementation, validation and applications." We have used the following example in a recent paper, also adding the GitHub address ( https://github.com/phicau/olaFlow ): The CFD model used in this work is olaFlow (Higuera et al., 2013a,b), developed within the OpenFOAM framework. We wanted to cite the wave generation and absorption capabilities plus the validation (Higuera et al., 2013a,b), but feel free to change the references. For example, if you would like to cite the porous model and applications it would be (Higuera et al., 2014a,b) or my thesis. For the moving-boundary wavemaker (Higuera et al., 2015)... I think you get the idea. Don't forget to give credit to Devolder et al. papers for that turbulence models! Best, Pablo |
|
May 24, 2018, 10:08 |
|
#35 |
New Member
Vania
Join Date: Dec 2016
Posts: 13
Rep Power: 9 |
Dear all,
I am new to OF and now testing the wavemakerFlume tutorial case, for a piston wavemaker, with olaFlow. I was wondering if : 1) it is possible (and how) to generate just only one single wave with the piston wavemaker , and not several continuous equal waves, and how to check the profile of the wave generated at some point along the channel; 2) in the case I am running, I got the error: --> FOAM FATAL ERROR: The time series is not long enough. and was wondering why could this be happening. I just changed the wave height and water depth values. I would appreciate very much any help you could provide with this. Thank you. |
|
May 28, 2018, 05:11 |
|
#36 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Hi Vania,
Yes, it is possible, you just need to edit wavemakerMovement.txt so that the wavemaker moves as you wish: one wave and then keep the position. You can sample free surface elevation, just check how it is done in the breakwater tutorial. That error indicates that the time series of movement that you provided is no long enough, your simulation time is longer. Select a shorter simulation time or give a longer time series. Best, Pablo |
|
June 1, 2018, 09:18 |
|
#37 |
New Member
Vania
Join Date: Dec 2016
Posts: 13
Rep Power: 9 |
Hello Pablo.
Thank you so much for your clarifications. I did not manage to check the plots for the breakwater case study yet. When running the case (SST and KEpsilon) I come across some errors in the Python scripts, namely concerning names and paths for files/folders I assume: bash runCaseKEpsilon blockMesh meshing... snappyHexMesh meshing... Preparing 0 folder... Setting the fields... Running... Simulation complete runCaseKEpsilon: line 44: [: OF_FLAVOUR: integer expression expected Sampling free surface... runCaseKEpsilon: line 59: postProcess: command not found Sampling pressure... runCaseKEpsilon: line 62: postProcess: command not found Processing and plotting variables Traceback (most recent call last): File "postSensVOF.py", line 21, in <module> a = os.listdir('./'+postPath) OSError: [Errno 2] No such file or directory: './postProcessing/sampleDictVOF' Traceback (most recent call last): File "postSensPres.py", line 21, in <module> a = os.listdir('./'+postPath) OSError: [Errno 2] No such file or directory: './postProcessing/sampleDictPres' The files sampleDictVOF and sampleDictPres are in the system folder. I tried to change the path but it does not work either. The errors are common to both SST and KEpsilon. Is anything missing in my installation? Also, is line 44 error relevant for this? Any suggestions to overcome this are very welcome. Thank you very much in advance. |
|
June 2, 2018, 05:35 |
|
#38 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
You found a glitch in the code that select the postprocessing tool for your version of OpenFOAM. I should be able to fix it quite easily. Which version are you using?
Please report what the terminal outputs for the following command: echo $WM_PROJECT_VERSION Pablo |
|
June 2, 2018, 07:30 |
|
#39 |
New Member
Vania
Join Date: Dec 2016
Posts: 13
Rep Power: 9 |
The command "echo $WM_PROJECT_VERSION" outputs "v1606+".
|
|
June 3, 2018, 22:13 |
|
#40 |
Senior Member
Pablo Higuera
Join Date: Jan 2011
Location: Auckland
Posts: 627
Rep Power: 19 |
Hi Vania,
all solved, there was a $ missing in line 44 of the script. All is working now. You can update the repository to the latest version following the update guide in https://sites.google.com/view/olaflowcfd/source-code Best, Pablo |
|
Tags |
olaflow, waves |
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Divergence detected in AMG solver: k when udf loaded | google9002 | Fluent UDF and Scheme Programming | 3 | November 8, 2019 00:34 |
udf problem | jane | Fluent UDF and Scheme Programming | 37 | February 20, 2018 05:17 |
UDF velocity profile | willroca | Fluent UDF and Scheme Programming | 2 | January 10, 2016 04:13 |
Error messages | atg | enGrid | 7 | August 30, 2013 12:16 |
Phase locked average in run time | panara | OpenFOAM | 2 | February 20, 2008 15:37 |