为何该非零旋度动态系统的轨迹为闭合轨道？

Great question—this is a neat example that highlights an important misconception: non-zero curl doesn't rule out closed trajectories entirely. Let's break down what's happening here.

First: Curl ≠ No Closed Trajectories #

First, let's clear up a common myth: A vector field having non-zero curl means it's not conservative (so not all closed paths have zero circulation), but it doesn't mean no closed paths can exist. Closed trajectories can still appear if the field has specific symmetries, conserved quantities, or if the initial conditions land the trajectory on a closed invariant set.

Looking at Your System's Specifics #

Your system is:
$$\dot{r}=F(r,\phi)=v \cos{\phi}$$
$$\dot{\phi}=G(r,\phi)=\frac{0.5 r^5 (\pi/2-\phi)}{v \cos{\phi}}-\frac{v \sin{\phi}}{r}$$

Let's focus on your initial condition: (r(0)=0.9, \phi(0)=\pi/2). At this point, (\cos\phi=0), so (\dot{r}=0) (radial velocity is zero) initially. The (\dot{\phi}) term needs careful handling here (since we have a 0 in the denominator), but numerically, we're looking at a point just below (\phi=\pi/2) (to avoid the singularity), where (\cos\phi) is a small positive number.

Key Observations for Closed Trajectories #

1. Conserved Quantity (First Integral):
   The most likely reason for closed trajectories here is the existence of a first integral (a function (H(r,\phi)) that's constant along trajectories, i.e., (\frac{dH}{dt}=0)). Even if the field has non-zero curl, a conserved quantity can trap trajectories on closed contour lines of (H).

   For your system, rearranging the ODEs to get the slope (\frac{d\phi}{dr}) gives:
   $$\frac{d\phi}{dr} = \frac{G(r,\phi)}{F(r,\phi)} = \frac{0.5 r^5 (\pi/2-\phi)}{v^2 \cos^2\phi} - \frac{\tan\phi}{r}$$
   While solving this exactly is messy, your numerical result of a closed trajectory implies there's some (H(r,\phi)) that stays constant—so the trajectory loops around a closed contour of (H).

2. Circulation vs. Pointwise Curl:
   Stokes' theorem tells us that the circulation around a closed path equals the integral of curl over the enclosed area. Just because curl is non-zero at some points doesn't mean the total circulation around your trajectory is non-zero. For your specific trajectory, the positive and negative contributions of curl over the enclosed area could cancel out, resulting in zero circulation—allowing the path to close.

3. Initial Condition Specificity:
   Your initial condition lands the trajectory in a region where the interplay between radial and angular motion is perfectly balanced. When (\phi) approaches (\pi/2) again, the radial velocity drops back to zero, and the system has enough "balance" in its motion to return to the initial (r) value, closing the loop.

Wrapping Up #

To sum it up: Non-zero curl means the field isn't conservative, but it doesn't prevent closed trajectories from existing. In your case, the system likely has a hidden conserved quantity that traps the trajectory on a closed contour, and the specific initial condition you chose falls exactly on one of these closed invariant sets.

内容的提问来源于stack exchange，提问作者Shengkai Li