TELEMAC-MASCARET Forum: Flow pattern of source points without velocity (1/1)

Welcome, Guest

TOPIC: Flow pattern of source points without velocity

Flow pattern of source points without velocity 5 years 8 months ago #33064

Binguo
OFFLINE
Junior Boarder
Posts: 34
Thank you received: 1

Dear all,

I have a question regarding the source points, hope someone could help me.

I put several source points in the vicinity of a solid boundary (see the circle in the figure below). Then discharge was add to sources without velocity. Yet the flow pattern around the source points in result beyond my understanding.

Why the water flow along the solid boundary rather than to the direction againt the wall? Is it right? Should I do some setting about the boundary? Thanks in advance.

Best Regards,
Bin

The administrator has disabled public write access.

Flow pattern of source points without velocity 5 years 8 months ago #33075

riadh	Hello On finite element Wall boundary condition is prescribed in a weak way. This means that impermeabilty of the wall is obtained in an integral way. This explains velocity pattern you get. It is not recommended to put sources close to solid boundary. If it is necessary prescribe them using velocity. You can also use dirac option (source option=2)
	The administrator has disabled public write access.

Flow pattern of source points without velocity 5 years 8 months ago #33099

Binguo
OFFLINE
Junior Boarder
Posts: 34
Thank you received: 1

Dear Riadh, Many thanks for your quick answer.

I have used the direct option but nothing change. This is weird because it should have some change for the next code in prosou.f and the area of my element is not equal to 1.

IF(OPTSOU.EQ.1) THEN
!             "NORMAL" VERSION
              SMH%R(IR)=SMH%R(IR)-DBUS(I)*UNSV2D%R(IR)
            ELSE
!             "DIRAC" VERSION
              SMH%R(IR)=SMH%R(IR)-DBUS(I)
            ENDIF

Sorry I still don't understand why the intergral of the boundary result in this problem. Do you mind to give me more detailes about this?

Regards,
Bin

The administrator has disabled public write access.