» Home  » Research  » CASA projects 
logo logo
  Home About us Research Education Vacancies Meetings Newsletter
Topics
« Research  
» PhD projects  
» A Platform for Numerical Computations ...  
  W.D. Drenth (2003)
» A methematical model analysing...  
  F. Hagman (2005)
» A variational model for diblock...  
  Y. van Gennip (2008)
» Acoustic liner - mean flow interaction  
  M. Darau (2012)
» Analysis of Eddy Currents...  
  J.M.B. Kroot (2005)
» Analysis of the Flow Instabilities in the Extrusion of Polymeric Melts  
  A.C.T. Aarts (1997)
» Anisotropic turbulence transport  
  R. Minero (2006)
» Applications of model order reduction for IC modeling  
  M.V. Ugryumova (2011)
» Aspects of Solving Non-Linear Boundary Value Problems Numerically  
  M.E. Kramer (1992)
» Automated focusing and astigmatism correction in electron microscopy  
  M.E. Rudnaya (2011)
» BEM simulation for glass parisons  
  K. Wang (2002)
» Capturing Detonation Waves for the Reactive Euler Equations  
  A.C. Berkenbosch (1995)
» Cartilaginous Tissues  
  A.J.H. Frijns (2000)
» Combustion associated noise...  
  M.L. Bondar (2007)
» Condition number of ...  
  W. Dijkstra (2008)
» Constitutive modeling of concentrated solutions of main-chain liquid crystalline polymers  
  O. Matveichuk (2013)
» Constraints in applied mathematics  
  R. Planqué (2005)
» Contour Dynamics  
  P.W.C. Vosbeek (1998)
» Contour Dynamics Simulations  
  R.M. Schoemaker (2003)
» Cooling machinery  
  I.A. Lyulina (2004)
» Coupled-Wave Analysis  
  M. van Kraaij (2011)
» Curved Jets of Viscous Fluid: Interactions with a Moving Wall  
  A. Hlod (2009)
» DEM simulations of toner  
  I.E.M. Severens (2005)
» Development of Maxwell's Solver  
  Shcherbakov E.
» Dynamic capillarity in porous media  
  Y. Fan (2012)
» Efficient simulation of flow and heat transfer in arbitrarily shaped pipes  
  P.I. Rosen Esquivel (2012)
» Electric Circuit Simulation  
  S.H.M.J. Houben (2003)
» Electrochemical Drilling  
  M.J. Noot (1997)
» Energy-conserving discretization methods for the incompressible Navier-Stokes equations  
  B. Sanderse (2013)
» Finite Antenna Arrays  
  Dave Bekers (2004)
» Flow Front Instabilities in ...  
  H.J.J Gramberg (2005)
» Fracture Mechanics  
  M.J. Patricio Dias (2008)
» Hele-Shaw and Stokes flow with a source or sink: Stability of spherical solutions  
  E. Vondenhoff (2009)
» High performance circuit  
  A. Verhoeven (2008)
» High-amplitude oscillatory gas ...  
  P. in 't panhuis (2009)
» Laminar flames  
  M.G. Graziadei (2004)
» Magnetohydrodynamic Waves and Instabilities in Rotating Tokamak Plasmas  
  J.W. Haverkort (2013)
» Mathematical Models ...  
  D. Bezanovic (2005)
» Microscopic Interpretation of Wasserstein Gradient Flows  
  D.R.M. Renger (2013)
» Mixed Finite Elements ...  
  K. Malakpoor (2007)
» Model Order Reduction for Coupled Systems using Low-rank Approximationsa  
  A. Lutowska (2012)
» Model Order Reduction for Multi-terminal Systems with Applications to Circuit Simulation  
  R. Ionutiu (2011)
» Model order reduction and sensitivity analysis  
  Z. Ilievski (2011)
» Modelling laser percussion drilling  
  J.C.J. Verhoeven (2008)
» Modelling of the Glass Press-Blow Process  
  S.M.A. Allaart-Bruin
» Multi-Scale Riemann-Finsler Geometry. Applications to ...  
  L. Astola (2010)
» Multi-Valued Geodesic Tractography for Diffusion Weighted Imaging  
  N. Sepasian (2011)
» Multibody Dynamics  
  P.M.E.J. Wijckmans (1996)
» Multibody systems  
  B. Tasic (2004)
» Multiscale Reaction-Diffusion Systems Describing Concrete Corrosion: Modeling and Analysis  
  T. Fatima (2013)
» Numerical Analysis of Viscous Flow  
  V. Nefedov (2001)
» Numerical Aspects of Laminar Flame Simulation  
  B. van 't Hof (1998)
» Numerical methods in combustion  
  M.J.H. Anthonissen (2001)
» Numerical optimization of the air film cooling  
  M. Sizov (2007)
» Numerical shape optimisation in blow moulding  
  J.A.W.M. Groot (2011)
» Numerical simulation of a three-stage Stirling-type pulse-tube refrigerator  
  M.A. Etaati (2011)
» Parabolic evolution equations for quasistationary free boundary problems in capillary fluid mechanics  
  G. Prokert (1997)
» Perturbation and Operator Methods for Solving Stokes Flow and Heat Flow Problems  
  T.C. Chandra (2002)
» Positive operator-valued measures and phase-space representations  
  R. Beukema (2003)
» Pressing of Glass in Bottle and Jar Manufacturing: Numerical Analysis and Computation  
  K.Y. Laevsky (2003)
» Radiative Heat Transfer in Glass: The Algebraic Ray Trace Method  
  B.J. van der Linden (2002)
» Representation and Manipulation of Images Based on Linear Functionals  
  B.J. Janssen (2009)
» Scattering from Finite Structures: An Extended Fourier Modal Method  
  M. Pisarenco (2011)
» Simulating Unsteady Conduit Flows with Smoothed Particle Hydrodynamics  
  Q. Hou (2012)
» Simulation and modelling underlying radio frequencies  
  P.J. Heres (2005)
» Solving Boundary Value Problems on Composite Grids ...  
  P.J.J. Ferket (1996)
» Stability and Evolution  
  G.J.M.Pieters (2004)
» Stability of Immersed Liquid Threads  
  A.Y. Gunawan (2004)
» Stream Function Approach for Determining Optimal Surface Currents  
  G.N. Peeren (2003)
» The Boundary Element Method: Errors and gridding for problems with hot spots  
  G. Kakuba (2011)
» The Rigorous Coupled-Wave Analysis  
  N.P. van der Aa (2007)
» Upscaling of Reactive Flows  
  K. Kumar (2012)
» Viscous Sintering  
  G.A.L. van de Vorst (1994)
» Visualisation and Simulation with Object-Oriented Networks  
  A.C. Telea (2000)

External links
» On-line dissertations

A Parallel Hierarchical-Element Method for Contour Dynamics Simulations

R.M. Schoemaker

Introduction

The HEM is a modified contour dynamics (CD) method in the sense that it accelerates highly complex 2D inviscid incompressible vortex flows for 1 processor systems [1]. The HEM has been parallelised on its turn for even faster computation, using OpenMP for up to 16 processors on an Origin 3800 shared-memory archictecture. CD is based on the fact that the evolution of a patch of uniform vorticity is fully determined by the evolution of its bounding contour [2].

A continuous vorticity distribution (red) is approximated by nested patches of uniform vorticity (blue). The contours are discretised by nodes and linear elements. The increasing complexity of the evolving contours demand node redistribution.

The HEM makes use of a hierarchical tree of levels with boxes and is based on the Fast Multipole Method, i.e. sources (vorticity in this case) are approximated by a computationally simpler description. The size of of the collected sources depend on the distance to the evaluation point.

Example of the velocity contributions from the far field to a certain node in the blue box. Faraway vorticity is approximated by large rings. These contributions are passed on to a dotted ring which lies above the blue box in the hierarchical tree and so on for each level. The light orange areas indicate regions for the optimal ring size at each level. The little dark orange area just around the blue box contain conventional interactions due to accuracy requirements.

Parallelisation Strategy

The HEM has been parallelised along the boxes one level at a time in the hierarchical tree. A fixed domain decomposition has been compared to a dynamical domain scheduling policy. Domain decomposition implies assigning processors to certain areas of a boxed level. In the static scheduling policy of OpenMP, processors have to wait for each other. The  dynamic scheduling permits `grazing' of the level, implying a load balancing effect for non-uniform node distributions.


Numerical Experiments and Discussion

In the load-imbalanced case (left) the domain is subdivided into P=2 fixed areas. In the load-balanced case (right) two processors 'graze' the domain with one box per processor and if a processor has finished one box it will resume on the next available box and will not wait for the other processors ($P = 2 is arbitrary). In the load-imbalanced case (left) the domain is subdivided into P=2 fixed areas. In the load-balanced case (right) two processors 'graze' the domain with one box per processor and if a processor has finished one box it will resume on the next available box and will not wait for the other processors (P = 2 is arbitrary)

References

1. P.W.C. Vosbeek, H.J.H. Clercx, and R.M.M. Mattheij, Acceleration of ContourDynamics Simulations with a Hierarchical-Element Method, J. Comp. Phys., 161, 287--311 (2000)
2. D.G. Dritschel, Contour Dynamics and Contour Surgery: Numerical algorithms for extendedhigh-resolution modelling of vortex dynamics in two-dimensional, inviscid, incompressible flows, Comp. Phys. Rep., 10, 77 (1989)

© Centre for Analysis, Scientific Computing and Applications. For questions please refer to the editor.
This page modified: Mon Nov 28 13:55:14 CET 2005