open TELEMAC-MASCARET The mathematically superior suite of solvers

just look at your screenshot and you will see that all lines for polygon4 are unaligned with the previous one!

I edited and now the line after YSOM3(49) is aligned, but still it gives same error. you can see the alignment (attached).

could you post the last error message?

Thanks for your reply.
I attached the error

Could you try to use your brain a little and read the error message instead of just posting and waiting for the answer !
The problem in is corpor!!
No longer in dragfo !!

As I already wrote previously:

correct!!
coding should start in column 7 or greater so that it can be read

Hi, I am using drag coefficient and porosity fortran subroutine to simulate a wall impact in the domain. The result without the wall looks ok. But when I run with the drag, it shows irregularity in free surface. The expected situation is that the drag will increase the free surface and decrease the velocity. But there is a suddenn upward or downward spike in free surface which can be an error. I used mesh size 0.25 where the drag is working. I attached some cross sectional value in here for showing the result. This is a 2d model. Do you have any idea how to fix it or what are the reasons behind this? Is it mesh size? Thanks.

Hello,

I am facing a similar problem, and thanks to previous post from this topic I rewrote user_dragfo.f to have two areas: I am trying to model a turbine made of 2 rotors in T2D in a uniform flow with no friction at the bottom.

However, my wake often remains dysymmetric:



uniform 1[m] mesh: geo_dragforce_back1_rotor0d2

This dysymmetry/wake is (strongly) mesh dependant:
- using a coarser mesh we/I get something else:




uniform 5[m] mesh: geo_dragforce_back5_rotor0d2

- Or using boxes - instead of uniform but still with a 1m cell-size in area of interest, wake is no longer dysymmetric:




mesh with 1[m] box: geo_dragforce_back10_box1_rotor0d2

Stuff I already tried and made no difference:
(- rewrite dragfo with the two distincts rotors)
- use as constant Unorm=U_inf=5
- tried ks=30, make no diff
- tried another turbulence model: Smagorinsky made no difference (as expected)
- ...

Parameters I don't want to change:
- friction coefficient (because it aims to reproduce a 2D-CFD model)


I have attached the case here: https://gitlab.com/phmusiedlak/public/-/raw/main/drag_force_2rotors_.zip?ref_type=heads

Any hints would be good :)

Have a nice day

Pierre-Henri

Hello,

I am facing a similar problem, and thanks to previous post from this topic I rewrote user_dragfo.f to have two areas: I am trying to model a turbine made of 2 rotors in T2D in a uniform flow with no friction at the bottom.

However, my wake often remains dysymmetric:



uniform 1[m] mesh: geo_dragforce_back1_rotor0d2

This dysymmetry/wake is (strongly) mesh dependant:
- using a coarser mesh we/I get something else:




uniform 5[m] mesh: geo_dragforce_back5_rotor0d2

- Or using boxes - instead of uniform but still with a 1m cell-size in area of interest, wake is no longer dysymmetric:




mesh with 1[m] box: geo_dragforce_back10_box1_rotor0d2

Stuff I already tried and made no difference:
(- rewrite dragfo with the two distincts rotors)
- use as constant Unorm=U_inf=5
- tried ks=30, make no diff
- tried another turbulence model: Smagorinsky made no difference (as expected)
- ...

Parameters I don't want to change:
- friction coefficient (because it aims to reproduce a 2D-CFD model)


I have attached the case here: https://gitlab.com/phmusiedlak/public/-/raw/main/drag_force_2rotors_.zip?ref_type=heads

Any hints would be good :)

Have a nice day

Pierre-Henri