Error Simulating

[SohailAusgur]

Hi,

I am having trouble running simulations with the DesignSPH macros. It come up as seen in the picture attached below. I tried it on my laptop and my PC both of which are powerful enough to be running simulations.

I update my graphics drivers aswell as it stated in the Error.

I checked all the parameters and they are all fine as i am doing this simulation from your Youtube channel. This is the simulation I am running: https://www.youtube.com/watch?v=8DKaDAHAFJk&t=312s

My parameters are the exact same as the one in the video so i do not know what is going on here.

Please advise on how to fix this Many Thanks

Error detailed:

Copyright (c) 2023 by Dr Jose M. Dominguez Alonso, Dr Alejandro Crespo, Prof. Moncho Gomez Gesteira, Prof. Benedict Rogers, Dr Georgios Fourtakas, Prof. Peter Stansby, Dr Renato Vacondio, Dr Corrado Altomare, Dr Angelo Tafuni, Dr Orlando Garcia Feal, Ivan Martinez Estevez, Dr Joseph O'Connor, Dr Aaron English 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 Lesser General Public License as published by the Free Software Foundation, either version 2.1 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 Lesser General Public License for more details. You should have received a copy of the GNU Lesser General Public License along with DualSPHysics. If not, see . Free download source files and documentation from DualSPHysics website (http://dual.sphysics.org). Full list of developers and contributors is available at https://dual.sphysics.org/developers Properly cite the use of DualSPHysics (Dominguez et al., 2022) (https://dual.sphysics.org/references). DualSPHysics5 v5.2.272 (12-06-2023) ==================================== Threads by host for parallel execution: 12 [Initialising JSphCpuSingle 25-02-2024 17:21:53] ProgramFile=".../dualsphysics/bin/DualSPHysics5.2_win64.exe" ExecutionDir=".../FreeCAD/Mod/DesignSPHysics" XmlFile=".../casewaves1/casewaves1_out/casewaves1.xml" OutputDir=".../dualsphcases/casewaves1/casewaves1_out" OutputDataDir=".../dualsphcases/casewaves1/casewaves1_out" XML-App: GenCase v5.0.278c (30-10-2023) XML-Vars (uservars + ctes): MResId=[0] CaseName=["casewaves1"] Data2D=[1] Data2DPosy=[0] H=[0.0254558] KernelSize=[0.0509117] B=[210214] Gamma=[7] Rhop0=[1000] Dp=[0.015] Gravity_x=[0] Gravity_y=[0] Gravity_z=[-9.81] MassFluid=[0.225] MassBound=[0.225] XML-Vars (parameters): TimeMax=[8] TimeOut=[0.025] *** WARNING: The use of DEM was disabled because there are no floating objects... **Basic case configuration is loaded **Special case configuration is loaded Loading initial state of particles... Loaded particles: 11111 MapRealPos(border)=(-0.106273,-0.00127279,-0.00127279)-(10.9963,0.00127279,0.706273) MapRealPos(final)=(-0.106273,-0.00127279,-0.00127279)-(10.9963,0.00127279,0.706273) **Initial state of particles is loaded **2D-Simulation parameters: CaseName="casewaves1" RunName="casewaves1" Simulate2DPosY=0 Symmetry=False SavePosDouble=False SvExtraParts="" SaveFtAce=False SvTimers=True Boundary="DBC" StepAlgorithm="Symplectic" Kernel="Wendland" Wendland.awen=859.5679 Wendland.bwen=-168847.2 Viscosity="Artificial" Visco=0.01 ViscoBoundFactor=0 DensityDiffusion="None" Shifting="None" RigidAlgorithm="None" FloatingCount=0 CaseNp=11111 CaseNbound=885 CaseNfixed=838 CaseNmoving=47 CaseNfloat=0 CaseNfluid=10226 PeriodicActive="None" Dp=0.015 KernelH=0.025456 (CoefficientH=1.2; H/Dp=1.69706) KernelSize=0.05091169 CteB=210214.3 Gamma=7 RhopZero=1000 Cs0=38.36013515030937 CFLnumber=0.2 DtIni=0.0006636015112436296 DtMin=3.318007605660313e-05 DtAllParticles=False MassFluid=0.225 MassBound=0.225 TimeMax=8 TimePart=0.025 Gravity=(0,0,-9.81) NpMinimum=885 RhopOut=True RhopOutMin=700 RhopOutMax=1300 WrnPartsOut=True CellMode="Full" ScellDiv=1 MapCells=(219,1,14) CellDomFixed=False **Requested CPU memory for 11239 particles: 1.7 MB. DomCells=(219,1,14) DomCellCode="1+14_8_9" **CellDiv: Requested cpu memory for 11794 particles: 0.4 MB. **CellDiv: Requested cpu memory for 3120 cells (CellMode=Full): 0.0 MB. RunMode="Pos-Double - OpenMP(Threads:12)" Particle summary: Fixed....: 838 id:(0-837) MKs:2 (10-11) Moving...: 47 id:(838-884) MKs:1 (12) Floating.: 0 Fluid....: 10226 id:(885-11110) MKs:1 (1) Total particles: 11111 (bound=885 (fx=838 mv=47 ft=0) fluid=10226) Total MK blocks: 4 (bound=3 (fx=2 mv=1 ft=0) fluid=1) Wave paddles configuration: **Saving motion paddle in .../casewaves1/casewaves1_out/WavePaddle_mkb0002.csv Paddle_0 (mkb:2) Type of paddle: Piston - Regular Movement direction: (1:0:0) [m] Start time: 0 [s] (end: 8 [s]) Gravity: 9.81 [m/s^2] Depth: 0.4 [m] InitialPhase: 0 times PI InitialRamp: 0 periods WaveOrder: 2nd WaveHeight: 0.5 [m] WavePeriod: 1.2 [s] WaveLength: 1.93625 [m] Relative depth (d/L): 0.206585 (Transitional water) GainStroke: 1 Stroke: 0.403308 Amplitude: 0.201654 [m] *** Warning: Waves breaking for achieved steepness (Miche 1944). *** Warning: Waves breaking for depth limitation (McCowan 1893). *** Warning: Madsen criterion is not satisfied (Madsen 1971). Particles: 838 - 884 (47) Allocated memory in CPU: 2226924 (2.12 MB) Part_0000 11111 particles successfully stored [WARNINGS] 1. The use of DEM was disabled because there are no floating objects... 2. Wave generation> Waves breaking for achieved steepness (Miche 1944). 3. Wave generation> Waves breaking for depth limitation (McCowan 1893). 4. Wave generation> Madsen criterion is not satisfied (Madsen 1971). [Initialising simulation (z29rvrds) 25-02-2024 17:21:53] PART PartTime TotalSteps Steps Time/Sec Finish time ========= ============ ============ ======= ========= =================== Part_0001 0.025025 202 202 38.95 25-02-2024 17:27:04 Part_0002 0.050066 416 214 39.86 25-02-2024 17:27:08 Part_0003 0.075052 627 211 40.14 25-02-2024 17:27:10 Part_0004 0.100054 836 209 40.40 25-02-2024 17:27:11 Part_0005 0.125040 1044 208 39.93 25-02-2024 17:27:11 Part_0006 0.150035 1245 201 38.55 25-02-2024 17:27:10 Part_0007 0.175041 1448 203 39.48 25-02-2024 17:27:10 Part_0008 0.200048 1645 197 37.93 25-02-2024 17:27:08 Part_0009 0.225057 1843 198 38.42 25-02-2024 17:27:07 Part_0010 0.250098 2038 195 37.88 25-02-2024 17:27:05 Part_0011 0.275033 2231 193 37.42 25-02-2024 17:27:04 Part_0012 0.300039 2424 193 38.39 25-02-2024 17:27:03 Part_0013 0.325043 2616 192 36.87 25-02-2024 17:27:02 Part_0014 0.350081 2808 192 38.37 25-02-2024 17:27:02 Part_0015 0.375063 3002 194 40.71 25-02-2024 17:27:03 Part_0016 0.400076 3195 193 38.13 25-02-2024 17:27:03 Part_0017 0.425060 3389 194 38.26 25-02-2024 17:27:02 Part_0018 0.450003 3585 196 38.78 25-02-2024 17:27:02 Part_0019 0.475121 3780 195 39.09 25-02-2024 17:27:03 Part_0020 0.500015 3975 195 38.68 25-02-2024 17:27:03 Part_0021 0.525046 4171 196 41.82 25-02-2024 17:27:04 Part_0022 0.550042 4367 196 41.49 25-02-2024 17:27:05 Part_0023 0.575030 4563 196 40.54 25-02-2024 17:27:05 Part_0024 0.600091 4759 196 38.47 25-02-2024 17:27:05 Part_0025 0.625099 4958 199 39.61 25-02-2024 17:27:05 Part_0026 0.650017 5155 197 38.98 25-02-2024 17:27:05 Part_0027 0.675031 5354 199 40.32 25-02-2024 17:27:05 Part_0028 0.700098 5556 202 39.73 25-02-2024 17:27:06 *** ERROR: Some boundary particle was excluded. *** TimeStep: 0.714696 (Nstep: 5671) Total boundary: 47 (fixed=0 moving=47 floating=0) Excluded for: position=47 rhop=0 velocity=0 Some boundary particle exceeded the -X limit (left limit) of the simulation domain. *** Exception (JSphCpuSingle::AbortBoundOut) at DualSPHysics/source/JSph.cpp:2549 Text: Fixed, moving or floating particles were excluded. Check VTK file Error_BoundaryOut.vtk with excluded particles. Finished execution (code=1).

[imestevez]

The error, as shown at the end of the details, says the following:

*** ERROR: Some boundary particle was excluded. ***
TimeStep: 0.714696 (Nstep: 5671)
Total boundary: 47 (fixed=0 moving=47 floating=0)
Excluded for: position=47 rhop=0 velocity=0
Some boundary particle exceeded the -X limit (left limit) of the simulation domain.

This means that the piston (moving object) exceeded the -X limit of the domain when generating waves. To fix that, increase the left side of the CaseLimits.

[SohailAusgur]

Hi,

I increased the Case limits about 10000mm in the -X direction and the same Issue Occurs.

---

From: Iván Martínez Estévez @.> Sent: 25 February 2024 17:47 To: DualSPHysics/DesignSPHysics @.> Cc: C19442692 Mohammad Ausgur @.>; Author @.> Subject: Re: [DualSPHysics/DesignSPHysics] Error Simulating (Issue #178)

The error, as shown at the end of the details, says the following:

ERROR: Some boundary particle was excluded. TimeStep: 0.714696 (Nstep: 5671) Total boundary: 47 (fixed=0 moving=47 floating=0) Excluded for: position=47 rhop=0 velocity=0 Some boundary particle exceeded the -X limit (left limit) of the simulation domain.

This means that the piston (moving object) exceeded the -X limit of the domain when generating waves. To fix that, increase the left side of the CaseLimits.

— Reply to this email directly, view it on GitHubhttps://github.com/DualSPHysics/DesignSPHysics/issues/178#issuecomment-1963011347, or unsubscribehttps://github.com/notifications/unsubscribe-auth/BGOE7GOORTS7X34KJC35BK3YVN2KBAVCNFSM6AAAAABDZAXPTSVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTSNRTGAYTCMZUG4. You are receiving this because you authored the thread.Message ID: @.***>

[imestevez]

The final domain size is adjusted to the particles created within the limits to optimise memory. So in your case, the modification in CaseLimits does not have effect. Therefore, you can try: 1) Increase the left bottom size a little until the piston movement does not exceed the bottom limit. or 2) Go to Execution Parameters->Simulation Doomain -> Increase the minimum X size by 10 % (If 10 % is not enough, then increase it a bit more). This option enforces the increase of the domain.

[SohailAusgur]

Hi,

I am increasing the Domain by 10% as advised. However, I am now getting another Error when running the GenCase:

See detail below:

Copyright (c) 2020 by Dr Jose M. Dominguez All rights reserved. GenCase is a free tool for DualSPHysics simulations. DualSPHysics project is an international collaboration between University of Vigo (Spain) and University of Manchester (UK). Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met: * Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer. * Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution. * Neither the name of the GenCase nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission. THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. GenCase v5.0.278c (30-10-2023) =============================== LoadDsphConfig> .../dualsphysics/bin/DsphConfig.xml OmpThreads: 12 [LoadXMLInit] List: MAIN Total operations: 18 [LoadXMLInit-End time:0.000s] [Draw] List of available variables: Gravity_x=[0] Gravity_y=[0] Gravity_z=[-9.81] Rhop0=[1000] MResId=[0] Dp=[0.015] PosMin_x=[-0.105] PosMin_y=[0] PosMin_z=[0] PosMax_x=[11.4] PosMax_y=[0] PosMax_z=[1] CaseName=["casewaves"] Data2D=[1] Data2DPosy=[0] H=[0.0254558] MassBound=[0.225] MassFluid=[0.225] Mk size: 1 bytes Case limits: (-0.105,0,0) - (11.4,0,1) Distance between points (Dp): 0.015 Reference position: (0,0,0) Case domain: (0,0,0)-(766,0,65) Domain points: 767 x 1 x 66 = 50622 Case parts: 1 x 66z Memory: 0.00Gb 0.05Mb (50622) > MAIN FileDbgDomains> C:/Users/sohau/Documents/dualsphcases/casewaves/casewaves_out/casewaves_dbg-fillbox.vtk domains:1 FileShapes> _Dp.vtk shapes:18 points:12 (time:0.001s) FileShapes> _Actual.vtk shapes:18 points:24 (time:0.001s)

[imestevez]

Try typing - 10.0% in that field

[SohailAusgur]

Hi,

I extended my Bottom cube in the -x direction and then increases the case Limits. This seemed to do the trick. However the -10.0% did not work.

Many Thanks Anyways.

Much Appreciated.

[TShorlar]

Okay, thank you.