Nonlinear PDEs arise when the equation depends nonlinearly on the unknown function or its derivatives, breaking the superposition principle that makes linear theory tractable. They are classified by degree of nonlinearity: semilinear (nonlinear in u but linear in highest derivatives, like -Δu = u³), quasilinear (coefficients of highest derivatives depend on u, like div(|∇u|^{p-2}∇u) = 0), and fully nonlinear (arbitrary dependence on highest derivatives, like det(D²u) = f). Each class requires different techniques, and fundamentally new phenomena emerge: blow-up in finite time, multiple solutions, pattern formation, and solitons.
Nonlinear PDEs are ubiquitous in science and mathematics: fluid mechanics (Navier-Stokes, Euler), general relativity (Einstein), quantum mechanics (nonlinear Schrodinger), materials science (Ginzburg-Landau), and geometry (Ricci flow, minimal surfaces) all involve nonlinear PDEs. The loss of the superposition principle means that the powerful linear machinery—Green's functions, Fourier transforms, spectral decompositions—cannot be applied directly. Instead, nonlinear PDE theory relies on a combination of a priori estimates, fixed-point theorems, variational methods, and comparison principles.
Semilinear equations, where the nonlinearity affects only the lower-order terms, are the most tractable. For -Δu = f(u) on a bounded domain, the linear theory for -Δ provides a framework (Green's functions, Sobolev regularity), and the nonlinearity f(u) is handled as a perturbation. The Leray-Schauder fixed-point theorem and the mountain pass theorem are key existence tools. The critical issue is the growth rate of f: subcritical nonlinearities (growing slower than u^{(n+2)/(n-2)} in the Sobolev-critical sense) allow compact embedding arguments, while supercritical nonlinearities can produce non-existence or blow-up.
Quasilinear equations, where the highest derivatives have coefficients depending on u or ∇u, present deeper challenges because the equation itself changes character as the solution evolves. The p-Laplacian div(|∇u|^{p-2}∇u) = 0 is degenerate elliptic (the coefficient |∇u|^{p-2} vanishes where ∇u = 0), and its solutions are C^{1,α} but generally not C^2. The regularity theory of De Giorgi-Nash-Moser provides Holder continuity for divergence-form quasilinear equations under very general structural conditions, using an iterative argument that bootstraps L^∞ bounds into Holder estimates.
Fully nonlinear equations like the Monge-Ampere equation det(D²u) = f and Hamilton-Jacobi equations H(x, Du) = 0 require entirely different frameworks—viscosity solutions for first-order equations, and the Evans-Krylov theory for second-order convex equations. The interplay between existence, uniqueness, regularity, and blow-up makes nonlinear PDE theory one of the richest and most active areas of modern mathematics. Many fundamental questions remain open: global regularity for Navier-Stokes (a Millennium Prize problem), singularity formation in the Euler equations, and the full regularity theory for the Monge-Ampere equation in non-convex settings.