1
$\begingroup$

Given

$$(\textrm{Pe}-1)u_{i+1} + 2u_i - (\textrm{Pe}-1)u_{i-1} = \frac{h^2}{\mu}$$

where $\textrm{Pe}$ stands for Péclet number, $h$ is the step size of the mesh and $\mu$ is the diffusion (constant)

How to solve such equation? I'm trying to solve it to understand why there are no oscillations in a finite element method approximation with the upwind schema.

UPDATE 1

The roots I got are: $r_1 = 1$ and $r_2 = \frac{1+\textrm{Pe}}{1-\textrm{Pe}}$

So I have the following guess for the general solution

$$u_i = C_1r_1^i + C_2r_2^i$$

To find the constants $C_1$ and $C_2$ I need to initial conditions. What are those? Are they related to the boundary conditions? I have $u(0) = 1$ and $u'(1) = 1$. How to proceed?

UPDATE 2

According to this text, I replace the roots in my guess for the solution to the homogeneous part.

$$u_i = C_1(1)^i + C_2\Big(\frac{1+\textrm{Pe}}{1-\textrm{Pe}}\Big)^i$$

Applying the initial conditions I have

$$u_0 = 1 = C_1(1)^0 + C_2\Big(\frac{1+\textrm{Pe}}{1-\textrm{Pe}}\Big)^0$$ $$1-C_1 = C2$$

for the other initial condition, what should I do, I mean I don't know how to proceed with that Neumann condition there.

UPDATE 3

Before applying the initial conditions to find the values of $C_1$ and $C_2$ I tried to solve for the particular case following the example on the information of this site

So I came up with using a guess like $f_i = ki$. Replacing on the original equation I have

$$(\textrm{Pe}-1)k(i+1) + 2ki - (\textrm{Pe}-1)k(i-1) = \frac{h^2}{\mu}$$

which after rearranging gives me

$$\big[(\textrm{Pe}-1)+2-(\textrm{Pe}-1)\big]ki + \big[(\textrm{Pe}-1)+(\textrm{Pe}+1)\big]k = \frac{h^2}{\mu}$$

finally $$k = \frac{h^2}{\mu \big[(\textrm{Pe}-1)+(\textrm{Pe}+1)\big]}$$

Hence, the complete (I think) equation is:

$$u_i = C_1(1)^i + C_2\Big(\frac{1+\textrm{Pe}}{1-\textrm{Pe}}\Big)^i+\frac{h^2 i}{\mu \big[(\textrm{Pe}-1)+(\textrm{Pe}+1)\big]}$$

Am I right? I just need to apply initial conditions. The $u(0) = 1$ is easy but the other one $u'(1) = 1$ how to apply that?

  • 0
    Solve for what? In Which Context? Are $A,B,C$ matrices? A little bit more explanation.2012-12-19
  • 0
    Your solution should be of the form $u(i)=u_h(i)+u_p(i)$, where $u_h$ and $u_p$ are the corresponding solutions of the homogeneous and particular solutions. $Here is a [related problem](http://math.stackexchange.com/questions/209640/solving-a-recurrence-relation-with-the-characteristic-equation/209652#209652).2012-12-19
  • 0
    This being a finite difference representation of a diffusion eq'n, I would set $u_0 = 1$ and $u_1 - u_0 = h$, or $u_1 = 1+h$. The initial conditions are different from boundary conditions: initial conditions deal with the behavior at $t=0$ while boundary conditions deal with the solution on a manifold. For a boundary-value problem, you would have to use a very different numerical representation of the equation.2012-12-19
  • 0
    @rlgordonma but this problem I'm trying to solve is the steady case. So there are no initial conditions. It's an elliptic 1D steady problem.2012-12-19
  • 0
    OK then, do this. You have mixed bounday conditions: Dirichlet on the left, Neumann on the right. For Dirichlet conditions, you have $u_0=1$. For Neumann, if there are N steps, you have $u_{N}-u_{N-1} = h$. You then need to apply the discrete equation successively into the interior of the region. To do this, you may be able to use the solution you derived in the initial value case to get the values of $u_{1}$, etc. Then take those new "boundary" values and propagate into the interior. Another way is to simply use a relaxation method (average until you reach equilibrium).2012-12-19
  • 0
    @rlgordonma could you put some formulas or explain in more detail I didn't get your point.2012-12-19
  • 0
    Oh! and please someone tell me (if possible) if this the way I should be analyzing to know why there are no oscillations on the numerical approximation of this problem.2012-12-19

1 Answers 1