open TELEMAC-MASCARET The mathematically superior suite of solvers

Thanks very much Christophe, that was really helpful.

Kind regards,

Niamh

Hi Christophe,

I have another question, is it possible in ARTEMIS to assign different reflection coefficients (Kr values) to the boundaries? In my model some of the walls bordering the channel are vertical while others are rock revetments which should absorb some of the wave energy rather then reflecting it.

Thanks for your help,

Niamh

Hi Niamh,

You can impose different reflexion coefficients to your boundaries using the fortran file borh.f, RP%R variable. It works as for phase definition on the boundary (loop on the boundary nodes).
As before, the test case btwi (file btwi-par.f) gives an example and should help you.

Regards,

Christophe

Hi Christophe,

I've got a question regarding the borh.f (niamh.f) file you attached for the download. I actually got the same problem as Niamh in terms of wave direction and was wondering if you have any particular advice or information regarding the phase and especially the wavenumber calculation-loop. This loop:

C --- WAVE LENGTH FROM INFINITY
C COMPUTE KINF FROM OMEGA AND DINF (explicit method)
DINF =5.2D0
X0 =0.
Y0 =0.

AUX1 =DINF*OMEGA**2/GRAV
AUX2 = 1.D0 + AUX1 *( 0.6522D0 +
& AUX1 *( 0.4622D0 +
& AUX1 *
& AUX1 *( 0.0864D0 +
& AUX1 *( 0.0675D0 ) )))
KINF = SQRT( AUX1*(AUX1 + 1.D0/AUX2) )
KINF = KINF/DINF

was just copy and pasted in niamh's file from the btwi example, as far I could see.
I tried to find any further explanation or documentation to understand what's behind this and how the wavenumber is calculated but I even didn't find a proper documentation about the variables. For instance, can you explain what is AUX1 and AUX2? Is this always the same loop no matter what the conditions (direction and HB) change? (except the predefined depth in DINF)? As far I can see it uses the given values in the both.f file, can you confirm this?

Any explanation would be appreciated, I'm trying to get a touch of Artemis and it feels a bit confusing at the moment.

Cheers, Tim

Hey guys,

I came another problem across. My simulation of monochromatic waves is fine after adopting Christophe's advices to my model. After the mono wave case I tried to run the mono-direction random wave with a jonswap spec. between 3. to 8. sec. I find it quite hard to interpret the waveheigt output compared to the free-surface elev.. Anyhow, after I ran my mono-directional random wave the output looked very strange, as I got strong gradients perpendicular to my actual direction so it looks like the waves are running rotated by 90° to my actual prescribed direction.

Does anyone has any idea how to correct define the wave direction, or is this the way it should look like...it feels strange anyway.
As additional information my wave should enter the domain at 155°. Attached you can find my three examples as screenshots.

Any help or idea how to improve the visualization or the simulation would be much appreciated.

Cheers Tim