search

Xiangyu-Hu / SPHinXsys

SPHinXsys provides C++ APIs for engineering simulation and optimization. It aims at complex systems driven by fluid, structure, multi-body dynamics and beyond. The multi-physics library is based on a unique and unified computational framework by which strong coupling has been achieved for all involved physics.
https://www.sphinxsys.org/
Apache License 2.0
259 stars 199 forks source link

Floating body test case seems not right #289

Closed Xiangyu-Hu closed 12 months ago

Xiangyu-Hu commented 12 months ago

The 2d floating body problem seems not right. Screenshot from 2023-05-18 16-08-23 There are quite some flow around the floating body.

ncsalis commented 12 months ago

I have updated the master and tested. In the mac system seems ok. floating_body

I will also check in linux system.

ncsalis commented 12 months ago

Also in linux system is working ok. floting_box

Maybe the bug is for the windows system. For curiosity have you tried the 3D floating body case? (3D cases kept failing in the windows CI test for the pull request and I have increased the dtwdistance)

ncsalis commented 12 months ago

pull request #267

Xiangyu-Hu commented 12 months ago

I am using Linux system. Let me check again.

Xiangyu-Hu commented 12 months ago

Which master did you do the above computation, and the resolution?

ncsalis commented 12 months ago

I updated the master with today merge before running. The resolution used is L/20.

Xiangyu-Hu commented 12 months ago

It is strange. I am using rhe same master on Linux. Same resolution as the master.

ncsalis commented 12 months ago

I have tried with lower resolution (L/10) and solution diverges. Tomorrow I will search for possible reasons

anim_l10

Xiangyu-Hu commented 12 months ago

I also use resolution L/20.

ncsalis commented 12 months ago

It is strange that does not work in different machines and/or with different resolution. Maybe tomorrow we can discuss it

Xiangyu-Hu commented 12 months ago

I have found that I am using Intel compiler.

Xiangyu-Hu commented 12 months ago

The problem may be due to some small number issue as difference compiler give quite different results.

Xiangyu-Hu commented 12 months ago

I mean probably origining from initial particle position.