Study guide
2.1 Differential & Integral Calculus in Engineering Applications
PRC Mechanical Engineer Licensure Examination MELE
Track your progress as you read
Finished sections get a check mark.
1 Chapter 1: Licensing Overview, Law & Professional Ethics
1.1 MELE Facts, RA 8495 & Board Organization 1.2 Code of Ethics for Mechanical Engineers & Professional Practice 1.3 Philippine Mechanical Engineering Code (PME Code) Governance & Legal Mandates
2 Chapter 2: Advanced Mathematics & Engineering Statistics
2.1 Differential & Integral Calculus in Engineering Applications 2.2 Differential Equations & Mathematical Modeling 2.3 Applied Engineering Statistics & Quality Control Probability
3 Chapter 3: Engineering Thermodynamics
3.1 Thermodynamic State Postulate, Properties & P-V-T Equations of State 3.2 First Law Analysis for Open & Closed Thermal Systems 3.3 Second Law, Entropy, Exergy & Reversibility Analysis
4 Chapter 4: Fluid Mechanics & Hydraulic Engineering
4.1 Fluid Statics, Manometry, Hydrostatic Force & Buoyancy 4.2 Fluid Kinematics, Dynamics & Bernoulli Equation 4.3 Pipe Flow, Friction Losses & Piping Network Analysis
5 Chapter 5: Heat Transfer & Thermal Systems
5.1 Conduction Heat Transfer & Thermal Insulation Design 5.2 Natural & Forced Convection Heat Transfer Correlations 5.3 Thermal Radiation & Industrial Heat Exchanger Design
6 Chapter 6: Mechanics & Strength of Materials
6.1 Statics & Dynamics of Rigid Bodies & Engineering Linkages 6.2 Normal, Shear & Triaxial Stress-Strain Relationships 6.3 Bending Stresses, Torsion & Column Buckling Analysis
7 Chapter 7: Engineering Economics & Financial Analysis
7.1 Time Value of Money & Cash Flow Equivalence Formulas 7.2 Asset Depreciation, Taxation & Economic Evaluation Methods 7.3 Break-Even Analysis, Payback Period & EOQ Inventory Models
8 Chapter 8: Engineering Materials & Metallurgy
8.1 Ferrous & Non-Ferrous Metals, Alloys & Polymers 8.2 Phase Diagrams, Iron-Carbon System & Heat Treatment 8.3 Mechanical Testing, Fatigue Failure & Fracture Mechanics
9 Chapter 9: Machine Elements & Stress Analysis
9.1 Failure Theories, Static & Fluctuating Stress Design 9.2 Power Transmission Shafts, Keys & Rigid/Flexible Couplings 9.3 Spur, Helical, Bevel & Worm Gear Design Calculations
10 Chapter 10: Power Transmission Components
10.1 Belts, V-Belts, Chains, Wire Ropes & Flywheels 10.2 Rolling Element & Hydrodynamic Journal Bearings 10.3 Mechanical Springs, Clutches & Industrial Brakes
11 Chapter 11: Manufacturing, Fasteners & GD&T
11.1 Threaded Fasteners, Rivets & Welded Joint Analysis 11.2 Manufacturing Processes, Machine Tools & Machining Calculations 11.3 Engineering Drawing, Limits, Fits, Tolerances & GD&T
12 Chapter 12: Power Plants, Boilers & Thermal Power Cycles
12.1 Steam Power Plants & Rankine Cycles 12.2 Steam Generators, Boilers & Combustion Engineering 12.3 Internal Combustion Engines: SI & CI Cycles 12.4 Gas Turbine Power Plants & Combined Cycles
13 Chapter 13: Renewable Energy, Fluid Machinery & Plant Systems
13.1 Hydroelectric, Geothermal & Renewable Energy Power Plants 13.2 Pumps & Hydraulic Machinery 13.3 Air Compressors, Fans & Blowers
14 Chapter 14: Refrigeration, HVAC & Building Mechanical Systems
14.1 Refrigeration Systems, Cycles & Refrigerants 14.2 Psychrometrics & Air Conditioning Systems 14.3 Piping/Ductwork Design, PME Code & Industrial Safety Management
2.1 Differential & Integral Calculus in Engineering Applications
Key Takeaways Explain This Topic AI
Derivatives measure instantaneous rates of change, enabling optimization of physical dimensions, velocity/acceleration relationships, and thermal gradients.
Definite integrals quantify cumulative physical quantities such as polytropic gas work (W = integral P dV), fluid forces, centroids, and area moments of inertia.
Partial derivatives facilitate multi-variable sensitivity analysis and 2D/3D gradient evaluations in heat transfer and fluid flow.
Area centroids and moments of inertia (Ix = integral y^2 dA) define structural neutral axis locations and section moduli for beam stress calculations.
Numerical integration methods like Simpson's 1/3 Rule provide accurate work and area approximations when experimental data points are discrete.
Last updated: July 2026
- 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 ) = d s d t v(t) = \frac{ds}{dt} v ( t ) = d t d s and linear acceleration a ( t ) = d v d t = d 2 s d t 2 = v d v d s a(t) = \frac{dv}{dt} = \frac{d^2s}{dt^2} = v \frac{dv}{ds} a ( t ) = d t d v = d t 2 d 2 s = v d s d v . Angular velocity ω ( t ) = d θ d t \omega(t) = \frac{d\theta}{dt} ω ( t ) = d t d θ and angular acceleration α ( t ) = d ω d t \alpha(t) = \frac{d\omega}{dt} α ( t ) = d t d ω .
Thermodynamic Heat Flux : Heat rate q ( t ) = d Q d t q(t) = \frac{dQ}{dt} q ( t ) = d t d Q , where power P P P is the derivative of work with respect to time ( P = d W d t = F ⃗ ⋅ v ⃗ P = \frac{dW}{dt} = \vec{F} \cdot \vec{v} P = d t d W = F ⋅ v ).
Fluid Dynamics : Velocity gradients d u d y \frac{du}{dy} d y d u governing shear stress τ = μ d u d y \tau = \mu \frac{du}{dy} τ = μ d y d u 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 0 0 \frac{0}{0} 0 0 or ∞ ∞ \frac{\infty}{\infty} ∞ ∞ frequently arise. L'Hôpital's Rule states: lim x → c f ( x ) g ( x ) = lim x → c f ′ ( x ) g ′ ( x ) \lim_{x \to c} \frac{f(x)}{g(x)} = \lim_{x \to c} \frac{f'(x)}{g'(x)} x → c lim g ( x ) f ( x ) = x → c lim g ′ ( x ) f ′ ( x )
provided g ′ ( x ) ≠ 0 g'(x) \neq 0 g ′ ( x ) = 0 near c c 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:
Express the objective function in terms of a single independent variable using geometric or physical constraint equations.
Calculate the first derivative d y d x \frac{dy}{dx} d x d y and locate critical points by setting d y d x = 0 \frac{dy}{dx} = 0 d x d y = 0 .
Perform the Second Derivative Test:
If d 2 y d x 2 > 0 \frac{d^2y}{dx^2} > 0 d x 2 d 2 y > 0 at the critical point, the function achieves a local minimum .
If d 2 y d x 2 < 0 \frac{d^2y}{dx^2} < 0 d x 2 d 2 y < 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 V_0 V 0 . Determine the optimal ratio of length h h h to diameter d d d ( h / d h/d h / d ) that minimizes the total surface area S S S (and hence material cost).
Step 1: Constraint & Objective Equations
Volume constraint: V = π r 2 h ⟹ h = V π r 2 V = \pi r^2 h \implies h = \frac{V}{\pi r^2} V = π r 2 h ⟹ h = π r 2 V
Surface Area objective: S = 2 π r 2 + 2 π r h S = 2\pi r^2 + 2\pi r h S = 2 π r 2 + 2 π r h
Step 2: Substitute Constraint into Objective S ( r ) = 2 π r 2 + 2 π r ( V π r 2 ) = 2 π r 2 + 2 V r S(r) = 2\pi r^2 + 2\pi r \left(\frac{V}{\pi r^2}\right) = 2\pi r^2 + \frac{2V}{r} S ( r ) = 2 π r 2 + 2 π r ( π r 2 V ) = 2 π r 2 + r 2 V
Step 3: Differentiate and Set to Zero d S d r = 4 π r − 2 V r 2 = 0 ⟹ 4 π r 3 = 2 V ⟹ r 3 = V 2 π \frac{dS}{dr} = 4\pi r - \frac{2V}{r^2} = 0 \implies 4\pi r^3 = 2V \implies r^3 = \frac{V}{2\pi} d r d S = 4 π r − r 2 2 V = 0 ⟹ 4 π r 3 = 2 V ⟹ r 3 = 2 π V
Step 4: Express Height in terms of Radius Substitute V = 2 π r 3 V = 2\pi r^3 V = 2 π r 3 back into the height equation: h = 2 π r 3 π r 2 = 2 r h = \frac{2\pi r^3}{\pi r^2} = 2r h = π r 2 2 π r 3 = 2 r
Since diameter d = 2 r d = 2r d = 2 r , we have: h = d ⟹ h d = 1.0 h = d \implies \frac{h}{d} = 1.0 h = d ⟹ d h = 1.0
Step 5: Verify Minimum via Second Derivative d 2 S d r 2 = 4 π + 4 V r 3 > 0 ( since r , V > 0 ) \frac{d^2S}{dr^2} = 4\pi + \frac{4V}{r^3} > 0 \quad (\text{since } r, V > 0) d r 2 d 2 S = 4 π + r 3 4 V > 0 ( since r , V > 0 )
Thus, a cylinder whose height equals its diameter ( h / d = 1.0 h/d = 1.0 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 ) h(P, T) h ( P , T ) depends on pressure and temperature.
Key Partial Derivative Relationships:
Partial Derivative : ∂ f ∂ x \frac{\partial f}{\partial x} ∂ x ∂ f differentiates f ( x , y ) f(x, y) f ( x , y ) holding y y y constant.
Total Differential : d f = ∂ f ∂ x d x + ∂ f ∂ y d y df = \frac{\partial f}{\partial x} dx + \frac{\partial f}{\partial y} dy df = ∂ x ∂ f d x + ∂ y ∂ f d y , used for error propagation in experimental measurements.
Gradient Vector : ∇ T = ∂ T ∂ x i ^ + ∂ T ∂ y j ^ + ∂ T ∂ z k ^ \nabla T = \frac{\partial T}{\partial x}\hat{i} + \frac{\partial T}{\partial y}\hat{j} + \frac{\partial T}{\partial z}\hat{k} ∇ T = ∂ x ∂ T i ^ + ∂ y ∂ T j ^ + ∂ z ∂ T k ^ , defining spatial heat flux via Fourier's Law q ⃗ ′ ′ = − k ∇ T \vec{q}'' = -k \nabla T q ′′ = − k ∇ T .
4. Integral Calculus & Physical Applications
Definite integration evaluates cumulative physical quantities over spatial domains or temporal intervals.
Physical Metric Integral Formula Units (SI)
Area under Curve A = ∫ a b [ f ( x ) − g ( x ) ] d x A = \int_{a}^{b} [f(x) - g(x)] dx A = ∫ a b [ f ( x ) − g ( x )] d x m 2 \text{m}^2 m 2
Polytropic Gas Work W = ∫ V 1 V 2 P d V W = \int_{V_1}^{V_2} P \, dV W = ∫ V 1 V 2 P d V J \text{J} J or kJ \text{kJ} kJ
Hydrostatic Force F = ∫ ρ g h d A F = \int \rho g h \, dA F = ∫ ρ g h d A N \text{N} N
Centroid Coordinate x ˉ = 1 A ∫ x d A , y ˉ = 1 A ∫ y d A \bar{x} = \frac{1}{A} \int x \, dA, \quad \bar{y} = \frac{1}{A} \int y \, dA x ˉ = A 1 ∫ x d A , y ˉ = A 1 ∫ y d A m \text{m} m
Area Moment of Inertia I x = ∫ y 2 d A , I y = ∫ x 2 d A I_x = \int y^2 \, dA, \quad I_y = \int x^2 \, dA I x = ∫ y 2 d A , I y = ∫ x 2 d A m 4 \text{m}^4 m 4
Worked Numerical Example: Polytropic Gas Expansion Work
A gas in a piston-cylinder assembly expands polytropically ( P V n = C P V^{n} = C P V n = C ) with exponent n = 1.3 n = 1.3 n = 1.3 . The initial state is P 1 = 500 kPa P_1 = 500 \text{ kPa} P 1 = 500 kPa and V 1 = 0.1 m 3 V_1 = 0.1 \text{ m}^3 V 1 = 0.1 m 3 . The final volume is V 2 = 0.3 m 3 V_2 = 0.3 \text{ m}^3 V 2 = 0.3 m 3 . Calculate the total boundary work done by the gas in kilojoules (kJ).
Step 1: Calculate Final Pressure P 2 P_2 P 2 P 1 V 1 n = P 2 V 2 n ⟹ P 2 = P 1 ( V 1 V 2 ) n = 500 ( 0.1 0.3 ) 1.3 = 500 ( 1 3 ) 1.3 = 119.72 kPa 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} P 1 V 1 n = P 2 V 2 n ⟹ P 2 = P 1 ( V 2 V 1 ) n = 500 ( 0.3 0.1 ) 1.3 = 500 ( 3 1 ) 1.3 = 119.72 kPa
Step 2: Apply Polytropic Work Integral Formula W = ∫ V 1 V 2 P d V = C ∫ V 1 V 2 V − n d V = P 1 V 1 − P 2 V 2 n − 1 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} W = ∫ V 1 V 2 P d V = C ∫ V 1 V 2 V − n d V = n − 1 P 1 V 1 − P 2 V 2
Step 3: Substitute Values W = ( 500 kPa ) ( 0.1 m 3 ) − ( 119.72 kPa ) ( 0.3 m 3 ) 1.3 − 1 W = \frac{(500 \text{ kPa})(0.1 \text{ m}^3) - (119.72 \text{ kPa})(0.3 \text{ m}^3)}{1.3 - 1} W = 1.3 − 1 ( 500 kPa ) ( 0.1 m 3 ) − ( 119.72 kPa ) ( 0.3 m 3 ) W = 50.00 − 35.916 0.3 = 14.084 kPa ⋅ m 3 0.3 = 46.95 kJ W = \frac{50.00 - 35.916}{0.3} = \frac{14.084 \text{ kPa}\cdot\text{m}^3}{0.3} = 46.95 \text{ kJ} W = 0.3 50.00 − 35.916 = 0.3 14.084 kPa ⋅ m 3 = 46.95 kJ
Thus, the total work output during expansion is 46.95 kJ 46.95 \text{ kJ} 46.95 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 σ = M y I \sigma = \frac{M y}{I} σ = I M y .
Parallel Axis Theorem : I x = I x c + A d y 2 I_x = I_{xc} + A d_y^2 I x = I x c + A d y 2 where I x c I_{xc} I x c is the moment of inertia about the centroidal axis, A A A is the cross-sectional area, and d y d_y d y is the perpendicular distance between parallel axes.
6. Numerical Integration (Trapezoidal & Simpson's Rule)
When an analytical function f ( x ) f(x) 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 n n subintervals of width h = b − a n h = \frac{b-a}{n} h = n b − a ): ∫ a b f ( x ) d x ≈ h 2 [ y 0 + 2 y 1 + 2 y 2 + ⋯ + 2 y n − 1 + y 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] ∫ a b f ( x ) d x ≈ 2 h [ y 0 + 2 y 1 + 2 y 2 + ⋯ + 2 y n − 1 + y n ]
Simpson's 1/3 Rule ( n n n MUST be an EVEN number of subintervals): ∫ a b f ( x ) d x ≈ h 3 [ y 0 + 4 ( y 1 + y 3 + … ) + 2 ( y 2 + y 4 + … ) + y n ] \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] ∫ a b f ( x ) d x ≈ 3 h [ y 0 + 4 ( y 1 + y 3 + … ) + 2 ( y 2 + y 4 + … ) + y n ]
Worked Numerical Example: Numerical Work Calculation
A machine engine test yields force F ( x ) F(x) F ( x ) at 0.5 m 0.5 \text{ m} 0.5 m intervals along a displacement path x = 0 x = 0 x = 0 to x = 2.0 m x = 2.0 \text{ m} x = 2.0 m as given below:
- x x x (m) 0.0 0.5 1.0 1.5 2.0
- F ( x ) F(x) F ( x ) (kN) 10.0 18.0 25.0 20.0 12.0
Calculate work W = ∫ 0 2.0 F ( x ) d x W = \int_0^{2.0} F(x) dx W = ∫ 0 2.0 F ( x ) d x using Simpson's 1/3 Rule.
Step 1: Identify Parameters
Subinterval step size h = 0.5 m h = 0.5 \text{ m} h = 0.5 m .
n = 4 n = 4 n = 4 subintervals (even, so Simpson's 1/3 Rule applies).
y 0 = 10.0 y_0 = 10.0 y 0 = 10.0 , y 1 = 18.0 y_1 = 18.0 y 1 = 18.0 , y 2 = 25.0 y_2 = 25.0 y 2 = 25.0 , y 3 = 20.0 y_3 = 20.0 y 3 = 20.0 , y 4 = 12.0 y_4 = 12.0 y 4 = 12.0 .
Step 2: Apply Simpson's Formula W ≈ 0.5 3 [ 10.0 + 4 ( 18.0 + 20.0 ) + 2 ( 25.0 ) + 12.0 ] W \approx \frac{0.5}{3} \left[ 10.0 + 4(18.0 + 20.0) + 2(25.0) + 12.0 \right] W ≈ 3 0.5 [ 10.0 + 4 ( 18.0 + 20.0 ) + 2 ( 25.0 ) + 12.0 ] W ≈ 0.5 3 [ 10.0 + 4 ( 38.0 ) + 50.0 + 12.0 ] = 0.5 3 [ 10.0 + 152.0 + 50.0 + 12.0 ] 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 ≈ 3 0.5 [ 10.0 + 4 ( 38.0 ) + 50.0 + 12.0 ] = 3 0.5 [ 10.0 + 152.0 + 50.0 + 12.0 ] W ≈ 0.5 3 [ 224.0 ] = 37.33 kN ⋅ m = 37.33 kJ W \approx \frac{0.5}{3} \left[ 224.0 \right] = 37.33 \text{ kN}\cdot\text{m} = 37.33 \text{ kJ} W ≈ 3 0.5 [ 224.0 ] = 37.33 kN ⋅ m = 37.33 kJ
MELE Exam Tips & Shortcuts
Radian Mode Warning : When taking derivatives or integrals of trigonometric expressions ( d ( sin ω t ) / d t d(\sin \omega t)/dt d ( sin ω t ) / d t ), ensure your calculator is set to RADIANS .
Polytropic Integrals : For isothermal processes ( n = 1 n = 1 n = 1 ), the polytropic formula P 1 V 1 − P 2 V 2 n − 1 \frac{P_1 V_1 - P_2 V_2}{n-1} n − 1 P 1 V 1 − P 2 V 2 yields division by zero! Use W = P 1 V 1 ln ( V 2 V 1 ) W = P_1 V_1 \ln\left(\frac{V_2}{V_1}\right) W = P 1 V 1 ln ( V 1 V 2 ) instead.
Simpson's Rule Condition : Always check if n n n is even (or if the number of data points is odd). If n n 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.