Efficient Natural Gradient Descent Methods for Large-Scale PDE-Based Optimization Problems