Scribd
Upload
251 views2 pages

Assignment4 PDF

This document provides instructions for Assignment 4 of the course ME 5107: Numerical Methods in Thermal Engineering. It includes 7 hand calculation problems and 3 programming problems involving numerical methods like bisection, Regula-Falsi, Newton's method, and solving systems of nonlinear equations. Students are asked to determine roots of equations and solve engineering problems related to heat transfer using iterative numerical techniques. Solutions are due by October 1st or 4th 2019 depending on which class batch the student is in.

Uploaded by

Vijay Vijay
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
251 views2 pages

Assignment4 PDF

This document provides instructions for Assignment 4 of the course ME 5107: Numerical Methods in Thermal Engineering. It includes 7 hand calculation problems and 3 programming problems involving numerical methods like bisection, Regula-Falsi, Newton's method, and solving systems of nonlinear equations. Students are asked to determine roots of equations and solve engineering problems related to heat transfer using iterative numerical techniques. Solutions are due by October 1st or 4th 2019 depending on which class batch the student is in.

Uploaded by

Vijay Vijay
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 2

ME 5107: Numerical Methods in Thermal Engineering Jul–Nov 2019

Assignment 4

General Instructions
• Solutions due by Submission due by 1st October 2019 for Batch I (Dr. Tiwari’s class); 4th October 2019 for
Batch II (Dr. Pallab’s class).
• Hand calculations can be done as homework and should be submitted as a report.

Hand calculations
1. The equation lnx = x2 − 1 has exactly√two real roots, α1 = 0.45 and α2 = 1. Determine for which initial
approximation x0 , the iteration xn+1 = 1 + lnxn converges to α1 and α2 .

2. Find the interval in which the smallest positive root of the following equations lies,
(a) tanx + tanhx = 0. Take the initial guess as (0, π/2).
(b) x3 − x − 4 = 0. Take the initial guess as (1, 2).
Determine the roots correct to two decimal places using bisection method.

3. For the following equations,


(a) x4 − x − 10 = 0, take the initial guess as (1, 2).
(b) x − e−x = 0, take the initial guess as (0, 1).
Determine the smallest positive root correct to three decimal places using Regula-Falsi and Secant methods.

4. Consider a function f (x) defined in [0, 2π] as



−1
 x ∈ [−4, 0]
f (x) = x − 1 x ∈ [0, 2]

1 x ∈ [2, 4]

Sketch the function f vs x. What is the appropriate method to be used to solve for the root of f (x) = 0 and
why? Perform five iterations using that method and report the solution at the end of each iteration.
5. A good approximation for n! is given by the function
 
1/2 x+1/2 −x 1 1
f (x) = (2π) x e 1+ +
12x 288x2

Though n! is defined for integer values of n, for large n, the approximation is very close. For this f (x), if
f (x) = 1000, calculate the value of x accurate up to 4 decimal places.
6. Determine the order of convergence of the iterative method

x0 f (xk ) − xk f (x0 )
xk+1 =
f (xk ) − f (x0 )

for finding a simple root of the equation f (x) = 0.

7. Calculate all the solutions of the system,

x2 + y 2 = 1.12 (1)
xy = 0.23 (2)

correct to three decimal places. Start with initial guess of x0 = 1 and y0 = 0.23.
ME 5107: Numerical Methods in Thermal Engineering Jul–Nov 2019

Figure 1: A schematic of the furnace indicating several heat transfer processes.

8. The hot combustion gases of a furnace are separated from the ambient air and its surroundings, as shown in
figure 1, which are at 25o C, by a brick wall 0.15 m thick. The brick has a thermal conductivity of 1.2 W/m.K and
a surface emissivity of 0.8. Under steady-state conditions an inner surface temperature of 352o C is measured.
Free convection heat transfer to the air adjoining the surface is characterised by a convection coefficient of h = 20
W/m2 . K. Calculate the brick outer surface temperature T2 using Newoton’s method by performing 5 iterations
with an initial guess of 1000 K. An energy balance on the outer surface of the furnace yields the following
equation:
T1 − T2
k = h(T2 − T∞ ) + σ(T24 − Tsur
4
) (3)
L
where σ = 5.67 × 10−8 W/m2 K4 is the Stefan-Boltzmann constant.

Programming
1. Show that the equation
 
π(x + 1)
f (x) = cos + 0.148x − 0.9062 = 0 (4)
8

has one root at the interval (−1, 0) and one in (0, 1). Calculate the negative root correct to six decimals after
10 iterations.

2. The equation

1 1
2e−x = + (5)
x+2 x+1
has two roots greater than −1. Calculate these roots correct to five decimal places.
3. Solve the system of non-linear equations given below with two starting guesses - x = [0.1, 1.2, 2.5]T and
x = [1, 0, 1]T . Do the two solutions converge to the same root? If not, why?

x+y+z =3 (6)
2 2 2
x +y +z =5 (7)
ex + xy − xz = 1 (8)

You might also like