TEP4280: Introduction to Computational Fluid Dynamics
# Model Equations
## Burger's Equation
$$\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} = \nu \frac{\partial ^2 u}{\partial x^2} $$
The inviscid version of the equation is
$$\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} = 0$$
The linear version of the Burger's Equation is often called the Convection–Diffusion equation
$$\frac{\partial u}{\partial t} + u_0 \frac{\partial u}{\partial x} = \alpha \frac{\partial ^2 u}{\partial x^2} $$
### Numerical Schemes
### Application
## Diffusion Equation
$$\frac{\partial u}{\partial t} = \alpha \frac{\partial ^2 u}{\partial x^2} $$
### Numerical Schemes
### Application
The diffusion equation is applied when a quantity $u$ is being distributed throughout a space or body with out moving the initial source of the quantity.
#### Heat Conduction
$$\frac{\partial T}{\partial t} = \alpha \frac{\partial ^2 T}{\partial x^2} $$
where $T$ is the temperature, and $alpha$ is the heat conduction constant.
#### Flow in porous media
$$\frac{\partial u}{\partial t} = c \frac{\partial ^2 u}{\partial x^2} $$
## Poisson Equation
$$\frac{\partial ^2 u}{\partial x^2} + \frac{\partial ^2 u}{\partial y^2} = f(x, y)$$
Setting $f(x, y) = 0 $ will give the Laplace equation.
### Numerical Schemes
### Application
## Wave Equation
$$\frac{\partial ^2 u}{\partial t^2} = \alpha_0 ^2 \frac{\partial ^2 u}{\partial x^2}$$
The linear version of the eqaution is
$$\frac{\partial u}{\partial t} + \alpha_0 \frac{\partial u}{\partial x} = 0$$
### Numerical Schemes
### Application
# Boundary Conditions
# Numerical Methods
## Euler
$$ y_{n+1} = y_n + h \cdot f(t_n, y_n)$$
### Errors
$$ e_{n+1} = \mathcal{O} (h^2)$$
$$ E_n = \mathcal{O} (h)$$
### Stability Region
$$ \{z \in \mathbb{C} |\quad |1 + z| \leq 1\}$$
## Implicitt Euler
$$ y_{n+1} = y_n + h \cdot f(t_{n+1}, y_{n+1})$$
### Errors
$$ E_n = \mathcal{O} (h)$$
### Stability Region
$$ \{z \in \mathbb{C} |\quad |\frac{1}{1- z}| \leq 1\}$$
## Trapezoidal
$$y_{n+1} = y_n + h\cdot \frac{1}{2}\cdot [f(x_n, y_n) +(f(x_{n+1}, y_{n+1})] $$
### Errors
$$ e_{n+1} = \mathcal{O} (h^3)$$
$$ E_n = \mathcal{O} (h^2)$$
### Stability region
$$ \{z \in \mathbb{C} | Re(z) < 0\}$$
## Heun
$$y_{n+1} = y_n + h\cdot \frac{1}{2}\cdot [f(x_n, y_n) +(f(x_{n+1}, y*_{n+1})] $$
where $y*_{n+1} = y_n + h\cdot f(x_n, y_n) $
### Errors
$$ e_{n+1} = \mathcal{O} (h^3)$$
$$ E_n = \mathcal{O} (h^2)$$
### Stability Region
$$ \{z \in \mathbb{C} |\quad |1 + z + \frac{z^2}{2}| \leq 1\}$$
## Runge-Kutta
$$ k_1 = f(t_n, y_n)$$
$$ k_2 = f(t_n + \frac{h}{2}, y_n + \frac{h}{2}\cdot k_1)$$
$$ k_3 = f(t_n + \frac{h}{2}, y_n + \frac{h}{2}\cdot k_2)$$
$$ k_2 = f(t_n + h, y_n + h\cdot k_3)$$
$$ y_{n+1} = y_n + \frac{h}{6} \cdot (k_1 + 2k_2 + 2k_3 + k4)$$
### Errors
$$ e_{n+1} = \mathcal{O} (h^5)$$
$$ E_n = \mathcal{O} (h^4)$$
### Stability Region
$$ \{z \in \mathbb{C} |\quad |1 + z + \frac{z^2}{2}+ \frac{z^3}{6}+ \frac{z^4}{24}| \leq 1\}$$
# Discretization methods
### FPM
### FVM
### FEM
### FTCS
### Spectral method
### Spectral element methods
# Termonology