{"id":436644,"date":"2025-07-11T22:14:53","date_gmt":"2025-07-11T22:14:53","guid":{"rendered":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/issues-with-modifying-official-ls-dyna-ex-with-s-ale-mesh-2\/"},"modified":"2025-07-14T09:25:13","modified_gmt":"2025-07-14T09:25:13","slug":"issues-with-modifying-official-ls-dyna-ex-with-s-ale-mesh-2","status":"closed","type":"topic","link":"https:\/\/innovationspace.ansys.com\/forum\/forums\/topic\/issues-with-modifying-official-ls-dyna-ex-with-s-ale-mesh-2\/","title":{"rendered":"Issues with modifying official LS-DYNA ex with S-ALE mesh"},"content":{"rendered":"<p>&lt;p&gt;Hello, I would greatly appreciate some help with this. I am modifying an official LS-DYNA example, Underwater F (https:\/\/lsdyna.ansys.com\/underwater-f\/), to instead have a square field that is constructed with an S-ALE mesh. Everything else stays the same, so the resulting simulations should look identical for the first ~30 time steps, until the shock wave reaches the boundary of the cutout in the example. However, they look different and the shock wave is not propagating&nbsp;outward into the water and it looks like there is a boundary effect.&lt;\/p&gt;&lt;p&gt;I created an S-ALE mesh and a box region at the bottom left corner (the origin). I filled the entire region with water and then the box region with the explosive material. I detonated the material at the origin at time 0.&lt;\/p&gt;&lt;p&gt;1. There seems to be some kind of boundary effect at the box region. I don&rsquo;t think I defined any boundaries and am wondering if there is a discontinuity or something between the two materials. However, I don&rsquo;t know how to fix this, particularly with this S-ALE method.&lt;\/p&gt;&lt;p&gt;2. Originally, I kept the element sizes the same, 1 element = 1 cm. However, for my simulation, the velocity contour did not show the shock wave, which was very strange, so I refined the mesh to have each element be 0.16 cm. Now I can see the shock wave, but it is strange that this refinement was not necessary for the official example but was for mine. Would anyone know why this might be?&lt;\/p&gt;&lt;p class=&#8221;MsoNormal&#8221;&gt;Please, I would greatly appreciate some help, thank you very much in advance.&lt;\/p&gt;&lt;p class=&#8221;MsoNormal&#8221;&gt;&nbsp;&lt;\/p&gt;&lt;p class=&#8221;MsoNormal&#8221;&gt;Here is my code:&nbsp;&lt;\/p&gt;&lt;div&gt;*KEYWORD&lt;\/div&gt;&lt;div&gt;*NODE&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp;1 0.000000000E+00 0.000000000E+00 0.000000000E+00&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp;2 0.000000000E+00 0.000000000E+00 0.000000000E+00&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp;3 0.100000000E+00 0.000000000E+00 0.000000000E+00&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp;4 0.000000000E+00 0.100000000E+00 0.000000000E+00&lt;\/div&gt;&lt;div&gt;*ALE_STRUCTURED_MESH_CONTROL_POINTS&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; 1001&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp;0.0&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 51&nbsp; &nbsp; &nbsp; &nbsp;8.0&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;201&nbsp; &nbsp; &nbsp; 32.0&lt;\/div&gt;&lt;div&gt;*ALE_STRUCTURED_MESH_CONTROL_POINTS&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; 1002&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp;0.0&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 51&nbsp; &nbsp; &nbsp; &nbsp;8.0&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;201&nbsp; &nbsp; &nbsp; 32.0&lt;\/div&gt;&lt;div&gt;*ALE_STRUCTURED_MESH_CONTROL_POINTS&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; 1003&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp;0.0&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&nbsp; &nbsp; &nbsp; &nbsp;1.0&lt;\/div&gt;&lt;div&gt;*DEFINE_COORDINATE_NODES&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp;234&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;3&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;4&lt;\/div&gt;&lt;div&gt;*ALE_STRUCTURED_MESH&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;7&nbsp; &nbsp; &nbsp; 5000&nbsp; &nbsp; &nbsp;10000&nbsp; &nbsp; &nbsp;10000&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; 1001&nbsp; &nbsp; &nbsp; 1002&nbsp; &nbsp; &nbsp; 1003&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp;234&nbsp; &nbsp; &nbsp;&lt;\/div&gt;&lt;div&gt;*DEFINE_BOX&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;9&nbsp; &nbsp; &nbsp; &nbsp;0.0&nbsp; &nbsp; &nbsp; &nbsp;8.0&nbsp; &nbsp; &nbsp; &nbsp;0.0&nbsp; &nbsp; &nbsp; &nbsp;8.0&nbsp; &nbsp; &nbsp; &nbsp;0.0&nbsp; &nbsp; &nbsp; &nbsp;1.0&lt;\/div&gt;&lt;div&gt;*PART&lt;\/div&gt;&lt;div&gt;&nbsp;&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&lt;\/div&gt;&lt;div&gt;*PART&lt;\/div&gt;&lt;div&gt;&nbsp;&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&lt;\/div&gt;&lt;div&gt;*MAT_NULL&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp;1.0&lt;\/div&gt;&lt;div&gt;*MAT_HIGH_EXPLOSIVE_BURN&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&nbsp; &nbsp; &nbsp; 1.63&nbsp; &nbsp; &nbsp;0.784&nbsp; &nbsp; &nbsp;0.260 0.000E+00&lt;\/div&gt;&lt;div&gt;*EOS_GRUNEISEN&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp;1.5&lt;\/div&gt;&lt;div&gt;&nbsp;&lt;\/div&gt;&lt;div&gt;*EOS_JWL&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&nbsp; &nbsp; &nbsp; 3.71 3.230E-02&nbsp; &nbsp; &nbsp; 4.15&nbsp; &nbsp; &nbsp;0.950&nbsp; &nbsp; &nbsp;0.300 4.300E-02&nbsp; &nbsp; &nbsp; 1.00&lt;\/div&gt;&lt;div&gt;*ALE_STRUCTURED_MULTI-MATERIAL_GROUP&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp;water&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; expl&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&lt;\/div&gt;&lt;div&gt;*ALE_STRUCTURED_MESH_VOLUME_FILLING&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;7&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;water&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp;ALL&lt;\/div&gt;&lt;div&gt;*ALE_STRUCTURED_MESH_VOLUME_FILLING&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;7&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; expl&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; BOXCOR&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;9&lt;\/div&gt;&lt;div&gt;*CONTROL_ALE&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&nbsp; &nbsp; &nbsp; -1.0&lt;\/div&gt;&lt;div&gt;&nbsp;&lt;\/div&gt;&lt;div&gt;*CONTROL_ENERGY&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;1&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;2&lt;\/div&gt;&lt;div&gt;*INITIAL_DETONATION&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; 5000 0.000E+00 0.000E+00 0.000E+00 0.000E+00&lt;\/div&gt;&lt;div&gt;*CONTROL_TERMINATION&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp;100.0&lt;\/div&gt;&lt;div&gt;*DATABASE_BINARY_D3PLOT&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp;1.0&lt;\/div&gt;&lt;div&gt;*DATABASE_GLSTAT&lt;\/div&gt;&lt;div&gt;&nbsp; &nbsp; &nbsp; &nbsp;0.1&lt;\/div&gt;&lt;div&gt;*END&lt;\/div&gt;&lt;div&gt;&nbsp;&lt;\/div&gt;&lt;div&gt;And here are the screenshots of the pressure contours:&lt;\/div&gt;&lt;div&gt;<img decoding=\"async\" src=\"https:\/\/innovationspace.ansys.com\/forum\/wp-content\/uploads\/sites\/2\/2025\/07\/11-07-2025-1752272045-1.jpg\" alt=\"\" \/>&lt;\/div&gt;&lt;div&gt;<img decoding=\"async\" src=\"https:\/\/innovationspace.ansys.com\/forum\/wp-content\/uploads\/sites\/2\/2025\/07\/11-07-2025-1752272055-2.jpg\" alt=\"\" \/>&lt;\/div&gt;&lt;div&gt;<img decoding=\"async\" src=\"https:\/\/innovationspace.ansys.com\/forum\/wp-content\/uploads\/sites\/2\/2025\/07\/11-07-2025-1752272065-3.jpg\" alt=\"\" \/>&lt;\/div&gt;&lt;div&gt;<img decoding=\"async\" src=\"https:\/\/innovationspace.ansys.com\/forum\/wp-content\/uploads\/sites\/2\/2025\/07\/11-07-2025-1752272075-4.jpg\" alt=\"\" \/>&lt;\/div&gt;<\/p>\n","protected":false},"template":"","class_list":["post-436644","topic","type-topic","status-closed","hentry","topic-tag-ale","topic-tag-Explosion-1","topic-tag-lsdyna","topic-tag-underwater-shockwave-1"],"aioseo_notices":[],"acf":[],"custom_fields":[{"0":{"_bbp_forum_id":["27792"],"_bbp_topic_id":["436644"],"_bbp_subscription":["572157","199"],"_bbp_author_ip":["71.179.15.164"],"_bbp_last_reply_id":["436699"],"_bbp_last_active_id":["436699"],"_bbp_last_active_time":["2025-07-14 09:25:05"],"_bbp_reply_count":["1"],"_bbp_reply_count_hidden":["0"],"_bbp_voice_count":["2"],"_bbp_engagement":["572157","199"],"_btv_view_count":["157"],"_bbp_topic_status":["unanswered"],"_bbp_notification_enabled":["572157"],"_bbp_status":["publish"]},"test":"e70e8cc42ca9c70bd83b8f98ebf17389fbd1db36"}],"_links":{"self":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/436644","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics"}],"about":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/types\/topic"}],"version-history":[{"count":0,"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/topics\/436644\/revisions"}],"wp:attachment":[{"href":"https:\/\/innovationspace.ansys.com\/forum\/wp-json\/wp\/v2\/media?parent=436644"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}