See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/369882349

CMP Lab Assignment week 7 Numerical integration

Presentation · April 2023

DOI: 10.13140/RG.2.2.20995.17445

CITATIONS READS

0 801

1 author:

Abhishek D. Patange
ABB
111 PUBLICATIONS 679 CITATIONS

SEE PROFILE

All content following this page was uploaded by Abhishek D. Patange on 08 April 2023.

The user has requested enhancement of the downloaded file.

Computational Methods &
Programming
[ME-21009]

Assignment for the 7th week

Submission deadline: 17th April 2023

Abhishek D. Patange, Ph.D.

College of Engineering Pune
Numerical integration assignment
No. Question Batch
Evaluate the following integral:
(a) by analytically;
(b) by single application trapezoidal rule (n=1)
(c) by multiple-application trapezoidal rule (i.e. n=2, n=4, n=6, n=8)
(d) by single application Simpson’s 1/3 rule (n=2)
(e) by multiple-application Simpson’s 1/3 rule (i.e. n=4, n=6, n=8, n=10)
1
(f ) by single application Simpson’s 3/8 rule (n=3)
(g) by multiple-application Simpson’s rule (i.e. n=6, n=9, n=12, n=15)
For each of the numerical estimates (b) through (g), determine the percent
relative error based on (a).
Estimate the ‘n’ required to get answer of trapezoidal, Simpson’s 1/3 rule,
Simpson’s 3/8 rule as same as analytical approach.

2 Suppose that the upward force of air resistance on a falling object is A

proportional to the square of the velocity. For this case, the velocity can be
computed as

where cd = a second-order drag coefficient. (a) If g = 9.8 m/s2, m = 68.1 kg

and cd = 0.25 kg/m, use analytical integration to determine how far the
object falls in 10 s. (b) Make the same evaluation, but evaluate the integral
with the multiple-segment Simpson’s 3/8 rule. Use a sufficiently high ‘n’ that
you get three significant digits of accuracy.

Suppose that the upward force of air resistance on a falling object is

proportional to the square of the velocity. For this case, the velocity can be
computed as

where cd = a second-order drag coefficient. (a) If g = 9.8 m/s2, m = 68.1 kg

and cd = 0.25 kg/m, use analytical integration to determine how far the
object falls in 10 s. (b) Make the same evaluation, but evaluate the integral
with the multiple-segment Simpson’s 1/3 rule. Use a sufficiently high ‘n’ that
you get three significant digits of accuracy.

Suppose that the upward force of air resistance on a falling object is

proportional to the square of the velocity. For this case, the velocity can be
computed as

where cd = a second-order drag coefficient. (a) If g = 9.8 m/s2, m = 68.1 kg

and cd = 0.25 kg/m, use analytical integration to determine how far the
object falls in 10 s. (b) Make the same evaluation, but evaluate the integral
with the multiple-segment trapezoidal rule. Use a sufficiently high ‘n’ that you
get three significant digits of accuracy.

Suppose that the upward force of air resistance on a falling object is

proportional to the square of the velocity. For this case, the velocity can be
computed as
D
where cd = a second-order drag coefficient. (a) If g = 9.8 m/s2, m = 68.1 kg
and cd = 0.25 kg/m, use analytical integration to determine how far the
object falls in 10 s. (b) Make the same evaluation, but evaluate the integral
with the multiple-segment Simpson’s 3/8 rule. Use a sufficiently high ‘n’ that
you get three significant digits of accuracy.

Suppose that the upward force of air resistance on a falling object is

proportional to the square of the velocity. For this case, the velocity can be
computed as

where cd = a second-order drag coefficient. (a) If g = 9.8 m/s2, m = 68.1 kg

and cd = 0.25 kg/m, use analytical integration to determine how far the
object falls in 10 s. (b) Make the same evaluation, but evaluate the integral
with the multiple-segment Simpson’s 1/3 rule. Use a sufficiently high ‘n’ that
you get three significant digits of accuracy.

Suppose that the upward force of air resistance on a falling object is

proportional to the square of the velocity. For this case, the velocity can be
computed as

where cd = a second-order drag coefficient. (a) If g = 9.8 m/s2, m = 68.1 kg

and cd = 0.25 kg/m, use analytical integration to determine how far the
object falls in 10 s. (b) Make the same evaluation, but evaluate the integral
with the multiple-segment trapezoidal rule. Use a sufficiently high ‘n’ that you
get three significant digits of accuracy.

Suppose that the upward force of air resistance on a falling object is

proportional to the square of the velocity. For this case, the velocity can be
computed as
G

where cd = a second-order drag coefficient. (a) If g = 9.8 m/s2, m = 68.1 kg

and cd = 0.25 kg/m, use analytical integration to determine how far the
object falls in 10 s. (b) Make the same evaluation, but evaluate the integral
with the multiple-segment Simpson’s 3/8 rule. Use a sufficiently high ‘n’ that
you get three significant digits of accuracy.

Suppose that the upward force of air resistance on a falling object is

proportional to the square of the velocity. For this case, the velocity can be
computed as

where cd = a second-order drag coefficient. (a) If g = 9.8 m/s2, m = 68.1 kg

and cd = 0.25 kg/m, use analytical integration to determine how far the
object falls in 10 s. (b) Make the same evaluation, but evaluate the integral
with the multiple-segment Simpson’s 1/3 rule. Use a sufficiently high ‘n’ that
you get three significant digits of accuracy.

Evaluate the integral of the following tabular data with (a) the trapezoidal
rule and (b) Both Simpson’s rules
3 A

Evaluate the integral of the following tabular data with (a) the trapezoidal
rule and (b) Both Simpson’s rules
4 B

Determine the mean value of the function between x = 2 and 10

5 C
using all 3 methods of numerical integration.

The following data was collected for a cross-section of a river (y = distance

6 D
from bank, H = depth and U = velocity):
Use numerical integration to compute the (a) average depth, (b) cross-
sectional area, (c) average velocity, and (d) the flow rate. Note that the
cross-sectional area (Ac) and the flow rate (Q) can be computed as

The outflow concentration from a reactor is measured at a number of times

over a 24-hr period:

The flow rate for the outflow in m3/s can be computed with the
following equation:

7 E

Use the best numerical integration method to determine the flow weighted
average concentration leaving the reactor over the 24-hr period,

An 11-m beam is subjected to a load, and the shear force follows the
equation

where V is the shear force and x is length in distance along the beam. We
know that V = dM/dx, and M is the bending moment. Integration yields the
8 relationship F

If Mo is zero and x = 11, calculate M using (a) analytical integration, (b)

multiple-application trapezoidal rule, and (c) multiple application
Simpson’s rules. For (b) and (c) use 1-m increments.
The work produced by a constant temperature, pressure volume
thermodynamic process can be computed as

9 G
where W is work, p is pressure, and V is volume. Using a combination of the
trapezoidal rule, Simpson’s 1/3 rule, and Simpson’s 3/8 rule, use the
following data to compute the work in kJ (kJ=kN ) m):

Determine the distance traveled for the following data:

10 (a) Use the trapezoidal rule, (b) the best combination of the trapezoidal and H
Simpson’s rules, and (c) analytically integrating second and third-order
polynomials determined by regression.

The total mass of a variable density rod is given by

where m = mass, ρ (x) = density, Ac(x) = cross-sectional area, x = distance

along the rod and L = the total length of the rod. The following data has been All
11 measured for a 10-m length rod. Determine the mass in kilograms to the batches
best possible accuracy.

A transportation engineering study requires that you determine the number

of cars that pass through an intersection traveling during morning rush hour.
You stand at the side of the road and count the number of cars that pass
every 4 minutes at several times as tabulated below. Use the best numerical All
12
method to determine batches
(a) the total number of cars that pass between 7:30 and 9:15, and
(b) the rate of cars going through the intersection per minute. (Hint: Be
careful with units.)
View publication stats