Scribd
Upload
12 views8 pages

Numerical Case Study (Project Fall 2024-2025)

Uploaded by

samer1234569g
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
12 views8 pages

Numerical Case Study (Project Fall 2024-2025)

Uploaded by

samer1234569g
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/ 8

The Hashemite University - Faculty of Engineering

Department of Allied Engineering sciences


Numerical Analysis Course # 110402303
Coordinator: Eng. Yousef Almashakbeh - Fall, 2024/2025

"List of Projects"
Note: See descriptions of the projects below this table.

Due Date: 14/12/2025.

Instructor will arrange and determine groups and the name of the project for each group.

1 Uniform beam under distributed load.


2 Temperature Distribution in Thin Rod
3 Newton divided difference polynomial and Numerical
Differential & Integration.
4 Spring Mass Problem
Case 1: Uniform beam under distributed load.
In the shown Figure, a uniform beam subject to a linearly increasing distributed load.
The deflection y (m) can be expressed by

Where E is the modulus of elasticity and I


is the moment of inertia (m4), L length of
beam.
Use the following parameters L=600 cm,
E=50,000 kN/cm2, I= 30.000 cm4, wo=2.5
kN/cm, to find the requirements

1) Develop MATLAB code to


determine the point of maximum
deflection by using numerical
method (bisection, false position
method,….). Hint(The value of x
where dy/dx=0).
2) Plot the point of maximum
deflection versus iteration number.
3) Plot the values of the relative
approximate error of the point of
maximum deflection (𝜖𝑎,𝑥) versus
iteration number.

4) Plot the following quantities versus distance along the beam


a) Displacement (y).
b) Slope 𝜃(𝑥) = 𝑑𝑦/𝑑𝑥.
c) Moment 𝑀(𝑥) = 𝐸𝐼𝑑2𝑦/𝑑𝑥2.
d) Shear 𝑉(𝑥) = 𝐸𝐼𝑑3𝑦/𝑑𝑥3.
e) Loading 𝑤(𝑥) = −𝐸𝐼𝑑4𝑦/𝑑𝑥4.
Case 2:Temperature Distrbution In Thin Rod
The differential equation derives from a heat balance for along thin rod as follow:

Where T=temperature (oC), x=distance along the rod (m), h'=a heat transfer coefficent between the
rod and the ambient air (m-2), and Ta=the temperature of the surrounding air (oC).This equation
can be transformed into a set of linear algebraic equations by using a finite dividd difference
approximation for the second derivative as follow:

Where Ti designates the temperature at node i. This approximation can be substituded into eq(1)
to give:

−𝑇𝑖−1 + (2 + ℎ′Δ𝑥2)𝑇𝑖 − 𝑇𝑖+1 = ℎ′Δ𝑥2𝑇𝑎 ................................................... 𝑒𝑞(3)

This equation can be written for each of the interior nodes of the rod resulting in a tridiagonal
system of equations. The first and the last nodes at the rods ends are fixed by boundaryconditions.
T(x=0)=200 oC, T(x=10)= 40 oC, h'=0.02 m-2, Ta=10 oC, Δx=2m.

If this system of linear equation can be solved using Gauss-Seidel and Jacobi iterations methods.
Please answer the following:
1. Develop a MATLAB program to solve the system above using both Gauss-Seidel and
Jacobi iteration methods after n of iterations.
2. Plot the values of each tempratures (T1,T2,T3,T4) versus iteration number in both
solutions. Compare the values you get from both solution methods.
3. Plot the values of the relative approximate error of each temperatures (𝜀𝑎) versus iteration
number in both solutions. Compare the values you get from both solution methods.
4. Use LU decomposition to find the exact solution {𝑇𝑖} and compare your exact values to
those you got from point 1.
5. Plot the temperatures you got from point 4 along the rod.
Case 3: Newton divided difference polynomial and Numerical Differantion &
Integration.

i 0 1 2 3 4 5 6
xi 1.4 1.6 1.8 2.0 2.2 2.4 2.6
yi 2.151 2.577 3.107 4.015 5.105 6.314 7.015

1. Develop a MATLAB code of Newton divided difference to estimate the interpolating


polynomial of the values given in the table above for x and y.
2. Develop a MATLAB code of numerical differention to estimate y'(1.8) and y''(1.8) as
accurate as possible.
2.6
3. Develop a MATLAB code to estimate∫1.4 𝑓(𝑥)𝑑𝑥 by using Simpson’s 1/3, Simpson’s
3/8 and Trapezodial Rule with end correction.
4. Plot the interpolating polynomial in part 1.
Case 4: Spring Mass Problem:
The mass m is attached to a spring of free length b and stiffness k. The coefficient of friction
between the mass and the horizontal rod is μ. The acceleration of the mass can be shown to be

If the mass is released from rest at x = b, its speed at x = 0 is given by

vo =

1. Compute v0 by numerical Simpson’s 1/3 and 3/8 integration and compare between them
with different step size, using the data m = 0.9 kg, b = 0.6 m, μ=0.3, k = 100 N/m, and g
= 9.81 m/s2,
2. Develop a MATLAB code to solve the equation for both methods.
3. Plot the acceleration of the mass versus x, and find the area under the curve by MATLAB
built-in function.
4. Can we find an exact solution???Try it.

Figure (1)
Numerical Analysis
Term Project Report Writing Format
Warning: Anything copied from the internet will not be accepted no matter how significant it may
be. All what you need is in textbooks. Read and rephrase.

The following items should be included in the report:

1. Title (Cover) Page (see template below):


In the first page of the report (not to be numbered), University, College, Department, Course Title and
Number, Project title, Student Name and Date have to be included. Student should use black color only.
2. Abstract ( maximum of two paragraphs):
In the abstract, students should introduce the problem, state the approach and scope of the work, and
summarize the findings/conclusions of the work.
3. Table of Contents:
In this table, the report contents with respective page numbers have to be enumerated.
4. List of Figures (if any):
In this section, all figures in the report are to be listed together with respective page numbers.
5. List of Tables (if any):
In this section, all tables in the report are to be listed together with respective page numbers.
6. Main Text Sections:
The Sections to be included in the report are as follows
− Section 1: Introduction
Includes Statement of the problem, objectives, and the numerical method used
− Section 2: Literature Review (if any): optional
Includes citation of relevant work and results.
− Section 3: Formulation of the problem in Numerical Sense (how the method is applied
and state the algorithm: solution steps)
− Section 4: Results
− Section 5: Conclusions

7. References: include at least one reference: Format as follows:

Books:
Abbas, A. (1999), “ Numerical Methods, " John Wiley and Sons, NY.

Journal Articles

Fulan, A.M.S. and Fulan, Y.N.M.(1999), “Drag reduction in pipe flow, “ Int. J. of
Multiphase Flow, vol. 56, no. 12, pp. 1234-1256.
8. Appendices (if any)
Appendix Number and Title (normally the computer programs are listed in
the Appendix)
Numerical Analysis
(# 303)
Term Project

Faculty of Engineering
The Hashemite University, Zarqa- Jordan

Project Title Goes here (in black)

(FirstName) (LastName) ( ID #) (Majors)

Section #:
Group #:

December, 2025

You might also like