- Practice · 5 marks
Second derivative central difference scheme
Short answer
Central formula: (f_{i+1}−2f_i+f_{i-1})/h², order O(h²). For f=x⁴ at x=1.2, exact f''=17.28. With h=0.2: 17.36; with h=0.1: 17.30. Error reduces fourfold on halving h, confirming second-order accuracy.
Derive the central difference formula for the second derivative and state its order. For , estimate using and , and compare with the exact value.
Derivation
Add the Taylor series of and :
Solve for :
It is second-order accurate, and it uses three points (stencil ). It is the basic formula for diffusion terms.
Numerical for at
Exact: .
h = 0.2: , ,
h = 0.1: ,
| Estimate | Error | |
|---|---|---|
| 0.2 | 17.36 | +0.08 |
| 0.1 | 17.30 | +0.02 |
| exact | 17.28 | 0 |
The error falls by a factor of 4 when is halved, confirming second-order accuracy. It also matches theory: for .
Answer: () and (); exact 17.28.
✓ Answer checked by Claude Opus 5.5 (high effort)
Related questions
- Derive forward, backward, central differences · Applied Computational Fluid Dynamics, chapter 3
- Upwind versus central scheme for convection · Applied Computational Fluid Dynamics, chapter 3
- Boundary derivative with one-sided differences · Applied Computational Fluid Dynamics, chapter 3
- Numerical check of difference schemes · Applied Computational Fluid Dynamics, chapter 3
- Order of accuracy, consistency and stability · Applied Computational Fluid Dynamics, chapter 3
- Saint-Venant, second-order explicit scheme · Computational Techniques in Civil Engineering, chapter 5 · asked 2×
- Second-order scheme, linear kinematic wave · Computational Techniques in Civil Engineering, chapter 5
- Derivative feedback in second order time response · Control Systems, chapter 7
- First and second order finite difference schemes · Computational Techniques in Civil Engineering, chapter 5 · asked 2×
All Applied Computational Fluid Dynamics chapter 3 questions →