Speaker
Prof.
Vladimir Gerdt
(Joint Institute for Nuclear Research)
Description
We apply modern computer algebra methods and software for the Lie symmetry analysis to the governing Cosserat system of twelve nonlinear partial differential equations (PDEs) describing dynamics of nearly one-dimensional flexible structures (rods, fibers, cables, etc.) and construct the general analytical solution to the kinematic part (six equations) of the system. The obtained solution depends on two arbitrary analytical vector functions, and we show that its knowledge helps to resolve stiffness of the governing Cosserat system. Based on this result we develope algorithms based on combinations of numerical (exponential integration) and analytical (solver of differential equations built-in Maple) treatments of the dynamical part of the governing PDE system. Our approach allows for larger step sizes compared to pure numerical solvers and at the same time combines efficiency and accuracy without sanctifying one for another. Based on this observation, we create a hybrid semi-analytical solver for highly viscous two-way coupled fluid-rod problems which allows for the interactive high-fidelity simulations of flagellated microswimmers, as a result of a substantial reduction of the numerical stiffness. Besides, we present in the talk experimental comparison of our method with the so-called $\alpha-$method, the best one among pure numerical methods developed for integrating Cosserat equations, and demonstrate superiority of our method.
Primary author
Prof.
Vladimir Gerdt
(Joint Institute for Nuclear Research)
Co-authors
Prof.
Andreas Weber
(Institute of Computer Science II, University of Bonn, Bonn, Germany)
Dr
Dmitry Lyakhov
(Visual Computing Center, King Abdullah University of Science and Technology, Thuwal, Saudi Arabia)
Prof.
Dominik Michels
(Visual Computing Center, King Abdullah University of Science and Technology, Thuwal, Saudi Arabia)