Study guide: Introduction to finite element methods

Loading [MathJax]/extensions/TeX/boldsymbol.js

$\newcommand{\uex}{{u_{\small\mbox{e}}}} \newcommand{\Aex}{{A_{\small\mbox{e}}}} \newcommand{\half}{\frac{1}{2}} \newcommand{\tp}{\thinspace .} \newcommand{\Oof}[1]{\mathcal{O}(#1)} \newcommand{\x}{\boldsymbol{x}} \newcommand{\X}{\boldsymbol{X}} \renewcommand{\u}{\boldsymbol{u}} \renewcommand{\v}{\boldsymbol{v}} \newcommand{\e}{\boldsymbol{e}} \newcommand{\f}{\boldsymbol{f}} \newcommand{\dfc}{\alpha} % diffusion coefficient \newcommand{\Ix}{\mathcal{I}_x} \newcommand{\Iy}{\mathcal{I}_y} \newcommand{\Iz}{\mathcal{I}_z} \newcommand{\If}{\mathcal{I}_s} % for FEM \newcommand{\Ifd}{{I_d}} % for FEM \newcommand{\Ifb}{{I_b}} % for FEM \newcommand{\sequencei}[1]{\left\{ {#1}_i \right\}_{i\in\If}} \newcommand{\basphi}{\varphi} \newcommand{\baspsi}{\psi} \newcommand{\refphi}{\tilde\basphi} \newcommand{\psib}{\boldsymbol{\psi}} \newcommand{\sinL}[1]{\sin\left((#1+1)\pi\frac{x}{L}\right)} \newcommand{\xno}[1]{x_{#1}} \newcommand{\Xno}[1]{X_{(#1)}} \newcommand{\xdno}[1]{\boldsymbol{x}_{#1}} \newcommand{\dX}{\, \mathrm{d}X} \newcommand{\dx}{\, \mathrm{d}x} \newcommand{\ds}{\, \mathrm{d}s}$

Method 2: Use all $\basphi_i$ and insert the Dirichlet condition in the linear system

Now $\baspsi_i=\basphi_i$ , $i=0,\ldots,N=N_n$

$\basphi_N(L)\neq 0$ , so $u'(L)\basphi_N(L)\neq 0$

However, the term $u'(L)\basphi_N(L)$ in $b_N$ will be erased when we insert the Dirichlet value in $b_N=D$

We can therefore forget about the term

$u'(L)\basphi_i(L)$ !

Boundary terms $u'\basphi_i$ at points $\xno{i}$ where Dirichlet values apply can always be forgotten.

$\begin{equation*} u(x) = \sum_{j=0}^{N=N_n} c_j\basphi_j(x) \end{equation*}$

$\begin{equation*} \sum_{j=0}^{N=N_n}\left( \int_0^L \basphi_i'(x)\basphi_j'(x) dx \right)c_j = \int_0^L f(x)\basphi_i(x) dx - C\basphi_i(0) \end{equation*}$

Assemble entries for $i=0,\ldots,N=N_n$ and then modify the last equation to $c_N=D$

Method 2: Use all \basphi_i \basphi_i and insert the Dirichlet condition in the linear system

Method 2: Use all $\basphi_i$ and insert the Dirichlet condition in the linear system