MATLAB Implementation of a Projection Method for 2D Micropolar Couette Flow
Description
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.
Files
Steps to reproduce
Ensure MATLAB is installed (version R2018a or later recommended). Download the two MATLAB files: run_micropolar_channel.m and MicropolarFlow2D_TopWall_Fixed.m. Place both files in the same directory. Open MATLAB and navigate to that directory. Run the main script by typing run_micropolar_channel in the Command Window. The simulation will execute and automatically generate two figures: figure1_main.png – Primary flow fields and convergence. figure2_additional.png – Vorticity and dissipation contours. To modify parameters (e.g., Reynolds number, micropolar parameters), edit the initial lines of run_micropolar_channel.m.
Institutions
- Samara State Technical University
- Ural'skij federal'nyj universitet imeni pervogo Prezidenta Rossii B N El'cina