chaos moving of fluid particles in 3d model

mtmtmtmtmtmt
edited March 2019 in DualSPHysics v4.2
hi i simulate a tank of water in 3d model in case DEM but when i run case and see my simulation the fluid particles have chaos movement & i dont know why!
i working in 3D model recently and i cant realize this problem. i use laminar+sps viscosity in my model.


Comments

  • Asalih3dAsalih3d
    It is because you have not activated delta-SPH scheme. Because of this you get some density fluctuations which lead to chaotic movement and finally unstability. If you want to know more read Dalrymple 2004.

    Kind regards
  • mtmtmtmtmtmt
    Dear @Asalih3d in my simulation DELTASPH=0.1 and when i use artificial viscosity this problem would be solved but when use laminar+sps and v=0.000001 problem occur. are you sure this problem relevant to deltasph term?
    Best Regards
  • Asalih3dAsalih3d
    In my case, delta-SPH has played a huge role in stabilizing SPS-LES simulations since they handle the artificial density increase due to SPS-LES. Another parameter would be the number of particles, what happens if you decrease and increase the amount of particles?

    A good test is always to try the artifical scheme and see if it is stable in that and then try fine tuning the laminar+sps scheme afterwards.

    Kind regards
  • AlexAlex
    Laminar viscosity + SPS turbulence model should be used only in the case of simulating a laminar flow and with a very high resolution... If not, it is recommend to use artificial viscosity where you can obtain more stable simulations.
    Alpha=0.01 is recommended for wave tanks.

    Regards
  • mtmtmtmtmtmt
    @Asalih3d Thanks For Your Advice I try your Idea.
    best Regards
  • mtmtmtmtmtmt
    edited March 2019
    Dear @Alex What do you mean about High Resolution? Because The Size Of Domain Relevent To Dp and When Create Small Domain actually use smaller dp. in my simulation dp=0.005 and size of domain x=0.5 m y=0.040 z=0.06? and when run 2D simulation with Laminar+SPS this chaotic movement very decrease!!! Why?

  • mtmtmtmtmtmt
    Dear @Alex I Increase The Y direction(double of last simulation -> y=0.080) & this problem is fix, is there any relation between this problem and resolution?
  • AlexAlex
    SPS means Sub-Particle Scale model so that this model only solves turbulence at a very small scale
  • mtmtmtmtmtmt
    @Alex Thanks Dear Professor
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 Thank...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('@Asalih3d Thank...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('@Asalih3d Thank...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('@Asalih3d Thank...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('@Asalih3d Thank...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('@Asalih3d Thank...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('@Asalih3d Thank...', '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/applications/vanilla/controllers/class.discussioncontroller.php(402): DiscussionController->index(1545, 'x', 'p1')
#17 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(862): DiscussionController->comment('2730')
#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 Thanks De...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('@Alex Thanks De...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('@Alex Thanks De...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('@Alex Thanks De...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('@Alex Thanks De...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('@Alex Thanks De...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('@Alex Thanks De...', '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), 9)
#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(1545, 'x', 'p1')
#17 /var/www/forums-dual-sphysics-org/library/core/class.dispatcher.php(862): DiscussionController->comment('2730')
#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