Science Score: 44.0%
This score indicates how likely this project is to be science-related based on various indicators:
-
✓CITATION.cff file
Found CITATION.cff file -
✓codemeta.json file
Found codemeta.json file -
✓.zenodo.json file
Found .zenodo.json file -
○DOI references
-
○Academic publication links
-
○Academic email domains
-
○Institutional organization owner
-
○JOSS paper metadata
-
○Scientific vocabulary similarity
Low similarity (11.8%) to scientific vocabulary
Repository
Stokes solver
Basic Info
- Host: GitHub
- Owner: skjelanger
- Language: Python
- Default Branch: main
- Size: 254 KB
Statistics
- Stars: 0
- Watchers: 1
- Forks: 0
- Open Issues: 0
- Releases: 1
Metadata Files
README.md
FEM Stokes Flow Solver with Periodic Boundaries and Adaptive Meshing
This project simulates 2D incompressible Stokes flow through a square domain with a circular obstacle. The solver is based on P2-P1 mixed finite elements, periodic boundary conditions, and is designed to assess permeability under varying geometries.
Mesh generation and classification are handled via Gmsh and Meshio, while linear systems are solved with Pypardiso for efficiency.
Features
- Quadratic P2 elements for velocity and linear P1 elements for pressure
- Support for periodic boundary conditions (top-bottom and left-right)
- Adaptive meshing around obstacles and domain walls
- Batch simulation support with logging and result saving
- Permeability calculation via Darcy's law
- Mesh and result visualization (pressure, velocity fields, vector plots)
Dependencies
Install the following Python packages (e.g., via pip or conda):
gmshmeshionumpyscipymatplotlibpypardiso(or switch to another solver manually if needed)
Running the Stokes Solver
To execute a simulation:
- Download stokes_solver.py and mesher.py.
- Install the required dependencies listed above.
- Run the main solver script:
bash python stokes_solver.py bash
This will:
- Generate a mesh using Gmsh.
- Assemble and solve the Stokes system using P2-P1 mixed finite elements.
- Apply periodic or inflow/outflow boundary conditions (as configured).
- Plot and save the velocity magnitude, velocity vectors, and pressure fields.
- Compute the permeability.
- Save the entire simulation state (solution + mesh + metadata) as a compressed file under simulations/.
Output
- Figures are saved to the plots/ folder as .png files.
- Full simulation results are saved as .pkl.gz files under simulations/.
You can modify mesh parameters and boundary settings directly in the main section at the bottom of stokes_solver.py.
Notes
Periodic boundary conditions are applied at the matrix level via slave-master DOF enforcement.
Meshing can be adapted by adjusting create_mesh() in mesher.py.
Pressure is anchored to one node to ensure solvability.
Acknowledgements
Created during the PCS911 course at Western Norway University of Applied Sciences
Owner
- Name: Skjelanger
- Login: skjelanger
- Kind: user
- Location: Bergen
- Repositories: 1
- Profile: https://github.com/skjelanger
Master degree in Theoretical Physics.
Citation (CITATION.cff)
cff-version: 1.2.0
message: "If you use this software, please cite it as below."
title: "Finite Element Solver in Python for 2D Stokes flow"
version: "1.0.0"
authors:
- family-names: "Skjelanger"
given-names: "Kristoffer"
date-released: 2025-05-29
repository-code: "https://github.com/skjelanger/stokes_solver"
license: "MIT"
keywords:
- stokes flow
- CFD
- numerical methods
- finite element
GitHub Events
Total
- Release event: 1
- Push event: 4
- Create event: 1
Last Year
- Release event: 1
- Push event: 4
- Create event: 1
Dependencies
- gmsh *
- matplotlib *
- meshio *
- numpy *
- scipy *