https://github.com/takah29/2d-fluid-simulator

2D incompressible fluid solver implemented in Taichi.
https://github.com/takah29/2d-fluid-simulator

cfd cfd-solvers cip fdm finite-difference-method fluid-dynamics fluid-simulation fluid-solver navier-stokes taichi vorticity-confinement

Last synced: about 2 months ago
JSON representation

2D incompressible fluid solver implemented in Taichi.

• Host: GitHub
• URL: https://github.com/takah29/2d-fluid-simulator
• Owner: takah29
• Created: 2022-04-08T15:39:23.000Z (over 2 years ago)
• Default Branch: main
• Last Pushed: 2024-05-10T19:44:18.000Z (8 months ago)
• Last Synced: 2024-08-02T11:23:38.795Z (5 months ago)
• Topics: cfd, cfd-solvers, cip, fdm, finite-difference-method, fluid-dynamics, fluid-simulation, fluid-solver, navier-stokes, taichi, vorticity-confinement
• Language: Python
• Homepage:
• Size: 37.5 MB
• Stars: 93
• Watchers: 0
• Forks: 16
• Open Issues: 0
• Metadata Files:
  • Readme: README.md

Awesome Lists containing this project

• awesome-taichi - 2D Fluid Simulator - A hard-core CFD (Computational Fluid Dynamics) solver powered by Taichi. (Applications / **Simulation**)
• awesome-taichi - 2D Fluid Simulator - A hard-core CFD (Computational Fluid Dynamics) solver powered by Taichi. (**Simulation**)

README

        
# 2D Fluid Simulator

![baundary_condition_2_dye](./images/bc2_res1600_cip_dye.jpg)

![baundary_condition_2_norm](./images/bc2_res1600_cip_norm.jpg)

## Features

- Finite Difference Method (MAC Method)

- Advection Scheme

  - Upwind Differencing

  - Kawamura-Kuwahara

  - CIP (Constrained Interpolation Profile)

- Flow Visualization

  - Norm (Velocity) and Pressure

  - Pressure

  - Vorticity

  - Dye

- Vorticity Confinement

## Requirements

- Python 3.11

- Taichi 1.7

GeForce GTX 1080 or higher recommended.

## Usage

- Boundary Condition 1, ReynoldsNumber = 1000, dt = 0.0005, VorticityConfinement is Disable

  ```bash

  python main.py -re 1000 -dt 0.0005 -vc 0.0

  ```

  Press `V` key switches the flow visualization method.

  `dt` is automatically determined even if not specified, but should be small for divergence.

- Boundary Condition 2, resolution = 800

  ```bash

  python main.py -bc 2 -res 800

  ```

  Boundary conditions can be specified from 1 to 6

- Boundary Condition 3, ReynoldsNumber = 10^8, resolution = 800, VorticityConfinement = 10

  ```bash

  python main.py -bc 3 -re 100000000 -res 800 -vc 10

  ```

- Help

  ```bash

  python main.py -h

  ```

- for CPU

  ```bash

  python main.py -dt 0.0005 -cpu

  ```

## Screenshots

### Flow Visualization

- Norm and Pressure

  ![norm_and_pressure](./images/bc5_res800_cip_norm.jpg)

- Pressure

  ![pressure](./images/bc5_res800_cip_pressure.jpg)

- Vorticity

  ![vorticity](./images/bc5_res800_cip_vorticity.jpg)

- Dye

  ![dye](./images/bc5_res800_cip_dye.jpg)

### Vorticity Confinement

- Disable

  ![no_vorticity_confinement](./images/bc3_res800_cip_dye_novc.jpg)

- Enable

  ![vorticity_confinement](./images/bc3_res800_cip_dye_vc.jpg)

## References

- [移流法](https://pbcglab.jp/cgi-bin/wiki/index.php?%E7%A7%BB%E6%B5%81%E6%B3%95)

- [2 次元 CIP 法による移流項の計算](https://i-ric.org/yasu/nbook2/04_Chapt04.html#cip)

- [GPU Gems Chapter 38. Fast Fluid Dynamics Simulation on the GPU

  ](https://developer.nvidia.com/gpugems/gpugems/part-vi-beyond-triangles/chapter-38-fast-fluid-dynamics-simulation-gpu)

- [Ronald Fedkiw, Jos Stam, Henrik Wann Jensen. Visual Simulation of Smoke.](https://web.stanford.edu/class/cs237d/smoke.pdf)