Calculating Artificial Viscosity

swaraj1112swaraj1112
in DualSPHysics v4.2
Is there any analytical formula to calculate artificial viscosity from kinematic viscosity and density ?

Comments

  • AlexAlex
    Some papers of Colagrossi et al defines a relationship... However that formula is only valid for some problems such as wave tanks (where a value of alpha=0.01 has been shown to provide good results), but it is not valid for dam-breaks or sloshing problems using DualSPHysics since alpha will vary with particle resolution.

    Regards
  • swaraj1112swaraj1112
    Ok @Alex. I want to simulate a Molten Metal Flow through a casting Mould and all i am getting is the artificial viscosity of water.
  • Asalih3dAsalih3d
    @swaraj1112

    Then you will have to find some physical experiments of these cases, use the artificial viscosity values until you hit the same experimental results. Then you can use this artificial viscosity you've found for your own case. Note that this is still an approximation, since you have not conducted an experiment your self (?) - if you have you can use that instead.

    As a final note, artificial viscosity, is not a property any fluid has - it should be understood as a tuning factor, which is able to stabilize SPH simulation without enormous numbers of particles and get results comparable to the experiments which have been conducted.

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('@swaraj1112 <br...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('@swaraj1112 <br...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('@swaraj1112 <br...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('@swaraj1112 \r\n\r...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('@swaraj1112 \r\n\r...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('@swaraj1112 \r\n\r...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('@swaraj1112 \r\n\r...', '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), 3)
#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('1477', 'calculating-art...', '')
#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