TheoremComplete

The Bendixson-Dulac Criterion

The Bendixson-Dulac criterion provides a sufficient condition for the non-existence of periodic orbits, complementing the Poincare-Bendixson existence theorem.


Statement

Theorem6.6Bendixson-Dulac theorem

Let DR2D \subseteq \mathbb{R}^2 be a simply connected domain and F=(f,g)C1(D)\mathbf{F} = (f, g) \in C^1(D). If there exists a C1C^1 function B:DRB: D \to \mathbb{R} (called a Dulac function) such that

(Bf)x+(Bg)y\frac{\partial(Bf)}{\partial x} + \frac{\partial(Bg)}{\partial y}

does not change sign in DD and is not identically zero on any open subset, then the system x=f(x,y)x' = f(x,y), y=g(x,y)y' = g(x,y) has no periodic orbits entirely contained in DD.


Proof

Proof

Suppose for contradiction that γ\gamma is a periodic orbit in DD. Let RR be the region enclosed by γ\gamma (which exists since DD is simply connected). By Green's theorem:

R[(Bf)x+(Bg)y]dA=γ(BfdyBgdx).\iint_R \left[\frac{\partial(Bf)}{\partial x} + \frac{\partial(Bg)}{\partial y}\right]dA = \oint_\gamma (Bf\,dy - Bg\,dx).

On γ\gamma: dx=fdtdx = f\,dt and dy=gdtdy = g\,dt, so

γ(BfdyBgdx)=γB(fggf)dt=0.\oint_\gamma (Bf\,dy - Bg\,dx) = \oint_\gamma B(fg - gf)\,dt = 0.

But the double integral is nonzero (since the integrand has one sign and is not identically zero). Contradiction. \blacksquare


Examples

ExampleSimple application

For x=y+x(x2+y21)x' = -y + x(x^2+y^2-1), y=x+y(x2+y21)y' = x + y(x^2+y^2-1): take B=1B = 1.

fx+gy=(3x2+y21)+(x2+3y21)=4(x2+y2)2.\frac{\partial f}{\partial x} + \frac{\partial g}{\partial y} = (3x^2+y^2-1) + (x^2+3y^2-1) = 4(x^2+y^2) - 2.

This changes sign (positive for x2+y2>1/2x^2+y^2 > 1/2, negative for x2+y2<1/2x^2+y^2 < 1/2), so B=1B = 1 does not work directly. A different Dulac function is needed.

However, in polar coordinates: r=r(r21)r' = r(r^2-1) has a limit cycle at r=1r = 1, so the criterion correctly cannot rule it out.

ExampleChoosing a Dulac function

For x=yx' = y, y=xy+x2y' = -x - y + x^2: try B(x,y)=e2yB(x,y) = e^{2y}.

(Be2yy)x+(Be2y(xy+x2))y=0+e2y[2(xy+x2)+(1)]=e2y(2x22x2y1).\frac{\partial(Be^{2y}\cdot y)}{\partial x} + \frac{\partial(Be^{2y}(-x-y+x^2))}{\partial y} = 0 + e^{2y}[2(-x-y+x^2) + (-1)] = e^{2y}(2x^2-2x-2y-1).

For the region where 2x22x2y1<02x^2-2x-2y-1 < 0, no periodic orbits exist.

RemarkGradient and Hamiltonian systems

Gradient systems x=V(x)\mathbf{x}' = -\nabla V(\mathbf{x}) never have periodic orbits (along any trajectory, V˙=V20\dot{V} = -|\nabla V|^2 \leq 0, so VV decreases and cannot return to its initial value). Hamiltonian systems x=H/yx' = \partial H/\partial y, y=H/xy' = -\partial H/\partial x satisfy div(F)=0\text{div}(\mathbf{F}) = 0 (Bendixson with B=1B=1), but this does not rule out periodic orbits — instead, Hamiltonian systems typically have families of periodic orbits.