Various studies in the biomechanics field have been conducted to analyze the aortic valvular disease and develop prosthetic heart valves. Optimization is an efficient solution to improve prosthetic heart valve performance. This thesis aims to optimize the bioprosthetic heart valve design to increase the coaptation area during the diastole and the orifice area during the systole. Large coaptation area prevents some possibilities of aortic regurgitation, and large orifice area increases the blood flow efficiency from the left ventricle to the aorta. The proposed solution for optimizing aortic heart valve design includes the surrogate management framework, which is a derivative-free pattern search optimization method, and isogeometric analysis, which is a method that tightly integrates geometric modeling and analysis. This combined approach significantly improves the optimization efficiency using high-fidelity analyses. In this study, the optimized design improves the coaptation area by 40% and the orifice area by 3% in comparison to the baseline design.