CFD101: 2D Lid Driven Cavity Flow

This repository provides MATLAB code for the lid-driven cavity flow where incompressible Navier Stokes equation is numerically solved using

2D Lid Driven Cavity Flow

This repo provides a MATLAB example code for the lid-driven cavity flow where incompressible Navier Stokes equation is numerically solved using a simple 2nd order finite difference scheme on a staggered grid system.

(Left: Re = 100, Right: Re = 10,000)

The arrow denotes the velocity field, and the contour denotes its magnitude.

Part 1: Getting Started with the Cavity Flow

• Click here for detailed documentation in English.
• 日本語のドキュメントはこちら から

The numerical scheme is kept primitive; the explicit treatment of viscous term (the solution diverges at low Reynolds number), and the time integration is Euler.

まずは単純な手法でキャビティ流れのシミュレーションを実施します。

Part 2: Implicit Scheme for the Viscous Terms

• Click here for detailed documentation in English.
• 日本語のドキュメントはこちら から

The implicit treatments for viscous terms are implemented, namely the Crank-Nicolson method. For better stability for non-linear terms, Adams-Bashforth, and 3 steps-Runge-Kutta is also implemented.

拡散項に対して陰解法を実装しました。対流項へアダムス・バッシュフォースを使用したもの、3段階のルンゲクッタ法の２つの時間発展を実装しています。

Part 3: Performance Comparison of the Implicit Methods

• Click here for detailed documentation in English.
• 日本語のドキュメントはこちら から

The implicit treatments for viscous terms results in solving the discretized Helmholtz equation at every time step. We compare the performance of five methods.

拡散項に対して陰解法を実装すると離散ヘルムホルツ方程式を毎時間ステップ解く必要が出てきます。直接法、反復法、離散サイン変換を使用した方法などの処理速度を検証します。

Part 4: Validation of the Numerical Scheme

• Click [here] for detailed documentation in English. (not ready)
• 日本語のドキュメントはこちら から

The results of spatial and temporal convergence tests are shown. Convergence tests are run using the method of manufactured solutions where the Navier-Stokes equations are forced so that the solution will be a prescribed time-dependent function.

ある外力項を加えた Navier-Stokes 方程式の数値解と解析解と比較することで、時間積分の精度（1次精度） と空間微分の精度（2次精度）を確認します。

Next to come

The plan is to allow arbitrary boundary conditions for more fun simulations.

Environment

• MATLAB R2019b-R2023b
• Signal Processing Toolbox if you use dct in solving Poisson eqn.

ToDo

1. Implement implicit treatment of viscous terms
2. Implement crank-Nicolson for the non-linear terms
3. Allow obstacles within the domain
4. Allow inflow from the wall
5. Make it to 3D

--

Copyright (c) 2020-2023, The MathWorks, Inc.