Modeling oxygen diffusion within the aortic heart valve leaflet tissue is challenging. This is because the aortic valve cusp is highly heterogeneous with varying physical properties with respect to location and direction. The conventional numerical solutions such as the finite element method, the finite difference method and other methods applicable to modeling the oxygen diffusion within the tissue contain error regardless of the nature of the given problem or modeling parameters. The major concern with numerical methods is computational cost such as the execution-time and the accuracy. This study focuses on modeling oxygen diffusion within the heart valve leaflet tissue using a novel inverse numerical technique. The Lagrange multiplier method and the least-squares algorithm are used to define a cost function which is discretized using the finite difference method and is optimized to reduce erroneous results. While our model proposes accuracy similar to that of an equivalent finite element approach, the execution-time of the solution is reduced. The numerical method developed in this study can have a broad application in tissue engineering/repair and regenerative medicine in which oxygenation of the cells residing within the scaffold/tissue is of particular importance.