File Exchange

image thumbnail

2D Unsteady Navier-Stokes

version 1.3 (446 KB) by Andrea La Spina
Solution of 2D Navier-Stokes equations


Updated 07 Aug 2018

View License

The "UNSTEADY_NAVIER_STOKES" script solves the 2D steady Navier-Stokes equations.
The space discretization is performed by means of the standard Galerkin approach.
Bilinear quadrangular elements are used for the pressure and biquadratic quadrangular elements are used for the velocity.
This pressure and velocity interpolation satisfies the so-called LBB condition, which ensures the solvability of the algebraic system.
For the time integration the theta-method has been implemented. According to the value of theta these schemes are obtained:
0 -> Forward Euler
1/2 -> Crank-Nicolson
3/4 -> Galerkin
1 -> Backward Euler
The FEM parameters such as the number of finite elements and the number of Gauss integration points can be easily chosen.
The functions and the examples are developed according with Chapter 6 "Viscous incompressible flows" of the book "Finite Element Methods for flow problems" of Jean Donea and Antonio Huerta.
If you like the file, please give feedbacks.

Comments and Ratings (2)

Liu Fusheng

very interesting!!!






MATLAB Release Compatibility
Created with R2014a
Compatible with any release
Platform Compatibility
Windows macOS Linux

Discover Live Editor

Create scripts with code, output, and formatted text in a single executable document.

Learn About Live Editor

2D Unsteady Navier-Stokes