MATLAB Code for 3D Steady Micropolar Couette Flow: Numerical Solver with Verification and Parametric Study
Description
This dataset provides a complete, self-contained MATLAB implementation for simulating three-dimensional steady micropolar Couette flow between parallel plates. The code solves the full nonlinear system of micropolar fluid equations—including conservation of linear momentum, angular momentum, and incompressibility—using a semi-implicit finite-difference scheme with upwind convection treatment and pressure correction via a Poisson solver. The model incorporates realistic boundary conditions: the upper plate moves with constant velocity while all other walls are stationary and no-slip. The package includes two main components: (1) a class-based solver (MicropolarCouette3D) that handles grid setup, time-stepping, boundary conditions, and physical property definitions; and (2) a driver script that performs numerical verification against the Newtonian limit (N=0) and conducts a detailed parametric study of the micropolarity parameter N, microstructure length scale m, and Reynolds number Re. Key outputs include velocity profiles, microrotation fields, energy dissipation distributions, convergence histories, and integral dissipation metrics. The analysis demonstrates how micropolar effects alter flow structure, enhance local rotation, and increase total energy dissipation compared to classical fluids. All visualizations are static, high-quality, and suitable for publication. This code is intended for researchers in theoretical and computational fluid dynamics, particularly those studying complex fluids with microstructure (e.g., liquid crystals, ferrofluids, blood analogs). It emphasizes physical fidelity, numerical robustness, and reproducibility—without analytical simplifications.
Files
Steps to reproduce
Ensure MATLAB R2023a or later is installed. Place MicropolarCouette3D.m and run_verification_and_analysis.m in the same directory. Run run_verification_and_analysis.m. The script will: Solve the Newtonian (N=0) and micropolar (N>0) cases for verification. Generate four verification plots (convergence, velocity profiles, microrotation magnitude, dissipation field). Perform parametric sweeps over N, m, and Re. Produce four analysis plots showing profile evolution and integral dissipation trends. Print quantitative results to the command window. All figures and data are generated in-memory; no external dependencies are required.
Institutions
- Samara State Technical University
- Ural'skij federal'nyj universitet imeni pervogo Prezidenta Rossii B N El'cina