1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
|
# Partial Differential Equation
Partial differential equations are defines when two or more partial derivatives are present in an equation. Due to the widespread application in engineering, we will be looking at **second-order equations** which can be expressed in the following general form:
$$
A\frac{\partial^2u}{\partial x^2} +
B\frac{\partial^2u}{\partial x \partial y} +
C\frac{\partial^2u}{\partial y^2} +
D
= 0
$$
where $A$, $B$ and $C$ are functions of both $x$ and $y$. $D$ is a a function of $x$, $y$, $u$, $\partial u / \partial x$, and $\partial u / \partial y$. Similar to our beloved quadratic formula, we can take the discriminant of the equation
$$
\Delta = B^2 - 4 AC
$$
Based on the discriminant we can categorize the equations into the following three categories:
| $\Delta$ | Category | Example |
| -------- | ---------- | ------------------------ |
| - | Elliptical | Laplace equation |
| 0 | Parabolic | Heat Conduction equation |
| + | Hyperbolic | Wave equation |
## Finite Difference Methods
### Elliptic Equations
- Used for steady-state, boundary value problems
- Examples where fields where these equations are used are: steady-state heat conduction, electrostatics and potential flow
Description of how the Laplace equations works
$$
\frac{\partial^2T}{\partial x^2}+\frac{\partial^2T}{\partial y^2}=0
$$
Finite-different solutions
- Laplacian Difference equations in dimension $x$ and $y$:
$$
\frac{\partial^2T}{\partial x^2}= \frac{T_{i+1,j}-2T_{i,j}+T_{i-1,j}}{\Delta x^2}
$$
$$
\frac{\partial^2T}{\partial y^2}= \frac{T_{i+1,j}-2T_{i,j}+T_{i-1,j}}{\Delta y^2}
$$
Boundary Conditions
Control-Volume approach
<img
style="display: block;
margin-left: auto;
margin-right: auto;
width: 50%;"
src="control_volume_approach.png"
alt="Two Different perspsectives for developing approximate solutions of PDE: (a) finite-difference or node and (b) control-volume based">
Computer Algorithms
### Parabolic Equations
- Used for unstead-state, initial + boundary conditions problems
For parabolic PDE equations we also consider the change in time as well as space.
Heat-conduction equation
Explanation of heat equation
$$
k\frac{\partial^2T}{\partial x^2}=\frac{\partial T}{\partial t}
$$
Explicit methods
Simple Implicit methods
Crank-Nicolson
ADI
### Hyperbolic Equations
MacCormack Method
## Finite-Element Method
General Approach
1. Discretization
2. Element Equations
3. Assembly
4. Boundary Conditions
5. Solutions
6. Postprocessing
### One-dimensional analysis
### Two-dimensional Analysis
# Problem 1: Finite-Element Solution of a Series of Springs
Problem 32.4 from Numerical Methods for Engineers 7th Edition Steven C. Chapra and Raymond P. Canale
A series of interconnected strings are connected to a fixed wall where the other is subject to a constant force F. Using the step-by-step procedure from above, determine the displacement of the springs.
## Problem 2: Finite
Solve the non-dimensional transient heat conduction equation in two dimensions, which represents the transient temperature distribution in an insulated plate
|
