A Stabilized Incompressible SPH Method by Relaxing the Density Invariance Condition

<table>Detect free surface numerically using the number of neighbouring particles for using <svg height="12.4375" id="M109" style="vertical-align:-0.1638pt" version="1.1" viewbox="0 0 36.262501 12.4375" width="36.262501" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.162)"><path d="M480 416q0 -21 -18 -41q-9 -11 -17 -7q-20 9 -42 9q-62 0 -140 -78q23 -69 88 -192q17 -31 27 -42t20 -11q16 0 62 46l17 -20q-64 -92 -119 -92q-35 0 -70 66q-41 73 -84 187q-36 -30 -62 -61q-27 -115 -35 -172q-41 -8 -78 -20l-6 6l140 612q7 28 0.5 34t-37.5 7l-34 1
l5 26q38 4 74 13.5t57 17t25 7.5q12 0 4 -32l-104 -443h2q35 38 97 93q39 35 65.5 56t62 41.5t58.5 20.5q19 0 30.5 -10t11.5 -22z" id="x1D458"></path></g><g transform="matrix(.017,-0,0,-.017,13.322,12.162)"><use xlink:href="#x3D"></use></g><g transform="matrix(.017,-0,0,-.017,28.025,12.162)"><path d="M285 378v-2q65 -13 102 -54.5t37 -97.5q0 -57 -30.5 -104.5t-74 -75t-85.5 -42t-72 -14.5q-31 0 -59.5 11t-40.5 23q-19 18 -16 36q1 16 23 33q13 10 24 0q58 -51 124 -51q55 0 88 40t33 112q0 64 -39 96.5t-88 32.5q-29 0 -64 -11l-6 29q77 25 118 57.5t41 84.5
q0 45 -26.5 69.5t-68.5 24.5q-67 0 -120 -79l-20 20l43 63q51 56 127 56h1q66 0 107 -37t41 -95q0 -42 -31 -71q-22 -23 -68 -54z" id="x33"></path></g>
</svg> and <svg height="16.450001" id="M110" style="vertical-align:-3.39066pt" version="1.1" viewbox="0 0 67.112503 16.450001" width="67.112503" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.162)"><path d="M499 88q-41 -46 -83.5 -73t-68.5 -27q-41 0 -19 98l52 225q10 43 8 59.5t-16 16.5q-37 0 -112 -65t-117 -135q-14 -42 -38 -179q-36 -7 -75 -20l-7 8l127 594q8 40 2.5 49t-35.5 9h-35l5 25q36 3 71.5 13t56.5 18t25 8q11 0 3 -38l-88 -417h2q77 102 195 167q44 24 79 24
q61 0 30 -130l-52 -220q-8 -33 4 -33q8 0 31 14.5t43 33.5z" id="x210E"></path></g><g transform="matrix(.017,-0,0,-.017,13.543,12.162)"><use xlink:href="#x3D"></use></g><g transform="matrix(.017,-0,0,-.017,28.246,12.162)"><path d="M384 0h-275v27q67 5 81.5 18.5t14.5 68.5v385q0 38 -7.5 47.5t-40.5 10.5l-48 2v24q85 15 178 52v-521q0 -55 14.5 -68.5t82.5 -18.5v-27z" id="x31"></path></g><g transform="matrix(.017,-0,0,-.017,36.406,12.162)"><use xlink:href="#x2E"></use></g><g transform="matrix(.017,-0,0,-.017,40.282,12.162)"><use xlink:href="#x32"></use></g><g transform="matrix(.017,-0,0,-.017,51.28,12.162)"><use xlink:href="#x1D451"></use></g>
<g transform="matrix(.012,-0,0,-.012,60.712,16.25)"><use xlink:href="#x30"></use></g>
</svg> for quintic spline kernel function at initial particle size <svg height="16.450001" id="M111" style="vertical-align:-3.39066pt" version="1.1" viewbox="0 0 87.737503 16.450001" width="87.737503" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink">
<g transform="matrix(.017,-0,0,-.017,.062,12.162)"><use xlink:href="#x1D451"></use></g>
<g transform="matrix(.012,-0,0,-.012,9.5,16.25)"><use xlink:href="#x30"></use></g>
<g transform="matrix(.017,-0,0,-.017,20.55,12.162)"><use xlink:href="#x3D"></use></g><g transform="matrix(.017,-0,0,-.017,35.237,12.162)"><use xlink:href="#x30"></use></g><g transform="matrix(.017,-0,0,-.017,43.396,12.162)"><use xlink:href="#x2E"></use></g><g transform="matrix(.017,-0,0,-.017,47.272,12.162)"><use xlink:href="#x30"></use></g><g transform="matrix(.017,-0,0,-.017,55.431,12.162)"><use xlink:href="#x30"></use></g><g transform="matrix(.017,-0,0,-.017,63.591,12.162)"><use xlink:href="#x35"></use></g><g transform="matrix(.017,-0,0,-.017,73.739,12.162)"><path d="M797 0h-214v26q46 5 58 17.5t12 63.5v167q0 116 -88 116q-48 0 -112 -52q3 -18 3 -42v-195q0 -46 11 -58t53 -17v-26h-214v26q48 5 59.5 17t11.5 60v172q0 115 -88 115q-47 0 -109 -52v-235q0 -47 10 -59.5t48 -17.5v-26h-211v26q50 5 62.5 17t12.5 60v207q0 48 -9 59.5
t-52 18.5v23q79 14 139 41v-81q3 2 29 20.5t44 29.5q45 28 81 28t66 -22t46 -61q25 20 83 57q39 26 76 26q59 0 93 -40t34 -113v-193q0 -49 10.5 -60.5t54.5 -16.5v-26z" id="x6D"></path></g>
</svg>.</table>

Journal of Applied Mathematics

fig5

Figure 5

Figure 5: A Stabilized Incompressible SPH Method by Relaxing the Density Invariance Condition 