Improving mixing is one of the important goals in flow control, e.g., to decrease concentration polarization in membrane systems to reduce fouling. As with many distributed parameter systems, fluid flow can be controlled using boundary value manipulation. Fluid manipulation using electro-osmosis is studied in this work, where several cylindrical electrodes are used to create multiple spatially non-uniform time-varying electric fields. The proposed approach converts the distributed parameter system into an infinite-dimensional system by spatial and spectral discretization. A virtual output variable is constructed to allow the optimization of a mixing objective function to be conducted using frequency response analysis, with consideration of the constraints of conservation of charge. The solution obtained in this paper is the input profile that provides the greatest achievable ratio of time-average dissipation function to time-average input energy satisfying the input constraints. (C) 2016 Elsevier Ltd. All rights reserved.