https://github.com/ingeborggjerde/stokes-paradox
https://github.com/ingeborggjerde/stokes-paradox
Last synced: about 1 month ago
JSON representation
- Host: GitHub
- URL: https://github.com/ingeborggjerde/stokes-paradox
- Owner: IngeborgGjerde
- Created: 2024-01-02T19:47:22.000Z (over 1 year ago)
- Default Branch: main
- Last Pushed: 2024-01-07T19:30:54.000Z (over 1 year ago)
- Last Synced: 2025-04-12T10:14:18.937Z (about 1 month ago)
- Language: Python
- Size: 15.6 KB
- Stars: 0
- Watchers: 1
- Forks: 0
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
Awesome Lists containing this project
README
# New insights on the Stokes paradox
This repo contains code used to generate the figures and results in [*New Insights on the Stokes Paradox for Flow in Unbounded Domains*]([https://arxiv.org/pdf/2301.00039.pdf](https://trebuchet.public.springernature.app/get_content/779d0e10-e068-464d-8230-0e72eadc58c4?utm_source=rct_congratemailt&utm_medium=email&utm_campaign=nonoa_20240106&utm_content=10.1140/epjp/s13360-023-04804-6)) by Ingeborg Gjerde and Ridgway Scott.
## Dependencies
Simulations were performed using `FEniCS`, with [`mshr`](https://bitbucket.org/fenics-project/mshr/) for meshing and `octave`/`matplotlib` for plotting.
## Install
To run the FEniCS scripts
```
conda create -n fenicsproject -c conda-forge fenics mshr
source activate fenicsproject
```
## Usage
For the results in Figure 2, we refer to the [repo](https://github.com/IngeborgGjerde/nitsche-method-for-navier-stokes-with-slip) on Nitsche's method for slip boundary conditions. A short example is given in `stokes_with_friction.py`.
The simulations for Figure 4 were generated using `Figure4/stokes_noslip.py` for the no-slip boundary conditions,
and `Figure4/stokes_freeslip.py` for the slip boundary conditions with zero friction. The results are collected and plotted in the octave scripts`Figure4/plot_results_stokes_freeslip.m` and `Figure4/plot_results_stokes_freeslip.m`, respectively.
Figure 5 was generated using `ploteffpee.m`. It requires an input "tol" which is an input to the code `effprimestopar.m` which computes the data for the individual curves in Figure 5.
The other input to `effprimestopar.m` is the value "b" indicated in Figure 5. The input "tol" is used to control the number of data points in the plot.