Fluid particles penetrate boundary

duanduan
in DualSPHysics v4.2
Hi, everyone!

I am new to dualsphysics and I am testing the dambreak case with Freecad following https://design.sphysics.org/wiki/doku.php?id=dam-break-example.

I just found a strange behavior that in the initial condition, the fluid particles penetrate the boundary, please see the attached figure, blue is boundary and red are the fluid particles.



To give more details, I changed the case to 2D and the boundary wall have thickness 10 mm while the dp is set 5 mm.

does anyone have any ideas to avoid it?

Regards,

Duan




Comments

  • bladomas1bladomas1
    It probably has to do something with the order of initialization. If you create the wall first and then the fluid, the fluid particles can replace the wall particles if the two initialization domains slightly overlap at the edges.
    To give more details, I changed the case to 2D and the boundary wall have thickness 10 mm while the dp is set 5 mm.
    Also I do not understand this concern. You set the wall to have a thickness of 10 mm, and set the particle spacing to 5mm, which means 2 particles will be "bound" and act as the wall.
  • duanduan
    @bladomas1

    Thank you for reply.

    What I did not understand is exactly why the fluid particles can replace the wall particles, because when I did the initialization, there is supposed to be no overlap, while in fact there is a overlap.

    I found that I could make a small distance (equals to dp) between the fluid and water to avoid the overlap. But I am not sure it is the best solution or not.
    You set the wall to have a thickness of 10 mm, and set the particle spacing to 5mm, which means 2 particles will be "bound" and act as the wall.
    This is another doubt. By setting the wall thickness = 10 mm, and set the particle spacing to 5mm, "2" particles are supposed to be "bound" and act as the wall. However as shown in the picture, there is "3" particles bound as the wall. Why?






  • Asalih3dAsalih3d
    There are three particles since, it starts from 0, then goes to 5 and then goes to 10 mm. Sometimes because of numerical error you would only get two because it might think that it starts at 0,001 or something. Getting three particles means that it has been done correctly in this case.

    The fluid particles can penetrate the boundary if the particle spacing is too big, which I believe is what happened here. If the particle spacing was 1 mm I don't think you would encounter a problem.

    Kind regards
  • bladomas1bladomas1
    @duan Ah, I made the same mistake as you. With thickness =10mm and dp=5mm, there are 2 gaps present. Kind of hard to explain. Lets say particle A is placed at point 1=at 0mm, then there is a gap of 5mm, and point 2 particle=at 5mm, another gap and last particle/point=at 10mm. This means 3 particles were placed, with 2 gaps in total.
  • duanduan
    @Asalih3d
    The fluid particles can penetrate the boundary if the particle spacing is too big, which I believe is what happened here. If the particle spacing was 1 mm I don't think you would encounter a problem.
    I tried with 1 mm, the problem remains at z=0 mm, see the attached figure.
    Since the wall is set to (-10 mm to 0 mm) and the water (0 mm to 300 mm) on Z direction, I think it is probably both fluid particles and wall particles starts from z=0 mm, so there is one layer particle overlap at z=0 mm, no matter how small the dp is. And if the water is set to (dp mm to 300 mm), the overlap is gone.


  • duanduan
    @bladomas1
    There are three particles since, it starts from 0, then goes to 5 and then goes to 10 mm.
    I think this is the right explaination.
  • Asalih3dAsalih3d
    Yes, seems like you just have to move the fluid a little bit higher up.

    Kind regards
  • TPouzolTPouzol
    Or as bladomas1 suggested in his first post, you should change the order of creation of the mks.
    Good day
Sign In or Register to comment.

Debug Trace

Notice 
Cannot assign an empty string to a string offset

#0 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1203): gdn_ErrorHandler(2, 'Cannot assign a...', '/var/www/forums...', 1203, Array)
#1 /var/www/forums-dual-sphysics-org/library/core/class.format.php(958): Gdn_Format::formatMentionsCallback('')
#2 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1292): Gdn_Format::tagContent('@bladomas1 <br ...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('@bladomas1 <br ...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('@bladomas1 <br ...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('@bladomas1 \r\n\r\n...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('@bladomas1 \r\n\r\n...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('@bladomas1 \r\n\r\n...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('@bladomas1 \r\n\r\n...', 'Html')
#9 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(170): formatBody(Object(stdClass))
#10 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/comments.php(19): writeComment(Object(stdClass), Object(DiscussionController), Object(Gdn_Session), 2)
#11 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/index.php(53): include('/var/www/forums...')
#12 /var/www/forums-dual-sphysics-org/library/core/class.controller.php(778): include('/var/www/forums...')
#13 /var/www/forums-dual-sphysics-org/library/core/class.controller.php(1382): Gdn_Controller->fetchView('', false, false)
#14 /var/www/forums-dual-sphysics-org/library/core/class.pluggable.php(217): Gdn_Controller->xRender()
#15 /var/www/forums-dual-sphysics-org/applications/vanilla/controllers/class.discussioncontroller.php(310): Gdn_Pluggable->__call('render', Array)
#16 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(862): DiscussionController->index('1549', 'fluid-particles...', '')
#17 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(279): Gdn_Dispatcher->dispatchController(Object(Gdn_Request), Array)
#18 /var/www/forums-dual-sphysics-org/index.php(29): Gdn_Dispatcher->dispatch()
#19 {main}
Notice 
Cannot assign an empty string to a string offset

#0 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1203): gdn_ErrorHandler(2, 'Cannot assign a...', '/var/www/forums...', 1203, Array)
#1 /var/www/forums-dual-sphysics-org/library/core/class.format.php(958): Gdn_Format::formatMentionsCallback('')
#2 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1292): Gdn_Format::tagContent('@duan Ah, I mad...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('@duan Ah, I mad...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('@duan Ah, I mad...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('@duan Ah, I mad...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('@duan Ah, I mad...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('@duan Ah, I mad...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('@duan Ah, I mad...', 'Html')
#9 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(170): formatBody(Object(stdClass))
#10 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/comments.php(19): writeComment(Object(stdClass), Object(DiscussionController), Object(Gdn_Session), 4)
#11 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/index.php(53): include('/var/www/forums...')
#12 /var/www/forums-dual-sphysics-org/library/core/class.controller.php(778): include('/var/www/forums...')
#13 /var/www/forums-dual-sphysics-org/library/core/class.controller.php(1382): Gdn_Controller->fetchView('', false, false)
#14 /var/www/forums-dual-sphysics-org/library/core/class.pluggable.php(217): Gdn_Controller->xRender()
#15 /var/www/forums-dual-sphysics-org/applications/vanilla/controllers/class.discussioncontroller.php(310): Gdn_Pluggable->__call('render', Array)
#16 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(862): DiscussionController->index('1549', 'fluid-particles...', '')
#17 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(279): Gdn_Dispatcher->dispatchController(Object(Gdn_Request), Array)
#18 /var/www/forums-dual-sphysics-org/index.php(29): Gdn_Dispatcher->dispatch()
#19 {main}
Notice 
Cannot assign an empty string to a string offset

#0 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1203): gdn_ErrorHandler(2, 'Cannot assign a...', '/var/www/forums...', 1203, Array)
#1 /var/www/forums-dual-sphysics-org/library/core/class.format.php(958): Gdn_Format::formatMentionsCallback('')
#2 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1292): Gdn_Format::tagContent('@Asalih3d \n<blo...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('@Asalih3d \n<blo...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('@Asalih3d \n<blo...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('@Asalih3d \r\n\r\n<...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('@Asalih3d \r\n\r\n<...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('@Asalih3d \r\n\r\n<...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('@Asalih3d \r\n\r\n<...', 'Html')
#9 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(170): formatBody(Object(stdClass))
#10 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/comments.php(19): writeComment(Object(stdClass), Object(DiscussionController), Object(Gdn_Session), 5)
#11 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/index.php(53): include('/var/www/forums...')
#12 /var/www/forums-dual-sphysics-org/library/core/class.controller.php(778): include('/var/www/forums...')
#13 /var/www/forums-dual-sphysics-org/library/core/class.controller.php(1382): Gdn_Controller->fetchView('', false, false)
#14 /var/www/forums-dual-sphysics-org/library/core/class.pluggable.php(217): Gdn_Controller->xRender()
#15 /var/www/forums-dual-sphysics-org/applications/vanilla/controllers/class.discussioncontroller.php(310): Gdn_Pluggable->__call('render', Array)
#16 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(862): DiscussionController->index('1549', 'fluid-particles...', '')
#17 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(279): Gdn_Dispatcher->dispatchController(Object(Gdn_Request), Array)
#18 /var/www/forums-dual-sphysics-org/index.php(29): Gdn_Dispatcher->dispatch()
#19 {main}
Notice 
Cannot assign an empty string to a string offset

#0 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1203): gdn_ErrorHandler(2, 'Cannot assign a...', '/var/www/forums...', 1203, Array)
#1 /var/www/forums-dual-sphysics-org/library/core/class.format.php(958): Gdn_Format::formatMentionsCallback('')
#2 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1292): Gdn_Format::tagContent('@bladomas1 \n<bl...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('@bladomas1 \n<bl...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('@bladomas1 \n<bl...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('@bladomas1 \r\n\r\n...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('@bladomas1 \r\n\r\n...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('@bladomas1 \r\n\r\n...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('@bladomas1 \r\n\r\n...', 'Html')
#9 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(170): formatBody(Object(stdClass))
#10 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/comments.php(19): writeComment(Object(stdClass), Object(DiscussionController), Object(Gdn_Session), 6)
#11 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/index.php(53): include('/var/www/forums...')
#12 /var/www/forums-dual-sphysics-org/library/core/class.controller.php(778): include('/var/www/forums...')
#13 /var/www/forums-dual-sphysics-org/library/core/class.controller.php(1382): Gdn_Controller->fetchView('', false, false)
#14 /var/www/forums-dual-sphysics-org/library/core/class.pluggable.php(217): Gdn_Controller->xRender()
#15 /var/www/forums-dual-sphysics-org/applications/vanilla/controllers/class.discussioncontroller.php(310): Gdn_Pluggable->__call('render', Array)
#16 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(862): DiscussionController->index('1549', 'fluid-particles...', '')
#17 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(279): Gdn_Dispatcher->dispatchController(Object(Gdn_Request), Array)
#18 /var/www/forums-dual-sphysics-org/index.php(29): Gdn_Dispatcher->dispatch()
#19 {main}

© Vanilla Keystone Theme 2025