Topics

• Introduction to Tensor Analysis: Definition of cartesian tensors in 3 dimensions. Transformation properties. Contraction of tensors in 3 dimensions.
• Numerical Analysis II:
  • Partial differential equation: Finite difference approximations to partial derivatives (O(h^2)). Solution of one dimensional heat conduction equation by explicit method. Qualitative idea of explicit and implicit methods. Laplace equation (2-d) using standard five point formula, Successive relaxation technique. Solution of 1-d Wave equation. Stability criterion — CFL condition (qualitative).
• Introduction to Scipy:
  • Interpolation Using scipy.interpolate.lagrange
  • Numerical Integration using: scipy.integrate.quad, scipy.integrate.trapz, scipy.integrate.simps
  • Solving first order and 2nd order ODE by scipy.integrate.odeint()
• Dirac-delta function:
  • Numerically handling improper integrals over infinite intervals over the range ∫_a^∞ correct upto given decimal place without using Scipy.
  • Evaluate 1⁄√2πσ²∫e^{-(x-2)2⁄2σ2} (x+3) dx for x = 1, 0.1, 0.01 and show that it tends to 5.
  • ∫_-∞^∞ exp[-(ax²+bx+c)]dx = √^π⁄^aexp(^b2⁄_4a+c).
  • Verifying that ∫_a-x1^a+x2 δ(x-a)f(x)dx = f(a) using different limiting representation of δ(x).
• Special Functions:
  • Use of special functions taken from scipy.special. Plotting and verification of the properties of special functions. Orthogonality relations and recursion relations. (Legendre and Hermite Only)
    • ∫_-∞^∞ P_n(μ)P_m(μ)dμ = δ_nm ²⁄_2n+1.
    • (1-x²)P_n′(x)+(n+1)xP_n(x)=(n+1)P_n+1(x).
• Solution of some basic PDEs:
  • Boundary value problems. Finite discrete method with fixed step sizes. Idea of stability. Application to simple physical problems. Dirichlet Boundary conditions only.
  • Laplace equation ∇²φ=0 on a square grid with specified potential at the boundaries.
  • Wave equation in 1+1 dimension: ∂_t²φ=λ∂_x²φ Vibration of a string with ends fixed with given initial configurations: φ(x,0) and ∂_tφ(x,0).
  • Heat equation in 1+1 dimension, ∂_tφ=α∂_x²φ with specified value of temperature at the boundaries with given initial temperature at the boundaries with given initial temperature profile.
• Fourier Series:
  Evaluate the Fourier coefficients of a given periodic function using scipy.integrate.quad(). Examples: square wave, triangular wave, saw-tooth wave. Plot to see a wave form from scipy.signal and the constructed series along with.
• Curve Fitting:
  Curve Fitting with user defined functions using scipy.optimize module

Study Materials

• Scientific Computing in python - Avijit Kar Gupta.
• Computational Physics problem solving with Computers - Landau, Paez, Bordeianu.
• Programming for Computation-Python - Svein Linge & Hans Petter Lantangen.
• Numerical Solution of Partial Differential Equations in Science and Engineering - L. Lapidus & G.F. Pinder.
• PDE (Schaum Series) - P. Duchateau & D.W. Zachmann.
• Spectral Methods in MATLAB - L.N.Trefethen.
• Limit Representation of Dirac Delta.