Computational Modeling of Individual Red Blood Cell Dynamics Using Discrete Flow Composition and Adaptive Time-Stepping Strategies

In this article, we present a finite element method for studying the dynamic behavior of deformable vesicles, which mimic red blood cells, in a non-Newtonian Casson fluid. The fluid membrane, represented by an implicit level-set function, adheres to the Canham-Helfrich model and maintains surface inextensibility constraint through penalty. We propose a two-step time integration scheme that incorporates higher-order accuracy by using an asymmetric composition of discrete flow based on the second-order backward difference formula, followed by a projection onto the real axis. Our framework incorporates variable time steps generated by an appropriate adaptation criterion. We validate our model through numerical simulations against existing experimental and numerical results in the case of purely Newtonian flow. Furthermore, we provide preliminary results demonstrating the influence of the non-Newtonian fluid model on membrane regimes.

Automated summary

In this article, a finite element method is presented to study the dynamic behavior of deformable vesicles in a non-Newtonian Casson fluid. The fluid membrane, represented by an implicit level-set function, follows the Canham-Helfrich model and maintains surface inextensibility through penalty. A two-step time integration scheme is proposed, incorporating higher-order accuracy and variable time steps generated by an appropriate adaptation criterion. The model is validated through numerical simulations and its influence on membrane regimes in a non-Newtonian fluid is demonstrated.