[chandanvarma07]

Hello everyone,

Inspired by the insightful teachings of Prof. Sandip Mazumder on YouTube, I have constructed a CFD code entirely from scratch using MATLAB for a classic problem: the lid-driven cavity flow. My goal is to implement a co-located mesh approach and harness the capabilities of the SIMPLE algorithm to ensure accurate pressure-velocity coupling. But after running the code, the residuals were blowing upto infinity.

In the early stages of debugging, I closely monitored the simulation's behavior. The initial iteration yielded a promising outcome, with both X and Y momentum well-preserved—a strong validation of my approach. Employing the TDMA method, I've attached an x momentum contour from the first iteration for your reference, exhibiting similarities to patterns observed in heat conduction cases solved by the ADI Method.

However, as I delved further into debugging, a significant challenge emerged: accurately calculating face velocities. It became evident that inaccuracies at this step were rippling through the simulation, affecting pressure source term calculations and subsequently impacting the pressure correction data. Notably, certain cells exhibited discrepancies between the sign of pressure source terms and correction data.

Through careful examination, I've begun to explore whether the absence of appropriate boundary conditions in face velocity calculation might contribute to these disparities. This has prompted a deeper consideration of boundary condition implementation and its potential influence on this critical simulation stage.

To provide a comprehensive view of my progress, I've attached the code I've been diligently working on. Your insights and expertise in reviewing this code would be invaluable in helping me address these challenges.

TIA