Sokolmish / heated_plate

Small research in numerical methods for partial differential equations solving

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Static heat distribution problem

General problem:

$$ \dfrac{\partial}{\partial x}\left(k(x, z)\dfrac{\partial u}{\partial x}\right) + \dfrac{\partial}{\partial z}\left(k(x, z)\dfrac{\partial u}{\partial z}\right) = -f(x, z) $$

Simplified problem:

$$ \dfrac{\partial^2 u}{\partial x^2} + \dfrac{\partial^2 u}{\partial z^2} = - \dfrac{f(x, z)}{k} $$

Finite differences methods:

  • Seidel
  • Thomas (tridiagonal) - faster and more accurate

Monte-carlo methods:

  • Fixed random walk - works fine
  • Float random walk - fast but innacurate near heating source
  • Semi-float random walk - works fine with step size $\approx 1$

Finite differences solution example

About

Small research in numerical methods for partial differential equations solving


Languages

Language:Jupyter Notebook 87.7%Language:Python 12.3%