Adaptive time stepping for explicit euler implementation of spherical and non-spherical particle speed up
Numerical implementation schemes of drag force effects on Lagrangianparticles can lead to instabilities or inefficiencies if static particle timestepping is used. Despite well known disadvantages, the programmingstructure of the underlying, C++ based, Lagrangian particle solver led to thechoice of a...
| Main Authors: | , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
MULTIPHYSICS
2016-04-01
|
| Series: | International Journal of Multiphysics |
| Online Access: | http://journal.multiphysics.org/index.php/IJM/article/view/103 |
| Summary: | Numerical implementation schemes of drag force effects on Lagrangianparticles can lead to instabilities or inefficiencies if static particle timestepping is used. Despite well known disadvantages, the programmingstructure of the underlying, C++ based, Lagrangian particle solver led to thechoice of an explicit EULER, temporal discretization scheme. To optimizethe functionality of the EULER scheme, this paper proposes a method ofadaptive time stepping, which adjusts the particle sub time step to the needof the individual particle. A user definable adjustment between numericalstability and calculation efficiency is sought and a simple time stepping ruleis presented. Furthermore a method to quantify numerical instability isdevised and the importance of the characteristic particle relaxation time asnumerical parameter is underlined. All derivations are being conducted for(non-)spherical particles and finally for a generalized drag forceimplementation. Important differences in spherical and non-sphericalparticle behaviour are pointed out. |
|---|---|
| ISSN: | 1750-9548 2048-3961 |