1
If you want to use the coefficient of the form
$$
\alpha = \begin{cases}
\alpha_1 & r < 0.5 \\
\alpha_2 & r > 0.5
\end{cases}
$$
it is useful to work in polar coordinates. You build two pieces of the solution
$$
u_1(r,\theta), \qquad r < 0.5
$$
and
$$
u_2(r,\theta), \qquad r > 0.5
$$
Then at the interface you ensure solution and flux ...
1
If you want to know the why and the how, you probably want to watch lecture 21.65 here:
https://www.math.colostate.edu/~bangerth/videos.html
Which you probably want to do after you've watched lecture 21.6.
Disclaimer: I'm the one in these videos.
Only top voted, non community-wiki answers of a minimum length are eligible
Related Tags
poisson × 113finite-element × 29
finite-difference × 26
boundary-conditions × 20
pde × 14
elliptic-pde × 14
matlab × 9
linear-solver × 9
convergence × 9
discretization × 9
python × 8
numerics × 8
iterative-method × 8
finite-volume × 8
fourier-analysis × 6
electromagnetism × 6
linear-algebra × 5
fluid-dynamics × 5
multigrid × 5
nonlinear-equations × 4
sparse-matrix × 4
conjugate-gradient × 4
spectral-method × 4
error-estimation × 3
fenics × 3