The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
17.1. FLOW PAST A RECTANGLE
135
is parallel to the xaxis, must be constant. Matching the constant to the value of
on the midline, we obtain = 0 along the boundary of the obstacle.
From the noslip co
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
12.3. PHASE PORTRAITS
12.3
85
Phase portraits
The phase space of a dynamical system consists of the independent dynamical variables. For example, the phase space of the damped, driven pendulum equatio
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
115
The third part of this course considers a problem in computational uid dynamics
(cfd). Namely, we consider the steady twodimensional ow past a rectangle or a circle.
116
Chapter 14
Derivation of
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
8.3. MULTIDIMENSIONAL INTERPOLATION
55
From the known value of , and using (8.6), (8.13), and (8.14), we have
= 3an
=3
2
1 hn 1
1
3hn
+ 2bn
(bn
1 hn 1
bn
1
1)
+ cn
h2 1
n
1
+ 2bn
1 hn 1
+
fn
hn
1
hn
3
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
Chapter 16
Stream function, vorticity
equations
16.1
Stream function
A streamline at time t is dened as the curve whose tangent is everywhere parallel to
the velocity vector. With dx along the tangent
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
10.2. PERIOD OF MOTION
65
the real parts and the imaginary parts of (10.6), and we nd
r
T
,
m
g sin .
r2 = g cos
r + 2r =
(10.7)
(10.8)
If the connector is a rigid rod, as we initially assumed, then
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
5.2. LU DECOMPOSITION
35
where L is a lower triangular matrix and U is an upper triangular matrix.
Using the same matrix A as in the last section, we show how this factorization
is realized. We have
0
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
13.5. PERIOD DOUBLING IN THE LOGISTIC MAP
105
A period2 cycle of the map f (x) necessarily corresponds to two distinct xed points
of the composite map (13.7). We denote these two xed points by x0 and
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
Chapter 1
IEEE arithmetic
1.1
Denitions
The real line is continuous, while computer numbers are discrete. We learn
here how numbers are represented in the computer, and how this can lead to
roundo er
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
Chapter 4
Dierential equations
We now discuss the numerical solution of ordinary dierential equations. We
will include the initial value problem and the boundary value problem.
4.1
Initial value probl
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
2.5. ORDER OF CONVERGENCE
15
or
n+1 = n +
f (xn )
.
f 0 (xn )
(2.1)
We use Taylor series to expand the functions f (xn ) and f 0 (xn ) about the root
r, using f (r) = 0. We have
1
r)f 0 (r) + (xn
2
f
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
Chapter 7
Iterative methods
7.1
Jacobi, GaussSeidel and SOR methods
Iterative methods are often used for solving a system of nonlinear equations.
Even for linear systems, iterative methods have some
The Hong Kong University of Science and Technology
Scientific Computing
MATH 164

Fall 2013
11.2. THE NONLINEAR PENDULUM
75
t, and reuse some of the other parameter names, with the understanding that the
damped, driven pendulum equation that we will now study numerically is dimensionless. We