[Phicau]

Hi Hao,

take a look here:

http://www.cfd-online.com/Forums/ope...servation.html

Best,

Pablo

[dneves]

Hi Pablo,

I wonder if you, or your colleagues, use any special turbulence model (or turbulence parameters) por wave breaking and propagation on a beach profile?

I'm trying to simulate some waves breaking on a beach profile using Keplison, Komega and Komega-sst but it is killing the wave...i'm trying to figure out quickly the parameters for this type of simulations.

Thank you for your attention,

Best regards,

Diogo

[Phicau]

Hi Diogo,

not really, we use the regular models and parameters by default. Take a look at Jacobsen's waves2foam paper, they did some changes to the k-omega model which may interest you.

As for the wave height numerical dissipation, it may be caused by a not-so-good mesh quality that makes turbulence grow in excess.

Best,

Pablo

[dneves]

Thank you very much for your quick response

I will follow your suggestions,

Best regards!

Diogo

[totalart]

Hello,

I encounter errors when I try to install IHFOAM on my Mac.
Do you know the reason?
Thank you in advance.




wmakeLnInclude: linking include files to ./lnInclude
Making dependency list for source file IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.C
Making dependency list for source file IH_Waves_InletVelocity/IH_Waves_InletVelocityFvPatchVectorField.C
SOURCE=IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.C ; clang++ -m64 -DdarwinIntel64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -Wno-overloaded-virtual -Wno-unused-comparison -Wno-deprecated-register -O3 -DNoRepository -ftemplate-depth-100 -DOFVERSION=230 -I/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude -I./IH_Waves_InletVelocity/velProfiles -I../common -I../common/checks -I../common/calculateWaterLevel -IlnInclude -I. -I/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude -I/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OSspecific/POSIX/lnInclude -fPIC -Ddarwin -std=c++11 -c $SOURCE -o Make/darwinIntel64ClangDPOpt/IH_Waves_InletAlphaFvPatchScalarField.o
In file included from IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.C:47:
In file included from IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.H:277:
../common/memberFun.H:32:24: error: conversion from 'tmp<Field<cmptType> >' to 'scalarField'
(aka 'Field<scalar>') is ambiguous
return patch().Cf().component(0); // patchX
^~~~~~~~~~~~~~~~~~~~~~~~~
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/tmp.H:124:16: note: candidate
function
inline operator const T&() const;
^
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/Field.H:220:9: note: candidate
constructor
Field(const tmp<Field<Type> >&);
^
In file included from IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.C:47:
In file included from IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.H:277:
../common/memberFun.H:38:24: error: conversion from 'tmp<Field<cmptType> >' to 'scalarField'
(aka 'Field<scalar>') is ambiguous
return patch().Cf().component(1); // patchY
^~~~~~~~~~~~~~~~~~~~~~~~~
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/tmp.H:124:16: note: candidate
function
inline operator const T&() const;
^
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/Field.H:220:9: note: candidate
constructor
Field(const tmp<Field<Type> >&);
^
In file included from IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.C:47:
In file included from IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.H:277:
../common/memberFun.H:273:31: error: conversion from 'tmp<vectorField>' to 'const vectorField'
(aka 'const Field<vector>') is ambiguous
const vectorField nVecCell = patch().nf();
^ ~~~~~~~~~~~~
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/tmp.H:124:16: note: candidate
function
inline operator const T&() const;
^
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/Field.H:220:9: note: candidate
constructor
Field(const tmp<Field<Type> >&);
^
IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.C:350:23: error: conversion from
'tmp<Field<double> >' to 'const scalarField' (aka 'const Field<scalar>') is ambiguous
const scalarField alphaCell =
^
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/tmp.H:124:16: note: candidate
function
inline operator const T&() const;
^
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/Field.H:220:9: note: candidate
constructor
Field(const tmp<Field<Type> >&);
^
IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.C:352:23: error: conversion from
'tmp<Field<Foam::Vector<double> > >' to 'const vectorField' (aka 'const Field<vector>') is
ambiguous
const vectorField UCell = U.boundaryField()[patchID].patchInternalField();
^ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/tmp.H:124:16: note: candidate
function
inline operator const T&() const;
^
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/Field.H:220:9: note: candidate
constructor
Field(const tmp<Field<Type> >&);
^
IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.C:355:23: error: conversion from
'tmp<vectorField>' to 'const vectorField' (aka 'const Field<vector>') is ambiguous
const vectorField nVecCell = patch().nf();
^ ~~~~~~~~~~~~
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/tmp.H:124:16: note: candidate
function
inline operator const T&() const;
^
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/Field.H:220:9: note: candidate
constructor
Field(const tmp<Field<Type> >&);
^
IH_Waves_InletAlpha/IH_Waves_InletAlphaFvPatchScalarField.C:359:23: error: conversion from
'tmp<Field<scalar> >' to 'const scalarField' (aka 'const Field<scalar>') is ambiguous
const scalarField patchHeight = patch().Cf().component(2)-cMin[2];
^ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/tmp.H:124:16: note: candidate
function
inline operator const T&() const;
^
/Users/totalart/OpenFOAM/OpenFOAM-2.3.0/src/OpenFOAM/lnInclude/Field.H:220:9: note: candidate
constructor
Field(const tmp<Field<Type> >&);
^
7 errors generated.
make: *** [Make/darwinIntel64ClangDPOpt/IH_Waves_InletAlphaFvPatchScalarField.o] Error 1
\n\nWave generation boundary conditions compilation failed

[Phicau]

Hi Youngkook,

I am sorry, but I do not work with OpenFOAM in Mac, so I cannot be of help here.

Best,

Pablo

[ashim]

Hi Pablo,

I am trying to use IHFOAM for wave generation in tank. As a starting point, I was running the baseWaveFlume tutorial in foam extended and OF2.3.1. It runs okay but when I have just altered the patch inlet and outlet and make the wave direction 180 deg, it is not generating wave in negative X-direction. it is just diverge. Am I making very simple mistake?

Ali

[Antoniorp]

Hello Pablo,

I want to simulate the case of bichromatic waves. I have the have height and period for each of the two waves. I used that information in IHWavesDict with the waveType "irregular" but I'm not sure i should be using this.

I wonder if I should use the waveType "wavemaker" with a time series of eta, calculated by the combination of the two components?

Also, I want to simulate another wave, in this case random, of which I only have the Jonswap spectrum. How can i translate the information in the spectrum to Ihfoam?

Thank you,

António

[Phicau]

Hi there,

sorry for the late answers, but here they go.

@Ashim Yes, it is a simple mistake, don't forget to change absDir as well to 180 (or to any value greater than 360 for an automatic calculation of the absorption direction). It is running flawlessly for me in OF2.1.1.

@Antonio You are right, you want to generate an irregular case with two components. The wavemaker type could also be applied, but it is more suitable for constant inflow conditions (e.g. uniform inlet velocity). If you want to simulate a Jonswap spectrum, you just need to discretize it in a number of components (H,T,phase) and feed the irregular type with them.

Best regards,

Pablo

[Antoniorp]

Hello Pablo,

Thank you very much for your answer, it was a great help!

António

[Rising Star]

hi Phicau

I succesfully compiled IhFoam on OF 2.3.0, I've received the tutorials, but they're written only up ti OF 2.2.0. What I've to do to make them working with my version of OF?

Thanks in advance

Best Regards

[Phicau]

Hi Rising Star,

you have to adapt those provided to fit the new standards of the new version (i.e. changes in the names of variables, in the way to define the numerics...)

I suggest that you take the damBreak tutorial of version 2.2.2 and look how it changes in version 2.3.0.

For the record, I recommend you not to work with version 2.3...

Best,

Pablo

[ashim]

Hi Pablo,

Thank you very much for your reply. I am trying to simulate wave generation in thin domain having 10 cell per height and 50 cell per wave length for wave height 0.120m, water depth 2 m and wave length 3.6m. The wave is generating in negative x direction using IHFoam in OF2.3.1. The wave generation does not look good at the beginning although it recovers after few seconds. The main problem is the waves seem to be reflected from outlet boundary and distort the waves in domain. I have uploaded files in dropbox. Here is the link. Could you please let know me which is causing this problem. I have also observed the same behavior for waveflume case also.

https://www.dropbox.com/sh/uz6wu4rfb...dV5uA3w-a?dl=0

Best Regards,
Ali

[rafa13]

Hi Pablo an IHfoamers,

I am trying to simulate a wave breake and i am struggling with the porosity dict,

I don understand witch valores of alpha and beta are used to calculate the resistence coefficientes (a and b). I am trying first to recreate the Dam break an I read your paper, where the results suggests to use alpha =10000 and beta= 3 (porosity=0.49;d50=0.0159m), but when i use this same valores i do not get a= 50 and b=2. So i am asking how to determine a and b ?

Or are a and b...the resistence coefficients( alpha and beta) of the forcheimer equation??

sorry if thats a noob question but i'm struggling for some days now.

thanks
RM

[Phicau]

Hi Rafael,

the values in porosityDict are alpha and beta in the Forcheimer drag. The 2014 papers were performed with a different porous media implementation, that is why the values are so different, but the results are as good.

You determine a and b based on your experimental or analytical results, by best fit; or at least based on previous experience.

Best,

Pablo

[rafa13]

Hi Pablo,

Thanks a lot for your help and for this really quick reply ,
now you suggest to calibrate the alpha an beta for my case, beginning with a = 50 and b=2 than changing this coefficients, using the same procedure as your paper?

Ok thats my next step

Rafael Marques

[me.ouda]

Hi Pablo,

I am trying to use the wave generation and absorption boundary conditions with other solver and I started by interFoam and I am planning to try also with interMixingFoam to simulate 3 phase flow because I am trying to simulate sediment transport under wave action. I have compiled ihFoam on Ubuntu 14.04 with OpenFoam 2.2.2. I am trying to run the baseWaveFlume tutorial with the same configurations except that I run it with interFoam. To do that I followed the instructions in the manual and included the required libraries to the controlDict. But when I tried to run the case I got the following error:

--> FOAM FATAL IO ERROR: 
keyword div((muEff*dev(T(grad(U))))) is undefined in dictionary "/home/meouda/OpenFOAM/meouda-2.2.2/run/tutorials/multiphase/IHFOAM/baseWaveFlume_2/system/fvSchemes.divSchemes"

file: /home/meouda/OpenFOAM/meouda-2.2.2/run/tutorials/multiphase/IHFOAM/baseWaveFlume_2/system/fvSchemes.divSchemes from line 32 to line 35.

    From function dictionary::lookupEntry(const word&, bool, bool) const
    in file db/dictionary/dictionary.C at line 402.

FOAM exiting

Simulation complete.

So could you please help me to solve this problem?

Best regards

[Phicau]

Hi,

the case is prepared to be run with IHFOAM, therefore, the schemes have porosity included. If you decide to run it with interFoam you need to define the scheme to use for the term:

div((muEff*dev(T(grad(U)))))

change your fvSchemes to include this one, with the same value as the one that looks similar but includes porosity.

Best,

Pablo

[me.ouda]

Hi Pablo again,

It is my mistake that I didn't search before asking, sorry. I've found the solution by adding the following line into the system/fvSchemes file in divSchemes:

div((muEff*dev(T(grad(U))))) Gauss linear;

The solution was reported by majkl here

But is this error means that I have problem in OpenFOAM or ihFoam compilation?

Best regards.

[me.ouda]

Originally Posted by Phicau

Hi,

the case is prepared to be run with IHFOAM, therefore, the schemes have porosity included. If you decide to run it with interFoam you need to define the scheme to use for the term:

div((muEff*dev(T(grad(U)))))

change your fvSchemes to include this one, with the same value as the one that looks similar but includes porosity.

Best,

Pablo

Seems that we reply at the same time

many thanks for your interest