Nonlinear finite elements/Homework 2

Problem 1: Classification of PDEs

Two commonly encountered PDEs in engineering are

$u_{,xx}=\alpha u_{,t}$

and

$u_{,xxxx}=\alpha u_{,tt}$

1. What physical processes do the two PDEs model?
2. Write out the PDEs in elementary partial differential notation (that is, using symbols of the form $\frac{\partial u}{\partial x}$).
3. Determine whether the PDEs are elliptic, hyperbolic, or parabolic. Show how you arrived at your classification.

Problem 2: Models of Physical Problems

The differential equation which can be written as

$\frac{d}{dx}(a(x)\frac{du(x)}{dx})+c(x)u(x)=f(x)\text{for}x\in (O,L).$

arises in a number of fields of physics and engineering.

Write out the ODEs for at least five of these physical situations and identify the meanings of the variables.

Problem 3: Weighted Residual Methods

Give examples from the scientific literature, with full references, where the collocation, subdomain, and least squares weighted residual methods have been used.

Explain why the standard finite element method has not been used in these examples.

One example for each method will suffice.

Problem 4: Least Squares Method

The residual over an element (for the general second-order ODE) is defined as

$R^e(x,c_1^e,c_2^e,\dots ,c_n^e):=\frac{d}{dx}\left(a\frac{d{u}_h^e}{dx}\right)+c{u}_h^e-f(x)$

where

$u_h^e(x)=\sum _{j=1}^n{c}_j^e{\phi }_j^e(x).$

In the least squares method, the residual $R^e$ is minimized in the following sense:

$\frac{\partial }{\partial c_i}\underset{x_a}{\overset{x_b}{\int }}[R^e(x,c_1^e,c_2^e,\dots ,c_n^e){]}^{2}dx=0,i=1,2,\dots ,n.$

1. What is the weighting function used in the least squares approach?
2. Develop a least-squares finite element model for the problem. You don't have to find the solution to a specific problem - just the components of $𝐊$, $𝐮$, and $𝐟$ for an element will do.
3. Discuss some functions (${\phi }_j^e$) that can be used to approximate $u$.

Problem 5: Heat Transfer

The differential equation that is used to model one-dimensional heat transfer is

$-\frac{d}{dx}\left(k\frac{dT}{dx}\right)=q\text{for}x\in (0,L)$

where $T$ is the temperature, $k$ is the thermal conductivity, and $q$ is the heat generation.

Let us assume that the boundary conditions are

$\begin{array}{cc}T(0)& =T_0\\ {[k\frac{dT}{dx}+\beta (T-T_{\mathrm{\infty }})+\hat{q}]|}_{x=L}& =0\end{array}$

where $T_0$ is a prescribed temperature, $\beta$ is a convective heat transfer coefficient, $T_{\mathrm{\infty }}$ is the ambient temperature, and $\hat{q}$ is a prescribed heat flux.

1. Formulate the finite element model for this problem over one element, that is, find expressions for the stiffness matrix and the load vector.
2. Use ANSYS (or any other software tool) to solve the problem of heat conduction in an insulated rod. Compare your solution with the exact solution (if available). Try at least three refinements of the mesh to confirm that your solution converges.
3. Write your own code to solve the problem in part 2 (use two elements). Compare your solution with the ANSYS (or other software) solution.

Use the following data:

$\begin{array}{cc}q& =0\\ \hat{q}& =0\\ L& =0.1\text{m}\\ k& =0.01\text{W/(m-K)}\\ \beta & =25{\text{W/(m}}^2\text{-K)}\\ T_0& =323\text{K}\\ T_{\mathrm{\infty }}& =278\text{K}\end{array}$

You can look at the 1-D heat conduction problem on the Introduction to Ansys page to see how ANSYS deals with heat conduction.