1 / 27

IV European Conference of Computational Mechanics

IV European Conference of Computational Mechanics Hrvoje Gotovac , Veljko Srzić, Tonći Radelja, Vedrana Kozulić University of Split, Department of Civil and Architectural Engineering, Croatia Explicit Adaptive Fup Collocation Method (EAFCM) for solving the parabolic problems

missy
Download Presentation

IV European Conference of Computational Mechanics

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. IV European Conference of Computational Mechanics Hrvoje Gotovac, Veljko Srzić, Tonći Radelja, Vedrana Kozulić University of Split, Department of Civil and Architectural Engineering, Croatia Explicit Adaptive Fup Collocation Method (EAFCM) for solving the parabolic problems Presentation ECCM, 21 May 2010, Paris, France.

  2. Presentation outline • General concept • Fup basis functions • Fup collocation transform (FCT) - space approximation • Explicit time integration for parabolic stiff problems • Numerical examples • Conclusions • Future directions

  3. 1. General concept • Developing adaptive numerical method which can deal with parabolic flow and transport stiff problems having wide range of space and temporal scales • Ability to handle multiple heterogeneity scales • Application target: unsaturated and multiphase flow, reactive transport and density driven flow in porous media, as well as structural mechanics problems

  4. Saturated – unsaturated flow

  5. Interaction between surface and subsurface flow

  6. Geothermal convective processes in porous media

  7. Typical physical and numerical problems • Description of wide range of space and temporal scales • Sharp gradients, fronts and narrow transition zones (‘fingering‘ and ‘layering’) • Artificial oscillations and numerical dispersion – advection dominated problems • Description of heterogeneity structure • Strong nonlinear and coupled system of equations

  8. Motivation for EAFCM • Multi-resolution and meshless approach • Continuous representation of variables and all its derivatives (fluxes) • Adaptive strategy • Method of lines (MOL) • Explicit formulation (no system of equations!!!) • Perfectly suited for parallel processing

  9. Flow chart of EAFCM

  10. 2. Fup basis functions • Atomic or Rbf class of functions • Function up(x) • Fourier transform of up(x) function • Function Fupn(x)

  11. Function Fup2(x)

  12. Function Fup2(x,y)

  13. Function Fupn(x) • Compact support • Linear combination of n+2 functions Fupn(x) exactly presents polynomial of order n • Good approximation properties • Universal vector space UP(x) • Vertexes of basis functions are suitable for collocation points • Fupn(x,y) is Cartesian product of Fupn(x) and Fupn(y)

  14. 3. Fup collocation transform (FCT) • Any function u(x) is presented by linear combination of Fup basis functions: j - level (from zero to maximum level J) jmin - resolution at the zero level k - location index in the current level - Fup coefficients - Fup basis functions n - order of the Fup basis function

  15. Spatial derivatives

  16. 4. Time numerical integration • Reduces to system of Ordinary Differential Equations (ODE) for adaptive grid and every time step (t – t+dt): • With appropriate initial conditions:

  17. Stabilized second-order Explicit Runge-Kutta method (SERK2) • Recently developed by Vaquero and Janssen (2009) • Extended stability domains along the negative real axis • Suitable for very large stiff parabolic ODE • Second – order method up to 320 stages • Public domain Fortran routine SERK2

  18. 5. Numerical examples • 1-D density driven flow problem • 2-D Henry salwater intrusion problem

  19. Mathematical model • Pressure-concentration formulation • Fluid mass balance: • Salt mass balance:

  20. Number of collocation points and compression coefficient

  21. 6. Conclusions • Development of mesh-free adaptive collocation algorithm that enables efficient modeling of all space and time scales • Main feature of the method is the space adaptation strategy and explicit time integration • No discretization and solving of huge system of equations • Continuous approximation of fluxes

  22. 7. Future directions • Multiresolution description of heterogeneity • Development of 3-D parallel EAFCM • Time subdomain integration • Description of complex domain with using other families of atomic basis functions • Further application to mentioned processes in porous media and other (multiphysics) problems

More Related