INTEC   05402
INSTITUTO DE DESARROLLO TECNOLOGICO PARA LA INDUSTRIA QUIMICA
Unidad Ejecutora - UE
artículos
Título:
An h-Adaptive Solution of the Spherical Blast Wave Problem
Autor/es:
GUSTAVO A. RÍOS RODRIGUEZ; STORTI M.; LÓPEZ, E.; SARRAF, SOFÍA
Revista:
INTERNATIONAL JOURNAL OF COMPUTATIONAL FLUID DYNAMICS
Editorial:
TAYLOR & FRANCIS LTD
Referencias:
Lugar: Londres; Año: 2011 vol. 25 p. 31 - 39
ISSN:
1061-8562
Resumen:
Shock waves and contact discontinuities usually appear in compressible
flows, requiring a fine mesh in order to achieve an acceptable accuracy
of the numerical solution. The usage of a mesh adaptation strategy is
convenient as uniform refinement of the whole mesh becomes prohibitive
in three-dimensional (3D) problems. An unsteady h-adaptive strategy for
unstructured finite element meshes is introduced. Non-conformity of the
refined mesh and a bounded decrease in the geometrical quality of the
elements are some features of the refinement algorithm. A 3D extension
of the well-known refinement constraint for 2D meshes is used to enforce
a smooth size transition among neighbour elements with different levels
of refinement. A density-based gradient indicator is used to track
discontinuities. The solution procedure is partially parallelised, i.e.
the inviscid flow equations are solved in parallel with a finite element
SUPG formulation with shock capturing terms while the adaptation of the
mesh is sequentially performed. Results are presented for a spherical
blast wave driven by a point-like explosion with an initial pressure
jump of 105 atmospheres. The adapted solution is compared to
that computed on a fixed mesh. Also, the results provided by the theory
of self-similar solutions are considered for the analysis. In this
particular problem, adapting the mesh to the solution accounts for
approximately 4% of the total simulation time and the refinement
algorithm scales almost linearly with the size of the problem.