EDP Sciences Journals List
Advanced Search
Subscriber Authentication Point
Home All issues Volume 1 Article
Free access article

Issue ESAIM: Proc.
Volume 1, 1996
Vortex flows an related numerical methods II
Page(s) 225 - 240
DOI http://dx.doi.org/10.1051/proc:1996039

ESAIM: Proc., 1996, Vol. 1, pp. 225-240
DOI: 10.1051/proc:1996039

Application of fast parallel and sequential tree codes to computing three-dimensional flows with the vortex element and boundary element methods

G.S. Winckelmans1, J.K. Salmon2, M.S. Warren3, A. Leonard4 and B. Jodoin5

1  Mechanical Engineering Dept. Université Catholique de Louvain Louvain-la-Neuve, B-1348, Belgium
2  Center for Advanced Computing Research California Institute of Technology Pasadena CA 91125, USA
3  Los Alamos National Laboratories Los Alamos, NM 87545, USA
4  Graduate Aeronautical Laboratories California Institute of Technology Pasadena CA 91125, USA
5  Chemical Eng. Dept. Université de Sherbrooke Sherbrooke (Qc), Canada


Abstract
A fast parallel oct-tree code originally developed for three-dimensional N-body gravitational simulations was modified into (1) a fast N-vortex code for viscous and inviscid vortex flow computations using the regularized vortex particle method (VEM), and (2) a fast N-panel code for solving boundary integral equations in potential flow aerodynamics using the boundary element method (BEM). The core of the fast tree code remains essentially unchanged between the different application codes: gravitation, VEM, BEM, etc. Only the modules that actually encode the physical model are changed. Particular attention is given to controlling the error introduced by the use of multipole expansions to represent the field produced by groups of elements, i.e., the tree code error. In particular, the acceptable error bound for use of any multipole expansion approximation is a run-time parameter. Program outputs include statistics on the errors for the field evaluation at all element locations. Problems in VEM and BEM involving N in the range 104 to over 106 are computed on parallel supercomputers. Problems with N in the range 103 to 105 are computed on workstations. Performance results are presented, together with sample computational results. For the VEM method, a high order particle redistribution scheme has been incorporated, in an efficient way, into the parallel tree code. It is applied, if necessary, to ensure that the core overlapping condition remains satisfied in long time computations. In addition, two different relaxation schemes have also been incorporated and partially tested. Such schemes are applied, if necessary, to ensure that the particle representation of the vorticity field remains a good representation of the true divergence free vorticity field in long time computations.



© EDP Sciences, ESAIM 1996


What is OpenURL?

The OpenURL standard is a protocol for transmission of metadata describing the resource that you wish to access. An OpenURL link contains article metadata and directs it to the OpenURL server of your choice. The OpenURL server can provide access to the resource and also offer complementary services (specific search engine, export of references...). The OpenURL link can be generated by different means.
  • If your librarian has set up your subscription with an OpenURL resolver, OpenURL links appear automatically on the abstract pages.
  • You can define your own OpenURL resolver with your EDPS Account. In this case your choice will be given priority over that of your library.
  • You can use an add-on for your browser (Firefox or I.E.) to display OpenURL links on a page (see http://www.openly.com/openurlref/). You should disable this module if you wish to use the OpenURL server that you or your library have defined.