CoMPhy Lab blogs

14 min read

📄 PDF Version

Download PDF

Lecture 2 -

PDF version

Download this page as PDF

From Newton’s Laws to Conservation Equations

Learning Goals

By the end of this lecture, you should understand:

  1. Newton’s second law for one particle: , and how it extends to many particles
  2. Microscopic fields from particles: how to convert a sum of point particles into continuous density fields using Dirac deltas
  3. The moving delta identity: (see detailed note)
  4. Microscopic continuity equation: how mass conservation emerges at the particle level
  5. Pairwise forces → divergence of stress: how pair forces become a flux tensor (see detailed note)
  6. Microscopic momentum balance:
  7. Coarse-graining → Cauchy stress: how ensemble averaging yields the Cauchy momentum equation
  8. Physical intuition: internal forces appear as stresses (forces per unit area) in the continuum description

Part A: Newton’s Second Law — Scalar to Many Particles

Single particle (warm-up)

Write on board:

Many particles (N particles, pairwise forces)

Write:

Important

We start here: every particle obeys Newton’s second law with external forces and pair interactions. The key assumption — which is exact in classical mechanics — is that the force one particle exerts on another is equal and opposite (, Newton’s third law). This symmetry is crucial: it will let us turn distributed forces into a stress tensor.

  1. : mass of particle
  2. : velocity of particle
  3. : external force (gravity, walls, etc.)
  4. : force exerted by particle on particle
  5. : action–reaction; internal forces cancel when you sum over all particles

Part B: Convert Particles to Fields Using Dirac Deltas

Big idea

Info

We want to describe the fluid at every point in space, not just at discrete particle locations. To do that, we represent each particle as a moving point mass — a Dirac delta spike at . Then we can define continuous fields: density, momentum, stress, etc.

Setup: Microscopic fields

Density:

Momentum:

where is a spike centered on particle .

TL;DR

Each particle contributes a delta spike to the density field. At the location , the spike has ‘height’ (or weight) equal to the particle’s mass. The momentum density does the same but weights by . These fields are exact — there’s no averaging yet. But they’re spiky and discontinuous; what we want eventually is a smooth continuum limit.

Part C: Microscopic Continuity — Mass Conservation at the Particle Level

The key identity: moving delta

Important

To convert particle equations into field equations, we need one crucial identity: how does a Dirac delta change when the particle moves? This is the moving delta identity (see detailed derivation).

This is the advection equation for a spike. The time rate of change at a fixed point is minus the spatial change (gradient) in the direction the spike moves. See detailed derivation for three ways to justify this.

Derive mass continuity

Start with the definition of :

Use the moving delta identity:

Factor out the gradient (it acts only on , not on ):

Recognize the sum as the momentum density :

Tip

That one identity—the moving delta—is the entire algebra. It converts a sum of moving spikes into a divergence of momentum density. This is mass conservation at the microscopic level: if mass flows out of a region at rate , the local density must decrease at rate . No approximation yet—this is exact for point particles.

Part D: Microscopic Momentum Balance — Where Stress Is Born

Idea

Now we take the time derivative of the momentum density. Newton’s second law says particles accelerate due to the net external forces. When we convert that into a field equation, the internal (pair) forces turn into a divergence of a tensor—and that tensor is the stress.

Step 1: Time derivative of momentum density

Use the product rule:

The first term is the acceleration field; the second is what needs conversion. See detailed note for the algebra, but the result is:

where is the outer product (dyad): .

Define the kinetic momentum-flux tensor:

So:

Step 2: Insert Newton’s second law

From Newton’s 2nd law:

So:

where is the external body force per unit mass (e.g., gravity), and is the external force per unit volume.

Step 3: Turn pairwise forces into a divergence of the virial tensor

This is the non-obvious step. See detailed derivation for the full proof, but the key result is:

where is the configurational (virial) momentum-flux tensor:

Note:

Each pair force is ‘smeared’ uniformly along the bond connecting the two particles. When you take the divergence, you recover the two point forces at the ends—equal and opposite, as required. The integral parameterizes the bond; it’s just a mathematical trick to convert point forces into a flux.

Important

Internal forces don’t vanish when summed over all particles, but they appear locally as a divergence of a stress-like tensor. This is the transition from particles to continuum: pair forces become components of the stress tensor.

Step 4: Complete microscopic momentum balance

Gather all terms:

Define the total microscopic momentum-flux tensor:

So:

Important

This is exact: no averaging, no approximation. Every term comes from Newton’s second law and the geometry of deltas. The tensor contains two pieces: kinetic (, from bulk motion and random thermal motion) and configurational (, from pair forces). Together, they form the microscopic stress.

Part E: Coarse-Grain to Get Cauchy’s Momentum Equation

Idea

The microscopic equation is spiky and exact, but a real fluid isn’t made of point particles—we have particles per cm. So we ensemble-average (or volume-average) over a mesoscopic region that contains many particles but is still small enough to resolve macroscopic gradients.

Define macroscopic fields via averaging

See detailed derivation for the complete mathematical treatment.

where the angle brackets denote ensemble (or volume) averaging.

Key observation: The momentum density has two components:

  1. Bulk momentum: (all particles drifting in the same direction)
  2. Random (thermal) momentum: from velocity fluctuations around the mean

Split the flux: convection + stress

Write the crucial definition:

where is the Cauchy stress tensor:

Breaking it down:

  • : bulk contribution to stress
  • : configurational (virial) contribution to stress (from pair forces)

Important

The stress tensor is defined as the part of the momentum flux that isn’t accounted for by the bulk drift. It captures how internal forces and thermal motion resist (or facilitate) deformation. For a gas at rest, it’s just pressure. For a flowing fluid under shear, it includes viscous stresses.

Derive Cauchy’s momentum equation

Average the microscopic balance:

Substitute and :

Use the continuity equation to simplify:

where is the material (convective) derivative.

TL;DR

This is the continuum momentum balance. It says: the rate of change of momentum in a moving fluid element equals the net stress applied to it plus external forces. It’s valid for any fluid, whether Newtonian (constant viscosity) or non-Newtonian. The specific form of is the constitutive relation, which depends on the fluid’s internal structure.

Part F: Why This Matters—Physical Intuition

The stress tensor is the flux of momentum

Draw on the board:

    → → → (fast fluid)
    __________
   |          |  interface (unit area, normal n)
   |__________|
    → → (slow fluid)

Imagine a small surface across the interface. Momentum flows across it in two ways: (1) bulk convection ( term), and (2) microscopic collisions creating a stress ( term). The total flux is the product force per unit area, i.e., traction = . If the fast fluid pushes on the slow fluid, the traction is positive; the stress tensor quantifies this.

Example: pressure (isotropic stress)

For a fluid at rest or in hydrostatic equilibrium:

where is pressure and is the identity. Then:

Pressure is the simplest stress: it acts equally in all directions. The gradient pushes fluid from high to low pressure. No viscosity, no shear—just pressure.

Example: viscous stress (Newtonian fluid)

For a Newtonian fluid with viscosity :

where is the viscous stress (from internal friction and momentum diffusion):

(The second term involves bulk viscosity ; for many fluids it’s negligible.)

The viscous stress is proportional to the strain rate . Faster shear means stronger viscous stress. This is the microscopic momentum diffusion we derived in Lecture 2: Viscosity—the kinetic theory. But here we’ve derived it from first principles: a consequence of particle interactions and collisions.

Summary

Common Misconceptions

MisconceptionCorrection
Dirac deltas are just very tall, thin spikes.They’re distributions: singular objects defined by their action on test functions. But physically, you can think of them as limits of smooth bumps.
The moving delta identity is obvious.It’s not! It requires careful treatment of the chain rule (or test functions). See the detailed note.
Internal forces cancel, so they don’t matter.They cancel globally, but they create local stresses and are essential for momentum balance in the fluid.
Stress is the same as pressure.Pressure is an isotropic stress (equal in all directions). Stress is more general: it can be anisotropic (different in different directions) and depends on the flow.
Coarse-graining is an approximation.It’s not an approximation—it’s a change of perspective. The microscopic equation is exact but spiky. Averaging converts it into a smooth, usable form.
The Cauchy equation only works for Newtonian fluids.No. It’s exact for any fluid. The assumption of Newtonian behavior (constant viscosity) comes later, when you specify the constitutive relation.

Thought Experiments

1. Why matters

Ask: “What if I summed all the internal forces on all particles and got a non-zero result? What would happen?”

Answer: The center of mass would accelerate without an external force—violating Newton’s first law. So action–reaction must hold. And when it does, internal forces create local stresses but don’t change the total momentum.”

2. Pressure vs. viscous stress

  1. Pressure: acts equally in all directions; no preferred axis.
  2. Viscous stress: depends on velocity gradients; has a preferred direction (the direction of shear).

Imagine honey being poured. Gravity provides a pressure gradient (accelerating the fluid downward), and viscosity creates shear stress (resisting the motion between fast-flowing and slow-flowing regions).

Homework / Conceptual Questions

  1. Coarse-graining in practice: explain what “coarse-graining” means and why it lets us treat soft matter as a continuum. Pick one system (foam, polymer solution, emulsion) and say when this breaks down (give a specific length scale or condition).
  2. Continuity in a nozzle: You squeeze shampoo through a pump and narrow outlet at steady rate. Using only words, explain why the fluid speeds up in the narrow section and name the conservation law you’re invoking. Give one case where this reasoning would fail.
  3. Momentum balance around a bend: Water flows steadily through a pipe that turns 90 degrees. Describe which forces turn the flow (pressure differences versus wall shear) and how the picture changes for very slow flows (creeping) versus very fast flows (inertial).
  4. Where does the work go? You slide two plates with glycerol between them at constant speed. Explain where the mechanical power you input ends up inside the fluid, name the phenomenon responsible (in words), and give a one-sentence microscopic picture for it.

References & Further Reading

  1. Chapman & Cowling (1970): The Mathematical Theory of Non-Uniform Gases. Section 2: kinetic theory and microscopic description.
  2. Bird, Stewart & Lightfoot (2007): Transport Phenomena. Chapter 1: detailed derivation of continuum equations from molecular dynamics.
  3. Irving & Kirkwood (1950): “The Statistical Mechanical Theory of Transport Processes.” Journal of Chemical Physics. (Original microscopic momentum-flux definition; slightly technical but canonical.)
  4. Doi & Edwards (1986): The Theory of Polymer Dynamics. Chapter 2: stress tensor and its microscopic origin.

Quick Assessment (Check your Understanding)

After this lecture, students should be able to:

  1. Newton’s second law for a single particle and for N particles.
  2. Microscopic fields (, ) using Dirac deltas and explain why deltas are the right tool.
  3. Explain the moving delta identity and use it to derive the microscopic continuity equation.
  4. Explain how pair forces (Newton’s third law) turn into a divergence of the virial tensor.
  5. Derive the microscopic momentum balance and identify the kinetic and configurational components of stress.
  6. Describe ensemble averaging and explain why it’s necessary to transition from spiky microscopic fields to smooth macroscopic fields.
  7. Derive the Cauchy momentum equation from the microscopic balance and identify the role of the Cauchy stress tensor.
  8. Explain the physical meaning of the stress tensor and give examples (pressure, viscous stress).
  9. Understand that the Cauchy equation is exact for any fluid; the Newtonian assumption comes later (constitutive relation).

Connection to Soft Matter

  1. These laws apply to flows of complex fluids (polymer solutions obey continuity and momentum conservation, with extra stresses)
  2. Apply to deformable solids (elastic energy is stored energy in the energy equation)
  3. Foundation for understanding:
    1. Fluid instabilities (next lecture)
    2. Rheology (later lectures)
    3. All continuum modeling of soft matter

Supporting Detailed Notes

These notes contain the mathematical details and alternative proofs that are too long for the main lecture but essential for deep understanding:

  1. Moving-Delta-identity: Three detailed proofs of (chain rule, distributions, 1D sanity check)
  2. Dyadic-conversion: The algebraic steps to show
  3. Virial-identity: Full derivation of and its application to pair forces

Back to main website

Home, Team, Research, Github

📝 Edit this page on GitHub

Graph View

Backlinks