MATLAB Code for Parametric Analysis of Generalized Couette Flow in Micropolar Fluids: Reynolds, Taylor and Micropolarity Parameter Study
Description
This MATLAB R2023b code performs a systematic parametric analysis of generalized Couette flow in a micropolar fluid with couple stresses. The study investigates the influence of three key dimensionless parameters: the micropolarity parameter s (ranging from 0.03 to 5.0), the Reynolds number Re (30 to 250), and the Taylor number Ta (10 to 250). The implementation is based on exact analytical solutions derived via the Lin–Sidorov–Aristov ansatz. The package includes comprehensive parameter scanning across extended ranges and generates eight types of publication-ready visualizations: three-dimensional surfaces of maximum velocity, heat maps of velocity and its gradients, derivative-based curvature analysis, full-field velocity isolines, comparative deviation metrics between classical and micropolar solutions, contour maps of maximum deviations, flow uniformity optimization analysis, and complex deviation surfaces in the multi-parameter space. This dataset provides a complete computational framework for researchers studying non-Newtonian and microstructured fluids in microfluidic, biomedical, and industrial applications.
Files
Steps to reproduce
This dataset contains a complete MATLAB implementation for reproducing the parametric analysis of generalized Couette flow in a micropolar fluid with couple stresses under different boundary conditions. To reproduce all results and figures, follow these steps: 1. Ensure MATLAB R2023b (or a compatible version) is installed. No additional toolboxes are required. 2. Download and place the file "couette_micropolar_parametric_analysis.m" in a dedicated working directory. 3. Open MATLAB and navigate to this directory. 4. Run the script by typing "couette_micropolar_parametric_analysis" in the Command Window and pressing Enter. 5. The script will automatically: - Compute exact analytical solutions for both boundary value problems (vanishing bending moments and vanishing shear stresses), - Perform a systematic scan over the micropolarity parameter (s ∈ [0.03, 5.0]), Reynolds number (Re ∈ [30, 250]), and Taylor number (Ta ∈ [10, 250]), - Generate 8 publication-ready figures illustrating velocity profiles, gradients, 3D surfaces, heat maps, deviation metrics, and flow quality optimization, - Display real-time progress messages in the Command Window. 6. Upon completion, the following output files will be saved in the same directory: - 8 PNG files (e.g., "Fig1_3D_Re_Ta_influence.png", "Fig2_parameter_heatmaps.png", etc.), - 8 corresponding FIG files for further editing in MATLAB, - A comprehensive log of parameter ranges and computational verification in the Command Window. The entire execution takes approximately 1–2 minutes on a standard laptop and requires no user interaction. The code is self-contained and uses only core MATLAB functions, ensuring full reproducibility.
Institutions
- Samara State Technical University
- Ural'skij federal'nyj universitet imeni pervogo Prezidenta Rossii B N El'cina