A data scientist and engineer detail building a Navier-Stokes solver entirely in Python, demystifying complex CFD software. The project, initially a personal challenge and course component, translates fundamental fluid dynamics equations into discretized code. Immersing a wing geometry onto a Cartesian grid using the immersed boundary method allows simulation of airflow around a bird's wing profile. The solver tackles the coupled pressure-velocity problem via a Poisson pressure equation and Jacobi iteration, updating velocity fields explicitly. Lift calculation from integrated pressure confirms the solver's physical plausibility, showing pressure dominance over drag by nearly 1000x in air.

While laminar and computationally intensive, the project serves as a foundational step towards understanding CFD principles and potential future enhancements like turbulence modeling or Finite Volume methods.