2.1 Differential & Integral Calculus in Engineering Applications

2.1 Differential & Integral Calculus in Engineering Applications

Core Calculus Concepts in Mechanical Engineering

Differential and integral calculus serve as the fundamental mathematical framework for describing mechanical phenomena involving continuous motion, variable force fields, rates of thermal transfer, and material geometry. On the PRC Mechanical Engineering Licensure Examination (MELE), calculus problems assess both direct mathematical operations and their physical applications in kinematics, structural design, fluid dynamics, and thermodynamics. Mastery of these topics requires rapid identification of rate-of-change models, optimization techniques, double/triple integration for physical centroids, and numerical integration algorithms.

---

1. Differential Calculus & Rates of Change

In mechanical systems, derivative functions represent instantaneous physical rates:

• Kinematics: Linear velocity $v(t) = \frac{ds}{dt}$ and linear acceleration $a(t) = \frac{dv}{dt} = \frac{d^2s}{dt^2} = v \frac{dv}{ds}$. Angular velocity $\omega(t) = \frac{d\theta}{dt}$ and angular acceleration $\alpha(t) = \frac{d\omega}{dt}$.
• Thermodynamic Heat Flux: Heat rate $q(t) = \frac{dQ}{dt}$, where power $P$ is the derivative of work with respect to time ($P = \frac{dW}{dt} = \vec{F} \cdot \vec{v}$).
• Fluid Dynamics: Velocity gradients $\frac{du}{dy}$ governing shear stress $\tau = \mu \frac{du}{dy}$ in Newtonian fluids.

Indeterminate Forms & L'Hôpital's Rule

When evaluating limits during fluid boundary layer analysis or asymptotic thermal decay, indeterminate forms such as $\frac{0}{0}$ or $\frac{\infty}{\infty}$ frequently arise. L'Hôpital's Rule states: $\lim_{x \to c} \frac{f(x)}{g(x)} = \lim_{x \to c} \frac{f'(x)}{g'(x)}$ provided $g'(x) \neq 0$ near $c$ and the limit of derivatives exists.

---

2. Engineering Optimization: Maxima & Minima

Designing mechanical components requires optimizing physical metrics such as minimizing material volume, weight, or fluid pressure drop, while maximizing volumetric capacity, structural stiffness, or heat dissipation.

Optimization Workflow:

1. Express the objective function in terms of a single independent variable using geometric or physical constraint equations.
2. Calculate the first derivative $\frac{dy}{dx}$ and locate critical points by setting $\frac{dy}{dx} = 0$.
3. Perform the Second Derivative Test:
   • If $\frac{d^2y}{dx^2} > 0$ at the critical point, the function achieves a local minimum.
   • If $\frac{d^2y}{dx^2} < 0$ at the critical point, the function achieves a local maximum.

Worked Numerical Example: Optimal Cylindrical Pressure Vessel

A cylindrical storage tank with closed flat ends is to be fabricated to hold a specified volume $V_0$. Determine the optimal ratio of length $h$ to diameter $d$ ($h/d$) that minimizes the total surface area $S$ (and hence material cost).

Step 1: Constraint & Objective Equations

• Volume constraint: $V = \pi r^2 h \implies h = \frac{V}{\pi r^2}$
• Surface Area objective: $S = 2\pi r^2 + 2\pi r h$

Step 2: Substitute Constraint into Objective $S(r) = 2\pi r^2 + 2\pi r \left(\frac{V}{\pi r^2}\right) = 2\pi r^2 + \frac{2V}{r}$

Step 3: Differentiate and Set to Zero $\frac{dS}{dr} = 4\pi r - \frac{2V}{r^2} = 0 \implies 4\pi r^3 = 2V \implies r^3 = \frac{V}{2\pi}$

Step 4: Express Height in terms of Radius Substitute $V = 2\pi r^3$ back into the height equation: $h = \frac{2\pi r^3}{\pi r^2} = 2r$

Since diameter $d = 2r$, we have: $h = d \implies \frac{h}{d} = 1.0$

Step 5: Verify Minimum via Second Derivative $\frac{d^2S}{dr^2} = 4\pi + \frac{4V}{r^3} > 0 \quad (\text{since } r, V > 0)$ Thus, a cylinder whose height equals its diameter ($h/d = 1.0$) minimizes the total material surface area for a given internal volume.

---

3. Partial Derivatives & Multi-Variable Calculus

Mechanical systems often depend on multiple independent state variables. For example, specific enthalpy $h(P, T)$ depends on pressure and temperature.

Key Partial Derivative Relationships:

• Partial Derivative: $\frac{\partial f}{\partial x}$ differentiates $f(x, y)$ holding $y$ constant.
• Total Differential: $df = \frac{\partial f}{\partial x} dx + \frac{\partial f}{\partial y} dy$, used for error propagation in experimental measurements.
• Gradient Vector: $\nabla T = \frac{\partial T}{\partial x}\hat{i} + \frac{\partial T}{\partial y}\hat{j} + \frac{\partial T}{\partial z}\hat{k}$, defining spatial heat flux via Fourier's Law $\vec{q}'' = -k \nabla T$.

---

4. Integral Calculus & Physical Applications

Definite integration evaluates cumulative physical quantities over spatial domains or temporal intervals.

Worked Numerical Example: Polytropic Gas Expansion Work

A gas in a piston-cylinder assembly expands polytropically ($P V^{n} = C$) with exponent $n = 1.3$. The initial state is $P_1 = 500 \text{ kPa}$ and $V_1 = 0.1 \text{ m}^3$. The final volume is $V_2 = 0.3 \text{ m}^3$. Calculate the total boundary work done by the gas in kilojoules (kJ).

Step 1: Calculate Final Pressure $P_2$ $P_1 V_1^n = P_2 V_2^n \implies P_2 = P_1 \left(\frac{V_1}{V_2}\right)^n = 500 \left(\frac{0.1}{0.3}\right)^{1.3} = 500 \left(\frac{1}{3}\right)^{1.3} = 119.72 \text{ kPa}$

Step 2: Apply Polytropic Work Integral Formula $W = \int_{V_1}^{V_2} P \, dV = C \int_{V_1}^{V_2} V^{-n} dV = \frac{P_1 V_1 - P_2 V_2}{n - 1}$

Step 3: Substitute Values $W = \frac{(500 \text{ kPa})(0.1 \text{ m}^3) - (119.72 \text{ kPa})(0.3 \text{ m}^3)}{1.3 - 1}$ $W = \frac{50.00 - 35.916}{0.3} = \frac{14.084 \text{ kPa}\cdot\text{m}^3}{0.3} = 46.95 \text{ kJ}$

Thus, the total work output during expansion is $46.95 \text{ kJ}$.

---

5. Centroids & Area Moments of Inertia

The centroid of a plane area represents its geometric center, which defines the neutral axis location for beam bending stress $\sigma = \frac{M y}{I}$.

• Parallel Axis Theorem: $I_x = I_{xc} + A d_y^2$ where $I_{xc}$ is the moment of inertia about the centroidal axis, $A$ is the cross-sectional area, and $d_y$ is the perpendicular distance between parallel axes.

---

6. Numerical Integration (Trapezoidal & Simpson's Rule)

When an analytical function $f(x)$ is unavailable and data is obtained from experimental testing (e.g., indicator diagrams of pressure vs volume), numerical quadrature methods are required.

Trapezoidal Rule ($n$ subintervals of width $h = \frac{b-a}{n}$):

$\int_{a}^{b} f(x) dx \approx \frac{h}{2} \left[ y_0 + 2y_1 + 2y_2 + \dots + 2y_{n-1} + y_n \right]$

Simpson's 1/3 Rule ($n$ MUST be an EVEN number of subintervals):

$\int_{a}^{b} f(x) dx \approx \frac{h}{3} \left[ y_0 + 4(y_1 + y_3 + \dots) + 2(y_2 + y_4 + \dots) + y_n \right]$

Worked Numerical Example: Numerical Work Calculation

A machine engine test yields force $F(x)$ at $0.5 \text{ m}$ intervals along a displacement path $x = 0$ to $x = 2.0 \text{ m}$ as given below:

Calculate work $W = \int_0^{2.0} F(x) dx$ using Simpson's 1/3 Rule.

Step 1: Identify Parameters

• Subinterval step size $h = 0.5 \text{ m}$.
• $n = 4$ subintervals (even, so Simpson's 1/3 Rule applies).
• $y_0 = 10.0$, $y_1 = 18.0$, $y_2 = 25.0$, $y_3 = 20.0$, $y_4 = 12.0$.

Step 2: Apply Simpson's Formula $W \approx \frac{0.5}{3} \left[ 10.0 + 4(18.0 + 20.0) + 2(25.0) + 12.0 \right]$ $W \approx \frac{0.5}{3} \left[ 10.0 + 4(38.0) + 50.0 + 12.0 \right] = \frac{0.5}{3} \left[ 10.0 + 152.0 + 50.0 + 12.0 \right]$ $W \approx \frac{0.5}{3} \left[ 224.0 \right] = 37.33 \text{ kN}\cdot\text{m} = 37.33 \text{ kJ}$

---

MELE Exam Tips & Shortcuts

• Radian Mode Warning: When taking derivatives or integrals of trigonometric expressions ($d(\sin \omega t)/dt$), ensure your calculator is set to RADIANS.
• Polytropic Integrals: For isothermal processes ($n = 1$), the polytropic formula $\frac{P_1 V_1 - P_2 V_2}{n-1}$ yields division by zero! Use $W = P_1 V_1 \ln\left(\frac{V_2}{V_1}\right)$ instead.
• Simpson's Rule Condition: Always check if $n$ is even (or if the number of data points is odd). If $n$ is odd, Simpson's 1/3 Rule cannot be applied across the entire domain without combining with the Trapezoidal Rule for the last panel.