New DualSPHysics v3.2
We are pleased to announce the release of new packages of DualSPHysics (v3.2) that includes the executable of the multi-phase code by Fourtakas et al. (2014) that deals with two-phase liquid-solid problems with flow-induced erosion of fully saturated sediment. The model predicts multi-phase sediment scouring phenomena that are induced by rapid flows creating shear forces at the surface of the sediment which causes the surface to yield and produce a shear layer of sediment suspended particles at the interface and finally sediment suspension in the fluid.
- This is a beta version, so support and stuff on the forum is limited and you should read the readme file ("MPhase_V3.0_Beta_README.txt") and the PDF ("Guide_CaseTwoPhases.pdf") included in the new version.
- A new testcase ("CASETWOPHASES") is also included to be executed with the new executable (“DualSPHysics_win64_MPhase”).
REFERENCE: Fourtakas, G., B.D. Rogers, and D.R. Laurence. 3-D SPH Modelling of Sediment Scouring Induced by Rapid Flows. in 9th International SPHERIC SPH Workshop. 2014. Paris, France.
- This is a beta version, so support and stuff on the forum is limited and you should read the readme file ("MPhase_V3.0_Beta_README.txt") and the PDF ("Guide_CaseTwoPhases.pdf") included in the new version.
- A new testcase ("CASETWOPHASES") is also included to be executed with the new executable (“DualSPHysics_win64_MPhase”).
REFERENCE: Fourtakas, G., B.D. Rogers, and D.R. Laurence. 3-D SPH Modelling of Sediment Scouring Induced by Rapid Flows. in 9th International SPHERIC SPH Workshop. 2014. Paris, France.
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('Hi All,<br /><b...', 'Gdn_Format::for...')
#3 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Html/HtmlEnhancer.php(59): Gdn_Format::mentions('Hi All,<br /><b...')
#4 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/Formats/HtmlFormat.php(74): Vanilla\Formatting\Html\HtmlEnhancer->enhance('Hi All,<br /><b...')
#5 /var/www/forums-dual-sphysics-org/library/Vanilla/Formatting/FormatService.php(34): Vanilla\Formatting\Formats\HtmlFormat->renderHtml('Hi All,\r\n\r\nI cu...')
#6 /var/www/forums-dual-sphysics-org/library/core/class.format.php(769): Vanilla\Formatting\FormatService->renderHTML('Hi All,\r\n\r\nI cu...', 'html')
#7 /var/www/forums-dual-sphysics-org/library/core/class.format.php(1479): Gdn_Format::html('Hi All,\r\n\r\nI cu...')
#8 /var/www/forums-dual-sphysics-org/applications/vanilla/views/discussion/helper_functions.php(24): Gdn_Format::to('Hi All,\r\n\r\nI cu...', '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), 10)
#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('566', 'new-dualsphysic...', '')
#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} |
Comments
Thank you for the update. I downloaded the new version and tried to run the two-phase case. However it stoped after the initialising step. I can see that for the raw case only 3506 particles so it seems not the problem of computational efficiency but more like something like a endless loop. Could you have a look at it? Thanks!
Yi Pan
Regards
Copyright (C) 2015 by Dr Georgios Fourtakas,
Dr Jose M. Dominguez, Dr Benedict Rogers, Dr Alejandro Crespo,
Prof. M. Gomez Gesteira, Dr Anxo Barreiro, Dr Ricardo Canelas,
Dr Renato Vacondio.
EPHYSLAB Environmental Physics Laboratory, Universidade de Vigo
School of Mechanical, Aerospace and Civil Engineering, University of Manchester
DualSPHysics is free software: you can redistribute it and/or
modify it under the terms of the GNU General Public License as
published by the Free Software Foundation, either version 3 of
the License, or (at your option) any later version.
DualSPHysics is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License,
along with DualSPHysics. If not, see .
DUALSPHYSICS for MultiPhase V3.0 (22-09-2015)
==============================================
[Select CUDA Device]
Device 0: "GeForce GTX 760"
Compute capability: 3.0
Multiprocessors: 6 (1152 cores)
Memory global: 2048 MB
Clock rate: 1.07 GHz
Run time limit on kernels: Yes
ECC support enabled: No
[GPU Hardware]
Device default: 0 "GeForce GTX 760"
Compute capability: 3.0
Memory global: 2048 MB
Memory shared: 49152 Bytes
[Initialising JSphGpuSingle v3.00 30-09-2015 16:29:50]
[Multi-phase configuration]
PhaseCount=2
Phase_0-> mkfluid:0 idbegin:800 count:594 rhop0:1090119490620926500000000000.000000 csound:72250710197730502000000000000.000000 gamma:0.000000
Phase_1-> mkfluid:1 idbegin:1394 count:2112 rhop0:198560776821242030000000000000.000000 csound:4592633630852250300000000000.000000 gamma:68872225043086548000000.000000
**Case configuration is loaded
Loading initial state of particles...
Loaded particles: 3506 (3506+0)
MapPos(border)=(0.002461,0.999961,0.002461)-(4.007539,1.000039,4.007539)
MapPos=(0.002461,0.999961,0.002461)-(4.007539,1.000039,24.032930)
**Initial state of particles is loaded
**2D-Simulation parameters:
CaseName="CaseTwoPhases"
RunName="CaseTwoPhases"
SvTimers=True
StepAlgorithm="Symplectic"
Kernel="Wendland"
Viscosity="SumSPS"
Visco=0.000100
ShepardSteps=0
DeltaSph="DBC"
DeltaSphValue=0.100000
CaseNp=3506
CaseNbound=800
CaseNfixed=800
CaseNmoving=0
CaseNfloat=0
PeriodicActive=0
Dx=0.030000
H=0.039032
CteB=554965.687500
Gamma=7.000000
Rhop0=1000.000000
Eps=0.000000
Cs0=3461172224.000000
CFLnumber=0.100000
DtIni=0.000000
DtMin=0.000000
MassFluid=0.900000
MassBound=0.450000
WendlandCte.awen=365.600677
WendlandCte.bwen=-46836.468750
TimeMax=2.000000
TimePart=0.010000
Gravity=(0.000000,0.000000,-9.810000)
PartOutMax=2706
RhopOut=True
RhopOutMin=100.000000
RhopOutMax=3000.000000
CellOrder="XYZ"
**Requested gpu memory for 3506 particles: 1.0 MB.
**CellDiv: Requested GPU memory for 3506 particles: 0.0 MB.
**CellDiv: Requested gpu memory for 16016 cells (CellMode=2H): 0.2 MB.
CellMode="2H"
Hdiv=1
MapCells=(52,1,308)
PtxasFile="../../EXECS/DualSPHysics_win64_MPhase_ptxasinfo"
**Without optimization of registers.
BsForcesBound=256 (? regs)
BsForcesFluid=256 (? regs)
BsForcesFluidCorr=256 (? regs)
RunMode="Single-Gpu"
Loading multiphase Constants to Constant memory
Finish Loading multiphase Constants to Constant memory
Allocated memory in CPU: 309420 (0.30 MB)
Allocated memory in GPU: 1365568 (1.30 MB)
Part_0000 3506 particles successfully stored
[Initialising simulation (u5h3rodd) 30-09-2015 16:29:50]
PART PartTime TotalSteps Steps Time/Seg Finish time
========= ============ ============ ======= ========= ===================
it just stop here
could you find the possible reason? 'Cause all other cases work well in my computer.
Thanks!
Lazycatyi we are looking into this
thanks for your patience
George
We found a peculiar situation related with my MVStudio configuration.
It was running in some machines but not in others.
We will update the files shortly and let you know!
Thanks for reporting the problem!
Regards,
George
The issues have been resolved.
We have updated the Package v3.2.
Please let us know if you have any more issues,
Happy multi-SPH_ing
George
Thank you for rapid update. I'm trying to use the code to deal with coastal processes with sediment. Could you email me two of your papers: "3-D SPH Modelling of Sediment Scouring Induced by Rapid Flows" and "Modelling sediment resuspension in industrial tanks using SPH"? If you have other papers that would be helpful to understand the code, could you also give me a copy?
Thank you very much! My email is panyi21hhu@gmail.com
Regards
Yi
I currently run a multi phase simulation. As I did before in v3.1 I copied my "_out"-Folder to somewhere else and created a previev (just executed the BoundaryVTK, PartVTK and Isosurface).
This works pretty good in v3.1 but in the new version all Executables crash.
I didn´t change so much in the ".bat"-File:
=========================================
@echo off
rem "name" and "dirout" are named according to the testcase
set name=Testcase1
set dirout=%name%_out
set boundaryvtk="../EXECS/BoundaryVTK_win64.exe"
rem "executables" are renamed and called from their directory
set gencase="../../EXECS/GenCase_win64.exe"
set dualsphysics="../../EXECS/DualSPHysics_win64_MPhase.exe"
set boundaryvtk="../../EXECS/BoundaryVTK_win64.exe"
set partvtk="../../EXECS/PartVTK_win64.exe"
set measuretool="../../EXECS/MeasureTool_win64.exe"
set isosurface="../../EXECS/IsoSurface_win64.exe"
rem "dirout" is created to store results or it is removed if it already exists
%boundaryvtk% -loadvtk %dirout%/%name%_Druckplatte_Actual.vtk -filexml %dirout%/%name%.xml -savevtkdata %dirout%/BoundaryMotion
%partvtk% -dirin %dirout% -savevtk %dirout%/PartFluid -onlymk:0 -vars:-all,+rhop,+vel,+mass,+press
%partvtk% -dirin %dirout% -savevtk %dirout%/PartSediment -onlymk:1 -vars:-all,+rhop,+vel,+mass,+press
if not "%ERRORLEVEL%" == "0" goto fail
%isosurface% -dirin %dirout% -filexml %dirout%/%name%.xml -onlymk:0 -saveiso %dirout%/IsosurfaceWater -onlymk:1 -saveiso %dirout%/IsosurfaceSoil
if not "%ERRORLEVEL%" == "0" goto fail
:success
echo All done
goto end
:fail
echo Execution aborted.
:end
pause
===========================================
Is it my mistake (maybe I am just to impatient
Best regards
Chris
Well if it all executables, then probably you went somewhere wrong,
where is the %dualsphysics% command line in your *.bat?
Check your script..
Cheers,
George
i removed it, because this *.bat is just for calculating the whole output after the simulation.
As is said it´s just for creating a preview. So I copy the current output-Folder (and XML, CAD-Files etc.) somewhere else and run the *.bat to create all desired outputs until the current state.
PartVTK and Isosurface are crashing both.
EDIT: PartVTK is crashing when it tries to extract the 3D Data:
PartVTK v1.41 (17-12-2013)
===========================
Threads for parallel execution: 8
LoadData> Einspuelschale_out/PartOut.bi2
LoadData> Einspuelschale_out/Part0000.bi2
LoadData> 3D data files...
Email me the *.bat file, I will take a look when I get some time, but it is very strange for your PartVTK to crash..
Cheers,
George
A few questions related to DPI and the new release:
a. is "density" and "body density" different from "PhaseRhop"? and why can't the first two be different for each fluid phase?
b. is "PhaseVisco" the fluid dynamic viscosity?
c. are the units SI? or is there no inherent set of units?
Thanks!
DPI was created well before the current multi-phase release of DualSPHysics was implemented, as such it refers only to the original single phase code. Saying that, as I understand it, there have been few changes to the form of the XML that the new GenCase accepts, with the majority being added as parameters. DPI supports free editing of non-standard parameters so it should still be useful to some regard, even if just to define a geometry.
I'll let somebody from the DualSPHysics multi-phase development team fill in the specifics though.
Also, I've sent you a message about some issues I got while accessing the "Geometry" tab on DPI. Any idea what's wrong?
Can anyone let me know about the possibility of generating say 1000 irregular waves in DualSPHysics. If it is possible, how long would be the computational time?
I want to do numerical simulation on wave runup and wave overtopping on a curved dike.
Thanks,
Suba
a 2D case can take hours but 3D case.... 1000 waves can take 2-3 days (on GPU)
Thanks for your reply. I would really wait for the new version.
Before that I would like to pre equip myself in DualSPHysics.
Is there is any pre requisite before I start exploring the new version?
Can you suggest me something in this regard.
Thanks,
Suba