- +1-530-264-8006
- info@tutorsglobe.com

18,76,764

Questions

Asked

21,311

Experts

9,67,568

Questions

Answered

Start Excelling in your courses, Ask an Expert and get answers for your homework and assignments!!

Submit Assignment
## Insulated Boundary Conditions and Differential Equations

Insulated Boundary Conditions::InsulationIn lots of the previous sections we have considered fixed boundary conditions, that is u(0) = a,u(L) = b. We executed these simply by assigning u

^{j}_{0}= a and u^{j}_{n}= b for all j.We as well considered variable boundary conditions such as u(0, t) = g

_{1}(t). For instance we mighthave u(0, t) = sin(t) which could stand for periodic heating and cooling of the end at x = 0.A third significant type of boundary condition is called the insulated boundary condition. It is thereforenamed for the reason that it mimics an insulator at the boundary. Physically the consequence of insulation is thatno heat flows across the boundary. This signifies that the temperature gradient is zero which impliesthat must require the mathematical boundary condition u′(L) = 0.

To utilize it in a program we should replace u′(L) = 0 by a discrete version. Evoke that in our discreteequations we typically have L = x

_{n}. Evoke from the section on numerical derivatives that there arethree different manners to replace a derivative by a difference equation left, right as well as central differences the three of them at x_{n}would be:u′(x

_{n}) ≈ (u_{n}− u_{n−1})/ h≈ (u_{n+1}− u_{n}) h ≈ (u_{n+1}− u_{n−1}) 2hIf x

_{n}is the last node of our grid, then it is clear that we cannot use the right or central difference, but are stuck with the first of these. Setting that expression to zero entails:u

_{n}= u_{n−1}.This restriction is able to be easily execute in a program simply by putting a statement u(n+1)=u(n) inside the loop that updates values of the profile. Nevertheless since this method replaces u′(L) = 0 by an expression that is only accurate to first order it is not extremely accurate and is usually avoided.

Instead we want to utilize the most accurate version the central difference. For that we must have:

u′(L) = u′(x

_{n}) = (u_{n+1}− u_{n−1})/ 2h = 0.or simply

u

_{n+1}= u_{n−1}.Nevertheless u

_{n+1}would represent u(x_{n+1}) and x_{n+1}would be L + h, which is outside the domain. This nevertheless isn’t an obstacle in a program. We can merely extend the grid to one more node x_{n+1}, and let u_{n+1}always equal u_{n−1}by copying u_{n−1 }into u_{n+1}whenever u_{n−1}changes. The point x_{n+1}is fictional however a computer doesn’t know the difference between fiction and reality!This idea is carried out in the calculation.

Example- of an insulated boundary condition using a fictional point x

_{n+1}with u_{n+1}= u_{n−1}.A manner to think of an insulated boundary that makes sense of the point L+h is to think of two bars joined end to end where you let the second bar be mirror image of the first bar. If you do this after that no heat will flow across the joint which is exactly the same effect as insulating.

Another practical manner to implement an insulated boundary is to let the grid points straddle the boundary. For illustration suppose we want to impose insulated boundary at the left end of a bar that is u′(0) = 0 then you could let the first two grid points be at x

_{0}= −h/2 and x_{1}= h/2. Then you can let,u

_{0}= u_{1}.This will once more force the central difference at x = 0 to be 0.

:Execution in a linear equation by eliminationConsider the BVP:

u

_{xx}= −1, u(0) = 5, u′(1) = 0. (38.1)This symbolizes the steady state temperature of a bar with a uniformly applied heat source, with one end held at a fixed temperature and the other end insulated.

If we utilize 4 equally spaced intervals then

m = 4, L = 1, ⇒ h = L/m=1/4

and

x

_{0}= 0, x_{1}= .25, x_{2 }= .5, x_{3}= .75, x_{4}= 1, and x_{5}= 1.25.The point x

_{5}= 1.25 is outside the region and therefore fictional and the boundary condition at x_{0}= 0 is implemented asu

_{0}= 5.For the insulated condition we will necessitate:

u

_{5}= u_{3}.This creates the central difference for u′(x

_{4}) be zero. We are able to write the differential equation as a difference equation:(u

_{i−1}− 2u_{i}+ u_{i+1})/ h^{2}= −1or

u

_{i−1}− 2u_{i}+ u_{i+1 }= −0.0625, i = 1, 2, 3, 4.For i = 1, recalling that u0 = 5 we have:

5 − 2u

_{1}+ u_{2}= −.0625 or − 2u_{1 }+ u_{2 }= −5.0625.For i = 2 and i = 3 we have:

u

_{1}− 2u_{2}+ u_{3}= −.0625 and u_{2}− 2u_{3}+ u_{4}= −.0625.For i = 4 we have:

u

_{3}− 2u_{4}+ u_{5 }= −.5.Note that we now have 5 unknowns in our problem: u

_{1}, . . . , u_{5}Nevertheless from the boundary condition u_{5}= u_{3}and thus we can eliminate u5 from our equations and write2u

_{3}− 2u_{4}= −.0625.Summarizing we are able to put the unknown quantities in a vector u = (u

_{1}, u_{2}, u_{3}, u_{4})′and write the equations as a matrix equation Au = b where:and b = (−5.0625,−.0625,−.0625. − .0625)′. Resolve this system and plot the results

> u = A\b

> u = [5 ; u]

> x = 0:.25:1

>plot(x,u,’d’)

the interpolate with a spline.

The precise solution of this BVP is:

U(x) = 5 + x − .5x2.

Use holds on and plots this function on the same graph to compare. >xx = 0:.01:1; > uu = 5 + xx - .5*xx.^2; must see that our approximate solution almost perfect!

:Insulated boundary conditions in time-dependent problemsTo execute the insulated boundary condition in an explicit difference equation with time we need to copy values from inside the region to fictional points just outside the region.

An example

The stable state temperature u(r) (given in polar coordinates) of a disk subjected to a radially symmetric heat load g(r) and cooled by conduction to the rim of the disk as well as radiation to its environment is determined by the boundary value problem

∂

^{2}u/∂r^{2}+(1/r)(∂u/∂r)= d(u^{4}− u^{4}b) − g(r) with u(R) = uRand u′(0) = 0.Here u

_{b}is the (fixed) background temperature as well as uRis the (fixed) temperature at the rim of the disk.It is an example of information flow for the explicit method with an insulated boundary condition. An open circles are the fictional points as well as the curved arrows represent copying the value.

The class website has a program myheatdisk.m that executes these equations for parameter values R = 5, d = .1, u

_{R}= u_{b}= 10 and g(r) = (r − 5)^{2}. Notice that the equations comprise a singularity (discontinuity) at r = 0. How does the program evade this problem? How does the program execute u_{R}= 10 and u′(0) = 0? Run the program.Latest technology based Matlab Programming Online Tutoring AssistanceTutors, at the

www.tutorsglobe.com, take pledge to provide full satisfaction and assurance inMatlab Programming helpviaonline tutoring. Students are getting 100% satisfaction byonline tutorsacross the globe. Here you can get homework help for Matlab Programming, project ideas and tutorials. We provide email basedMatlab Programming help. You can join us to ask queries 24x7 with live, experienced and qualified online tutors specialized in Matlab Programming. ThroughOnline Tutoring, you would be able to complete your homework or assignments at your home. Tutors at theTutorsGlobeare committed to provide the best qualityonline tutoringassistance forMatlab Programming Homework helpandassignment helpservices. They use their experience, as they have solved thousands of the Matlab Programming assignments, which may help you to solve your complex issues of Matlab Programming.TutorsGlobeassure for the best quality compliance to your homework. Compromise with quality is not in our dictionary. If we feel that we are not able to provide thehomework helpas per the deadline or given instruction by the student, we refund the money of the student without any delay.