Arbitrary Unstructured Grids (arbitrary + unstructured_grid)

Distribution by Scientific Domains


Selected Abstracts


An accurate gradient and Hessian reconstruction method for cell-centered finite volume discretizations on general unstructured grids

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 9 2010
Lee J. Betchen
Abstract In this paper, a novel reconstruction of the gradient and Hessian tensors on an arbitrary unstructured grid, developed for implementation in a cell-centered finite volume framework, is presented. The reconstruction, based on the application of Gauss' theorem, provides a fully second-order accurate estimate of the gradient, along with a first-order estimate of the Hessian tensor. The reconstruction is implemented through the construction of coefficient matrices for the gradient components and independent components of the Hessian tensor, resulting in a linear system for the gradient and Hessian fields, which may be solved to an arbitrary precision by employing one of the many methods available for the efficient inversion of large sparse matrices. Numerical experiments are conducted to demonstrate the accuracy, robustness, and computational efficiency of the reconstruction by comparison with other common methods. Copyright © 2009 John Wiley & Sons, Ltd. [source]


Voronoi cell finite difference method for the diffusion operator on arbitrary unstructured grids

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING, Issue 1 2003
N. SukumarArticle first published online: 11 MAR 200
Abstract Voronoi cells and the notion of natural neighbours are used to develop a finite difference method for the diffusion operator on arbitrary unstructured grids. Natural neighbours are based on the Voronoi diagram, which partitions space into closest-point regions. The Sibson and the Laplace (non-Sibsonian) interpolants which are based on natural neighbours have shown promise within a Galerkin framework for the solution of partial differential equations. In this paper, we focus on the Laplace interpolant with a two-fold objective: first, to unify the previous developments related to the Laplace interpolant and to indicate its ties to some well-known numerical methods; and secondly to propose a Voronoi cell finite difference scheme for the diffusion operator on arbitrary unstructured grids. A conservation law in integral form is discretized on Voronoi cells to derive a finite difference scheme for the diffusion operator on irregular grids. The proposed scheme can also be viewed as a point collocation technique. A detailed study on consistency is conducted, and the satisfaction of the discrete maximum principle (stability) is established. Owing to symmetry of the Laplace weight, a symmetric positive-definite stiffness matrix is realized which permits the use of efficient linear solvers. On a regular (rectangular or hexagonal) grid, the difference scheme reduces to the classical finite difference method. Numerical examples for the Poisson equation with Dirichlet boundary conditions are presented to demonstrate the accuracy and convergence of the finite difference scheme. Copyright © 2003 John Wiley & Sons, Ltd. [source]


Numerical simulation of a single bubble by compressible two-phase fluids

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 6 2010
Siegfried Müller
Abstract The present work deals with the numerical investigation of a collapsing bubble in a liquid,gas fluid, which is modeled as a single compressible medium. The medium is characterized by the stiffened gas law using different material parameters for the two phases. For the discretization of the stiffened gas model, the approach of Saurel and Abgrall is employed where the flow equations, here the Euler equations, for the conserved quantities are approximated by a finite volume scheme, and an upwind discretization is used for the non-conservative transport equations of the pressure law coefficients. The original first-order discretization is extended to higher order applying second-order ENO reconstruction to the primitive variables. The derivation of the non-conservative upwind discretization for the phase indicator, here the gas fraction, is presented for arbitrary unstructured grids. The efficiency of the numerical scheme is significantly improved by employing local grid adaptation. For this purpose, multiscale-based grid adaptation is used in combination with a multilevel time stepping strategy to avoid small time steps for coarse cells. The resulting numerical scheme is then applied to the numerical investigation of the 2-D axisymmetric collapse of a gas bubble in a free flow field and near to a rigid wall. The numerical investigation predicts physical features such as bubble collapse, bubble splitting and the formation of a liquid jet that can be observed in experiments with laser-induced cavitation bubbles. Opposite to the experiments, the computations reveal insight to the state inside the bubble clearly indicating that these features are caused by the acceleration of the gas due to shock wave focusing and reflection as well as wave interaction processes. While incompressible models have been used to provide useful predictions on the change of the bubble shape of a collapsing bubble near a solid boundary, we wish to study the effects of shock wave emissions into the ambient liquid on the bubble collapse, a phenomenon that may not be captured using an incompressible fluid model. Copyright © 2009 John Wiley & Sons, Ltd. [source]


Employment of the second-moment turbulence closure on arbitrary unstructured grids

INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS, Issue 4 2004
B. BasaraArticle first published online: 12 JAN 200
Abstract The paper presents a finite-volume calculation procedure using a second-moment turbulence closure. The proposed method is based on a collocated variable arrangement and especially adopted for unstructured grids consisting of ,polyhedral' calculation volumes. An inclusion of 23k in the pressure is analysed and the impact of such an approach on the employment of the constant static pressure boundary is addressed. It is shown that this approach allows a removal of a standard but cumbersome velocity,pressure ,Reynolds stress coupling procedure known as an extension of Rhie-Chow method (AIAA J. 1983; 21: 1525,1532) for the Reynolds stresses. A novel wall treatment for the Reynolds-stress equations and ,polyhedral' calculation volumes is presented. Important issues related to treatments of diffusion terms in momentum and Reynolds-stress equations are also discussed and a new approach is proposed. Special interpolation practices implemented in a deferred-correction fashion and related to all equations, are explained in detail. Computational results are compared with available experimental data for four very different applications: the flow in a two-dimensional 180o turned U-bend, the vortex shedding flow around a square cylinder, the flow around Ahmed Body and in-cylinder engine flow. Additionally, the performance of the methodology is assessed by applying it to different computational grids. For all test cases, predictions with the second-moment closure are compared to those of the k,,model. The second-moment turbulence closure always achieves closer agreement with the measurements. A moderate increase in computing time is required for the calculations with the second-moment closure. Copyright © 2004 John Wiley & Sons, Ltd. [source]