Application of grid convergence index to shock wave validated with LS-DYNA and ProsAir

Abstract

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. RESUMENLa 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.Palabras clave: ındice de convergencia de malla (GCI), equivalente de TNT, LS-DYNA, ProsAir.How to ...

Accuracy of CFD Simulations on Indoor Air Ventilation: Application of Grid Convergence Index on Underfloor Air Distribution (UFAD) System Design

“…Users can see how big of an error that might be making and where it might fall inside the error band. While low GCI values indicate that the computational solution is close to the asymptotic on the mesh size, this also demonstrates how much the solution will vary when the mesh size is refined [30], [31]- [32]. If the solution is already good enough, making the mesh smaller will not make it much different.…”

“…Figure 5 shows a scatter plot of measured versus modelled values for horizontal, along with three velocity components, Ux. Again, all data points (32) in the vertical profiles are included. Figure 6 compares the graph of Chen's [10] simulation and the present simulation of RNG k-ɛ and STD k-ɛ model.…”

Numerical Study of the Opening Angle Incidence in Michell-Banki Turbine’s Performance without Guide Blades

“…In numerical studies, the estimation of the discretization error is a transcendental part since it determines the fit and agreement of the numerical results with the experimental data [40]. Computer software packages can solve a set of differential or integral equations that describe the physical phenomenon that governs the analyzed situation.…”

Prediction of blast loading on protruded structures using machine learning methods

“…The three-dimensional blast simulation program ProSAir, a compressible CFD package, was used to collect maximum overpressure and impulse values at various locations of the protruded structural façade. This package has been validated against experimental data and is widely used for blast load modeling (Alterman et al, 2019; Castedo et al, 2019; Remennikov and Rose, 2005). Detailed information on the underlying computational equations and assumptions for ProSAir can be found on the Cranfield University website (ProSAir computational blast loading tool, 2022).…”