Advanced Numerical Project
Advanced Numerical Project
For each of the following four differential equations, state whether is a PDE or an ODE,
linear, or non – linear, homogeneous, or non – homogeneous, and state what order the
equation is.
𝛛𝛓 𝛛𝛓 𝛛𝟐 𝛓
+ + + 𝛑𝛙 + 𝛓 ∗ (𝛓 + 𝐗) = 𝛓
𝛛𝐗 𝛛𝛙 𝛛𝛑𝟐
𝐚′′ 𝐚
+ 𝐬𝐢𝐧(𝐚𝟐 ) + 𝐛 = 𝐜𝐨𝐬(𝐛) +
𝐛 𝐛
𝐘
𝐘ðð + 𝐘𝛕 =
𝛔𝟐
∆𝐪𝐞 + ∆ ∗ 𝐘 𝟐 + ∆ 𝐥𝐨𝐠(𝐯𝐪𝟑 )𝟐 = 𝟎
Question 2
Classify each of the following four 2nd order PDEs as elliptic, parabolic, or hyperbolic.
x 2
−3( −1)
ux (t, 0) = 0, ux (t, 1) = 0, u(0, x) = e 0.2 , ut (0, x) =0
Use COMSOL with 50 elements and a time step of 0.1 to solve the given equation.
Question 4
For the following rectangular domain, solve Laplace’s equation using 5x7 mesh points
(including boundaries) as shown. The boundary conditions are given.
b) Start with a guess of 0 at all points at which values are unknown and use Gauss –
Seidel iteration to find the solution accurate at all locations within five significant
digits. At which point is convergence the most rapid?
c) Repeat ‘b’ using SOR. Use relaxation factor of 1.23. Compare the number of
iterations to reach the desired tolerance with the Gauss – Seidel method.
Question 5
a) Solve the system using Explicit Euler with h = 0.00001, 0.0001, and 0.01 for 0 ≤
x ≤ 0.5.
b) Solve the system using Implicit Euler with h = 0.00001, 0.0001, and 0.01 for 0 ≤
x ≤ 0.5.
c) Plot the solutions (y1 vs. x) and (y2 vs. x) in a and b on the same graph and make
your comments on the behavior of the system, accuracy and stability.