Unlike macroscopic swimmers, microswimmers operate in a low-Reynolds-number regime dominated by viscous forces. This paper investigates the controllability of a magnetic microswimmer composed of a spherical magnetic head and an elastic, non-magnetic flagellum. The swimmer evolves in a Stokes flow and is modeled using the resistive force theory. We prove that, under planar motion, the system is not small-time locally controllable and numerically identify regions that remain inaccessible. Nevertheless, simulations show that trajectory tracking can still be achieved via Bayesian optimization, though it requires large-amplitude transverse deformations.

Trajectory tracking for microswimmers remains a key challenge in microrobotics, where low-Reynolds-number dynamics make control design particularly complex. In this work, we formulate the trajectory tracking problem as an optimal control problem and solve it using a combination of B-spline parametrization with Bayesian optimization, allowing the treatment of high computational costs without requiring complex gradient computations. Applied to a flagellated magnetic swimmer, the proposed method reproduces a variety of target trajectories, including biologically inspired paths observed in experimental studies. We further evaluate the approach on a three-sphere swimmer model, demonstrating that it can adapt to and partially compensate for wall-induced hydrodynamic effects. The proposed optimization strategy can be applied consistently across models of different fidelity, from low-dimensional ODE-based models to high-fidelity PDE-based simulations, showing its robustness and generality. These results highlight the potential of Bayesian optimization as a versatile tool for optimal control strategies in microscale locomotion under complex fluid-structure interactions.

Understanding and optimizing the design of helical microswimmers is crucial for advancing their application in various fields. This study presents an innovative approach combining Free-Form Deformation with Bayesian Optimization to enhance the shape of these swimmers. Our method facilitates the computation of generic swimmer shapes that achieve optimal average speed and efficiency. Applied to both monoflagellated and biflagellated swimmers, our optimization framework has led to the identification of new optimal shapes. These shapes are compared with biological counterparts, highlighting a diverse range of swimmers, including both pushers and pullers.