One or more methods of designing microchannel fluid flow networks in a fuel cell bipolar plate includes executing one or more programs on one or more computing devices having one or more processors to optimize the spatially varying orientations of homogenized anisotropic porous media by iteratively executing a gradient-based algorithm that incorporates objective functions of reaction uniformity and flow resistance, and then generate, in response to the homogenized anisotropic porous media optimization, one or more microchannel fluid flow networks by dehomogenizing the optimized anisotropic porous media.