---
title: "Short note: Smoothed particle hydrodynamics"
subject: "Computational Techniques in Civil Engineering"
chapter: 1
key: "semester_8/computational_techniques_in_civil_engineering:short-note-sph"
source: https://www.bctengineer.com/q/semester_8/computational_techniques_in_civil_engineering/short-note-sph
site: BCT Engineering
---

# Short note: Smoothed particle hydrodynamics

Computational Techniques in Civil Engineering · Chapter 1: Introduction

- Asked in: 2078 Kartik, 3 marks
- Video: [SPH](https://www.youtube.com/watch?v=lQ3a77WrNAU) — Civil Engineering UGent
- Answer checked by Claude Opus 5.5 (high effort) on 2026-10-14
- Page: https://www.bctengineer.com/q/semester_8/computational_techniques_in_civil_engineering/short-note-sph
- Whole chapter: https://www.bctengineer.com/md/semester_8/computational_techniques_in_civil_engineering/chapter-1

## Question

Write a short note on smoothed particle hydrodynamics.

## Short answer

SPH is a mesh-free Lagrangian particle method representing the medium by particles with mass, density and velocity. Field values are interpolated by a smoothing kernel over neighbouring particles, and governing equations become ordinary differential equations integrated in time. It handles large deformation and free surfaces, but has costly neighbour search and boundary treatment difficulties.

## Answer

**Smoothed particle hydrodynamics (SPH)** is a mesh-free, Lagrangian particle method originally developed for astrophysics (Lucy, and Gingold and Monaghan, 1977) and now used for fluids and solids with large deformation. The medium is represented by particles carrying mass, density, velocity and pressure; there is no grid.

**Principle:** the value of a field $f$ at a point is the weighted sum of the values at neighbouring particles within a smoothing length $h$:
$$
f(x_i)\approx\sum_j \frac{m_j}{\rho_j}\,f_j\,W(x_i-x_j,h)
$$
where $W$ is the smoothing (kernel) function. Gradients are obtained by differentiating $W$, so the governing equations (continuity and momentum) become ordinary differential equations for each particle, which are integrated in time.

**Advantages:** no mesh distortion; free surfaces and large deformation are handled easily; mass is conserved automatically (each particle carries a fixed mass).

**Disadvantages:** costly neighbour search; boundary treatment and tensile instability are difficult; lower accuracy than FEM for small deformation.

**Applications:** dam break, wave impact, flood flow, landslide and debris flow, soil flow.

## Check yourself

The answer's key points as multiple-choice questions, with the correct option and why.

1. What is SPH?
   A. A boundary integral method
   B. A mesh-free Lagrangian particle method
   C. A grid-based finite difference method
   D. An assembly of rigid blocks
   Answer: B. SPH represents the medium by particles carrying mass, density, velocity and pressure, with no grid. (Key point: SPH nature)

2. How is a field value at a point obtained in SPH?
   A. As a weighted sum of neighbouring particle values within the smoothing length h using a kernel W
   B. By solving [K]{U} = {F}
   C. By Green's function on the boundary
   D. By differences on a grid
   Answer: A. f(xi) is approximated by sum over neighbours of (mj/rhoj) fj W. (Key point: SPH kernel)

3. Why is SPH good for dam-break and debris flow?
   A. It is more accurate than FEM for small deformation
   B. It needs only boundary discretisation
   C. It has cheap neighbour search
   D. No mesh to distort, so free surfaces and large deformation are handled easily
   Answer: D. Being mesh-free, SPH avoids mesh distortion; neighbour search is actually costly. (Key point: SPH advantage)