MATLAB Implementation of a Projection Method for 2D Micropolar Couette Flow

This dataset contains MATLAB code for simulating the flow of a micropolar fluid in a two-dimensional rectangular channel. The upper wall of the channel moves with a constant velocity, inducing shear-driven flow, while the lower and side walls remain stationary. The simulation implements a finite-difference-based numerical solver that couples the linear momentum equations with the angular momentum equation governing micropolar fluid behavior. Key features of the code include:

Implementation of the micropolar fluid model, which accounts for microscopic rotational degrees of freedom and couple stresses.

Use of the SIMPLE-like algorithm for pressure-velocity coupling.

Upwind differencing for convective terms to enhance numerical stability.

Dynamic relaxation and adaptive time-stepping for convergence control.

Calculation of derived quantities such as vorticity and dissipation.

The code is designed for research and educational purposes in the field of non-Newtonian fluid mechanics. It provides a clear framework for studying the effects of micropolar parameters (such as the coupling number N and characteristic length m) on flow patterns, velocity profiles, and microrotation fields. The results are visualized through vector plots, streamlines, contour maps, and convergence histories.

Typical applications include the analysis of complex fluids like liquid crystals, polymeric suspensions, blood flow, and ferrofluids. The dataset is accompanied by example results and is ready to run in MATLAB without requiring additional toolboxes.