Problem of fluid particles splash during simulation

dongdong
in DualSPHysics v4.2
Hi guys:
I was recently troubled with a problem of particles splashing during simulation. There are some floatings moving in the water and generate waves , the phenomenon in the previous section is correct ,but some particles will splash and into the air, such as follows:


Is this problem relative to any parameters?

Comments

  • Asalih3dAsalih3d
    Hi Dong

    Could you please describe more in detail what we are looking at? It can be hard to provide any useful insight without understanding your situation enough.

    I will note that I have beforehand experienced "splashing" and in my case it was realistic since the water got disconnected from the main volume and starting mixing with air.
  • dongdong
    @Asalih3d

    In this simulation a stream of water flows on the floatings, and some fluid particles fly away. Is this realistic?
  • AlexAlex
    Try to use a higher viscosity value
    OR
    Try to use double precision
  • dongdong
    @Alex
    Now I have used double precision and I use laminar+SPS , the value is 1e-6，just like water.

  • AlexAlex
    Try to use artificial viscosity so that you can control better those instabilities.
    Laminar viscosity is a good example of programming a ACTUAL and PHYSICAL viscosity, however most of the times is not numerically stable, so I suggest to move to artificial viscosity.
    Other viscosity treatments such as k-eps should work in a better way than the laminar formulation we have implemented

    Regards
  • dongdong
    @Alex
    If I use artificial viscosity and set the value 0.01 or 0.001. When I compute the force with ComputeForces4_linux64 if I can get the right result?
  • AlexAlex
    Artificial viscosity will affect the velocity propagation of the dam (of volume of fluid) so that impact time should vary, however the force magnitude should be corrected.

    Note that the force = pressure terms + viscous forces, and generally pressure terms are much higher (at least that seems in your problem)

    Regards
  • dongdong
    @Alex
    But how to correct the force if I use artifical viscosity? Is there any reference or principle I can use?
  • AlexAlex
    You do not need to correct the force. Force will be computed using that value of viscosity.
    Please try it and tell us if it works or not and we can help you again.
  • dongdong
    Thanks all
    I have solved the problem by changing the RigidAlgrithm. I use the dem so I must set RigidAlgrithm=2(DEM).
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('@Asalih3d <br /...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('@Asalih3d <br /...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('@Asalih3d <br /...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('@Asalih3d \r\n<im...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('@Asalih3d \r\n<im...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('@Asalih3d \r\n<im...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('@Asalih3d \r\n<im...', '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/applications/vanilla/controllers/class.discussioncontroller.php(402): DiscussionController->index(1472, 'x', 'p1')
#17 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(862): DiscussionController->comment('2448')
#18 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(279): Gdn_Dispatcher->dispatchController(Object(Gdn_Request), Array)
#19 /var/www/forums-dual-sphysics-org/index.php(29): Gdn_Dispatcher->dispatch()
#20 {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('@Alex <br />Now...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('@Alex <br />Now...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('@Alex <br />Now...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('@Alex \r\nNow I h...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('@Alex \r\nNow I h...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('@Alex \r\nNow I h...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('@Alex \r\nNow I h...', '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/applications/vanilla/controllers/class.discussioncontroller.php(402): DiscussionController->index(1472, 'x', 'p1')
#17 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(862): DiscussionController->comment('2448')
#18 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(279): Gdn_Dispatcher->dispatchController(Object(Gdn_Request), Array)
#19 /var/www/forums-dual-sphysics-org/index.php(29): Gdn_Dispatcher->dispatch()
#20 {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('@Alex <br />If ...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('@Alex <br />If ...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('@Alex <br />If ...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('@Alex \r\nIf I us...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('@Alex \r\nIf I us...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('@Alex \r\nIf I us...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('@Alex \r\nIf I us...', '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/applications/vanilla/controllers/class.discussioncontroller.php(402): DiscussionController->index(1472, 'x', 'p1')
#17 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(862): DiscussionController->comment('2448')
#18 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(279): Gdn_Dispatcher->dispatchController(Object(Gdn_Request), Array)
#19 /var/www/forums-dual-sphysics-org/index.php(29): Gdn_Dispatcher->dispatch()
#20 {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('@Alex <br />But...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('@Alex <br />But...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('@Alex <br />But...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('@Alex \r\nBut how...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('@Alex \r\nBut how...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('@Alex \r\nBut how...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('@Alex \r\nBut how...', '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), 8)
#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/applications/vanilla/controllers/class.discussioncontroller.php(402): DiscussionController->index(1472, 'x', 'p1')
#17 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(862): DiscussionController->comment('2448')
#18 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(279): Gdn_Dispatcher->dispatchController(Object(Gdn_Request), Array)
#19 /var/www/forums-dual-sphysics-org/index.php(29): Gdn_Dispatcher->dispatch()
#20 {main}

© Vanilla Keystone Theme 2025