3D Micropolar Fluid Flow Solver: Numerical Simulation of a Lid-Driven Cavity using MATLAB

This dataset contains a complete MATLAB implementation for the numerical simulation of a three-dimensional, steady-state flow of an incompressible micropolar fluid in a closed cubic cavity with a moving upper lid. The code solves the full system of Eringen's micropolar hydrodynamics equations using a finite difference method on a uniform staggered grid combined with a projection scheme for pressure-velocity coupling. The solver features explicit time integration, second-order central differencing for spatial derivatives, and a Successive Over-Relaxation (SOR) method for the Poisson pressure equation. The simulation accounts for key physical parameters: Reynolds number (Re), and micropolarity parameters (N, m). The package includes functions for applying boundary conditions (no-slip for velocity, homogeneous Neumann for microrotation on stationary walls, and zero microrotation on the moving lid), calculating energy dissipation, and generating publication-quality visualizations of velocity, pressure, microrotation, and dissipation fields. The code is designed for reproducibility and includes detailed convergence monitoring, grid independence verification, and physical consistency checks (e.g., mass conservation). This work provides a foundational numerical framework for researchers and students in computational fluid dynamics studying the influence of internal microstructure on macroscopic flow characteristics in confined domains.