[openfoamdjxych]

Originally Posted by Phicau

Hi Yuqian,
The Dev branch was added for OF7 / OF+20xx
Full forward compatibility will happen when I release the new version of the code, for which I still have no fixed date.

Best,
Pablo

Hi Pablo,
I have tried the Dev branch in my openfoamv2206,but when I allwmake,I get some mistakes,when I use olaflow without the Dev branch,allwmake can work well,can you give me some suggestions, thank you！

[goodstah]

Does anyone know or can explain a couple options in waveDict as well as for wave absorption BC:

1. rereadU + rereadAlpha: I've played around with this and am getting unexpected behavior with reread* in waveDict -- what would be the appropriate case for using this option set to "true"? If I have a simulation using wave theory "tveta" with a long list of time, velocity and free surface hts, and I stopped my sim at say 50 s, and I wanted to restart the sim where I left off, would I set "reread*" to true? or is it for another use?
2. nEdgeMin + nEdgeMax: I tried using these parameters to stabilize a sim, and played around with employing them at all four corners of the rectangular domain including the inlet BC (at both corners of the inlet) as well as the inlet side of the lateral BC's, and got some really unexpected behavior. My inlet wave maker had absorption turned on, so was thinking the "nEdge*" features might behave like with the outlet, but unfortunately not. My question is, are these parameters only intended for use with the "waveAbsorption2D(3D)Velocity" BC, and not so much as an absorption element included in the waveDict? Has anyone successfully employed the nEdgeMin/Max at all four corners of a domain?

Mahalo
Chris

[Phicau]

Hi pcwang,
As you mention, it is highly likely that there is something wrong with OpenFOAM in your workstation. I would recommend re-compiling OpenFOAM first and then olaFlow to see if you resolve it.

Hi openfoamdjxych,
Please note that OF+20xx refers to versions 2006 and 2012. Compatibility with newer versions is in the making with a lot of other code changes.

Hi Chris,
Happy to explain those.
1) No. Reread is only used when restarting the case AND changing the wave generation conditions. For example: you run a case with current, wait for it to achieve a steady state and stop it. Next, you add waves on top of the current in waveDict. When you restart you need to set both rereadU and rereadAlpha to true so that olaFlow knows that you have changed waveDict. Otherwise it will just assume that you continue the simulation with the previous conditions (which should be fully stored in the alpha.water and U files).
2) nEdgeMin and nEdgeMax tell the absorption boundary condition that water can only go out at that particular vertical transects. They are handy in cases in which absorption occurs at a sharp corner (see the example in my papers), in which if water is pumped into the domain, it may reach the top; or when using waveAbsorption3DVelocity theory at both sides of a corner (because this technique uses velocity measurements too). Please take a look at the example in the irreg45degTank tutorial. Note how it is only used at the corners of the 3D absorption to 3D absorption boundaries, and not at the wave generation boundary.

Best,
Pablo

[philiplu]

hello, ola-users and Phicau,
just FYI:
As I tried once, it looks like, the current "ola" works also fine under OF2112
if I'm wrong, sorry in advance.

best

philip

[Phicau]

The code in the main branch of olaFlow is compatible with the latest versions of OpenFOAM (10 and 2206). There has been quite some changes in the file naming conventions in version 10, so I will update the tutorials shortly.

The code in the Dev branch of olaFlow is only compatible up to versions 7 and 20xx. Moving forward this branch will be fully integrated in the main branch and updated to work with the latest OpenFOAM versions.

Pablo

[philiplu]

hello, Phicau,
many thanx for the info and "ola"
best
Philip

[Stan Zhou]

Hi all,

I wonder how to extract the time history of pressure at a certain point on a moving object (such as the floating object in the tutorial "waveFloatingObject") when post-processing a dynamic mesh case?



The OpenFOAM version I'm using is OpenFOAM 8.


Thanks in advance!

[goodstah]

This is sort of a Hail Mary at this point, but curios if anyone has insights or recommendations for models of longer run times, like 120s+? I have been running some wave models using the tveta wave boundary, which has been working great. But I'm running into an implacable issue with instabilities at the inlet at between 29 and 42 seconds, where the pressure blows up and alpha.water grows well beyond 1. I have experimented with using different numbers of nPaddles, nEdgeMin/Max, and even different wave generation theories such as cnoidal and stream function, and even lowering wave heights to extremely low heights such as 0.1 m (0.01% of domain ht), and yet the model still blows up after 42 s, and always at the inlet, either at the water surface of sea floor, with extreme negative pressure and alpha.water > 1, starting at one cell (not always the same cell either). Just wondering if anyone has any advice or recommendations or best management practices for long duration wave models?
Many thanks

[Biheng]

Hi all,

I am a fresh user of olaFlow. I have a problem with running olaFlow tutorial cases. I ran the cases after compiling the code following the tutorial with OpenFOAM-10. However, I get an error when I run the final olaFlow command. The error is following:

blockMesh meshing...
Preparing 0 folder...
Setting the fields...
Running...


--> FOAM FATAL ERROR:
fvMeshStitcher for changing mesh not found

Valid fvMeshStitchers are :
1(stationary)

From function static Foam::autoPtr<Foam::fvMeshStitcher> Foam::fvMeshStitcher::New(Foam::fvMesh&, bool)
in file fvMesh/fvMeshStitchers/fvMeshStitcher/fvMeshStitcherNew.C at line 86.

FOAM exiting

Simulation complete.

Could you guys give me some suggestions about how to fix that?
Many thanks.

[OussCous]

Hello Pablo,

I'm having the same problem as Callum, I really want to implement a power law current profile along with waves.
With simpleFoam/pimple solvers, i used codedFixedValue (or GroovyBC) with no problem but i can't figure out how to implement this profile with olaflow.
Based on the conversation you had with Callum, i see that i have to change the ***.h files. unfortunetly, the discussion ended in private.

Originally Posted by Phicau

Hi Callum,
The first implementation looks alright, however for the second you might need to implement it slightly different, with += instead of = to prevent the wave values being overwritten by the current.
We'll continue our communication via email.

.
.
.

Best,
Pablo

Could you please send me via email also the correct way to change the ***.h files in order to have a correct wave+current profile interaction.
my adress : oussema.mrzuk@gmail.com

Thank you very much in advance.
Kind regards.

[goodstah]

Originally Posted by goodstah

...running into an implacable issue with instabilities at the inlet at between 29 and 42 seconds, where the pressure blows up and alpha.water grows well beyond 1.

Problem Solved. If having similar issues with long runs, check out the use of "limitFields," but recommend to use judiciously...

[Stan Zhou]

Dear Pablo,
I was trying to compile the olaIsoFlow but some errors happened. The isoAdvector has been successfully compiled and the directory used during the compilation of olaIsoFlow is correct. The olaIsoFlow compilation log is attached. I'll be appreciated if you could help me solve this problem.

Best regards,
Stan

[lucia_995]

Hi, I try to reproduce irregular waves in olaFlow, but I don't know all wave profile but only Ts,Hs. How can i have waveHeights, wavePeriods, wavePhases to put in waveDict??

[goodstah]

Originally Posted by lucia_995

Hi, I try to reproduce irregular waves in olaFlow, but I don't know all wave profile but only Ts,Hs. How can i have waveHeights, wavePeriods, wavePhases to put in waveDict??

Hi Lucia, not sure if I understand completely, but it sounds like you just don't know what to plug in for wave phase (your H and T are good right)? If that's the case, you can pretty much use any value between 0 and 2pi, maybe with some form of random number generator if you have many bands to populate. If starting from scratch, you can just take your wave spectra and breakdown into a 9-band table (or other # of bands) of H/T/Dir, and then use Matlab or excel to generate random values for phase (0 < phase < 2pi) unless you know exactly where you want your wave crests to be in the time domain. Hope this helps.
Aloha
Chris

[lucia_995]

No, I want this file :

/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.3 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object waveDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

waveType irregular;

genAbs 1;

absDir 0.0;

nPaddles 1;

secondOrder 1; // 0/1, true/false, on/off

waveHeights
439
(
0.000001
0.000001
0.000002
0.000003
0.000005
0.000008
0.000011
0.000016
0.000022
0.000030
0.000041
0.000055
0.000071
0.000092
0.000117
0.000147
0.000183
0.000224
0.000272
0.000326
0.000388
0.000457
0.000534
0.000619
0.000711
0.000812
0.000921
0.001038
0.001163
0.001295
0.001434
0.001580
0.001732
0.001890
0.002054
0.002222
0.002394
0.002571
0.002751
0.002933
0.003118
0.003304
0.003491
0.003678
0.003865
0.004052
0.004238
0.004422
0.004605
0.004786
0.004966
0.005143
0.005318
0.005490
0.005661
0.005830
0.005999
0.006166
0.006333
0.006501
0.006670
0.006842
0.007018
0.007200
0.007388
0.007584
0.007789
0.008006
0.008237
0.008481
0.008741
0.009019
0.009314
0.009628
0.009962
0.010314
0.010683
0.011068
0.011465
0.011870
0.012280
0.012687
0.013086
0.013468
0.013827
0.014153
0.014437
0.014673
0.014854
0.014975
0.015033
0.015030
0.014984
0.014898
0.014773
0.014612
0.014418
0.014195
0.013946
0.013675
0.013387
0.013086
0.012776
0.012459
0.012141
0.011824
0.011511
0.011203
0.010904
0.010615
0.010337
0.010072
0.009819
0.009580
0.009354
0.009140
0.008940
0.008753
0.008578
0.008415
0.008262
0.008120
0.007989
0.007866
0.007751
0.007645
0.007545
0.007452
0.007364
0.007282
0.007204
0.007131
0.007061
0.006994
0.006931
0.006870
0.006811
0.006754
0.006699
0.006645
0.006592
0.006541
0.006490
0.006441
0.006392
0.006343
0.006296
0.006248
0.006201
0.006155
0.006109
0.006063
0.006018
0.005973
0.005928
0.005883
0.005839
0.005795
0.005751
0.005707
0.005664
0.005621
0.005579
0.005536
0.005494
0.005452
0.005411
0.005369
0.005328
0.005287
0.005247
0.005207
0.005167
0.005127
0.005088
0.005049
0.005010
0.004972
0.004934
0.004896
0.004858
0.004821
0.004784
0.004748
0.004711
0.004675
0.004640
0.004604
0.004569
0.004534
0.004500
0.004465
0.004431
0.004398
0.004364
0.004331
0.004298
0.004266
0.004233
0.004201
0.004170
0.004138
0.004107
0.004076
0.004045
0.004015
0.003985
0.003955
0.003926
0.003896
0.003867
0.003839
0.003810
0.003782
0.003754
0.003726
0.003699
0.003671
0.003644
0.003618
0.003591
0.003565
0.003539
0.003513
0.003487
0.003462
0.003437
0.003412
0.003387
0.003363
0.003339
0.003315
0.003291
0.003267
0.003244
0.003221
0.003198
0.003175
0.003153
0.003130
0.003108
0.003086
0.003065
0.003043
0.003022
0.003001
0.002980
0.002959
0.002938
0.002918
0.002898
0.002878
0.002858
0.002838
0.002819
0.002799
0.002780
0.002761
0.002742
0.002724
0.002705
0.002687
0.002669
0.002651
0.002633
0.002615
0.002598
0.002580
0.002563
0.002546
0.002529
0.002512
0.002496
0.002479
0.002463
0.002447
0.002431
0.002415
0.002399
0.002383
0.002368
0.002352
0.002337
0.002322
0.002307
0.002292
0.002278
0.002263
0.002249
0.002234
0.002220
0.002206
0.002192
0.002178
0.002164
0.002151
0.002137
0.002124
0.002110
0.002097
0.002084
0.002071
0.002058
0.002046
0.002033
0.002020
0.002008
0.001996
0.001983
0.001971
0.001959
0.001947
0.001935
0.001924
0.001912
0.001901
0.001889
0.001878
0.001867
0.001855
0.001844
0.001833
0.001822
0.001812
0.001801
0.001790
0.001780
0.001769
0.001759
0.001749
0.001738
0.001728
0.001718
0.001708
0.001698
0.001688
0.001679
0.001669
0.001659
0.001650
0.001640
0.001631
0.001622
0.001613
0.001603
0.001594
0.001585
0.001576
0.001567
0.001559
0.001550
0.001541
0.001533
0.001524
0.001515
0.001507
0.001499
0.001491
0.001482
0.001474
0.001466
0.001458
0.001450
0.001442
0.001434
0.001426
0.001419
0.001411
0.001403
0.001396
0.001388
0.001381
0.001373
0.001366
0.001359
0.001352
0.001344
0.001337
0.001330
0.001323
0.001316
0.001309
0.001302
0.001295
0.001289
0.001282
0.001275
0.001269
0.001262
0.001255
0.001249
0.001242
0.001236
0.001230
0.001223
0.001217
0.001211
0.001205
0.001198
0.001192
0.001186
0.001180
0.001174
0.001169
0.001163
0.001157
0.001151
0.001145
0.001140
0.001134
0.001128
0.001123
0.001117
0.001111
0.001106
0.001100
0.001095
0.001090
0.001084
0.001079
0.001074
0.001069
0.001063
0.001058
0.001053
0.001048
0.001043
0.001038
0.001033
0.001028
0.001023
0.001018
0.001013
0.001008
0.001004
0.000999
0.000994
0.000989
0.000985
);
wavePeriods
439
(
6.215700
6.207500
6.085600
6.022800
5.955000
5.911700
5.805300
5.743000
5.673700
5.623100
5.595600
5.546800
5.457600
5.401100
5.346800
5.294000
5.234700
5.180400
5.160500
5.101700
5.052300
5.012700
4.964600
4.919400
4.843600
4.819900
4.781000
4.741800
4.660700
4.651800
4.600000
4.556600
4.531400
4.492400
4.453200
4.416200
4.384300
4.347200
4.309200
4.247000
4.238000
4.200400
4.164100
4.112100
4.092500
4.071100
4.026100
3.996400
3.953500
3.932900
3.905500
3.871700
3.841600
3.835200
3.798400
3.777200
3.751300
3.702700
3.697000
3.659300
3.623300
3.601500
3.576500
3.554600
3.525100
3.506800
3.481900
3.462100
3.432600
3.425400
3.403100
3.377400
3.347200
3.322900
3.305600
3.300900
3.265700
3.260100
3.224900
3.217000
3.186700
3.163600
3.158600
3.138100
3.107600
3.094000
3.081900
3.059600
3.047800
3.024200
3.004100
2.992900
2.965800
2.962300
2.944300
2.924800
2.909300
2.882200
2.881000
2.850600
2.848400
2.827800
2.811700
2.799900
2.777700
2.772400
2.745500
2.732300
2.719900
2.711400
2.693700
2.688500
2.668500
2.648600
2.638400
2.630800
2.611500
2.604400
2.587800
2.573400
2.564900
2.551400
2.540000
2.522700
2.508500
2.501700
2.484600
2.482700
2.470400
2.457000
2.446300
2.426200
2.424300
2.411800
2.394900
2.392000
2.377200
2.370300
2.353700
2.342500
2.333800
2.319300
2.312000
2.303100
2.291700
2.282200
2.273400
2.260700
2.257400
2.239800
2.236300
2.220000
2.216000
2.209000
2.192300
2.189400
2.173600
2.168900
2.160200
2.151300
2.140700
2.131000
2.127900
2.114400
2.108500
2.099300
2.093900
2.080100
2.070500
2.063900
2.061300
2.049900
2.044800
2.034900
2.026700
2.017500
2.009100
2.004200
1.996600
1.983600
1.976300
1.968300
1.961600
1.954600
1.946700
1.940700
1.938800
1.924800
1.923400
1.917000
1.905700
1.902600
1.893600
1.883700
1.882800
1.875800
1.869200
1.856200
1.850400
1.844800
1.836500
1.832900
1.825300
1.820900
1.816200
1.808600
1.799800
1.796500
1.790500
1.785900
1.776500
1.772500
1.762100
1.756500
1.755500
1.747400
1.743600
1.734500
1.730000
1.721800
1.718400
1.715500
1.706600
1.703100
1.695200
1.693100
1.685500
1.679900
1.674800
1.668300
1.661400
1.657500
1.653500
1.646900
1.642900
1.637400
1.631300
1.629100
1.624200
1.616000
1.609900
1.605400
1.601400
1.595600
1.593100
1.588900
1.583400
1.580000
1.572100
1.569200
1.563100
1.557200
1.553500
1.548700
1.546700
1.542200
1.535900
1.532100
1.525900
1.522200
1.519200
1.515800
1.511400
1.503700
1.499600
1.495500
1.491900
1.488100
1.481800
1.479500
1.477200
1.471900
1.468400
1.462600
1.458000
1.454600
1.452500
1.444700
1.441300
1.437500
1.436200
1.430400
1.425700
1.422900
1.417400
1.413500
1.412600
1.406200
1.402000
1.399200
1.397500
1.393700
1.387100
1.384200
1.381800
1.378800
1.375300
1.370200
1.366500
1.363800
1.360400
1.355800
1.353300
1.350700
1.344800
1.344100
1.337300
1.336500
1.331800
1.327400
1.323700
1.322800
1.317100
1.315700
1.311800
1.310100
1.306200
1.300900
1.297300
1.295200
1.291700
1.288700
1.285600
1.283000
1.279500
1.275400
1.274100
1.269900
1.267200
1.265000
1.260800
1.257900
1.255100
1.251000
1.247700
1.245200
1.241800
1.239900
1.236400
1.235600
1.231800
1.228500
1.226000
1.223700
1.218900
1.217600
1.215500
1.212300
1.209400
1.207100
1.202700
1.200600
1.196500
1.195700
1.190700
1.189500
1.186700
1.182900
1.181700
1.177400
1.176100
1.173600
1.169600
1.166900
1.164900
1.163600
1.159700
1.158600
1.154300
1.152600
1.150600
1.146800
1.146100
1.142400
1.139800
1.136700
1.134400
1.131400
1.130900
1.127700
1.124200
1.122100
1.119200
1.118200
1.115900
1.112400
1.111500
1.109200
1.106600
1.102700
1.100500
1.099000
1.097100
1.094000
1.092000
1.088700
1.087000
1.084800
1.083100
1.081400
1.077800
1.077500
1.075100
1.071500
1.069100
1.066800
1.065400
1.062400
1.061400
1.057900
1.057600
1.054100
1.051600
1.049900
1.047600
1.046300
1.044700
1.042800
1.039700
1.037200
1.036300
1.034600
1.031100
1.029100
1.028200
1.025800
1.022600
1.021800
1.020200
1.018500
1.014700
1.013200
1.011300
1.009300
1.007400
1.006100
1.004000
1.002100
0.999270
);
wavePhases
439
(
6.031400
0.448110
4.733100
4.063200
3.980900
1.721500
5.597100
5.439500
6.122800
5.028700
1.612400
0.477510
3.855700
3.846200
3.726400
3.558500
4.287000
4.558200
0.712850
1.750500
1.738700
0.577700
0.632330
0.419590
4.546600
1.599900
0.838260
0.225270
5.969000
1.028400
2.669700
3.167900
0.889190
0.889450
1.014900
0.892000
0.018819
0.102550
0.471020
5.315500
0.669370
1.316700
1.831400
5.503400
3.029900
0.962790
3.674100
3.421600
6.053300
4.097600
3.665200
4.696100
5.042300
0.197510
2.192400
0.732790
0.450080
5.616100
0.708680
3.558300
6.163900
5.355700
5.443400
4.793300
6.229300
4.799400
5.175700
4.276300
6.136000
1.860900
1.833900
2.879300
5.351700
6.241900
5.112600
0.233210
4.638800
0.084133
4.781800
0.977080
4.385500
5.597900
0.968030
1.448400
5.430600
3.778500
1.649600
3.107400
0.990880
3.069200
4.046200
1.852400
5.574100
0.587980
1.044600
2.194000
1.895500
6.186200
0.369800
6.236400
0.827080
2.902700
3.271900
1.881900
4.989800
0.937820
6.255100
5.708000
4.882600
2.415800
4.060100
0.127070
2.992200
5.895200
4.417300
1.703700
4.636400
1.787100
3.579900
4.356000
2.292800
2.779500
2.228800
4.804200
5.876700
3.172300
5.873000
0.605880
0.903330
1.863200
1.384500
6.209600
0.973740
1.717900
4.998600
0.400700
2.606000
0.351300
3.830800
4.167900
3.107100
5.624200
3.766500
2.999200
3.767500
3.458900
2.695600
4.518100
0.337510
5.493900
1.503400
5.948000
2.376900
0.779420
5.763600
1.449600
6.083700
3.065600
2.803900
2.853000
4.066500
4.774200
0.699350
4.291400
2.320500
2.875100
0.555910
4.665700
5.731800
4.484000
0.189900
2.780400
0.482500
2.003700
2.175200
3.308500
3.821400
1.497800
1.446800
5.876200
5.752100
6.207400
5.594800
5.248900
5.785200
4.755000
0.128760
6.139800
1.050000
0.422220
4.314200
0.863630
2.735400
5.499300
0.115410
0.330230
0.196630
6.162500
5.434700
4.486600
6.263300
3.465800
4.632600
2.697000
1.051300
2.400600
4.920500
2.004400
1.837200
0.261430
3.666500
1.548500
6.207300
5.781200
0.644490
3.023600
0.772700
4.315500
3.010300
5.705400
3.170400
0.109400
3.857500
1.517300
4.162700
0.235650
2.704400
2.913500
2.584300
3.928300
5.804400
4.129500
2.622700
4.240500
2.909100
3.242400
4.480100
0.952260
0.629100
4.671400
5.986300
5.426100
4.297000
5.377400
2.300900
1.549000
2.405700
0.610010
4.643400
2.291200
4.133700
5.749000
4.475000
4.817300
1.271700
1.126800
3.548200
2.621300
4.981000
3.938700
1.845000
0.385440
0.376060
5.288800
4.931100
4.657800
3.651700
2.909100
6.019100
3.187600
0.304740
2.089000
1.071900
3.618200
4.392000
3.318700
0.319570
6.257400
5.257600
5.087000
0.770860
3.778700
5.105200
3.264900
5.994400
6.004300
1.284600
5.512200
6.277500
4.547400
1.184300
1.292800
6.144300
4.882700
2.724800
1.517700
1.272700
3.857800
4.044700
2.638500
2.363000
4.197900
2.307100
0.727230
5.007200
0.026421
6.144800
1.406000
3.720700
5.649900
6.227200
1.626200
6.026000
2.407600
3.490400
0.540870
1.681300
5.541400
6.244100
4.143100
4.630600
4.241700
4.222000
3.074900
3.694900
5.777200
2.100900
4.196500
3.498900
1.768900
4.074800
3.822600
3.256400
5.549200
6.144500
5.109700
6.181400
3.840500
5.064000
0.516210
2.309400
3.342000
2.458200
1.099400
5.460400
2.069800
0.322110
1.167000
1.376000
0.276440
3.827400
2.458500
5.496900
0.918240
6.257600
2.674000
3.034300
5.424200
2.152300
5.792700
2.720600
2.319000
5.534600
5.768300
4.256500
1.085100
4.357900
0.747840
4.938800
2.853700
1.639400
4.752700
0.261380
3.130800
3.438900
4.929100
4.501700
5.996900
0.800730
2.949200
5.605300
4.722100
6.144900
2.364700
2.136100
5.125500
1.270600
1.000700
1.823800
5.804100
5.411900
3.285700
2.217000
4.394200
3.564000
6.263900
4.947100
4.531100
3.018000
1.538800
5.455600
0.085704
0.488350
4.396900
4.994100
5.250200
3.138800
5.608700
1.992500
6.067300
0.713250
4.563500
5.772800
4.359900
4.955000
2.669500
1.241400
0.556330
3.856400
5.122600
1.628100
0.496190
4.874000
4.825700
1.180500
2.304600
6.193300
2.296700
1.249900
0.314150
5.936600
4.482000
4.280600
4.209200
4.161600
2.180200
2.563700
2.506400
5.539000
);
waveDirs
439
{
0
};

// ************************************************** *********************** //




but I don't know how to produce it.
I know only significant waveheight and peak period of my wave.
Thank you very much.

Lucia.

[Jeggi]

Is there any tutorial on to use olaFlow with adaptive mesh refinement? I've tried setting it up myself, but I keep getting errors after the first iteration.
I've attached my current attempt at setting it up.

[Mahdi_Kh]

Hi everyone,
I have been running the compilation of Olaflow for like 16 hours on my laptop, and I am not sure it should take such time, or I have made a mistake during the process. Can you tell me how long the usual time is needed for compilation.
Thank you.

[AUbuntu]

Hi all,

Could anyone let me know if a fully developed current flow profile (e.g. log-low equation) can be combined with wave in olaFlow at the inlet?

In the waveDict as much as I see, we can only add the mean value for uCurrent!

if yes, could you please let me know how!

Kind regards,
Amir

[mdelisle]

Hello all!

I'm currently trying to simulate specific waves from a field experiment using the wavemaker (tveta) in olaFlow. However, I'm noticing discrepancies between the input eta and velocity and the modelled eta and velocity at the offshore boundary. The offshore boundary is at ~8m depth which stretches Pablo's recommendation of using the wavemaker for shallow water conditions, however I don't believe that this would cause the discrepancy in eta and velocity at the offshore boundary that I am seeing. The initial water surface level is set to the first value in the waveDict.

Below are plots showing the eta and velocity results for several model runs at the offshore boundary of my domain. The different colored lines represent different model runs where I have linearly subtracted a specific value from the grey line (experimental data I want to replicate). Model run 1 (bright green line) uses the grey line (experimental) input directly (for both eta and velocity), but as you can tell, produces significantly higher eta and velocity. Model run 7 (red line) subtracts 0.91 m from the grey line as the eta input (no change to velocity input from the grey line) and reproduces the grey line well, though still overpredicts the velocity. I've also included a snapshot of the model domain below for reference.

gauges.png

domain.png

Do you have any idea where I might be going wrong in using the wavemaker to force the model? I've attached my waveDict and 0/U files (for model run 1; input is the grey lines for both eta and velocity).

Thanks in advance for any suggestions!

U FILE:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.x |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type waveVelocity;
waveDictName waveDict;
value uniform (0 0 0);
/*
type fixedValue;
value uniform (0 0 0);
*/
}
atmosphere
{
type pressureInletOutletVelocity;
value uniform (0 0 0);
}
bottom
{
type fixedValue;
value uniform (0 0 0);
}
slope
{
type fixedValue;
value uniform (0 0 0);
}
defaultFaces
{
type empty;
}
outlet
{
type waveAbsorption2DVelocity;
nPaddles 1;
nEdgeMin 0;
nEdgeMax 0;
absorptionDir 400;
value uniform (0 0 0);
}

// ************************************************** *********************** //

WAVEDICT
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.3.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object waveDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

waveType wavemaker;

waveTheory tveta;

genAbs 1;

absDir 0.0;

nPaddles 1;

timeSeries
300001
(
0
0.001
0.002
0.003
0.004
0.005
EXCLUDING 0.006 - 299.997 FOR BREVITY
299.998
299.999
300
);

paddleEta
300001
(
2.3803
2.3804
2.3806
2.3807
2.3808
...
1.7956
1.7955
1.7955
1.7954
);

paddleVelocity
300001
(
0.178
0.1782
0.1784
0.1786
0.1788
...
-0.1933
-0.1933
-0.1933
-0.1932
);

[abas.rahmani86]

Dear olaFlower,

I'm modelling a wave flume with an intermediate depth condition, using the second-order Stokes wave theory. I was checking the maximum horizontal velocity that faced a local peak near the bottom. Firstly, I used boundary layer mesh and eliminated the WallFunctions, but the problem was not solved. I checked the velocity profile next to the inlet BC (Figure 1) and far from BC (Figure 2), and the problem is in both of them. I also faced the same difficulty with shallow water depth conditions under the Cnoidal wave theory. One can conclude that the problem is due to wave generation BC. If you have a solution, please let me know.

Best regards
Ali
Figure 1:
https://drive.google.com/file/d/1G5L...ew?usp=sharing
Figure 2:
https://drive.google.com/file/d/1eSL...ew?usp=sharing