First Order Partial Differential Equations Examples

First Order Partial Differential Equations Examples: A Practical Walkthrough

You open a textbook on partial differential equations and immediately hit a wall. That's where most people check out — and it's a shame, because first order PDEs are genuinely useful once you see them in action. That said, pages of symbols, boundary conditions, and what looks like abstract algebra with no clear purpose. The trick is working through enough first order partial differential equations examples to build intuition. That's what this article is for.

Whether you're a student grinding through a math course, an engineer modeling fluid flow, or someone who stumbled into this topic out of curiosity, this guide breaks down what first order PDEs actually are, why they matter, and — most importantly — gives you concrete examples you can follow step by step.

What Are First Order Partial Differential Equations

A partial differential equation is any equation that involves an unknown function of multiple variables and its partial derivatives. When the highest derivative appearing in the equation is first order — meaning no second derivatives, no third derivatives — you've got a first order PDE.

The general form looks like this:

a(x, y) * ∂u/∂x + b(x, y) * ∂u/∂y = c(x, y, u)

Here, u is the unknown function of x and y, and a, b, and c are given functions. That's it. So no second derivatives, no integrals lurking in the background. Just first partial derivatives of u with respect to each independent variable.

But don't let the simplicity fool you. Even this stripped-down form produces a rich variety of behaviors depending on what a, b, and c look like and what boundary conditions you impose.

The Difference Between ODEs and PDEs

An ordinary differential equation deals with a function of one variable and its derivatives. Consider this: a PDE deals with functions of multiple variables. That extra dimension changes everything. Solutions to PDEs are typically functions themselves — surfaces, fields, distributions — rather than curves. And the methods for solving them are fundamentally different from what you use with ODEs.

Why Studying Examples Matters

Here's the honest truth: reading about PDEs without working through examples is like reading about swimming without ever getting in the water. The notation becomes familiar, but the actual problem-solving skill doesn't develop.

When you sit down with first order partial differential equations examples, you start to notice patterns. The same technique keeps showing up under different disguises. Different physical situations reduce to the same mathematical structure once you strip away the context. That recognition — that's what turns a student into someone who can actually solve these things.

Examples also expose you to the common traps. Sign errors, misapplied boundary conditions, confusing characteristics with solution curves — these mistakes are predictable, and the best way to learn to avoid them is by making them yourself on practice problems.

How First Order PDEs Work

The General Structure and What Each Piece Means

Let's look at the general first order PDE again:

a(x, y) * ∂u/∂x + b(x, y) * ∂u/∂y = c(x, y, u)

Each term plays a specific role. The left side is a directional derivative of u in the direction defined by the vector (a, b). Think of it this way: you're asking, "how does u change as you move in a particular direction through the (x, y) plane?" The right side, c, tells you what that rate of change equals.

The function u(x, y) is what you're solving for. It could represent temperature, concentration, pressure, density — whatever physical quantity depends on position in a 2D domain.

The Method of Characteristics

The primary technique for solving first order PDEs is the method of characteristics. It converts the PDE into a system of ordinary differential equations along special curves called characteristic curves.

The idea is deceptively simple. Instead of trying to solve the PDE over the entire (x, y) plane at once, you trace paths — the characteristics — along which the PDE reduces to an ODE. Solve the ODE along each path, then piece the solutions together.

The characteristic equations are:

dx/dt = a(x, y), dy/dt = b(x, y), du/dt = c(x, y, u)

You solve this system of ODEs with appropriate initial conditions, and the solution u(x, y) emerges from that process.

Common First Order PDE Examples

This is the heart of the article. Let's go through several canonical examples, each illustrating a different flavor of first order PDE and a different solving approach.

The Linear Transport Equation

The transport equation is the simplest nontrivial first order PDE and the one you'll encounter first in most courses. It looks like:

∂u/∂t + c * ∂u/∂x = 0

Here, c is a constant representing the speed of transport. Day to day, the variable t is usually time and x is a spatial coordinate. This equation models a quantity u being carried along without changing shape — think of a wave moving through a medium, or a pollutant drifting in a river.

The characteristic curves are straight lines: x - ct = constant. Along each of these lines, u stays constant. So if you know the initial condition u(x, 0) = f(x), the solution is simply u(x, t) = f(x - ct).

That's remarkably clean. Here's the thing — the initial profile f just slides to the right at speed c without distorting. This is the linear transport equation — "linear" because u and its derivatives appear without being multiplied together or composed with nonlinear functions.

For more on this topic, read our article on total surface area of right circular cylinder or check out how do you find the height of an obtuse triangle.

The Inviscid Burgers' Equation

Now make it nonlinear. The inviscid Burgers' equation is:

∂u/∂t + u * ∂u/∂x = 0

The key difference: the "speed" of transport now depends on u itself. Here's the thing — where u is large, the wave moves fast. Where u is small, it moves slow. This creates a feedback loop that eventually causes the solution to steepen and form a shock — a discontinuity where u jumps abruptly.

The method of characteristics still applies, but now the characteristic curves are not straight lines. They're curves whose slopes depend on the value of u, and they can cross. When characteristics cross, the classical solution breaks down, and you need to introduce a shock condition to continue the solution past that point.

This example is a gateway to nonlinear PDE theory. It shows that even first order equations can produce genuinely complex behavior — shocks, rarefaction waves, multivalued solutions that need to be resolved by physical reasoning.

Clairaut's Equation

Clairaut's equation has a special form:

u = x * ∂u/∂x + y * ∂u/∂y + f(∂u/∂x, ∂u/∂y)

where f is some given function of the two partial derivatives. This looks unusual because the unknown function u appears on the left and its

Clairaut’s equation is noteworthy because it admits both a family of straight‑line solutions and, generically, a singular envelope that cannot be obtained by fixing the constants in the general solution. Writing (p = u_x) and (q = u_y), the equation reads

[ u = xp + yq + f(p,q). ]

Treating (p) and (q) as parameters, the characteristic system derived from the Lagrange‑Charpit equations becomes

[ \frac{dx}{dt}=x+f_p,\qquad \frac{dy}{dt}=y+f_q,\qquad \frac{du}{dt}=p\bigl(x+f_p\bigr)+q\bigl(y+f_q\bigr),\qquad \frac{dp}{dt}=0,\qquad \frac{dq}{dt}=0, ]

so that (p) and (q) remain constant along each characteristic. Integrating the first two equations gives

[ x(t)=\bigl(x_0+f_p\bigr)e^{t}-f_p,\qquad y(t)=\bigl(y_0+f_q\bigr)e^{t}-f_q, ]

and substituting these into the expression for (u) yields the general solution

[ u(x,y)=xf(p,q)+yf(q,p)+f(p,q), ]

where ((p,q)) are arbitrary constants. Eliminating (p,q) from the pair

[ x+f_p=0,\qquad y+f_q=0 ]

produces the singular solution, which is the envelope of the family of straight lines and often corresponds to a physically relevant branch (for instance, the caustic in geometric optics).

Hamilton–Jacobi and Eikonal Equations

Another important class of first‑order PDEs appears in optimal control and wave propagation:

[ u_t + H\bigl(x,\nabla u\bigr)=0, ]

known as the Hamilton–Jacobi equation. When the Hamiltonian takes the form (H(x,p)=|p|) (or a variable‑speed version (c(x)|p|)), the equation reduces to the eikonal equation

[ |\nabla u| = \frac{1}{c(x)}. ]

Its characteristics are the rays of geometric optics; solving the characteristic ODEs yields the travel‑time function (u(x)) as the minimum‑action principles. The method again furnishes a complete integral, and singularities (caustics) arise when neighboring characteristics intersect, signalling the breakdown of the classical solution and the need for viscosity‑type selections.

Conservation Laws and Riemann Problems

Many physical models lead to first‑order PDEs in conservation form

[ u_t + \nabla!\cdot!F(u)=0, ]

with a flux (F:\mathbb{R}^m\to\mathbb{R}\to\mathbb{R}^R). Consider this: the one‑dimensional scalar case (u_t+F(u)_x=0) is tackled by the Rankine‑Hugoniot jump condition at discontinuities and the entropy condition to select the physically admissible shock. The Riemann problem—piecewise constant initial data with a single jump—produces elementary waves (shocks, rarefactions, contact discontinuities) whose interaction builds up general solutions via the wave‑front tracking or Godunov schemes.

Conclusion

First‑order PDEs, though seemingly simple, encapsulate a rich variety of phenomena ranging from passive transport and nonlinear steepening to wavefront propagation and optimal control. The unifying tool is the method of characteristics, which converts the PDE into a system of ordinary differential equations along curves in the independent‑variable space. When these characteristics remain distinct, the solution is obtained by propagating initial data unchanged (or with a known factor).

...and entropy conditions, or, in the context of Hamilton–Jacobi theory, viscosity solutions. These selection principles not only restore well-posedness but also reflect the underlying physics: shocks model irreversible dissipation, while viscosity solutions capture the asymptotic behavior of diffusive regularizations.

Beyond the classical examples, first-order PDEs permeate diverse domains. In economics, Hamilton–Jacobi equations arise in optimal investment strategies; in aerodynamics, conservation laws govern compressible flow; and in image processing, the eikonal equation drives contour evolution algorithms. Mastery of the method of characteristics, coupled with an appreciation for its limitations, equips analysts to figure out both smooth regimes and the singular frontiers where classical differentiability fails.

Thus, while the formalism may appear elementary, the theory of first-order PDEs offers a profound lens through which to view and predict the evolution of complex systems, blending geometric intuition with rigorous analytical tools to illuminate the interplay between continuity and discontinuity in nature’s tapestry.