Published
Application of grid convergence index to shock wave validated with LS-DYNA and ProsAir
Aplicación del índice de convergencia de malla frente a onda de choque validado con LS-DYNA y ProsAir
Keywords:
Grid convergence index (GCI), TNT equivalent, LS-DYNA, ProsAir. (en)índice de convergencia de malla (GCI), equivalente de TNT, LS-DYNA, ProsAir (es)
The discretization error is not always calculated, even though it is essential for the studies of computational solid mechanics. However, it is well known that an error committed by the mesh used can be as large as the measured variable, which greatly invalidates the results obtained. The grid convergence index (GCI) method makes possible to determine on a solid basis, the order of convergence and the asymptotic solution. This method seems to be a suitable estimator despite further research is needed in the context of blast situations and finite element (FE) calculations. For this purpose, field trials were performed consisting in the detonation of a spherical hanging load of homemade explosive. The pressure generated by the shock wave was measured in different positions at two distances. With these data, a TNT equivalent has been obtained and used to calculate the shock propagation with the solvers LS-DYNA and ProsAir. This work aims to verify the GCI method by comparing its results with field data along with the simulations carried out. The comparison also seeks to validate the methodology used to obtain the TNT equivalent.
This research shows that the GCI gives good results for both solvers despite the complexity of the physical problem. Besides, LS-DYNA displays better correlation with the experimental data than the ProsAir results, with an error of less than 10% in all values.
La estimación del error de discretización no siempre se calcula, aunque es algo fundamental para el estudio de la mecánica computacional de los sólidos. Sin embargo, es bien sabido que el error cometido por la malla utilizada puede ser del mismo orden que la variable medida, lo que llega a invalidar los resultados obtenidos. El método del índice de convergencia de la malla (GCI) permite determinar sobre una base sólida el orden de convergencia y la solución asintótica, por lo que parece ser un buen estimador, a pesar de que es necesario seguir investigando en el contexto de las situaciones de ondas de choque (explosivos) y de los cálculos de elementos finitos (FE). Para este fin, se realizaron pruebas de campo consistentes en la detonación de una carga esférica colgada de explosivo casero. La presión generada por la onda de choque se midió en diferentes posiciones a dos distancias. Con estos datos, se obtuvo un equivalente de TNT que se utilizó para calcular la propagación del choque con los programas LS-DYNA y ProsAir. Este trabajo pretende verificar el método GCI comparando sus resultados con los datos de campo junto con las simulaciones realizadas. También, la comparación busca validar la metodología empleada para la obtención del equivalente TNT.
La investigación muestra que el GCI da buenos resultados para ambos programas a pesar de la complejidad del problema físico. Además, el LS-DYNA produce una mejor correlación con los datos experimentales que los aportados por el ProsAir, con todos los valores por debajo del 10 % de error.
References
Alañón, A., Cerro-Prada, E., Vázquez-Gallo, M. J., and Santos, A. P. (2018). Mesh size effect on finite-element modeling of blast-loaded reinforced concrete slab. Engineering with Computers, 34(4), 649-658. DOI: 10.1007/s00366-017-0564-4
Anderson, A. E., Ellis, B. J., and Weiss, J.A. (2007) Verification, validation and sensitivity studies in computational biomechanics. Computer Methods in Biomechanics and Biomedical Engineering, 10(3), 171-184. DOI: 10.1080/10255840601160484
Chiquito, M., Castedo, R., López, L. M., Santos, A. P., Mancilla, J. M., and Yenes, J. I. (2019). Blast Wave Characteristics and TNT Equivalent of Improvised Explosive Device at Small scaled Distances. Defence Science Journal, 69(4), 328-335. DOI: 10.14429/dsj.69.13637
Forth, S. (2012). ProSAir User Guide. Cranfield, UK: Crandfield University, Applied Mathematics and Scientific Computing Department.
Huang, Y., Willford, M. R., and Schwer, L. E. (2012, June). Validation of LS-DYNAr MMALE with blast experiments. Paper presented at the 12th International LS-DYNA Users Conference, Detroit, DYNAmore. https://www.dynalook.com/conferences/12th-international-ls-dyna-conference/blast-impact20-c.pdf/view
Jin, Y., and Shaw, B. D. (2010). Computational modeling of n-heptane droplet combustion in air–diluent environments under reduced-gravity. International Journal of Heat and Mass Transfer, 53(25-26), 5782-5791. DOI: 10.1016/j.ijheatmasstransfer.2010.08.005
Kwasniewski, L. (2013). Application of grid convergence index in ´ FE computation. Bulletin of the Polish Academy of Sciences: Technical Sciences, 61(1), 123-128. DOI: 10.2478/bpasts2013-0010
Lapoujade, V., Van Dorsselaer, N., Kevorkian, S., and Cheval, K. (2010, June). A study of mapping technique for air blast modelling. Paper presented at the 11th International LS-DYNA Users Conference. Detroit, DYNAmore. https: //www.dynalook.com/conferences/international-conf-2010/BlastImpact-1-3.pdf/view
Lee, E., Finger, M., and Collins, W. (1973). JWL equation of state coefficients for high explosives [Technical report No. UCID16189]. Livermore, CA: Lawrence Livermore National Lab (LLNL). DOI: 10.2172/4479737
Liang, Y., and Tao, L. (2017). Interaction of vortex shedding processes on flow over a deep-draft semisubmersible. Ocean Engineering, 141, 427-449. DOI: 10.1016/j.oceaneng.2017.06.056
LSTC. (2019). LS-DYNA Keyword User’s Manual (Volume II: Material Models). Livermore, CA: Livermore Software Technology Corporation. http://ftp.lstc.com/anonymous/outgoing/jday/manuals/LS-DYNA_Manual_Volume_II_R11_Ver2.pdf
Mobaraki, B., and Vaghefi, M. (2015). Numerical study of the depth and cross-sectional shape of tunnel under surface explosion. Tunnelling and Underground Space Technology, 47, 114-122. DOI: 10.1016/j.tust.2015.01.003
Ndebele, B. B., and Skews, B. W. (2019). The interaction of a cylindrical shock wave segment with a converging– diverging duct. Shock Waves, 29, 817-831. DOI: 10.1007/s00193-018-00888-7
Rebelo, H. B., and Cismasiu, C. (2017, May). A Comparison between three air blast simulation techniques in LS-DYNA. Paper presented at the 11th European LS-DYNA Conference, Salzburg, DYNAmore. https://www.dynalook.com/conferences/11th-european-ls-dyna-conference/air-blast-2/a-comparison-between-three-air-blast-simulation-techniques-in-ls-dyna/view
Rigby, S. E., Tyas, A., Bennett, T., Fay, S. D., Clarke, S. D., and Warren, J. A. (2014). A numerical investigation of blast loading and clearing on small targets. International Journal of Protective Structures, 5(3), 253-274. DOI: 10.1260/2041-4196.5.3.253
Rigby, S. E., Fay, S. D., Tyas, A., Warren, J. A., and Clarke, S. D. (2015). Angle of incidence effects on far-field positive and negative phase blast parameters. International Journal of Protective Structures, 6(1), 23-42. DOI: 10.1260%2F2041-4196.6.1.23
Roache, P. J. (1994). Perspective: a method for uniform reporting of grid refinement studies. Journal of Fluids Engineering, 116(3), 405-413. DOI: 10.1115/1.2910291
Rose, T. A. (2001). An approach to the evaluation of blast loads on finite and semi-infinite structures (Ph.D. thesis, Cranfield University). http://hdl.handle.net/1826/4262
Schwer, L. E. (2008). Is your mesh refined enough? Estimating discretization error using GCI. Paper presented at the German LS-DYNA Forum, Bamberg, DYNAmore. https://www.dynamore.de/de/download/papers/forum08/dokumente/I-I-03.pdf
Swisdak, M. M. (1975). Explosion effects and properties. Part I. Explosion effects in air [Technical report No. NSWC/WOL/TR-75-116]. White Oak, MD: Naval Surface Weapons Center, White Oak Lab Silver Spring. https://apps.dtic.mil/dtic/tr/fulltext/u2/a018544.pdf
Trajkovski, J., Kunc, R., Perenda, J., and Prebil, I. (2014). Minimum mesh design criteria for blast wave development and structural response-MMALE method. Latin American Journal of Solids and Structures, 11(11), 1999-2017. DOI: 10.1590/S1679-78252014001100006
How to Cite
License
Copyright (c) 2019 Ricardo Castedo, Carlos Reifarth, Anastasio P Santos, Jorge J Losada, Lina M López, Maria Chiquito, Juan M Mancilla
This work is licensed under a Creative Commons Attribution 4.0 International License.
The authors or holders of the copyright for each article hereby confer exclusive, limited and free authorization on the Universidad Nacional de Colombia's journal Ingeniería e Investigación concerning the aforementioned article which, once it has been evaluated and approved, will be submitted for publication, in line with the following items:
1. The version which has been corrected according to the evaluators' suggestions will be remitted and it will be made clear whether the aforementioned article is an unedited document regarding which the rights to be authorized are held and total responsibility will be assumed by the authors for the content of the work being submitted to Ingeniería e Investigación, the Universidad Nacional de Colombia and third-parties;
2. The authorization conferred on the journal will come into force from the date on which it is included in the respective volume and issue of Ingeniería e Investigación in the Open Journal Systems and on the journal's main page (https://revistas.unal.edu.co/index.php/ingeinv), as well as in different databases and indices in which the publication is indexed;
3. The authors authorize the Universidad Nacional de Colombia's journal Ingeniería e Investigación to publish the document in whatever required format (printed, digital, electronic or whatsoever known or yet to be discovered form) and authorize Ingeniería e Investigación to include the work in any indices and/or search engines deemed necessary for promoting its diffusion;
4. The authors accept that such authorization is given free of charge and they, therefore, waive any right to receive remuneration from the publication, distribution, public communication and any use whatsoever referred to in the terms of this authorization.
