- +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 Assignment2015 © Tutors Globe. All rights reserved.

## Finite Difference Method-Nonlinear ODE

Finite Difference Method – Nonlinear ODE::Heat conduction with radiationIf we once more consider the heat in a metal bar of length L but this time consider the effect of radiation as well as conduction then the steady state equation has the form

u

_{xx}− d(u^{4}− u^{4}b) = −g(x),Where u

_{b}is temperature of the background d incorporates a coefficient of radiation as well as g(x) is the heat source.If we again restore the continuous problem by its discrete approximation then we get:

(u

_{i+1}− 2u_{i}+ u_{i−1})/ h^{2}− d(u^{4}_{i}− u^{4}_{b}) = −g_{i}= −g(x_{i}).This equation is nonlinear in the unknowns therefore we no longer have a system of linear equations to solve however a system of nonlinear equations. One way to resolve these equations would be by the multivariable Newton method. In its place we introduce another interactive method.

:Relaxation Method for Nonlinear Finite DifferencesWe are able to rewrite equation as:

u

_{i+1}− 2u_{i}+ u_{i−1}= h^{2}d(u^{4}_{i}− u^{4}_{b}) − h^{2}g_{i}.From this we are able to solve for uiin terms of the other quantities:

2u

_{i}= u_{i+1}+ u_{i−1}− h^{2}d(u^{4}_{i}− u^{4}_{b}) + h^{2}g_{i}.After that we add uito both sides of the equation to obtain:

3u

_{i}= u_{i+1}+ u_{i}+ u_{i−1}− h^{2}d(u^{4}_{i}− u^{4}_{b}) + h^{2}g_{i},And afterwards divide by 3 to get:

u

_{i}=1/3(u_{i+1}+ u_{i}+ u_{i−1}) −h^{2}/3(d(u^{4}_{i}− u^{4}_{b}) + g_{i})Now for the major idea. We will start with an initial guess for the value of uifor each i which we can represent as a vector u

^{0}. Subsequently we will use the above equation to get better estimates, u^{1}, u^{2}, . . . , and expect that they converge to the correct answer.If we let:

u

^{j}= (u^{j}_{0}, u^{j}_{1}, u^{j}_{2}, . . . , u^{j}_{n−1}, u^{j}_{n})Denote the jth approximation afterwards we can obtain that j + 1st estimate from the formula:

Notice that g

_{i}and u_{b}don’t change. In the resulting equation we have u_{i}at every successive step depending on its previous value as well as the equation itself.:Implementing the Relaxation MethodIn the following program we resolve the finite difference equations with the boundary conditions:

u(0) = 0 and u(L) = 0. We let L = 4, n = 4, d = 1, and g(x) = sin(πx/4).

Notice that the vector u forever contains the current estimate of the values of u.

% mynonlinheat (lacks comments)L = 4; n = 4;

h = L/n; hh = h^2/3;

u0 = 0; uL = 0; ub = .5; ub4 = ub^4;

x = 0:h:L; g = sin(pi*x/4); u = zeros(1,n+1);

steps = 4;

u(1)=u0; u(n+1)=uL;

for j = 1:steps

u(2:n) = (u(3:n+1)+u(2:n)+u(1:n-1))/3 + hh*(-u(2:n).^4+ub4+g(2:n));

end

plot(x,u)

If you run this program with the specified n and steps the result will not seem reasonable.

We are able to plot the initial guess by adding the command plot(x,u) right before the for loop. We can as well plot successive iterations by moving the last plot(x,u) before the end. Now we can experiment as well as see if the iteration is converging. Try a variety of values of steps and n to produce a good plot. You will notice so as to this method converges quite slowly. In particular as we raise n we need to increase steps like n2 that is if n is large then steps needs to be really large.

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.