Article ID | Journal | Published Year | Pages | File Type |
---|---|---|---|---|
6933912 | Journal of Computational Physics | 2013 | 18 Pages |
Abstract
In multifluid problems with surface tension the fluid pressure and its derivative are discontinuous at fluid interfaces. We present a Cartesian grid embedded boundary method for numerically resolving these discontinuities in which we use Neumann-Neumann preconditioned iterative substructuring to solve the governing equations. We validate this method by computing several well-known Poisson problems with discontinuous coefficients, and we compare its performance to an approach based on simple iteration. By analogy with the conjugate gradient method, we hypothesize that the scaling of the Neumann-Neumann preconditioned iterative substructuring is O(h-Dlnh-1) where h is the cell size and D=2,3 is the dimensionality of the problem. In contrast, we show that the simple iterative procedure scales like O(h-(D+1)) and is slower by a factor of 4000 for a small (i.e., 64Ã64 cell) model calculation with physical parameters corresponding to a 1.5Â mm air bubble in water. We present an analytical model to explain the scaling of this iterative procedure.
Keywords
Related Topics
Physical Sciences and Engineering
Computer Science
Computer Science Applications
Authors
Gregory H. Miller, Elbridge Gerry Puckett,