Homework 4
Chris Nkinthorn
Finite Element Methods
Problem 1
Prompt: In section 2.3 of the text (pg. 60-64) the strong form and weak form of the equations for linear heat conduction are given. Using a Galerkin based method of weighted residuals, show the construction of the given weak form starting from the strong form, stated as:
Given $f : \Omega \rightarrow \mathbb{R}, q : \Gamma{o} \rightarrow \mathbb{R} \text { and } h : \Gamma{h} \rightarrow \mathbb{R}, \text { find } u : \overline{\Omega} \rightarrow \mathbb{R}$ such that
With the functions from candidate spaces: $u \in \delta$ and $w \in \mathcal{V}$.
Beginning with the Heat Equation $[\ref{eqnHeat}]$, we cannot solve in the strong form presented. Instead, we minimize a weighting function $w$ in the domain $\Omega$, where the closure of $\Omega$ is $\Gamma$. Distributing the weight function to produce two integral terms. Solving the first term using integration by parts (IBP) produces
Algebraic rearrangement with knowns on the right and the unknown on the left:
Problem 2
$\textcolor{red}{\texttt{To Be Graded }}$ Exercise 1 on page 68 of the text book. This exercise is a multidimensional analog of the the one contained in Sec. 1.8. Let
One side of $\Gamma{int}$ is arbitrarily designated to be the “+ side” and the other is the “- side”. Let $n^+$ and $n^-$ be unit normals to $\Gamma{int}$, with the relationship:
Consider the weak formulation$[\ref{eqn2BC}]$ and assume that $w$ and $u$ are smooth on the element interiors but may experience discontinuities in gradient across element boundaries. Restating $[\ref{eqn2BC}]$, the equation is made homogenous by moving the left hand side (LHS) term to the right hand side (RHS) of the equation.
where IBP was
so substitution into
This statement is local, so on the global mesh is the sum of all elements from $1$ to $n_{el}$:
The normal flux relationship of Eqn $[\ref{eqnJump}]$ is applied into the middle term. Decomposing the $\Gamma = \cancelto{0}{\Gamma_g} + \Gamma_h$, produces
Problem 3
Exercise 1 on page 71 of the text book
Let
where
${\bf d}^e$ is the element temperature vector. Show that the heat flux vector at point $x \in \Omega^e$ can be calculated from the formulation
Not sure how this expression was constructed, to move the $B$ matrix as a summation of element nodes. It looks like collocation because of evaluation at specific locations but I’m not sure how to prove it.
Problem 4
Exercise 2 on page 71 of the textbook, Consider the strong statement where the replacement on $\Gamma_h$ is
where $\lambda\geq0$ is a function of $x \in\Gamma_h$. To generalize, continue from $[\ref{eqnLambda}]$, where the additional term creates:
Assuming that the original weak form expression was also positive definite, then the additional contribution to $k_{ab}^e$ is based purely on $\lambda$. To prove that $\bf K$ is positive definite, then we would need to show that $\bf c^T\cdot K \cdot c \geq 0$ and $0$ only when $\bf c$ is the 0 vector, $\bf 0 $.
Last updated