Numerical Methods for Engineers
Numerical Methods for Engineers
7th Edition
ISBN: 9780073397924
Author: Steven C. Chapra Dr., Raymond P. Canale
Publisher: McGraw-Hill Education
bartleby

Concept explainers

Differential Equations
Have you ever observed the movement of the satellites and rockets or how the planets go around the sun in their orbits? How will you determine their motion? Definitely they follow some scientific laws and these kinds of problems are framed as differentia…
bartleby

Videos

Textbook Question
Book Icon
Chapter 24, Problem 44P

A jet fighter's position on an aircraft carrier's runway was timed during landing:

t, s 0 0.52 1.04 1.75 2.37 3.25 3.83
x, m 153 185 208 249 261 271 273

where x is the distance from the end of the carrier. Estimate (a) velocity ( d x / d t ) and (b) acceleration ( d v / d t ) using numerical differentiation.

(a)

Expert Solution
Check Mark
To determine

To calculate: The velocity of the jet fighter landing on an aircraft carrier’s runaway using Numerical Integration.

Answer to Problem 44P

Solution:

The velocities at different points of time of a Jet fighter are given below,

t,s 0 0.52 1.04 1.75 2.37 3.25 3.83
x,m 153 185 210 249 261 271 273
v,m/s 68.26923 54.80769 50.97398 35.93855 16.05180 6.59273 0.303817

Explanation of Solution

Given Information:

The position of a Jet fighter on an aircraft carrier’s runway is,

t,s 0 0.52 1.04 1.75 2.37 3.25 3.83
x,m 153 185 210 249 261 271 273

Formula Used:

Second-order Lagrange interpolating polynomial.

f′(t)=f(ti−1)2t−ti−ti+1(ti−1−ti)(ti−1−ti+1)+f(ti)2t−ti−1−ti+1(ti−ti−1)(ti−ti+1)+f(ti+1)2t−ti−1−ti(ti+1−ti−1)(ti+1−ti)…… (1)

Here, t is the value at which the derivative is to be estimated.

Calculation:

Consider the function f(x) to be position (x) as shown in table above. So, f′(x) will represent the acceleration (v).

Apply second-order Lagrange interpolating polynomial to each set of three adjacent points for calculation of velocities at different points.

Calculate f′(0).

Substitute t=0 and i=1 in equation(1).

f′(0)=f(0)2(0)−t1−t2(t0−t1)(t0−t2)+f(t1)2(0)−t0−t2(t1−t0)(t1−t2)+f(t2)2(0)−t0−t1(t2−t0)(t2−t1)

Substitute the value of function from above table.

f′(0)=153(2(0)−0.52−1.04(0−0.52)(0−1.04))+185(2(0)−0−1.04(0.52−0)(0.52−1.04))+210(2(0)−0−0.52(1.04−0)(1.04−0.52))=68.26923

Calculate f′(0.52).

Substitute t=0.52 and i=2 in equation(1).

f′(0.52)=f(0.52)2(0.52)−t2−t3(t1−t2)(t1−t3)+f(t2)2(0.52)−t1−t3(t2−t1)(t2−t3)+f(t3)2(0.52)−t1−t2(t3−t1)(t3−t2)

Substitute the value of function from above table.

f′(0.52)=153(2(0.52)−0.52−1.04(0−0.52)(0−1.04))+185(2(0.52)−0−1.04(0.52−0)(0.52−1.04))+210(2(0.52)−0−0.52(1.04−0)(1.04−0.52))=54.80769

Calculate f′(1.04).

Substitute t=1.04 and i=3 in equation(1).

f′(1.04)=f(1.04)2(1.04)−t3−t4(t2−t3)(t2−t4)+f(t3)2(1.04)−t2−t4(t3−t2)(t3−t4)+f(t4)2(1.04)−t2−t3(t4−t2)(t4−t3)

Substitute the value of function from above table.

f′(1.04)=185(2(1.04)−1.04−1.75(0.52−1.04)(0.52−1.75))+210(2(1.04)−0.52−1.75(1.04−0.52)(1.04−1.75))+249(2(1.04)−0.52−1.04(1.75−0.52)(1.75−1.04))=50.97398

Calculate f′(1.75).

Substitute t=1.75 and i=1 in equation(1).

f′(1.75)=f(1.75)2(1.75)−t4−t5(t3−t4)(t3−t5)+f(t4)2(1.75)−t3−t5(t4−t3)(t4−t5)+f(t5)2(1.75)−t4−t3(t5−t3)(t5−t4)

Substitute the value of function from above table.

f′(1.75)=210(2(1.75)−1.75−2.37(1.04−1.75)(1.04−2.37))+249(2(1.75)−1.04−2.37(1.75−1.04)(1.75−2.37))+261(2(1.75)−1.04−1.75(2.37−1.04)(2.37−1.75))=35.93855

Calculate f′(2.37).

Substitute t=2.37 and i=5 in equation(1).

f′(2.37)=f(2.37)2(2.37)−t5−t6(t4−t5)(t4−t6)+f(t5)2(2.37)−t4−t6(t5−t6)(t5−t4)+f(t6)2(2.37)−t4−t5(t6−t4)(t6−t5)

Substitute the value of function from above table.

f′(2.37)=249(2(2.37)−2.37−3.25(1.75−2.37)(1.75−3.25))+261(2(2.37)−1.75−3.25(2.37−1.75)(2.37−3.25))+271(2(2.37)−1.75−2.37(3.25−1.75)(3.25−2.37))=16.05180

Calculate f′(3.25).

Substitute t=3.25 and i=6 in equation(1).

f′(3.25)=f(3.25)2(3.25)−t6−t7(t5−t6)(t5−t7)+f(t6)2(3.25)−t6−t7(t6−t5)(t6−t7)+f(t7)2(3.25)−t5−t6(t7−t6)(t7−t5)

Substitute the value of function from above table.

f′(3.25)=261(2(3.25)−3.25−3.83(2.37−3.25)(2.37−3.83))+271(2(3.25)−2.37−3.83(3.25−2.37)(3.25−3.83))+273(2(3.25)−2.37−3.25(3.83−2.37)(3.83−3.25))=6.59273

Calculate f′(3.83).

Substitute t=3.83 and i=7 in equation(1).

f′(3.83)=f(3.83)2(3.83)−t7−t8(t6−t7)(t6−t8)+f(t7)2(3.83)−t6−t8(t7−t6)(t7−t8)+f(t8)2(3.83)−t6−t8(t8−t6)(t8−t7)

Substitute the value of function from above table.

f′(3.83)=261(2(3.83)−3.25−3.83(2.37−3.25)(2.37−3.83))+271(2(3.83)−2.37−3.83(3.25−2.37)(3.25−3.83))+273(2(3.83)−2.37−3.25(3.83−2.37)(3.83−3.25))=0.303817

The velocities at different points of time of a Jet fighter are given below,

t,s 0 0.52 1.04 1.75 2.37 3.25 3.83
x,m 153 185 210 249 261 271 273
v,m/s 68.26923 54.80769 50.97398 35.93855 16.05180 6.59273 0.303817

(b)

Expert Solution
Check Mark
To determine

To calculate: The acceleration of the jet fighter landing on an aircraft carrier’s runaway using Numerical Integration.

Answer to Problem 44P

Solution:

The complete table of acceleration for different values of velocity is given below,

t,s x,m v=dxdt a=dvdt
0 153 68.26923 -35.14510
0.52 185 54.80769 -16.63004
1.04 210 50.97398 -13.20841
1.75 249 35.93855 -26.99478
2.37 261 16.05180 -23.26046
3.25 271 6.59273 -10.80560
3.83 273 0.303817 -10.88029

Explanation of Solution

Given Information:

The velocities at different points of time of a Jet fighter are given below as calculated in Part (a).

t,s 0 0.52 1.04 1.75 2.37 3.25 3.83
x,m 153 185 210 249 261 271 273
v,m/s 68.26923 54.80769 50.97398 35.93855 16.05180 6.59273 0.303817

Formula Used:

Second-order Lagrange interpolating polynomial.

f′(t)=f(ti−1)2t−ti−ti+1(ti−1−ti)(ti−1−ti+1)+f(ti)2t−ti−1−ti+1(ti−ti−1)(ti−ti+1)+f(ti+1)2t−ti−1−ti(ti+1−ti−1)(ti+1−ti)…… (1)

Here, t is the value at which the derivative is to be estimated.

Calculation:

Consider the function f(x) to be velocity (v) as shown in table above. So, f′(x) will represent the acceleration (a).

Calculate the acceleration (a) at different points with respect to the velocities (v) given in the above table with the help of Equation (1).

Calculate f′(0).

Substitute t=0 and i=1 in equation(1).

f′(0)=f(0)2(0)−t1−t2(t0−t1)(t0−t2)+f(t1)2(0)−t0−t2(t1−t0)(t1−t2)+f(t2)2(0)−t0−t1(t2−t0)(t2−t1)

Substitute the value of function from above table.

f′(0)=68.26923(2(0)−0.52−1.04(0−0.52)(0−1.04))+54.80769(2(0)−0−1.04(0.52−0)(0.52−1.04))+50.97398(2(0)−0−0.52(1.04−0)(1.04−0.52))=−35.14510

Calculate f′(0.52).

Substitute t=0.52 and i=2 in equation(1).

f′(0.52)=f(0.52)2(0.52)−t2−t3(t1−t2)(t1−t3)+f(t2)2(0.52)−t1−t3(t2−t1)(t2−t3)+f(t2)2(0.52)−t1−t2(t3−t1)(t3−t2)

Substitute the value of function from above table.

f′(0.52)=68.26923(2(0.52)−0.52−1.04(0−0.52)(0−1.04))+54.80769(2(0.52)−0−1.04(0.52−0)(0.52−1.04))+50.97398(2(0.52)−0−0.52(1.04−0)(1.04−0.52))=−16.63004

Calculate f′(1.04).

Substitute t=1.04 and i=3 in equation(1).

f′(1.04)=f(1.04)2(1.04)−t3−t4(t2−t3)(t2−t4)+f(t1)2(1.04)−t2−t4(t3−t2)(t3−t4)+f(t4)2(1.04)−t2−t3(t4−t2)(t4−t3)

Substitute the value of function from above table.

f′(1.04)=54.80769(2(1.04)−1.04−1.75(0.52−1.04)(0.52−1.75))+50.97398(2(1.04)−0.52−1.75(1.04−0.52)(1.04−1.75))+35.93855(2(1.04)−0.52−1.04(1.75−0.52)(1.75−1.04))=−13.20841

Calculate f′(1.75).

Substitute t=1.75 and i=4 in equation(1).

f′(1.75)=f(1.75)2(1.75)−t4−t5(t3−t4)(t3−t5)+f(t1)2(1.75)−t3−t5(t4−t3)(t4−t5)+f(t2)2(1.75)−t4−t3(t5−t3)(t5−t4)

Substitute the value of function from above table.

f′(1.75)=50.97398(2(1.75)−1.75−2.37(1.04−1.75)(1.04−2.37))+35.93855(2(1.75)−1.04−2.37(1.75−1.04)(1.75−2.37))+16.05180(2(1.75)−1.04−1.75(2.37−1.04)(2.37−1.75))=−26.99478

Calculate f′(2.37).

Substitute t=2.37 and i=1 in equation(1).

f′(2.37)=f(2.37)2(2.37)−t1−t2(t0−t1)(t0−t2)+f(t1)2(2.37)−t0−t2(t1−t0)(t1−t2)+f(t2)2(2.37)−t0−t1(t2−t0)(t2−t1)

Substitute the value of function from above table.

f′(2.37)=35.93855(2(2.37)−2.37−3.25(1.75−2.37)(1.75−3.25))+16.05180(2(2.37)−1.75−3.25(2.37−1.75)(2.37−3.25))+6.59273(2(2.37)−1.75−2.37(3.25−1.75)(3.25−2.37))=−23.26046

Calculate f′(3.25).

Substitute t=3.25 and i=1 in equation(1).

f′(3.25)=f(3.25)2(3.25)−t6−t7(t5−t6)(t5−t7)+f(t6)2(3.25)−t6−t7(t6−t5)(t6−t7)+f(t7)2(3.25)−t5−t6(t7−t6)(t7−t5)

Substitute the value of function from above table.

f′(3.25)=16.05180(2(3.25)−3.25−3.83(2.37−3.25)(2.37−3.83))+6.59273(2(3.25)−2.37−3.83(3.25−2.37)(3.25−3.83))+0.303817(2(3.25)−2.37−3.25(3.83−2.37)(3.83−3.25))=−10.80560

Calculate f′(3.83).

Substitute t=3.83 and i=7 in equation(1).

f′(3.83)=f(3.83)2(3.83)−t7−t8(t6−t7)(t6−t8)+f(t7)2(3.83)−t6−t8(t7−t6)(t7−t8)+f(t8)2(3.83)−t6−t8(t8−t6)(t8−t7)

Substitute the value of function from above table.

f′(3.83)=16.05180(2(3.83)−3.25−3.83(2.37−3.25)(2.37−3.83))+6.59273(2(3.83)−2.37−3.83(3.25−2.37)(3.25−3.83))+0.303817(2(3.83)−2.37−3.25(3.83−2.37)(3.83−3.25))=−10.88029

Want to see more full solutions like this?

Subscribe now to access step-by-step solutions to millions of textbook problems written by subject matter experts!
Previous
Chapter 24, Problem 43P
Next
Chapter 24, Problem 45P
Students have asked these similar questions
The speed of sound traveling through the sea is a function of temperature, salinity, and pressure. It is modeled by the function C = 1449.2 + 4.67T -0.055T² +0.00029T³ + (1.34-0.01T) (S-35) + 0.016D where C is the speed of sound (in meters per second), T is the temperature (in degrees Celsius), S is the salinity (the concentration of salts in parts per thousand, which means grams of dissolved solids per 1000 grams of water), and D is the depth below the sea surface, in meters. Calculate the following (Check attached image d3) when T = 10 °C, S = 35 parts per thousand, and D = 100 m. Explain the physical meaning of these derivatives.
the following data is obtained for the velocity of a vehicle during a crash test: t(ms)      v(mph) 0             30 10          29.5 20         28 30         23 40          10 50           5 60           2 70           0.5 80           0 if the vehicle's weight is 2,400 lb, approximate the force f in lb acting on the vehicle 40ms after the crash. the force can be calculated by  F=m(dv/dt) , and the mass of the car m is 2400/32.2 slug. note that 1 ms=10^-3 s and 1 mile = 5,280 ft. convert t to secons and mph to fps. use 1st order forward, backward, and centered difference and apply Richardson's approximation. use h1=0.02 and h2-0.01s.Show detailed solution. Thanks!
Problem statement : A simple flow-measurement device for streams and channels is a notch, of angle α, cut into the side of a dam, as shown in Fig. P5.69. The volume flow Q depends only on α, the acceleration of gravity g, and the height δ of the upstream water surface above the notch vertex. Tests of a model notch, of angle α = 55°, yield the following flow rate data: δ, cm: 10 20 30 40 Q, m3/h: 8 47 126 263 Part (A) - Find a dimensionless correlation for the data. (b) Use the model data to predict the flow rate of a prototype notch, also of angle α = 55°, when the upstream height δ is 3.2 m.

Chapter 24 Solutions

Numerical Methods for Engineers

Ch. 24 - Perform the same computation as Sec. 24.1, but...Ch. 24 - 24.2 Repeat Prob. 24.1, but use Romberg...Ch. 24 - 24.3 Repeat Prob. 24.1, but use a two- and a...Ch. 24 - 24.4 Integration provides a means to compute how...Ch. 24 - Use numerical integration to compute how much mass...Ch. 24 - 24.6 Fick’s first diffusion law states...Ch. 24 - The following data were collected when a large oil...Ch. 24 - 24.8 You are interested in measuring the fluid...Ch. 24 - Prob. 10PCh. 24 - 24.11 Glaucoma is the second leading cause of...
Ch. 24 - One of your colleagues has designed a new...Ch. 24 - Video an giography is used to measure blood flow...Ch. 24 - 24.14 Perform the same computation as in Sec....Ch. 24 - Perform the same computation as in Sec. 24.2, but...Ch. 24 - 24.16 As in Sec. 24.2, compute F using the...Ch. 24 - Stream cross-sectional areas (A) are required for...Ch. 24 - 24.18 As described in Prob. 24.17, the...Ch. 24 - 24.21 A transportation engineering study requires...Ch. 24 - 24.22 A wind force distributed against the side of...Ch. 24 - 24.23 Water exerts pressure on the upstream ...Ch. 24 - 24.24 To estimate the size of a new dam, you have...Ch. 24 - The data listed in the following table gives...Ch. 24 - The heat flux q is the quantity of heat flowing...Ch. 24 - 24.27 The horizontal surface area of a lake at a...Ch. 24 - 24.28 Perform the same computation as in Sec....Ch. 24 - 24.29 Repeat Prob. 24.28, but use five...Ch. 24 - Repeat Prob. 24.28, but use Romberg integration to...Ch. 24 - Faradays law characterizes the voltage drop across...Ch. 24 - 24.32 Based on Faraday’s law (Prob. 24.31), use...Ch. 24 - Suppose that the current through a resistor is...Ch. 24 - If a capacitor initially holds no charge, the...Ch. 24 - 24.35 Perform the same computation as in Sec....Ch. 24 - 24.36 Repeat Prob. 24.35, but use (a) Simpson’s ...Ch. 24 - 24.37 Compute work as described in Sec. 24.4, but...Ch. 24 - As was done in Sec. 24.4, determine the work...Ch. 24 - 24.39 The work done on an object is equal to the...Ch. 24 - The rate of cooling of a body (Fig. P24.40) can be...Ch. 24 - 24.41 A rod subject to an axial load (Fig....Ch. 24 - If the velocity distribution of a fluid flowing...Ch. 24 - 24.43 Using the following data, calculate the work...Ch. 24 - 24.44 A jet fighter’s position on an aircraft...Ch. 24 - 24.45 Employ the multiple-application Simpson’s...Ch. 24 - The upward velocity of a rocket can be computed by...Ch. 24 - Referring to the data from Problem 20.61, find the...Ch. 24 - Fully developed flow moving through a 40-cm...Ch. 24 - Fully developed flow of a Bingham plasticfluid...Ch. 24 - 24.50 The enthalpy of a real gas is a ...Ch. 24 - Given the data below, find the isothermal work...Ch. 24 - 24.52 The Rosin-Rammler-Bennet (RRB) equation is...Ch. 24 - For fluid flow over a surface, the heat flux to...Ch. 24 - The pressure gradient for laminar flow through a...Ch. 24 - 24.55 Velocity data for air are collected at...
Knowledge Booster
Background pattern image
Advanced Math
Learn more about
Need a deep-dive on the concept behind this application? Look no further. Learn more about this topic, advanced-math and related others by exploring similar questions and additional content below.
Similar questions
SEE MORE QUESTIONS
Recommended textbooks for you
Text book image
Elements Of Electromagnetics
Mechanical Engineering
ISBN:9780190698614
Author:Sadiku, Matthew N. O.
Publisher:Oxford University Press
Text book image
Mechanics of Materials (10th Edition)
Mechanical Engineering
ISBN:9780134319650
Author:Russell C. Hibbeler
Publisher:PEARSON
Text book image
Thermodynamics: An Engineering Approach
Mechanical Engineering
ISBN:9781259822674
Author:Yunus A. Cengel Dr., Michael A. Boles
Publisher:McGraw-Hill Education
Text book image
Control Systems Engineering
Mechanical Engineering
ISBN:9781118170519
Author:Norman S. Nise
Publisher:WILEY
Text book image
Mechanics of Materials (MindTap Course List)
Mechanical Engineering
ISBN:9781337093347
Author:Barry J. Goodno, James M. Gere
Publisher:Cengage Learning
Text book image
Engineering Mechanics: Statics
Mechanical Engineering
ISBN:9781118807330
Author:James L. Meriam, L. G. Kraige, J. N. Bolton
Publisher:WILEY
SEE MORE TEXTBOOKS
01 - What Is A Differential Equation in Calculus? Learn to Solve Ordinary Differential Equations.; Author: Math and Science;https://www.youtube.com/watch?v=K80YEHQpx9g;License: Standard YouTube License, CC-BY
Higher Order Differential Equation with constant coefficient (GATE) (Part 1) l GATE 2018; Author: GATE Lectures by Dishank;https://www.youtube.com/watch?v=ODxP7BbqAjA;License: Standard YouTube License, CC-BY
Solution of Differential Equations and Initial Value Problems; Author: Jefril Amboy;https://www.youtube.com/watch?v=Q68sk7XS-dc;License: Standard YouTube License, CC-BY