View source

You do not have permission to edit this page, for the following reason:

You can view and copy the source of this page:

:The matrix form of the residual formulation can be written as: <math>\mathbf{K} \cdot \mathbf{a} = \mathbf{f} </math> with: ::<math> \begin{cases} K_{ij} & = \int_{\Omega} W_i(x) \left[ \frac{d}{dx} k \frac{d}{dx} \left( N_j (x) \right) \right] d\Omega \\ f_i & = \frac{}{} \int_{\Omega} W_i(x) Q(x) d\Omega \end{cases} </math> :For simplicity reasons, we select as basis functions the Fourier Series <math>(N_i(x)=\sin \frac {\pi x i}{l} \,)</math>, that automatically accomplish the boundary conditions <math>(N_i(0)=0 \,</math> and <math>N_i(x_L)=0 \,)</math>. If the material is constant, we can write: ::<math> \frac{d}{dx} k \frac{d}{dx} \left ( \sin \frac {\pi x i}{l} \right) = \left ( - k \left (\frac{\pi i}{l} \right )^2 \sin \frac {\pi x i}{l} \right) = - k \left (\frac{\pi i}{l} \right )^2 \sin \frac {\pi x i}{l} </math> ::<math> \begin{cases} K_{ij} & = - \int_{\Omega} W_i(x) k \left (\frac{\pi j}{l} \right )^2 \sin \frac {\pi x j}{l} d\Omega \\ f_i & = - \frac{}{} \int_{\Omega} W_i(x) Q(x) d\Omega \end{cases} \Rightarrow \begin{cases} K_{ij} & = \left (\frac{\pi j}{l} \right )^2 \int_{\Omega} W_i(x) k \sin \frac {\pi x j}{l} d\Omega \\ f_i & = \frac{}{} \int_{\Omega} W_i(x) Q(x) d\Omega \end{cases} </math> :{| style="background:Silver; color:green" |- ! sources ! medium ! boundary conditions ! shape functions |- style="background:Beige; color:green" | :''constant source'' | :''homogeneous'' | :''Dirichlet'' | :''Fourier Series'' |- style="background:Ivory; color:green" | :<math> Q(x) = \begin{cases} Q_0 & x_0 \leqslant x \leqslant x_f \\ 0 & x_f \leqslant x \leqslant x_L \end{cases} </math> | :<math>k(x) = cte \frac{ }{ }</math> | :<math>\varphi (x_0)= \varphi_0 </math> :<math>\varphi (x_L)= \varphi_L </math> | :<math>N_i(x)=\sin \frac {\pi x i}{l} \,</math> |} :By using different weighting functions, a new range of methods can be classified: == Collocation Method &nbsp; &nbsp; <math>(\frac{}{} W_i(x)=\delta(x - x_i),</math> <math>\frac{}{} i=1, 2,... , n)</math> == {| ! General expression ! k and Q constant, using Fourier series |- | ::<math> \begin{cases} K_{ij} & = \frac{d}{dx} k \frac{d}{dx} N_j (x) \bigg|_{x=x_i} \\ f_i & = Q(x_i) \, \end{cases} </math> | ::<math> \begin{cases} K_{ij} & = k \left (\frac{\pi j}{l} \right )^2 \sin \frac {\pi x_i j}{l} \\ f_i & = Q(x_i) \, \end{cases} </math> |} <br> :Number of points = '''1''' ::<math> \begin{cases} x_1 & = \frac{l}{2} \\ Q(x_1) & = \frac{Q}{2} \\ K_{11} & = k \left (\frac{\pi}{l} \right )^2 \sin \frac {\pi x_1}{l} \\ f_1 & = Q(x_1) \, \end{cases} </math> :Number of points = '''2''' ::<math> \begin{cases} x_1 & = \frac{l}{4} \\ x_2 & = \frac{3 l}{4} \\ Q(x_1) & = Q \\ Q(x_2) & = 0 \\ K_{ij} & = k \left (\frac{\pi j}{l} \right )^2 \sin \frac {\pi x_i j}{l} \\ f_i & = Q(x_i) \, \end{cases} </math> :Number of points = '''n''' ::<math> \begin{cases} x_i & = \frac{i l}{n+1} \\ Q(x_i) & = \begin{cases} Q & if ~ x < \frac{l}{2} \\ 0 & if ~ x > \frac{l}{2} \end{cases} \\ K_{ij} & = k \left (\frac{\pi j}{l} \right )^2 \sin \frac {\pi x_i j}{l} \\ f_i & = Q(x_i) \, \end{cases} </math> <br> :[[Image:collocationWRM.jpg]] :''To modify the parameters, edit this'' [[Matlab code Collocation Method of the Resolution of the Poisson's equation with the WRM using global Shape Functions|'''Matlab code''']] == Subdomain Method &nbsp; &nbsp; <math>(W_i(x)=1 \quad \forall x \in \Omega_i \quad and \quad W_i(x)=0 \quad \forall x \notin \Omega_i )</math> == {| ! General expression ! k and Q constant, using Fourier series |- | ::<math> \begin{cases} K_{ij} & = \int_{\Omega_i} \left[ \frac{d}{dx} k \frac{d}{dx} \left( N_j (x) \right) \right] d\Omega \\ f_i & = \frac{}{} \int_{\Omega_i} Q(x) d\Omega \end{cases} </math> | ::<math> \begin{cases} K_{ij} & = k \left (\frac{\pi j}{l} \right )^2 \int_{\Omega_i} \sin \frac {\pi x j}{l} d\Omega = - k \left (\frac{\pi j}{l} \right ) \cos \frac {\pi x j}{l} \Big|_{\Omega_i} \\ f_i & = Q(x) \Big|_{\Omega_i} \int_{\Omega_i} d\Omega \, \end{cases} </math> |} <br> :Number of points = '''1''' ::<math> \begin{cases} x_1 & = \frac{l}{2} \\ Q(x_1) & = \frac{Q}{2} \\ K_{11} & = 2 k \left (\frac{\pi}{l} \right ) \\ f_1 & = Q(x_1) \, \end{cases} </math> :Number of points = '''2''' ::<math> \begin{cases} x_1 & = x_0 \\ x_2 & = \frac{l}{2} \\ x_3 & = x_L \\ Q(x_1) & = Q \\ Q(x_2) & = 0 \\ K_{ij} & = - k \left (\frac{\pi j}{l} \right ) \left ( \cos \frac {\pi x_{i+1} j}{l} - \cos \frac {\pi x_{i} j}{l}\right )\\ f_i & = Q(x_i) (x_{i+1} - x_i)\, \end{cases} </math> :Number of points = '''n''' ::<math> \begin{cases} x_1 & = x_0 \\ x_{i+1} & = \frac{i l}{n} \\ Q(x_i) & = \begin{cases} Q & if ~ x < \frac{l}{2} \\ 0 & if ~ x > \frac{l}{2} \end{cases} \\ K_{ij} & = - k \left (\frac{\pi j}{l} \right ) \left ( \cos \frac {\pi x_{i+1} j}{l} - \cos \frac {\pi x_{i} j}{l}\right )\\ f_i & = Q(x_i) (x_{i+1} - x_i) \, \end{cases} </math> <br> :[[Image:subdomainWRM.jpg]] :''To modify the parameters, edit this'' [[Matlab code Subdomain Method of the Resolution of the Poisson's equation with the WRM using global Shape Functions|'''Matlab code''']] == Galerkin Method &nbsp; &nbsp; <math>(\frac{}{} W_i(x) \equiv N_i(x))</math> == {| ! General expression ! k and Q constant, using Fourier series |- | ::<math> \begin{cases} K_{ij} & = \int_{\Omega} N_i(x) \left[ \frac{d}{dx} k \frac{d}{dx} \left( N_j (x) \right) \right] d\Omega \\ f_i & = \frac{}{} \int_{\Omega} N_i(x) Q(x) d\Omega \end{cases} </math> | ::<math> \begin{cases} K_{ij} & = \int_{\Omega} \sin \frac {\pi x i}{l} k \left (\frac{\pi j}{l} \right )^2 \sin \frac {\pi x j}{l} d\Omega \\ & = k \left (\frac{\pi j}{l} \right )^2 \int_{\Omega} \sin \frac {\pi x i}{l} \sin \frac {\pi x j}{l} d\Omega \\ f_i & = \frac{}{} \int_{\Omega} \sin \frac {\pi x i}{l} Q(x) d\Omega \end{cases} </math> |} <br> :Number of points = '''1''' ::<math> \begin{cases} x_1 & = \frac{l}{2} \\ Q(x_1) & = Q \, \\ K_{11} & = k \left ( \frac{\pi^2}{2 l} \right) \\ f_1 & = \left ( \frac{Q(x_1) l}{\pi} \right)\, \end{cases} </math> :Number of points = '''2''' ::<math> \begin{cases} x_1 & = x_0 \, \\ x_2 & = \frac{l}{2} \\ x_3 & = x_L \, \\ Q(x_1) & = Q \\ Q(x_2) & = 0 \\ K_{ij} & = \begin{cases} k \left (\frac{\pi j}{l} \right )^2 \frac {l}{2} & if ~ i = j \\ k \left (\frac{\pi j}{l} \right )^2 \frac {l}{2 \pi} \left( \frac {sin(\pi (i-j))}{i-j} - \frac {sin(\pi (i+j))}{i+j} \right) & if ~ i <> j \end{cases} \\ f_i & = Q l \frac{1 - cos(\frac{i \pi}{2})}{i \pi} \, \end{cases} </math> :Number of points = '''n''' ::<math> \begin{cases} x_1 & = x_0 \, \\ x_{i+1} & = \frac{i l}{n} \\ Q(x_i) & = \begin{cases} Q & if x < \frac{l}{2} \\ 0 & if x > \frac{l}{2} \end{cases} \\ K_{ij} & = \begin{cases} k \left (\frac{\pi j}{l} \right )^2 \frac {l}{2} & if ~ i = j \\ k \left (\frac{\pi j}{l} \right )^2 \frac {l}{2 \pi} \left( \frac {sin(\pi (i-j))}{i-j} - \frac {sin(\pi (i+j))}{i+j} \right) & if ~ i <> j \end{cases} \\ f_i & = Q l \frac{1 - cos(\frac{i \pi}{2})}{i \pi} \, \end{cases} </math> <br> :[[Image:galerkinWRM.jpg]] :''To modify the parameters, edit this'' [[Matlab code Galerkin Method of the Resolution of the Poisson's equation with the WRM using global Shape Functions|'''Matlab code''']] == Least Squares Method &nbsp; &nbsp; <math>(W_i(x)=A(\varphi) \quad and \quad \overline{W_i}(x)=B(\varphi))</math> == {| ! General expression ! k and Q constant, using Fourier series |- | ::<math> \begin{cases} K_{ij} & = \int_{\Omega} 2 \left[ \frac{d^2 N_i(x)}{dx^2} \right] \left[ \frac{d}{dx} k \frac{d}{dx} \left( N_j (x) \right) \right] d\Omega \\ f_i & = \frac{}{} \int_{\Omega} 2 \left[ \frac{d^2 N_i(x)}{dx^2} \right] Q(x) d\Omega \end{cases} </math> | ::<math> \begin{cases} K_{ij} & = - \int_{\Omega} 2 \left( \frac{\pi i}{l} \right)^2 N_i(x) \left[ \frac{d}{dx} k \frac{d}{dx} \left( N_j (x) \right) \right] d\Omega \\ & = - 2 \left( \frac{\pi i}{l} \right)^2 \int_{\Omega} N_i(x) \left[ \frac{d}{dx} k \frac{d}{dx} \left( N_j (x) \right) \right] d\Omega \\ f_i & = - \int_{\Omega} 2 \left( \frac{\pi i}{l} \right)^2 \frac{}{} N_i(x) Q(x) d\Omega \\ & = - 2 \left( \frac{\pi i}{l} \right)^2 \frac{}{} \int_{\Omega} N_i(x) Q(x) d\Omega \end{cases} </math> |} :Note that, '''in this case and without generality''', due to the selected basis functions, the obtained solution is the same that the obtained one using the Galerkin method scaled with a factor <math>2 \left( \frac{\pi i}{l} \right)^2</math>. <br> :[[Image:LeastSquaresWRM.jpg]] :''To modify the parameters, edit this'' [[Matlab code Least Squares Method of the Resolution of the Poisson's equation with the WRM using global Shape Functions|'''Matlab code''']] ---- :[[Resolution of the Poisson's equation with the WRM using global Shape Functions | ''Back to the Resolution of the Poisson's equation with the Weighted Residual Method using global Shape Functions'']] :[[Poisson's_Equation#Weighted_Residual_Method_.28WRM.29 | ''Back to the Weighted Residual Method applied to the Poisson's equation'']] [[Category:Theory]]

Return to Specific Examples of the Resolution of the Poisson's equation with the WRM using global Shape Functions.