FDM Heat Equation FTCS Scheme with NBC

  • #1
Adyssa
203
3

Homework Statement



I have a quiz question that I'm struggling with. We've been working on using a FTCS scheme with two Essential Boundary Conditions, and now I have a problem with one EBC (ie static) and a Natural Boundary Condition (ie a derivative). The condensed problem statement:

[itex]u_{t} = 2u_{xx}, 0 < x < 1, t > 0[/itex]
[itex]u_{x}(0,t) = 0.5[/itex] <- NBC
[itex]u(1,t) = 3[/itex] <- EBC

and grid spacing [itex]Δx = 0.2, Δt = 0.01[/itex]

Construct a Finite Difference Approximation to the PDE using FTCS scheme. Approximate the derivative boundary condition using a central difference approximation at x = 0.

Homework Equations



As above. [itex][/itex]

The Attempt at a Solution



The basic scheme I have been using is this:

[itex]U^{n+1}_{i} = sU^{n}_{i-1} + (1-2s)U^{n}_{i}+sU^{n}_{i+1}[/itex], where [itex]s = \kappa \frac{Δt}{(Δx)^{2}} = 2 \frac{0.01}{0.04} = 0.5[/itex]

With regards to the derivative boundary condition, I believe that because [itex]u_{x}(0,t) = 0.5[/itex] is a positive gradient, that there is heating (not cooling) at the boundary. In my notes I have the following scheme for a Neumann Condition (heat flux at boundary):

Second order central dierence approximation to the boundary condition. Need grid point
outside the boundary to satisfy the boundary condition.


[itex]k \frac{U^{n}_{j+1} - U^{n}_{j-1}}{2Δx} = Fl[/itex] where [itex]k[/itex] is the conductivity and [itex]Fl[/itex] is the heat flux.

If this is even the correct expression, I'm not sure how to integrate it into my FTCS scheme in order to take into account the derivative boundary condition. In class we did some algebraic rearrangement of the equation to get [itex]U^{n}_{j-1}[/itex] by itself (LHS) and then substituted the RHS into the scheme, which seems appropriate:

[itex]U^{n+1}_{i} = \frac{2sΔx.Fl}{k} + (1-2s)U^{n}_{i}+2sU^{n}_{i+1}[/itex]

but we never worked the solution. I'm not sure which values to use for [itex]k[/itex] and [itex]Fl[/itex]. I think one of them will be 0.5? Also, it appears to me that I have now defined a constant heat flux at each point [itex]Δx[/itex] as I compute along the spatial dimension for [itex]U^{n}_{i}, (i=0,1,2,3,4,5) (n = time step)[/itex], but this is just speculation, I'm not a physics/eng major and this is all very new to me.
 
Last edited:
Physics news on Phys.org
  • #2
You basically need another grid line at x = -h. Then you have the usual equation relating u at -h, 0, and +h, plus you approximate the derivative with [u(+h, t) - u(-h, t)]/2h = a, where a is whatever constant the x-derivative must be at x = 0. So you get 2xN additional equations.
 
Back
Top