2D - FSI Leakage with Erosion Problem

79 views
Skip to first unread message

Furkan Keleş

unread,
Dec 16, 2021, 12:31:47 PM12/16/21
to LS-DYNA2
Dear all,

I'm conducting a 2D-Axisymmetrical FSI simulation involving explosives. During the simulation, some regions of the lagrange parts are subjected to very high FSI loads.

Usually the simulation runs fine with no leakage problems. However, in some cases where the FSI loads are higher, lagrange elements coupled to fluids erode and that causes leakage. I have no leakage problem as long as no erosion occurs.

I'm using Constrained_Lagrange_in_Solid keyword with CTYPE=5 and Johnson-Cook damage parameters for lagrange parts.

How can I solve this leakage issue without removing erosion in lagrange parts? 

Thank you,

Furkan.

l...@schwer.net

unread,
Dec 16, 2021, 1:27:13 PM12/16/21
to Furkan Keleş, LS-DYNA2

Consider a further reduction in the time step via TSSFAC on Control Timestep keyword. CTYPE=5 is supposed to handle eroded elements and FSI, but maybe it needs a little help?            --len

--
You received this message because you are subscribed to the Google Groups "LS-DYNA2" group.
To unsubscribe from this group and stop receiving emails from it, send an email to ls-dyna2+u...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/ls-dyna2/f2d3eed3-25f5-4f5e-b071-557d8873cfc8n%40googlegroups.com.

James M. Kennedy

unread,
Dec 16, 2021, 5:20:21 PM12/16/21
to Furkan Keleş, LS-DYNA2

Dear Furkan,

 

See if these presentations are of some help:

 

https://ftp.lstc.com/anonymous/outgoing/support/FAQ/2d_general_condensed

 

A 2D Multi-Material Arbitrary Lagrange Euler (MM-ALE) code was implemented in LS-DYNA. Like the 3D MM-ALE, already available, each 2D computational cycle is divided into two steps. First a multi-material version of the two-dimensional shell formulations solves the physical problem on quadrangle meshes during the Lagrangian step. The 2D shell formulations are plane strain and are-weighted axisymmetric. An advection cycle adapted to the 2D shell approaches follows to control the mesh motion. 2D ALE data of the last cycle can be mapped on 3D ALE mesh:

 

Aquelet, N., and Souli, M., "2D to 3D ALE Mapping", 10th International LS-DYNA User's Conference, Dearborn, Michigan, June, 2008.

 

http://www.dynalook.com/international-conf-2008/FluidStructure-3.pdf

 

Pure MMALE was tested on impacts and explosives studies from defense and spatial fields. A high velocity impact, a long rod penetration, and explosively formed projectile, a shaped charge jet and an air blast were modeled using 2D axisymmetric models and compared with experimental data:

 

Van Dorsselaer, N., and Lapoujade, V., "A Contribution to New ALE 2D Method Validation", 11th International LS-DYNA Users Conference, Dearborn, Michigan, June, 2010.

 

http://www.dynalook.com/international-conf-2010/Simulation-2-5.pdf

 

A numerical simulation was conducted to model explosive detonation and blast wave propagation in the open air field. The mesh size and boundary conditions as well as size

of air domain were the sensitive variables which may significantly affect the predicted pressure wave magnitude and rising time in blast simulations. The current approach focused on determining the optimal key parameters to predict the blast wave accurately. A 2D to 3D mapping was performed to save computational time:

 

Kalra, A., Zhu, F., Yang, K.H., and King, A.I., "Key Parameters in Blast Loading Using 2D to 3D ALE Mapping Technique", 13th International LS-DYNA Users Conference, Dearborn, Michigan, June, 2014.

 

http://www.dynalook.com/13th-international-ls-dyna-conference/blast/key-parameters-in-blast-modeling-using-2d-to-3d-ale-mapping-technique

 

LS-DYNA incorporated modeling of detonation process via the arbitrary-Lagrangian-Eulerian (ALE) elements or by using *LOAD_BLAST keyword. The pressure loading developed due to detonation of an explosive can be coupled with a Lagrangian part to model the fluid structure interaction (FSI). This was done with the *CONSTRAINED_ LAGRANGE_IN_SOLID (*CLIS) keyword using penalty type coupling (CTYPE 5). The concrete slab was modeled with *MAT_159, the continuous surface cap model (CSCM) constitutive law. Preliminary investigations were carried out for the blast pressures developed by the detonation process in LS-DYNA. The keyword parameters and mesh sensitivity were examined both in 2D and 3D modeling:

 

Dua, A., Braimah, B., and Matsagar, V., "Understanding the Response of Reinforced Concrete Slabs due to Contact Explosion of TNT", 6th International Conference on Engineering Mechanics and Materials (CSCE2017), Vancouver, British Columbia, Canada, May-June, 2017.

 

https://www.xcdsystem.com/csce/proceedings2017/EMM/FinalPaper_522.pdf

 

------------------------------------------

 

For further ALE help in 2D regarding the details, Nicolas Aquelet has generously shared the following.

 

Here are a few input decks and power point presentations, plus the KEYWORD description entries. The first input deck (slam.tar.gz) is an example of coupling. It is a wedge of beams entering water at a constant velocity.

 

input decks:

ftp://ftp.lstc.com/outgoing/aquelet/ale2d/slam.tar.gz

ftp://ftp.lstc.com/outgoing/aquelet/ale2d/explo.tar.gz

ftp://ftp.lstc.com/outgoing/aquelet/ale2d/airblastreflection.tar.gz

ftp://ftp.lstc.com/outgoing/aquelet/ale2d/underwaterexplo.bounmapping.tar.gz

 

power points (the 2008 LS-DYNA Conference presentation plus a set of class type like notes that are similar to the Conference presentation):

ftp://ftp.lstc.com/outgoing/aquelet/ale2d/dyna08.zip

ftp://ftp.lstc.com/outgoing/aquelet/ale2d/class09.zip

 

KEYWORD descriptions (most of the mapping applications only require the first 2 entries):

ftp://ftp.lstc.com/outgoing/aquelet/ale2d/SECTION_ALE2D.doc

ftp://ftp.lstc.com/outgoing/aquelet/ale2d/INITIAL_ALE_MAPPING.doc

ftp://ftp.lstc.com/outgoing/aquelet/ale2d/BOUNDARY_ALE_MAPPING.doc

 

------------------------------------------

 

Sincerely,

James M. Kennedy

KBS2 Inc.

December 16, 2021

 

From: ls-d...@googlegroups.com [mailto:ls-d...@googlegroups.com] On Behalf Of Furkan Keles
Sent: Thursday, December 16, 2021 11:32 AM
To: LS-DYNA2 <ls-d...@googlegroups.com>
Subject: [LS-DYNA2] 2D - FSI Leakage with Erosion Problem

 

Dear all,

--

James M. Kennedy

unread,
Dec 16, 2021, 5:46:04 PM12/16/21
to Furkan Keleş, LS-DYNA2

p.s. Perhaps some better links for the following presentation:

 

Dua, A., Braimah, B., and Matsagar, V., "Understanding the Response of Reinforced Concrete Slabs due to Contact Explosion of TNT", 6th International Conference on Engineering Mechanics and Materials (CSCE2017), Vancouver, British Columbia, Canada, May-June, 2017

 

James M. Kennedy

unread,
Dec 16, 2021, 6:52:16 PM12/16/21
to Furkan Keleş, LS-DYNA2

Dear Furkan,

 

Perhaps this presentation would be of some interest:

 

During the past two years, the author has worked on a new fluid-structure interaction (FSI)

package dedicated to be used with S-ALE SOLVER. The objective is to shorten the running

time, stop leakage, and make the input deck user friendly. In this paper this new FSI package

with its keyword – *ALE_STRUCTURED_FSI is introduced.

 

Chen, H., "ALE_STRUCTURED_FSI The New S-ALE Solver", 16th International LS-DYNA

Users Conference, Dearborn, Michigan, June, 2020.

 

https://www.dynalook.com/conferences/16th-international-ls-dyna-conference/fsi-ale-computing-technology-t8-1/t8-1-b-fsi-ale-081.pdf

 

Sincerely,

James M. Kennedy

KBS2 Inc.

December 16, 2021

 

From: James M. Kennedy [mailto:j...@kbs2.com]
Sent: Thursday, December 16, 2021 4:46 PM
To: 'Furkan Keleş' <kmuc...@gmail.com>; 'LS-DYNA2' <ls-d...@googlegroups.com>
Subject: RE: [LS-DYNA2] 2D - FSI Leakage with Erosion Problem

 

p.s. Perhaps some better links for the following presentation:

 

Dua, A., Braimah, B., and Matsagar, V., "Understanding the Response of Reinforced Concrete Slabs due to Contact Explosion of TNT", 6th International Conference on Engineering Mechanics and Materials (CSCE2017), Vancouver, British Columbia, Canada, May-June, 2017

 

https://docplayer.net/203103353-Understanding-the-response-of-reinforced-concrete-slabs-due-to-contact-explosion-of-tnt.html

 

https://ir.library.carleton.ca/pub/24336

 

https://www.csce.ca/elf/apps/CONFERENCEVIEWER/conferences/2018/pdfs/Paper_DM13_0608011503.pdf

 

 

From: ls-d...@googlegroups.com [mailto:ls-d...@googlegroups.com] On Behalf Of James M. Kennedy
Sent: Thursday, December 16, 2021 4:20 PM
To: 'Furkan Keleş' <
kmuc...@gmail.com>; 'LS-DYNA2' <ls-d...@googlegroups.com

>

Furkan Keleş

unread,
Dec 17, 2021, 1:06:23 AM12/17/21
to James M. Kennedy, LS-DYNA2
Dear Len & James,

Thanks for all the links, I'm sure they are going to be very useful in solving my problem.

James M. Kennedy <j...@kbs2.com>, 17 Ara 2021 Cum, 02:52 tarihinde şunu yazdı:

James M. Kennedy

unread,
Dec 17, 2021, 2:31:04 PM12/17/21
to Furkan Keleş, LS-DYNA2

Dear Furkan,

 

LS-DYNA wave-structure interaction simulations using ALE (standard) and S-ALE  models.

 

------------------

 

ALE (standard) - Intermediate example presentation

 

https://www.dynaexamples.com/ale/intermediate-example

 

Intermediate : Wave-Structure Interaction

 

A container is partially filled with water subjected to gravity. As the wave propagates, it hits a rigid column placed roughly in the middle of the container. This model was developed to make/allow comparisons with an available S-ALE model. To make a direct comparison, using identical cell discretization (0.010 meters) one can modify the S-ALE model (wavehitcol.k) below by making a simple parameter change (node numbering) in the three *ALE_STRUCTURED_MESH_CONTROL_POINTS keyword entries.

 

Roughly 330,000 elements are in the ALE and S-ALE models with a cell discretization of 0.010 meters (m), while roughly 2,590,000 elements are in the S-ALE model (wavehitcol.k) with a  cell discretization of 0.005 meters (m).

 

The Murnaghan equation of state has become available as a new keyword: *EOS_MURNAGHAN. This equation of state can be used to enforce quasi-static incompressibility while allowing for a reasonable time step in explicit calculations.  Any ALE model part is therefore represented by a *MAT_NULL card, describing the density and viscosity of the material, and an *EOS_MURNAGHAN card, describing the pressure-density relationship.

 

------------------

 

S-ALE - Intermediate example presentation

 

http://ftp.lstc.com/anonymous/outgoing/hao/sale/models/wavehitcol/

 

Intermediate : Wave-Structure Interaction

 

A container is partially filled with water subjected to gravity. As the wave propagates, it hits a rigid column placed roughly in the middle of the container. The total contact forces are recorded for the entire simulation, and compared to experimental data [1].

 

[1} Gomez-Gesteira, M., Crespo, A.J.C., Rogers, B.D. Dalrymple, R.A., Dominguez, J.M. and Barreiro, A., “SPHysics - Development of a Free-Surface Fluid Solver - Part 2: Efficiency and Test Cases,” Computers & Geosciences, Vol. 48, pp. 300–307, November, 2012.

 

https://www.sciencedirect.com/science/article/pii/S0098300412000842

 

[2] Chen, H., "LS-DYNA Structured ALE (S-ALE) Solver", 14th International LS-DYNA Users Conference, Dearborn, Michigan, June, 2016.

 

http://www.dynalook.com/14th-international-ls-dyna-conference/simulation/ls-dyna-r-structured-ale-s-ale-solver

 

[3] Chen, H, and Do, I., "Recent Developments in LS-DYNA S-ALE", 15th International LS-DYNA Users Conference, Dearborn, Michigan, June, 2018.

 

https://www.dynalook.com/15th-international-ls-dyna-conference/fsi-ale/recent-developments-in-ls-dyna-r-s-ale

 

[4] Chen, H., "ALE_STRUCTURED_FSI The New S-ALE Solver", 16th International LS-DYNA Users Conference, Dearborn, Michigan, June, 2020.

 

https://www.dynalook.com/conferences/16th-international-ls-dyna-conference/fsi-ale-computing-technology-t8-1/t8-1-b-fsi-ale-081.pdf

 

-------------------

 

Sincerely,

James M. Kennedy

KBS2 Inc.

December 17, 2021

 

 

From: ls-d...@googlegroups.com [mailto:ls-d...@googlegroups.com] On Behalf Of Furkan Keles
Sent: Thursday, December 16, 2021 11:32 AM
To: LS-DYNA2 <ls-d...@googlegroups.com>
Subject: [LS-DYNA2] 2D - FSI Leakage with Erosion Problem

 

Dear all,

--

Reply all
Reply to author
Forward
0 new messages