Adaptive Grid Refinement for Free-Surface Hydrodynamic Flows
Résumé
Free-surface water flows feature the creation of gravity waves due to the interaction of the surface with immersed bodies, as well as the propagation of these waves and their interaction with other flow phenomena. For Reynolds-Averaged Navier Stokes (RANS) simulation of such flows, grids are needed that are refined in the vicinity of the water surface in order to resolve the surface model (for example, a level set or volume-of-fluid equation). Furthermore, fine grids are needed below the surface for the accurate resolution of the pressure-velocity fields which drive the surface deformation. The goal of the research presented in this paper is to create such refined meshes with automatic adaptive mesh refinement. The work is based on the ISIS-CFD unstructured finite-volume RANS solver which our group develops. This solver contains a surface-capturing model for the water surface, where the surface is represented as a discontinuity in the volume fraction for the water ; the volume fraction is computed with a linear convection equation. The automatic grid refinement method developed for ISIS-CFD [1] performs anisotropic refinement of unstructured hexahedral grids by division of the cells. Like the solver, it is fully parallel. The first part of the paper concerns the appropriate choice of a refinement criterion, which indicates where the grid should be refined. We propose combined refinement criteria , consisting of refinement around the surface in order to resolve correctly the volume fraction equation, combined with pressure Hessian-based refinement below the surface to capture the pressure and velocity fields associated with the waves. The computation of the Hessian matrix is complicated by the variations in cell size that appear in unstruc-tured hexahedral meshes, which lead to large errors in the computed second derivatives. A double application of Gauss's theorem, followed by smoothing of the obtained first and second derivatives to eliminate the irregularities, is proposed. The difficulty of eliminating numerical errors while preserving the effects of physical pressure peaks is discussed.
Origine : Fichiers produits par l'(les) auteur(s)
Loading...